gamma函数的一阶导数值公式(阶乘也很有趣:从阶乘到伽玛函数到非整数的阶乘)
阶乘运算(Factorial)
任何大于等于1 的自然数n 阶乘:
也即
下表给出了一些自然数的阶乘值:
https://en.wikipedia.org/wiki/Factorial
100!是一个158位的整数
100!这么大的数到底怎么算出来的呢?
阶乘的计算
直接求阶乘,需要经过大量的乘法运算,位数太多,计算机也无法表示出来。此时,往往采用对数方法,将阶乘的乘法运算化为加法运算。如
编写一段Python语言代码求等式右边的值:
import math
digit_num =0.0
for i in range(100):
digit_num += math.log10(i+1)
print(digit_num)
运【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的行得到
157.97000365471575(近似值),即
这说明100!是一个158位的数。根据对数函数与指数函数的关系,可以反求出阶乘值:
以前不理解对数意义的朋友这里可以体会到对数的强大威力了吧?
另,阶乘有一个有趣的近似公式:
斯特林()公式- Stirlings approximation
斯特林公式与阶乘曲线对比
我们实际验证一下斯特林公式的误差。将n=100代入上述公式,得到100!≈
9.32484762526934324776475612717870232347097456474180622928179581533688495555540466030862391627555227673250【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的661579827505817302017886487207720230946742094857267442225508190492286520310411195040966964294345297084311638093420567576481015234062861600852667351728186398316114266209416847362850304098552423112683442073070730677904381912557360138125732653622702291187198097261154385694104026076300353130469579563925【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的66366745658132452941877904052886947223641749037779513877635612354880691524914259437590327045612488757528210… × 10^157与我们用对数求得的值之间的误差大约为0.08329%,即万分之8.3,相当精准吧!!!
阶乘的延拓
可以将点(n, n!)即(0, 0!), (1,1!), (2,2!), (3,3!),…在平面坐标系上表示出来。
n!, n=0..4
n!, n=0..6
n!, n=0..10
我们能不能找到一条数学曲线,能够穿越上述所有点(n,n!)呢?找到这样一条曲线的过程就【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的是数学上的解析延拓,从整数域解析延拓到实数域。
伽玛函数
人类恰恰找到了这样一个函数,即伽玛函数(Gamma Function)。伽玛函数的定义如下:
伽玛函数是一个用定积分公式定义的函数,所以求伽玛函数变成了求定积分。不难求得:
进而
伽玛函数与实数域阶层的关系
这些结论我就不做证明了,一方面这些知识可以很便捷地索到,另一也是更重要的方面是,毕竟我的目标不是吓唬大家和显摆自己的学问,而是希望尽可能充分地向大家分享、呈现数学的奥妙、美丽和魅力。
从该等式可以看出,阶乘不就是伽玛函数从实数域降维到整数域的降维函数吗?反之,伽玛函数不正是阶乘序列在从整数域向实数域的延拓吗?
伽玛函数衍生出的一个常数,即为弗朗桑-【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的罗宾逊常数(Fransén–Robinson Constant):
问题:伽玛函数是阶乘运算的唯一解析拓展函数吗?
答案是否定的,因为满足这样的拓展函数有无数个。如如下函数在横坐标为整数时的值也等于对应的阶乘值:
实数域的阶乘函数
因为
也就是说
用如下方式来表示这个阶乘函数:
该阶乘函数有如下递推性质(从小到大,算正数的阶乘时用到):
从上面的递推公式,我们可以得到新的递推公式(从大到小,算负数的阶乘时用到):
我们试着求一下几个非整数实数的阶乘函数值:
根据这个值可以推出其他一系列值:
这π(x)就是实数域阶乘函数的一个合理定义公式。阶乘函数y=π(x)=x!的曲线如下图:
我们发现上述阶乘函数在负整数处不连续,【资 ;源 之.家.】ziyuanzhij. cn 每日持.续更新.可.实操.的即不收敛,与我们计算的结果相符:
最后再求两个特殊的阶乘
其实阶乘还可以延拓到复数域,如
复函数
曲线图如下
(cosx+isinx)!的曲线图
人家在何许?云外一声鸡