幂法与反幂法.doc

幂法与反幂法 1 功能 幂法是一种计算矩阵主特征值(矩阵按模最大的特征值)及对应特征向量的迭代方法, 特别是用于大型稀疏矩阵。 反幂法用来计算矩阵按模最小的特征值及其特征向量，也可用来计算对应与一个给定近似特征值的特征向量。 2算法描述 2.1 幂法 （1）取初始向量u（例如取u=(1,1,…1)）,置精度要求，置k=1. （2）计算 v=Au,m=max(v), u= v/ m （3）若| m= m|，则停止计算（m作为绝对值最大特征值，u作为相应的特征向量）否则置k=k+1，转（2） 2.2 反幂法 （1）取初始向量u（例如取u=(1,1,…1)）,置精度要求，置k=1. （2）对A作LU分解，即A=LU （3）解线性方程组 Ly=u,Uv=y （4）计算 m=max(v), u= v/ m 若|m=m|，则停止计算（1/m作为绝对值最小特征值，u作为相应的特征向量）;否则置k=k+1，转（3）. 3 Matlab程序的实现 3.1 幂法 function[m,u]=pow(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; while k=N 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 输入： A=[7 3 -2;3 4 -1;-2 -1 3]; [m,u]=pow(A,1e-6) Enter 结果： m = 9.6056 u =1.0000 0.6056 -0.3944 4.2 反幂法 function[m ,u]=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); while k=N v=invA*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 ; 输入： A=[6 2 1;2 3 1;1 1 3]; [m,u]=pow_inv(A,1e-6) Enter 结果： m = 1.7733 u = -0.3473 1.0000 -0.5321 4 Matlab程序调用示例 随机产生一对称矩阵，对不同的原点位移和初值(至少取3个)分别使用幂法求计算矩阵的主特征值及主特征向量，用反幂法求计算矩阵的按模最小特征值及特征向量。 要求1）比较不同的原点位移和初值说明收敛性 2）给出迭代结果，生成DOC文件。 3）程序清单，生成M文件。wo 解答： A=rand(5) %随机产生5*5矩阵 求随机矩阵 A = 0.7094 0.1626 0.5853 0.6991 0.1493 0.7547 0.1190 0.2238 0.8909 0.2575 0.2760 0.4984 0.7513 0.9593 0.8407 0.6797 0.9597 0.2551 0.5472 0.2543 0.6551 0.3404 0.5060 0.1386 0.8143 B=A+A %A矩阵和A的转置相加，得到随机对称矩阵B B = 1.4187 0.9173 0.8613 1.3788 0.8044 0.9173 0.2380 0.7222 1.8506 0.5979 0.8613 0.7222 1.5025 1.2144 1.3467 1.3788 1.8506 1.2144 1.0944 0.3929 0.8044 0.5979 1.3467 0.3929 1.6286 B= 5 C语