幂法-反幂法求解矩阵最大最小特征值及其对应特征向量.pdfVIP

数值计算解矩阵的按模最大最小特征值及对应的特征向量 一.幂法 1. 幂法简介： 当矩阵 A 满足一定条件时，在工程中可用幂法计算其主特征值 (按模最大 ) 及其特征向量。矩阵 A 需要满足的条件为： (1) | 1 | | 2 | ... | n | 0 , i 为A 的特征值 (2) 存在 n 个线性无关的特征向量，设为 x , x ,..., x 1 2 n 1.1 计算过程 ： n 对任意向量 x ( 0) ，有x (0 ) u , i i i 不全为 0，则有 i 1 x(k 1) Ax(k) ... Ak 1x( 0 ) n n k 1 k 1 A αu αλ u i i i i i i 1 i 1 k 1 2 k 1 n k 1 λ u ( ) a u ( ) a u 1 1 1 2 2 n n 1 1 k 1 u 1 1 1 可见，当 | 2 |越小时，收敛越快；且当 k 充分大时，有 1 (k 1) k 1 (k 1) x 1 1u1 x (k 1) (k ) k (k ) 1 ，对应的特征向量即是 x 。 x u x 1 1 1 2 算法实现 (1).输入矩阵 A ，初始向量 x，误差限 ,最大迭代次数 N (k ) ( k ) x (2).k 1, 0; y ( k ) max(abs(x ) (3).计算 x Ay, max(x); 若 输出 否则 转 (4). | | , , y, , (5) (5).若 k N , 置k k 1, ,转3 , 否则输出失败信息 ,停机 . 3 matlab 程序代码 1 / 17 function [t,y]=lpowerA,x0,eps,N) % t 为所求特征值， y