UA MATH564 概率论 计算至少有一个发生的概率:容斥原理与庞加莱公式
UA MATH564 概率論 計算至少有一個發(fā)生的概率:容斥原理與龐加萊公式
- 事件的并的Poincare公式
- 事件的交的Poincare公式
上一講介紹了P(?i=1nAi)P(\bigcup_{i=1}^{n} A_i)P(?i=1n?Ai?)的上界(Gunias不等式)和下界(鐘開萊-艾爾迪希不等式),這一講介紹準(zhǔn)確計算P(?i=1nAi)P(\bigcup_{i=1}^{n} A_i)P(?i=1n?Ai?)的方法——Poincare公式。
先考慮一個簡單情況,當(dāng)n=2n=2n=2時,A1∪A2A_1 \cup A_2A1?∪A2?可以很容易寫成不相交的三個子事件的并:
A1∪A2=(A1?A2)∪(A1∩A2)∪(A2?A1)A_1 \cup A_2 = (A_1-A_2) \cup (A_1\cap A_2) \cup (A_2 - A_1)A1?∪A2?=(A1??A2?)∪(A1?∩A2?)∪(A2??A1?)
根據(jù)概率的有限可加性,
P(A1∪A2)=P(A1?A2)+P(A1∩A2)+P(A2?A1)=P(A1)+P(A2)?P(A1∩A2)P(A_1 \cup A_2) = P(A_1-A_2) + P(A_1\cap A_2) + P(A_2 - A_1) \\ = P(A_1) + P(A_2)- P(A_1\cap A_2) P(A1?∪A2?)=P(A1??A2?)+P(A1?∩A2?)+P(A2??A1?)=P(A1?)+P(A2?)?P(A1?∩A2?)
這個表達(dá)式類似兩個集合的容斥原理。與容斥原理類似,我們可以嘗試把這個結(jié)果推廣到nnn個事件,獲得準(zhǔn)確計算P(?i=1nAi)P(\bigcup_{i=1}^{n} A_i)P(?i=1n?Ai?)的方法。
事件的并的Poincare公式
引入一個記號:
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??)
下面給出Poincare公式的形式:
P(?i=1nAi)=∑m=1n(?1)m+1SmP(\bigcup_{i=1}^{n} A_i) = \sum_{m=1}^n (-1)^{m+1}S_mP(i=1?n?Ai?)=m=1∑n?(?1)m+1Sm?
證明
其實早在Boole不等式的證明中,我們就已經(jīng)窺得Poincare公式的端倪。在Boole不等式的證明中,我們構(gòu)造了一個集列{Bi}i=1n\{B_i\}_{i=1}^n{Bi?}i=1n?,其中B1=A1,B2=A2?A1,?,Bn=An??i=1n?1AnB_1 = A_1,B_2 = A_2-A_1,\cdots,B_n = A_n - \bigcup_{i=1}^{n-1}A_nB1?=A1?,B2?=A2??A1?,?,Bn?=An???i=1n?1?An?,顯然Bi∩Bj=?,?i≠jB_i \cap B_j = \phi, \forall i \ne jBi?∩Bj?=?,?i?=j,并且
?i=1nBi=?i=1nAi\bigcup_{i=1}^n B_i = \bigcup_{i=1}^{n} A_ii=1?n?Bi?=i=1?n?Ai?
現(xiàn)在我們利用這個結(jié)構(gòu)做計算
P(?i=1nAi)=P(?i=1nBi)=P(?i=1n?1Bi)+P(Bn)=P(?i=1n?1Ai)+P(An??i=1n?1An)=P(?i=1n?1Ai)+P(An)?P(An∩?i=1n?1Ai)P(\bigcup_{i=1}^{n} A_i) = P(\bigcup_{i=1}^n B_i) = P(\bigcup_{i=1}^{n-1} B_i)+P(B_n) \\= P(\bigcup_{i=1}^{n-1} A_i) + P(A_n - \bigcup_{i=1}^{n-1}A_n) = P(\bigcup_{i=1}^{n-1} A_i) + P(A_n) - P(A_n \cap \bigcup_{i=1}^{n-1} A_i)P(i=1?n?Ai?)=P(i=1?n?Bi?)=P(i=1?n?1?Bi?)+P(Bn?)=P(i=1?n?1?Ai?)+P(An??i=1?n?1?An?)=P(i=1?n?1?Ai?)+P(An?)?P(An?∩i=1?n?1?Ai?)
根據(jù)這個遞推關(guān)系,
P(?i=1nAi)=∑i=1nP(Ai)?∑j=1n?1P(Aj+1∩?i=1jAi)P(\bigcup_{i=1}^{n} A_i) = \sum_{i=1}^n P(A_i) - \sum_{j=1}^{n-1} P(A_{j+1} \cap \bigcup_{i=1}^{j} A_i)P(i=1?n?Ai?)=i=1∑n?P(Ai?)?j=1∑n?1?P(Aj+1?∩i=1?j?Ai?)
下面考慮P(Aj+1∩?i=1jAi)P(A_{j+1} \cap \bigcup_{i=1}^{j} A_i)P(Aj+1?∩?i=1j?Ai?),同樣依據(jù)上面那個遞推關(guān)系:
P(Aj+1∩?i=1jAi)=P(?i=1jAi∩Aj+1)=P(Aj∩Aj+1)?P(Aj∩Aj+1∩?k=1j?1Ak)P(A_{j+1} \cap \bigcup_{i=1}^{j} A_i) = P(\bigcup_{i=1}^{j} A_i\cap A_{j+1}) = P(A_j \cap A_{j+1}) - P(A_{j}\cap A_{j+1} \cap \bigcup_{k=1}^{j-1} A_k)P(Aj+1?∩i=1?j?Ai?)=P(i=1?j?Ai?∩Aj+1?)=P(Aj?∩Aj+1?)?P(Aj?∩Aj+1?∩k=1?j?1?Ak?)
因此
P(?i=1nAi)=∑i=1nP(Ai)?∑1≤i<j≤nP(Ai∩Aj)+∑j=1n?1P(Aj∩Aj+1∩?k=1j?1Ak)P(\bigcup_{i=1}^{n} A_i) = \sum_{i=1}^n P(A_i) - \sum_{1 \le i < j \le n}P(A_i \cap A_j) + \sum_{j=1}^{n-1} P(A_{j}\cap A_{j+1} \cap \bigcup_{k=1}^{j-1} A_k)P(i=1?n?Ai?)=i=1∑n?P(Ai?)?1≤i<j≤n∑?P(Ai?∩Aj?)+j=1∑n?1?P(Aj?∩Aj+1?∩k=1?j?1?Ak?)
接下來就是對P(Aj∩Aj+1∩?k=1j?1Ak)P(A_{j}\cap A_{j+1} \cap \bigcup_{k=1}^{j-1} A_k)P(Aj?∩Aj+1?∩?k=1j?1?Ak?)用那個遞推公式,把所有三三相交的概率分離出來,一直做下去做到nnnnnn相交的概率即可!
證畢
事件的交的Poincare公式
基于事件的并的Poincare公式與de Morgan律,我們可以推出事件的交的Poincare公式。
P(?i=1nAi)=1?P(?i=1nAˉi)=1?∑m=1n(?1)m+1∑1≤i1<i2<?<im≤nP(Aˉi1∩Aˉi2∩?∩Aˉim)=1?∑m=1n(?1)m+1∑1≤i1<i2<?<im≤n[1?P(Ai1∪Ai2∪?∪Aim)]=1?∑m=1n(?1)m+1Cnm+∑m=1n(?1)m+1∑1≤i1<i2<?<im≤nP(Ai1∪Ai2∪?∪Aim)=∑m=1n(?1)m+1∑1≤i1<i2<?<im≤nP(Ai1∪Ai2∪?∪Aim)P(\bigcap_{i=1}^{n} A_i) = 1 -P(\bigcup_{i=1}^{n} \bar{A}_i) = 1-\sum_{m=1}^n (-1)^{m+1}\sum_{1 \le i_1 < i_2 < \cdots < i_m \le n} P(\bar{A}_{i_1} \cap \bar{A}_{i_2} \cap \cdots \cap \bar{A}_{i_m}) \\ = 1-\sum_{m=1}^n (-1)^{m+1}\sum_{1 \le i_1 < i_2 < \cdots < i_m \le n} [1-P(A_{i_1} \cup A_{i_2} \cup \cdots \cup A_{i_m})] \\ = 1-\sum_{m=1}^n (-1)^{m+1}C_n^m + \sum_{m=1}^n (-1)^{m+1}\sum_{1 \le i_1 < i_2 < \cdots < i_m \le n}P(A_{i_1} \cup A_{i_2} \cup \cdots \cup A_{i_m}) \\ = \sum_{m=1}^n (-1)^{m+1}\sum_{1 \le i_1 < i_2 < \cdots < i_m \le n}P(A_{i_1} \cup A_{i_2} \cup \cdots \cup A_{i_m})P(i=1?n?Ai?)=1?P(i=1?n?Aˉi?)=1?m=1∑n?(?1)m+11≤i1?<i2?<?<im?≤n∑?P(Aˉi1??∩Aˉi2??∩?∩Aˉim??)=1?m=1∑n?(?1)m+11≤i1?<i2?<?<im?≤n∑?[1?P(Ai1??∪Ai2??∪?∪Aim??)]=1?m=1∑n?(?1)m+1Cnm?+m=1∑n?(?1)m+11≤i1?<i2?<?<im?≤n∑?P(Ai1??∪Ai2??∪?∪Aim??)=m=1∑n?(?1)m+11≤i1?<i2?<?<im?≤n∑?P(Ai1??∪Ai2??∪?∪Aim??)
定義記號
S~m=∑1≤i1<i2<?<im≤nP(Ai1∪Ai2∪?∪Aim)\tilde{S}_m = \sum_{1 \le i_1 < i_2 < \cdots < i_m \le n}P(A_{i_1} \cup A_{i_2} \cup \cdots \cup A_{i_m})S~m?=1≤i1?<i2?<?<im?≤n∑?P(Ai1??∪Ai2??∪?∪Aim??)
事件的交的Poincare公式可以寫成
P(?i=1nAi)=∑m=1n(?1)m+1S~mP(\bigcap_{i=1}^{n} A_i) = \sum_{m=1}^n (-1)^{m+1}\tilde{S}_mP(i=1?n?Ai?)=m=1∑n?(?1)m+1S~m?
總結(jié)
以上是生活随笔為你收集整理的UA MATH564 概率论 计算至少有一个发生的概率:容斥原理与庞加莱公式的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH564 概率论 计算至少有
- 下一篇: UA MATH566 统计理论 QE练习