maxtr1.m

本程序使用实现了矩阵的反幂法,希望有所帮助

function [mt,my]=maxtr1(a,ep)
n=length(a);
v0=ones(n,1);
k=1;
v=inv(a)*v0;
while norm(v-v0,'inf')>ep
    k=k+1;
    q=v;
    [c,I]=max(abs(v));
    m=sign(v(I));
    u=v/(m*c);
    v=inv(a)*u;
    v0=q;
end
mt=max(abs(v))
my=u