poweriteration.m

来自「实现兰乔斯迭代」· M 代码 · 共 27 行

function [V]=poweriteration(A,kMax,err)
%------200612009佟跃--------
%此程序实现幂迭代算法,输出A的最大特征值的特征向量
tic;
switch nargin
    case 4
        err = eps;
    case 3
        kMax = 10e6;
        err = eps;
    case 2  
      kMax = 10e6;
        err = eps;
    case 1
        KMax=10e6;
        err=eps;
    otherwise
end
[m,n]=size(A);
v0=rand(m,1);
v0=v0/norm(v0);
for k=1:kMax
    w=A*v0;
    v0=w/norm(w);
    end
V=v0;
t=toc;