全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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 endendendA2、区域划分找出各个元素相对应的可达集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,通过查找对角线为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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024建筑工程服务合同范本
- 2023年光导纤维导光环项目风险评价报告
- 幕墙工程施工施工安排
- 2024外汇买卖合同法
- 2024单位汽车租赁合同官方
- 2023年高性能铜镍合金带、线材项目风险分析及评价报告
- 2024房屋个人拆除协议书
- 2024工程建设招标投标合同模板
- 2024建设工程设计合同
- 2023年铌酸钾晶体(KNBO3)项目安全评价报告
- 人民医院申请购置医疗设备可行性报告
- 2023年高考考点考务工作总结
- GB 5009.3-2016食品安全国家标准食品中水分的测定
- 2023年华英学校小升初数学试卷(含答案)
- 三通两平台建设方案
- 成本计划及保证措施
- 从十四五规划看数字化转型课件
- 最新NXT操作员培训资料课件
- 《雷雨》(完美版)(共28张)课件
- 《发展汉语(第二版)中级综合(Ⅰ)》第2课课件
- 心脏介入诊疗技术操作规范及流程
评论
0/150
提交评论