UA MATH564 概率论 计算至少有一个发生的概率:Waring公式
UA MATH564 概率論 計算至少有一個發生的概率:Waring公式
- 基于組合學的證明方法
- 基于概率論的證明方法
在推出Poincare公式后,計算nnn個事件中至少有一個事件發生的概率的問題就徹底解決了。但數學家們又開始糾結了,Poincare公式只能夠計算至少有一個發生的概率,那有沒有能計算恰好有mmm個事件同時發生的概率的方法呢?
假設用BmB_mBm?表示A1,?,AnA_1,\cdots,A_nA1?,?,An?中正好有mmm個發生的事件,則
P(Bm)=∑k=mn(?1)k?mCkmSkP(B_m) = \sum_{k=m}^n (-1)^{k-m}C_k^mS_kP(Bm?)=k=m∑n?(?1)k?mCkm?Sk?
這個公式叫做Waring公式。其中SkS_kSk?就是上一講定義的記號Sm=∑1≤i1<i2<?<im≤nP(Ai1∩Ai2∩?∩Aim)S_m = \sum_{1 \le i_1 < i_2 < \cdots < i_m \le n} P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_m})Sm?=1≤i1?<i2?<?<im?≤n∑?P(Ai1??∩Ai2??∩?∩Aim??)
基于組合學的證明方法
Feller的Introduction to Probability and its application中給出了基于組合學的證明方法。用反證法。我們要證明的公式等價于:假設EEE是一個恰好包含在lll個事件中的子事件,當且僅當l=ml=ml=m時,P(E)P(E)P(E)才是P(Bm)P(B_m)P(Bm?)的一部分。顯然這時P(E)P(E)P(E)也是S1,S2,?,SlS_1,S_2,\cdots,S_lS1?,S2?,?,Sl?的一部分,并且P(E)P(E)P(E)只會出現在SlS_lSl?的一項中。如果l<ml <ml<m,則P(E)P(E)P(E)不是P(Bm)P(B_m)P(Bm?)的一部分。如果l>ml>ml>m,則P(E)P(E)P(E)會出現在Sk,k≥mS_k,k\ge mSk?,k≥m的ClkC_l^kClk?項中,因此對Sk,k≥mS_k,k\ge mSk?,k≥m的貢獻是ClkCkmC_l^k C_k^mClk?Ckm?根據容斥原理,P(E)P(E)P(E)對P(Bm)P(B_m)P(Bm?)的貢獻為
ClmCmm?Clm+1Cm+1m+Clm+2Cm+2m+1??=Clm(Cl?m0?Cl?m1+?)=Clm(1?1)n?m=0C_l^mC_m^m - C_l^{m+1}C_{m+1}^m + C_l^{m+2}C_{m+2}^{m+1} - \cdots \\ = C_l^m(C_{l-m}^0 - C_{l-m}^1 + \cdots) = C_l^m(1-1)^{n-m} = 0Clm?Cmm??Clm+1?Cm+1m?+Clm+2?Cm+2m+1???=Clm?(Cl?m0??Cl?m1?+?)=Clm?(1?1)n?m=0
這就說明了當且僅當l=ml=ml=m時,P(E)P(E)P(E)才是P(Bm)P(B_m)P(Bm?)的一部分,因此上述公式成立。
基于概率論的證明方法
施利亞耶夫的概率論習題集提到了基于事件的Indicator function與期望的證明方法。我們記恰好發生的mmm個事件為i1,i2,?,imi_1,i_2,\cdots,i_mi1?,i2?,?,im?,記沒有發生的事件為j1,j2,?,jn?mj_1,j_2,\cdots,j_{n-m}j1?,j2?,?,jn?m?,記隨機變量X=IAX = I_AX=IA?,則恰好有mmm個事件發生的隨機變量是
∑Xi1Xi2?Xim(1?Xj1)?(1?Xjn?m)\sum X_{i_1}X_{i_2}\cdots X_{i_m} (1-X_{j_1})\cdots (1-X_{j_{n-m}})∑Xi1??Xi2???Xim??(1?Xj1??)?(1?Xjn?m??)
上式是對任意可能的i1,?,imi_1,\cdots,i_mi1?,?,im?的選取的求和。事件的Indicator function的期望等于事件發生的概率,因此
P(Bm)=E∑Xi1Xi2?Xim(1?Xj1)?(1?Xjn?m)P(B_m) = E\sum X_{i_1}X_{i_2}\cdots X_{i_m} (1-X_{j_1})\cdots (1-X_{j_{n-m}})P(Bm?)=E∑Xi1??Xi2???Xim??(1?Xj1??)?(1?Xjn?m??)
注意(1?Xj1)?(1?Xjn?m)(1-X_{j_1})\cdots (1-X_{j_{n-m}})(1?Xj1??)?(1?Xjn?m??)的展開式中一共有2n?m2^{n-m}2n?m項,其中一項為111,因此概率中有一項為
E∑Xi1Xi2?Xim=∑P(Ai1∩Ai2∩?∩Aim)E\sum X_{i_1}X_{i_2}\cdots X_{i_m} = \sum P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_m})E∑Xi1??Xi2???Xim??=∑P(Ai1??∩Ai2??∩?∩Aim??)
其余展開式期望的計算與此類似,事件的Indicator function乘積的期望等于事件的交發生的概率。下面研究(1?Xj1)?(1?Xjn?m)(1-X_{j_1})\cdots (1-X_{j_{n-m}})(1?Xj1??)?(1?Xjn?m??)的展開式。
引理1 用示性函數表示的容斥原理為:Xi=IAiX_i = I_{A_i}Xi?=IAi??
I?i=1nAi=1?∏i=1n(1?Xi)=∑i=1nXi?∑1≤i1<i2≤nXi1Xi2+?+(?1)m+1∑1≤i1<?<im≤nXi1?Xim?+(?1)nX1?XnI_{\bigcup_{i=1}^n A_i} =1 - \prod_{i=1}^n (1-X_i) = \sum_{i=1}^n X_i - \sum_{1 \le i_1 < i_2 \le n} X_{i_1}X_{i_2} + \cdots \\+ (-1)^{m+1}\sum_{1 \le i_1 < \cdots < i_m \le n} X_{i_1}\cdots X_{i_m} \cdots +(-1)^n X_1\cdots X_nI?i=1n?Ai??=1?i=1∏n?(1?Xi?)=i=1∑n?Xi??1≤i1?<i2?≤n∑?Xi1??Xi2??+?+(?1)m+11≤i1?<?<im?≤n∑?Xi1???Xim???+(?1)nX1??Xn?
基于這個引理把(1?Xj1)?(1?Xjn?m)(1-X_{j_1})\cdots (1-X_{j_{n-m}})(1?Xj1??)?(1?Xjn?m??)展開,然后把期望改寫成概率就可以得到Waring公式了。
總結
以上是生活随笔為你收集整理的UA MATH564 概率论 计算至少有一个发生的概率:Waring公式的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: html 360登录自动填写,汇总:如何
- 下一篇: 如何永久删除360SANDBOX?问题解