MATLAB-矩阵计算-反幂法程序

/反幂法程序，函数名：pow_inv.m/

function[m,u,index]=pow_inv(A,ep,N)

%A为矩阵；ep为精度要求；N为最大迭代次数；m为绝对值最小的特征值；u为对应最小特征值的特征向量。

N=100;

ep=1e-6;

n=length(A);u=ones(n,1);

index=0;k=0;m1=0;

invA=inv(A);

whilek<=N

v=invAu;[vmax,i]=max(abs(v));

m=v(i);u=v/m;

ifabs(m-m1)
index=1;break;

end

m1m;k=k+1;

end

m=1/m;