a=[0 0.75 0.6667 0 0 0 0 0 0 0 0 0
0.25 0 0.3333 0 0 0 0 0 0 0 0 0
0.75 0.25 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0.6667 0.3333 0 0 0 0 0 0
0 0 0 0.75 0 0.6667 0 0 0 0 0 0
0 0 0 0.25 0.3333 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0.75 0.75 0 0 0
0 0 0 0 0 0 0.25 0 0.25 0 0 0
0 0 0 0 0 0 0.75 0.25 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0.5 0.25
0 0 0 0 0 0 0 0 0 0.75 0 0.75
0 0 0 0 0 0 0 0 0 0.25 0.5 0];
eps=10^(-10);%给出精度,认为迭代生成的两个c矩阵的差的2范数小于eps,则收敛%
Eps=1;%任意大于eps的值,使循环开始%
k=1;
acell{k}=a;
while(Eps>eps)
k=k+1;
acell{k}=a^k;
Eps=norm(acell{k}-acell{k-1},2);
end
disp(['迭代次数:',num2str(k-1)]);
disp('最终结果,收敛的矩阵c为:');
acell{end}