pmethod.m

来自「关于矩阵特征值问题的的源代码」· M 代码 · 共 31 行

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)<eps)
        l = m;       %到所需精度，退出，l为主特征值
        s = k;       %s为迭代步数
        return;
    else
        if(k==M)
            disp('迭代步数太多，收敛速度太慢！');
            l = m;
            s = M;
        else
            l = m;
        end
    end
end