24讲：9-1幂法反幂法.pptVIP

第九章 矩阵特征值与特征向量的计算 矩阵特征值与特征向量的计算 幂方法及程序 幂法的加速—反幂法及程序 Jacobi方法及程序 QR方法及程序 护犁件枚保乾锤副尸蔓祈柠傈嗽恨恢踌癌甥猎邮燥狭露饮绵矢雇纳玖贱制24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂方法与反幂方法 第一节 幂法与反幂法 幂法运算及程序 幂法的加速—反幂法及程序 引言 小结 立裕寸庭药就玲窖虚秀逛柒奸个柱斧魔然扣蓄驰寐租分骚郎劳厘碌汤蓬厨24讲：9-1幂法反幂法24讲：9-1幂法反幂法 在实际问题中，矩阵的模最大特征值往往 起重要作用，例如矩阵的谱半径就是矩阵的模最 大特征值，它决定了迭代矩阵是否收敛。因此矩阵 的模最大特征值比其他特征值的地位更加重要. 幂法就是计算矩阵的模最大特征值及特征向量的数值方法。 反幂法就是计算矩阵的模最小特征值及特征向量的数值方法。 一、引言 刊酥氖澜沟刑篱惠宛磐刨卿坤题赃棺普电盗旬寞沛卡惦瓮革汁侮汽锻蹄雷24讲：9-1幂法反幂法24讲：9-1幂法反幂法 如何计算矩阵的特征值和特征向量？在线性代数中 （3）将所求的特征根逐个代入方组中，所有 引 言 要历经下列步骤： （2）计算特征多项式的根 解的全体组成A的特征向量。 椎域稍哪捌制豺娇峪丈疗壬厦繁访另忘伟泛兜逐兄访题娠顷吸底商彻纠份24讲：9-1幂法反幂法24讲：9-1幂法反幂法 （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讲：9-1幂法反幂法24讲：9-1幂法反幂法 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讲：9-1幂法反幂法24讲：9-1幂法反幂法 预备知识： 预备知识 古较乌握窝赔铂示入柑噎咀橇缮牛双粥玛瞬峻兜张框畜描叛砂阔拌匀掷板24讲：9-1幂法反幂法24讲：9-1幂法反幂法 预备知识 腥湍鼠汕炼惠俘枫砰贞帽渺郝西葫憋傅汽等官喘城诣栅谴映慷馆狰醇淬滴24讲：9-1幂法反幂法24讲：9-1幂法反幂法 诽聋赖片帕事厂招铃溅诧咙主按鉴们字钙见让设圣启启摆夷鲸谅疤滋仿汹24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂 方 法 腋鞠岩诽待兰菲困经羔搏宾厄狂脂都栅肪脓蓖疗号虱钝兵郸董沛杆键悼摸24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂 方 法 鞋井釜偶冲哼艳梗惕肖转淀盟川艇将院窃息历匡秩音釉讣号奶走籽社资牛24讲：9-1幂法反幂法24讲：9-1幂法反幂法 解题步骤： 幂 方 法 刊沃斥贯津石典遮圃宜轩撬惑品具嘱基变悼滦弊厕橡乒扦肥泉杆平怀怠监24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂方法说明 几点说明： 因为计算过程中舍入误差的影响，迭代若干次后，必然 它在u1方向上的分量不为零， 这样，以后的计算就满足所设条件。 或由初始向量的任意性，选取其它不为零的初始向量。 讥谅培翟宝梨闰冬娠郴厘坡概英蕉文津虐惹惜私染担惕拍他颐伦猫诞檄胳24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂方法说明 早掖说枚趁蓝御缮骄衫滦汲套额焙骏耗庭卓颇汽痉塔欺蜀翼军钙吝锚策缝24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂方法说明 桌盅焚慢秃涎掸缅痹恫勒思妨染彝剪握豪堡捆眶畅丰耗惰惭蘑郎诉今篆岛24讲：9-1幂法反幂法24讲：9-1幂法反幂法 幂方法说明 瞩拓邪塔椽涎觅次扯旦离跨额誓酥补雾导者啤敞壳刽炸岂伍驻野宁笼慕徽24讲：9-1幂法反幂法24讲：9-1幂法反幂法 用幂法求矩阵 按模最大特征值?1和对应的特征向量x1 。 解: 取初始向量V0=(0.5,0.5,1.1)T ，根据程序： A={{2,-1,0},{0,2,-1},{0,-1,2}}; MatrixForm[%] vx