系统工程ISMmatlab代码

1、系统工程作业本次作业在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

2、 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

3、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 求综合影响矩阵中的行和、列和、中心度

4、、原因度：在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)

5、=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

6、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如下：计算系统的可达性集合和先行集合，以及二者的共同集合，如下表所示：iL（si）F(si)L(si)F(si)111-13121,22,7231,33,7341,44-9,11,1

7、3451,4-6,10-125,11,135,1161,4,65,6,11,13671-4,7-97781,4,87,8,13891,4,97,9,139101,10,125,10,11,1310111,4-6,10-125,11,135,11121,125,10-1312131,4-6,8-131313由上表可知: T=7,13，由于L(s3)L(s7),所以s3、s7属于同一个区域中，即M不可分区。 s1是第一级要素。划去s1对应的行列，所得表如下：iL（si）F(si)L(si)F(si)222,72333,73444-9,11,13454-6,10-125,11,135,1164,65

8、,6,11,13672-4,7-97784,87,8,13894,97,9,1391010,125,10,11,1310114-6,10-125,11,135,1112125,10-1312134-6,8-131313由上表可知：s2、s3、s4、s12是第二级要素。同理可得：s6、s8、s9、s10是第三级要素；s5、s7、s11是第四级要素；s13是第五级要素。按照上面所划分的要素等级，对可达矩阵M进行重新排列，得到如下：第三题：对于上述问题研究的系统，利用DEMATEL方法判断各要素之间的综合影响关系，并求出各要素的原因度和中心度。计算规范化直接影响矩阵：继续算得综合影响矩阵各要素之间的

9、综合影响矩阵要 素s1s2s3s4s5s6s7s8s9s10s11s12s13s10000000000000s20.1667000000000000s30.1667000000000000s40.1667000000000000s50.2397000.17620.02860.02860000.02860.17140.00480s60.1944000.1667000000000s70.28700.16670.16670.0556000016670.16670000s80.1944000.1667000000000s90.1944000.1677000000000s100.194400000000

10、000.16670s110.2714000.05710.17140.17140000.17140.02860.02860s120.1667000000000000s130.3362000.11270.17140.171400.16670.16670.17140.02860,02860各系统要素的原因度以及中心度要素行和列和原因度中心度s102.5787-2.57872.5787s20.16670.166700.3334s30.16670.166700.3334s40.16670.9027-0.7361.0694s50.67790.37140.30651.0493s60.36110.3714-0

11、.01030.7325s71.009401.00941.0094s80.36110.33340.02770.6945s90.36210.33340.02870.6955s100.36110.3714-0.01030.7325s110.89990.22860.67131.1285s120.16670.2287-0.0620.3954s131.353701.35371.3537从上表的结果中可以得到：各要素在系统中的重要程度依次是：s1、s13、s11、s4、s5、s7、s6、s10、s9、s8、s12、s2、s3；系统中原因要素是：s1、s4、s6、s10、s12；系统中结果要素是：s5、s7、

12、s8、s9、s11、s13在该系统中，s2、s3的原因度最低，其中心度也最低，可以考虑删除该要素，达到减少要素的目的；第四题：求解问题二中的系统结构模型。观察矩阵，不难发现，s5与s11的相应的行列元素完全相同，可以把二者当成一个系统看待，从而可缩减相应的行和列，得到缩减矩阵如下：在矩阵中，先找出一级与二级之间的关系，再找出二级与三级之间的关系，直到把各级找完为止。从中知=1，说明节点s2与处于第一级的节点s1有关，即s2s1，以此类推，则有：第一级：s1第二级：s2s1，s3s1，s4s1，s12s1；第三级：s6s4，s8s4，s9s4，s10s12；第四级：s5s6，s5s10，s7s2

13、，s7s3，s7s8，s7s9；第五级：s13s5。根据上述画出结构矩阵如下：从而，可以绘制出系统的多级递阶有向结构图，如下图所示：第五题：探讨ISM方法与DEMATEL方法的区别与联系。区别：DEMATEL方法是由美国学者提出的一种运用图论和矩阵论原来进行系统因素分析的方法，它借助系统中各因素之间的逻辑关系构建直接影响矩阵，计算各因素对其他因素的影响程度以及被影响程度，从而计算各因素的中心度和原因度。根据因素所对应的中心度和原因度，得出该因素所属种类（原因因素还是结果因素），也可根据中心度和原因度的取值掉正系统的结构图，使系统的结构更加合理。ISM方法是现代系统工程学中广泛应用的一种分析和揭

