matlab上机试题(8)

if abs(x1-x)<0.0001 break; elseif j>1000

disp('迭代不收敛或者迭代次数太少'); end end

631 3

A 32110。

5．用逆幂迭代法求 最接近于11的特征值和特征向量，准确到 111

解答：以下所得结果是最小特征值对应的参数结果

function [m,u,index,k]=pow_inv(A,ep,it_max) %求矩阵最小特征值的反幂法，其中A为矩阵； %ep为精度要求，缺省为1e-5；

%it_max为最大迭代次数，缺省为100；m为绝对值最大的特征值； %index，当index=1时，迭代成功，当index=0时，迭代失败 if nargin<3 it_max=1000;end if nargin<2 ep=1e-5; end

n=length(A);index=0;k=0;m1=0;m0=0;

%修改移位参数，原点移位法加速收敛，为0时，即为反幂法? I=eye(n); T=A-m0*I; invT=inv(T); u=ones(n,1)

while k<=it_max v=invT*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 m=1/m; m=m+m0;

以下所得结果是最大特征值对应的参数结果 function [m,u,index,k]=pow(A,ep,it_max)

% 求矩阵最大特征值的幂法，其中A为矩阵； % ep为精度要求，缺省为1e-5；

% it_max为最大迭代次数，缺省为100；