伽马函数
$$\Gamma(x)=\int_0^{+\infty}t^{x-1}e^{-t}dt$$
作为阶乘的推广,会让很多初学者感到困惑,对于笔者来说也不例外。一个最自然的问题就是:这般复杂的推广公式是如何得到的?

在cos.name的文章《神奇的伽马函数》中,有比较详细地对伽马函数的历史介绍,笔者细读之后也获益匪浅。但美中不足的是,笔者还是没能从中找到引出伽马函数的一种“自然”的办法。所谓“自然”,并不是说最简单的,而是根据一些基本的性质和定义,直接把伽马函数的表达式反解出来。它的过程和运算也许并不简单,但是思想应当是直接而简洁的。当然,我们不能苛求历史上伽马函数以这种方式诞生,但是作为事后探索是有益的,有助于我们了解伽马函数的特性。于是笔者尝试了以下途径,得到了一些结果,可是也得到了一些困惑。

在实数的世界里:干不了!

伽马函数作为阶乘函数的推广,自然要满足阶乘函数所具有的性质,比如阶乘函数满足
$$n!=n\times (n-1)!\tag{1}$$
现在我们定义伽马函数$\Gamma(x)$,在$x$是正整数的时候有
$$\Gamma(x)=(x-1)!\tag{2}$$
于是根据$(1)$有
$$\Gamma(x+1)=x\Gamma(x)\tag{3}$$
上式开始只对$x$为正整数时成立,而作为推广,我们要求$(3)$对于所有的实数都成立(这正是伽马函数的性质),而为了求出$\Gamma(x)$的具体表达式,我们就得从$(3)$中解出$\Gamma(x)$,这是函数方程的求解问题。根据《算符的艺术:差分、微分与伯努利数》中的思想,记$D\equiv \frac{d}{dx}$,我们就可以将$(3)$写成
$$e^D\Gamma(x)=x\Gamma(x)\tag{4}$$
利用傅里叶变换来求解它,两边作傅里叶变换,其中
$$D\to i\omega,\quad x\to i\frac{d}{d\omega}$$
从而$(4)$变为
$$e^{i\omega}\mathcal{F}[\Gamma(x)]=i\frac{d}{d\omega}\mathcal{F}[\Gamma(x)]$$
这方程是很容易解的,解为
$$\mathcal{F}[\Gamma(x)]=C\exp\left(-e^{i\omega}\right)$$
最后作逆变换得到
$$\Gamma(x)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} C\exp\left(-e^{i\omega}\right)e^{i\omega x}d\omega\tag{5}$$
似乎挺完美的,但是请注意,积分$(5)$不会得到一个良好的函数(除非$C=0$),比如说$x$是整数的时候,被积函数是一个周期函数!非零的周期函数自然不会积出一个有意义的结果来。严格来说,上述积分并不存在,它的结果只能算是一个广义函数。因此,我们走这条路失败了。事实上,这表明$\Gamma(x)$的傅里叶变换不是一个普通的函数,因此不能像研究普通函数那样研究$\Gamma(x)$的傅里叶变换,所以也就就没有从它的傅里叶变换反过来求伽马函数的做法了。

在复数的世界里:柳暗花明

数学中往往有这样的现象:在某个范围内解决某个特殊的问题显得特别困难,但是把问题推广了,把范围扩大了,一次性解决更多的问题,反而显得更加简洁,如前不久我们谈到过的《正弦级数和余弦级数》就是一例。而在这篇文章中,我们把$\Gamma(x)$看成是阶乘函数对一般实数的推广,但是前面我们没有办法从$(3)$式中反解出$\Gamma(x)$来。有一个奇妙的技巧,让我们可以从中得到有意义的结果,那就是干脆把$\Gamma(x)$延拓到复数域中去,把它认为是复数域中的函数$\Gamma(z)$!这样做能够让我们从$(3)$式中反解出$\Gamma(x)$来,经历了实数世界中的“山重水复疑无路”之后,复数的世界给我们“柳暗花明又一村”的感觉。

设$x$是一个实数变量,那么根据$(2)$式,对于纯虚数$xi$,我们有
$$\Gamma\left((x-i)i\right)=\Gamma(x i+1)=x i\Gamma(x i)\tag{6}$$
同样设$D\equiv\frac{d}{dx}$,那么类似$(4)$,有
$$e^{-i D}\Gamma(x i)=x i\Gamma(x i)\tag{7}$$
同样利用傅里叶变换,此时得到
$$e^{\omega}\mathcal{F}[\Gamma(x i)]=-\frac{d}{d\omega}\mathcal{F}[\Gamma(x i)]$$
解得
$$\mathcal{F}[\Gamma(x i)]=C\exp\left(-e^{\omega}\right)$$
作逆变换得
$$\Gamma(x i)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} C\exp\left(-e^{\omega}\right)e^{i\omega x}d\omega\tag{8}$$
那么
$$\Gamma(z)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} C\exp\left(-e^{\omega}\right)e^{\omega z}d\omega\tag{9}$$
注意不同于$(5)$式,$(9)$式是一个简单的实积分,它可以继续化简
$$\begin{aligned}\Gamma(z)=&\frac{1}{2\pi}\int_{-\infty}^{+\infty} C\exp\left(-e^{\omega}\right)e^{\omega z}d\omega\\
=&\frac{C}{2\pi}\int_{-\infty}^{+\infty} \exp\left(-e^{\omega}\right)e^{\omega (z-1)}e^{\omega}d\omega\\
=&\frac{C}{2\pi}\int_{0}^{+\infty} \exp\left(-e^{\omega}\right)e^{\omega (z-1)}de^{\omega}\\
=&\frac{C}{2\pi}\int_{0}^{+\infty} e^{-t} t^{z-1}dt\end{aligned}$$
最后根据$\Gamma(1)=0!=1$,确定$\frac{C}{2\pi}=1$,从而
$$\Gamma(z)=\int_{0}^{+\infty} e^{-t} t^{z-1}dt\tag{10}$$
$(10)$式正好是我们要寻求的伽马函数的表达式。

这表明,伽马函数作为复自变量的虚部的一个实函数(也就是说,固定自变量的实部,虚部为变量),才具有意义普通函数性质的傅里叶变换。

简述

本文是笔者试探式的结果,通过求解函数方程$(3)$来反解伽马函数,我们看到了傅里叶变换在其中的作用,但也看到了某些不足,一些步骤需要一些特定的技巧才能进一步求解下去。当然,更重要的是不断尝试,不断探索,才能真正去体味到其中的乐趣。本文的后半部分,也就是在复数中求解的思路,也是经过长时间的尝试无果,在今天突然想到的。总的来说,多多尝试,总会有收获的。

转载到请包括本文地址:https://www.kexue.fm/archives/3108

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

如果您需要引用本文,请参考:

苏剑林. (Dec. 08, 2014). 《伽马函数的傅里叶变换之路 》[Blog post]. Retrieved from https://www.kexue.fm/archives/3108

@online{kexuefm-3108,
        title={伽马函数的傅里叶变换之路},
        author={苏剑林},
        year={2014},
        month={Dec},
        url={\url{https://www.kexue.fm/archives/3108}},
}