主題
Search

Erf

DOWNLOAD Mathematica Notebook下載 Wolfram 筆記本
Erf

erf(z) 是在積分正態分佈（它是高斯函式的歸一化形式）中遇到的“誤差函式”。它是由以下定義的整函式

erf(z)=2/(sqrt(pi))int_0^ze^(-t^2)dt.
(1)

請注意，一些作者（例如，Whittaker 和 Watson 1990，第 341 頁）將 erf(z) 定義為沒有前導因子 2/sqrt(pi)。

Erf 在 Wolfram 語言 中實現為Erf[z]。 給出 erf(z_1)-erf(z_0) 的雙引數形式也實現為Erf[z0, z1]。

Erf 滿足以下恆等式

erf(z)=1-erfc(z)
(2)
=(2z)/(sqrt(pi))_1F_1(1/2;3/2;-z^2)
(3)
=(2ze^(-z^2))/(sqrt(pi))_1F_1(1;3/2;z^2),
(4)

其中 erfc(z) 是 erfc，互補誤差函式，而 _1F_1(a;b;z) 是第一類合流超幾何函式。對於 z>0，

erf(z)=pi^(-1/2)gamma(1/2,z^2),
(5)

其中 gamma(a,x) 是不完全伽瑪函式。

Erf 也可以定義為麥克勞林級數

erf(z)=2/(sqrt(pi))sum_(n=0)^(infty)((-1)^nz^(2n+1))/(n!(2n+1))
(6)
=2/(sqrt(pi))(z-1/3z^3+1/(10)z^5-1/(42)z^7+1/(216)z^9+...)
(7)

(OEIS A007680)。 同樣地，

erf^2(z)=4/pi(z^2-2/3z^4+(14)/(45)z^6-4/(35)z^8+(166)/(4725)z^(10)+...)
(8)

(OEIS A103979 和 A103980)。

對於 x<<1，erf(x) 可以從以下公式計算

erf(x)=1/(sqrt(pi))e^(-x^2)sum_(n=0)^(infty)((2x)^(2n+1))/((2n+1)!!)
(9)
=2/(sqrt(pi))e^(-x^2)[x+(2x^3)/(1·3)+(4x^5)/(1·3·5)+...]
(10)

(OEIS A000079 和 A001147；Acton 1990)。

對於 x>>1，

erf(x)=2/(sqrt(pi))(int_0^inftye^(-t^2)dt-int_x^inftye^(-t^2)dt)
(11)
=1-2/(sqrt(pi))int_x^inftye^(-t^2)dt.
(12)

使用分部積分得到

int_x^inftye^(-t^2)dt=-1/2int_x^infty1/td(e^(-t^2))
(13)
=-1/2[(e^(-t^2))/t]_x^infty-1/2int_x^infty(e^(-t^2)dt)/(t^2)
(14)
=(e^(-x^2))/(2x)+1/4int_x^infty1/(t^3)d(e^(-t^2))
(15)
=(e^(-x^2))/(2x)-(e^(-x^2))/(4x^3)-...,
(16)

因此

erf(x)=1-(e^(-x^2))/(sqrt(pi)x)(1-1/(2x^2)-...)
(17)

繼續此過程得到漸近級數

erf(x)∼1-(e^(-x^2))/(sqrt(pi))sum_(n=0)^(infty)((-1)^n(2n-1)!!)/(2^n)x^(-(2n+1))
(18)
∼1-(e^(-x^2))/(sqrt(pi))(x^(-1)-1/2x^(-3)+3/4x^(-5)-(15)/8x^(-7)
(19)
+(105)/(16)x^(-9)+...)
(20)

(OEIS A001147 和 A000079)。

Erf 的值如下

erf(0)=0
(21)
erf(infty)=1.
(22)

它是奇函式

erf(-z)=-erf(z),
(23)

並滿足

erf(z)+erfc(z)=1.
(24)

Erf 可以用第一類合流超幾何函式 M 表示為

erf(z)=(2z)/(sqrt(pi))M(1/2,3/2,-z^2)
(25)
=(2z)/(sqrt(pi))e^(-z^2)M(1,3/2,z^2).
(26)

它的導數是

(d^n)/(dz^n)erf(z)=(-1)^(n-1)2/(sqrt(pi))H_(n-1)(z)e^(-z^2),
(27)

其中 H_n 是埃爾米特多項式。一階導數是

d/(dz)erf(z)=2/(sqrt(pi))e^(-z^2),
(28)

積分是

interf(z)dz=zerf(z)+(e^(-z^2))/(sqrt(pi)).
(29)
ErfReImAbs
最小值 最大值
實部
虛部 Powered by webMathematica

Erf 也可以擴充套件到複平面，如上圖所示。

一個簡單的涉及 erf 的積分，Wolfram 語言 無法計算，由下式給出

int_0^pe^(-x^2)erf(p-x)dx=1/2sqrt(pi)[erf(1/2sqrt(2)p)]^2
(30)

