Bootstrap

信道容量计算题

目录

T1

(二进无噪信道)设 X = Y = { 0 , 1 } , Q = I 2 X=\mathcal{Y}=\{0,1\},Q=I_2 X=Y={0,1},Q=I2 ,试求信道容量。

解:

I ( p ; Q ) = p ( 0 ) [ Q ( 0 ∣ 0 ) log ⁡ Q ( 0 ∣ 0 ) p ( 0 ) Q ( 0 ∣ 0 ) + p ( 1 ) Q ( 0 ∣ 1 ) + Q ( 1 ∣ 0 ) log ⁡ Q ( 1 ∣ 0 ) p ( 0 ) Q ( 1 ∣ 0 ) + p ( 1 ) Q ( 1 ∣ 1 ) ] + p ( 1 ) [ Q ( 0 ∣ 1 ) log ⁡ Q ( 0 ∣ 1 ) p ( 0 ) Q ( 0 ∣ 0 ) + p ( 1 ) Q ( 0 ∣ 1 ) + Q ( 1 ∣ 1 ) log ⁡ Q ( 1 ∣ 0 ) p ( 0 ) Q ( 1 ∣ 0 ) + p ( 1 ) Q ( 1 ∣ 1 ) ] = p ( 0 ) log ⁡ 1 p ( 0 ) + p ( 1 ) log ⁡ 1 p ( 1 ) = H ( X ) \begin{aligned}&I(p;Q)=p(0)[Q(0|0)\log\frac{Q(0|0)}{p(0)Q(0|0)+p(1)Q(0|1)}\\&+Q(1|0)\log\frac{Q(1|0)}{p(0)Q(1|0)+p(1)Q(1|1)}]\\&+p(1)[Q(0|1)\log\frac{Q(0|1)}{p(0)Q(0|0)+p(1)Q(0|1)}\\&+Q(1|1)\log\frac{Q(1|0)}{p(0)Q(1|0)+p(1)Q(1|1)}]\\&=p(0)\log\frac1{p(0)}+p(1)\log\frac1{p(1)}=H(X)\end{aligned} I(p;Q)=p(0)[Q(0∣0)logp(0)Q(0∣0)+p(1)Q(0∣1)Q(0∣0)+Q(1∣0)logp(0)Q(1∣0)+p(1)Q(1∣1)Q(1∣0)]+p(1)[Q(0∣1)logp(0)Q(0∣0)+p(1)Q(0∣1)Q(0∣1)+Q(1∣1)logp(0)Q(1∣0)+p(1)Q(1∣1)Q(1∣0)]=p(0)logp(0)1+p(1)logp(1)1=H(X)
因此 C = max ⁡ n ( x ) H ( X ) = log ⁡ ∣ ∣ X ∣ ∣ = 1 C=\max_{n(x)}H(X)=\log||\mathcal{X}||=1 C=maxn(x)H(X)=log∣∣X∣∣=1bit,当且仅当 p p p为均匀分布时取得。



T2

(对称信道)设 X = Y = { 0 , 1 } \mathcal{X} = \mathcal{Y} = \{ 0, 1\} X=Y={0,1}, Q = [ 1 − ε ε ε 1 − ε ] Q= \begin{bmatrix} 1- \varepsilon & \varepsilon \\ \varepsilon & 1- \varepsilon \end{bmatrix} Q=[1εεε1ε],求信道容量。
解:
计算可得

( p ; Q ) = − [ p ( 0 ) ( 1 − ε ) + p ( 1 ) ε ] log ⁡ [ p ( 0 ) ( 1 − ε ) + p ( 1 ) ε ] − [ p ( 0 ) ε + p ( 1 ) ( 1 − ε ) ] log ⁡ [ p ( 0 ) ε + p ( 1 ) ( 1 − ε ) ] − ε log ⁡ ε − ( 1 − ε ) log ⁡ ( 1 − ε ) = H ( Y ) − h ( ε ) ≤ log ⁡ ∣ ∣ Y ∣ ∣ − h ( ε ) = 1 − h ( ε ) . \begin{aligned}&(p;Q)=-[p(0)(1-\varepsilon)+p(1)\varepsilon]\log[p(0)(1-\varepsilon)+p(1)\varepsilon]\\&-[p(0)\varepsilon+p(1)(1-\varepsilon)]\log[p(0)\varepsilon+p(1)(1-\varepsilon)]\\&-\varepsilon\log\varepsilon-(1-\varepsilon)\log(1-\varepsilon)\\&=H(Y)-h(\varepsilon)\leq\log||\mathcal{Y}||-h(\varepsilon)=1-h(\varepsilon).\end{aligned} (p;Q)=[p(0)(1ε)+p(1)ε]log[p(0)(1ε)+p(1)ε][p(0)ε+p(1)(1ε)]log[p(0)ε+p(1)(1ε)]εlogε(1ε)log(1ε)=H(Y)h(ε)log∣∣Y∣∣h(ε)=1h(ε).

