メインメニューを開く

ニュートン法

方程式系を数値計算によって解くための反復法による求根アルゴリズムの1つ

数値解析の分野において、ニュートン法(ニュートンほう、: Newton's method)またはニュートン・ラフソン法: Newton-Raphson method)は、方程式系を数値計算によって解くための反復法による求根アルゴリズムの1つである。対象とする方程式系に対する条件は、領域における微分可能性と2次微分に関する符号だけであり、線型性などは特に要求しない。収束の速さも2次収束なので古くから数値計算で使用されていた。名称はアイザック・ニュートンジョゼフ・ラフソンに由来する。

導入編集

 
ニュートン法の一手順の概念図 (青い線が関数 f のグラフで、その接線を赤で示した). xn よりも xn+1 のほうが、 f(x)=0 の解 x についてのよりよい近似を与えている.

この方法の考え方は以下のようである:まず初めに、予想される真の解に近いと思われる値をひとつとる。次に、そこでグラフの接線を考え、その x 切片を計算する。このx切片の値は、予想される真の解により近いものとなるのが一般である。以後、この値に対してそこでグラフの接線を考え、同じ操作を繰り返していく。

上の考え方は次のように定式化される。 ここでは、考える問題を f: RR, xRとして

 

となる x を求めることに限定する。このとき、x の付近に適当な値 x0 をとり、次の漸化式によって、x収束する数列を得ることができる場合が多い。

 

例として、2 を計算で求める場合に、

 

とおき、f(x) = 0 の解を求めることを考える。

 

であるので、(1) の式は

 

と書き表せる。たとえば x0 = 1 とおくと、この数列は √2 に収束するが、その収束の仕方は

x0 = 1
x1 = 1.5
x2 = 1.416666…
x3 = 1.4142156862745…
x4 = 1.4142135623746899… (下線部は√2と一致する部分)

となる。

また、x0 = −1 とおくと、この数列は −√2 に収束する。

理論編集

局所収束定理編集

初期値 を解の十分近くに選ぶことを要求した上で、

 

の解を考える(解の存在を仮定する)。    でのテーラー展開をすると

 

このとき、(右辺)=0の解は、(左辺)=0の根の  での多項式次数一次の近似となっている。 右辺の解は

 

次に、この近似値が、 より根に近づいている ということに関する意味を考える。 上式を、次のような離散力学系として考える。

 , ただし  

この力学系において、  となる は明らかに固定点である。

したがって 、つまり が沈点(アトラクター、安定固定点)であり、 与えられた初期条件 が、このアトラクターの吸引領域に属していれば   -極限( )は  となる に収束する。

 が沈点である保証は、常に担保されてはいない。 例えばx軸の漸近線や関数 極値近傍では固定点が不安定になる事が知られている。 [1]

たとえば を、適当な近傍の点 で展開すると なら、二次の余剰項 として

 

よって2次で収束する。

半局所収束定理編集

前節では解の存在を仮定した上で初期値 を解の十分近くに選ぶことを要求した。これに対して、解の存在を仮定せず、初期値 がある条件を満たすときに解の存在と反復の収束を示す定理を半局所収束定理(Semi-Local Convergence Theorem)という。1次元の場合での半局所収束定理はコーシーによって1829年に示され、 次元ユークリッド空間での場合はファイン(Fine)によって1916年に示された。その後、バナッハ空間での半局所収束定理がカントロビッチ(Kantorovich)によって1948年に示され、現代ではニュートン-カントロビッチの定理と呼ばれている[2]。この定理にはいくつかの変種が知られており、[3]にまとめられている。

高次元の場合編集

ニュートン法は、接線を一次近似式、接線のx切片を一次近似式の零点と考えることにより、より高次元の関数の場合に一般化できる。 対象となる関数を f: RmRm, xRm とし、

 

なる点 x を求めるには次のようにする。(f が同じ次元の空間の間の関数であることに注意せよ。)

まず、今 x0Rm が既知であるとする。x0における f(x) の一次近似式

 

を考える。ただし、∂f(x0) は、m × mヤコビ行列である。

この一次近似式の零点を求める。ヤコビ行列∂f(x0) が正則行列であるとして、

 

を解いて、

 

となる。

コンピュータで計算を行う場合 ∂f(x0)-1f(x0) を直接求めることは困難なので、

 

という方程式をガウスの消去法などの解法によって線形方程式系を解き r を求め、x = x0 - r によって x を求める。

ここで求めた xx0よりも f(x) = 0 の解に近いことが見込まれる。そこで、今求めた xx1 として、再び同様の計算を繰り返す。計算を繰り返すことによって xnf(x) = 0 の解に近づいていく。

逆行列を求めることを避けるために共役勾配法を用いることがある。

注意編集

ニュートン法により近似値を求めようとする場合にはヤコビ行列が陽に分からなければ計算できない。しかし、関数 f によってはヤコビ行列が陽に分からない場合もある。この場合にはヤコビ行列を必要としない準ニュートン法を用いる。

また、f (x*) = 0 を満たす真の解 x* において導関数が f '(x*) = 0 であった場合、収束の速さは1次に退化する[4]

改良編集

ニュートン法の考え方を少し改良することにより、(1) の代わりに次の式を用いる方法を得る。

 

この方法は、場合によっては従来の方法より速い[5]

平野の変形ニュートン法編集

ニュートン法の局所的な2次収束性を保ちつつ、不安定な振る舞いを抑えるように工夫した方法として平野の変形ニュートン法が知られている[6]

簡易ニュートン法編集

ニュートン法における導関数の計算を簡略化したものが簡易ニュートン法である:

 

簡易ニュートン法に対する半局所収束定理は占部の定理として知られる。占部の定理は元々数学的帰納法を使って示されたが、その後、バナッハの不動点定理を使った別証明が与えられた[7]

クラフチック法編集

簡易ニュートン法に対する半局所収束定理の条件をより容易に評価するために開発された簡易ニュートン法の変種がクラフチック(Krawczyk)法である[8][9]

区間ニュートン法編集

ニュートン法に対する半局所収束定理の条件をより容易に評価するために開発されたニュートン法の変種が区間ニュートン法である[10]

 -ニュートン法編集

ニュートン法において微分の代わりに微分のq-類似 -微分)を使う反復を -ニュートン法という[11]:

 

 -ニュートン法に対する半局所収束定理が -ニュートン-カントロビッチの定理である[12]

