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(ε)=1−h(ε).
其中 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−ε。口