主题
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, 高斯积分, 高斯函数, 莱布尼茨积分法则, 正态分布

使用 Wolfram|Alpha 探索

参考文献

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.

在 Wolfram|Alpha 上被引用

高斯积分

引用为

Weisstein, Eric W. "高斯积分。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/GaussianIntegral.html

主题分类