pmethod.m
function [l,v,s]=pmethod(A,x0,eps)
if nargin==2
eps = 1.0e-6;
end
v = x0; %v为主特征向量
M = 5000; %迭代步数限制
m = 0;
l = 0;
for(k=1:M)
y = A*v;
m = max(y); %m为按模最大的分量
v = y/m;
if(abs(m - l)
s = k; %s为迭代步数
return;
else
if(k==M)
disp('迭代步数太多，收敛速度太慢！');
l = m;
s = M;
else
l = m;
end
end
end
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024