不完全伽玛函数

在數學中，**上不完全 Γ 函數**和**下不完全 Γ 函數**是 $\Gamma $ 函數的推廣。它們的定義分別如下：

　
: $\Gamma (s,x)=\int _{x}^{\infty }t^{s-1}\,e^{-t}\,{\rm {d}}t.\,\!\qquad \gamma (s,x)=\int _{0}^{x}t^{s-1}\,e^{-t}\,{\rm {d}}t.\,\quad \Re (s)>0,x\in \mathbb {R} _{0}^{+}$

通過解析延拓可以將定義域拓展到 **C**×**C** （除去可數個奇點外），詳見下文。

## 記號

如無特別說明，在本文中，以 x 表示非負實數，以 z 表示任意複數。

## 基本性質

在上面定義的約束條件下，通過分部積分，可以計算得遞歸關係：

　
: $\Gamma (s,x)=(s-1)\Gamma (s-1,x)+x^{s-1}e^{-x}$

以及反過來：

　
: $\gamma (s,x)=(s-1)\gamma (s-1,x)-x^{s-1}e^{-x}$

因為正常的 $\Gamma $ 函數定義為：

　
: $\Gamma (s)=\int _{0}^{\infty }t^{s-1}\,e^{-t}\,{\rm {d}}t$

故有：

　
: $\Gamma (s)=\Gamma (s,0)$

以及

　
: $\gamma (s,x)+\Gamma (s,x)=\Gamma (s).$

## 解析延拓

### 下不完全伽瑪函數

#### 解析延拓的方法

在最原始的定義式中，積分是沿著實軸進行的，故要求在 γ(s,z) 中，

　
: $z\in \mathbb {R} _{0}^{+},\Re (s)>0$

運用上一小節裡面的關係式，可以用下式來進行解析延拓：

　
: $\gamma (s,z)=\sum _{k=0}^{\infty }{\frac {z^{s}e^{-z}z^{k}}{s(s+1)\dots (s+k)}}=z^{s}\,\Gamma (s)\,e^{-z}\sum _{k=0}^{\infty }{\frac {z^{k}}{\Gamma (s+k+1)}}$

由比值審斂法可知，右邊的級數的收斂半徑是無窮大。

由魏爾斯特拉斯原理，下式中的函數，有時記作 $\gamma ^{*}$ ，是關於 s 和 z 的整函數。

　
: $\gamma ^{*}(s,z):=e^{-z}\sum _{k=0}^{\infty }{\frac {z^{k}}{\Gamma (s+k+1)}}$

因此下面的分解式

　
: $\gamma (s,z)=z^{s}\,\Gamma (s)\,\gamma ^{*}(s,z)$

的確給出了下不完全伽瑪函數的一個解析延拓。其中前兩個因子給出了下不完全伽瑪函數的奇點（即 z=0 或 s 為非正整數），而後面的因子則給出了下不完全伽瑪函數的零點。

#### 多值性

下不完全伽瑪函數的多值性來自於因子 zs 的多值性。如無特別說明，本文限制 z 的輻角絕對值小於 π。

#### 積分表達式

在選定了 zs 的單值分支之後，下不完全伽瑪函數的積分定義式可以自然地拓展到 z 為任意複數的情形，只是此時該積分應該理解為複平面上的路徑積分，且積分路徑需避開單值分支間的割線。需注意的是此時仍然要求 s 的實部大於 0，否則積分不收斂。

#### z→∞ 時的極限

s 為正實數的情形，有定義式有：

　
: $\lim _{z\rightarrow +\infty }\gamma (s,z)=\Gamma (s)$

s 為複數且不為非正整數的情形，可以證明：

　
: $\lim _{|z|\rightarrow +\infty }\gamma (s,z)=\Gamma (s),\quad |\arg z|<{\frac {\pi }{2}}-\varepsilon $

後面的條件相當於要求 z 的實部為正值且輻角取主值。

#### 總結

根據上面的討論，下不完全伽瑪函數有下列性質：

* 當 s 為正整數時是 z 的整函數；
* 當 s 不是整數時是 z 的多值全純函數，z=0 是其枝點；
* 對於確定的 z，指定主分支後，下不完全伽瑪函數是 s 的亞純函數，非正整數是其一階極點。

### 上不完全伽瑪函數

#### 解析延拓的方法

