近似数値微分によるニュートンラプソン近似の欠点


17

関数あり、ようなを見つけたいとします。Newton-Raphson法を使用する場合があります。しかし、これには微分関数を知っている必要があります。の分析式が利用できない場合があります。たとえば、は、実験値のデータベースを参照する複雑なコンピューターコードによって定義される場合があります。X 、F X 0 F 'X F Ffxf(x)0f(x)ff

しかし、が複雑な場合でも、小さな数字を選択してを計算することにより、特定のを近似できます。。、F 'ε F 'F A + ε - F ff(a)aϵf(a)f(a+ϵ)f(a)ϵ

このアプローチには明らかに不利な点があると聞きましたが、それが何であるかはわかりません。ウィキペディアは、「この近似を使用すると、ニュートンの方法よりも収束が遅い割線法のようなものになる」と示唆しています。

誰かがこれについて詳しく説明し、この手法の問題を特に説明するリファレンスを提供できますか?


5
secantメソッドは、導関数の計算が高価な場合の優れた代替手段です。3段割線は一般に2段ニュートンとほぼ同等であり、より安価です。

1
(提案しているように)有限差分によって導関数を数値的に計算するときは常に、関数内のノイズが増幅されるため、イプシロンを慎重に選択する必要があります。1つの可能性は、解に近づいたときに、2分割法に切り替えることです。fが局所的に単調である限り、収束することが保証されています。
マイクダンラベイ

2
アンドレが述べたように、2点の数値微分は、あなたが提案するように、再開されたセカント法と同等です。ただし、収束を高速化するために、いわゆるイリノイアルゴリズムをお勧めします。これは、Secantメソッドに近いもので、ステップごとに1つのポイントのみを使用します(2つではなく)。間違った位置の方法。
ペドロ

の次元は何ですか?次元が高いほど、デリバティブはより価値があります。ヤコビアンフリーのニュートンクリロフは、明示的な微分を必要としないオプションです(ただし、条件の悪いシステムでは前提条件が重要です)。バツ
ジェドブラウン

回答:


12

表記のために、(つまり、ベクトルを入力として受け取り、同じサイズのベクトルを出力するベクトル値関数)を想定します。懸念事項は2つあります。計算コストと数値精度です。fRnRn

誘導体算出(ヤコビ行列をJ X 、又はF X T、または任意のお好み)有限差分を用いては必要としようとしているn個の機能評価を。定義から直接浮動小数点演算を使用して導関数を計算できる場合、差分商を計算する必要がありますDfバツJバツfバツTn

Dfバツe=リムε0fバツ+εefバツε

それぞれについて、のスパースパターンあなたが知っている(または検出することができます)ので、あなたが(カーティス・パウエルレイドのような)「差分スマート有限」の任意の並べ替えをしないと仮定すると、Dの Fnが大きい場合、多くの関数の評価になる可能性があります。D fの分析式がある場合、それを計算する方が安くなる可能性があります。また、自動(アルゴリズムとしても知られている)微分法を使用して、関数評価のコストの約3〜5倍でD fを計算することもできます。=1nDfnDfDf

数値的な懸念もあります。明らかに、コンピューターでは、ゼロになるスカラーの制限をとることができないため、を近似すると、εを実際に「小さい」ものとして計算します。Dfε

Dfバツefバツ+εefバツε

ここで、は近似値であることを意味し、本当に良い近似値であることを願っています。あなたが選ぶ場合ので、浮動小数点演算では、この近似を計算することは難しいですε大きすぎる、あなたの近似が悪いかもしれないが、あなたが選ぶ場合はεが小さすぎ、大きな丸め誤差がある可能性があります。これらの効果は、表面的な詳細における数値微分に関するウィキペディアの記事で説明されています。より詳細なリファレンスは記事内で見つけることができます。εε

ヤコビ行列の誤差が大きすぎない場合、ニュートンラプソンの反復は収束します。詳細な理論解析については、Nick Highamによる精度と安定性の数値アルゴリズムの第25章、またはそれが基づいているFrançoiseTisseurの論文を参照してください。Df

通常、ライブラリはこれらのアルゴリズムの詳細を処理し、通常、Newton-Raphsonアルゴリズム(またはその変形)のライブラリ実装は非常にうまく収束しますが、多くの場合、欠点のために問題を引き起こす問題があります上記。スカラーの場合、その堅牢性と実際の良好な収束率のために、ブレントの方法を使用します。n=1

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.