ISM模型的matlab实现_第1页
ISM模型的matlab实现_第2页
ISM模型的matlab实现_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、ISM模型的matlab实现1、可达矩阵的求解A：邻接矩阵E=zeros（A）；通过布尔运算求出A2、A3直至An=An-1,停止运算，此时的An便是要求的可达矩阵。具体实施：通过将上次循环得出来的A赋给E，并判断新一轮循环得出的A是否与E相等，决定是否终止循环。循环终止意味着此时的A就是最终的可达矩阵。具体程序如下：n=input(请输入矩阵维数:);A=input(请输入邻接矩阵:);E=zeros(n);B=A;while(norm(A-E)0) E=A;for i=1:n for j=1:n for k=1:n if A(i,k)&B(k,j) A(i,j)=1; end end en

2、dendendA2、区域划分找出各个元素相对应的可达集P、先行集Q以及两者的交集S求解P：找出每一行中元素为1对应的列求解Q：找出每一列中元素为1对应的行求解S：套用Matlab本身自带的函数S=intersect(P,Q); 或者编写M文件函数，如下所示：具体程序如下：for i=1:nP=find(A(i,:);Q=find(A(:,i);S=intersect(P,Q); P; Q; S;3、级别划分因为S是P与Q的交集，所以只需判断P与S的长度是否相等便可进行级别划分。M=zeros（n）；r=1；r为第几级求出每个元素的P、Q、S，再将相应的对角线元素A（i，i）赋予0，通过查找对角

3、线为0的元素位置将所对应的行列均赋值0。每次循环r自增1，再进行循环，求出每一级的元素，直至A=M。具体程序如下：r=1;M=zeros(n);while(isequal(A,M)for i=1:nP=find(A(i,:);Q=find(A(:,i);S=intersect(P,Q); P; Q; S;if(isempty(P)&isempty(Q)&( length(P)=length(S) disp(第r级:) r disp(元素为) i A(i,i)=0; endendfor i=1:nif A(i,i)=0A(i,:)=0;A(:,i)=0;endend r=r+1;end4.课本例题程序验证请输入矩阵维数:7请输入邻接矩阵: 1 1 1 0 1 1 0 0 1 0 0 1 1 0 0 1 1 0 1 1 0 0 1 1 1 1 1 0 0 0 0 0 1 0 0 0 1 0 0 1 1 0 0 1 1 0 1 1 1第r级:r = 1,元素为i =5第r级:r =2,元素为i =2第r级:r = 2,元素为i =6第r级:r =3,

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论