(M. R. D'Orsogna，私人通訊，2004 年 5 月 9 日)。更復雜的積分包括

int_0^infty(e^(-(p+x)y))/(pi(p+x))sin(asqrt(x))dx=-sinh(asqrt(p)) +(e^(-asqrt(p)))/2erf(a/(2sqrt(y))-sqrt(py))+(e^(asqrt(p)))/2erf(a/(2sqrt(y))+sqrt(py)) int_0^infty(sqrt(x)e^(-(p+x)y))/(pi(p+x))cos(asqrt(x))dx=(e^(-[py+a^2/(4y)]))/(sqrt(piy))+sqrt(p)[-cosh(asqrt(p))-(e^(-asqrt(p)))/2erf(a/(2sqrt(y))-sqrt(py))+(e^(asqrt(p)))/2erf(a/(2sqrt(y))+sqrt(py))]
(31)

(M. R. D'Orsogna，私人通訊，2005 年 12 月 15 日)。

Erf 有連分數

int_0^xe^(-t^2)dt=1/2sqrt(pi)erf(x)
(32)
=1/2sqrt(pi)-(1/2e^(-x^2))/(x+1/(2x+2/(x+3/(2x+4/(x+...)))))
(33)

(Wall 1948，第 357 頁)，最早由拉普拉斯於 1805 年和勒讓德於 1826 年提出 (Olds 1963，第 139 頁)，由雅可比證明，並由拉馬努金重新發現 (Watson 1928；Hardy 1999，第 8-9 頁)。

涉及 erf(x) 的定積分包括 涉及 erf(x) 的定積分包括

int_0^inftye^(-px^2)erf(ax)erf(bx)dx=1/(sqrt(pip))tan^(-1)((ab)/(sqrt(p(a^2+b^2+p))))
(34)
int_0^inftyxe^(-px^2)erf(ax)erf(bx)erf(cx)dx=1/(pip)[a/(sqrt(a^2+p))tan^(-1)((bc)/(sqrt((a^2+b^2+c^2+p)(a^2+p))))+b/(sqrt(b^2+p))tan^(-1)((ac)/(sqrt((a^2+b^2+c^2+p)(b^2+p))))+c/(sqrt(c^2+p))tan^(-1)((ab)/(sqrt((a^2+b^2+c^2+p)(c^2+p))))]
(35)
int_0^inftye^(-x)erf(sqrt(x))dx=1/2sqrt(2)
(36)
int_0^inftye^(-x)erf^2(sqrt(x))dx=(2sqrt(2)cot^(-1)(sqrt(2)))/pi
(37)
int_0^inftye^(-x)erf^3(sqrt(x))dx=(3sqrt(2)cot^(-1)(2sqrt(2)))/pi.
(38)

其中前兩個出現在 Prudnikov et al. （1990，第 123 頁，方程式 2.8.19.8 和 2.8.19.11）中，其中 R[p]>0，|arg(a)|,|argb|,|argc|<pi/4。

複數對 erf(x) 的推廣定義為

w(z)=e^(-z^2)erfc(-iz)
(39)
=e^(-z^2)(1+(2i)/(sqrt(pi))int_0^ze^(t^2)dt).
(40)

僅在上半平面 I[z]>0 中有效的積分表示由下式給出

w(z)=i/piint_(-infty)^infty(e^(-t^2))/(z-t)dt
(41)
=(2iz)/piint_0^infty(e^(-t^2))/(z^2-t^2)dt.
(42)

另請參閱

道森積分, Erfc, Erfi, 菲涅爾積分, 高斯函式, 高斯積分, 反誤差函式, 正態分佈函式, 歐文 T 函式, 機率積分 在 課堂中探索此主題

相關的 Wolfram 網站

http://functions.wolfram.com/GammaBetaErf/Erf/, http://functions.wolfram.com/GammaBetaErf/Erf2/

使用 探索

參考文獻

Abramowitz, M. and Stegun, I. A. (Eds.). "誤差函式和菲涅爾積分。" 第 7 章，數學函式手冊，包含公式、圖表和數學表格，第 9 次印刷。 New York: Dover, pp. 297-309, 1972.Acton, F. S. 有效的數值方法，第 2 次印刷。 Washington, DC: Math. Assoc. Amer., p. 16, 1990.Arfken, G. 物理學家的數學方法，第 3 版。 Orlando, FL: Academic Press, pp. 568-569, 1985.Hardy, G. H. 拉馬努金：關於他的生活和工作提出的主題的十二次講座，第 3 版。 New York: Chelsea, 1999.Havil, J. 伽瑪：探索尤拉常數。 Princeton, NJ: Princeton University Press, p. 105, 2003.Olds, C. D. 連分數。 New York: Random House, 1963.Prudnikov, A. P.; Brychkov, Yu. A.; and Marichev, O. I. 積分與級數，第 2 卷：特殊函式。 New York: Gordon and Breach, 1990.Sloane, N. J. A. Sequences A000079/M1129, A001147/M3002, A007680/M2861, A103979, A103980 in "整數序列線上百科全書."Spanier, J. and Oldham, K. B. "誤差函式 erf(x) 及其補函式 erfc(x)。" 第 40 章，函式圖集。 Washington, DC: Hemisphere, pp. 385-393, 1987.Wall, H. S. 連分數解析理論。 New York: Chelsea, 1948.Watson, G. N. "拉馬努金提出的定理（四）：關於近似積分和級數求和的定理。" J. London Math. Soc. 3, 282-289, 1928.Whittaker, E. T. and Robinson, G. "誤差函式。" §92 in 觀測演算：數值數學專著，第 4 版。 New York: Dover, pp. 179-182, 1967.Whittaker, E. T. and Watson, G. N. 現代分析教程，第 4 版。 Cambridge, England: Cambridge University Press, 1990.

在 中被引用

Erf

引用為

Weisstein, Eric W. "Erf." 來自 網路資源。 https://mathworld.tw/Erf.html

主題分類