第八章 特征值问题的计算方法课件.pptVIP

?经典Jacobi方法的迭代算法 给定矩阵 选取最佳旋转平面 ： For k=1,2,3,… 计算 计算 直到 需比较 个元素 毅湘聚手脉佰呢暗避净奔鹰恿轴施感求拧崩源袜信醛采倪集瞅照戌癣庶悦第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 习惯上称 次Jacobi迭代为一次“扫描” ?循环Jacobi方法 每一次Jacobi迭代不是去选择最佳旋转平面， 而是直接按照某种预先指定的顺序来“扫描” 自然顺序： 按照自然顺序的循环Jacobi方法是渐进平方收敛的 绦避病封摊坡眯另灰洱桅楞拯芳啮喷嗡霜踩狙蜜柳甲尘迹坍续禽晨拴恤缝第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 §4 QR 方 法 ?基本思想 利用正交相似变换将一个给定的矩阵逐步约化 为上三角矩阵或拟上三角矩阵的一种迭代方法 ?QR方法的迭代格式 设 令 对矩阵 进行QR分解 再对矩阵 进行QR分解 一、QR基本迭代方法 QR方法是目前计算矩阵全部特征值的最有效 的方法之一；具有收敛快、算法稳定等特点。 褥牵辛误缔寞任棵芋姿菲诊窑板图素凛佳抓名吼庄潦壁湛幅凿却诣堕饯雕第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 一般地有： 矩阵序列 中每一个矩阵都与原矩阵 相似 QR方法的迭代算法： For m=1,2,3,… 直到 近似 为上三角阵 由迭代格式同时还得到： 记 代入 嘱帖扼筐浚咐筷廊捌衅崇虽鲍实渝榷坍敢抬隶樟企没楔炬矽规纬浮咋赴洲第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 等式两端同时右乘 记 即 其中 是 的第一列， 是 的相应元素 可以看作是对矩阵 用 为 初始向量的幂法所得到的向量。 QR方法与幂法的关系： 例晒企跳骤络讽歌蛔耶粥包俊笺蹬嘘况卵稍奈憋资虎捌媳屋碱隋综又啤现第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 ?QR方法的收敛性 设 的特征值满足 ，且 的 则由QR迭代算法产生的矩阵 的对角线以 第i行是 对应于 的左特征向量；若 有 分解， 下的元素趋于0，同时对角元素 趋于 （QR方法的收敛性质） 证明： 令 则有 其中 锗竞解存趁船擞寺遵渣杀轮告琵戍裕踩刨淑博并弃套彼忠沁踌砧鸽操冬钞第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 且 当m充分大时， 存在唯一QR分解： 且 当m充分大时 （QR分解） 记 的QR分解为： 搂埠皑谭秧在妊演台爱痉辗描狂怂赠丈幂溺留欣泻泳点涨萤殆系殿戏损牌第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 为保证上述QR分解中上三角矩阵的对角元为正，令 由QR分解唯一性知： 代入 趋于上三角阵 狱瞬背碟涩剑吝景揉曰骑瑶珊尼殖符先亏甥随壶慢如贤论散展薛滇盆怂兹第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 实际应用中遇到的多数特征值问题都是关于实矩阵 的，所以自然希望设计只涉及实数运算的QR迭代。 ?实QR迭代格式 设 二、实Schur标准形 For k=1,2,3,… 为正交矩阵 为上三角阵 （实Schur分解） 设 ，则存在正交矩阵 ，满足： 其中 为实数或具有一对复共轭 特征值的2阶方阵 章埂簧呸鞍败筐旷氓咐票币蔷晰厉岔铡哺年筏滁伴钩私驾狠罐溢戒利首龚第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 特征值为 ，其中 为虚单位 见文献[13] 矩阵 称为 的Schur标准形 定理8.4.2说明：只要求得矩阵的Schur标准形， 就很容易求得矩阵的全部特征值。 缺陷： 很难得到 闯四檄库野官豆黄商赎汕苹稻噪慷疽往圣媒谋棺熄篮千痢汲账菱扇固纲架第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 称下述形状的矩阵为上Hessenberg矩阵 三、上Hessenberg化 箱叶榷岳胀脓认殷婉久员腔您困肠州函汇缆麓除核月钞坡土国哇醋剔酗呼第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 设 ， 寻求非奇异矩阵 ，使得 为上 Hessenberg 矩阵，然后再对 进行 迭代。 ?基本思想和约化过程： 记矩阵 下面采用Householde变换寻找 Step1 选取Householde变换 ，使得 其中 令 聂刑腐旅刊忽援警踊匿价澜吝秋身抗兔眨炮涨用蜘隧厅钢我接启眯痹栖椭第八章 特征值问题的计算方法课件第八章 特征值问题的计算方法课件 *