第八篇 章 特征值问题 矩阵理论课件.pptVIP

第八章 特征值问题;特征值问题是线性代数的研究重点，在理论和应用上都非常重要。 理论上 ，矩阵的特征值就是线性算子的谱。因此可以从泛函分析里找到理论的支撑和生长点。 应用上，常微分方程(ODE)和偏微分方程(PDE)中许多问题都可以转化为矩阵特征值问题。 矩阵特征值问题的算法也是高性能计算机的主要计算任务之一。可大致分为求解稠密、中小型矩阵全部特征值的变换类方法和求解稀疏、大型矩阵部分特征值的投影类方法。;§1、特征值的估计;工程计算中，求解特征值问题;由于矩阵的特征多项式的系数是矩阵元素的连续函数，而多项式的根都是其系数的连续函数，因此矩阵的特征值作为特征多项式的零点都连续地依赖于矩阵的元素。因此矩阵元素的连续变化时，必有对应特征值的连续变化。;二、盖尔（Gerschgorin）定理;定义3 对方阵 ，称;定理4 （Gerschgorin）对方阵 ;（1）的证明：;例 5 矩阵;Wilkson在名著《代数特征值问题》中，先将矩阵变换成Jordan标准型，再用Gerschgorin定理和对角相似变换，对不同特征值结构的特征值问题进行了细致的扰动分析。;例 6 矩阵;三个行Gerschgorin圆分别收缩为：;%ex801.m j=sqrt(-1); A=[20 5 0.8; 4 10 1;1 2 10*j] gersch(A) %调用ATLAST程序库函数gersch hold on v=[1 1 0.5]; D=diag(v); B=inv(D)*A*D %通过对角相似变换隔离Gerschgorin圆 gersch(B) hold off;定义7 对方阵 ，如果;定理8 （Levy--Desplanques） 严格对角占优矩阵是可逆矩阵。 ;根据定理8，严格对角占优矩阵 没有零特征值，而;因此矩阵 严格对角占优，根据定理8;除了Gerschgorin区域外，Brauer发现优美的Cassini椭圆形也可以用于估计特征值。;三、特征值的界;（1）;（1）的证明：;（2）的证明：;在上述证明中，当且仅当 是正规矩阵时，上三角阵 为对角矩阵，即 因此等号都成立。;（1）;（1）的证明：;例 12 矩阵;§2、多项式特征值问题;标准特征值问题 (SEP) ，即为;这里系数矩阵为 矩阵，将矩阵元素即 多项式的次数推广到 次，可得多项式特征值问题(PEP)：;二、广义特征值问题;例 1 (广义特征值问题的直接变换法);例 2 (广义特征值问题的位???求逆法);例 3 (广义特征值问题的Cayley变换法);例 4 (广义特征值问题的Cholesky分解法);设 是广义特征值问题 的特征对， 是标准特征值问题 的特征对。因为矩阵 是Hermite阵，所以有完备的标准正交特征向量系，如果 都是单位特征向量，那么;又由于;这样，对任意 ，可得;例 5 (广义特征值问题的同时合同对角化);证明：由于 是Hermite正定矩阵，所以有;因此;例 6 (基于广义Schur分解的QZ算法****);虽然矩阵束 也存在与Jordan标准型类似的Weierstrass标准型及Kronecher标准型，但也同Jordan标准型一样，在数值计算上存在困难。从数值观点看，更吸引人的是Moler和Stewart (1973)提出的广义Schur分解：;我们知道，在矩阵 的Schur分解 中，上三角 的对角元就是 的特征值。;求解GEP 的QZ算法分三步：;§3、Rayleigh商和广义Rayleigh商;一、 Rayleigh商;定义1 设 是Hermite矩阵，称;单位球面 是闭集，又因为 是 的连续函数，因此根据多元函数的最值定理， 在 上存在最大值和最小值。由于特征值与Rayleigh商的关系，我们自然猜测，这些最值与矩阵 的特征值有什么关系？;可以证明,Rayleigh商是Hermite特征值问题的“最佳逼近”，即对任意 ，有;由于 是Hermite矩阵，因此设 有谱分解;所以;其余特征值是否有类似结论呢？;类似地，假定 ，即;综上，一般地，我们有;遗憾的是，