14、示系统结构的方法。ISM作为一种分析系统结构的方法，将系统要素之间复杂、凌乱的关系分解成清晰的、多层级的结构形式。以ISM法得到的系统结构是一种宏观的定性结构，它揭示系统的几何学的定性结构，而不是对其结构做出精确的代数描述，或给出数量上、统计上的性质。也就是说。ISM法在分析系统结构时，是以系统的各要素为研究对象，以各要素间是否存在某种关系（从属关系、并列关系、因果关系等）为构造模型的唯一依据，即如果两要素间存在着某种关系，在模型中它们之间就会有连线相接，否则两个要素就是独立的。可见，DEMATEL方法以要素间的影响程度为依据，以分析各要素的影响度、被影响度、中心度和原因度为核心，它不仅反映了

15、系统各要素间的相互影响关系及相应的影响程度，而且反映了各要素在系统中的重要程度；ISM方法是以分析系统要素间的联系（形成关系）为基础，以分析系统层级化结构为核心，并可以进一步进行系统的序列化、聚类分析等，因此它不需要具体因素间的数量化关心。联系：从数学模型可以看出，DEMATEL方法和ISM方法并不是两个彼此相互独立的方法，两个基本矩阵的转化是这两种方法相结合进行系统分析的基础。DEMATEL方法中的直接影响矩阵可以转化为ISM方法中的邻接矩阵。DEMATEL方法中的综合影响矩阵可以转化为ISM方法中的可达矩阵。DEMATEL方法与ISM方法各有侧重，但是研究的每一方面都是系统分析的重要组成部

16、分。第六题：食物链的分析求得A的可达矩阵M如下：计算系统的可达性集合和先行集合，以及二者的共同集合，如下表所示：iL（si）F(si)L(si)F(si)11-121122,11,121,2233,10-121,3344,11,121,4455-121,5566-121,5,6677-121,5-7788,10,121,5-8899,11,121,5-7,991010,121,3,5-8,101011111-7,9,111112121-10,1212由上表可知: T=1，由于L(s1)L(si),所以s1与其他要素属于同一个区域中，即M不可分区。 s11、s12是第一级要素。划去s12对应的行

17、列，所得表如下：iL（si）F(si)L(si)F(si)11-1011221,2233,101,33441,4455-101,5566-101,5,6677-101,5-7788,101,5-88991,5-7,9910101,3,5-8,1010由上表可知：s2、s4、s9、s10是第二级要素。同理可得：s3、s8是第三级要素； s7是第四级要素；s6是第五级要素；s5是第六级要素；s1是第七级要素。按照上面所划分的要素等级，对可达矩阵M进行重新排列，得到如下：对于食物链的例子，用DEMATEL方法，算的综合影响矩阵T如下：要 素s1s2s3s4s5s6s7s8s9s10s11s12s10

18、0.25000.25000.25000.25000.06250.07810.08200.09770.08300.21190.2327s200000000000.25000.2500s30000000000.25000.25000.3125s400000000000.25000.2500s5000000.25000.31250.32810.39060.08200.09770.1182s60000000.25000.06250.31250.01560.07810.0820s700000000.25000.25000.06250.06250.0781s80000000000.250000.0625s

19、900000000000.25000.2500s10000000000000.2500s11000000000000s12000000000000各系统要素的原因度以及中心度要素行和列和原因度中心度s11.8479 0.0000 1.8479 1.8479 s20.5000 0.2500 0.2500 0.7500 s30.8125 0.2500 0.5625 1.0625 s40.5000 0.2500 0.2500 0.7500 s51.5791 0.2500 1.3291 1.8291 s60.8008 0.3125 0.4883 1.1133 s70.7031 0.6406 0.0625 1.3438 s80.3125 0.7227 -0.4102