線形システムls_lanzcos2_ls.eq2.gif(ただし,Aは対称行列)について,初期近似値ls_lanzcos2_ls.eq3.gifが与えられたとき, mステップ目の近似値ls_lanzcos2_ls.eq4.gifはFOMと同様に以下のようになる.

ls_lanzcos2_ls.eq5.gif

ここで,ls_lanzcos2_ls.eq6.gifが三重対角行列(tridiagonal matrix)になるのでls_lanzcos2_ls.eq7.gifと置き換えている. また,ls_lanzcos2_ls.eq8.gifである.

この式からLanczos法を使った線形システムの解法アルゴリズムは以下となる.

ls_lanzcos2_ls.eq9.gifを計算
ls_lanzcos2_ls.eq10.gifを設定
for(j = 1,2,...,m){
  ls_lanzcos2_ls.eq11.gif
  ls_lanzcos2_ls.eq12.gif
  ls_lanzcos2_ls.eq13.gif
  ls_lanzcos2_ls.eq14.gif
  if(ls_lanzcos2_ls.eq15.gif) 反復終了
  ls_lanzcos2_ls.eq16.gif
}
ls_lanzcos2_ls.eq17.gifで構成される三重対角行列ls_lanzcos2_ls.eq7.gifを設定
ls_lanzcos2_ls.eq18.gif
ls_lanzcos2_ls.eq19.gif


添付ファイル: filels_lanzcos2_ls.eq12.gif 504件 [詳細] filels_lanzcos2_ls.eq14.gif 495件 [詳細] filels_lanzcos2_ls.eq13.gif 520件 [詳細] filels_lanzcos2_ls.eq10.gif 476件 [詳細] filels_lanzcos2_ls.eq4.gif 478件 [詳細] filels_lanzcos2_ls.eq3.gif 461件 [詳細] filels_lanzcos2_ls.eq18.gif 436件 [詳細] filels_lanzcos2_ls.eq19.gif 446件 [詳細] filels_lanzcos2_ls.eq6.gif 438件 [詳細] filels_lanzcos2_ls.eq11.gif 428件 [詳細] filels_lanzcos2_ls.eq17.gif 464件 [詳細] filels_lanzcos2_ls.eq16.gif 436件 [詳細] filels_lanzcos2_ls.eq1.gif 427件 [詳細] filels_lanzcos2_ls.eq15.gif 430件 [詳細] filels_lanzcos2_ls.eq9.gif 425件 [詳細] filels_lanzcos2_ls.eq7.gif 459件 [詳細] filels_lanzcos2_ls.eq2.gif 446件 [詳細] filels_lanzcos2_ls.eq5.gif 449件 [詳細] filels_lanzcos2_ls.eq8.gif 424件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2012-07-12 (木) 17:03:57 (3308d)