線形システム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.eq5.gif 537件 [詳細] filels_lanzcos2_ls.eq6.gif 527件 [詳細] filels_lanzcos2_ls.eq7.gif 542件 [詳細] filels_lanzcos2_ls.eq8.gif 510件 [詳細] filels_lanzcos2_ls.eq9.gif 509件 [詳細] filels_lanzcos2_ls.eq1.gif 512件 [詳細] filels_lanzcos2_ls.eq10.gif 589件 [詳細] filels_lanzcos2_ls.eq11.gif 515件 [詳細] filels_lanzcos2_ls.eq12.gif 617件 [詳細] filels_lanzcos2_ls.eq13.gif 634件 [詳細] filels_lanzcos2_ls.eq14.gif 599件 [詳細] filels_lanzcos2_ls.eq15.gif 519件 [詳細] filels_lanzcos2_ls.eq16.gif 528件 [詳細] filels_lanzcos2_ls.eq17.gif 550件 [詳細] filels_lanzcos2_ls.eq18.gif 521件 [詳細] filels_lanzcos2_ls.eq19.gif 535件 [詳細] filels_lanzcos2_ls.eq2.gif 527件 [詳細] filels_lanzcos2_ls.eq3.gif 543件 [詳細] filels_lanzcos2_ls.eq4.gif 594件 [詳細]

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