跳轉到

9.1 GAMMA函數#

本章導言#

階乘只定義在整數上；Gamma 函數把它推廣到複數——複分析中最重要的特殊函數。Euler 1729/1730 的積分、Weierstrass 的乘積、函數方程 \(\Gamma(z+1) = z\Gamma(z)\) 、餘元公式 \(\Gamma(z)\Gamma(1-z) = \frac{\pi}{\sin \pi z}\) 、Stirling 漸近—— \(\Gamma\) 是「一個函數、五副面孔」的教科書典範，也是 Bohr–Mollerup 定理（凸性唯一性）的主角。本節講述其五種表示、關鍵恆等式、漸近分析，與 Legendre、Gauss、Weierstrass 的接力史。

定義與五副面孔#

一、Euler 積分（1729/1738 出版）：

\[ \Gamma(z) = \int_0^\infty e^{-t} t^{z-1}\,dt, \qquad \operatorname{Re} z > 0. \]

二、函數方程：分部積分

\[ \Gamma(z+1) = \int_0^\infty e^{-t}t^{z}\,dt = \left[-e^{-t}t^z\right]_0^\infty + z\int_0^\infty e^{-t}t^{z-1}dt = z\,\Gamma(z). \]

配合 \(\Gamma(1) = 1\) ：** \(\Gamma(n + 1) = n!\) ——階乘的複變推廣。注意**： \(\Gamma\) 在 \(z\) 處的值對應 \((z-1)!\) （平移一位的歷史約定）。

三、Gauss 乘積形式（1812）：

\[ \Gamma(z) = \lim_{n\to\infty}\frac{n!\,n^z}{z(z+1)\cdots(z+n)}. \]

四、Weierstrass 乘積（8.1 篂）：

\[ \frac{1}{\Gamma(z)} = z\,e^{\gamma z}\prod_{n=1}^\infty\left(1 + \frac zn\right)e^{-z/n}, \qquad \gamma = \lim\left(\sum_{k\le n}\frac1k - \log n\right) = 0.5772\ldots \]

五、Hankel 圍道（7.3 篂場景二）：

\[ \frac{1}{\Gamma(z)} = \frac{1}{2\pi i}\int_{\mathcal H} e^t t^{-z}\,dt \]

——繞負實軸的鑰匙孔圍道，全平面的定義，割線跳躍產生 \(\frac{1}{\Gamma}\) 的整函數性。

解析性質（由乘積面孔讀出）： \(\Gamma\) 半純，極點 \(z = 0, -1, -2, \ldots\) ，留數 \(\frac{(-1)^n}{n!}\) ；無零點（ \(1/\Gamma\) 整）。五副面孔的分工：積分管計算、方程管遞推、Gauss 管極限、乘積管延拓、Hankel 管圍道積分。

關鍵恆等式#

餘元公式（reflection formula，Euler 1749/Legendre 1809）：

\[ \Gamma(z)\,\Gamma(1 - z) = \frac{\pi}{\sin \pi z}. \]

證明思路（由 8.1 篂正弦乘積）： \(\Gamma(z)\Gamma(1-z) = \int\!\!\int e^{-(s+t)}s^{z-1}t^{-z}ds\,dt\) ，換元 \(u = s + t, v = \frac{s}{s+t}\) 化為 Beta 積分 \(B(z, 1-z) = \int_0^\infty \frac{v^{z-1}}{1+v}dv = \frac{\pi}{\sin \pi z}\) （後者由圍道積分/部分分數——8.2 篂技術）。餘元公式 = \(\zeta\) 函數方程的姊妹（8.3 篂函數方程中的 \(\sin\frac{\pi s}{2}\Gamma(1-s)\) 因子正是它）。

倍角公式（Legendre 1809）：

\[ \Gamma(z)\,\Gamma\left(z + \frac12\right) = 2^{1 - 2z}\sqrt\pi\,\Gamma(2z). \]

Beta–Gamma： \(B(x, y) = \int_0^1 t^{x-1}(1-t)^{y-1}dt = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}\) ——二項係數的複變推廣： \(\binom{\alpha}{n} = \frac{\Gamma(\alpha+1)}{n!\,\Gamma(\alpha - n + 1)}\) 對任意複 \(\alpha\) （4.1 篂二項級數的係數語言）。

