ニュートン法を連立非線形方程式に一般化する.

rf_newton_md.eq1.gif

ベクトル表記では,

rf_newton_md.eq2.gif

ここで,

rf_newton_md.eq3.gif

である.

rf_newton_md.eq4.gifステップ目の近似値をrf_newton_md.eq5.gifとし, rf_newton_md.eq5.gifの周りで上式をテイラー展開する.

rf_newton_md.eq6.gif

ここで,rf_newton_md.eq7.gifrf_newton_md.eq8.gifrf_newton_md.eq9.gifの要素とするヤコビ行列である. 2次以上の項を無視すると,連立非線形方程式は以下となる.

rf_newton_md.eq10.gif

rf_newton_md.eq11.gifとすると,

rf_newton_md.eq12.gif

が得られる. この式はrf_newton_md.eq13.gifを未知数とした線形連立方程式であり, LU分解などで解くことで,rf_newton_md.eq13.gifが得られる. そして,以下の式でrf_newton_md.eq14.gifを計算する.

rf_newton_md.eq15.gif

例)2元連立非線形方程式の場合

rf_newton_md.eq16.gif

そして,

rf_newton_md.eq17.gif

である.よって,rf_newton_md.eq18.gifに関する式は以下となる.

rf_newton_md.eq19.gif

rf_newton_md.eq13.gifについて解くと,

rf_newton_md.eq20.gif

これらを用いてrf_newton_md.eq21.gifを更新する.

rf_newton_md.eq22.gif

添付ファイル: filerf_newton_md.eq16.gif 884件 [詳細] filerf_newton_md.eq17.gif 891件 [詳細] filerf_newton_md.eq18.gif 822件 [詳細] filerf_newton_md.eq19.gif 755件 [詳細] filerf_newton_md.eq2.gif 742件 [詳細] filerf_newton_md.eq20.gif 761件 [詳細] filerf_newton_md.eq21.gif 840件 [詳細] filerf_newton_md.eq22.gif 905件 [詳細] filerf_newton_md.eq3.gif 793件 [詳細] filerf_newton_md.eq4.gif 699件 [詳細] filerf_newton_md.eq5.gif 757件 [詳細] filerf_newton_md.eq6.gif 753件 [詳細] filerf_newton_md.eq7.gif 780件 [詳細] filerf_newton_md.eq8.gif 778件 [詳細] filerf_newton_md.eq9.gif 705件 [詳細] filerf_newton_md.eq1.gif 754件 [詳細] filerf_newton_md.eq10.gif 734件 [詳細] filerf_newton_md.eq11.gif 695件 [詳細] filerf_newton_md.eq12.gif 735件 [詳細] filerf_newton_md.eq13.gif 792件 [詳細] filerf_newton_md.eq14.gif 749件 [詳細] filerf_newton_md.eq15.gif 763件 [詳細]

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