系统工程ismmatlab代码

系统工程作业本次作业在 matlab 中的代码如下： 求可达矩阵 M：在 matlab 中输入代码如下：A=0 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 1 0 0 0 0 0 0 1 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 1 1 0 0 0 0 1 1 0 0 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 1 0;1 0 0 0 1 1 0 0 0 1 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 1 1 0 1 1 1 0 0 0;I=eye(13);M=(A+I)13;M(M=0)=1;M 求综合影响矩阵 T：在 matlab 中输入代码如下：A=0 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 1 0 0 ;1 0 0 1 0 0 0 0 0 0 0 0 0 ;1 1 1 0 0 0 0 1 1 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 1 0 ;1 0 0 0 1 1 0 0 0 1 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 1 1 0 1 1 1 0 0 0 ;d=0;s=0;for i=1:13for j=1:13s=s+A(i,j);if dsd=s;endends=0;endG=A/d;I=eye(13);T=G/(I-G)T 求综合影响矩阵中的行和、列和、中心度、原因度：在 matlab 中代码如下：x=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13for j=1:13x(i,1)=x(i,1)+T(i,j);endendx（行和）y=0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0;for i=1:13for j=1:13y(1,i)=y(1,i)+T(j,i);endendy（列和）z=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13z(i,1)=x(i,1)+y(1,i);endz（中心度）c=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13c(i,1)=x(i,1)-y(1,i);endc（原因度）食物链的例子代码基本相同，就是矩阵输入不同，再次则不以此写出，只需将 A 矩阵替换为 0 1 1 1 1 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 0 0 0 0 1 1 1;0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 1 1 1 1 0 0 0;0 0 0 0 0 0 1 0 1 0 0 0;0 0 0 0 0 0 0 1 1 0 0 0;0 0 0 0 0 0 0 0 0 1 0 0;0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 0 0 0 0 0 0 1;0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0;即可。其他代码一致。具体题目及其答案如下第二题：一个系统的邻接矩阵为 A，求 A 的可达矩阵，并对其进行分解。求得 A 的可达矩阵 M 如下：计算系统的可达性集合和先行集合，以及二者的共同集合，如下表所示：i L（si） F(si) L(si)F(si)1321098765432113209876543A 01010000101 0100101 00000 32098765432132098765431M 11010001011 000 1011 00001 1 1-13 12 1,2 2,7 23 1,3 3,7 34 1,4 4-9,11,13 45 1,4-6,10-12 5,11,13 5,116 1,4,6 5,6,11,13 67 1-4,7-9 7 78 1,4,8 7,8,13 89 1,4,9 7,9,13 9101,10,125,10,11,13 1011 1,4-6,10-12 5,11,13 5,1112 1,12 5,10-13 1213 1,4-6,8-13 13 13由上表可知: T=7,13，由于 L(s3)L(s7),所以 s3、 s7 属于同一个区域中，即 M 不可分区。 s1 是第一级要素。划去 s1 对应的行列，所得表如下：i L（si） F(si) L(si)F(si)2 2 2,7 23 3 3,7 34 4 4-9,11,13 45 4-6,10-12 5,11,13 5,116 4,6 5,6,11,13 67 2-4,7-9 7 78 4,8 7,8,13 89 4,9 7,9,13 910 10,12 5,10,11,13 1011 4-6,10-12 5,11,13 5,1112 12 5,10-13 1213 4-6,8-13 13 13由上表可知：s2、s3、s4、 s12 是第二级要素。同理可得：s6、s 8、s 9、s10 是第三级要素；s5、s7、s 11 是第四级要素；s13 是第五级要素。按照上面所划分的要素等级，对可达矩阵 M 进行重新排列，得到 如下：M第三题：对于上述问题研究的系统，利用 DEMATEL 方法判断各要素之间的综合影响关系，并求出各要素的原因度和中心度。计算规范化直接影响矩阵：继续算得综合影响矩阵各要素之间的综合影响矩阵要 素s1 s2 s3 s4 s5 s6 s7s8 s9 s10 s11 s12 s13s10 0 0 0 0 0 0 0 0 0 0 0 0s20.16670 0 0 0 0 0 0 0 0 0 0 0137510986124321137510986124M 10111000110000100000 AaGjii61mx13113i)(ITs30.16670 0 0 0 0 0 0 0 0 0 0 0s40.16670 0 0 0 0 0 0 0 0 0 0 0s50.23970 0 0.17620.02860.02860 0 0 0.02860.17140.00480s60.19440 0 0.16670 0 0 0 0 0 0 0 0s70.28700.16670.16670.05560 0 0 016670.16670 0 0 0s80.19440 0 0.16670 0 0 0 0 0 0 0 0s90.19440 0 0.16770 0 0 0 0 0 0 0 0s100.19440 0 0 0 0 0 0 0 0 0 0.16670s110.27140 0 0.05710.17140.17140 0 0 0.17140.02860.02860s120.16670 0 0 0 0 0 0 0 0 0 0 0s130.33620 0 0.11270.17140.17140 0.16670.16670.17140.02860,02860各系统要素的原因度以及中心度要素 行和 列和 原因度 中心度s1 0 2.5787 -2.5787 2.5787s2 0.1667 0.1667 0 0.3334s3 0.1667 0.1667 0 0.3334s4 0.1667 0.9027 -0.736 1.0694s5 0.6779 0.3714 0.3065 1.0493s6 0.3611 0.3714 -0.0103 0.7325s7 1.0094 0 1.0094 1.0094s8 0.3611 0.3334 0.0277 0.6945s9 0.3621 0.3334 0.0287 0.6955s10 0.3611 0.3714 -0.0103 0.7325s11 0.8999 0.2286 0.6713 1.1285s12 0.1667 0.2287 -0.062 0.3954s13 1.3537 0 1.3537 1.3537从上表的结果中可以得到： 各要素在系统中的重要程度依次是：s1、s13、s11 、s4、s5、s7、s 6、s10、s9、s8、s12、s 2、s3； 系统中原因要素是：s1、s4、s 6、s10、s12；系统中结果要素是：s5、s 7、s 8、s9、s11、s13 在该系统中，s2、s 3 的原因度最低，其中心度也最低，可以考虑删除该要素，达到减少要素的目的；第四题：求解问题二中的系统结构模型。观察矩阵 ，不难发现，s5 与 s11 的相应的行列元素完全相同，可以把二者当成一个M系统看待，从而可缩减相应的行和列，得到缩减矩阵 如下：Ma在矩阵 中，先找出一级与二级之间的关系，再找出二级与三级之间的关系，直到把各M级找完为止。从 中知 =1，说明节点 s2 与处于第一级的节点 s1 有关，即21ms2s1，以此类推，则有：第一级：s113751098612432113750986124Ma 10110001100 00100 3759862432137509861243IMa 0111000010000000第二级：s2s 1，s 3s1，s4s1，s 12s1；第三级：s6s 4，s 8s4，s9s4，s 10s12；第四级：s5s 6，s 5s10，s7s2，s 7s3，s7 s8，s7s 9；第五级：s13s 5。根据上述画出结构矩阵 如下：A从而，可以绘制出系统的多级递阶有向结构图，如下图所示：137509861243113750986124A 010010010第五题：探讨 ISM 方法与 DEMATEL 方法的区别与联系。区别：DEMATEL 方法是由美国学者提出的一种运用图论和矩阵论原来进行系统因素分析的方法，它借助系统中各因素之间的逻辑关系构建直接影响矩阵，计算各因素对其他因素的影响程度以及被影响程度，从而计算各因素的中心度和原因度。根据因素所对应的中心度和原因度，得出该因素所属种类（原因因素还是结果因素） ，也可根据中心度和原因度的取值掉正系统的结构图，使系统的结构更加合理。ISM 方法是现代系统工程学中广泛应用的一种分析和揭示系统结构的方法。ISM 作为一种分析系统结构的方法，将系统要素之间复杂、凌乱的关系分解成清晰的、多层级的结构形式。以 ISM 法得到的系统结构是一种宏观的定性结构，它揭示系统的几何学的定性结构，而不是对其结构做出精确的代数描述，或给出数量上、统计上的性质。也就是说。ISM 法在分析系统结构时，是以系统的各要素为研究对象，以各要素间是否存在某种关系（从属关系、并列关系、因果关系等）为构造模型的唯一依据，即如果两要素间存在着某种关系，在模型中它们之间就会有连线相接，否则两个要素就是独立的。可见，DEMATEL 方法以要素间的影响程度为依据，以分析各要素的影响度、被影响度、中心度和原因度为核心，它不仅反映了系统各要素间的相互影响关系及相应的影响程度，而且反映了各要素在系统中的重要程度；ISM 方法是以分析系统要素间的联系（形成关系）为基础，以分析系统层级化结构为核心，并可以进一步进行系统的序列化、聚类分析等，因此它不需要具体因素间的数量化关心。联系：从数学模型可以看出，DEMATEL 方法和 ISM 方法并不是两个彼此相互独立的方法，两个基本矩阵的转化是这两种方法相结合进行系统分析的基础。 DEMATEL 方法中的直接影响矩阵可以转化为 ISM 方法中的邻接矩阵。 DEMATEL 方法中的综合影响矩阵可以转化为 ISM 方法中的可达矩阵。 DEMATEL 方法与 ISM 方法各有侧重，但是研究的每一方面都是系统分析的重要组成部分。第六题：食物链的分析 120987654321A 120987654 0001000101100001求得 A 的可达矩阵 M 如下：计算系统的可达性集合和先行集合，以及二者的共同集合，如下表所示：i L（si） F(si) L(si)F(si)1 1-12 1 12 2,11,12 1,2 23 3,10-12 1,3 34 4,11,12 1,4 45 5-12 1,5 56 6-12 1,5,6 67 7-12 1,5-7 78 8,10,12 1,5-8 89 9,11,12 1,5-7,9 910 10,12 1,3,5-8,10 1011 11 1-7,9,11 1112 12 1-10,12 12由上表可知: T=1，由于 L(s1)L(si),所以 s1 与其他要素属于同一个区域中，即 M 不可分区。 s11、s12 是第一级要素。120987654321M120987654 1000110001011000011划去 s12 对应的行列，所得表如下：i L（si） F(si) L(si)F(si)1 1-10 1 12 2 1,2 23 3,10 1,3 34 4 1,4 45 5-10 1,5 56 6-10 1,5,6 67 7-10 1,5-7 78 8,10 1,5-8 89 9 1,5-7,9 910 10 1,3,5-8,10 10由上表可知：s2、s4、s9、 s10 是第二级要素。同理可得：s3、s 8 是第三级要素； s7 是第四级要素；s6 是第五级要素；s 5 是第六级要素；s1 是第七级要素。按照上面所划分的要素等级，对可达矩阵 M 进行重新排列，得到 如下：M对于食物链的例子，用 DEMATEL 方法，算的综合156783109421M1567831094 11100010001000影响矩阵 T 如下：要 素s1s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12s1 0 0.25000.25000.25000.25000.06250.07810.08200.09770.08300.21190.2327s2 0 0 0 0 0 0 0 0 0 0 0.25000.2500s3 0 0 0 0 0 0 0 0 0 0.25000.25000.3125s4 0 0 0 0 0 0 0 0 0 0 0.25000.2500s5 0 0 0 0 0 0.25000.31250.32810.39060.08200.09770.1182s6 0 0 0 0 0 0 0.25000.06250.31250.01560.07810.0820s7 0 0 0 0 0 0 0 0.25000.25000.06250.06250.0781s8 0 0 0 0 0 0 0 0 0 0.25000 0.0625s9 0 0 0 0 0 0 0 0 0 0 0.25000.2500s100 0 0 0 0 0 0 0 0 0 0 0.2500s110 0 0 0 0 0 0 0 0 0 0 0s120 0 0 0 0 0 0 0 0 0 0 0各系统要素的原因度以及中心度要素 行和 列和 原因度 中心度s1 1.8479 0.0000 1.8479 1.8479 s2 0.5000 0.2500 0.2500 0.7500 s3 0.8125 0.2500 0.5625 1.0625 s4 0.5000 0.2500 0.2500 0.7500 s5 1.5791 0.2500 1.3291 1.8291 s6 0.8008 0.3125 0.4883 1.1133 s7 0.7031 0.6406 0.0625 1.3438 s8 0.3125 0.7227 -0.4102 1.0352 s9