日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

UA MATH636 信息论5 信道编码定理

發布時間:2025/4/14 编程问答 18 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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(YX),其中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而言,它有下面三條性質

  • ∑j=1lyPij=1\sum_{j=1}^{l_y} P_{ij} = 1j=1ly??Pij?=1,這條性質就是概率的歸一性
  • ?i,i′\forall i,i'?i,i, ?σ∈Ply\exists \sigma \in P_{l_y}?σPly??, Pi′.=σ(Pi.)P_{i'.} = \sigma(P_{i.})Pi.?=σ(Pi.?)
  • ?j,j′\forall j,j'?j,j, ?τ∈Plx\exists \tau \in P_{l_x}?τPlx??, P.j=τ(P.j)P_{.j} = \tau(P_{.j})P.j?=τ(P.j?),性質2和3定義了對稱性
  • 其中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(YX)
    先計算H(Y∣X)H(Y|X)H(YX),根據性質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(YX)=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)直接看出的信道容量的兩個性質:

  • 0≤C≤min?(log?2lx,log?2ly)0\le C \le \min (\log_2 l_x, \log_2 l_y)0Cmin(log2?lx?,log2?ly?)
  • CCC具有唯一性
  • 事實上找信道容量并不是一個容易的事情,對于更復雜的信道以及網絡信道,找信道容量一般是要做比較復雜的優化的,證明完這個定理后會介紹計算信道容量的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:WXnnn表示使用的信道數目。則信源編碼的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}XD
    信道可以用p(Y∣X)p(Y|X)p(YX)表示,不妨假設XXXYYYnnn維字符,即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(YX)=i=1n?p(YiXi)
    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)?=iX=Xn(i))=YYn?p(YXn(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=1M?λ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 CRC。證明信道編碼定理需要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)的一個集合,并且滿足以下性質:

  • ∣?log?2p(xn)n?H(X)∣<?|-\frac{\log_2 p(x^n)}{n} - H(X)|<\epsilon?nlog2?p(xn)??H(X)<?,Marginal typical on x
  • ∣?log?2p(yn)n?H(Y)∣<?|-\frac{\log_2 p(y^n)}{n} - H(Y)|<\epsilon?nlog2?p(yn)??H(Y)<?,Marginal typical on y
  • ∣?log?2p(xn,yn)n?H(X,Y)∣<?|-\frac{\log_2 p(x^n,y^n)}{n} - H(X,Y)|<\epsilon?nlog2?p(xn,yn)??H(X,Y)<?,Joint typical on (x,y)
  • 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((xn,yn)∈A?(n))→1,asn→∞P((x^n,y^n) \in A_{\epsilon}^{(n)}) \to 1, as\ n \to \inftyP((xn,yn)A?(n)?)1,as?n
  • ∣A?(n))∣≤2n(H(X,Y)+?)|A_{\epsilon}^{(n)})| \le 2^{n(H(X,Y)+\epsilon)}A?(n)?)2n(H(X,Y)+?),當nnn足夠大時,∣A?(n))∣≥(1??)2n(H(X,Y)??)|A_{\epsilon}^{(n)})| \ge (1-\epsilon) 2^{n(H(X,Y)-\epsilon)}A?(n)?)(1??)2n(H(X,Y)??)
  • P((x~n,y~n)∈A?(n))P((\tilde{x}^n,\tilde{y}^n) \in A_{\epsilon}^{(n)})P((x~n,y~?n)A?(n)?)的上界為2?n(I(X;Y)?3?)2^{-n(I(X;Y)-3\epsilon)}2?n(I(X;Y)?3?),下界為(1??)2?n(I(X;Y)+3?)(1-\epsilon)2^{-n(I(X;Y)+3\epsilon)}(1??)2?n(I(X;Y)+3?),當nnn足夠大時下界成立。
  • 前兩個性質可以用之前的證明直接推廣就可以了,方法比較平凡。第三個性質之前沒有遇到過,這里給一個簡單證明:
    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 信道编码定理的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。