主題
Search

對稱 LQ 方法


共軛梯度法可以看作是 Lanczos 方法對於正定對稱系統的特殊變體。最小殘差法和對稱 LQ 方法 (SYMMLQ) 是可以應用於對稱不定系統的變體。

共軛梯度法中的向量序列對應於類似於係數矩陣的三對角矩陣的分解。因此,如果矩陣是不定的,則演算法可能會發生崩潰,對應於零主元。此外,對於不定矩陣,共軛梯度法的最小化性質不再明確定義。MINRES 和 SYMMLQ 方法是 CG 方法的變體,它們避免了 LU 分解,並且不會遭受崩潰。SYMMLQ 求解投影系統,但不最小化任何東西(它保持殘差與所有先前的殘差正交)。

A 不是正定的,而是對稱的時,我們仍然可以透過三項遞推關係為 Krylov 子空間構造一個正交基。消除共軛梯度法方程中的搜尋方向得到一個遞推式

 Ar^((i))=r^((i+1))t_(i+1,i)+r^((i))t_(i,i)+r^((i-1))t_(i-1,i),
(1)

可以寫成矩陣形式為

 AR_i=R_(i+1)T^__i,
(2)

其中 T^__i 是一個 (i+1)×i 三對角矩陣

在這種情況下,我們遇到了 (·,·)_(A) 不再定義內積的問題。但是,我們仍然可以透過獲得以下方式嘗試最小化 2-範數中的殘差

 x^((i)) in {r^((0)),Ar^((0)),...,A^(i-1)r^((0))},    x^((i))=R_iy^_
(3)

從而最小化

|Ax^((i))-b|_2=|AR_iy^_-b|_2
(4)
=|R_(i+1)T^__iy-b|_2.
(5)

現在我們利用以下事實,如果

 D_(i+1)=diag(|r^((0))|_2,|r^((1))|_2,...,|r^((i))|_2),
(6)

那麼 R_(i+1)D_(i+1)^(-1) 是關於當前 Krylov 子空間的正交變換

 |Ax^((i))-b|_2=|D_(i+1)T^__iy-|r^((0))|_2e^((1))|_2,
(7)

並且這個最終表示式可以簡單地看作是最小范數最小二乘問題。

一種方法是求解系統 T_iy=|r^((0))|_2e^((1)),就像在共軛梯度法中一樣(T_iT^__i 的上 i×i 部分)。然而,與共軛梯度法不同,我們不能依賴 Cholesky 分解的存在(因為 A 不是正定的)。另一種方法是透過 LQ 分解來分解 T_i。這導致簡單的遞推關係,並且所得方法被稱為 SYMMLQ (Paige and Saunders 1975)。


另請參閱

雙共軛梯度法, 切比雪夫迭代, 正規方程上的共軛梯度法 共軛梯度法, 共軛梯度平方方法, 靈活的廣義最小殘差法, 廣義最小殘差法, 線性方程組, 最小殘差法, 非平穩迭代法, 預處理器, 準最小殘差法 平穩迭代法, 無轉置準最小殘差法

此條目由 Noel Black 和 Shirley Moore 貢獻,改編自 Barrett 等人 (1994) (作者連結)

使用 探索

參考文獻

Barrett, R.; Berry, M.; Chan, T. F.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; 和 van der Vorst, H. 線性系統求解模板:迭代方法的構建模組,第 2 版。 Philadelphia, PA: SIAM, 1994. http://www.netlib.org/linalg/html_templates/Templates.html.Paige, C.; Parlett, B.; 和 van der Vorst, H. "來自 Krylov 子空間的近似解和特徵值界限。" Numer. Lin. Alg. Appl. 29, 115-134, 1995.Paige, C. 和 Saunders, M. "稀疏不定線性方程組的解。" SIAM J. Numer. Anal. 12, 617-629, 1975.

在 上被引用

對稱 LQ 方法

引用為

Black, NoelMoore, Shirley. "對稱 LQ 方法。" 來自 Web 資源,由 Eric W. Weisstein 建立。 https://mathworld.tw/SymmetricLQMethod.html

主題分類