幂法求矩阵的主特征值及主特征向量
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
幂法求矩阵的主特征值及主特征向量
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) l = m; %到所需精度，退出，l为主特征值 s = k; %s为迭代步数 return; else if(k==M) disp('迭代步数太多，收敛速度太慢！'); l = m; s = M; else l = m; end end end