當 z 為正實數，s 為實部大於 0 的複數時，有定義顯然有：

　
: $\Gamma (s,z)=\Gamma (s)-\gamma (s,z)$

由於伽瑪函數和下不完全伽瑪函數關於 s,z 都至少是亞純函數，上式可以自然地作解析延拓，並以此作為上不完全伽瑪函數的定義。下不完全伽瑪函數的多值性自然地導致上不完全伽瑪函數的多值性，下面的討論基於主分支。

進一步地，由黎曼可去奇點原理，由於等號右邊在 s 取非正整數時的鄰域內有界，故作為 s 的函數，非正整數是上不完全伽瑪函數的可去奇點，可以通過對等號右邊取極限來定義非正整數時上不完全伽瑪函數的值。下面以 s=0 為例來說明這種極限過程，其它情況可以類推得到。

事實上，在下不完全伽瑪函數的積分表達式中，將指數函數用其泰勒展開式代換，得到：

　
: $\gamma (s,x)=\int _{0}^{x}t^{s-1}e^{-t}\operatorname {d} t=\int _{0}^{x}\sum _{k=0}^{\infty }(-1)^{k}\,{\frac {t^{s+k-1}}{k!}}\operatorname {d} t=\sum _{k=0}^{\infty }(-1)^{k}\,{\frac {x^{s+k}}{k!(s+k)}}=x^{s}\,\sum _{k=0}^{\infty }{\frac {(-x)^{k}}{k!(s+k)}}$

即，

　
: $\gamma ^{*}(s,x)=\sum _{k=0}^{\infty }{\frac {(-x)^{k}}{k!\,\Gamma (s)(s+k)}}.$

上式實際上給出了 γ(s,z) 的一個級數表示，給定 s 後，由比值審斂法知上式的收斂半徑為無窮大。下面的討論將 x 換成 z（z≠0）。

　
: $\gamma (s,z)-{\frac {1}{s}}=-{\frac {1}{s}}+z^{s}\,\sum _{k=0}^{\infty }{\frac {(-z)^{k}}{k!(s+k)}}={\frac {z^{s}-1}{s}}+z^{s}\,\sum _{k=1}^{\infty }{\frac {(-z)^{k}}{k!(s+k)}},\quad \Re (s)>-1,\,s\neq 0$

上式等號右邊第二項當 s=0 時有良好的定義，第一項在 s→0 時的極限存在，故等號右邊於 s→0 時的極限存在，並可以用它來定義等號左邊的值。

　
: $\lim _{s\rightarrow 0}\gamma (s,z)-{\frac {1}{s}}=\ln z+\sum _{k=1}^{\infty }{\frac {(-z)^{k}}{k\cdot k!}}$

另一方面，由伽瑪函數的魏爾斯特拉斯無窮乘積表示有：

　
: $\lim _{s\rightarrow 0}\Gamma (s)-{\frac {1}{s}}=-\gamma $

γ 是歐拉 - 馬歇羅尼常數。於是，可以定義

　
: $\Gamma (0,z)=\lim _{s\rightarrow 0}\left[\Gamma (s)-{\frac {1}{s}}-\left(\gamma (s,z)-{\frac {1}{s}}\right)\right]=-\gamma -\ln z-\sum _{k=1}^{\infty }{\frac {(-z)^{k}}{k\,(k!)}}$

這樣就證明了 s=0 的確是上不完全伽瑪函數的可去奇點。

#### 總結

上不完全伽瑪函數的其它解析性質可以由下不完全伽瑪函數和（完全）伽瑪函數的解析性質得到。結果如下：

* 當 s 是正整數時，是 z 的整函數；
* 當 s 不是整數時，是 z 的多值全純函數，z=0 是其枝點；
* 選定單值分支後，對 z≠0，是 s 的整函數；
* 當 s 的實部大於零且 z=0 時，等於（完全）伽瑪函數 Γ(s)；

注意最後一條對一般的 s 並不成立。特別地，當 s 為負實數且不為整數時，Γ(s) 是實數，而 Γ(s,0) 沒有定義。

## 特殊值

　
: $\Gamma (s,0)=\Gamma (s),\quad \Re (s)>0$

