matlab中可达矩阵的计算

邻接矩阵A 单位矩阵I

A+I的表示

可达矩阵计算程序代码为：n=size(S,1);

p=S;

for i=2:n

p=p+S^i;

end

p(p~=0)=1;

P即为所求的可达矩阵