%% **********计算图的可达矩阵********** % 输入参数： A 图的邻接矩阵 // % 输出： P 图的可达矩阵 // % 功能：计算可达矩阵 // %% function P = dgraf(A) n = size(A,1); P = A; %算法思想：Bn = A + A^2 + A^3+...+A^n for i =2:n P = P+A^i; end P(P~=0)=1;
