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

歡迎訪問 生活随笔!

生活随笔

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

编程问答

UA MATH636 信息论7 高斯信道简介

發布時間:2025/4/14 编程问答 28 豆豆
生活随笔 收集整理的這篇文章主要介紹了 UA MATH636 信息论7 高斯信道简介 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

UA MATH636 信息論7 高斯信道簡介

  • 微分熵的性質
  • Gaussian Channel簡介

微分熵的性質

討論Gauss信道之前,先給出幾條微分熵的有用的性質:

  • Gibbs不等式:D(p∣∣q)≥0D(p||q) \ge 0D(pq)0
  • 互信息的非負性:I(X;Y)≥0I(X;Y) \ge 0I(X;Y)0
  • Conditioning Reducing Entropy:h(X)≥h(X∣Y)h(X) \ge h(X|Y)h(X)h(XY)
  • Chain Rule:h(X,Y)=h(X)+h(X∣Y)h(X,Y) = h(X)+h(X|Y)h(X,Y)=h(X)+h(XY)
  • 根據3和4:h(X,Y)≤h(X)+h(Y)h(X,Y) \le h(X)+h(Y)h(X,Y)h(X)+h(Y)
  • 性質1-5與離散型隨機變量的熵的性質一樣,證明也基本一樣。

  • h(X+c)=h(X)h(X+c)=h(X)h(X+c)=h(X)ccc是一個常數
  • h(aX)=h(X)+ln?∣a∣h(aX)=h(X)+\ln |a|h(aX)=h(X)+lnaaaa是一個常數
  • 性質6-7給出了隨機變量的線性運算的微分熵,性質6比較簡單,因為X+cX+cX+cXXX同分布。考慮性質7:
    首先定義Y=aXY=aXY=aX,則fY(y)=1∣a∣fX(ya)f_Y(y)=\frac{1}{|a|}f_X(\frac{y}{a})fY?(y)=a1?fX?(ay?)
    根據定義計算h(Y)h(Y)h(Y)
    h(Y)=?∫fY(y)ln?fY(y)dy=?∫1∣a∣fX(ya)ln?1∣a∣fX(ya)dy=ln?∣a∣∫fY(y)dy?∫fX(ya)ln?fX(ya)dya=h(X)+ln?∣a∣h(Y) =- \int f_Y(y)\ln f_Y(y)dy = -\int \frac{1}{|a|}f_X(\frac{y}{a}) \ln \frac{1}{|a|}f_X(\frac{y}{a})dy \\ = \ln |a| \int f_Y(y)dy -\int f_X(\frac{y}{a}) \ln f_X(\frac{y}{a})d \frac{y}{a} =h(X)+\ln |a|h(Y)=?fY?(y)lnfY?(y)dy=?a1?fX?(ay?)lna1?fX?(ay?)dy=lnafY?(y)dy?fX?(ay?)lnfX?(ay?)day?=h(X)+lna

  • h(AX+B)=h(X)+ln?∣det?(A)∣h(AX+B)=h(X)+\ln |\det(A)|h(AX+B)=h(X)+lndet(A)
  • 性質8是性質6、7到多元的推廣。證明也比較簡單,概率論那個系列給出了多遠隨機變量的變換的方法的,在上面的計算過程中套用一下就好。

  • h(X)≤12ln?((2πe)ndet?(Σ))h(X) \le \frac{1}{2} \ln \left((2\pi e)^n \det(\Sigma)\right)h(X)21?ln((2πe)ndet(Σ)),其中EX=0,Cov(X)=ΣEX=0,Cov(X)=\SigmaEX=0,Cov(X)=Σ
  • 注意到這個上界,它正好是N(0,Σ)N(0,\Sigma)N(0,Σ)的微分熵。證明一下這個性質:考慮XXX的概率密度為g(x)g(x)g(x),記XG~N(0,Σ)X_G\sim N(0,\Sigma)XG?N(0,Σ),分布記為ΦG(x)\Phi_G(x)ΦG?(x),則
    ΦG(x)=(2π)?n/2det?(Σ)?1/2exp?(?12xTΣ?1x)ln?ΦG(x)=?12xTΣ?1x?12ln?((2π)ndet?(Σ))\Phi_G(x) = (2\pi)^{-n/2} \det(\Sigma)^{-1/2} \exp \left( -\frac{1}{2} x^T \Sigma^{-1} x\right) \\ \ln \Phi_G(x) = -\frac{1}{2} x^T \Sigma^{-1} x - \frac{1}{2} \ln ((2\pi )^n \det(\Sigma))ΦG?(x)=(2π)?n/2det(Σ)?1/2exp(?21?xTΣ?1x)lnΦG?(x)=?21?xTΣ?1x?21?ln((2π)ndet(Σ))
    根據Gibbs不等式,
    0≤D(g∣∣ΦG)=∫g(x)ln?g(x)ΦG(x)dx=?h(X)?∫g(x)ln?ΦG(x)dx=?h(X)?∫g(x)(?12xTΣ?1x?12ln?((2π)ndet?(Σ)))dx0 \le D(g || \Phi_G) = \int g(x) \ln \frac{g(x)}{\Phi_G(x)} dx \\ = -h(X) - \int g(x) \ln \Phi_G(x)dx \\ = -h(X)-\int g(x) \left( -\frac{1}{2} x^T \Sigma^{-1} x - \frac{1}{2} \ln ((2\pi )^n \det(\Sigma))\right) dx0D(gΦG?)=g(x)lnΦG?(x)g(x)?dx=?h(X)?g(x)lnΦG?(x)dx=?h(X)?g(x)(?21?xTΣ?1x?21?ln((2π)ndet(Σ)))dx
    其中
    ∫12ln?((2π)ndet?(Σ))g(x)dx=12ln?((2π)ndet?(Σ))∫12xTΣ?1xg(x)dx=E[12XTΣ?1X]\int \frac{1}{2} \ln ((2\pi )^n \det(\Sigma))g(x)dx = \frac{1}{2} \ln ((2\pi )^n \det(\Sigma)) \\ \int \frac{1}{2} x^T \Sigma^{-1} x g(x) dx=E[\frac{1}{2} X^T \Sigma^{-1} X]21?ln((2π)ndet(Σ))g(x)dx=21?ln((2π)ndet(Σ))21?xTΣ?1xg(x)dx=E[21?XTΣ?1X]
    用上一講計算Gauss隨機變量的微分熵的方法,可以得到
    0≤?h(X)?12ln?((2πe)ndet?(Σ))0 \le -h(X)-\frac{1}{2} \ln ((2\pi e)^n \det(\Sigma)) 0?h(X)?21?ln((2πe)ndet(Σ))
    因此上界成立。

    Gaussian Channel簡介

    之前討論信道編碼的時候已經給出了信道的含義與功能,現在假設噪聲信道的模型是:
    Yi=Xi+Zi,Zi~N(0,N)Y_i = X_i + Z_i,Z_i \sim \mathcal{N}(0,N)Yi?=Xi?+Zi?,Zi?N(0,N)
    通常假設ZiZ_iZi?XiX_iXi?是獨立的,ZiZ_iZi?又被稱為Gaussian Noise,因此這個模型定義的信道又叫Gaussian Channel。
    關于信道的輸入,通常有兩種類型的約束,存在約束的高斯信道叫input constrained Gaussian channel:

  • Power constrained input
    假設一串隨機信號表示為{X1,?,Xn}\{X_1,\cdots,X_n\}{X1?,?,Xn?},則這串隨機信號的平均能量一般定義為∑E[Xi2]/n\sum E[X_i^2]/nE[Xi2?]/n,假設PPP為能量上界,則
    ∑i=1nE[Xi2]n≤P\frac{\sum_{i=1}^n E[X_i^2]}{n} \le Pni=1n?E[Xi2?]?P
    因為隨機信號獨立同分布,左邊這一項其實就是EX2EX^2EX2
  • 定理1 這個約束下高斯信道的信道容量是
    C=12ln?(1+PN)C = \frac{1}{2}\ln \left(1 + \frac{P}{N} \right)C=21?ln(1+NP?)
    其中P/NP/NP/N是信噪比(signal to noise ratio, SNR)。
    證明
    C=max?f(x):EX2≤PI(X;Y)C = \max_{f(x):EX^2 \le P} I(X;Y)C=f(x):EX2Pmax?I(X;Y)
    計算
    I(X;Y)=h(Y)?h(Y∣X)=h(Y)?h(X+Z∣X)=h(Y)?h(Z∣X)=h(Y)?h(Z)=h(Y)?12ln?(2πeN)I(X;Y)=h(Y)-h(Y|X)=h(Y)-h(X+Z|X) \\ = h(Y)-h(Z|X)=h(Y)-h(Z)=h(Y) - \frac{1}{2} \ln (2 \pi e N)I(X;Y)=h(Y)?h(YX)=h(Y)?h(X+ZX)=h(Y)?h(ZX)=h(Y)?h(Z)=h(Y)?21?ln(2πeN)
    根據性質9,
    h(Y)≤12ln?(2πeVar(Y))h(Y) \le \frac{1}{2} \ln (2 \pi e Var(Y))h(Y)21?ln(2πeVar(Y))
    計算
    E[Y2]=E[(X+Z)2]=EX2+EZ2+2EXEZ≤P+NE[Y^2] = E[(X+Z)^2] = EX^2 + EZ^2 + 2EXEZ \le P+ NE[Y2]=E[(X+Z)2]=EX2+EZ2+2EXEZP+N
    這里用到了約束EX2≤PEX^2 \le PEX2P。所以
    h(Y)≤12ln?(2πe(P+N))h(Y) \le \frac{1}{2} \ln (2 \pi e (P+N))h(Y)21?ln(2πe(P+N))
    當且僅當Y~N(0,P+N)Y \sim N(0,P+N)YN(0,P+N)是取等。從而
    I(X;Y)≤12ln?(2πe(P+N))?12ln?(2πeN)=12ln?(P+NN)=12ln?(1+PN)I(X;Y) \le \frac{1}{2} \ln (2 \pi e (P+N))-\frac{1}{2} \ln (2 \pi e N) \\ = \frac{1}{2} \ln \left( \frac{P+N}{N} \right) = \frac{1}{2}\ln \left(1 + \frac{P}{N} \right)I(X;Y)21?ln(2πe(P+N))?21?ln(2πeN)=21?ln(NP+N?)=21?ln(1+NP?)
    當且僅當X~N(0,P)X \sim N(0,P)XN(0,P)時取等。

    總結

    以上是生活随笔為你收集整理的UA MATH636 信息论7 高斯信道简介的全部內容,希望文章能夠幫你解決所遇到的問題。

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