跳转到内容

Γ函数

本页使用了标题或全文手工转换
维基百科,自由的百科全书
(重定向自Gamma function)
Γ函數在实数域上的函數圖形。

伽瑪函數(英語:gamma function;記作Γ)是數學中將階乘函數解析延拓至實數域或複數域後,所得的函數。對於定義域中的複數z,伽瑪函數滿足遞迴關係: Γ(z+1)=zΓ(z). 更進一步,對於正整數n,其函數值滿足Γ(n)=(n−1)!。

形式上,伽瑪函數可定義為以下瑕積分: Γ(z)=∫0∞tz−1e−tdt. 不過,積分結果只有當z的實部ℜ(z)>0時收斂。因此,須再經由遞迴關係進行解析延拓,以求出其餘複數的函數值(除了零與負整數)。

定義

[编辑]

Γ函數可以通过欧拉(Euler)第二类积分定義:

Γ(z)=∫0∞tz−1e−tdt

对复数z,我们要求Re(z)>0。

Γ函數还可以通过对e−t做泰勒展开,解析延拓到整个复平面: Γ(z)=∫1∞tz−1etdt+∑n=0∞(−1)nn!1n+z

这样定义的Γ函數在全平面除了z=0,−1,−2,…以外的地方解析。

Γ函數也可以用无穷乘积的方式表示:

Γ(z)=1z∏n=1∞(1+zn)−1(1+1n)z

这说明Γ(z)是亚纯函数,而1Γ(z)是全纯函数。

历史動機

[编辑]

Γ函數本身可以被看作是一個下列插值問題的解:

『找到一個光滑曲線連接那些由 y=(x−1)! 所給定的點(x,y),並要求x要為正整數』

由前幾個的階乘清楚地表明這樣的曲線是可以被畫出來的,但是我們更希望有一個精確的公式去描述這個曲線,並讓階乘的操作不會依賴於x值的大小。而最簡單的階乘公式 x!=1×2×⋯×x 不能直接應用在x值為分数的時候,因為它被限定在x值為正整數而已。相對而言,并不存在一個有限的關於加總、乘積、冪次、指數函數或是對數函數可以表達 x!,但是是有一個普遍的公式藉由微積分的積分與極限去表達階乘的,而 Γ函數就是那個公式。[1]

階乘有無限多種的連續擴張方式將定義域擴張到非整數:可以通過任何一組孤立點畫出無限多的曲線。Γ函數是實務上最好的一個選擇,因為是解析的(除了非正整數點),而且它可以被定義成很多種等價形式。然而,它並不是唯一一個擴張階乘意義的解析函數,只要給予任何解析函數,其在正整數上為零,像是 ksin⁡(mπx),會給出其他函數有著階乘性質。

無窮乘積

[编辑]

Γ函數可以用無窮乘積表示:

Γ(z)=limn→∞n!nz∏k=0n(z+k)−1
Γ(z)=e−γzz∏n=1∞(1+zn)−1ezn

其中γ是欧拉-马歇罗尼常数。

Γ積分

[编辑]
1=∫0∞xα−1λαe−λxΓ(α)dx

⟹Γ(α)λα=∫0∞xα−1e−λxdx

递推公式

[编辑]

Γ函数的递推公式为: Γ(x+1)=xΓ(x),

对于正整数n,有

Γ(n+1)=n!,

可以说Γ函数是階乘的推廣。

递推公式的推导

[编辑]

Γ(n+1)=∫0∞e−xxn+1−1dx=∫0∞e−xxndx

我们用分部积分法来计算这个积分:

∫0∞e−xxndx=[−xnex]0∞+n∫0∞e−xxn−1dx

当x=0时,−0ne0=01=0。当x趋于无穷大时,根据洛必达法则,有:

limx→∞−xnex=limx→∞−n!⋅0ex=0.

因此第一项[−xnex]0∞变成了零,所以:

Γ(n+1)=n∫0∞xn−1exdx

等式的右面正好是nΓ(n), 因此,递推公式为:

Γ(n+1)=nΓ(n).

重要性质

[编辑]
Γ(z)Γ(1−z)=πsin⁡πz(0<Re(z)<1).
由此可知当 z=12时,Γ(12)=π.
  • 乘法定理:
Γ(z)Γ(z+12)=21−2zπΓ(2z)。
Γ(z)Γ(z+1m)Γ(z+2m)⋯Γ(z+m−1m)=(2π)m−12m12−mzΓ(mz).
  • 此外:
Γ(n+12)=(2n)!πn!4n.
  • 使用乘法定理推導的關係:
