積分の精度を推定する方法は?


11

コンピュータグラフィックスで非常に一般的な状況は、一部のピクセルの色が一部の実数値関数の積分に等しいことです。多くの場合、関数は複雑すぎて分析的に解くことができないため、数値近似を行います。しかし、この関数は計算に非常にコストがかかることも多いため、計算できるサンプルの数には大きな制約があります。(たとえば、100万サンプルを取得して、そのままにしておくことはできません。)

次に、一般的に、推定積分が「十分に正確」になるまで、ランダムに選択されたポイントで関数を評価します。これは私の実際の質問に私をもたらします:積分の「精度」をどのように推定しますか?


具体的には、があります。これは、いくつかの複雑で低速なコンピューターアルゴリズムによって実装されます。見積もりたいf:RR

k=abf(x) dx

任意のxに対してを計算できますが、コストがかかります。したがって、ランダムにいくつかのx値を選択し、kの推定が許容できるほど正確になったときに停止します。もちろん、これを行うには、現在の見積もりが実際にどれほど正確であるかを知る必要があります。f(x)xxk

この種の問題にどの統計ツールが適切であるかさえわかりません。しかし、私がfについてまったく何も知らなければ、問題は解決できないようです。たとえば、f x を1000回計算し、それが常にゼロの場合、推定積分はゼロになります。しかし、については何も知りませんfは、それがあることはまだ可能だF xは= 1 000 000をお使いの推定値は恐ろしく間違っているので、あなたは、サンプルに起こった点を除いてどこでも!ff(x)ff(x)=1,000,000

ff


編集: OK、これは多くの応答を生成したようで、これは良いことです。それぞれに個別に返信するのではなく、ここで追加の背景を記入してみます。

ffff

fff

f

また、「モンテカルロ」が出現した回数を考えると、それがこの種の統合の専門用語だと思いますか?


ff

2
通常、既知の関数を介して統合すると、モンテカルロ統合よりもはるかに優れた機能を実行できます。モンテカルロは割合で真の値に収束します1/NN1/N(lnN)n/Nn

1
f

1
ff

1
f

回答:


2

簡単にするために、すべてのx in [a、b]についてf(x)> = 0と仮定すると、[a、b]のすべてのxについてf(x)<MとなるMがわかります。[a、b]上のfの積分Iは、幅baおよび高さMの長方形で囲むことができます。fの積分は、関数fに該当する長方形にM(ba)を掛けた比率です。長方形の点をランダムに選択し、曲線の下にある場合は成功として、そうでない場合は失敗としてカウントすると、ベルヌーイ試行が設定されます。内部のポイントのサンプル部分は二項比率であり、したがって、平均pと分散p(1-p)/ nを持ちます。ここで、nは取得したポイントの数です。したがって、pの信頼区間を構築できます。また、I = p M(ba)なので、推定のI ^ = p ^ M(ba)、Var(I ^)= Mなので、Iの信頼区間です。222222


3
0Mf

1
@Macro fについて何も知らなければ、固定有限点セットでの評価に基づいて積分の推定値の統計的精度についてどのように言えるかわかりません。私の仮定はかなり最小限です。fが区間[a、b]で制限されている場合、fの上限として使用できる十分な大きさのMが必要です。
Michael R. Chernick

M

2
それは仮定です。決定的な答えに到達するために、できるだけ少ない仮定をしていると言うために、「ミミラル」という用語を使用しました。
Michael R. Chernick

f

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