ガウス積分

函数 ƒ(x) = exp(−x²) のグラフおよび、函数と x-軸との間の面積 (= √π)

ガウス積分（がうす-せきぶん、""英: Gaussian integral""）あるいはオイラー-ポワソン積分（オイラーポワソンせきぶん、英: Euler–Poisson integral）はガウス関数 exp(−x²) の実数直線全体に亙る（両側無限）積分

$\int_{-\infty}^\infty e^{-x^2}dx = \sqrt{\pi}$

である。名称は、数学・物理学者のカール・フリードリヒ・ガウスに由来する。

この積分の応用は広い。例えば、変数の微小変化に伴う正規分布の正規化定数の計算に用いられる。積分の上の限界を有限な値に替えることで、誤差函数や正規分布の累積分布函数とも深く関連する。

誤差函数を表す初等函数は存在しないけれども、リッシュのアルゴリズムによって微分積分学の道具立てを用いてガウス積分の値が解析的に求まることが証明できる。つまり、$\textstyle\int e^{-x^2}\,dx$ に対する初等的な不定積分は存在しないけれど、定積分 $\textstyle\int_{-\infty}^\infty e^{-x^2}\,dx$ は評価することができるのである。

ガウス積分は物理学で非常に頻繁に現れ、またガウス積分の様々な一般化が場の量子論に現れる。

積分値の計算

極座標を用いて

ガウス積分を求める標準的な方法として、以下のアイデアはポワソンまで遡れる：

平面 R² 上の函数 exp(−x² − y²) = exp(−r²) を考え、これを二通りの方法で計算する。

1. 一つは直交座標系に関する二重積分として計算し、その値は求める値の平方になることを確かめる。
2. いま一つは極座標系に関する二重積分（いわゆるバウムクーヘン積分）として計算し、その値が π となることを確かめる。

広義積分が現れることに注意して、これら二つの計算を比較して積分の値が求まる。即ち、面積要素 dA が xy-直交座標系では dA = dx dy, rθ-極座標系では dA = r dr dθ で与えられることに注意すれば、

$\begin{align} \int_{\mathbb{R}^2} \exp(-(x^2+y^2))\,dA = \int_{-\infty}^\infty \int_{-\infty}^\infty e^{-(x^2+y^2)}\,dx\,dy = \left ( \int_{-\infty}^\infty \exp(-t^2)\,dt \right )^2, \end{align}$

および

$\begin{align} \int_{\mathbb{R}^2} \exp(-(x^2+y^2))\,dA &= \int_0^{2\pi} \int_0^{\infin} \exp(-r^2)\,r\,dr\,d\theta\\ &= 2\pi \int_0^\infty r\exp(-r^2)\,dr= \pi \int_{-\infty}^0 e^s\,ds = \pi \end{align}$

と計算できる。後者では s = −r² なる置換を行って、ds = −2r dr となることを用いている。さてこれらの結果から

$\left ( \int_{-\infty}^\infty \exp(-x^2)\,dx \right )^2=\pi$

であり、符号を考慮して

$\int_{-\infty}^\infty \exp(-x^2)\,dx=\sqrt{\pi}$

を得る。

上記の考察において、広義二重積分や二つの式を等しいとおいたことに対する正当性を再考しておこう。まずは近似函数

$I(a)=\int_{-a}^a \exp(-x^2)\,dx$

を考える。求めるガウス積分が絶対収斂ならば、それはコーシー主値、即ち

$\lim_{a\to\infty} I(a)$

なる極限によって求められることになる。これを見るには、

$\int_{-\infty}^\infty |\exp(-x^2)|\,dx < \int_{-\infty}^{-1} -x \exp(-x^2)\,dx + \int_{-1}^1 \exp(-x^2)\,dx+ \int_{1}^{\infty} x\exp(-x^2)\,dx<\infty$

が成り立つという事実を確かめればよい。故に I(a) の平方をとれば

$I(a)^2 = \int_{-a}^a \int_{-a}^a \exp(-(x^2+y^2))\,dx\,dy$

と書くことができて、フビニの定理により、これは xy-座標平面における面積分

$\int \exp(-(x^2+y^2))\,dA$

に等しいことが確かめられる。ただし、積分域は {(−a, a), (a, a), (a, −a), (−a, −a)} を頂点集合とする正方形である。

指数函数は全実数に対して正の値を取るから、上記の積分域の内接円上での積分は I(a)² よりも小さく、同様に外接円上での積分は I(a)² よりも大きい。これら二つの円板上での積分は、直交座標系から極座標系へ

$x = r \cos \theta,\quad y = r \sin\theta,\quad dx\,dy = r dr\,d\theta$

なる標準的な変換でうつれば容易に計算できるから、積分を実行して

$\pi (1-e^{-a^2}) < I(a)^2 < \pi (1 - e^{-2a^2})$

なる評価を得ることができる。a → ∞ なる極限をとれば、挟み撃ちの原理によって等式