下面一組關係式都能夠由積分表達式直接得出，其中第三、四式中的函數是指數積分，第五、六式中的函數分別是余誤差函數和誤差函數：

　
: $\Gamma (1,z)=e^{-z}$

　
: $\gamma (1,z)=1-e^{-z}$

　
: $\Gamma (0,x)=-{\rm {Ei}}(-x),\quad x>0,$

　
: $\Gamma (s,z)=z^{s}\,{\rm {E}}_{1-s}(z),\quad s\in \mathbb {R} ^{+}$

　
: $\Gamma \left({\tfrac {1}{2}},x\right)={\sqrt {\pi }}\,{\rm {erfc}}\left({\sqrt {x}}\right),$

　
: $\gamma \left({\tfrac {1}{2}},x\right)={\sqrt {\pi }}\,{\rm {erf}}\left({\sqrt {x}}\right).$

## 與合流超幾何函數的關係

上不完全伽瑪函數和下不完全伽瑪函數都可以用合流超幾何函數表示，詳見合流超幾何函數一文。

## 導函數

由不完全伽瑪函數的積分表達式顯然有：

　
: ${\frac {\partial \gamma (s,z)}{\partial z}}=-{\frac {\partial \Gamma (s,z)}{\partial z}}=z^{s-1}e^{-z}$

儘管積分表達式為保證收斂性要加上 s 的實部大於零的限制，但上式並沒有這種限制。這可以通過 s 的實部大於零時的對應表達式兩邊作解析延拓證明。

另一方面，不完全伽瑪函數對 s 的偏導數是 Meijer G- 函數的特例，事實上，定義

　
: $T(m,s,z)=G_{m-1,\,m}^{\,m,\,0}\!\left(\left.{\begin{matrix}0,0,\dots ,0\\s-1,-1,\dots ,-1\end{matrix}}\;\right|\,x\right)$

則有：

　
: ${\frac {\partial \Gamma (s,z)}{\partial s}}=\Gamma (s,z)\ln z+z\,T(3,s,z)$

　
: ${\frac {\partial ^{2}\Gamma (s,z)}{\partial s^{2}}}=\Gamma (s,z)\ln ^{2}z+2z[\ln z\,T(3,s,z)+T(4,s,z)]$

更一般地：

　
: ${\frac {\partial ^{m}\Gamma (s,z)}{\partial s^{m}}}=\Gamma (s,z)\ln ^{m}z+mz\,\sum _{n=0}^{m-1}(m-1)_{n}\ln ^{m-n-1}z\,T(3+n,s,z)$

式中 (m-1)n 是下降階乘冪的 Pochhammer 記號。

事實上

　
: $\int _{x}^{\infty }{\frac {t^{s-1}\ln ^{m}t}{e^{t}}}{\rm {d}}t={\frac {\partial ^{m}}{\partial s^{m}}}\int _{x}^{\infty }{\frac {t^{s-1}}{e^{t}}}{\rm {d}}t={\frac {\partial ^{m}}{\partial s^{m}}}\Gamma (s,x)$

利用上式和 Mellin 變換的性質，並作解析延拓，就可以得到上不完全伽瑪函數對參變量的高階偏導數的表達式。

## 注

## 參考資料

