mifa.m

来自「牛顿迭代法」· M 代码 · 共 21 行

function [m, u, index]=mifa(A, ep, it_max)
% 求矩阵特最大模征值的幂法，其中
%    A      --- 矩阵
%    ep     --- 精度要求
%    it_max --- 最大迭代次数
%    m      --- 最大特征值
%    u      --- 对应最大特征值的特征向量
%   index   --- index=1, 迭代成功，index=0, 迭代失败

if nargin <3  it_max=100; end
if nargin <2  ep=1e-5;    end
n=length(A); u=ones(n,1); index=0; k=0; m1=0;
while k<=it_max
   v=A*u; [vmax,i]=max(abs(v)); m=v(i); u=v/m
   if abs(m-m1)<ep
       index=1; break;
   end
   m1=m 
   k=k+1;
end