数学中的概率分析之伯努利大数定律

三、伯努利大数定律

现在我们来介绍伯努利《推测术》的最重要部分――包含了如今我们称之为伯努利大数

定律的第4部分。回到本章开始那个缶中抽球的模型：缶中有a 白球，b 黑球，p ＝a

a b +。有放回地从缶中抽球N 次，记录得抽到白球得次数为X ，以X

N 去估计p 。这个估计法现今仍

是数理统计学中最基本的方法之一。此处的条件是，每次抽取时都要保证缶中a ＋b 个球的每一个有同等机会被抽出。这一点在实践中并不见得容易。例如，产生中奖号码时用了复杂的装置。在实际工作中，统计学家有时用一种叫做“随机数表”的工具。这时一本大书，各页按行、列排列着数字0，1，…,9,它们是用据信是“充分随机”的方法产生的。在使用时，“随机地”翻到其中一页并“随机”点到一个位置，以其处地数字决定抽出地对象。

伯努利企图证明的是：用X

N 估计p 可以达到事实上的确定性――他称为道德确定性。其

确切含义是：任意给定两个数0ε>和0η>，总可以取足够大的抽取次数N ，使事件

X p N ε⎧⎫−>⎨⎬⎩⎭的概率不超过η。这意思是很显然：X

p N ε−>表明估计误差未达到指定的

接近程度ε，但这种情况发生的可能性可以随心所欲地小（代价是加大N ）

。为忠实于伯努利地表达形式，应指出两点：一是伯努利把ε限定为1

()a b −+，虽然其证明对一般ε也有效。他作这一限定与所有缶子模型的特殊性有关：必要时把缶中的白、黑球分别改为ra 和个，

则p 不改变，rb 1

()a b −+改为1

ra rb +，只须r 取足够大，可使此数任意小。其次，伯努利要证

的是：对任给c>0，只须抽取次数N 足够大，可使

X X P p cP p N

N εε⎛⎞⎛−≤>−>⎜⎟⎜⎝⎠⎝⎞⎟

⎠. (8) 这与前面所说是一回事。因为由上式得

1

(1)X P p c N ε−⎛⎞

−><+⎜⎟⎝⎠, (9)

取c 充分大可使它小于η。另外要指出的是：伯努利使用的这个缶子模型使被估计的p 值只能取有理数，因而似乎有损于其结果的普遍性，但其证明对任意的p 成立，故这一细节并不重要。

伯努利上述对事实上确定性的数学理解，即（8）式，有一个很值得赞赏之点，即他在概率论的发展刚起步的阶段，就给出了问题的一个适当的提法。因为，既然我们想要证明的是

当N 充分大时，X

N 和p 可以任意接近，则一个看来直截了当的提法是

lim

N X p N →∞

=, (10)

而这不可能实现。因为原则上不能排除“每次抽到白球”的可能性，这时X

N 总为1，不能

收敛于p<1。或者退一步：要求（10）式成立的概率为1，这个结论是对的，但直到1909年才由波莱尔证明，其难度也比伯努利的提法大得多。设想如当时伯努利就采用这个提法，他也许不一定能在有生之年完成这一工作。波莱尔得结论比伯努利强，故现今把它们得结论分别称为强大数律和弱大数律。

如今具有概率论初步知识的人都知道，伯努利大数律是契比谢夫不等式的简单推论。但在伯努利时代尚无方差概念，更不用说这一不等式了。伯努利用的是直接估计概率的方法，大意如下：令

()

0A P Np X Np N ε=<<+，

((1))A P Np kN X Np k N k εε=+<≤++，k=1,2,……

只须证明：当N 充分大时有（注3）

, (11)

()012A c A A >++⋅⋅⋅这就解决了X>Np 的一边。对X

附带指出：可以把伯努利的结论（9）引申一点点：如果我们知道缶中球的总数a+b ，或者更广一些，知道a ＋b 不超过某已知数M ，则可以把（3）式改进为：可以找到p 的一个估计

（不是ˆ()p X X

N ），使当N 充分大时有 。 （12）

1

ˆ(())(1)P p

X p c −≠<+

但如不给定a+b 的界限，则找不到这样的估计量（注4）。

ˆ()p X 伯努利当初提出的目标，比单纯证明（9）式要高：（9）式只肯定了当取N 充分大时，用X

N

估计p 可达到任意指定的精度ε，而可靠度不小于1

1(1)

c −−+。伯努利希望弄清楚到底需要

N 多大。解决了这个问题，在实用上就可以根据所需的精度和可靠度，去规划所须观测次数N 。他证明了以下的结果：定义

＝不小于1m

[]

log (1)log(1)log c b a −+−a 的最小整数， （13）

＝不小于2m

[]

log (1)log(1)log c a b −+−b 的最小整数， （14）

111()()(11m a b b a b m N a +++−=

+)

， （15） 222()()(11m a b a a b m N b )+++−=

+。 （16）

则取

能满足（9）式。伯努利给了若干数字例子，其一为：a ＝30，b ＝20

（p ＝max(,)

12N N N =35）

，1

50ε=

，c ＝1000。用上述结果算出所需的次数N 为25550。可以与由契比谢夫不等式计算的结果作一比较。按此不等式，有（注5）

23

113260015505055X P N N N −⎛⎞⎛⎞⎛⎞⎛⎞−−>≤=

⎜⎟⎜⎟⎜⎟⎜⎟⎝⎠

⎝⎠⎝⎠⎝⎠。 为使此值不超过

1

1(1)1001c −+=

，N 至少应为600600，这比伯努利给出的值大20多倍。

这反映了一个事实：伯努利在证明（9）式中所作的概率估值，比契比谢夫不等式所作出的要精细得多。虽然如此，25550这个数仍嫌过大。美国统计史学者斯蒂格勒认为，伯努利之所以久未发布其研究成果，与他对一点的不满意有关。因为在伯努利时代，一个中等城市的规模尚不过几千人，25550简直可算时“天文数字”。不过，后世的学者所看重的不在这些地方。如今大家都公认由伯努利工作发端的大数定律已成为整个数理统计学的基础。人们也对伯努利工作的哲学意义给予很高的评价。如斯蒂格勒指出：伯努利证明了数学家不仅可以后验地认识世界，还可以用数学去估量他们的知识的限度。伯努利在结束《推测术》时就其