隱式QR法求實(shí)矩陣的全部特征值matlab實(shí)現(xiàn)

隱式QR法求實(shí)矩陣的全部特征值matlab實(shí)現(xiàn)

要求：用matlab編寫(xiě)通用子程序，利用隱式QR法求實(shí)矩陣的全部特征值和特征向量。

思想：隱式QR法實(shí)質(zhì)上就是將一個(gè)矩陣 Schur化，之后求解特征值就比較方便。而隱式QR法還需要用到household變換，以及上hessenberg變換。最后使用QR迭代，達(dá)到Schur化的結(jié)果。

步驟：

1.將矩陣A上hessenberg化(算法6.4.1)，送而得到一個(gè)上hessenberg形矩陣H；

2.可約性判定，也就是判斷次對(duì)角線元素是否非零，如果次對(duì)角線元素非零，則不可約。

3.Schur化，也就是通過(guò)QR迭代，將矩陣H變化成為某些次對(duì)角線元素變成0，同時(shí)還要滿足，這些元素之間間隔最大為1，那么，所得到的最重的矩陣H就是一個(gè)Schur形矩陣。

4.假如兩個(gè)等于0的次對(duì)角線元素間隔為0，那么該元素的上面一個(gè)元素，也就是H的對(duì)角線上的元素，即為其中一個(gè)特征值；假如兩個(gè)等于0的次對(duì)角線元素間隔為1，那么在這兩個(gè)元素之間就形成了一個(gè)2*2的矩陣，可以求解一個(gè)一元二次方程來(lái)得到兩個(gè)共軛的特征值。

實(shí)驗(yàn)代碼：詳見(jiàn)附錄2

實(shí)驗(yàn)結(jié)果：(代碼相見(jiàn)附錄2)

(i)設(shè)矩陣A如下：

9.1 4.2 3.2 6.13.02.64.0 5.34.71.6 1.79.4x 4.93.56.2

求x=0.9, 1.0, 1.1時(shí)的特征值和特征向量。

總結(jié)

以上是生活随笔為你收集整理的豪斯荷尔德变换 matlab,隐式QR法求实矩阵的全部特征值matlab实现的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。