ISM的计算机实现.docxVIP

ISM技术的计算机实现 陈林艺 1302011006 辛全超 1302011066 学员二旅五营 摘要：本文将解析结构模型（ISM）的传统算法利用MATLAB语言进行实现，应用矩阵的运算法则，把系统中的要素及要素之间的上下位关系用邻接矩阵表示，并做出系统的关系图，使复杂的系统分解成多级递阶结构形式。 关键词：邻接矩阵、可达矩阵、级别划分、ISM 引言： ISM法(Interpretative Structural Modeling Method,)是一种分析系统结构的方法。解析结构模型可将系统单元之间复杂、凌乱的关系分解成清晰的、多级递阶的结构形式。凡系统必有结构，系统的结构决定系统功能；破坏结构，就会完全破坏系统的总体功能，这说明了系统结构的普遍性与重要性。总之，要研究一个由大量单元组成的、各单元之间又存在这相互关系的系统，就必须了解系统的结构。 ２０世纪７０年代以来解析结构模型（ISM）和其他结构模型在社会经济系统中得到了广泛的应用，解析结构模型属于静态的定性模型，建立系统的递阶结构模型，是ISM技术的核心内容。如果用传统的矩阵算法来划分系统层次，其运算量之大是让任何人都难以承受的。因此采用一种简便的算法设计开发出一个ISM法系统的程序是很有必要的。 ISM实现的算法思想 系统中的要素及要素之间的上下位关系可用对应的矩阵表示，这种矩阵称为关系矩阵，其中最直接的一种是邻接矩阵，用来表示各要素之间的连接关系。由于邻接矩阵是布尔矩阵，所以矩阵元素按布尔运算（逻辑和、逻辑与、逻辑乘）的法则进行运算，将邻接矩阵转化为可达矩阵，再将矩阵进行级间分解，然后再通过人-机结合，使复杂的系统分解成多级递阶结构形式。 ISM技术算法的功能设计 算法以MATLAB软件运行后，在kedaxingjuzhen.m文件中将系统的关系矩阵输入，例如：A = [1 0 0 0 0 0 0;1 1 0 0 0 0 0;0 0 1 1 1 1 0;0 0 0 1 1 1 0;0 0 0 0 1 0 0;0 0 0 1 1 1 0;1 1 0 0 0 0 1];运行Guanxitu.m文件后，在Command Window界面中分别输入：A（邻接矩阵）、R（可达矩阵）、L(级别划分矩阵)。在.m文件中将系统的关系矩阵输入，例如：rel=[1 0 0 0 0 0 0;1 1 0 0 0 0 0;0 0 1 1 1 1 0;0 0 0 1 1 1 0;0 0 0 0 1 0 0;0 0 0 1 1 1 0;1 1 0 0 0 0 1];运行后在Figure 1对话框中将会出现矩阵的关系图。 ISM技术的算法代码 Guanxitu.m文件代码 % ISM 解释型结构模型 根据邻接矩阵求可达矩阵，进行级划分的算法 % step 1 求可达矩阵 % A 是邻接矩阵，可以由用户输入 A = [1 0 0 0 0 0 0;1 1 0 0 0 0 0;0 0 1 1 1 1 0;0 0 0 1 1 1 0;0 0 0 0 1 0 0;0 0 0 1 1 1 0;1 1 0 0 0 0 1]; N=size(A,1); % N 是矩阵的阶数，是所有要素的数目 r = [N:N]; r=A; for(n=1:1:N) for(i=1:1:N) for(j=1:1:N) sum = 0.0; for(k=1:1:N) sum=sum+r(i,k)*A(k,j);%此处是采用普通的矩阵乘法，但是下面可以根据sum的值是否大于或等于1 %来判断连通型，从而等价得到r(i,j)的实际值。 end if(sum = 1) r(i,j)= 1; else r(i,j)=0; end end end end R=r; % R 是可达矩阵 %step 2 级别划分 L=[N:N]; % L 是二维数组，存储级划分的结果，下面首先初始化为0 for(i=1:1:N) for(j=1:1:N) L(i,j)=0; end end for(p=1:1:N) % p 是存储层次级数 l 的变量，也控制了总循环的次数 k=1; % k 是记录每一级(p)内要素，在L内存储下标的变量,形式为L(p,k) %下面是利用二重循环，求解p级内的要素 for(i=1:1:N) sign=0; % 是标志变量，初始为0，如果对某要素考察后，其值仍为0，则表明该要素是顶点要素。否则，不是顶点要素 sum=0; % sum 是计数器，用于判断当前考察的要素所在的矩阵行向量是否是值全为0的向量。因为本程序的算法是???样的， % 即如果已经发现某要素是某级的顶点，则在求下一级顶点时，需要去除上面各级的要素。在本程序是通过变通的方法实现 % 相同的目的，即所有已经是顶点的要素的所在行和列的值全