Γ(1/6)=Γ(1/3)2/π∗22/3∗sin⁡(π/3).
Γ(5/6)=1/Γ(1/3)2∗π3∗24/3/3.
Γ(1/10)=Γ(1/5)∗Γ(2/5)/π∗24/5∗sin⁡(2∗π/5).
Γ(3/10)=Γ(1/5)/Γ(2/5)∗π/23/5/sin⁡(3∗π/10).
Γ(7/10)=Γ(2/5)/Γ(1/5)∗π∗23/5.
Γ(9/10)=1/(Γ(1/5)∗Γ(2/5))∗π3/24/5/(sin⁡(π/10)∗sin⁡(2∗π/5)).

[2]

此式可用來協助計算t分布機率密度函數、卡方分布機率密度函數、F分布機率密度函數等的累計機率。

  • 極限性質

對任何實數α

limn→∞Γ(n+α)Γ(n)nα=1,α∈𝐑

斯特靈公式

[编辑]
Γ函數與斯特靈公式
xy圖例Γ(z+1) (-0.81): 4.85Γ(z+1) (-0.8): 4.56Γ(z+1) (-0.79): 4.27Γ(z+1) (-0.77): 3.99Γ(z+1) (-0.75): 3.7Γ(z+1) (-0.74): 3.42Γ(z+1) (-0.71): 3.13Γ(z+1) (-0.69): 2.86Γ(z+1) (-0.65): 2.58Γ(z+1) (-0.62): 2.31Γ(z+1) (-0.57): 2.05Γ(z+1) (-0.51): 1.81Γ(z+1) (-0.44): 1.59Γ(z+1) (-0.36): 1.4Γ(z+1) (-0.27): 1.25Γ(z+1) (-0.17): 1.13Γ(z+1) (-0.07): 1.05Γ(z+1) (0.03): 0.98Γ(z+1) (0.14): 0.94Γ(z+1) (0.24): 0.91Γ(z+1) (0.35): 0.89Γ(z+1) (0.46): 0.89Γ(z+1) (0.56): 0.89Γ(z+1) (0.67): 0.9Γ(z+1) (0.77): 0.92Γ(z+1) (0.88): 0.96Γ(z+1) (0.99): 0.99Γ(z+1) (1.09): 1.04Γ(z+1) (1.19): 1.1Γ(z+1) (1.3): 1.17Γ(z+1) (1.4): 1.24Γ(z+1) (1.5): 1.33Γ(z+1) (1.6): 1.43Γ(z+1) (1.7): 1.54Γ(z+1) (1.8): 1.67Γ(z+1) (1.89): 1.81Γ(z+1) (1.98): 1.96Γ(z+1) (2.07): 2.13Γ(z+1) (2.15): 2.31Γ(z+1) (2.23): 2.5Γ(z+1) (2.31): 2.7Γ(z+1) (2.38): 2.91Γ(z+1) (2.45): 3.13Γ(z+1) (2.51): 3.37Γ(z+1) (2.57): 3.6Γ(z+1) (2.63): 3.85Γ(z+1) (2.68): 4.1Γ(z+1) (2.74): 4.35Γ(z+1) (2.78): 4.61Γ(z+1) (2.83): 4.87(2πx)(x/e)x (0.02): 0.33(2πx)(x/e)x (0.06): 0.49(2πx)(x/e)x (0.12): 0.59(2πx)(x/e)x (0.18): 0.66(2πx)(x/e)x (0.25): 0.69(2πx)(x/e)x (0.32): 0.72(2πx)(x/e)x (0.39): 0.74(2πx)(x/e)x (0.46): 0.75(2πx)(x/e)x (0.54): 0.77(2πx)(x/e)x (0.61): 0.79(2πx)(x/e)x (0.68): 0.8(2πx)(x/e)x (0.75): 0.83(2πx)(x/e)x (0.82): 0.85(2πx)(x/e)x (0.89): 0.87(2πx)(x/e)x (0.96): 0.9(2πx)(x/e)x (1.03): 0.93(2πx)(x/e)x (1.1): 0.97(2πx)(x/e)x (1.16): 1.01(2πx)(x/e)x (1.23): 1.05(2πx)(x/e)x (1.3): 1.1(2πx)(x/e)x (1.37): 1.15(2πx)(x/e)x (1.44): 1.2(2πx)(x/e)x (1.51): 1.26(2πx)(x/e)x (1.57): 1.33(2πx)(x/e)x (1.64): 1.4(2πx)(x/e)x (1.7): 1.47(2πx)(x/e)x (1.77): 1.56(2πx)(x/e)x (1.83): 1.64(2πx)(x/e)x (1.89): 1.74(2πx)(x/e)x (1.95): 1.84(2πx)(x/e)x (2.01): 1.94(2πx)(x/e)x (2.07): 2.05(2πx)(x/e)x (2.13): 2.17(2πx)(x/e)x (2.18): 2.29(2πx)(x/e)x (2.23): 2.42(2πx)(x/e)x (2.29): 2.55(2πx)(x/e)x (2.33): 2.69(2πx)(x/e)x (2.38): 2.83(2πx)(x/e)x (2.43): 2.97(2πx)(x/e)x (2.47): 3.12(2πx)(x/e)x (2.52): 3.27(2πx)(x/e)x (2.56): 3.43(2πx)(x/e)x (2.6): 3.59(2πx)(x/e)x (2.64): 3.75(2πx)(x/e)x (2.67): 3.92(2πx)(x/e)x (2.71): 4.08(2πx)(x/e)x (2.74): 4.25(2πx)(x/e)x (2.78): 4.42(2πx)(x/e)x (2.81): 4.6(2πx)(x/e)x (2.84): 4.77
Γ(z+1)(藍色)、2πz(ze)z(橘色),數字越大2πz(ze)z,會越趨近Γ(z+1)。但2πz(ze)z會在負值則會因為出現虛數而無法使用。