$\left ( \int_{-\infty}^\infty \exp(-x^2)\,dx \right )^2=\pi$

が正当化できる。

直交座標を用いて

ガウス積分を計算する別な方法として、以下は Laplace (1812) にまで遡れる。

$y=xs,\quad dy=x\, ds$

と置くと、y を ±∞ へ近づけるとき s の極限は x の符号で決まるから、exp(−x²) が偶函数ゆえに実数全体に亙る積分が正の実数全体に亙る積分の2倍となること、つまり

$\int_{-\infty}^\infty e^{-x^2} \, dx=2\int_0^\infty e^{-x^2} \, dx$

であることを利用すれば計算が簡単になる。即ち、積分範囲を x ≥ 0 に限れば、変数 y と s とは同じ極限を持ち、

$I^2 =4\int_0^\infty \int_0^\infty e^{-(x^2 +y^2)} dy\, dx$

が成り立つ。故に

$\begin{align} \frac{I^2}{4} &=\int_0^\infty \left( \int_0^\infty e^{-(x^2 +y^2)} \, dy\right) \, dx=\int_0^\infty \left( \int_0^\infty e^{-x^2(1+s^2)} x\, ds\right) dx \\[5pt] &=\int_0^\infty \left( \int_0^\infty e^{-x^2(1+s^2 )} x\, dx\right) \, ds=\frac{1}{2} \int_0^\infty \frac{ds}{1+s^2} \\[5pt] &=\left. \frac{1}{2} \arctan s\, \right|_0^\infty =\frac{\pi}{4} \end{align}$

となり、所期の I = √π を得る。

ガンマ関数との関係

被積分関数が偶関数ゆえ

$\int_{-\infty}^\infty \exp (-x^2 )\, dx= 2\int_0^\infty \exp (-x^2 )\, dx$

が成り立ち、これに変数変換 x = t^1/2 を行えばオイラー積分

$2\int_0^\infty \exp(-x^2 )\, dx = 2\int_0^\infty \frac{1}{2} \exp(-t)\, t^{-1/2} \, dt=\Gamma \left( \frac{1}{2} \right) =\sqrt{\pi}$

が得られる。ここで γ はガンマ関数。この式は、半整数値の階乗が √π の有理数倍となる理由を示している。より一般に、

$\int_0^\infty \exp (-ax^b )\, dx = \frac{1}{b} a^{-\frac{1}{b}} \, \Gamma \left(\frac{1}{b} \right)$

が成り立つ。

一般化

ガウス関数の積分

詳細は「ガウス函数の積分」を参照

勝手なガウス函数の積分は

$\int_{-\infty}^{\infty}\exp\left(\frac{-(x+b)^2}{c^2}\right)dx = c\sqrt{\pi}$

あるいは

$\int_{-\infty}^{\infty}\exp(-x^2 + bx + c)dx = \sqrt{\pi}\,\exp\left(\frac{b^2}{4}+c\right)$

で与えられる。

多変数化

詳細は「多変数正規分布」を参照

A = (α_ij) が正定値対称（従って可逆な）共変行列（二階共変テンソル）ならば

$\int_{-\infty}^\infty \exp\left(-\frac{1}{2}\sum_{i,j=1}^{n}\alpha_{ij} x_i x_j \right)\,dx = \sqrt{\frac{(2\pi)^n}{\det A}}$

が成り立つ。ここで積分は Rⁿ 全体でとる。この事実は多変数正規分布の研究に応用される。また、

$\int x^{k_1}\cdots x^{k_{2N}}\,\exp\left(-\frac{1}{2}\sum_{i,j=1}^n \alpha_{ij} x^i x^j\right)dx = \sqrt{\frac{(2\pi)^n}{\det A}}\;\frac{1}{2^N N!}\sum_{\sigma \in S_{2N}}(A^{-1})^{k_{\sigma(1)}k_{\sigma(2)}}\cdots (A^{-1})^{k_{\sigma(2N-1)}k_{\sigma(2N)}}$

が成立する。ここで、σ は {1, ..., 2N} の置換であり、右辺に現れる余分な因子は A⁻¹ のN 個のコピーを {1, ..., 2N} の組合せ対 (combinatorial pairing) の全体に亘って加えた和を意味する。

あるいはまた、A⁻ = (β_ij) として

$\int f(\vec{x})\exp(-\frac{1}{2}\sum_{i,j=1}^n \alpha_{ij}x^i x^j)dx = \sqrt{(2\pi)^n\over \det A}\;\left.\exp\left({1\over 2}\sum_{i,j=1}^n \beta_{ij}{\partial \over \partial x^i}{\partial \over \partial x^j}\right)f(\vec{x})\right|_{\vec{x}=0}$

がいくつかの解析関数 f に対して成立する。f は増加具合が適当に制限されているとかあるいはほかの技術的な判定条件を満足する必要がある。これは特定の関数に対してはうまく行くがそうでないものもある。たとえば多項式ならば成立する。また微分作用素変数の指数関数 exp は冪級数として理解され、あらたな微分作用素を定めるものである。