* Incomplete Gamma function. [2013-03-10]. （原始內容存檔於 2009-09-02）.  §6.5.
* Allasia, Giampietro; Besenghi, Renata. Numerical calculation of incomplete gamma functions by the trapezoidal rule. Numer. Math. 1986, **50** (4): 419–428. doi:10.1007/BF0139666.
* Amore, Paolo. Asymptotic and exact series representations for the incomplete Gamma function. Europhys. Lett. 2005, **71** (1): 1–7. MR 2170316.
* G. Arfken and H. Weber. _Mathematical Methods for Physicists_. Harcourt/Academic Press, 2000. _(See Chapter 10.)_
* DiDonato, Armido R.; Morris, Jr., Alfred H. Computation of the incomplete gamma function ratios and their inverse. ACM Transactions on Mathematical Software (TOMS). 1986 年 12 月， **12** (4): 377–393. doi:10.1145/22721.23109.
* Barakat, Richard. Evaluation of the Incomplete Gamma Function of Imaginary Argument by Chebyshev Polynomials. Math. Comp. 1961, **15** (73): 7–11. MR 0128058.
* Carsky, Petr; Polasek, Martin. Incomplete Gamma F\_m(x) functions for real and complex arguments. J. Comput. Phys. 1998, **143** (1): 259–265. MR 1624704. doi:10.1006/jcph.1998.5975.
* Chaudhry, M. Aslam; Zubair, S. M. On the decomposition of generalized incomplete Gamma functions with applications to Fourier transforms. J. Comput. Appl. Math. 1995, **59** (101): 253–284. MR 1346414.
* DiDonato, Armido R.; Morris, Jr., Alfred H. ALGORITHM 654: FORTRAN subroutines for computing the incomplete gamma function ratios and their inverse. ACM Transactions on Mathematical Software (TOMS). 1987 年 9 月， **13** (3): 318–319. doi:10.1145/29380.214348.
* Früchtl, H.; Otto, P. A new algorithm for the evaluation of the incomplete Gamma Function on vector computers. ACM Trans. Math. Softw. 1994, **20** (4): 436–446.
* Gautschi, Walter. The incomplete gamma function since Tricomi. Atti Convegni Lincei. 1998, **147**: 203–237. MR 1737497.
* Gautschi, Walter. A Note on the recursive calculation of Incomplete Gamma Functions. ACM Trans. Math. Softw. 1999, **25** (1): 101–107. MR 1697463.
* Gradshteyn, I.S.; Ryzhik, I.M. Tables of Integrals, Series, and Products 4th. New York: Academic Press. 1980. ISBN 0-12-294760-6.  _（見 8.35 章）_
* Jones, William B.; Thron, W. J. On the computation of incomplete gamma functions in the complex domain **12–13**: 401–417. 1985. MR 0793971.
* Hazewinkel, Michiel (編), Incomplete gamma-function, 数学百科全书，Springer, 2001, ISBN 978-1-55608-010-4
* Mathar, Richard J. Numerical representation of the incomplete gamma function of complex-valued argument. Numerical Algorithms. 2004, **36** (3): 247–264. MR 2091195. doi:10.1023/B:NUMA.0000040063.91709.5.
* Miller, Allen R.; Moskowitz, Ira S. On certain Generalized incomplete Gamma functions. J. Comput. Appl. Math. 1998, **91** (2): 179–190.
* Paris, R. B., 不完全 Γ 函數，Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (編), NIST Handbook of Mathematical Functions, Cambridge University Press, 2010, ISBN 978-0521192255, MR2723248
* Paris, R. B. A uniform asymptotic expansion for the incomplete gamma function. J. Comput. Appl. Math. 2002, **148** (2): 323–339. MR 1936142. doi:10.1016/S0377-0427(02)00553-8.
* Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP. Section 6.2. Incomplete Gamma Function and Error Function. Numerical Recipes: The Art of Scientific Computing 3rd. New York: Cambridge University Press. 2007 [2013-03-10]. ISBN978-0-521-88068-8. （原始內容存檔於 2021-04-15）.
* Takenaga, Roy. On the Evaluation of the Incomplete Gamma Function. Math. Comp. 1966, **20** (96): 606–610. MR 0203911.
* Temme, Nico. Uniform Asymptotic Expansions of the Incomplete Gamma Functions and the Incomplete Beta Function. Math. Comp. 1975, **29** (132): 1109–1114. MR 0387674.
* Terras, Riho. The determination of incomplete Gamma Functions through analytic integration. J. Comp. Phys. 1979, **31**: 146–151. MR 0531128.
* Tricomi, Francesco G. Sulla funzione gamma incompleta. Ann. Mat. Pura Applic. 1950, **31**: 263–279. MR 0047834. doi:10.1007/BF02428264.
* Tricomi, F. G. Asympotitsche Eigenschaften der unvollst. Gammafunktion. Math. Zeitsch. 1950, **53** (2): 136–148. MR 0045253.
* van Deun, Joris; Cools, Ronald. A stable recurrence for the incomplete gamma function with imaginary second argument. Numer. Math. 2006, **104**: 445–456. MR 2249673. doi:10.1007/s00211-006-0026-1.
* Winitzki, Serge. Computing the incomplete gamma function to arbitrary precision. Lect. Not. Comp. Sci. 2003, **2667**: 790–798. MR 2110953.
* 埃里克・韋斯坦因. Incomplete Gamma Function. MathWorld.

## 參見

* Gamma 函數

## 外部連結

* http://functions.wolfram.com/GammaBetaErf/Gamma2/ （頁面存檔備份，存於網際網路檔案館） wolfram 上關於上不完全 $\Gamma $ 函數的一些恆等式