矩阵特征值乘幂法.docVIP

矩阵特征值乘幂法 题目 求四阶矩阵的按模最大的特征值及相应的特征向量 精确至第5位有次数字 方法 幂法 程序 pow.M的程序代码 function[m,u,flag] pow A,delta,max1 if nargin 3 max1 100;end if nargin 2 delta 1e-5;end n length A ;u ones n,1 ; flag Fail!;k 0;m1 0; while k max1 v A*u;[vmax,i] max abs v ; m v i ;u v/m; if abs m-m1 delta flag OK!;break; end m1 m;k k+1; end k 结果 format compact A [2 1 0 0;1 2 1 0;0 1 2 1;0 0 1 2] A 2 1 0 0 1 2 1 0 0 1 2 1 0 0 1 2 [m,u,flag] pow A k 12 m 3.6180 u 0.6180 1.0000 1.0000 0.6180 flag OK! eig A ans 0.3820 1.3820 2.6180 3.6180 矩阵的按模最大的特征值是3.6180，相应的特征向量是 ，迭代次数为12次。 拓展 1、应用幂法求矩阵的主特征值的收敛速度主要由比值来决定，但当接近于1时，收敛可能很慢。这时补救方法可以采用加速收敛。 A [2 1 0 0;1 2 1 0;0 1 2 1;0 0 1 2]; [m,u,flag] pow A k 12 m 3.6180 u 0.6180 1.0000 1.0000 0.6180 flag OK! eig A ans 0.3820 1.3820 2.6180 3.6180 2.6180/3.6180 ans 0.7236 A [2 1 0 0;1 2 1 0;0 1 2 1;0 0 1 2.2]; [m,u,flag] pow A k 22 m 3.6491 u 0.5816 0.9591 1.0000 0.6901 flag OK! eig A ans 0.4067 1.4510 2.6932 3.6491 2.6932/3.6491 ans 0.7380 A [2 1 0 0;1 2 1 0;0 1 2 1;0 0 1 2.5]; [m,u,flag] pow A k 26 m 3.7119 u 0.5181 0.8868 1.0000 0.8251 flag OK! eig A ans 0.4353 1.5404 2.8124 3.7119 2.8124/3.7119 ans 0.7577 A [2 1 0 0;1 2 1 0;0 1 2 1;0 0 1 2.9]; [m,u,flag] pow A k 34 m 3.8384 u 0.3944 0.7250 0.9383 1.0000 flag OK! eig A ans 0.4625 1.6336 2.9656 3.8383 2.9656/3.8383 ans 0.7726 由以上结果可以知道，只是微小的改动矩阵中的一个数，在逐渐增大的同时，迭代次数迅速增大。可想而知，当接近1时，收敛速度也会很慢。 2、程序的缺陷 有些矩阵如，返回的特征值不是最大的，这与所选取的初始向量有关系，如果改动初始向量的取值，则能得出正确的结果。此程序中，所选的初始向量为分量值全为1的维列向量，对于有些矩阵就不能得出正确结果。