UA MATH636 信息论7 高斯信道简介
UA MATH636 信息論7 高斯信道簡介
- 微分熵的性質(zhì)
- Gaussian Channel簡介
微分熵的性質(zhì)
討論Gauss信道之前,先給出幾條微分熵的有用的性質(zhì):
性質(zhì)1-5與離散型隨機(jī)變量的熵的性質(zhì)一樣,證明也基本一樣。
性質(zhì)6-7給出了隨機(jī)變量的線性運(yùn)算的微分熵,性質(zhì)6比較簡單,因?yàn)?span id="ze8trgl8bvbq" class="katex--inline">X+cX+cX+c與XXX同分布。考慮性質(zhì)7:
首先定義Y=aXY=aXY=aX,則fY(y)=1∣a∣fX(ya)f_Y(y)=\frac{1}{|a|}f_X(\frac{y}{a})fY?(y)=∣a∣1?fX?(ay?)
根據(jù)定義計算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=?∫∣a∣1?fX?(ay?)ln∣a∣1?fX?(ay?)dy=ln∣a∣∫fY?(y)dy?∫fX?(ay?)lnfX?(ay?)day?=h(X)+ln∣a∣
性質(zhì)8是性質(zhì)6、7到多元的推廣。證明也比較簡單,概率論那個系列給出了多遠(yuǎn)隨機(jī)變量的變換的方法的,在上面的計算過程中套用一下就好。
注意到這個上界,它正好是N(0,Σ)N(0,\Sigma)N(0,Σ)的微分熵。證明一下這個性質(zhì):考慮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(Σ))
根據(jù)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) dx0≤D(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隨機(jī)變量的微分熵的方法,可以得到
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簡介
之前討論信道編碼的時候已經(jīng)給出了信道的含義與功能,現(xiàn)在假設(shè)噪聲信道的模型是:
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)
通常假設(shè)ZiZ_iZi?與XiX_iXi?是獨(dú)立的,ZiZ_iZi?又被稱為Gaussian Noise,因此這個模型定義的信道又叫Gaussian Channel。
關(guān)于信道的輸入,通常有兩種類型的約束,存在約束的高斯信道叫input constrained Gaussian channel:
假設(shè)一串隨機(jī)信號表示為{X1,?,Xn}\{X_1,\cdots,X_n\}{X1?,?,Xn?},則這串隨機(jī)信號的平均能量一般定義為∑E[Xi2]/n\sum E[X_i^2]/n∑E[Xi2?]/n,假設(shè)PPP為能量上界,則
∑i=1nE[Xi2]n≤P\frac{\sum_{i=1}^n E[X_i^2]}{n} \le Pn∑i=1n?E[Xi2?]?≤P
因?yàn)殡S機(jī)信號獨(dú)立同分布,左邊這一項(xiàng)其實(shí)就是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):EX2≤Pmax?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(Y∣X)=h(Y)?h(X+Z∣X)=h(Y)?h(Z∣X)=h(Y)?h(Z)=h(Y)?21?ln(2πeN)
根據(jù)性質(zhì)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+2EXEZ≤P+N
這里用到了約束EX2≤PEX^2 \le PEX2≤P。所以
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))
當(dāng)且僅當(dāng)Y~N(0,P+N)Y \sim N(0,P+N)Y~N(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?)
當(dāng)且僅當(dāng)X~N(0,P)X \sim N(0,P)X~N(0,P)時取等。
總結(jié)
以上是生活随笔為你收集整理的UA MATH636 信息论7 高斯信道简介的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH564 概率论V 中心极限
- 下一篇: UA MATH565C 随机微分方程I