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

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 887件 [詳細] filerf_newton_md.eq17.gif 894件 [詳細] filerf_newton_md.eq18.gif 825件 [詳細] filerf_newton_md.eq19.gif 758件 [詳細] filerf_newton_md.eq2.gif 747件 [詳細] filerf_newton_md.eq20.gif 764件 [詳細] filerf_newton_md.eq21.gif 843件 [詳細] filerf_newton_md.eq22.gif 909件 [詳細] filerf_newton_md.eq3.gif 796件 [詳細] filerf_newton_md.eq4.gif 704件 [詳細] filerf_newton_md.eq5.gif 761件 [詳細] filerf_newton_md.eq6.gif 756件 [詳細] filerf_newton_md.eq7.gif 783件 [詳細] filerf_newton_md.eq8.gif 782件 [詳細] filerf_newton_md.eq9.gif 708件 [詳細] filerf_newton_md.eq1.gif 758件 [詳細] filerf_newton_md.eq10.gif 737件 [詳細] filerf_newton_md.eq11.gif 701件 [詳細] filerf_newton_md.eq12.gif 740件 [詳細] filerf_newton_md.eq13.gif 795件 [詳細] filerf_newton_md.eq14.gif 752件 [詳細] filerf_newton_md.eq15.gif 766件 [詳細]

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