幂法,反幂法求解矩阵最大最小特征值和对应的特征向量.docxVIP

WORD 格式整理 数值计算解矩阵的按模最大最小特征值及对应的特征向量 一. 幂法 幂法简介： 当矩阵 A 满足一定条件时，在工程中可用幂法计算其主特征值 ( 按模最大 ) 及其特征向量。矩阵 A 需要满足的条件为： (1) | 1 | | 2 | ... | n | 0 , i 为A的特征值 (2) 存在 n 个线性无关的特征向量，设为 x1, x2 ,..., xn 1.1 计算过程 ： 对任意向量 x(0 )，有 x( 0) n i ui , i 不全为 0，则有 i 1 x(k 1) Ax(k) ... Ak 1 x( 0 ) n A k 1 n k 1 αu αλ u i i 1 i i i 1 i i k 1 u ( 2 ) k 1 a u ( n ) k 1 a u λ n 1 1 1 2 2 n 1 1 k 1 u 1 1 1 可见，当 | 2 | 越小时，收敛越快；且当 k 充分大时，有 1 x (k 1 ) k 1 x(k 1 ) 1 1u1 ，对应的特征向量即是 x(k 1 ) x (k ) k 1u1 x (k )1 。 1 算法实现 专业知识分享 WORD 格式整理 输入矩阵 ，初始向量 ，误差限 , 最大迭代次数 N (1). A x (2).k 1, 0; y( k) x( k ) max(abs( x( k ) ) 计算 x Ay, max(x); (3). 若 | | , 输出 , y, 否则 转 (5) (4). , 若 k N , 置 k k 1, , 转 3 , 否则输出失败信息 , 停机 . (5). 3 matlab 程序代码 function [t,y]=lpowerA,x0,eps,N) % t 为所求特征值， y是对应特征向量 k=1; z=0; % z 相当于 y=x0./max(abs(x0)); % 规范化初始向量 x=A*y; % 迭代格式 b=max(x); % b 相当于 if abs(z-b)eps % 判断第一次迭代后是否满足要求 t=max(x); return ; end while abs(z-b)eps kN k=k+1; z=b; y=x./max(abs(x)); x=A*y; b=max(x); 专业知识分享 WORD 格式整理 end [m,index]=max(abs(x)); % 这两步保证取出来的按模最大特征值 t=x(index); % 是原值，而非其绝对值。 end 举例验证 选取一个矩阵 A，代入程序，得到结果，并与 eig(A) 的得到结果比较，再计 算 A*y-t*y ，验证 y 是否是对应的特征向量。结果如下： 专业知识分享 WORD 格式整理 结果正确，表明算法和代码正确， 然后利用此程序计算 15 阶 Hilb 矩阵，与 eig(A) 的得到结果比较，再计算 A*y-t*y ，验证 y 是否是对应的特征向量。设 置初始向量为 x0=ones(15,1) ，结果显示如下 专业知识分享 WORD 格式整理 可见，结果正确。得到了 15阶 Hilb 矩阵的按模最大特征值和对应的特征向 量。 二. 反幂法 反幂法简介及其理论 在工程计算中，可以利用反幂法计算矩阵按模最小特征值及其对应特征向 量。其基本理论如下，与幂法基本相同： 由 Ax x x A 1 ，则 A 1 1 -1 的特征值互为倒数， ( x) x x ，可知，A 和 A -1 求 A 按模最小特征值即求 A 的按模最大特征值 , 取倒数即为 A 的按模最小特征值 所 以 算 法 基 本 相 同 ， 区 别 就 是 在 计 算 专业知识分享 WORD 格式整理 x( k 1)时，不是令 x( k 1) Ay(k ) , 而是 x( k 1) A-1 y(k )具体计算时，变换为 Ax( k 1) y( k ) ; 对 A做LU 分解，来计算 x( k 1) 算法实现 (1).输入矩阵 A,初始向量 x,误差限 , 最大迭代次数 N , (2).置 k 1, 0 0, y x , max( abs( x)) (3).作三角分解 A LU (4).解方程组 LUx y( Lz y,Ux z), (5). max( x), (6).若 | 0 | ,输出 1 , y,停机 ,否则转 (7), (7).若 k N , 置 k k 1, 0, y x , 转 (4); max(abs( x)) 否则输出失败信息 ,停机 . 3 matlab 程序代码 function [s,y]=invpower(A,x0,eps,n) % s 为按模最小特征值， y是对应特征 向量 k=1;