四章矩阵特征问题求解简.ppt

2．雅克比方法 引例：考虑矩阵 当适当选取 ? 让 时，有 其中 ?1=a11cos2?? a12sin2? + a22sin2?， ?2=a11sin2?+ a11sin2? + a22cos2? 定义：称矩阵 为平面旋转矩阵 第i行 第i列 第j列 第j行 平面旋转矩阵P 的性质： 1) P是正交矩阵，即 PPT=E; 2) A是实对称矩阵，则 B=PTAP 仍是实对称矩阵，且A与B具有相同的特征值； 3) 若 B=PTAP ，则||A ||F= || B ||F，其中 即经旋转变换后，对角线上元素的平方和将增加. 4)取 时，有 雅可比 (Jacobi) 旋转法的方法步 1.选取矩阵A的非对角线元素中绝对值最大的元素 aij= aji (ij ) ； 2.确定平面旋转矩阵P1=Pij(?),其中 ?与 aij 同号。作变换A1= P1TAP1=(aij(1)) 3.重复2的工作：作变换Ak= PkTAPk=(aij(k)), Pk=Pij(?), 则有如下关系： 将? 限制在下列范围内 如果 这样得到的矩阵序列{Ak}是收敛的， 且Ak ??（k ? ?）aii(k) ??i （ k ? ?）; 在迭代过程中，若max| aij(k) | ? ，结束运算。 2.取 P=P1 P2 …Pk中的列向量即为?i所对应的特征向量的近似值 则1.取?i ? aii(k) ， 例：用Jacobi方法计算矩阵A的特征值和相应的特征向量 特征值的精确解 ? 1=20.96915893998615 , ?2= -0.9339707868005558 , ?3＝ 0 .465911846814403 特征向量: 直接从三角函数关系式计算sin? 和cos?，记 则 当 时，有下面三角恒等式： * 计算方法与数值计算 University of Shanghai for Science and Technology College of Science 上 海 理 工 大 学 理 学 院 §1 引言 矩阵特征问题的求解 物理、力学及工程技术领域中的许多问题在数学上都可归结为求方阵特征值和特征向量的问题，如振动问题：桥梁的振动、机械振动、电磁振荡、地震引起的建筑物的振动等等；物理学中的临界点、临界值的确定；微分系统中的稳定性研究。这些常见的问题都与方阵的特征值和特征向量有关。 求方阵的特征值和特征向量是数学、物理学、力学、电磁电工学以及工程技术所面临的一个共同问题。 定义：矩阵 A?R n?n，若有数?和非零向量x ?R n 使得A x＝ ? x，则 称? 为A的特征值， x称为?对应的特征向量。 在线性代数中的求法：解特征多项式| ?E- A |=0. 利用线性代数中的上述方法计算特征值和特征向量是十分困难的; 我们的目的是寻求一种适合计算机运行的近似解法,且简单、可行、有效。 复习相关理论 定义 设矩阵A, B?R n?n，若有可逆阵P，使 定理1 若矩阵A, B?R n?n且相似，则 B =P -1 AP 则称A与B相似。 （1）A与B的特征值完全相同； （2）若x是B的特征向量，则Px便为A的特征向量。 定理2： 设A?R n?n具有完全的特征向量系，即存在n个线性无关 其中?i为A的特征值，P的各列为相应于?i的特征向量。 的特征向量构成Rn的一组基底，则经相似变换可化A为对角阵。 即 矩阵A与对角阵?相似的充要条件是A 有n 个线性无关的特征向量 即有可逆阵P，使 定理3 ：A?R n?n，?1, …, ?n为A的特征值，则 （2）A的行列式值等于全体特征值之积，即 （1）A的迹数等于特征值之和，即 定义： 设A?R n?n为对称矩阵，对任意x?R n，x≠0，则称比值 矩阵 A关于向量 x 的瑞雷(Rayleigh)商 瑞雷(L.Rayleigh)——英国数学家，1842－1919 定理4 设A?R n?n为对称矩阵，其特征值?1≥?2≥…≥?n，则 1） 对任意 x?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的特征值。 1 乘幂法 定理6 设A ?Rn?n有n个线性无关的特征向量，若?1, ?2, …, ?n为A的n个特征值且满足 对任取初始向量x(0) ?Rn， (x(