clear;clc A=[ 0 0 0 0 0 1
1 0 0 0 0 0
0 1 0 0 0 0
0 1 1 0 0 0
0 0 1 0 0 0
0 0 1 0 1 0]; %有向图邻接矩阵
[n,n] = size(A);
sn = sum(A,2); %按列求矩阵A各行的和
d = 0.85; %阻尼因子
D = zeros(1,n);
for j=1:n
if sn(j)==0
A(j,:)=ones(1,n)/n;
else
A(j,:)=A(j,:)/sn(j);
end
end
% A = (1-d)/n*ones(n)+d*A';
p = ones(n,1)/n;
z = zeros(n,1);
count=0;
while max(abs(p-z))>0.0001
z=p;
% p=A*p;