主題
Search

高斯積分

DOWNLOAD Mathematica Notebook下載 Wolfram Notebook

高斯積分，也稱為機率積分，與 erf 函式密切相關，是一維高斯函式在 (-infty,infty) 上的積分。它可以使用組合兩個一維高斯函式的技巧來計算

int_(-infty)^inftye^(-x^2)dx=sqrt((int_(-infty)^inftye^(-x^2)dx)(int_(-infty)^inftye^(-x^2)dx))
(1)
=sqrt((int_(-infty)^inftye^(-y^2)dy)(int_(-infty)^inftye^(-x^2)dx))
(2)
=sqrt(int_(-infty)^inftyint_(-infty)^inftye^(-(x^2+y^2))dydx).
(3)

這裡，利用了積分中的變數是一個啞變數的事實，它在最後被積分掉，因此可以從 x 重新命名為 y。切換到 極座標 後得到

int_(-infty)^inftye^(-x^2)dx=sqrt(int_0^(2pi)int_0^inftye^(-r^2)rdrdtheta)
(4)
=sqrt(2pi[-1/2e^(-r^2)]_0^infty)
(5)
=sqrt(pi).
(6)

也存在一個簡單的證明這個恆等式的方法，不需要轉換到 極座標 (Nicholas and Yates 1950)。

從 0 到有限上限 a 的積分可以用連分數表示

int_0^ae^(-t^2)dt=1/2sqrt(pi)erf(a)
(7)
=1/2sqrt(pi)-(e^(-a^2))/(2a+)1/(a+)2/(2a+)3/(a+)4/(2a+...),
(8)

其中 erfx 是 erf (誤差函式)，正如拉普拉斯首次提出的，雅可比證明的，以及拉馬努金重新發現的 (Watson 1928; Hardy 1999, pp. 8-9)。

一般形式的積分類 of the form

I_n(a)=int_0^inftye^(-ax^2)x^ndx
(9)

可以透過設定以下公式解析求解

x=a^(-1/2)y
(10)
dx=a^(-1/2)dy
(11)
y^2=ax^2.
(12)

然後

I_n(a)=a^(-1/2)int_0^inftye^(-y^2)(a^(-1/2)y)^ndy
(13)
=a^(-(n+1)/2)int_0^inftye^(-y^2)y^ndy.
(14)

對於 n=0，這只是通常的高斯積分，所以

I_0(a)=(sqrt(pi))/2a^(-1/2)=1/2sqrt(pi/a).
(15)

對於 n=1，被積函式可以透過積分法積分，

I_1(a)=a^(-1)int_0^inftye^(-y^2)ydy=a^(-1)[-1/2e^(-y^2)]_0^infty=1/2a^(-1).
(16)

要計算 I_n(a)，對於 n>1，使用恆等式

-partial/(partiala)I_(n-2)(a)=-partial/(partiala)int_0^inftye^(-ax^2)x^(n-2)dx
(17)
=-int_0^infty-x^2e^(-ax^2)x^(n-2)dx
(18)
=int_0^inftye^(-ax^2)x^ndx
(19)
=I_n(a).
(20)

對於 n=2s 偶數，

I_n(a)=(-partial/(partiala))I_(n-2)(a)
(21)
=(-partial/(partiala))^2I_(n-4)
(22)
=...=(-partial/(partiala))^(n/2)I_0(a)
(23)
=(partial^(n/2))/(partiala^(n/2))I_0(a)
(24)
=(sqrt(pi))/2(partial^(n/2))/(partiala^(n/2))a^(-1/2),
(25)

所以

int_0^inftyx^(2s)e^(-ax^2)dx=((s-1/2)!)/(2a^(s+1/2))
(26)
=((2s-1)!!)/(2^(s+1)a^s)sqrt(pi/a),
(27)

其中 n!! 是雙階乘。如果 n=2s+1 是 奇數，那麼

I_n(a)=(-partial/(partiala))I_(n-2)(a)
(28)
=(-partial/(partiala))^2I_(n-4)(a)
(29)
=...=(-partial/(partiala))^((n-1)/2)I_1(a)
(30)
=(partial^((n-1)/2))/(partiala^((n-1)/2))I_1(a)
(31)
=1/2(partial^((n-1)/2))/(partiala^((n-1)/2))a^(-1),
(32)

所以

int_0^inftyx^(2s+1)e^(-ax^2)dx=(s!)/(2a^(s+1)).
(33)

因此解是

int_0^inftye^(-ax^2)x^ndx={((n-1)!!)/(2^(n/2+1)a^(n/2))sqrt(pi/a) for n even; ([1/2(n-1)]!)/(2a^((n+1)/2)) for n odd.
(34)

因此，前幾個值是

I_0(a)=1/2sqrt(pi/a)
(35)
I_1(a)=1/(2a)
(36)
I_2(a)=1/(4a)sqrt(pi/a)
(37)
I_3(a)=1/(2a^2)
(38)
I_4(a)=3/(8a^2)sqrt(pi/a)
(39)
I_5(a)=1/(a^3)
(40)
I_6(a)=(15)/(16a^3)sqrt(pi/a).
(41)

一個相關的，通常有用的積分是

H_n(a)=1/(sqrt(pi))int_(-infty)^inftye^(-ax^2)x^ndx,
(42)

這可以簡單地由下式給出

H_n(a)={(2I_n(a))/(sqrt(pi)) for n even; 0 for n odd.
(43)

更一般的積分 x^ne^(-ax^2+bx) 具有以下閉合形式，

int_(-infty)^inftyx^ne^(-ax^2+bx)dx=i^(-n)a^(-(n+1)/2)sqrt(pi)e^(b^2/(4a))U(-1/2n;1/2;-b^2/4a)
(44)
=sqrt(pi/a)e^(b^2/(4a))sum_(k=0)^(|_n/2_|)(n!)/(k!(n-2k)!)((2b)^(n-2k))/((4a)^(n-k))
(45)
=sqrt(pi/a)e^(b^2/(4a))sum_(k=0)^(|_n/2_|)(n; 2k)(2k-1)!!(2a)^(k-n)b^(n-2k)
(46)

對於整數 n>0 (F. Pilolli, 私人通訊)。對於 (45) 和 (46)，a,b in C-{0} ( 穿孔平面 )，R[a]>0，且 (-1)!!=1。 這裡，U(a;b;x) 是第二類合流超幾何函式，(n; k) 是 二項式係數。

另請參閱

Erf, 高斯積分, 高斯函式, 萊布尼茨積分法則, 正態分佈

使用 探索

參考文獻

Guitton, E. "Démonstration de la formule." Nouv. Ann. Math. 65, 237-239, 1906.Hardy, G. H. Ramanujan: Twelve Lectures on Subjects Suggested by His Life and Work, 3rd ed. New York: Chelsea, 1999.Nicholas, C. B. and Yates, R. C. "The Probability Integral." Amer. Math. Monthly 57, 412-413, 1950.Papoulis, A. Probability, Random Variables, and Stochastic Processes, 2nd ed. New York: McGraw-Hill, pp. 147-148, 1984.Watson, G. N. "Theorems Stated by Ramanujan (IV): Theorems on Approximate Integration and Summation of Series." J. London Math. Soc. 3, 282-289, 1928.

在 上被引用

高斯積分

引用為

Weisstein, Eric W. "高斯積分。" 來自 Web 資源。 https://mathworld.tw/GaussianIntegral.html

主題分類