多項式で近似するのが難しい連続関数の例


16

教育目的のために、多項式で近似するのが「難しい」単一変数の連続関数が必要です。つまり、この関数をうまくフィットさせるには、べき級数の非常に高いべき乗が必要です。私は、学生にべき級数で達成できることの「限界」を示すつもりです。

私は「うるさい」と何かを紡ぎ上げについて考えたが、代わりに自分を転がり、私はちょうど人々がそれらに多少同様に、近似/補間アルゴリズムをテストするために使用することを標準「難しい機能」のようなものがあるかどうかを疑問に思って最適化テスト機能数多く持っています素朴なアルゴリズムが簡単にスタックするローカルミニマム。

この質問の形式が適切でない場合はおApびします。非数学者に慈悲を与えてください。

回答:


14

絶対値関数を単に表示しないのはなぜですか?

ルジャンドル多項式展開などの近似は機能しますが、かなりひどいです:

多項式による絶対値関数の逐次近似

テイラー展開はもちろんここでは完全に役に立たず、常に線形関数のみを与え、常に減少または常に増加します(展開する点が負か正かによって異なります)。


補間できます| x | Chebyshev補間を使用して、nbviewer.jupyter.org / github / cpraveen / na / blob / master /…を参照してください。これは非常に高速に収束します。たとえば、コードのN = 2 * iをN = 15 + iに変更して、より大きな度合いをテストできます。これは拡張方法ではありませんが、依然として多項式に基づいています。
-cpraveen

バツ=0L2

チェビシェフ補間のように不均一な点を持っていることはまったく問題ありません。次数が約20の場合、投稿で表示するルジャンドルよりもはるかに正確な近似が得られます。エラーをより定量的に測定します。| x |のチェビシェフ級数近似を行うこともできます。これはルジャンドル展開よりも正確です。
cpraveen

バツ|バツ|

完璧な方法はありません。問題は、多項式を近似するのが難しい関数です。したがって、多項式を含むすべての可能な方法を調べて、それらのどれも良い仕事をしていないと結論付けなければなりません。ルジャンドルは| x |を近似する最良の方法ではありません そのため、多項式は| x |に対してあまりにも悪いという誤った印象を与えます。チェビシェフを使用すると、収束とルジャンドルよりもはるかに優れた近似が得られ、ルジャンドルほどひどく振動しませんが、関数が十分に滑らかではないx = 0付近でゆっくり収束します。
cpraveen


6

近似は、近似される関数によって困難になるだけでなく、近似が「良好な適合」になる間隔によって困難になります。そして、あなたは「良い適合」のための尺度を定義する必要があります。つまり、あなたが許容したい最大(絶対または相対)エラーは何ですか?

expバツ[010]バツ[02π]ここに画像の説明を入力してくださいここに画像の説明を入力してください


テイラー展開は関数を近似する良い方法ではないという点を強調するために、私のコースでそのような例を示します。
cpraveen

6

多項式は、関数近似で驚くほど効果的です[1]。少なくともリプシッツ連続性がある場合、チェビシェフ近似は収束します。もちろん、収束は遅くなる可能性があり、それが非平滑関数を処理するために支払う代償です。

今日、コンピューターは多くの数値分析の本が書かれた時代よりもはるかに高速であり、巧妙なアルゴリズムによってさらに速度が向上しているため、より多くの用語を使用することは以前ほど悪くないかもしれません。

ワイエルシュトラスモンスター機能のような病理学的例は、理論的な観点からは興味深いものですが、ほとんどの実際のアプリケーションコンテキストを代表するものではありません。

|バツ|バツ=0

多項式による近似の難しさを教えることは重要ですが、これらの問題に対処できる誤差推定と適応アルゴリズムを構築できることを生徒に伝えることも重要です。

[1] https://people.maths.ox.ac.uk/trefethen/mythspaper.pdf

[2] http://www.chebfun.org


Lloyd Trefethenによる「神話論文」をリンクしてくれた+1、トピックIMOの非常に良い調査、ありがとう。
ラリックスデシドゥア

2

fバツ=1バツ2+1

1バツ2+1=1バツ2+バツ4バツ6+バツ8バツ10+バツ12

1<バツ<1バツ=0バツ=2


0

y=snバツ?周期関数は、多項式を有限間隔に制限するオプションがない限り、多項式を使用して近似するのが難しいはずです。


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