概率
麦克斯韦说:“我们这个世界的真正逻辑寓于概率的计算之中。”这句话强调了概率的空前的重要性。但仔细想想,会发现“概率”究竟是什么并不是一个容易回答的问题。我们平时在生活中经常使用“机会”或“可能性”这样的词,通常我们是在做一种猜测。为什么我们要猜测呢?因为我们希望在不完全掌握所有信息时做出我们必须做出的决定。事实上,我们的任何一个判断本质上都是一种猜测。所谓“概率”,就是指我们在大量重复某个观察时对其中出现的某个特定结果的最有可能分数的猜测(估计),这时候我们就把发生的次数除以总次数称为“概率”。例如当我们重复了100000次抛硬币的实验,观察到有将近50000次实验中正面朝上,我们就说正面朝上的概率是1/2。
由此可见,我们只能对可重复的观察谈概率。概率有赖于我们的知识以及进行估计的能力,因此它并不是完全客观的!很有可能,当我们所掌握的知识发生变化后,对事物的概率的估计就会变得完全不同。
随机事件
概率空间
为了严格地描述“概率是什么”,人们抽象出了概率空间的概念。概率空间由三个要素组成:样本集Ω,事件集F,概率函数Pr。我们可以把它写成三元组(Ω,F,Pr)。样本集是所有可能结果的集合(我们总是默认Ω是可数集);F=2Ω,它是Ω的所有子集的集合;Pr是F→R的函数,对应一个事件发生的概率。
只有样本集中的每个元素是现实世界中产生的结果,并且一旦随机事件进行了,一定只会出现样本集中的一个元素的结果。而“事件”是我们对“可能的结果”的一种描述,一个“事件”可能描述了许许多多种可能发生的结果。概率函数就是要描述自然语言描述的那些可能结果合在一起发生的可能性。比如,我们尝试用概率空间(Ω,F,Pr)来描述“掷一次骰子”这件事。那么样本集只能是Ω={1,2,3,4,5,6}。而我们的事件可以是“掷出偶数”,它对应着A={2,4,6}∈F,对应的概率Pr(A)=1/2。
事件是集合,因此我们可以用集合的符号来表示来表示事件:例如A∪B表示两个事件的并集,A∩B表示两个事件的交集。我们还可以定义A=Ω∖A,表示A的补事件。
我们发现,概率函数的取值不能是任意的。首先,它的取值必须在[0,1]之间,并且满足Pr(∅)=0,Pr(Ω)=1。由此可见Pr(A)+Pr(A)=1。我们还必须规定,对于两两不相交的有限或可数无穷序列A1,A2,⋯,必须满足Pr(i≥1⋃Ai)=i≥1∑Pr(Ai)。这些事实称为“概率公理”,它描述了符合我们日常经验的概率的性质。必须满足这些条件才能称为“概率空间”。
Ω中的每一个元素本身就可以单独构成一个事件(集合),这些集合两两不交;同时,任意一个事件都是Ω的一个子集,因此一定是这些集合中的某一些的无交并。所以我们就把这些单独构成的集合称为“基本事件”,任何一个事件都是一系列基本事件的并,求任何一个事件的概率只需要把对应的基本事件的概率累加。所以,只要直到了所有基本事件的概率,理论上就可以计算出任何事件的概率。
容斥原理
对于两个有交集的事件A,B,如何求出Pr(A∪B)?直觉就是Pr(A∪B)=Pr(A)+Pr(B)−Pr(A∩B)。但我们必须用公理来推导它。根据公理,我们只会计算无交集的两个集合的并的概率,所以我们做这样的的拆分:Pr(A)=Pr(A∖(A∩B))+Pr(A∩B),Pr(B)=Pr(B∖(A∩B))+Pr(A∩B)。而Pr(A∪B)=Pr(A∖(A∩B))+Pr(B∖(A∩B))+Pr(A∩B)。这样我们就解出来Pr(A∪B)=Pr(A)+Pr(B)−Pr(A∩B)。
用同样的方法,类比集合中我们采用的算贡献的想法,可以证明概率的容斥原理Pr[i∈[n]⋂Ai]=J⊆[n]∑(−1)∣J∣⋅Pr(j∈J⋂Aj)。(证明每个基本事件在等式两边的贡献相等)
Union Bound
从二元情形的容斥原理Pr(A∪B)=Pr(A)+Pr(B)−Pr(A∩B)中可以看出,如果我们丢掉−Pr(A∩B)这一项,因为概率函数一定是非负的,我们就会直接得到Pr(A∪B)≤Pr(A)+Pr(B),推广到n元就是Pr(i∈I⋃Ai)≤i∈I∑Pr(Ai),这个结论称为“Union Bound”。
这个结论虽然推导很简单,但却在应用中有着巨大的价值。
为了看到Union Bound的巨大价值,我们来看一个复杂的例子。
我们发现一个事实,6个人之间要么至少有三个人互相之间都认识,要么至少有三个人之间互相都不认识。用一般的数学语言来描述,在6个节点的完全图(K6)上,每条边要么染成红色要么染成蓝色(代表认识与不认识),那么对于任何一种染色方案,“所有边都是红色的K3子图”和“所有边都是蓝色的K3子图”都至少存在一个。
我们可以证明(枚举),在5个点的完全图上这一点是做不到的。也即要做到这一点至少需要6个节点,在5个点的图上至少存在一种染色方案使得红色和蓝色的K3子图都不存在。
由此我们定义Ramsey数R(r,s),他表示使得完全图KN的边任意红蓝染色始终存在“红色的Kr子图”或“蓝色的Ks子图”中的至少一个的最小的N。这是Well-Defined的,因为我们已经看到太小的N是不行的,而显然N足够大时一定是可行的,所以它被限制在一个闭区间里,所以最小值一定存在。上面这个例子可以描述为R(3,3)=6。
我们可以证明(r−1r+s−2)是R(r,s)的一个上界,这并不需要用到概率的知识:我们首先证明R(r,s)≤R(r−1,s)+R(r,s−1)。令M=R(r−1,s)+R(r,s−1),只需证明KM中一定存在红色Kr或蓝色Ks。对于任意某个特定的染色,考虑与某个点(比如点1)相邻的M−1条边的染色情况,由于M=R(r−1,s)+R(r,s−1),因此根据鸽巢原理,这M−1条边中“至少有R(r−1,s)条红边”与“至少有R(r,s−1)”条蓝边中总有一个成立,不然总边数就小于等于M−2,矛盾。如果是前者满足,那么把这些边连出去的点收集在一起就有R(r−1,s)个点了,根据定义,这些点构成的完全子图里一定要么存在一个红色Kr−1或者蓝色Ks。如果存在蓝色Ks则证毕,否则对于得到的Kr−1附加上原先固定的点1就一定能得到一个红色Kr,因此无论如何都成立。“至少有R(r,s−1)”的情况也是完全相同的。所以我们证明了这个递推不等式。对r+s归纳,那么R(r,s)≤R(r−1,s)+R(r,s−1)≤(r−2r+s−3)+(r−1r+s−3)=(r−1r+s−2),得证。
下面我们证明下界。这里我们只能给出r=s时的下界R(s,s)>2s/2。只需证明对于M=2s/2,存在一种染色方案使得KM中不存在任何同色的Ks子图。这里我们要用到概率方法——转而证明如果给KM随机着色,出现“不存在同色Ks子图”这一事件的概率大于0,这等价于“存在同色Ks子图”这一事件的概率小于1。在这里,概率空间中样本集就是所有可能的染色方案(每条边可以选择染红或染蓝,共有2∣E∣种)。“存在同色Ks子图”这一事件对应着样本集中的一个子集,即那些符合我们条件的染色方案。这是不容易计算的。然而,如果我们枚举所有的大小为s的子图,那么每一个子图的“同色”也是一个事件。并且我们发现对于我们枚举的所有的子图,这些事件的并集恰好就是“存在同色Ks子图”这一事件(如果一个子图同色,那么它在“存在同色Ks”里有贡献;如果它在“存在同色Ks”里有贡献,那么一定在适当的时候被枚举到)。这些事件之间可能有交集,因此精确计算他们的并集的概率必须要容斥,但我们可以用Union Bound给出上界,而对于每个事件的概率是容易计算的——只需保证枚举到的点集里的所有边同色,边的个数为(2s),因此每个事件的概率就是[21](2s)×2(乘以2是因为有两种可能的“同色”)。因此总的概率就可以被放缩为(sM)×21−(2s),代入M=2s/2,s!(2s/2−s)!2s/2!×2×22s(s−1)1,它是小于1的!这样证明就结束了。
当我们要证明“可能出现不存在的情况”时,只需证明“存在的概率小于1”。为此,我们可以枚举所有可能情况,这些可能情况的并集恰好构成“存在”这一事件。我们不直接计算这个并集的概率,而是计算每个可能情况的概率之和,这个和通过Union Bound给出了“存在的概率”一个上界,如果这个上界都已经小于1,那么我们想证明的概率就一定小于1。这就是证明“存在”的一种概率方法。
独立事件
对于概率空间(Ω,F,Pr)中的两个事件A,B∈F,定义两个事件是独立的当且仅当Pr(A∩B)=Pr(A)⋅Pr(B)。对于多个事件A1,⋯,An,定义它们“互相独立”当且仅当对于任何I⊆[n],Pr(i∈I⋂Ai)=i∈I∏Pr(Ai)。
特别需要注意的是,“互相独立”的定义与“两两独立”有所不同。它们并不是一回事。“互相独立”的要求更高。我们将会在随机变量一节中给出两两独立但并不“互相独立”的反例。
“独立”这个概念根据字面意思来理解,就是两个事件之间“没有关联”。比如就“在1,2,3,4中选数”这个样本空间中,“选到1或2”这个事件对应着A={1,2},“选到1或3”这个事件对应着B={1,3},那么发现Pr(A∩B)=Pr(A)Pr(B),所以这两个事件根据我们的定义就是独立的!
我们来看一个多项式恒等判定的例子。如果两个多项式P(x),Q(x)是黑箱,那么如何判定它们是否恒等?假设它们的次数都不超过d次,那么假如我们选取d个不同的实数xi,对于每个都满足P(xi)=Q(xi),那么对于多项式T(x)=P(x)−Q(x),x1,⋯,xd就是T(x)的d个根。如果我们再选一个xd+1,依然满足T(xd+1)=0,那说明T有d+1个根——这只可能是T(x)≡0,因为代数基本定理告诉我们一个d次多项式最多只能有d个不同的实数根。
由此我们可以采用概率方法来判定两个多项式是否恒等。我们选取一个样本集A,里面是一些不同的实数(整数或许更加方便)。均匀随机(Uniformly At Random, u.a.r.)从A中选取一个数字a,如果P(a)=Q(a)就返回二者恒等,否则返回二者不等。这个算法怎么样呢?如果P,Q确实恒等,那么我们的算法始终返回的是正确结果;如果P,Q不等, 我们的算法返回“恒等”的概率是多少?由于P,Q不等, 能满足P(xi)=Q(xi)的xi在实数范围内就一定不超过d个,因此我们抽取的a恰好满足的概率不超过∣A∣d。因此如果选择大小为100d的样本集,算法的出错概率就达到了≤1001。
由于∣A∣的增大是有代价的,使得∣A∣越大并不一定越是好事。如何进一步提高算法的正确性呢?如果我们多次选择a,比如选t次,只有当每一次都满足P(ai)=Q(ai)才返回“恒等”,否则只要有一个不相等就返回“不恒等”,这样的算法在多项式明明不恒等却返回恒等的概率是多少?此时我们的样本空间变成了At,其中的基本事件为“一种取t次的方法”。我们考虑“第i次选到的数x满足P(x)=Q(x)”这一事件发生的概率,记为Pr[P(ai)=Q(ai)]。我们可以验证P(ai)=Q(ai)与P(aj)=Q(aj)是独立的,因为Pr[P(ai)=Q(ai)∩P(aj)=Q(aj)]=∣A∣2∣D∣2,其中D是A中P(x)−Q(x)的根的集合;而Pr[P(ai)=Q(ai)]=Pr[P(aj)=Q(aj)]=∣A∣∣D∣,这样我们就验证了独立性,并且可以用类似的方法验证第1次到第n次每一次满足P(ai)=Q(ai)所有这些事件是“两两独立的”,因此直接得到Pr[i∈[t]⋂P(ai)=Q(ai)]=i=1∏tPr[P(ai)=Q(ai)]=∣A∣t∣D∣t,当∣A∣取100d时,可放缩为100t1。
条件概率
对于概率空间(Ω,F,Pr)中的两个事件A,B∈F,定义A事件在B事件发生下的条件概率为Pr(A∣B)=Pr(B)Pr(A∩B),记为A⊥B。这是符合我们的直观的,当我们讨论“B发生的条件下”时,我们把讨论区域限制在了B内,因此B就“好像是整个空间”了——所以我们除掉B发生的概率。
两个互相独立的事件其中一个在另一个的条件下发生的概率应该就是它本身发生的概率,代入验证这确实成立:由于Pr(A∩B)=Pr(A)⋅Pr(B)恰好得到Pr(A∣B)=Pr(B)Pr(A∩B)=Pr(B)Pr(A)Pr(B)=Pr(A)。这也是我们理解独立的一种方法,两事件独立等价于其中一个在另一个发生下的条件概率就是它自己,这表明A,B是“没有关系”的。
我们可以把条件概率的等式变形为Pr(A∩B)=Pr(A∣B)⋅Pr(B)。于是我们能够得到一个求解事件的交集的概率的链式法则:Pr(A1∩⋯∩An) =Pr(An∣A1∩⋯∩An−1)⋅Pr(A1∩⋯∩An−1) =⋯ =Pr(A1)⋅Pr(A2∣A1)⋅Pr(A3∣A1∩A2)⋯∩Pr(An−1∣A1∩⋯∩An−2)⋅Pr(An∣A1∩⋯∩An−1)。它在直观上似乎更好理解,因为所有事件同时发生的概率可以写作“第一个事件发生的概率乘上在第一个事件发生的条件下第二个事件发生的概率乘上第一个事件和第二个事件都发生的条件下第三个事件发生的概率……”
Birthday Paradox
23个人里有两个人生日在同一天的概率达到50%,60个人里有两个人生日在同一天的概率达到99%。如何计算这个“生日悖论”的概率?
抽象为数学问题,我们可以这样描述:把m个小球随机放进n个箱子里,存在一个箱子里有至少两个小球的概率是多大?我们的概率空间是所有从[m]到[n]的映射f,要计算的是“映射不为单射”这一事件的概率。为此,我们再考虑一系列事件Ai,其中Ai表示“f(i)=f(j)对j∈[i−1]恒成立”这一事件,它对应了一系列映射,直观上表示的是如果我们把丢球的过程拆分成一个个按顺序丢球的过程,这个事件代表第i个球在丢的时候箱子里还是空的。于是“映射不为单射”这一事件等价于所有Ai都发生,因此我们要算的就是Pr[A1∩A2∩⋯∩Am]。利用链式法则,我们可以把它转化为条件概率:Pr[Ai∣A1∩⋯∩Ai−1]是容易计算的,由于前i−1个球都对应了不同的箱子,留下的空箱子只有n−(i−1)个,因此它的值就是nn−i+1。于是求得Pr[A1∩A2∩⋯∩Am]=i=1∏mnn−i+1,它可以写作(1−n1)⋯(1−nm−1),根据1+x≤ex恒成立,1−nk<e−nk,因此原式放缩为e−n1[1+2+⋯+(m−1)]=e−2nm(m−1)。
取n=365,则当m=23时值恰好在0.5左右。当m=60时它小于0.01,因此生日重复的概率高达99%。
全概率公式
对于A,B∈F,恒成立Pr(A)=Pr(A∩B)+Pr(A∩B)。它可以推广,因为关键在于B,B构成了全集Ω的一个划分。如果全集能被划分为n个部分B1,⋯,Bn,那么
∀A,Pr(A)=i=1∑nPr(A∩Bi)。
给定三个n×n 的矩阵A,B,C,验证是否成立AB=C。
如果采用随机算法,我们可以随机一个长度为n的01向量x,验证是否成立ABx=Cx,如果成立返回等于,否则返回不等于。这和验证多项式是否恒等的例子很相似,我们来分析这个算法出错的概率——同样的,只有可能在AB=C而ABx=Cx时算法出错。由于移项可得AB−C=0,(AB−C)x=0,我们的问题等价于验证对于给定的一个矩阵D且D=0,要求u.a.r随机一个01向量x时Dx=0的概率。
注意,D是固定的。我们的随机在于向量x的随机,因此我们的样本空间就是{0,1}n,我们想求出Pr[Dx=0]。我们做一点放缩,如果只要求满足i=1∑nD1ixi=0,即只令Dx的第一维坐标取0,那么满足要求的x相比于Dx=0更多了,换言之事件Dx=0被包含于事件i=1∑nD1ixi=0,因此肯定有Pr[Dx=0]≤Pr[i=1∑nD1ixi=0]。这里我们不妨假设存在一个D1k=0(不然这一行D全0,我们会得到概率1,放缩放得太大了!所以我们不能取全0行来做这个放缩,这时我们得换一行),那么i=1∑nD1ixi=0可以移项得xk=−D1k1i=k∑D1ixi。如果除了xk以外都固定,只有xk随机,那我们发现等式右边就是一个固定的值,因此0和1中最多只有一个值是答案,所以概率小于等于1/2。这是我们的直观想法,把“局部固定”这件事严格化就恰好用到了全概率公式——对于每个固定的x1,⋯,xk−1,xk+1,⋯,xn,让它附带上xk=0和xk=1,就构成了样本空间的一个划分Ej。把xk=−D1k1i=k∑D1ixi记为事件G,于是Pr[G]=j∑Pr[G∩Ej],转化为条件概率=j∑Pr[G∣Ej]⋅Pr[Ej]。易知Pr[Ej]是常数,由于是划分求和后为1。而此时Pr[G∣Ej]这个条件概率就转化成我们显然能理解的了,有Pr[G∣Ej]≤1/2。因此j∑Pr[G∣Ej]⋅Pr[Ej]≤21j∑Pr[Ej],由于Ej是划分,j∑Pr[Ej]=1。综上Pr[G]≤1/2。因此Pr[Dx=0]≤1/2。
随机变量
随机变量是从样本空间到实数的映射Ω→R。比如“投两个骰子”的样本集有36个基本事件,此时“点数之和”是一个随机变量,(2,5)这个样本映射到实数7。
形式化地,我们用X=a来表示随机变量X的取值为a这一事件,它其实表示的是集合{s∈Ω∣X(s)=a}。从映射的角度,我们要在样本空间里找到a对应的原像的所有那些事件。更一般的,对于实数中的集合I类似的也可以定义事件X∈I,它对应着{s∈Ω∣X(s)∈I}。
我们特别地定义一个记号1[A]表示一个随机变量(函数),其中A是一个事件。这个随机变量作用在一个基本事件上,1[A](ω)=1当且仅当ω∈A,否则函数值为0。它直观上用来指示事件A是否会发生。
根据基本事件的定义,Pr(X=a)=X(s)=a,s∈Ω∑Pr(s)。由此可见我们其实是在用随机变量描述事件!因此我们可以把事件的性质应用到随机变量上。比如我们可以定义随机变量的独立性:X与Y独立当且仅当对于任意的a,b始终成立Pr((X=a)∩(Y=b))=Pr(X=a)⋅Pr(Y=b),记为X⊥Y。同样地,对于多个随机变量X1,⋯,Xn定义它们“互相独立”当且仅当对于任何I⊆[n],Pr[i∈I⋂(Xi=ai)]=i∈I∏Pr(Xi=ai)。
现在我们给出两两独立但不互相独立的反例:假设我们抛两次硬币,抛到正面取1,抛到反面取0。因此基本事件就是一个长度为2的01串。设随机变量X(ω)是从基本事件到其第一位的数字的映射(表示第一次的结果),Y(ω)是到第二位数字的映射,而定义Z(ω)=X(ω)⊕Y(ω)(异或),我们可以验证它们两两独立,但当I=[3]时Pr[X=0∩Y=0∩Z=0]=1/4,而Pr[X=0]=1/2,Pr[Y=0]=1/2,Pr[Z=0]=1/2,因此不成立,所以它们“两两独立”但是并不是“互相独立”的。
期望的线性性
期望是随机变量的一个重要特征:随机变量的期望定义为其概率的加权平均值。E[X]=i∑iPr(X=i)。它也可以按照基本事件来给出E[X]=ω∑X(ω)Pr(ω)。
一个最重要的定理称为“期望的线性性”,它指出当“期望”这样一个映射作用在若干个随机变量的“和”上时,等价于分别作用于每个随机变量再对它们求和:对于有限个随机变量X1,⋯,Xn始终满足E[i=1∑nXi]=i=1∑nE[Xi]。证明就是根据定义:E[X+Y]=i∑j∑(i+j)Pr((X=i)∩(Y=j)),提出固定的变量得到i∑ij∑Pr((X=i)∩(Y=j))+j∑ji∑Pr((X=i)∩(Y=j)),而后面的那个和式恰好就是全概率公式,因此直接可以化简得到i∑iPr(X=i)+j∑jPr(Y=j),这就是E[X]+E[Y]。这就证明了期望的线性性。从更高的角度看,期望的线性性其实来自于我们的求和的性质,这个性质在连续世界里就是“积分的线性性”。
重要的是,期望的线性性没有对随机变量提出任何的要求,即使两个随机变量不是“独立”的,它们的期望也始终符合“线性性”这一特征。这对我们计算期望提供极大的方便。
例如,我们求掷两次骰子的期望点数之和,一种方法是令X这个随机变量为点数之和,求得E=i=2∑12i⋅61⋅61⋅min{i−1,13−i}=7,而如果分别用X1和X2表示单独掷一次骰子的期望(实际上是只看某一次结果的随机变量),容易算出E0=i=1∑66i=7/2,两者相加就能直接得到7。
再例如,我们想知道一个给定的长度为k的01串在长度为n的01随机串里出现的期望次数。 利用期望的线性性,我们只需求出每个固定位置上出现该串的概率,再把它们相加即可。(这里我们用到了“某事件是否发生的期望就是它的概率”这一事实,严格化书写就是E[1[A]]=1⋅Pr[A]+0⋅Pr[A]=Pr[A])
再看一个有趣的例子。假设在我们在数轴上随机游走,从原点出发,每秒钟等概率向左或向右移动1,问t秒以后期望走到哪里。那么根据期望的线性性,每一步都期望移动距离0,因此t秒之后还是期望移动距离0,因为我们没有理由偏爱任何一个方向,我们不可能破坏对称性。事实上,对称性让这个问题变得无趣了,任意一个走到x的方案都会有一个走到−x的方案与它抵消。 一个更有意思的问题是问位移绝对值的期望,那样正负相消就不起效果了——我们从数学上来讨论他,记Zt=i=1∑tXi,其中Xi是表示第i步往右走还是往左走的随机变量。我们用一点数学技巧来解决这个问题,把∣Zt∣写作i=1∑t(∣Zi∣−∣Zi−1∣),于是E[∣Zt∣]=i=1∑tE[∣Zi∣−∣Zi−1∣],这是容易计算的,因为E[∣Zi∣−∣Zi−1∣]=2(∣Zi−1+1∣−∣Zi−1∣)+(∣Zi−1−1∣−∣Zi−1∣),画出函数图像注意到当Zi−1≥1或≤−1时函数值恒为0,只有在[−1,1]形成了一个45度角的突起。而由于任何一个Z都是整数,因此这个函数只有可能在Z=0时取到1,其余都为0。这样就得到E[∣Zi∣−∣Zi−1∣]=1[Zi−1=0]。因此E=i=1∑t1[Zi−1=0]。容易发现我们只有可能在偶数步的时候回到原点,因此也等价于E=i=0∑T/21[Z2i=0]。这是组合问题,我们必须走相同个数的左边和右边才会回到原点,因此写出E=i=0∑T/2(i2i)22i1,用Stirling公式代换阶乘,每一项可以近似为πi1,于是E≈1+i=1∑T/2πi1=O(T)(用积分来估计)。
当然,为了完善“线性性”这一概念,我们还应当证明E[cX]=cE[X]:这是容易的,E[cX]=i∑iPr(cX=i)=ci∑i/cPr(X=i/c)=cE[X]。
必须要指出的是,期望的线性性只对“有限个变量求和”成立。当变量有无穷个时它会出问题。我们再次从连续世界来看它,期望本身是一个积分,而对无穷多个变量求和本质上也是一个积分,因此期望的线性性就是一个“重积分能否化为累次积分的问题”,我们知道这是不一定能做到的。
条件期望
根据条件概率我们可以定义条件期望:E[X∣A]=i∑iPr(X=i∣A)或者E[X∣A]=ω∑X(ω)Pr(ω∣A)。
全期望公式
根据全概率公式Pr(A)=i=1∑nPr(A∩Bi)可以验证全期望公式E[X]=i=1∑nE[X∣Ai]Pr[Ai]。
根据定义证明:E[X]=i=1∑n(Pr[Ai]x∑xPr(X=x∣Ai))=i=1∑n(x∑xPr(X=x∣Ai)Pr[Ai]) =i=1∑nx∑xPr(X=x∩Ai)=x∑x(i=1∑nPr(X=x∩Ai))=x∑xPr(X=x)。
Markov不等式
概率与期望之间的大小关系有一个直接的必然联系:设随机变量X≥0,则∀a>0,Pr[X≥a]≤aE[X]。这称为Markov不等式。
根据期望定义,E[X]=x∑xPr[X=x]=x<a∑xPr[X=x]+x≥a∑xPr[X=x] ≥x≥a∑xPr[X=x]≥x≥a∑aPr[X=x]≥aPr[X≥a]。证毕。
二项分布
如果事件表示某件事的成功与否,成功概率为p,变量X=1,否则失败,变量X=0,那么这个变量X就称为伯努利随机变量。独立重复n次,记成功的次数为Y,那么称Y为二项随机变量,满足二项分布:
Pr(Y=i)=(in)pi(1−p)n−i
由此可以求出二项随机变量的期望:
E[Y]=i=0∑ni(in)pi(1−p)n−i=i=1∑n(i−1)!(n−i)!n!⋅pi(1−p)n−i=npi=1∑n(i−1)!(n−i)!(n−1)!pi−1(1−p)n−i=np(p+1−p)n−1=np
事实上还有更简单的证明方法,二项随机变量Y可以被拆解成n个基本事件的随机变量Xi,其中Xi表示第i次实验是否成功。于是就有关系Y=i=1∑nXi。两边同时取期望,E[Y]=E[i=1∑nXi],根据期望的线性性,=i=1∑nE[Xi]=i=1∑np=np。
几何分布
二项分布可以理解为丢n次硬币统计正面朝上的次数,那么几何分布就可以相应地理解为不停丢直到第一次出现正面的次数。因此Pr(X=n)=(1−p)n−1p。
如何算几何分布的期望?从直觉上看,丢一个1/10概率出正面的硬币大约需要10次,因此我们可以期望E[X]=p1。它的计算就是暴力求解等差乘等比数列的和。