共役勾配法の計算手順において,残差ls_krylov.eq1.gifを用いた. この残差はどこから来たのかを考える. まず,ls_krylov.eq2.gif

ls_krylov.eq3.gif

とかける.これを反復式とすると,

ls_krylov.eq4.gif

これで残差ベクトルがでてきた. さらに残差ベクトル間の関係を調べる.

ls_krylov.eq5.gif

よって,

ls_krylov.eq6.gif

同様にls_krylov.eq7.gifに関しても,

ls_krylov.eq8.gif

これらの式から,ls_krylov.eq9.gifls_krylov.eq10.gifls_krylov.eq11.gif の線形結合で表されることがわかる.これを式にすると,

ls_krylov.eq12.gif
ls_krylov.eq13.gif

ここで,ls_krylov.eq14.gifはベクトルls_krylov.eq15.gifの線形結合の集合で表される部分空間であり, 上式のような部分空間をクリロフ部分空間と呼ぶ.

ls_krylov.eq16.gif

ls_krylov.eq17.gifの次元nは近似解を求めるための反復ごとに増えていく. そして,クリロフ部分空間内の任意の点はls_krylov.eq18.gifと書ける. ここでls_krylov.eq19.gifは次元がm-1以下の多項式を表している. つまり,クリロフ部分空間内の点はAに関するm-1次以下の多項式とls_krylov.eq20.gifの積の形で書き表せる.

ls_krylov.eq21.gifからクリロフ部分空間ls_krylov.eq17.gifの中を探索することで解を得る非定常な反復解法のことを クリロフ部分空間法と呼ぶ.共役勾配法もクリロフ部分空間法のひとつである (ヤコビ反復やガウス・ザイデルは定常な反復解法).

クリロフ部分空間法としては他に,

  • 双共役勾配法(Bi-Conjugate Gradient method : BiCG法)
  • 安定化双共役勾配法(Bi-Conjugate Gradient STABilized method : BiCGSTAB法)
  • 自乗共役勾配法(Conjugate Gradiate Squared method : CGS法)
  • 共役残差法(Conjugate Residual method : CR法)
  • 一般化共役残差法(Generalized Conjugate Residual method : GCR法)
  • 一般化最小残差法(Generalized Minimal RESidual method : GMRES法) などが提案されている.

添付ファイル: filels_krylov.eq2.gif 586件 [詳細] filels_krylov.eq20.gif 589件 [詳細] filels_krylov.eq21.gif 570件 [詳細] filels_krylov.eq3.gif 590件 [詳細] filels_krylov.eq4.gif 644件 [詳細] filels_krylov.eq5.gif 585件 [詳細] filels_krylov.eq6.gif 686件 [詳細] filels_krylov.eq7.gif 562件 [詳細] filels_krylov.eq8.gif 581件 [詳細] filels_krylov.eq9.gif 733件 [詳細] filels_krylov.eq11.gif 587件 [詳細] filels_krylov.eq1.gif 621件 [詳細] filels_krylov.eq10.gif 607件 [詳細] filels_krylov.eq12.gif 556件 [詳細] filels_krylov.eq13.gif 594件 [詳細] filels_krylov.eq14.gif 975件 [詳細] filels_krylov.eq15.gif 664件 [詳細] filels_krylov.eq16.gif 596件 [詳細] filels_krylov.eq17.gif 696件 [詳細] filels_krylov.eq18.gif 613件 [詳細] filels_krylov.eq19.gif 745件 [詳細]

トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2024-03-08 (金) 18:06:04