上海工程技术大学 《系统工程》课程论文
for(j=1:1:N) sum=0;
for(k=1:1:N)
sum=sum+r(i,k)*A(k,j); end
if(sum>=1) r(i,j)=1; else
r(i,j)=0; end end end end R=r;
L=[N:N]; for(i=1:1:N) for(j=1:1:N) L(i,j)=0; end end
for(p=1:1:N) k=1;
for(i=1:1:N) sign=0; sum=0; for(j=1:1:N)
sum=sum+r(i,j);
if(r(i,j)==1&(r(i,j)~=r(j,i))) sign=1 break end end
if(sum~=0)
if(sign==0) L(p,k)=i; K=k+1; end end end
for(g=1:1:N)
if(L(p,g)~=0) for(h=1:1:N)
r(L(p,g),h)=0;
15
上海工程技术大学 《系统工程》课程论文
r(h,L(p,g))=0; end end end end clc
disp('邻接矩阵为:'); A
disp('可达矩阵为:'); R
附上截图:
16
上海工程技术大学 《系统工程》课程论文
17