第九章nbsp;nbsp;nbsp;nbsp;nbsp;矩阵特征问题的求解.pptVIP

* * 第九章????? 矩阵特征问题的求解 §1 引言 定义1 设矩阵A, B?R n?n，若有可逆阵P，使 则称A与B相似。 定理1 若矩阵A, B?R n?n且相似，则 （1）A与B的特征值完全相同； （2）若X是B的特征向量，则PX便为A的特征向量。 定理2： 设A?R n?n具有完全的特征向量系，即存在n个线性无关 其中?i为A的特征值，P的各列为相应于?i的特征向量。 的特征向量构成Rn的一组基底，则经相似变换可化A为 对角阵，即有可逆阵P，使 定理3 ：A?R n?n，?1, …, ?n为A的特征值，则 （2）A的行列式值等于全体特征值之积，即 （1）A的迹数等于特征值之和，即 定理4 设A?R n?n为对称矩阵，其特征值?1≥?2≥…≥?n，则 （1）对任意A?R n，X≠0， （2） （3） 定理5 （Gerschgorin圆盘定理） 设A?R n?n，则 表示以aii为中心，以 半径为的复平面上的n个圆盘。 （2）如果矩阵A的m个圆盘组成的并集S（连通的）与其余 （1）A的每一个特征值必属于下述某个圆盘之中， n – m个圆盘不连接，则S内恰包含m个A的特征值。 §2 乘幂法与反幂法 9.2.1 乘幂法 定理6 设A ?Rn?n有完全特征向量系，若?1, ?2,…, ?n为A的n个特征值且满足 对任取初始向量X(0) ?Rn，对乘幂公式 确定的迭代序列{Xk}，有下述结论： （1）当 时，对i = 1, 2, …, n 收敛速度取决于 的程度，r 1收敛快，r ? 1收敛慢， 且X(k)（当k充分大时）为相应于?1的特征向量的近似值。 （2）当 时 a）若?1 = ?2，则主特征值?1及相应特征向量的求法同（1）； 收敛速度取决于 的程度。向量 、 c）若 ，则连续迭代两次，计算出X(k+1)，X(k+2)， 分别为主特征值?1、?2相应的特征向量的近似值。 然后对j = 1, 2, …, n 解方程 b）若?1 = -?2，对i = 1, 2, …, n 求出 、 后，由公式 解出主特征值?1、?2。此时收敛速度取决于 的程度。 向量 、 分别为相应于?1，?2 的特征向量的近似值。 规范化乘幂法 令max(x)表示向量X分量中绝对值最大者。即如果有某i0，使 则 max (X) = xi 对任取初始向量X(0)，记 则 一般地，若已知X(k)，称公式 定理7 设A?Rn?n具有完全特征向量系，?1, ?2, …, ?n为A 则对任初始向量X(0)，由规范化的乘幂法公式确定的向量序列 (1) （2）Y(k)为相应于主特征值?1的特征向量近似值 的n个特征值，且满足 Y(k)，X(k)满足 例2 用规范化乘幂法计算矩阵A的主特征值及相应特征向量 取max(x(12)) = 6.00083作为主特征值?1的近似值，与真值?1 = 6相比，有较好的近似程度，相应于?1的特征向量的近似值取为Y(2) = (1, 0.714315, -0.249896)T。 (1, 0.714315, -0.249896) 6.00083 12 (6.00083, 4.286485, -1.499584) (1, 0.714345, -0.239792) 6.00167 11 (6.00167, 4.287265, -1.499172) (1, 0.714405, -0.249586) 6.00335 10 (6.00335, 4.28825, -1.496354) (1, 0.714525, -0.249177) 6.00671 9 (6.00671, 4.291945, -1.496732) (1, 0.714765, -0.248366) 6.013458 8 (6.013458, 4.298211, -1.483536) (1, 0.715247, -0.246768) 6.027024 7 (6.027024, 4.310808, -1.487278) (1, 0.716216, -0.243639) 6.054548 6 (6.054548, 4.336336, -1.475122) (1, 0.718182, -0.237561)