powermethod.m

来自「Numerical analysis （Gauss、Bisection、jaco」· M 代码 · 共 19 行

function [s,u,c]=powermethod(a,u,N)
if nargin<3
    return
end
i=0;
while i<N
    u=a*u;
    j=1;
    c=u(j);
    for j=1:length(u)-1
        if abs(u(j))<abs(u(j+1))
            c=u(j+1);
        end
    end
    u=u/c;
    s=a*u;
    i=i+1;
end
return