主題
Search

共軛梯度平方方法


雙共軛梯度方法 中,殘差向量 r^((i)) 可以被視為 r^((0)) 和關於 Ai 次多項式的乘積,即:

 r^((i))=P_i(A)r^((0)).
(1)

這個相同的多項式滿足

 r^~^((i))=P_i(A^(T))r^~^((0))
(2)

因此

rho_i=(r^~^((i)),r^((i)))
(3)
=(P_i(A^(T))r^~^((0)),P_i(A)r^((0)))
(4)
=(r^~^((0)),P_i^2(A)r^((0))).
(5)

這表明,如果 P_i(A)r^((0)) 縮小為一個較小的向量 r^((i)),那麼將這個“收縮”運算元應用兩次並計算 P_i^2(A)r^((0)) 可能是有利的。迭代係數仍然可以從這些向量中恢復(如上所示),並且事實證明很容易找到 x 的相應近似值。這種方法是共軛梯度平方(CGS)方法(Sonneveld 1989)。

通常人們觀察到 CGS 的收斂速度大約是 雙共軛梯度方法 的兩倍,這與同一個“收縮”運算元應用兩次的觀察結果相符。然而,沒有理由認為收縮運算元,即使它真的減少了初始殘差 r^((0)),也應該減少一次減少後的向量 r^((k))=P_k(A)r^((0))。CGS 經常出現高度不規則的收斂行為就證明了這一點。應該意識到,對當前解的區域性校正可能非常大,以至於會發生抵消效應。這可能導致解的精度不如更新後的殘差所暗示的那麼高(van der Vorst 1992)。如果初始猜測接近解,則該方法傾向於發散。

CGS 每次迭代所需的操作次數與 雙共軛梯度方法 大致相同,但不涉及與 A^(T) 的計算。因此,在與 A^(T) 的計算不切實際的情況下,CGS 可能很有吸引力。


另請參閱

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

此條目由 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. 線性系統求解模板:迭代方法的構建塊,第二版 Philadelphia, PA: SIAM, 1994. http://www.netlib.org/linalg/html_templates/Templates.html.Sonneveld, P. "CGS:非對稱線性系統的快速 Lanczos 型求解器。" SIAM 科學與統計計算雜誌 10, 36-52, 1989.van der Vorst, H. "Bi-CGSTAB:Bi-CG 的快速平滑收斂變體,用於求解非對稱線性系統。" SIAM 科學與統計計算雜誌 13, 631-644, 1992.

在 中被引用

共軛梯度平方方法

請引用為

Black, NoelMoore, Shirley. "共軛梯度平方方法." 來自 Web 資源,由 Eric W. Weisstein 建立. https://mathworld.tw/ConjugateGradientSquaredMethod.html

主題分類