1 概要
ガウス求積法は、定積分を少数の標本点で高精度に近似する数値積分法である。被積分関数の値を、あらかじめ選ばれた節点で評価し、それらに重みを付けて和を取ることで積分値を推定する。
この方法の大きな利点は、同じ節点数で比較した場合に、比較的高い次数の多項式まで正確に扱える点にある。計算資源が限られる状況でも有効であり、解析解が得にくい問題に対して広く利用されてきた。
1.1 定積分の数値近似
定積分の数値近似では、連続的な積分を有限個の計算可能な演算に置き換える。区分求積法のように区間を細かく分ける方法もあれば、少数の点を精密に選んで近似する方法もある。ガウス求積法は後者に属し、関数の形に応じて効率のよい近似を与える。
1.2 ガウス求積法の基本的な考え方
ガウス求積法では、積分区間と重み関数に適した直交多項式を用い、その零点を節点として採用する。重みは節点に対応して決まり、適切に構成すると、限られた点数でも高い精度が得られる。
この考え方は、単に点を増やすよりも、点の配置を最適化することが重要であるという発想に基づく。とくに、多項式を対象としたときに非常に良好な性質を示す。
1.3 他の数値積分法との比較
一般的な台形則やシンプソン則は、区間分割に基づくため実装が容易である一方、同じ計算量で得られる精度はガウス求積法に及ばないことが多い。ガウス型の公式は前処理として節点と重みの計算を要するが、一度求めれば繰り返しの積分評価に向く。
ただし、被積分関数に急な変化や特異性がある場合には、別の手法のほうが扱いやすいこともある。したがって、用途に応じた使い分けが重要である。
2 理論的基礎
ガウス求積法の理論は、直交多項式の性質と密接に結び付いている。節点の選択、重みの算出、誤差の構造は、いずれもこの枠組みの中で理解できる。
2.1 直交多項式
直交多項式は、ある内積のもとで互いに直交する多項式列である。ガウス求積法では、積分区間と重み関数に対応する直交系が中心的役割を果たす。
2.1.1 重み関数
重み関数は、積分で各点の寄与をどのように評価するかを定める関数である。これにより、ルジャンドル型、エルミート型、ラゲール型など、異なるガウス求積法が区別される。
重み関数の選択は、扱う積分の自然な領域や関数の減衰特性に合致している必要がある。適合がよいほど、少ない節点で安定した近似が得られる。
2.1.2 零点と節点
ガウス求積法では、直交多項式の零点が節点として使われる。零点は区間内に分布し、公式全体の精度を決める核心部分となる。
この選択により、節点が偏りにくく、重みとの組み合わせで高次の多項式まで正確に積分できる。零点の数値計算は、実用上も重要な課題である。
2.2 最適性の性質
ガウス求積法は、与えられた節点数に対して最良の精度を実現するように構成される。この性質が、他の求積法と比べた際の大きな特徴である。
2.2.1 多項式の完全性
ある求積公式が特定次数以下の多項式をすべて厳密に積分できるとき、その範囲で完全性をもつという。ガウス求積法は、この完全性を最大限に高めるよう作られている。
そのため、多項式で近似できる関数に対しては、非常に効率よく高精度を達成しやすい。実際の計算では、この性質が高い信頼性につながる。
2.2.2 最大次数の精度
節点数が同じなら、ガウス求積法は可能な限り高い次数まで厳密さを保つ。これは、重み付き積分の条件を満たすよう節点と重みが最適化されているためである。
この最適性は、単純な補間や等間隔の分割では得にくい。少数点で高精度を要する場面では、特に価値が大きい。
2.3 誤差解析
求積公式の実用性は、近似値だけでなく誤差の見積もりによっても左右される。ガウス求積法では、理論誤差の形が比較的明瞭であり、解析的な評価が行いやすい。
2.3.1 剰余項
剰余項は、真の積分値と求積近似の差を表す。ガウス求積法では、この差が高次の導関数や直交多項式の性質に依存する形で表されることが多い。
関数が十分滑らかであれば、誤差は急速に小さくなる傾向がある。逆に、非滑らかな点がある場合には、理論通りの改善が得られにくい。
2.3.2 誤差の評価
誤差評価では、関数の滑らかさ、区間の長さ、重み関数との整合性が考慮される。実際の数値計算では、理論誤差だけでなく丸め誤差や節点計算の不安定性も問題となる。
そのため、誤差評価は単なる公式の導出にとどまらず、実装の信頼性を支える重要な要素となる。
3 主な種類
ガウス求積法には、重み関数や積分区間の違いに応じた複数の型がある。代表的なものは、いずれも特定の直交多項式に対応している。
3.1 ガウス・ルジャンドル求積
ガウス・ルジャンドル求積は、区間が有限で重み関数が1の積分に用いられる。最も基本的なガウス型公式として知られ、標準区間の積分を扱う際の代表例である。
多くの教科書や数値計算ソフトで標準的に実装されており、他の方法の基準としても参照される。
3.2 ガウス・エルミート求積
ガウス・エルミート求積は、無限区間でガウス形の重みをもつ積分に適する。エルミート多項式の零点を節点として用いる。
正規分布に関係する積分や、指数減衰を伴う関数の処理に適しているため、確率計算や物理分野で有用である。
3.3 ガウス・ラゲール求積
ガウス・ラゲール求積は、半無限区間で指数減衰を含む積分に向く。ラゲール多項式に基づき、正の半直線上の積分を効率的に近似できる。
減衰が自然に組み込まれているため、積分核に指数因子が現れる問題で扱いやすい。
3.4 その他のガウス型求積
ガウス求積法には、基本三種以外にも多様な派生形がある。区間の形や重みの特性に応じて、専用の公式が設計される。
3.4.1 ガウス・チェビシェフ求積
ガウス・チェビシェフ求積は、チェビシェフ多項式に関連する重みをもつ。解析的な性質がよく、特定の積分に対して簡潔な形を取ることがある。
高速フーリエ変換と結び付く扱いもあり、近似計算の文脈で便利である。
3.4.2 ガウス・ジャコビ求積
ガウス・ジャコビ求積は、区間端点付近に特定の重み特性を持つ積分に用いられる。ジャコビ多項式を基礎とし、一般化された重み関数を扱える。
ルジャンドル型やチェビシェフ型を含む広い枠組みとして理解でき、端点近傍の振る舞いを反映しやすい。
4 計算法と応用
実際の利用では、理論だけでなく節点と重みをどう求めるかが重要である。加えて、対象積分を適切な形へ変換する工夫が精度と効率を左右する。
4.1 節点と重みの計算方法
節点と重みの計算は、ガウス求積法の実装上の中心課題である。小規模な場合は既知の表を使えるが、一般には数値的に求める必要がある。
4.1.1 直交多項式の零点の利用
節点は直交多項式の零点として与えられるため、まずその零点を数値的に求める。反復法や漸化式を利用して、高精度の近似値を得る方法が用いられる。
零点が精密に計算できれば、対応する重みも安定して決定できる。これは高次の公式ほど重要になる。
4.1.2 固有値問題による計算法
直交多項式の再帰関係は、三重対角行列の固有値問題として表せる。これを利用すると、節点を固有値、重みを固有ベクトルの成分から求められる。
この方法は数値的に安定で、大規模な節点数にも対応しやすい。現代の計算環境では広く用いられている。
4.2 変数変換と区間変換
標準形のガウス求積をそのまま使えない場合、変数変換によって積分区間や重みを整える。有限区間を標準区間に写す手法や、無限区間を扱いやすい形に変える手法が代表的である。
この前処理により、異なる問題を共通の求積公式で処理できる。変換の選び方は、数値安定性に直接影響する。
4.3 実際の応用分野
ガウス求積法は、解析解が困難な積分を含む多くの分野で使われる。高精度と計算効率の両立が求められる場面で特に有効である。
4.3.1 科学計算
科学計算では、微分方程式の離散化、量子力学、相関関数の評価などに現れる積分の近似に利用される。高精度が必要な反復計算との相性がよい。
4.3.2 工学解析
工学解析では、構造解析、熱伝導、電磁場計算などで積分評価が現れる。有限要素法や境界要素法の内部でも、要素積分の手段として重要である。
4.3.3 統計と確率計算
統計と確率計算では、期待値やモーメント、分布関数に関する積分の計算に用いられる。特に、重み付きの積分や正規分布に関連する式で有効である。
5 関連事項
ガウス求積法は数値積分全般の一部であり、求積公式や高次元化の議論とも接続している。周辺分野を理解すると、方法の位置付けがより明確になる。
5.1 数値積分全般
数値積分全般は、解析的に求めにくい積分を近似計算で解く分野である。区分的な近似、補間に基づく方法、適応的手法などが含まれる。
ガウス求積法は、その中でも少数点高精度型の代表例として位置付けられる。
5.2 求積公式
求積公式は、積分を有限和に置き換える一般的な式である。節点、重み、適用区間の組み合わせによって多様な形式が存在する。
ガウス求積法は、求積公式の中でも最適性が高い群に属する。
5.3 高次元積分への拡張
高次元積分では、単純に一変数の方法を繰り返すと計算量が急増する。そこで、テンソル積型の構成や、次元削減の工夫が必要になる。
ガウス求積法も高次元へ拡張できるが、実用上は次元の呪いが課題となる。そのため、問題構造に応じた変形や近似が併用されることが多い。