UA MATH636 信息论5 信道编码定理
UA MATH636 信息論5 信道編碼定理
- 信道編碼問題
- 信道容量的正式定義
- 信道編碼定理
- Joint Typical Set
- Joint AEP
上一篇簡介里面介紹了通訊的過程,并用下面的流程圖來表示了這個過程,并且引入了一些基本概念。這一講先對上一篇的例1做一個推廣,然后介紹一下分析信道編碼的工具。
WXYVSourceEncoderNoisy ChannelDecoderReceiver例1續:Symmetric Channels
假設矩陣PPP表示p(Y∣X)p(Y|X)p(Y∣X),其中Pij=P(Yj∣Xi)P_{ij} = P(Y_j|X_i)Pij?=P(Yj?∣Xi?),定義Pi.P_{i.}Pi.?表示矩陣PPP的第iii行,P.jP_{.j}P.j?表示矩陣PPP的第jjj列。已知Y={Y1,Y2,?,YlY}Y=\{Y_1,Y_2,\cdots,Y_{l_Y}\}Y={Y1?,Y2?,?,YlY??},X={X1,X2,?,XlX}X=\{X_1,X_2,\cdots,X_{l_X}\}X={X1?,X2?,?,XlX??}。之所以提出這類信道是因為它的容量比較簡單好算。對于矩陣PPP而言,它有下面三條性質
其中PlyP_{l_y}Ply??與PlxP_{l_x}Plx??分別為lyl_yly?與lxl_xlx?階的置換運算的集合。這類信道的容量為
C=log?2ly?H(Pi.)C = \log_2 l_y - H(P_{i.})C=log2?ly??H(Pi.?)
其中H(Pi.)H(P_{i.})H(Pi.?)表示任一行的熵。
證明:
I(X;Y)=H(Y)?H(Y∣X)I(X;Y) = H(Y) - H(Y|X)I(X;Y)=H(Y)?H(Y∣X)
先計算H(Y∣X)H(Y|X)H(Y∣X),根據性質2
H(Y∣X)=∑ip(X=Xi)H(Pi.)=H(Pi.)H(Y|X) = \sum_i p(X=X_i)H(P_{i.})=H(P_{i.})H(Y∣X)=i∑?p(X=Xi?)H(Pi.?)=H(Pi.?)
比較直觀地,根據熵的有界性,H(Y)≤log?2lyH(Y) \le \log_2 l_yH(Y)≤log2?ly?,當且僅當YYY為均勻分布時成立,要讓YYY為均勻分布,根據性質1,只需XXX為均勻分布即可。
注意到這個證明不需要性質3,如果定義weakly symmetric channel(只滿足性質1,2),則同樣有
C=log?2ly?H(Pi.)C = \log_2 l_y - H(P_{i.})C=log2?ly??H(Pi.?)
在正式定義信道編碼問題之前,先給出可以從I(X;Y)I(X;Y)I(X;Y)直接看出的信道容量的兩個性質:
事實上找信道容量并不是一個容易的事情,對于更復雜的信道以及網絡信道,找信道容量一般是要做比較復雜的優化的,證明完這個定理后會介紹計算信道容量的Blahut-Arimoto算法,并比較這個算法與EM算法的區別。
信道編碼問題
信道容量的正式定義
假設W∈{1,2,?,M}W \in \{1,2,\cdots,M\}W∈{1,2,?,M},即信源可能發送MMM種可能的信息,編碼器記作Xn:W→X\mathcal{X}^n:W\to XXn:W→X,nnn表示使用的信道數目。則信源編碼的codebook可以記作D={Xn(1),?,Xn(M)}\mathcal{D} =\{\mathcal{X}^n(1),\cdots,\mathcal{X}^n(M)\}D={Xn(1),?,Xn(M)},從而X∈DX \in \mathcal{D}X∈D。
信道可以用p(Y∣X)p(Y|X)p(Y∣X)表示,不妨假設XXX與YYY是nnn維字符,即X=(X1,?,Xn)X = (X^1,\cdots,X^n)X=(X1,?,Xn),Y=(Y1,?,Yn)Y=(Y^1,\cdots,Y^n)Y=(Y1,?,Yn),假設信道是無記憶的,即
p(Y∣X)=∏i=1np(Yi∣Xi)p(Y|X) = \prod_{i=1}^n p(Y^i|X^i)p(Y∣X)=i=1∏n?p(Yi∣Xi)
記YYY的值域為Yn\mathcal{Y}^nYn。解碼器可以表示為g:Yn→{1,2,?,M}g:\mathcal{Y}^n \to \{1,2,\cdots,M\}g:Yn→{1,2,?,M},ggg為解碼函數。信道編碼想要研究的對象之一就是這組通訊編碼(可以記為(M,n)(M,n)(M,n)-code)的錯誤率λi\lambda_iλi?:
λi=p(g(Y)≠i∣X=Xn(i))=∑Y∈Ynp(Y∣Xn(i))I(g(Y)≠i)\lambda_i = p(g(Y) \ne i | X = \mathcal{X}^n(i)) \\ = \sum_{Y \in \mathcal{Y}^n} p(Y|\mathcal{X}^n(i))I(g(Y) \ne i)λi?=p(g(Y)?=i∣X=Xn(i))=Y∈Yn∑?p(Y∣Xn(i))I(g(Y)?=i)
由此可以定義最大錯誤率
λ(n)=max?iλi\lambda^{(n)} = \max_i \lambda_iλ(n)=imax?λi?
以及平均錯誤率
pe(n)=1M∑i=1Mλip_e^{(n)} = \frac{1}{M} \sum_{i=1}^M \lambda_ipe(n)?=M1?i=1∑M?λi?
定義(M,n)(M,n)(M,n)-code的傳輸率為
R=log?2MnR = \frac{\log_2 M}{n}R=nlog2?M?
通常給定一個傳輸率,使用nnn個channel可以傳輸的信號為
M=ceiling(2nR)M = ceiling(2^{nR})M=ceiling(2nR)
實際上是M=2nRM = 2^{nR}M=2nR,但一般要取整數。如果λ(n)→0\lambda^{(n)} \to 0λ(n)→0,稱傳輸率RRR對于(M,n)(M,n)(M,n)-code為可實現的(achievable)。現在可以給出信道容量的正式定義:
max?{R:Risachievable}\max\{R:R\ is\ achievable\}max{R:R?is?achievable}
信道編碼定理
信道編碼定理是通訊理論中一個非常重要的定理,但它的表述非常簡單:所有小于CCC的傳輸率是可實現的。這里的CCC就是我們之前定義的
C=max?p(X)I(X;Y)C = \max_{p(X)} I(X;Y)C=p(X)max?I(X;Y)
因此這個表述指的是
C=max?p(X)I(X;Y)=max?{R:Risachievable}C = \max_{p(X)} I(X;Y) = \max\{R:R\ is\ achievable\}C=p(X)max?I(X;Y)=max{R:R?is?achievable}
這個表述的逆命題同樣成立,對于任何(M,n)(M,n)(M,n)-code,如果λ(n)→0\lambda^{(n)} \to 0λ(n)→0,則R≤CR \le CR≤C。證明信道編碼定理需要random coding的思想以及Joint typical set,證明逆命題需要Fano不等式以及數據處理不等式,在證明之前先介紹一下需要的工具。
Joint Typical Set
從聯合分布p(x,y)p(x,y)p(x,y)中生成隨機樣本(xn,yn)={xi,yi}i=1n(x^n,y^n) = \{x_i,y_i\}_{i=1}^n(xn,yn)={xi?,yi?}i=1n?,記Joint Typical Set為A?(n)A_{\epsilon}^{(n)}A?(n)?,它是(xn,yn)(x^n,y^n)(xn,yn)的一個集合,并且滿足以下性質:
Joint AEP
之前講信源編碼的時候就講過AEP了,這里推廣到joint case:記(x~n,y~n)(\tilde{x}^n,\tilde{y}^n)(x~n,y~?n)是從p(x)p(y)p(x)p(y)p(x)p(y)中生成的隨機樣本
前兩個性質可以用之前的證明直接推廣就可以了,方法比較平凡。第三個性質之前沒有遇到過,這里給一個簡單證明:
P((x~n,y~n)∈A?(n))=∑(x~n,y~n)∈A?(n)p(x~n,y~n)=∑(x~n,y~n)∈A?(n)p(x~n)p(y~n)≤∑(x~n,y~n)∈A?(n)2?n(H(X)??)2?n(H(Y)??)=∣A?(n)∣2?n(H(X)+H(Y)?2?)P((\tilde{x}^n,\tilde{y}^n) \in A_{\epsilon}^{(n)}) = \sum_{(\tilde{x}^n,\tilde{y}^n) \in A_{\epsilon}^{(n)}} p(\tilde{x}^n,\tilde{y}^n) = \sum_{(\tilde{x}^n,\tilde{y}^n) \in A_{\epsilon}^{(n)}} p(\tilde{x}^n)p(\tilde{y}^n) \\ \le \sum_{(\tilde{x}^n,\tilde{y}^n) \in A_{\epsilon}^{(n)}} 2^{-n(H(X)-\epsilon)} 2^{-n(H(Y)-\epsilon)} = |A_{\epsilon}^{(n)}| 2^{-n(H(X)+H(Y)-2\epsilon)} P((x~n,y~?n)∈A?(n)?)=(x~n,y~?n)∈A?(n)?∑?p(x~n,y~?n)=(x~n,y~?n)∈A?(n)?∑?p(x~n)p(y~?n)≤(x~n,y~?n)∈A?(n)?∑?2?n(H(X)??)2?n(H(Y)??)=∣A?(n)?∣2?n(H(X)+H(Y)?2?)
將性質2帶入,可以得到上界。類似的過程可以證明下界。
總結
以上是生活随笔為你收集整理的UA MATH636 信息论5 信道编码定理的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH636 信息论5 信道编码
- 下一篇: UA MATH636 信息论5 信道编码