関数の相対誤差が与えられた場合の導関数の相対誤差の境界


8

関数仮定相対誤差に結合するように計算することができることである、すなわち及び 、それぞれ計算され、正確な値であり及びR f x = f x 1 + r f f f | r | RfRf(x)=f(x)(1+r)fff|r|R

次の微分近似の相対誤差を、一般的なに対してとでしたいR fhRf

f(x)f(x+h)f(xh)2hf(x)f(x+h)2f(x)+f(xh)h2

ラルストンとラビノウィッツでは、境界はそれぞれRh4Rh2与えられます。しかし、これは証明されておらず、リチャードソン外挿についての説明の一部として言及されています。

その証明に関するアイデアはありますか?


5
指定した式には、エラーに両方の項が含まれていませんの評価が不正確であることと、打ち切りエラーが原因でエラーが発生する可能性があります(が大きすぎます)極端な場合(関数の正確な評価)、リストした数式は0エラーになりますが、説明する切り捨てエラーはまだあります。打ち切り誤差と不正確な関数評価による誤差の式を導き出し、合計誤差がどのように変化するかを確認することは良い練習です。fhR=0h
Brian Borchers 2013年

1
この境界について、より具体的なステートメント、または少なくとも正確な参照(定理とページ)を提供できると助かります。
クリスチャンクラソン2013年

回答:


1

この定理は、質問の作成者によって誤って解釈されたか、参照された本に誤りがあります。次のカウンターの例を考えてみます。

f(x)=100+x
h=0.01
R=0.01

で各機能の評価における絶対誤差は、我々は持っているので、 最悪のシナリオでは、2つのエラー項は同じ符号を持ち、キャンセルされません。したがって、微分近似の相対誤差はにもなり、よりはるかに大きくなります。x=01000.01=1

f(0)=f(h)f(h)2h=(100+0.01)±1((1000.01)±1)0.02
f(0)=0.020.02±10.02±10.02
100R/h=1

形式関数を選択することにより、導関数近似の相対誤差は常に増やすだけで増やすことができるため、私が知る限り、一般的な相対誤差には制限がありません。。ff(x)=n+xn

一方、依存する範囲を計算できます。十分に小さいとの絶対誤差の限界は次のとおりです 証明: ここで、周りでをテイラー展開し、と両方が小さいため、または以上の次項を無視します。同様に、 したがって、 fhR

f(x)=f(x+h)f(xh)2h±f(x)Rh
f(x+h)=(f(x)+hf(x))(1±R)
f(x+h)=f(x)+hf(x)±Rf(x)
fxhRh2hR
f(xh)=f(x)hf(x)±Rf(x)
f(x)=2hf(x)±Rf(x)±Rf(x)2h
f(x)=f(x)±Rf(x)h
ここで、エラーが増加する最悪のシナリオをもう一度検討します。

したがって、相対誤差の範囲は依存し、 f(x)

f(x)=f(x)(1±Rf(x)hf(x))

同様に、

f(x)=f(x)(1±4Rf(x)h2f(x))

もちろん場合、の打ち切り誤差は項に比べて支配的になりますhRO(h2)O(R/h)
SimonSciComp

ああ、私はこの賞金を完全に忘れていました!とにかく、これはかなりうまくカバーしているようです。
David Z

-1

あなたの直接の質問に答えるには(そしてブライアンボーチャーズのコメントの切り捨てエラーには対応しない):

の定義により、その相対誤差であり、定義はそれを明示的に示していませんが、は定数ではないため、はです。f|(ff)/f|Rr|f(x+h)f(xh)|2R

これは、がになる相対誤差に直接つながり、同様にがになる相対誤差につながります。fR/hf4R/h2


2
これにより、値に 限界が生じますではなく、上のエラー|f|f |ff|
クリスチャンクラソン2013年
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.