§3 Bernoulli试验和直线上的随机游动
文章目录
1.Bernoulli 概型
定义2.3.1(Bernoulli)试验
只有两种可能结果的实验被称为 Bernoulli试验。在这类问题中,我们可将事件域取为:
F={∅,A,A,Ω}
并称出现 A 为“成功”,出现 A 为“失败”。
在某些情况中,试验的结果不止两个。但是,如果我们可以认为:出现且仅仅出现某一个特定的结果时认为其“成功”,其余任何结果都视为“失败”,则此时我们也可以把问题视为 Bernoulli 试验。
在 Bernoulli 试验中,首先需要给出下面概率:
P(A)=p, P(A)=q
显然, p,q⩾0,且 p+q=1.
现在,我们考虑重复进行 n 次独立的 Bernoulli 试验。
定义2.3.2(n 重Bernoulli试验)
在每次试验中事件 A 和事件 A 出现的概率都保持不变。我们称这样的试验为 n重Bernoulli试验 ,记为 En.
对于 n 重Bernoulli试验,我们有下列的四个约定:
- 每次试验最多出现两个可能结果之一: A 或 A ;
- A 在每次实验中出现的概率 p 保持不变;
- 各次实验相互独立;
- 总共进行 n 次试验。
下面,我们给出 n 重Bernoulli 试验的概率空间:
n 重Bernoulli试验 En 的样本点形如:
(A^1,A^1,⋯,A^n)
其中, A^i 是 Ai 或 Ai ,分别表示第 i 次试验中出现 A 或 A.显见,这样的样本点总共有 2n 个,这是一个有限样本空间。
为书写方便起见,我们进行如下约定:
(A1,A2,⋯,An−1,An)
表示前 n−1 次试验均出现事件 A ,而第 n 次试验出现事件 A,简记为
A1A2⋯An−1An.
为给定 n 重Bernoulli试验的样本点的概率,我们主要考察在它中 A 或 A 出现的次数。若其中有 l 个 A ,从而有 n−l 个 A ,则利用试验的独立性可知:
P(A^1A^1⋯A^n)=P(A^1)P(A^2)⋯P(A^n)
=plqn−l.
一般事件的概率由它所含样本点的概率求和得到。这样一来,我们已对 n 重Bernnoulli试验给定了概率空间。
Bernoulli试验是一种重要的概率模型。它是“在相同条件下进行重复试验”的一种数学模型,尤其在讨论某事件出现的频率时常用这种模型。
2. Bernoulli概型中的一些分布
2.1 Bernoulli 分布
若我们只进行一次Bernoulli试验,则这种概率分布称为 Bernoulli分布,这是最简单的情况。
2.2 二项分布
我们记 n 重Bernoulli试验中事件 A 出现 k 次的概率,记之为 b(k;n,p):
若以 Bk 记 n 重Bernoulli试验中事件 A 恰好出现 k 次这一事件,而用 Ai 表示第 i 次试验中出现事件 A,以 Ai 表示第 i 次试验中出现 A ,则
Bk=A1A2⋯AkAk+1⋯An+⋯+A1A2⋯An−kAn−k+1⋯An.
右边的每一项表示在某 k 次试验中出现事件 A ,而在另外 n−k 次实验中出现 A,这种项共有 (kn) 个,而且两两互不相容。
显见右边各项所对应事件的概率均为 pkqn−k ,我们利用概率的可加性即得:
P(Bk)=(kn)pkqn−k
即:
b(k;n,p)=(kn)pkqn−k, k=0,1,2,⋯,n
注意到 b(k;n,p), k=0,1,2,⋯,n 是二项式 (q+ps)n 展开式中 sk 项的系数,因此上式称为 二项分布 。
特别地:
k=0∑nb(k;n,p)=k=0∑n(kn)pkqn−k=(p+q)n=1.
2.3 几何分布
下面, 我们讨论在 Bernoulli 试验中"首次成功"出现在第 k 次试验的概率. 要使"首次成功"出现在第 k 次试验时, 必须且只需在前 k−1 次试验中均出现事件 A, 而在第 k 次时出现 $A$, 因此该事件 (我们记为 Wk ) 可表示为
Wk=P(A1)P(A2)⋯P(Ak−1)P(Ak)=qk−1p.
记
g(k;p)=qk−1p, k=1,2,⋯
g(k;p) 是几何级数的一般项, 故称上式为 几何分布. 此处有;
k=1∑∞g(k;p)=k=1∑∞qk−1p=p1−q1=1.
几何分布给出了等待事件 A 出现共需要试验 k 次的概率. 这类概率在许多问题中都有涉及.
2.4 Pascal分布
下面, 我们讨论更为复杂的情况,也就是 Pascal 分布. 可以视其为几何分布的一种推广.
考虑 Bernoulli 试验,让我们考察要经过多长时间才会出现第 r 次成功:
若第 r 次成功发生在第 ζ 次试验,则必然有 ζ⩾r.
下面以 Ck 表示"第 r 次成功发生在第 k 次试验" 这一事件,并且记其概率为 f(k;r,p).
Ck 发生当且仅当前面的 k−1 次试验中有 r−1 次成功, k−r 次失败,而第 k 次试验的结果为成功,这两个事件的概率分别为:(r−1k−1)pr−1qk−r 和 p. 利用试验的独立性,得到:
P(Ck)=(r−1k−1)pr−1qk−r⋅p=P(Ck)=(r−1k−1)prqk−r.
即:
f(k;r,p)=(r−1k−1)prqk−r, k=r,r+1,⋯
注意到:
k=r∑∞f(k;r,p)=k=r∑∞(r−1k−1)prqk−r=l=0∑∞(r−1r+l−1)prql=l=0∑∞(lr+l−1)prql=l=0∑∞(l−r)(−1)lprql=pr(1−q)−r=1.
我们称 f(k;r,p) 为 Pascal分布. 特别地,当 r=1 时,我们得到几何分布.
3. 直线上的随机游动
定义2.3.3 (随机游动)
我们考虑位于 x 轴上的一个质点,并限制它只能位于整数点. 在时刻 t=0 时, 它处于初始位置 a a∈Z , 以后每个单位时间他总受到某个外力的随机作用使位置发生变化,分别以概率 p 和概率 q=1−p 向正方向/负方向移动一个单位.
在这个问题中, 我们所感兴趣的是质点在时刻 t=n 时的位置. 我们称用这种方式描述的质点运动为 随机游动.
定义2.3.4(无限制随机游动和有吸收壁的随机游动)
若质点可在整个数轴的整数点上游动,则称这种随机游动为 无限制随机游动.
若在某点 d 设有一个吸收壁, 质点一旦到达该点则被吸收从而不再游动, 因而游动过程结束,称这种随机游动为 在 d 点有吸收壁的随机游动.
此外, 我们还可以相应地考虑带有 “反射壁” 及 “弹性壁” 的随机游动, 在一个随机游动中还可以存在不止一个壁.
当 p=q=21 时,称这样的随机游动为 对称的, 此时质点向两方向移动的可能性相等.
在自然科学中, 我们可以将大量的问题都归结为随机游动问题. 例如: 随机游动模型可作为布朗运动的初步近似. 概率论中的一些古典问题也可引导到随机游动问题. 实际上, 随机游动可以视为 Bernoulli试验的一种描述法.
下面,我们简介随机游动的两个最简单的模型:
3.1 无限制随机游动
假定质点在时刻 0 从原点出发, 以 Sn 记它在时刻 t=n 时的位置. 为了使质点在时刻 t=n 时位于 k (k 也可以是负整数), 质点必须而且只需在前 n 次游动中向右游动的次数比向左游动的次数多 k 次.
若以 x 记它在前 n 次游动中向右游动的次数, y 记向左游动的次数,则
{x+y=nx−y=k
即 x=2n+k, 因为 x 是整数, 故 k 必须和 n 具有相同的奇偶性.
事件 {Sn=k} 发生相当要求在前 n 次游动中有 2n+k 次向右, 2n−k 次向左. 利用二项分布即得:
P{Sn=k}=(2n+kn)q2n−kp2n+k.
当 k 和 n 奇偶性相反时, 概率为 0.
3.2 两端带有吸收壁的随机游动
假定质点在时刻 t=0 时, 位于 x=a, 而在 x=0 和 x=a+b 处各有一个吸收壁,我们求质点在 x=0 被吸收或在 x=a+b 被吸收的概率. 使用的是差分方程法:
4. 推广的 Bernoulli 试验与多项分布
二项分布可以被很容易地推广到 n 次重复独立试验且每次试验可能有若干个结果的情形. 将每次试验的可能结果记为 A1,A2,⋯,Ar, 而 P(Ai)=pi, i=1,2,⋯,r, 且
p1+p2+⋯+pr=1, pi⩾0.
当 r=2 时, 我们即得到 Bernoulli 试验.
在这种推广的 Bernoulli 试验中,不难推出在 n 次试验中 Ai 出现 ki 次 (i=1,2,⋯,r) 的概率为;
k1!k2!⋯kr!n!p1k1p2k2⋯prkr.
此处 ki⩾0, 且 k1+k2+⋯+kr=n.
上式称为 多项分布, 因其为(p1+p2+⋯+pr)n 展开式的一般项,且知:
ki⩾0 k1+k2+⋯+kr=1∑k1!k2!⋯kr!n!p1k1p2k2⋯prkr.
显然,多项分布是二项分布的推广, 二项分布的许多结论均和多项分布的场合平行. 以后, 我们只详细讨论二项分布的有关问题.