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

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 885件 [詳細] filerf_newton_md.eq17.gif 892件 [詳細] filerf_newton_md.eq18.gif 823件 [詳細] filerf_newton_md.eq19.gif 756件 [詳細] filerf_newton_md.eq2.gif 745件 [詳細] filerf_newton_md.eq20.gif 762件 [詳細] filerf_newton_md.eq21.gif 842件 [詳細] filerf_newton_md.eq22.gif 907件 [詳細] filerf_newton_md.eq3.gif 794件 [詳細] filerf_newton_md.eq4.gif 701件 [詳細] filerf_newton_md.eq5.gif 758件 [詳細] filerf_newton_md.eq6.gif 753件 [詳細] filerf_newton_md.eq7.gif 780件 [詳細] filerf_newton_md.eq8.gif 779件 [詳細] filerf_newton_md.eq9.gif 705件 [詳細] filerf_newton_md.eq1.gif 756件 [詳細] filerf_newton_md.eq10.gif 735件 [詳細] filerf_newton_md.eq11.gif 698件 [詳細] filerf_newton_md.eq12.gif 738件 [詳細] filerf_newton_md.eq13.gif 793件 [詳細] filerf_newton_md.eq14.gif 750件 [詳細] filerf_newton_md.eq15.gif 764件 [詳細]

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