さらに無限次元への一般化としての汎関数積分には厳密な定義は無く、多くの場合それは計算的でさえないが、ガウス汎関数積分を有限次元の場合の類似物として「定義」することができる。もちろん問題はあって、単純に有限次元の場合の式を無限次元の場合に適用しようとすれば (2π)^∞ は無限大に発散してしまうし、汎函数行列式 (functional determinant) も一般には無限大となりうる。これらのことを考慮して比

$\begin{align} &\frac{\int f(x_1)\cdots f(x_{2N}) \exp(-\frac{1}{2}\iint A(x_{2N+1},x_{2N+2}) f(x_{2N+1}) f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2})\,\mathcal{D}f}{\int \exp(-\frac{1}{2}\iint A(x_{2N+1},x_{2N+2}) f(x_{2N+1}) f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2})\,\mathcal{D}f}\\[5pt] &\qquad =\frac{1}{2^N N!}\sum_{\sigma \in S_{2N}}A^{-1}(x_{\sigma(1)},x_{\sigma(2)})\cdots A^{-1}(x_{\sigma(2N-1)},x_{\sigma(2N)})\end{align}$

のみを考えることにするならばガウス汎関数積分を扱うことができるという意味である。ドヴィット記法 (deWitt notation) を使えば、この等式は有限次元の場合と同じ形に書くことができる。

一次の項を持つ多変数ガウス積分

A = (α_ij) をやはり正定値対称行列として

$\int \exp\left(-\frac{1}{2}\sum_{i,j=1}^{n} \alpha_{ij} x_i x_j+\sum_{i=1}^{n} b_i x_i\right)d^nx = \sqrt{\frac{(2\pi)^n}{\det{A}}}\exp\left(\frac{1}{2}{}^t\mathbf{b}A^{-1}\mathbf{b}\right)$

が成り立つ。ただし、b = (b_i) で ^t は行列の転置とする。

被積分函数の多項式倍

同様の積分として、

$\begin{align} \int_0^\infty x^{2n} \exp\left(-\frac{x^2}{a^2}\right)dx &= \sqrt{\pi} \frac{(2n-1)!!}{2^{n+1}} a^{2n+1} =\sqrt{\pi}\frac{\left(2n\right)!}{n!}\left(\frac{a}{2}\right)^{2n+1},\\ \int_0^\infty x^{2n+1} \exp\left(-\frac{x^2}{a^2}\right)dx &= \frac{n!}{2} a^{2n+2} \end{align}$

が成立する。これらを導出するには積分記号下での微分法を用いるのが簡便である：

$\begin{align} & \int_{-\infty}^\infty x^{2n} \exp(-\alpha x^2)\,dx = \left(-1\right)^n\int_{-\infty}^\infty \frac{\partial^n}{\partial \alpha^n} e^{-\alpha x^2}\,dx = \left(-1\right)^n\frac{\partial^n}{\partial \alpha^n} \int_{-\infty}^\infty \exp(-\alpha x^2)\,dx \\[8pt] &= \sqrt{\pi} (-1)^n\frac{\partial^n}{\partial \alpha^n}\alpha^{-1/2} = \sqrt{\frac{\pi}{\alpha}}\frac{(2n-1)!!}{\left(2\alpha\right)^n}. \end{align}$

冪指数が高階多項式の場合

被積分函数の冪指数がもっと別の偶数次多項式に変わった場合も、級数解は容易に計算することができる。例えば四次多項式を冪指数とする指数函数の積分は

$\begin{align} &\int_{-\infty}^{\infty} \exp(\alpha x^4 + \beta x^3 + \gamma x^2 + \delta x + \epsilon)\,dx\\[5pt] &\qquad = \frac{1}{2}\,\exp(\epsilon)\sum_{n,m,p=0 \atop n+p=0 \bmod 2}^{\infty} \frac{\beta^n}{n!}\;\frac{\gamma^m}{m!}\;\frac{\delta^p}{p!}\; \frac{\Gamma(\frac{3n+2m+p+1}4)}{\exp(\frac{3n+2m+p+1}4\,\log(-\alpha))} \end{align}$

で表される。n + p = 0 mod 2 であることが要求されるのは、−∞ から 0 までの積分が各項に (−1)^n+p/2 なる因子として寄与し、0 から +∞ までの積分が各項に 1/2 の因子として寄与することによる。これらの積分は場の量子論に属する話題である。

関連項目

• 正規分布
• 誤差函数

参考文献

• Weisstein, Eric W., "Gaussian Integral" - MathWorld.（英語）
• David Griffiths. Introduction to Quantum Mechanics. 2nd Edition back cover.
• Abramowitz, M. and Stegun, I. A. Handbook of Mathematical Functions, Dover Publications, Inc. New York