斯特靈公式能用以估計Γ(z)函数的增長速度。公式為:

Γ(z+1)∼2πz(ze)z,

其中e約等於2.718281828459。

特殊值

[编辑]
Γ(−32)=43π≈2.363271801207Γ(−12)=−2π≈−3.544907701811Γ(12)=π≈1.772453850906Γ(1)=0!=1Γ(32)=12π≈0.886226925453Γ(2)=1!=1Γ(52)=34π≈1.329340388179Γ(3)=2!=2Γ(72)=158π≈3.323350970448Γ(4)=3!=6

连分数表示

伽马函数也可以在复数域表示为两个连分数之和[3]:

Γ(z)=e−12+0−z+1z−12+2−z+2z−22+4−z+3z−32+6−z+4z−42+8−z+5z−52+10−z+⋱+e−1z+0−z+0z+1+1z+2−z+1z+3+2z+4−z+2z+5+3z+6−⋱

导数

[编辑]
Γ函數的微分
xy圖例Γ(x) (-1.98): 30.23Γ(x) (-1.98): 22.79Γ(x) (-1.97): 18.29Γ(x) (-1.95): 10.96Γ(x) (-1.87): 4.31Γ(x) (-1.54): 2.32Γ(x) (-1.21): 4.67Γ(x) (-1.09): 11.24Γ(x) (-1.05): 18.12Γ(x) (-1.04): 24.83Γ(x) (-1.03): 31.37Γ(x) (-1.03): 37.17Γ(x) (-1.02): 45.56Γ(x) (-1.01): 50Γ(x) (-1): 50Γ(x) (-1): 50Γ(x) (-1): 50Γ(x) (-1): -20Γ(x) (-1): -20Γ(x) (-1): -20Γ(x) (-1): -20Γ(x) (-0.99): -20Γ(x) (-0.99): -20Γ(x) (-0.93): -15.12Γ(x) (-0.87): -8.21Γ(x) (-0.61): -3.75Γ(x) (-0.27): -4.68Γ(x) (-0.1): -10.72Γ(x) (-0.06): -17.65Γ(x) (0.02): 50Γ(x) (0.03): 39.18Γ(x) (0.03): 31.1Γ(x) (0.04): 23.82Γ(x) (0.06): 17.15Γ(x) (0.09): 10.11Γ(x) (0.24): 3.8Γ(x) (0.57): 1.55Γ(x) (0.93): 1.05Γ(x) (1.29): 0.9Γ(x) (1.65): 0.9Γ(x) (2): 1Γ(x) (2.36): 1.21Γ(x) (2.72): 1.57Γ(x) (3.07): 2.14Γ(x) (3.43): 3.08Γ(x) (3.78): 4.58Γ(x) (4.12): 6.95Γ(x) (4.43): 10.49Γ(x) (4.69): 15.25Γ(x) (4.91): 20.93(Γ(x))′ (-1.85): -19.81(Γ(x))′ (-1.82): -14.08(Γ(x))′ (-1.78): -8.45(Γ(x))′ (-1.68): -3.05(Γ(x))′ (-1.51): 1.47(Γ(x))′ (-1.36): 6.28(Γ(x))′ (-1.27): 11.75(Γ(x))′ (-1.23): 17.42(Γ(x))′ (-1.2): 23.12(Γ(x))′ (-1.18): 28.71(Γ(x))′ (-1.17): 34.37(Γ(x))′ (-1.16): 40.05(Γ(x))′ (-1.14): 46.48(Γ(x))′ (-1.05): 50(Γ(x))′ (-0.86): 50(Γ(x))′ (-0.85): 41.95(Γ(x))′ (-0.84): 36.42(Γ(x))′ (-0.82): 30.68(Γ(x))′ (-0.81): 24.91(Γ(x))′ (-0.78): 19.27(Γ(x))′ (-0.75): 13.56(Γ(x))′ (-0.69): 7.98(Γ(x))′ (-0.58): 2.67(Γ(x))′ (-0.44): -2.23(Γ(x))′ (-0.32): -7.48(Γ(x))′ (-0.26): -13.04(Γ(x))′ (-0.22): -18.79(Γ(x))′ (0.03): -20(Γ(x))′ (0.23): -17.94(Γ(x))′ (0.28): -12.32(Γ(x))′ (0.37): -6.9(Γ(x))′ (0.56): -2.74(Γ(x))′ (0.83): -1.03(Γ(x))′ (1.11): -0.39(Γ(x))′ (1.4): -0.06(Γ(x))′ (1.68): 0.18(Γ(x))′ (1.97): 0.4(Γ(x))′ (2.26): 0.66(Γ(x))′ (2.54): 0.99(Γ(x))′ (2.83): 1.46(Γ(x))′ (3.11): 2.15(Γ(x))′ (3.39): 3.16(Γ(x))′ (3.67): 4.66(Γ(x))′ (3.93): 6.84(Γ(x))′ (4.18): 9.84(Γ(x))′ (4.39): 13.65(Γ(x))′ (4.57): 18.09(Γ(x))′ (4.72): 22.97(Γ(x))′ (4.85): 28.12(Γ(x))′ (4.95): 33.44
Γ函數(藍色)、Γ函數的微分(橘色),其中,大於50與小於-20的部分被截掉。

