- 1、本文档共43页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第八章 矩阵特征值问题计算 第八章 代数特征值问题 两种特殊情况 幂法小结 二、幂法的加速 三、反幂法 反幂法的一个应用 因为幂法的收敛速度是线性的,而且依赖于比值 ,当比值接近于1时,幂法收敛很慢。幂法 加速有多种,介绍两种。 反幂法是计算矩阵按模最小的特征值及特征向量的方法,也是修正特征值、求相应特征向量的最有效的方法。 * 数值分析 数值分析 * 工程实践中有多种振动问题,如桥梁 或建筑物的振动,机械机件、飞机机翼的振动,及 一些稳定性分析和相关分析可转 化为求矩阵特征值与特征向量的问题。 但高次多项式求根精度低 , 一般不作为求解方法. 目前的方法是针对矩阵不同的特点给出不同的有效方法. 第一节 特征值的估计和数值稳定性 第二节 幂法和反幂法 第三节 求矩阵全部特征值的QR方法 定理1 设 为 的特征值, , 则 (1) 为 的特征值( 为常数 ); (2) 为 的特征值,即 (3) 为 的特征值; 第一节 特征值的估计和数值稳定性 定理2 (1)设 可对角化,即存在非奇异矩阵 使 的充要条件是 具有 个线性无关的特征向量. (2) 如果 有 个 不同的特征值 则对应的特征向量 线性无关. 而 的列 向量 为 的对应于 的特征向量. 定理3 设 为对称矩阵,则: (1) 的特征值均为实数; (2) 有 个线性无关的特征向量; (3)存在一个正交矩阵 使 且 为 特征值, 记 称为矩阵 的瑞利(Rayleigh)商. 定理4 设 为对称矩阵(其特征值次序记为 则 证明 只证 1. 由于 为实对称矩阵,可将 对应的特征 向量 正交规范化,则有 (1.3) 设 为 中任一向量,则有展开式 于是 从而1成立. 结论1说明瑞利商必位于 和 之间. 8.1.2 特征值估计与扰动 称复平面上以 为圆心,以 为半径的所有圆盘为 的格什戈林(Gerschgorin)圆盘. 定义1 设 .令 (1) (2) 集合 . 定理5 (格什戈林圆盘定理) (1) 设 , 则 的每一个特征值必属于下述某个圆盘之中 (1.4) 或者说, 的特征值都在复平面上 个圆盘的并集中. (2) 如果 有 个圆盘组成一个连通的并集 , 且 与余下 个圆盘是分离的,则 内恰包含 的 个 特征值. 特别地,如果 的一个圆盘 是与其他圆盘分离的 (即孤立圆盘),则 中精确地包含 的一个特征值. 证明 只就(1)给出证明. 设 为 的特征值,即 记 考虑 的第 个方程, 即 或 于是 即 这说明, 的每一个特征值必位于 的一个圆盘中, 并且相应的特征值 一定位于第 个圆盘中. 其中 是对应特征向量 绝对值最大的分量的下标. 利用相似矩阵性质,有时可以获得 的特征值进一步 的估计,即适当选取非奇异对角阵 并做相似变换 . 适当选取 可使某些圆盘半径及连通性发生变化. 例1 估计矩阵 特征值的范围. 解 的3个圆盘为 由定理5,可知 的3个特征值位于3个圆盘的并集中, 由于 是孤立圆盘,所以 内恰好包含 的一个特征值 ( 为实特征值),即 的其他两个特征值 包含在 的并集中. 现选取对角阵 做相似变换 的3个圆盘为 这样,3个圆盘都成为了孤立圆盘,每一个圆盘都包 含 的一个特征值
文档评论(0)