mifa.m

来自「计算特征值和特征向量的幂法和反幂法以及计算矩阵的minres 算法的matlab」· M 代码 · 共 18 行

function [p,v,k,toll]=mifa(A)
[m,n]=size(A);
v=[1;ones(m-1,1)];
p=v'*A*v;
k=1;
tol(1)=1;
while tol(k)>1e-6
    v1=A*v;
    v2=norm(v1);
    v=v1/v2;
    p=v'*A*v;
    k=k+1;
    a=A*v-p*v;
    tol(k)=norm(a)/norm(A*v);
end
[m,n]=size(tol);
toll=tol(2:n);
k=k-1;