脚注編集

  1. ^ Newton's Method(Wolfram Math World), http://mathworld.wolfram.com/NewtonsMethod.html 2010年6月8日閲覧。 
  2. ^ 数値解析入門(増訂版)、山本哲郎、サイエンス社、2003年。
  3. ^ J. Ortega & W. Rheinboldt, Iterative Solution of Nonlinear Equations in Several Variables, Academic Press (1970).
  4. ^ 小澤一文 『Cで学ぶ数値計算アルゴリズム』 共立出版、2008年、49頁。ISBN 978-4-320-12221-5 
  5. ^ 長島隆廣「ニュートン近似の改良」数学セミナー, 第19巻, 第5号, 通巻222号, p.112, 1980年5月, 日本評論社.
  6. ^ 室田一雄「平野の変形Newton法の大域的収束性」数理解析研究所講究録 (1981), 422: 104-118
  7. ^ 非線形解析入門、大石進一コロナ社、1997年。
  8. ^ 精度保証付き数値計算、大石進一コロナ社、2000年。
  9. ^ 大石進一 et. al. (2018), 精度保証付き数値計算の基礎, コロナ社.
  10. ^ Interval Newton Method
  11. ^ Rajković, P., Stanković, M., & Marinković D, S. (2002). Mean value theorems in g-calculus. Matematički vesnik, 54(3-4), 171-178.
  12. ^ Rajković, P. M., Marinković, S. D., & Stanković, M. S. (2005). On q-Newton–Kantorovich method for solving systems of equations. Applied Mathematics and Computation, 168(2), 1432-1448.

関連項目編集

外部リンク編集