其中 h h h为第一章中的二进熵函数。显然等号成立当且仅当输出分布 Y = μ = ( 1 2 , 1 2 ) Y=\mu=(\frac12,\frac12) Y=μ=(21,21)是均匀分布。

由此,反解输入分布:设 ξ Q = μ \xi Q=\mu ξQ=μ
则当 ε ≠ 1 2 \varepsilon\neq\frac{1}{2} ε=21时, Q Q Q可逆, ξ \xi ξ有唯一解,用Cramer 法则解之可得 ξ = ( 1 2 , 1 2 ) \xi=(\frac12,\frac12) ξ=(21,21) 也是均匀分布。
ε = 1 2 \varepsilon=\frac12 ε=21 ∣ Q ∣ = 0 |Q|=0 Q=0,不可逆。此时解得: ξ 1 + ξ 2 = 1 \xi_1+\xi_2=1 ξ1+ξ2=1 ,即只要 ξ \xi ξ是一个输入分布即可,也就是任何输入分布均可达到信道容量。 □ \square

ξ Q = μ \xi Q=\mu ξQ=μ有时不一定有解,这时需采取其他方法。参看T3



T3

(对称删除信道) 设 X = { 0 , 1 } , Y = { 0 , 1 , 2 } , Q = [ 1 − ε 0 ε 0 1 − ε ε ] \mathcal{X}=\{0,1\},\mathcal{Y}=\{0,1,2\},Q=\begin{bmatrix}1-\varepsilon&0&\varepsilon\\0&1-\varepsilon&\varepsilon\end{bmatrix} X={0,1},Y={0,1,2},Q=[1ε001εεε]。试求信道容量。

解:
类似T2,可得 I ( p ; Q ) = H ( Y ) − h ( ε ) I(p;Q)=H(Y)-h(\varepsilon) I(p;Q)=H(Y)h(ε) ,但是这次不一定存在使 Y Y Y服从均匀分布的初始分布了,原因是由于一般没有(这里利用的是线性方程组有解当且仅当系数矩阵的秩等于增广矩阵的秩)

r [ 1 − ε 0 1 3 0 1 − ε 1 3 ε ε 1 3 ] = 2   , r\begin{bmatrix}1-\varepsilon&0&\frac13\\\\0&1-\varepsilon&\frac13\\\\\varepsilon&\varepsilon&\frac13\end{bmatrix}=2\:, r 1ε0ε01εε313131 =2,

于是只能利用定义进行计算了,事实上

H ( Y ) = − ( 1 − ε ) p ( 0 ) log ⁡ ( 1 − ε ) p ( 0 ) − ( 1 − ε ) p ( 1 ) log ⁡ ( 1 − ε ) p ( 1 ) − ε log ⁡ ε ≤ − ( 1 − ε ) log ⁡ 1 − ε 2 ε log ⁡ ε \begin{aligned}H(Y)&=-(1-\varepsilon)p(0)\log(1-\varepsilon)p(0)-(1-\varepsilon)p(1)\log(1-\varepsilon)p(1)-\varepsilon\log\varepsilon\\ &\leq-(1-\varepsilon)\log\frac{1-\varepsilon}2\varepsilon\log\varepsilon\end{aligned} H(Y)=(1ε)p(0)log(1ε)p(0)(1ε)p(1)log(1ε)p(1)εlogε(1ε)log21εεlogε

其中等号成立当且仅当 X X X服从均匀分布,其推导与二进熵函数类似。于是

C = − ( 1 − ε ) log ⁡ 1 − ε 2 − ε log ⁡ ε + ε log ⁡ ε + ( 1 − ε ) log ⁡ ( 1 − ε ) = 1 − ε 。口 C=-(1-\varepsilon)\log\frac{1-\varepsilon}2-\varepsilon\log\varepsilon+\varepsilon\log\varepsilon+(1-\varepsilon)\log(1-\varepsilon)=1-\varepsilon\text{。口} C=(1ε)log21εεlogε+εlogε+(1ε)log(1ε)=1ε。口



;