適応型求積法

適応型求積法は、被積分関数が困難な積分区間を自動的に細分化し、局所誤差推定を用いて、可能な限り少ない関数評価で要求された精度を達成します。

Definition

適応型求積法とは、局所近似誤差の推定値を用いて、積分領域をどこでどのように細分化するかを決定し、所定の全体誤差許容範囲を効率的に達成する数値積分戦略の総称です。

Scope

このトピックでは、異なる次数または細分化レベルの規則を比較することによる局所誤差推定、再帰的な区間二分法（適応型シンプソン法および適応型ガウス・クロンロッド法）、大域的誤差予算と停止基準、特異点および急峻な特徴の処理、QUADPACKライブラリなどの実用的な自動積分器の設計について扱います。

Core questions

• 正確な積分値を知らずに、求積推定値の局所誤差はどのように計算されるのでしょうか？
• 再帰的な細分化は、被積分関数が最も変化する箇所にどのように努力を集中させるのでしょうか？
• 無駄な作業を避けつつ、要求された許容誤差を確実に達成する停止基準は何でしょうか？
• 可積分特異点と不連続点はどのように検出され、ロバストに処理されるのでしょうか？

Key theories

局所誤差推定と細分化

部分区間における粗い推定値とより細かい（または高次の）推定値を比較することで、局所誤差の推定値が得られます。この推定値がその部分区間に割り当てられた許容誤差の割合を超えた場合、部分区間は分割され、手順は再帰的に繰り返されます。そうでない場合、その寄与は受け入れられます。

大域適応型戦略

部分区間を独立して扱うのではなく、大域適応型積分器は推定誤差によって順序付けられた部分区間のキューを保持し、常に最悪のものを細分化します。これにより、局所化された特異点が効率的に処理され、QUADPACKルーチンの基盤となっています。

Mechanisms

各部分区間において、積分器は埋め込み規則ペア（例えば、ガウス・クロンロッドペア、または異なる細分化レベルでの2つのシンプソン推定値）を評価し、その差によって局所誤差を推定します。局所適応型手法は、推定誤差が大きすぎる部分区間を二分することで再帰的に処理します。大域適応型手法は、推定誤差をキーとする部分区間の優先度キューを維持し、合計誤差推定値が許容範囲を満たすまで、現在の最悪の違反者を繰り返し細分化します。端点特異点や振動する被積分関数に対処するために、外挿法や特殊な重み処理が追加されます。

Clinical relevance

適応型求積法は、科学ソフトウェアの汎用積分ルーチンが、ユーザーが被積分関数を分析することなく、ユーザー指定の精度で結果を提供するために依拠するものです。これは、固定規則では対処できないピーク、境界層挙動、または可積分特異点を持つ被積分関数にとって不可欠であり、広く使用されている数値および統計パッケージの自動積分器の基盤となっています。

History

自動誤差制御積分は1970年代から1980年代初頭にかけて成熟し、QUADPACKパッケージ（1983年）で頂点に達しました。その適応型ガウス・クロンロッドルーチンと外挿法は事実上の標準となり、後に多くの数値および統計ソフトウェアシステムで採用、移植、または再実装されました。

Key figures

• Robert Piessens
• Philip J. Davis
• Philip Rabinowitz

Seminal works

• davis1984
• piessens1983

Frequently asked questions

適応型積分器は、答えを知らないのにどのように誤差を知るのでしょうか？

同じ部分区間上で異なる精度の2つの近似（例えば、高次規則と低次規則）を比較することで、局所誤差を推定します。それらの差が誤差を近似し、真の積分値が不明であっても、どこを細分化すべきかを導きます。

適応型求積法はどのような場合に困難に直面しますか？

サンプル点では滑らかだがその間に隠れた特徴を持つ被積分関数、強く振動する被積分関数、または非可積分特異点によって誤った結果を導くことがあります。その場合、特殊な規則、変換、または振動積分法が必要となります。