前章では、2元の非線型連立方程式のニュートン法での計算方法を示した。こ
こでは、それを一般化する。ここで示す方法は、複素数解にも適用できる。
N元の非線型連立方程式は、
と書くことができる。未知数は、
![$\displaystyle \boldsymbol{X}=(x_1,x_2,x_3,\cdots+x_N)$](img38.png) |
(12) |
とベクトルで表現する。すると、
![$ i$](img18.png)
番目の方程式は、
![$ f_i(\boldsymbol{X})$](img39.png)
と書き表
されるので、表現が簡単になる。これを、先ほどと同じようにテイラー展開す
ると
![$\displaystyle f_i(\boldsymbol{X}+\Delta\boldsymbol{X})= f_i(\boldsymbol{X})+ \f...
... \cdots \frac{\partial f_i}{\partial x_N}\Delta x_N+ O(\Delta \boldsymbol{X}^2)$](img40.png) |
(13) |
となる。
![$ i=1,2,3,\cdots,N$](img41.png)
の全てににおいて、
![$ f_i(\boldsymbol{X}+\Delta\boldsymbol{X})=0$](img42.png)
になるように、
![$ \Delta\boldsymbol{X}=(\Delta x_1,
\Delta x_2, \Delta x_3,\cdots ,\Delta x_N)$](img43.png)
を選ぶ。そのように選ぶため
には、2次以降の高次の項を無視すると
![$\displaystyle \begin{pmatrix}\frac{\partial f_1}{\partial x_1} & \frac{\partial...
...bol{X}) \\ -f_3(\boldsymbol{X}) \\ \ldots \\ -f_N(\boldsymbol{X}) \end{pmatrix}$](img44.png) |
(14) |
の線形であるN元1次連立方程式が成り立つ。これを解いて、
![$ \Delta\boldsymbol{X}=(\Delta x_1,
\Delta x_2, \Delta x_3,\cdots ,\Delta x_N)$](img43.png)
を求める。そうすると、より
真の解に近い
![$ \boldsymbol{X}^{new}$](img45.png)
は、
![$ \boldsymbol{X}^{new}=\boldsymbol{X}^{old}+\Delta
\boldsymbol{X}$](img46.png)
と計算できる。しつこいようであるが、成分で書き表すと
となる
非線型の連立方程式を線形の連立方程式で計算しているわけである。解きやす
い式になった分、反復計算が必要となっている。
ホームページ:
Yamamoto's laboratory著者:
山本昌志
Yamamoto Masashi
平成16年9月13日