ガウス求積

ガウス求積法は、求積規則のノードと重みの両方を選択することで、その多項式の厳密次数を最大化し、わずかn回の関数評価で2n-1次の多項式を正確に積分します。

Definition

ガウス求積法は、重み関数に関連付けられた直交多項式の根をノードとし、与えられたノード数に対して可能な限り最大の厳密次数を達成するために、それらの重みとともに選択される求積規則のファミリーです。

Scope

このトピックでは、直交多項式の根からのガウス規則の構築、ガウス-ルジャンドル規則および重み付き変種（ガウス-チェビシェフ、ガウス-エルミート、ガウス-ラゲール）、ノードと重みを計算するためのゴラブ-ウェルシュ固有値アルゴリズム、および実用的な誤差推定に使用されるガウス-クロンロッド拡張について説明します。

Core questions

• 直交多項式の根にノードを配置すると、固定ノード規則と比較して厳密次数がどのように2倍になるのでしょうか？
• 与えられた重み関数に対して、ノードと重みはどのように正確に計算されるのでしょうか？
• 重み付きガウス規則は、特異な重み関数や無限領域の重み関数を持つ積分をどのように処理するのでしょうか？
• ガウス-クロンロッドペアなどを通じて、信頼性の高い誤差推定はどのように得られるのでしょうか？

Key theories

厳密次数の最大化

n点求積規則は、2n-1次までの多項式に対して厳密であり得、この最大値は、ノードが重み関数に対するn次直交多項式の根であり、すべての重みが正である場合に正確に達成されます。

ゴラブ-ウェルシュアルゴリズム

ガウス規則のノードと重みは、直交多項式の漸化係数から形成される対称三重対角ヤコビ行列の固有値と、固有ベクトルの最初の成分の二乗として得られ、求積法の構築を固有値計算に変換します。

Mechanisms

直交多項式は、対称な三重対角ヤコビ行列を構成する係数を持つ3項漸化式を満たします。ゴラブ-ウェルシュアルゴリズムは、その固有値（求積ノード）を計算し、固有ベクトルの最初の成分を使用して重みを安定して回復します。重み関数を、特異点を持つもの、または半直線や全直線上でサポートされるものに変更すると、ガウス-チェビシェフ、ガウス-ラゲール、またはガウス-エルミート規則が得られ、解析的に困難な挙動を吸収します。ガウス-クロンロッド規則は、ガウスノードを再利用し、インターレースノードを追加することで、わずかな追加コストで高次の推定値、ひいては誤差推定値を得ることができます。

Clinical relevance

ガウス求積法は、有限要素解析における要素積分および剛性積分の評価、統計学および不確実性定量化における確率重み関数に対するモーメントおよび期待値の計算、そして高価な被積分関数評価の数を最小限に抑えることが最も重要である物理学および工学全体における滑らかな積分の高精度評価において、主要なツールとして用いられています。

History

ガウスは1814年に最適な求積法を導出し、ヤコビはそれを直交多項式に結びつけました。現代の計算処理は、1969年のゴラブ-ウェルシュアルゴリズムによって確立され、これによりノードと重みが日常的に計算可能となり、ガウス規則が標準的な数値ライブラリに組み込まれるようになりました。

Key figures

• Carl Friedrich Gauss
• Carl Gustav Jacob Jacobi
• Gene H. Golub
• Walter Gautschi

Seminal works

• davis1984
• gautschi2004

Frequently asked questions

n個の点で2n-1次の多項式を正確に積分できるのはなぜですか？

n個のノードとn個の重みの両方が自由パラメータであるため、2n個の自由度があり、2n個の基底多項式（0次から2n-1次）の積分に一致させるのに十分です。ノードを直交多項式の根に配置することで、まさにこれが達成されます。

ガウス規則の精度は実際にどのように確認されますか？

一般的なアプローチはガウス-クロンロッドペアであり、これはガウス規則に追加のノードを加えて高次の推定値を生成します。2つの推定値の差は、適応積分器によって使用される実用的な誤差推定値として機能します。