Γ 函數在实数域上的函數圖形 。
此條目介紹的是擴展自階乘
x ! 的數學函數
Γ 。关于統計學分布,请见「
伽瑪分布 」。
伽瑪函數 (英語:gamma function ;記作Γ )是數學 中將階乘函數 解析延拓 至實數域 或複數域 後,所得的函數。對於定義域 中的複數z ,伽瑪函數滿足遞迴關係 :
Γ ( z + 1 ) = z Γ ( z ) .
更進一步,對於正整數n ,其函數值滿足Γ ( n ) = ( n − 1 ) ! 。
形式上,伽瑪函數可定義為以下瑕積分 :
Γ ( z ) = ∫ 0 ∞ t z − 1 e − t d t .
不過,積分結果只有當z 的實部 ℜ ( z ) > 0 時收斂。因此,須再經由遞迴關係進行解析延拓,以求出其餘複數的函數值(除了零 與負整數 )。
Γ 函數可以通过欧拉 (Euler)第二类积分定義:
Γ ( z ) = ∫ 0 ∞ t z − 1 e − t d t
对复数 z ,我们要求R e ( z ) > 0 。
Γ 函數还可以通过对e − t 做泰勒展开 ,解析延拓 到整个复平面 :
Γ ( z ) = ∫ 1 ∞ t z − 1 e t d t + ∑ n = 0 ∞ ( − 1 ) n n ! 1 n + z
这样定义的Γ 函數在全平面除了z = 0 , − 1 , − 2 , … 以外的地方解析。
Γ 函數也可以用无穷乘积 的方式表示:
Γ ( z ) = 1 z ∏ n = 1 ∞ ( 1 + z n ) − 1 ( 1 + 1 n ) z
这说明Γ ( z ) 是亚纯函数,而1 Γ ( z ) 是全纯函数。
Γ函數本身可以被看作是一個下列插值問題的解:
『找到一個光滑曲線連接那些由 y = ( x − 1 ) ! 所給定的點( x , y ) ,並要求x 要為正整數』
由前幾個的階乘 清楚地表明這樣的曲線是可以被畫出來的,但是我們更希望有一個精確的公式去描述這個曲線,並讓階乘的操作不會依賴於x 值的大小。而最簡單的階乘公式 x ! = 1 × 2 × ⋯ × x 不能直接應用在x 值為分数 的時候,因為它被限定在x 值為正整數而已。相對而言,并不存在一個有限的關於加總、乘積、冪次、指數函數或是對數函數可以表達 x ! ,但是是有一個普遍的公式藉由微積分的積分與極限去表達階乘的,而 Γ函數就是那個公式。[ 1]
階乘有無限多種的連續擴張方式將定義域擴張到非整數:可以通過任何一組孤立點畫出無限多的曲線。Γ函數是實務上最好的一個選擇,因為是解析的 (除了非正整數點),而且它可以被定義成很多種等價形式。然而,它並不是唯一一個擴張階乘意義的解析函數,只要給予任何解析函數,其在正整數上為零,像是 k sin ( m π x ) ,會給出其他函數有著階乘性質。
Γ 函數可以用無窮乘積 表示:
Γ ( z ) = lim n → ∞ n ! n z ∏ k = 0 n ( z + k ) − 1
Γ ( z ) = e − γ z z ∏ n = 1 ∞ ( 1 + z n ) − 1 e z n
其中γ 是欧拉-马歇罗尼常数 。
1 = ∫ 0 ∞ x α − 1 λ α e − λ x Γ ( α ) d x
⟹ Γ ( α ) λ α = ∫ 0 ∞ x α − 1 e − λ x d x
Γ 函数的递推公式为:
Γ ( x + 1 ) = x Γ ( x ) ,
对于正整数 n ,有
Γ ( n + 1 ) = n ! ,
可以说Γ 函数是階乘 的推廣。
Γ ( n + 1 ) = ∫ 0 ∞ e − x x n + 1 − 1 d x = ∫ 0 ∞ e − x x n d x
我们用分部积分法 来计算这个积分:
∫ 0 ∞ e − x x n d x = [ − x n e x ] 0 ∞ + n ∫ 0 ∞ e − x x n − 1 d x
当x = 0 时,− 0 n e 0 = 0 1 = 0 。当x 趋于无穷大 时,根据洛必达法则 ,有:
lim x → ∞ − x n e x = lim x → ∞ − n ! ⋅ 0 e x = 0 .
因此第一项[ − x n e x ] 0 ∞ 变成了零,所以:
Γ ( n + 1 ) = n ∫ 0 ∞ x n − 1 e x d x
等式的右面正好是n Γ ( n ) , 因此,递推公式 为:
Γ ( n + 1 ) = n Γ ( n ) .
Γ ( z ) Γ ( 1 − z ) = π sin π z ( 0 < R e ( z ) < 1 ) .
由此可知当 z = 1 2 时,Γ ( 1 2 ) = π .
Γ ( z ) Γ ( z + 1 2 ) = 2 1 − 2 z π Γ ( 2 z ) 。
Γ ( z ) Γ ( z + 1 m ) Γ ( z + 2 m ) ⋯ Γ ( z + m − 1 m ) = ( 2 π ) m − 1 2 m 1 2 − m z Γ ( m z ) .
Γ ( n + 1 2 ) = ( 2 n ) ! π n ! 4 n .
Γ ( 1 / 6 ) = Γ ( 1 / 3 ) 2 / π ∗ 2 2 / 3 ∗ sin ( π / 3 ) .
Γ ( 5 / 6 ) = 1 / Γ ( 1 / 3 ) 2 ∗ π 3 ∗ 2 4 / 3 / 3 .
Γ ( 1 / 10 ) = Γ ( 1 / 5 ) ∗ Γ ( 2 / 5 ) / π ∗ 2 4 / 5 ∗ sin ( 2 ∗ π / 5 ) .
Γ ( 3 / 10 ) = Γ ( 1 / 5 ) / Γ ( 2 / 5 ) ∗ π / 2 3 / 5 / sin ( 3 ∗ π / 10 ) .
Γ ( 7 / 10 ) = Γ ( 2 / 5 ) / Γ ( 1 / 5 ) ∗ π ∗ 2 3 / 5 .
Γ ( 9 / 10 ) = 1 / ( Γ ( 1 / 5 ) ∗ Γ ( 2 / 5 ) ) ∗ π 3 / 2 4 / 5 / ( sin ( π / 10 ) ∗ sin ( 2 ∗ π / 5 ) ) .
[ 2]
此式可用來協助計算t分布 機率密度函數、卡方分布 機率密度函數、F分布 機率密度函數等的累計機率。
對任何實數α
lim n → ∞ Γ ( n + α ) Γ ( n ) n α = 1 , α ∈ 𝐑
Γ函數與斯特靈公式
x y 圖例 Γ ( 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 ( z e ) z (橘色),數字越大2 π z ( z e ) z , 會越趨近Γ ( z + 1 ) 。但2 π z ( z e ) z 會在負值則會因為出現虛數而無法使用。
斯特靈公式 能用以估計Γ ( z ) 函数的增長速度。公式為:
Γ ( z + 1 ) ∼ 2 π z ( z e ) z ,
其中e 約等於2.718281828459。
Γ ( − 3 2 ) = 4 3 π ≈ 2.363 271 801 207 Γ ( − 1 2 ) = − 2 π ≈ − 3.544 907 701 811 Γ ( 1 2 ) = π ≈ 1.772 453 850 906 Γ ( 1 ) = 0 ! = 1 Γ ( 3 2 ) = 1 2 π ≈ 0.886 226 925 453 Γ ( 2 ) = 1 ! = 1 Γ ( 5 2 ) = 3 4 π ≈ 1.329 340 388 179 Γ ( 3 ) = 2 ! = 2 Γ ( 7 2 ) = 15 8 π ≈ 3.323 350 970 448 Γ ( 4 ) = 3 ! = 6
连分数表示
伽马函数也可以在复数域表示为两个连分数 之和[ 3] :
Γ ( z ) = e − 1 2 + 0 − z + 1 z − 1 2 + 2 − z + 2 z − 2 2 + 4 − z + 3 z − 3 2 + 6 − z + 4 z − 4 2 + 8 − z + 5 z − 5 2 + 10 − z + ⋱ + e − 1 z + 0 − z + 0 z + 1 + 1 z + 2 − z + 1 z + 3 + 2 z + 4 − z + 2 z + 5 + 3 z + 6 − ⋱
Γ函數的微分
x y 圖例 Γ ( 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 ,有
d n d z n Γ ( z ) = ∫ 0 ∞ t z − 1 e − t ( ln t ) n d t
於是,對任何正整數 m
Γ ′ ( m + 1 ) = m ! ( − γ + ∑ k = 1 m 1 k )
其中γ是歐拉-馬歇羅尼常數 。
Γ ( x + i y ) = { ∫ 1 ∞ t x − 1 e t cos ( y ln t ) d t + ∑ k = 0 ∞ ( − 1 ) k k ! [ k + x ( k + x ) 2 + y 2 ] } + i { ∫ 1 ∞ t x − 1 e t sin ( y ln t ) d t − ∑ k = 0 ∞ ( − 1 ) k k ! [ y ( k + x ) 2 + y 2 ] }
Γ函數的絕對值函數圖形
注意到在Γ 函數的積分定義中若取z 為實部大於零之複數 、則積分存在,而且在右半複平面上定義一個全純函數 。利用函數方程
Γ ( z ) Γ ( 1 − z ) = π sin π z ( 0 < R e ( z ) < 1 )
並注意到函數sin ( π z ) 在整個複平面上有解析延拓,我們可以在R e ( z ) < 1 時設
Γ ( z ) = π Γ ( 1 − z ) sin π z
從而將Γ 函數延拓為整個複平面上的亞純函數 ,它在z = 0 , − 1 , − 2 , − 3 ⋯ 有單極點 ,留數為
R e s ( Γ , − n ) = ( − 1 ) n n ! .
許多程式語言或試算表軟體有提供Γ函数或對數的Γ函数,例如EXCEL。而對數的Γ函数還要再取一次自然指數才能獲得Γ函数值。例如在EXCEL中,可使用GAMMALN函数,再用EXP[GAMMALN(X)],即可求得任意實数的伽玛函数的值。
例如在EXCEL中:EXP[GAMMALN(4/3)]=0.89297951156925
而在沒有提供Γ函数的程式環境中,也能夠過泰勒級數或斯特靈公式等方式來近似,例如Robert H. Windschitl在2002年提出的方法,其在十進制可獲得有效數字八位數的精確度[ 4] ,已足以填滿單精度浮點數 的二進制有效數字24位:
Γ ( z ) ≈ 2 π z ( z e z sinh 1 z + 1 810 z 6 ) z