推广阶乘函数到实数域(伽马函数连乘积)
はるもと
2023年01月29日 17:55

阶乘函数,相信我们不会很陌生. 对于高中阶段使用的阶乘函数(通常记为 n!)的定义域为全体自然数 %5Cmathbb%7BN%7D%5Ccup%5C%7B0%5C%7D,并“规定”了 0!%3D1.

这样的阶乘函数定义为:

n!%3Dn%5Ctimes%20(n-1)%5Ctimes%20%5Ccdots%20%5Ctimes%203%5Ctimes%202%5Ctimes%201

由定义可知阶乘函数有以下性质:

n!%3Dn(n-1)!(n%5Cgeq1)

我们将尝试由此出发. 根据阶乘的定义,我们知道对于一个 k%5Cin%5Cmathbb%7BN%5E%2B%7D,有:

(n%2Bk)!%3Dn!(n%2B1)%5Ccdots(n%2Bk)

写为连乘积形式,有:

(n%2Bk)!%3Dn!%5Cprod_%7Bi%3D1%7D%5E%7Bk%7D(n%2Bi)

同理,我们也有:

(n%2Bk)!%3Dk!%5Cprod_%7Bi%3D1%7D%5E%7Bn%7D(k%2Bi)

两式右边相等,即有:

n!%5Cprod_%7Bi%3D1%7D%5E%7Bk%7D(n%2Bi)%3Dk!%5Cprod_%7Bi%3D1%7D%5E%7Bn%7D(k%2Bi)

我们为什么需要这样的一个包含了两个正整数变量的式子呢?其实是因为我们希望能借助 n%5Cto%5Cinfty 时极限的力量. 如果这样能够帮助化简,那么我们如果能把剩下的 k 给分离出来(指不在巨算符的上下界),那么把 k! 给分离出来再换 k 为 x,想必就能得到关于 x! 的表达式了吧.

可惜的是,左右两边一大串的乘法还是比较难以处理的,如果我们能有方法换成加法,那么或许有更多熟悉的极限技巧可以运用. 那么有没有一种函数能变乘为加呢?有,就是对数.

我们在这里为了方便,使用自然对数.

%5Cln(n!)%2B%5Csum_%7Bi%3D1%7D%5E%7Bk%7D%5Cln(n%2Bi)%3D%5Cln(k!)%2B%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%5Cln(k%2Bi)

当 n%5Cto%5Cinfty 时,对于一个正数 m%5Cln(n%2Bm)-%5Cln(n)%3D%5Cln%5Cleft(1%2B%5Cfrac%7Bm%7D%7Bn%7D%5Cright)%3D%5Cln(1)%3D0,可以视为

%5Cln(n)%3D%5Cln(n%2Bm)

此时我们可以利用这个性质来化简得到:

%5Cln(n!)%2Bk%5Cln(n)%3D%5Cln(k!)%2B%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%5Cln(k%2Bi)

这样,就可以把巨算符上的 k 给成功地拿下来了. 进一步整理得:

%5Cln(k!)%3D%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%5Cln(i)-%5Csum_%7Bi%3D1%7D%5En%5Cln(k%2Bi)%2B%5Cln(n%5Ek)

%5Cln(k!)%3D%5Csum_%7Bi%3D1%7D%5En%5Cln%5Cleft(%5Cfrac%7Bi%7D%7Bk%2Bi%7D%5Cright)%2B%5Cln(n%5Ek)

由此,已经能看见希望了!

%5Cexp(%5Cln(k!))%3Dk!%3D%5Cexp%5Cleft(%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%5Cleft(%5Cfrac%7Bi%7D%7Bk%2Bi%7D%5Cright)%2B%5Cln(n%5Ek)%5Cright)%3Dn%5Ek%5Cprod_%7Bi%3D1%7D%5E%7Bk%7D%5Cleft(%5Cfrac%7Bi%7D%7Bk%2Bi%7D%5Cright)

进而把 k 换成 x,就能得到我们的终极目标:

%5Cboxed%7Bx!%3D%5Clim_%7Bn%5Cto%5Cinfty%7Dn%5Ex%5Cprod_%7Bk%3D1%7D%5E%7Bn%7D%5Cleft(%5Cfrac%7Bk%7D%7Bx%2Bk%7D%5Cright)%7D

我们可以验证一下其是否符合我们先前的规定以及性质:

0!%3D%5Clim_%7Bn%5Cto%5Cinfty%7Dn%5E0%5Cprod_%7Bk%3D1%7D%5E%7Bn%7D%5Cleft(%5Cfrac%7Bk%7D%7B0%2Bk%7D%5Cright)%3D1

%5Cfrac%7Bx!%7D%7B(x-1)!%7D%3D%5Clim_%7Bn%5Cto%5Cinfty%7D%5Cfrac%7Bn%5Ex%5Cprod_%7Bk%3D1%7D%5E%7Bn%7D%5Cleft(%5Cfrac%7Bk%7D%7Bx%2Bk%7D%5Cright)%7D%7Bn%5E%7Bx-1%7D%5Cprod_%7Bk%3D1%7D%5E%7Bn%7D%5Cleft(%5Cfrac%7Bk%7D%7Bx%2Bk-1%7D%5Cright)%7D%3D%5Clim_%7Bn%5Cto%5Cinfty%7Dn%5Cfrac%7B%5Cfrac%7B1%7D%7B(x%2B1)%5Ccdots(x%2Bn)%7D%7D%7B%5Cfrac%7B1%7D%7Bx(x%2B1)%5Ccdots(x%2Bn-1)%7D%7D%3D%5Clim_%7Bn%5Cto%5Cinfty%7Dn%5Cfrac%7Bx%7D%7Bx%2Bn%7D%3Dx

是符合的!至此,我们成功了.