#中阶乘怎么表示_你不知道的阶乘与gamma函数

你不知道的阶乘

阶乘对于有数学基础的人来说都不陌生,简单理解就是数的累乘。10的阶乘10!=10*9*8*7*6*5*4*3*2*1。但是我们有没有思考过,如分数的阶乘是如何运算的?有没有方法估算一个数的阶乘?

96f6423c9810faad3ab67dcd866534eb.png

其实,1/2的阶乘等于π的平方根的一半,本文通过两种方式来求该等式。每一种都很巧妙,看了都能开拓数学思维。一种是利用极限与多项式,一种通过gamma函数。下面我们来证明这个等式。

首先,我们来证明沃利斯公式。

e6f8d613ae0554f13f084a0acc2cc406.png

正弦函数sin x有无穷多个零点0,±π,±2π,±3π,···。一个多项式如果有若干个零点x1,x2,x3,x4 ···,xn,那么该多项式一定可以表示为

39409c48e9b5387e2c1f3f3dcf332ba5.png

可以将sin x大胆展开得

b5d7f7b35aee04a93ef6b145b96e7934.png

将x=π/2代入上式可以得出

bd579d49771cb996f98b8f46cdb308cd.png

最后得到

ab8b5ae02b00b431b7afc29a22fc847a.png

上式就是沃利斯公式。这种证明方法不是特别严格,沃利斯通过求圆弧下的面积同样证明了沃利斯公式,见参考链接《神奇的伽玛函数上》。

证明了沃利斯公式,接下来估算n!,得到n阶乘的一般形式,再求1/2的阶乘,欧拉采用无穷乘积给出了n!的一个插值公式。

2c6853c4c7b5ab689dde8aac2968de0d.png

改成极限形式为

3cee6f24b9d94ecbd8db714dd8dde05f.png

整理式子得

12759a4b02cf5846ec2d5934741c8568.png

则n的阶乘的该插值公式得到证明,可以看出该公式也适用于n是分数的情况,将n=1/2代入得

cbf9309d812036445722eb06ecf3c2e5.png

可以惊奇的发现根号内的式子与沃利斯公式形式几乎一样,只少乘了最前面的因子2。将沃利斯公式代入上式得

1155a70be5ebd0b56f506b95eef18042.png

这样我们就求出了1/2的阶乘的值。

阶乘与gamma函数

gamma函数的一般形式为

66dd098b421e3ad11724bff25dc05153.png

利用分部积分法,可以得出

c7952cff992ca80a180ac50158d33567.png
d5e866afae3ef8d53c03b18e63d39972.png

所以可以得到

d2b33a6432c91bb9b9031e1de064dd3f.png

那么gamma函数的一般形式是如何得出的呢?欧拉通过n的阶乘推导出了gamma函数的一般形式。由于得出了1/2的阶乘的结果中有π的存在,因此欧拉自然联想到阶乘的计算会与积分有关,提出了以下一般的积分形式:

35abd19c47a1a378b9e73e1fa920663c.png

此处n为正整数,e为正实数,利用分部积分法得

7579d4780b2e9d772fbb9e1d649d8c36.png

通过重复迭代上面的公式得

931adb9be968629ca3e562f6a4cf68fb.png

则可以得到求n的阶乘的式子

a0d579334a3876890301f7a5f282deb6.png

现在已经成功的将n的阶乘表示成积分的形式,但是由于n为整数,式子中的非积分部分无法推广分数的情况,因此要继续简化式子。

要让一个量从一个数学等式中消失,数学家惯用的做法就是让这个量取一个极端的值。这里让e趋于无穷。取e = f/g得

30e6c5aa229804b0c760d3163b80fe26.png

然后令f趋于1,g趋于0。左边显然趋于n的阶乘,右边还需要简化计算,令x等于t的h次方,其中h=g/(f+g),得

c30571eec574b07528c83bfe466e41b5.png

当f趋于1,g趋于0时,h显然趋于0,利用洛必达法则,可以得到常用极限

80f4a6c72784bc5c75234ced5c4dfe33.png

对两边同时取极限,见证奇迹的时刻

fff532cb6618e4a0ba39f6bd6e52f333.png

通过n的阶乘推导出了gamma函数的一般形式。

下面我们利用gamma函数来求1/2的阶乘。

由于gamma函数的阶乘形式只满足n为正整数的情况,因为我们要通过gamma函数的一般形式来算1/2的阶乘。我们将n=-1/2代入gamma函数一般表达式得

938efeb6499dd1729b624c7d6ce30a1f.png

仔细观察函数部分与某一种分布的概率密度函数神似,就是正态分布。正态分布的概率密度函数为

344a8d360bedcce8fe26765adab8fb32.png

一般形式的正态分布概率密度函数参数较多,不便于我们观察,我们取正态分布的标准形式,即μ取0,σ取1,则正态分布的概率密度函数变成

79f9ef640292c6172a446f96449ca905.png

关于y轴对称,则大于0的部分积分为1/2,可以得到

a86d3dc04951a46c5a88b6bb661d1b30.png

令x = √(2t),代入上式得

61d80b293b8eda3290526dfbf3a88ce6.png

惊奇的发现积分部分包含我们求的(1/2),于是得

835731a19a14fb5fc2b4dfc0c4966bdf.png

现在我们已经求到(1/2),求1/2的阶乘,也就是(3/2)就很简单了,由于

2e1e2db4f6bddbd2456758b9c20304e7.png

我们可以得出

0fc966a6cf9087e49755e178f37d2c88.png

本文只是简单介绍了gamma函数,接下来的文章会介绍与gamma函数相关的gamma分布与beta分布。

参考链接:

http://www.flickering.cn/数学之美/2014/06/神奇的伽玛函数上/

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值