特殊值： \(\Gamma\left(\frac12\right) = \sqrt\pi\) （Gauss 積分 \(\int e^{-t^2}dt\) ）、 \(\Gamma\left(\frac{n}{2}\right)\) 族（ \(n\) 維球體積公式 \(V_n = \frac{\pi^{n/2}r^n}{\Gamma(\frac n2 + 1)}\) ——** \(\Gamma\) 是高維幾何的計量器）、 \(\Gamma\left(\frac13\right) \approx 2.6789\) （無初等形式——超越性**：Chudnovsky 1976 證 \(\Gamma(\frac1q)\) 的代數無關結果族）。

Bohr–MOLLERUP 定理：唯一性的「性格」#

問題： \(\Gamma(z+1) = z\Gamma(z)\) 、 \(\Gamma(1) = 1\) 的解有無窮多（乘以 \(p(z)\) 型 1-週期函數： \(\Gamma(z)p(z)\) 同樣滿足）——什麼額外條件鎖定 \(\Gamma\) ？

Bohr–Mollerup 定理（1922）： \(f(x) = \Gamma(x)\) 是唯一在 \((0, \infty)\) 上滿足

  1. \(f(1) = 1\) ；
  2. \(f(x+1) = xf(x)\) ；
  3. \(\log f\) 凸

的函數。

證明思想：由 2 與 3， \(\log f\) 的差商單調： \(\frac{\log f(x) - \log f(n)}{x - n}\) 夾在整數點差商之間；取 \(n \to \infty\) （用 Gauss 極限形式）把 \(f\) 完全夾死。∎

意義： \(\Gamma\) 不是「隨便一個」階乘推廣——它是唯一「對數凸」的。這條定理使「 \(\Gamma\) 是階乘的正確推廣」從信仰變成定理。衍生應用：Binet 公式、Karamata 型不等式、機率（Gamma 分布的對數凸性）。Harald Bohr（物理學家 Niels Bohr 之弟、丹麥足球國腳——1908 奧運銀牌！）與數學家 Mollerup 的這條定理，是「函數的特徵化」方法的經典：用少數性質唯一鎖定對象。

STIRLING 漸近#

Stirling 公式（1730 實數版、複變推廣）：

\[ \Gamma(z) \sim \sqrt{2\pi}\,z^{z - \frac12}e^{-z}\left(1 + \frac{1}{12z} + \frac{1}{288z^2} - \cdots\right) \qquad (|\arg z| < \pi - \delta). \]

