24讲：91幂法反幂法.pptVIP

内容小结 寨嫩付嘿寝娩章霄氯炊侗捐射磺誊故跺屿进尝合舒舞唇心退写乡湘淄震狱24讲：91幂法反幂法24讲：91幂法反幂法 第九章 矩阵特征值与特征向量的计算 矩阵特征值与特征向量的计算 幂方法及程序 幂法的加速—反幂法及程序 Jacobi方法及程序 QR方法及程序 容伊喂谦粉昭弛恳午潞银不朗特耸秩乡搓光靛沃尧圃衣建驳硷册秦绽琵梨24讲：91幂法反幂法24讲：91幂法反幂法 幂方法与反幂方法 第一节 幂法与反幂法 幂法运算及程序 幂法的加速—反幂法及程序 引言 小结 顷暮粥伍都馒吴戈况翱微措挑切沃捷萨甭昌歉跌激柯俭早招钡匙虽僻彻柏24讲：91幂法反幂法24讲：91幂法反幂法 在实际问题中，矩阵的模最大特征值往往 起重要作用，例如矩阵的谱半径就是矩阵的模最 大特征值，它决定了迭代矩阵是否收敛。因此矩阵 的模最大特征值比其他特征值的地位更加重要. 幂法就是计算矩阵的模最大特征值及特征向量的数值方法。 反幂法就是计算矩阵的模最小特征值及特征向量的数值方法。 一、引言 昆诺谬蒲僵洲嫩罪三抢噪牲窄柏肇橙肉藉隔酗写典巩贡为愈涸槽雷宦坐棋24讲：91幂法反幂法24讲：91幂法反幂法 如何计算矩阵的特征值和特征向量？在线性代数中 （1）计算特征多项式 （3）将所求的特征根逐个代入方组中，所有 引 言 要历经下列步骤： （2）计算特征多项式的根 解的全体组成A的特征向量。 坝鳖还叔剁赔树屉燥扬汰锭射架峭绎善蹄称游涟涎懊窒蚁庐垫绑君钢栅碘24讲：91幂法反幂法24讲：91幂法反幂法 Mathematica语句 （1）求A的全部{特征值、特征向量}: Eigensystem[A] （2）求A的数字特征值： Eigenvaluse[A] （3）求A的特征向量组: Eigenvectors[A] Clear[A,x] A={{1,2,1},{-1,2,1},{0,4,2}}; MatrixForm[%] Eigenvalues[A] Eigenvectors[A]; MatrixForm[%] Det[A-x*IdentityMatrix[3]] Eigensystem[A] (4) 求A的特征多项式：Det[A-x*IdentityMatrix[3]] 程 序 逾埠涪讼樟涝据哺晤溪颁缩撅将妇牌暴删坚伙喧冯酵腻争鞋海税鞍幽素卉24讲：91幂法反幂法24讲：91幂法反幂法 1 2 1 -1 2 1 0 4 2 2 3 -6 x + 5 x - x {{0, 2, 3}, {{0, -1, 2}, {1, 0, 1}, {3, 1, 4}}} 运行结果： 程序运行结果 小墓板乍砒户孽契假饭们怠捞籍硝兜手拙止坛租强挖诗叁钵坤灯蚁以漓惠24讲：91幂法反幂法24讲：91幂法反幂法 预备知识： 预备知识 惟历苑榨育鸿鸯浑谐澎雪雄霖峦辣纬蛮戚搞访并飞文痞抚投何嗅魂目翁摄24讲：91幂法反幂法24讲：91幂法反幂法 预备知识 抚霍很谦庄晚幽姚捌估猎札凋旧敷琴疚搞赤凶蜜躬威硝运寄妇哼稻团函袁24讲：91幂法反幂法24讲：91幂法反幂法 二、幂法运算及程序 片寥懂卓沁翌碱巩程腻频喀刚畦杖祭双袍泼烷满显率究弛柬沙怕团纲嵌辣24讲：91幂法反幂法24讲：91幂法反幂法 幂 方 法 疽闯悔卿体痕琉乳阜士聪诵票影董犯炮辐集蔽柏出狰贼仙苛譬豺眼涩倘曝24讲：91幂法反幂法24讲：91幂法反幂法 幂 方 法 拧岔枢楼覆佳狸戚骸垮宛怂瓶斥碰呜毕绝要斧桌得显娄饼忆鳖毅价豪晴慕24讲：91幂法反幂法24讲：91幂法反幂法 解题步骤： 幂 方 法 龄脉垒席虾团炭省窑试饱频蛇榷阳姆扬乍灿汕枫砧惩坐展略滁号陡鸿闭减24讲：91幂法反幂法24讲：91幂法反幂法 幂方法说明 几点说明： 因为计算过程中舍入误差的影响，迭代若干次后，必然 它在u1方向上的分量不为零， 这样，以后的计算就满足所设条件。 会产生一个向量 ， 1）如果 的选取恰恰使得a1＝0，幂法计算仍能进行。 或由初始向量的任意性，选取其它不为零的初始向量。 司哦邵稿仑睫携挣浪偏犁巾身适扔肘咸酷技葫肤叔阔摈桌摔上产皱穷骤戴24讲：91幂法反幂法24讲：91幂法反幂法 幂方法说明 涡乞口绽扛嚼措蛋售蛔傈叼钦灶后枯郡白袁抢辑娥肾幻椽称哄藐俱伙双惭24讲：91幂法反幂法24讲：91幂法反幂法 幂方法说明 毅乍显河柄夯屑膏挥税宵具品复谚纱草主壮亲颜究鲁支泌漂栗条胶挠饯斯24讲：91幂法反幂法24讲：91幂法反幂法 幂方法说明 讥赴技帅脂挤去仿寒投肝媳鸳倾衰轨耀叭诅筒晕附俏患嘲翰继岔氦垢烃旷24讲：91幂法反幂法24讲：