對任何複數z,滿足 Re(z) > 0,有

dndznΓ(z)=∫0∞tz−1e−t(ln⁡t)ndt

於是,對任何正整數 m

Γ′(m+1)=m!(−γ+∑k=1m1k)

其中γ是歐拉-馬歇羅尼常數。

复数值

[编辑]
Γ(x+iy)={∫1∞tx−1etcos⁡(yln⁡t)dt+∑k=0∞(−1)kk![k+x(k+x)2+y2]}+i{∫1∞tx−1etsin⁡(yln⁡t)dt−∑k=0∞(−1)kk![y(k+x)2+y2]}

解析延拓

[编辑]
Γ函數的絕對值函數圖形

注意到在Γ函數的積分定義中若取z為實部大於零之複數、則積分存在,而且在右半複平面上定義一個全純函數。利用函數方程

Γ(z)Γ(1−z)=πsin⁡πz(0<Re(z)<1)

並注意到函數sin⁡(πz)在整個複平面上有解析延拓,我們可以在Re(z)<1時設

Γ(z)=πΓ(1−z)sin⁡πz

從而將Γ函數延拓為整個複平面上的亞純函數,它在z=0,−1,−2,−3⋯有單極點,留數為

Res(Γ,−n)=(−1)nn!.

程式實現

[编辑]

許多程式語言或試算表軟體有提供Γ函数或對數的Γ函数,例如EXCEL。而對數的Γ函数還要再取一次自然指數才能獲得Γ函数值。例如在EXCEL中,可使用GAMMALN函数,再用EXP[GAMMALN(X)],即可求得任意實数的伽玛函数的值。

  • 例如在EXCEL中:EXP[GAMMALN(4/3)]=0.89297951156925

而在沒有提供Γ函数的程式環境中,也能夠過泰勒級數或斯特靈公式等方式來近似,例如Robert H. Windschitl在2002年提出的方法,其在十進制可獲得有效數字八位數的精確度[4],已足以填滿單精度浮點數的二進制有效數字24位:

Γ(z)≈2πz(zezsinh⁡1z+1810z6)z

参见

[编辑]

參考文獻

[编辑]
  1. ^ P. J., Davis. Leonhard Euler's Integral: A Historical Profile of the Gamma Function. American Mathematical Monthly. 1959 [2023-01-01]. doi:10.2307/2309786. (原始内容存档于2023-01-01). 
  2. ^ Mada, L. Relations of the Gamma function. R code on Github. Code publicly available on Github [Personal Research]. 2020-04-24 [2020-04-24]. (原始内容存档于2021-04-02). Relations of the Gamma function 
  3. ^ Exponential integral E: Continued fraction representations. [2023-01-01]. (原始内容存档于2022-11-09). 
  4. ^ Viktor T. Toth. "Programmable Calculators: Calculators and the Gamma Function". 2006 [2018-11-18]. (原始内容存档于2007-02-23). 

外部链接

[编辑]