證明技術： \(\Gamma(z) = \int_0^\infty e^{-t}t^{z-1}dt\) 換元 \(t = zu\) ： \(= z^z e^{-z}\int_0^\infty e^{-z(u - \log u - 1)}du\) ——被積函數在 \(u = 1\) 有鞍點（ \(\phi(u) = u - \log u - 1\) 、 \(\phi'(1) = 0\) 、 \(\phi''(1) = 1 > 0\) ），最陡下降法（5.2 篂場景一）展開：

\[ \Gamma(z) \sim z^{z-\frac12}e^{-z}\sqrt{2\pi}\sum \frac{a_n}{z^n}. \]

注意：右側級數發散（ \(a_n \sim (-1)^n (2n)!/(…)\) 增長）——漸近級數（4.2 篂場景四）：截斷給出越來越好的近似（誤差 < 首個捨棄項）。Stirling 的階 \(\Gamma(z) \sim z^{z-\frac12}e^{-z}\) 給出： \(\Gamma\) 的階（order） \(\rho = 1\) （ \(1/\Gamma\) 的階 1）、** \(n! \sim \sqrt{2\pi n}(n/e)^n\) **——機率（Stirling 逼近）、組合學、資訊論（熵近似 \(\log n! = n\log n - n + O(\log n)\) ）的日常工具。

歷史：Abraham de Moivre 1730 已得 \(n! \sim C\sqrt n (n/e)^n\) （常數 \(C\) 未知）；Stirling 1730 確定 \(C = \sqrt{2\pi}\) （透過 Wallis 乘積！）—— \(\pi\) 又一次出現在階乘與機率的深處：常態分布的歸一化 \(\frac{1}{\sqrt{2\pi}}\) 與 Stirling 常數同源（1.3 篂場景四的機率連線閉環）。

一個具體例題#

例題 9.1.1：(a) 用餘元公式計算 \(\Gamma\left(-\frac12\right)\) ；(b) 用倍角公式計算 \(\Gamma\left(\frac{9}{2}\right)\) ；(c) 用 Stirling 估計 \(10!\) ，給出誤差界。

解：

(a) 餘元公式取 \(z = \frac12\) ： \(\Gamma\left(\frac12\right)^2 = \frac{\pi}{\sin\frac\pi2} = \pi\) ，故 \(\Gamma\left(\frac12\right) = \sqrt\pi\) 。再取 \(z = -\frac12\) ： \(\Gamma\left(-\frac12\right)\Gamma\left(\frac32\right) = \frac{\pi}{\sin(-\frac\pi2)} = -\pi\) 。而 \(\Gamma\left(\frac32\right) = \frac12\Gamma\left(\frac12\right) = \frac{\sqrt\pi}{2}\) 。故

\[ \Gamma\left(-\frac12\right) = \frac{-\pi}{\sqrt\pi/2} = -2\sqrt\pi. \]

核對函數方程： \(\Gamma\left(\frac12\right) = \left(-\frac12\right)\Gamma\left(-\frac12\right) = \left(-\frac12\right)(-2\sqrt\pi) = \sqrt\pi\) ✓。

(b) 反覆倍角或遞推： \(\Gamma\left(\frac92\right) = \frac72\cdot\frac52\cdot\frac32\cdot\frac12\,\Gamma\left(\frac12\right) = \frac{7\cdot5\cdot3\cdot1}{2^4}\sqrt\pi = \frac{105\sqrt\pi}{16}\) 。一般規律： \(\Gamma\left(n + \frac12\right) = \frac{(2n)!}{4^n\,n!}\sqrt\pi\) ——半整數階乘的封閉式（量子力學諧振子、氫原子的徑向波函數歸一化常數來源）。

(c) Stirling： \(10! \approx \sqrt{2\pi\cdot10}\left(\frac{10}{e}\right)^{10}\left(1 + \frac{1}{120}\right)\) 。數值： \(\sqrt{20\pi} \approx 7.9267\) ； \(\left(\frac{10}{e}\right)^{10} = 10^{10}e^{-10} = 10^{10} \times 4.53999\times10^{-5} = 453999\) ；乘積 \(= 7.9267 \times 453999 \approx 3598699\) ；修正 \(\times 1.00833\) ： \(\approx 3628669\) 。真值 \(10! = 3628800\) ；誤差 \(\approx 131\) （相對 \(3.6\times10^{-5}\) ）✓（理論誤差界：下一項 \(\frac{1}{288\cdot100} = 3.47\times10^{-5}\) 的量級——誤差落在預測區間）。兩項 Stirling 在 \(n = 10\) 已達萬分之三精度——漸近級數的實用威力。

歷史背景：三位主角與一段遺憾#

  • Euler（1729–1738）：21 歲的 Euler 在 Daniel Bernoulli 通信中提出積分定義（1729.10.13 致 Goldbach 的信），1738 發表。** \(\Gamma\) 的積分是 18 歲天才的第一擊**。
  • Legendre（1809–1811）：命名「gamma 函數」、倍角公式、廣義 Beta、以及（遺憾的）記號選擇： \(\Gamma(z)\) vs Gauss 的 \(\Pi(z) = \Gamma(z+1)\) 。Gauss 的記號（ \(\Pi(z) = z!\) 直覺）更自然但未流行——記號的偶然性影響兩百年的教學。
  • Gauss（1812）：《關於超幾何級數的一般研究》以無窮乘積極限定義（第三面孔），系統研究 \(\Gamma\) 的超越性質。
  • Weierstrass（1856）： \(1/\Gamma\) 乘積與 \(\gamma\) 常數的嚴格化。
  • Stirling（1730）/ de Moivre（1730–33）：漸近——機率論的需求（二項分布的近似）驅動。
  • Bohr–Mollerup（1922）：特徵化定理——二十世紀的收尾。

** \(\gamma\) （Euler–Mascheroni 常數）的補注： \(\gamma = 0.5772156649\ldots\) ——它是有理還是無理至今未知（與 RH 同級的未解問題之一，雖然形式更簡單）；已算到 \(10^{12}\) 位。「階乘的推廣」連帶產生的常數，其最基本問題仍未解**——數學的深淵常藏在最樸素的定義旁。

歷史啟示#

  • 推廣的正確性需要「特徵化定理」： \(\Gamma\) 的多種定義（積分/方程/乘積）+ 唯一性定理（Bohr–Mollerup）= 「正確推廣」的完整論證。新定義的合法性 = 等價性 + 特徵性。
  • 漸近比精確更常用： \(\Gamma\) 的「真值」少有人用，Stirling 前兩項每天都在用。分析學的實用層在漸近——精確定義存在是為了漸近展開的合法性。
  • 記號的歷史包袱： \(\Gamma(z) = (z-1)!\) 的平移錯位是歷史偶然；讀者應自行建立「 \(\Gamma(z+1) = z!\) 」的語感防呆。

與其他章節及領域的關聯#

  • 8.1 篂： \(1/\Gamma\) 乘積是 Weierstrass 定理的元應用。
  • 8.3 篂： \(\zeta\) 函數方程中的 \(\Gamma(1-s)\) 因子。
  • 4.1 篂：二項級數係數 = \(\Gamma\) 比值。
  • 5.2 篂：Hankel 圍道與最陡下降（Stirling 證明）。
  • 6.3 篂：SC 積分的 Beta 函數。
  • 9.2 篂：橢圓函數的母函數（ \(\Gamma\) 常數進入橢圓積分的週期）。
  • 機率論：Gamma 分布、Chi 方分布（ \(\Gamma(\frac n2)\) 歸一化）、Dirichlet 分布（多元 Beta）、Beta 分布——統計推斷的母函數。
  • 物理：諧振子、氫原子（半整數 \(\Gamma\) ）、路徑積分的測度（Gauss 積分 + Stirling）、弦論的 Veneziano 振幅 \(B\left(-\alpha(s), -\alpha(t)\right) = \frac{\Gamma(-\alpha s)\Gamma(-\alpha t)}{\Gamma(-\alpha s - \alpha t)}\) ——1968 年弦論的誕生公式就是 Beta 函數！
  • 數值：lgamma 函數（對數 Gamma，防溢出）是所有科學計算庫的標配。

練習#

  1. 由 Euler 積分證明 \(\Gamma(z + 1) = z\Gamma(z)\) ，並用極限過程說明 \(\operatorname{Re} z > 0\) 外如何由方程遞推延拓。
  2. 計算 \(\Gamma\left(\frac{5}{2}\right)\) 、 \(\Gamma\left(-\frac{3}{2}\right)\) ，並驗證它們與函數方程自洽。
  3. 用餘元公式證明 \(\int_0^\infty \frac{t^{x-1}}{1+t}dt = \frac{\pi}{\sin\pi x}\) （ \(0 < x < 1\) ）。（提示： \(t = \frac{u}{1-u}\) 化為 Beta。）
  4. 用 Stirling 前兩項估計 \(\binom{100}{50}\) （用 \(\Gamma\) 表示後兩次 Stirling），給出誤差量級。
  5. 思考題：Bohr–Mollerup 的「對數凸」條件為什麼是「正確的」條件？若改用「凸」或「解析」會發生什麼？（提示： \(\Gamma(z)p(z)\) 的 \(p\) 可選什麼？）

進階研讀： \(\Gamma\) 的五個深化場景#

場景一：複 Stirling 與 Stokes 線。複變 Stirling 在扇形 \(|\arg z| < \pi - \delta\) 有效；跨過 \(\arg z = \pm\pi\) 時漸近展開的指數項 \(\frac{1}{12z}\) 附近跳變（Stokes 現象，7.1 篂場景四）—— \(1/\Gamma\) 在負實軸的對數跳躍由 Hankel 圍道讀出： \(\log\Gamma(-x \pm i0) = \log\Gamma(x) \mp i\pi x\) 型。** \(\Gamma\) 是 Stokes 現象的最乾淨標本**。

場景二：Barnes G 函數與多重 \(\Gamma\) 。 \(G(z+1) = \Gamma(z)G(z)\) （「超階乘」）、Barnes 雙 \(\Gamma\) 、 \(q\) -Gamma（ \(q\) -階乘）： \(\Gamma\) 的譜系——每個「遞推方程 + 極點配置 + 凸性特徵」組合定義一族函數。現代應用：量子場論的等效作用量（Conformal Field Theory 的結構常數含 \(G\) ）、隨機矩陣的 Selberg 積分（Selberg 積分 = 多元 Beta，1944—— \(\Gamma\) 家族的高維集成，隨機矩陣譜的歸一化母體）。

場景三：WKB 與譜漸近。 \(n\) 維諧振子的能階 \(E_n = n + \frac12\) 中的 \(\frac12\) （零點能）來自 \(\Gamma\) 的半整數值；Stirling 給Weyl 律（高激發態能階密度 \(\sim \frac{1}{2\pi}\) 型）的 \(\Gamma\) 證明。「量子譜的漸近計數 = \(\Gamma\) + Stirling」——從 Bohr–Sommerfeld（1915）到 Gutzwiller 跡公式（1971）的譜漸近傳統。

場景四：Veneziano 振幅與弦論。1968 年 Veneziano 發現強子散射振幅 = Euler Beta 函數——弦論從一個 Gamma 函數公式誕生（後來理解為弦的開弦一階樹級振幅；多弦推廣 = 多元 Beta）。 \(\Gamma\) 的極點（ \(s, t\) 平面的束縛態塔）= 弦的無窮激發譜——「 \(\Gamma\) 的極點結構 = 粒子譜」：8.1 篂「極點決定函數」哲學的物理實現。

場景五： \(p\) -adic 與算術 \(\Gamma\) 。Morita \(p\) -adic \(\Gamma\) （1975）： \(\Gamma_p\) 滿足 \(p\) -進域上的類比性質； \(\Gamma\) 的算術化連接 Iwasawa 理論與 \(\zeta\) 的 \(p\) -進類比（Kubota–Leopoldt \(p\) -adic \(\zeta\) ，1964）——數論在 \(p\) -進世界重建 \(\Gamma/\zeta\) 的整套架構。一個 1729 年的積分函數，在二十一世紀仍是數論與物理的通用接口。

名詞速查表（本節出現的核心術語）#

  • Euler 積分： \(\Gamma(z) = \int_0^\infty e^{-t}t^{z-1}dt\) （ \(\operatorname{Re} z > 0\) ）。
  • 函數方程： \(\Gamma(z+1) = z\Gamma(z)\) 。
  • 餘元公式： \(\Gamma(z)\Gamma(1-z) = \frac{\pi}{\sin\pi z}\) 。
  • Beta 函數： \(B(x,y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}\) 。
  • Weierstrass 乘積： \(\frac1{\Gamma(z)} = ze^{\gamma z}\prod(1+\frac zn)e^{-z/n}\) 。
  • Euler–Mascheroni 常數： \(\gamma \approx 0.5772\) （有理性未解）。
  • Stirling 公式： \(\Gamma(z) \sim \sqrt{2\pi}z^{z-\frac12}e^{-z}\) 。
  • Bohr–Mollerup 定理：對數凸特徵化（1922）。
  • Hankel 圍道： \(1/\Gamma\) 的全平面表示。
  • Veneziano 振幅：Beta 函數的弦論應用（1968）。

年表速查（本節涉及的關鍵年份）#

  • 1729：Euler 致 Goldbach 信，積分定義。
  • 1730–1733：de Moivre、Stirling 漸近與 \(\sqrt{2\pi}\) 。
  • 1738：Euler 發表積分表示。
  • 1749：Euler 餘元公式的原始形式。
  • 1812：Gauss 乘積定義與超幾何論文。
  • 1809–1811：Legendre 命名、倍角公式。
  • 1856：Weierstrass \(1/\Gamma\) 乘積。
  • 1908：Bohr 足球銀牌（同年輩的趣聞）。
  • 1922：Bohr–Mollerup 定理。
  • 1944：Selberg 積分。
  • 1968：Veneziano 振幅（弦論誕生）。
  • 1975：Morita \(p\) -adic \(\Gamma\) 。

深化問答三則#

問：為什麼 \(\Gamma\) 被稱為「特殊函數之母」？ 答：三個理由。結構上： \(\Gamma\) 出現在幾乎所有其他特殊函數的表示中（Bessel 的級數、超幾何的參數、zeta 的函數方程）。歷史上：它是第一個「非初等但系統研究」的函數，其理論（延拓、漸近、乘積、方程）成為特殊函數論的方法模板。應用上：機率、物理、組合的常數歸一化幾乎都經 \(\Gamma\) 。

問： \(\Gamma\) 的極點在負整數——這「應該」嗎？ 答：由函數方程必然： \(\Gamma(0) = \frac{\Gamma(1)}{0}\) 型發散是遞推的宿命。階乘往負方向「爆炸」是遞推結構的幾何（留數 \((-1)^n/n!\) 衰減——極點越深越弱）。函數方程 + 單極點 = 半純函數的「最小病態」實現；Bohr–Mollerup 保證這已是唯一選擇。

問：Stirling 級數發散，為什麼還是「定理」？ 答：因為陳述是漸近的： \(\Gamma(z) = \sqrt{2\pi}z^{z-\frac12}e^{-z}\left(\sum_{k=0}^{N}\frac{a_k}{z^k} + O(z^{-N-1})\right)\) 對每個固定 \(N\) 成立——不是「級數收斂於函數」，而是「函數的餘項可逐階控制」。漸近展開在 \(N\) 增大時先好後壞（最優截斷在誤差最小處）——「發散但有用」的精確語義由餘項界給出（5.2 篂最陡下降法的推論）。

（本節完）