為什麼矩陣QR分解迭代的主對角會最終收斂到特徵值?

時間 2021-06-08 23:58:14

1樓:sy sy

並不很熟悉數值代數,試著來強答一下。。

給定階方陣,QR演算法的迭代步驟為

將證明如下的收斂定理:

設的特徵值滿足|\lambda_2|>...>|\lambda_n|>0" eeimg="1"/>,則可對角化,。進一步,假設有LU分解。

則時,對角線下的元素趨於0,對角線元素趨於特徵值。

證明需要以下四個引理。

引理1(Cholesky分解)

正定,當且僅當存在唯一的對角元均正的上三角陣,。

引理2(Cholesky分解的連續性)

Cholesky分解誘導了「正定陣對角元正的上三角陣」的連續對映:。

引理3(QR分解)

可逆,當且僅當存在唯一的酉陣,唯一的對角元正的上三角陣,。

引理4(QR分解的連續性)

QR分解誘導了「可逆陣酉陣對角元正的上三角陣」的連續對映:。

引理1和3是標準的矩陣論結果,證明略掉。如果引理2對,則通過的連續性,可以得到的連續性,從而引理4成立。只需證明引理2:

階數,,連續。

假設連續性在階成立。考慮分塊

對比左右兩邊各分塊,可得

由此可見是的連續函式,由歸納假設是的連續函式,證畢。

收斂定理的證明:

一方面,記。反覆使用QR演算法的迭代公式,可得

另一方面,已知,。通過插入可逆的對角陣,可以使對角元全為1。再令的QR分解為,於是

可逆,令其QR分解為,於是

是上三角陣,通過左乘么模可逆對角陣,可以讓其對角元全為正,因此有

對比兩個分解式,由QR分解的唯一性,有

將其帶入分解式,可得

注意到因此,由QR分解的連續性,我們有

將其帶入分解式,注意到是對角線全1的上三角陣,因此時,對角線下的元素趨於0,對角線元素趨於特徵值,證畢。

實對稱矩陣相似對角化為什麼要弄成正交矩陣?

不開心的Lair 一般矩陣相似對角化的時候,只需要求出線性無關的特徵向量,將特徵向量組成乙個矩陣就可以對矩陣進行正交化,而對實對稱矩陣需要進行正交單位化的操作,是因為這裡是用正交矩陣來對實對稱矩陣進行對角化。然而對於實對稱矩陣特有的乙個性質就是不同特徵值對應的特徵向量相互正交。我們知道特徵向量乘乙個...

N 階七對角矩陣的特徵值是什麼?

這雖然是乙個數學問題,但我作為學物理的看到這個問題首先考慮的是他的物理意義.題主想要考慮的是如下形式的矩陣 在固體物理裡,這個矩陣對應著乙個有著次次次近臨躍遷的緊束縛模型.根據固體物理的知識,立刻知道這個系統的能量本徵值是 簡單起見,我們取無關輕重的 a 0.其中 k 是動量.剩下唯一的問題就是在這...

主成分分析PCA演算法 為什麼要對資料矩陣進行均值化?

我叫平沢唯 只有居中,才能使均方誤差最小 以下來自位維基百科 多維資料處理之主成分分析 PCA 我叫平沢唯 部落格園 劉理理 PCA可以從SVD演變過來,SVD其實就是在單位圓上找乙個向量,讓資料X在其上的投影的長度平方和最大化,非中心化時中心偏移了 0,0 這時候也可以找那個最大化投影,但不能準確...