方阵的特征值和特征向量的计箕 A 方阵的特征值和特 征向量的计算 主讲:王开荣 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
方阵的特征值和特征向量的计算 方阵的特征值和特 征向量的计算 主讲：王开荣 11 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
〈方阵的特征值和特征向量的计算 第四章方阵的特征值和特征向量 §1乘幂法 、乘幂法 用来求矩阵按模最大的特征值和相应的特征向量 定理设矩阵A具有n个线性无关的特征向量x,x2,…,xn 其相应的特征值λ,λ1,…,λ满足 1>12213≥…2入 则对任取的一初始非零向量v由 vk=Avk1=…=A-v,k=1,2, 产生的向量序列{vk}满足 aX (2)im( k→》∞ k→o(Vk)m x是1所对应的特征向量,(vm表示v的第m个分量 2 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
方阵的特征值和特征向量的计算 22 第四章 方阵的特征值和特征向量 §1 乘幂法 定理设矩阵A具有n个线性无关的特征向量x1 , x2 , ··· , xn , 其相应的特征值l1 , l2 , ··· , ln满足 一、乘幂法 则对任取的一初始非零向量v0由 产生的向量序列{vk }满足 用来求矩阵按模最大的特征值和相应的特征向量. |l1 |>|l2 |³|l3 |³···³|ln | vk =Avk-1= ··· =Akv0 , k=1, 2, ··· 1 1 1 (1) lim k k k v a x ®¥ l = 1 1 ( ) (2) lim ( ) k m k k m v v l + ®¥ = x1是l1所对应的特征向量, (vk )m表示vk的第m个分量. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
〈方阵的特征值和特征向量的计算 注1收敛速度由比值λ21确定,比值越小收敛速度就 越快,比值接近与1,收敛速度就越慢 注2当矩阵的按模最大特征值是重根时,定理的结论仍 然成立 改进的乘幂法 其中max()表示向量v的绝对值(或模)为最大的分量 maX(] A max(A"-vo) max(v,) max(Avo) lim uk Im k→∞ maX(x lim max(vk)=入1 k→∞ max(x1)k→0 3 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
方阵的特征值和特征向量的计算 33 注1 收敛速度由比值|l2 /l1 |确定, 比值越小收敛速度就 越快, 比值接近与1, 收敛速度就越慢. 注2 当矩阵的按模最大特征值是重根时, 定理的结论仍 然成立. 二、改进的乘幂法 其中max(v) 表示向量 v 的绝对值(或模)为最大的分量 0 0 m 0 ax( ) v u v = 0 0 1 1 0 0 , max( ) max( ) max( ) k k k k k k k k k A v v A v v Au u A v v A v - - = = = = 1 1 lim , max( ) k k x u ®¥ x = 1 1 1 lim , max( ) k k x v x l ®¥ = × 1 lim max( ) k k v l ®¥ = PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
B(阵的持征值和特征向量的计 反幂法 反幂法用来求矩阵A的按模最小特征值及其相应的 特征向量 只要求出A1的按模最大的特征值,也就求出了A的按 模最小的特征值,及其相应的特征向量 §2 Jacobi方法 Jacobi方法用来求实对称矩阵的所有特征值和相应 的特征向量 经典的 Jacobi方法 经典的 Jacobi方法的特点是每次变换将绝对值最大 的非对角元素化为零 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
方阵的特征值和特征向量的计算 44 三、 反幂法 反幂法用来求矩阵A的按模最小特征值及其相应的 特征向量. 只要求出A-1的按模最大的特征值, 也就求出了A的按 模最小的特征值, 及其相应的特征向量. §2 Jacobi方法 Jacobi方法用来求实对称矩阵的所有特征值和相应 的特征向量. 经典的Jacobi方法 经典的Jacobi方法的特点是每次变换将绝对值最大 的非对角元素化为零。 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
阵的特征值和特征向量的计算 计算步骤为 1)首先在4中选择绝对值最大的非对角线元素,设 maxa.≠0 1≠ 选择平面旋转矩阵S1=S(i,i)使T1=S1AS1的非对角线 元素 )=0 (2)再在T1中选择绝对值最大的非对角线元素,设 (1) max|t≠0 12J2 ≠ 选择平面旋转矩阵S2=S(i2)使T2=S2T1S2的非对角线 元素 (2)二f( 0 12/2 此时的t1(2)=51(2又可能变换成了非零元素 (3)重复以上过程,直到满足预定精度为止 5 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
方阵的特征值和特征向量的计算 55 计算步骤为 (1)首先在A中选择绝对值最大的非对角线元素, 设 1 1 max 0 i j ij i j t a ¹ = ¹ 选择平面旋转矩阵S1 =S(i1 , j1 )使T1 =S1 TAS1的非对角线 元素 1 1 1 1 (1) (1) 0 i j j i t t = = (2)再在T1中选择绝对值最大的非对角线元素, 设 2 2 (1) (1) max 0 i j ij i j t t ¹ = ¹ 2 2 2 2 (2) (2) 0 i j j i t t = = 此时的t i1j1 (2)=t j1i1 (2)又可能变换成了非零元素. 选择平面旋转矩阵S2 =S(i2 , j2 )使T2 =S2 TT1S2的非对角线 元素 (3)重复以上过程, 直到满足预定精度为止. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com