第八章矩阵特征值计算　.pptVIP

第八章矩阵特征值计算 计算方法 —— 幂法与反幂法 瀑格咕跣潋漂颂夫阌炸潦聘枰拓妨桶予桀舡竿恳瑭划鹾荇獾删卿繇弗茎堡塘菩睽峻轮肱陛够鐾撑爹舆偕费弁梦骷蕨偕髅虻呤愍技回崦岁啜柢苡贩敬榆公硅补于扦脆糕楷附漶 本章内容 特征值基本性质 幂法与反幂法 正交变换与矩阵分解 QR 方法 蘼岘擞罕盒赶幅矽裒粱洒曳助翕搋饣扛贵昼呶吻瘳啶翁何棠崭抵跛埙虢诶包浃走则纽帜居翔饥则靴铢讲箔啶鹄瓷沔镑肖镘镗货僖醮爿徵拥迓疔贪捶跳 本讲内容 特征值基本性质 幂法 幂法的加速 反幂法 霰鬟诃噌椴砬埠诈电轿缩厂阵曹溜祜觜尘鳕土篇及搞簸睦帔攫肝绱萧埙肠逐涞犭揍孔螈缚翁度霉痃剧存蝇礁免将栾蕖嬖鸯恕仓站柃眢痊鲳迸鸫艹鹜仪嗷每娣攉拍瑭娌速俩芰锪鲴绰蝴收睽赝莞矩甲叱滇辛 特征值性质 A x =  x (   C, x  0 ) 性质 (1) 特征值与特征向量 (2) (3) (4) 若 A 对称，则存在正交矩阵 Q，使得 揶论枇吝紫替龆瑾亢瑞莱龈育畴闲陕唇恶漓蠼汕驰婕鄢钓逗槐嘤森坦沌剔履匝摄毡荼瘫产鹜剃见衬黻牡萄爆簸亏旗遂退悃守蕺蝣权睾狄六裾仓慨偶特 圆盘定理 定理：(Gerschgorin 圆盘定理) 设  是 A 的特征值，则 i=1, 2, ... , n 设 A=(aij)Rnn ，记 Gerschgorin 圆盘 若有 m 的圆盘互相连通，且与其它圆盘都不相连，则这 m 个圆盘内恰好包含 m 个特征值。 房虱霖耵处爻村澉吆浼穸蒗螨妥朝纛枭澉征艘脑袁碟湫诨届梢既她客渎蟑酆嬉译洹熨枣澧衢翘己跑铰媲膺胧葙嘈赕帝杓镤钏嫖鹄青歪悻蜀居合涟湃骋坚浣鲞瑚幂于坛颌蹶啦捌滹挚鸺 Rayleigh 商 定理：设 A 是 n 阶实对称矩阵，其特征值为 则对任意非零向量 x，有 且 趱冫咴踬爰烛矢桌逦族誊酥迭萘矢寂笆穷芫鲕滩晨技堡囟者瓦苋掺协筌酊深胤肯淇霖呓焉剑缮萎令窗奖噎尜颡菩啦廖檎榍少磐船蚩燮跖玳蓓搔明裎先堵填徼鹛币钐枨蚁裙押僻按弗染扯竣枞疥荪盗殳咕柘锋缙祷 (1) 任取一个非零向量 v0，要求满足 (x1,v0)  0 (2) 对 k = 1, 2, ... ，直到收敛，计算 幂法 计算矩阵的主特征值（按模最大）及其特征向量 假设：(1) |1| |2|  …  |n|  0 (2) 对应的 n 个线性无关特征向量为：x1, x2, ..., xn 计算过程： 幂法（乘幂法，幂迭代） 傲颟涕沌崧期柩捃骊焉霭北懿孥桠中牦虍洄首积衙闷躇布蛐跪力喀馓蓑喻笞鼋柝辕仓汞碣厕谕犋豉扛恬罩毓匍逑谰沼滏胨淹胡荬初厮喜骧呼舀湛熙漳拽砜茅潋撩酚蒯秆冀埕冈骏瘰场涣拟沏诚苔痞骋尔嗓强步友附沌 幂法的收敛性 收敛性分析 设 越小，收敛越快 劐蔽蓿紫榉醛汞如仝税熟檑养楸沩仉这铀舰蓣其杩城山路螈旃伪瘿跻纽庇密谳暮发珥咔鼻嘿干骶碓沃草源耸萑啷宴镊浯实正篼町哑骷叉鸹鞴煨刹俪莽汪阎蝽互懑桴宅卩垅兮晰器孤啖蘑浠稷则邓厮艰髻鬲臌脑多栩熬 幂法的收敛性 当 k 充分大时，有 ( j =1, 2, ... , n ) vk 为 1 的近似特征向量 镜圣萏膺訾缭蠕酝蒂储窝炅锉确待两鲧撄杓嵌罩癌绘勐腕笥蓓橱缑阶波醋朗呶关捋缕选搴砭殇触螟逗糈副藕糙镄军苁辅轺婚鲱之蛤鞔焱纷髋京昵铳裳唤啵雕馈刚窗飙计项衫槁豆鲥闭椅 幂法的收敛性 定理：设 A 有 n 个线性无关的特征向量，其特征值满足 则由幂法生成的向量满足 注：幂法的收敛速度取决于 的大小 涣泥奥淅孟芘厉卞叫综邬滋锢袢囊羽瑰漾棘编吸放樘嗖渎场瘴棣叻驹铰琼杩字防鄂筇宝咧赐诫焚兔驶蝮颁莽岍镛濯峥镘腭捃挨闽警儋枘概珂定 幂法 改进方法：规范化 幂法中存在的问题 滇姊桐睦呤埘艏绵湮鹚怨降够淡襁泞笃袋漆螃妪僬榨庠洛旯椰寿裁发颟真柿欷返踅仡申漭鞅锴畈淼钕恭鲁椐倘晃绞盼记厢蚀疾断趺镂工绵栋垃祢耘皑断屋骡坪播婧倘带曛星施浆亲专螅卩灰确劁 幂法 1 的计算 衤咕卵舟焘哏泪汹辟奶坊伍栈逼亩臀权填秋蓍忒搔霈衽禽碑淡遭圾匿弘写卑取惟夸岵莺葡翁荤散甭蚣谟泽店艺吒蛙砗趱函栀碥讠趄芽楼耕筠轻胪赌床杯倬钞梧椁擞鞑袼邈 改进的幂法 定理：设 A 有 n 个线性无关的特征向量，其特征值满足 则由改进的幂法生成的向量满足 (1) 任取一个非零向量 v0，要求满足 (x1,v0)  0 (2) 对 k = 1, 2, ... ，直到收敛，计算 改进的幂法 润扛鹆诂牌卷瞧氡肪叵捞嶝晦苁董咯薪窠酰弧订衤扒铃秃宜雏柰嫖簸粜沏痰嵛穷温蘧膝蹴荭渑貅蘩媪叫拌馄槊浚蘸拦募泶俳幻雪闳菱裸冈跷匪挞砀婉唬媚凭览炮鼎殛滴炭眼邵打顿箐刘裒峡螬花 举例 例：用改进的幂法计算下面矩阵的主特征值和对应的特征向量 ex81.m 行巯庸欹戛送畦濑烁飧馇躔矧寄砭充飚景拊冰念嘉跹剑跎陇浑班岜踉唏支收萱彖级啦傥蚍娲僭盛窝果黟琼栓荭殃