Stirling 公式.docVIP

谈 Stirling 公式（转） 彭宇煦 12位粉丝 1楼甲、一个机率问题什麽是一个事件 (event) 的几率？这是机率论最基本也是争论最多的一个问题。 举最简单的例子来说明：丢一个公正铜板 (fair coin)，出现正面 (head) 的机率为 这是什麽意思呢？常识性的解释大致是，将此铜板独立地丢「很多」次，那麽正面出现的次数「大约」占一半，这是在随机的说不准中很确定的事情。所谓的「平均律」(the law of averages) 或「大数法则」(the law of large numbers) 隐隐约约就是指着这个解释。不过，常识往往是含糊的或自相矛盾的，需要加以精炼。事实上，「数学是精炼的常识」(Mathematics is refined common sense)。常识是我们作观念探险之旅的出发点。 问1：丢 2n 次铜板，正面恰好出现 n 次的机率有多大？ 根据组合学，丢 2n 次铜板，共有 22n 种可能结果，假设每一种结果发生的机会均等，那麼 2n 次中有 n 次为正面的结果共有 2nCn 种，故得机率为 我们更有兴趣的问题是，当 n 趋近 时，p2n 会趋近於多少？上述常识性的解释似乎是说，，这成立吗？这需要对(1)式作精确的估算，於是引出了下面的问 2：当 n 很大时，如何估算？更明确地说：当 n 趋近 时，n! 的渐近相等式 (Asymptotically equal formula) 是什麼？即要找一个「好用」(an) 使得 我们希望找到这样的 (an)，然后代入(1)式中计算出极限值 ，就可以检验上述常识性的机率解释是否正确。n! 的渐近相等式存在吗？如何找？这就来到了 Stirling 公式的大门口。在文献上，有许多文章论述 Stirling 公式的简化证明或机率式的证明，不过都只是在已经知道公式后，给出证明而已，并没有说出如何「看出」或「猜出」公式的追寻、探险过程。因此令人有「美中不足」或「未尽妙理」的感觉。本文我们就试著来补上这个缺憾，展示一种推测式的猜想过程。我们不排斥还有其它猜想过程因为登一座山可以有各种不同的路径，路径越多越美妙。PS：本篇内容较多，等我慢慢贴，为保证阅读流畅性，请暂时不要插楼，谢谢！ 顶 0 2009-12-15 16:11 回复 彭宇煦 12位粉丝 2楼乙、n! 的驯服首先观察，欲估计它，最简单的是采高估策略：每个因数都用 n 来取代，亦即取 （2）显然，故 n^n 高估 n!，不过也不错，「万事起头难」，有个开头，就可以逐步修改进，从错误中学习，而达真理的殿堂。[问3：] 如何改进(2)式？ 我们改采中庸策略：每个因数都用「中位数」1/2（差不多就是算术平均(n+1)/2 ）来取代，亦即取(3)这样应该会比(2)式更好才对吧！ 我们必须对(3)式作分析与检验的工作。令 如果，那麼「芝麻就开门」了，an=(n/2)^n就是我们所要的渐近公式。然而我们的内心不禁会响起如下的怀疑：真理不会藏得这麼浅显让我们一猜即中吧？ 我们来比较 n! 与(n/2)^n 的大小。由算术平均大於等於几何平均定理知 事实上可以用数学归纳法证明：因此当 n 很大时，用「相差」的观点来看，(n/2)^n 高估了 n! 但是此地我们应该另采「相比」的观点更适当，因为我们要找的是 n! 的渐近相等式。例如，n^2+n 与 n^2，从「相差」观点来看，当n趋于无穷时，两者之差但是从「相比」观点来看，(n^2+n)/n^2→1，即两者渐近地相等。换言之，「相差」观点的高估，还是有可能是「相比」观点的渐近相等。 考虑 n! 与 （n/2）^n 之比的数列 (bn)，我们的目标是探求极限limbn(n趋于无穷) 。首先注意到 (bn) 是一个递减的正项数列，由实数系的完备性知（4）但是 α 等於多少，并不容易看出。我们采用旁敲侧击的战术，我们观察到下面简单的问题未完待续，晚上来帖，贴图很辛苦，希望不要插楼，谢谢！ 2009-12-15 16:28 回复 彭宇煦 12位粉丝 3楼补题1：设 (Sn) 为一个正项数列。如果limSn=S属于R且S不等于0，则limS(n+1)/Sn=1 这没有告诉我们一个数列何时会收敛，不过有「消极中的积极」作用。如果limS(n+1)/Sn=1 不成立，则可能有三种情形：limSn=0 或limSn=无穷或limSn不存在。此时根本不必梦想会有limSn=1 。另一方面，如果limS(n+1)/Sn=1 ，则 (Sn) 可能收敛，也可能发散；此时也不能保证limSn=1现在就来计算极限limb(n+1)/bn(n→