信息论笔记

还是复习。

第一章 绪论

不提

第二章 离散信息的度量

习题 2.6
在某地篮球联赛的每个赛季中,最终只有 A、B 两队进入决赛争夺冠军。决赛采用 7 场 4 胜制,首先赢得 4 场胜利的球队获得冠军,并结束比赛。把产生冠军的事件 x 用 A、B 两队各场次比赛结果表示,作为信源 X 产生的随机事件,例如:AAAA 表示事件“A 队前 4 场获得冠军”;ABBAAA 表示事件“ A 队在第 1、4、5、6 场取胜获得冠军(而 B 队在第 2、3 场取胜)”……假设两队在每场比赛中的取胜机会均等,每场比赛只有“A 胜”或“B 胜”两种结果,并且各场比赛的结果是彼此独立的。

(1)求信源的熵 H(X)。

事件 数目 单事件的概率 概率的和
赛 4 场获冠军(AAAA) 1 1/16 1/16
赛 5 场获冠军(前 4 场 B 胜 1 场) 4 1/32 1/8
赛 6 场获冠军(前 5 场 B 胜 2 场) 10 1/64 5/32
赛 7 场获冠军(前 6 场 B 胜 3 场) 20 1/128 5/32
总概率 1/2

H(X)=2×[(1/16)log216+4×(1/32)log232+10×(1/64)log264+20×(1/128)log2128]=186/32=5.8125 比特/符号

注意到上面乘 2 是因为 A 队和 B 队的概率是相等的。

(2)求事件“两队打满 7 场”所提供的信息量。
“两队打满7场”事件数为 40,所求概率为 40×(1/128)=5/16 ,事件“两队打满7场”所提供的信息量为
I1=log2(5/16)=1.6781 bit
(3)列出 A 队前 3 场都失利的所有情形,求“A 队前 3 场都失利”所提供的信息量。
如表所示:

事件 概率 结果
BBBB 1/16 B 夺冠
BBBAB 1/32 B 夺冠
BBBAAB 1/64 B 夺冠
BBBAAAB 1/128 B 夺冠
BBBAAAA 1/128 A 夺冠

概率和就是 1/8,所以有

I2=log28=3 bit
(4)求事件“A 队在前 3 场获胜的条件下又取得冠军”所提供的信息量。
A队在前 3 场失利的条件下又取得冠军的条件概率为
1/1281/16+1/32+1/64+2×(1/128)=116
“A队在前 3 场失利的条件下又取得冠军”的信息量为
I3=log216=4 bit

习题 2.7
已知随机变量 X,Y 的联合概率分布为 PXY(ak,bj) ,满足:
PX(a1)=1/2, PX(a2)=PX(a3)=1/4,
PY(b1)=2/3, PY(b2)=PY(b3)=1/6。试求能使 H(XY) 取最大值的 XY 的联合概率分布。

因为 H(XY)H(X)+H(Y)当且仅当 X,Y 独立时,等号成立。所以当 X,Y 独立时,H(XY) 取最大值。因此 XY 的联合概率分布为

p(a1b1)=13,p(a1b2)=112,p(a1b3)=112,p(a2b1)=16,p(a2b2)=124,p(a2b3)=124,p(a3b1)=16,p(a3b2)=124,p(a3b3)=124.

习题 2.16
两随机变量 X,Y ,联合概率 p(x,y) 如下:

y=0 y=1
x=0 1/8 3/8
x=1 3/8 1/8

Z=XY(一般乘积),试计算:

(1) H(X),H(Y),H(Z),H(XZ),H(YZ),H(XYZ);
(2) H(X|Y),H(Y|X),H(X|Z),H(Z|X),H(Y|Z),H(Z|Y),H(Z|YZ),H(Y|XZ),H(Z|XY);
(3) I(X;Y),I(X;Z),I(Y;Z),I(X;Y|Z),I(Y;Z|X),I(X;Z|Y)
p(xz)=yp(xyz);p(yz)=xp(xyz);p(xyz)=p(xy)p(z|xy)
其概率取值容易计算如下:

PX(0)=1/2,即 P(X=0)=1/2,下同
PX(1)=1/2

PY(0)=1/2
PY(1)=1/2

由表知,x=0y=0 都能让 Z 为 0,故
PZ(0)=1/8+3/8+3/8=7/8
PZ(1)=1/8

(1)
主要利用自信息的平均值为熵,即
H(X)=Ep(x)[logp(x)]
注意 log 默认是以 2 为底的
比如 H(XZ),先列出所有可能,知 P(X=0,Z=0)=1/2P(X=1,Z=0)=3/8P(X=1,Z=1)=1/8,所以 H(XZ)=(1/2)log(1/2)(3/8)log(3/8)(1/8)log(1/8)=1.406 比特/符号
H(X)=H(1/2)=1 比特/符号H(Y)=H(1/2)=1 比特/符号H(Z)=H(1/8)=0.544 比特/符号H(XZ)=H(1/2,3/8,1/8)=1.406 比特/符号H(YZ)=H(XZ)=1.406 比特/符号H(XYZ)=H(1/8,3/8,3/8,1/8)=1.811 比特/符号

(2)
主要利用熵的可加性,即 H(XY)=H(X)+H(Y|X)
H(X|Y)=H(XY)H(Y)=1.8111=0.811 比特/符号H(Y|X)=H(XY)H(X)=1.8111=0.811 比特/符号H(X|Z)=H(XZ)H(Z)=1.4060.544=0.862 比特/符号H(Z|X)=H(XZ)H(X)=1.4061=0.406 比特/符号H(Y|Z)=H(YZ)H(Z)=1.4060.544=0.862 比特/符号H(Z|Y)=H(YZ)H(Y)=1.4061=0.406 比特/符号H(X|YZ)=H(XYZ)H(YZ)=1.8111.406=0.405 比特/符号H(Y|XZ)=H(XYZ)H(XZ)=1.8111.406=0.405 比特/符号H(Z|XY)=H(XYZ)H(XY)=1.811H(1/8,3/8,3/8,1/8)=0

(3)
主要利用平均互信息与熵的关系,即
I(X;Y)=H(X)H(X|Y)=H(Y)H(Y|X)=H(X)+H(Y)H(XY)
I(X;Y)=H(X)H(X|Y)=10.811=0.189 比特/符号I(X;Z)=H(X)H(X|Z)=10.862=0.138 比特/符号I(Y;Z)=H(Y)H(Y|Z)=10.862=0.138 比特/符号I(X;Y|Z)=H(X|Z)H(X|YZ)=0.8620.405=0.457 比特/符号I(Y;Z|X)=H(Y|X)H(Y|XZ)=0.8110.405=0.406 比特/符号I(X;Z|Y)=H(X|Y)H(X|YZ)=0.8110.405=0.406 比特/符号

第三章 离散信源

重点公式:πT=πTP

习题 3.9
x1,x2,,xn1,xn, 为平稳序列(未必是马氏链),那么下面的论断哪些是正确的?对正确的进行证明,对错误的举出反例。(提示:下面论断至少有一个是错的)
(1) H(Xn|X0)=H(Xn|X0);
正确,由于平稳性有
H(Xn|X0)=H(X0Xn)=H(X0Xn)H(Xn)=H(XnX0)H(X0)=H(Xn|X0)
(2) H(Xn|X0)H(Xn1|X0);
错误,若 x0,,xn1,xn 构成马氏链,则
I(X0;Xn1)I(X0;Xn)H(Xn1)H(Xn1|X0)H(Xn)H(Xn|X0)H(Xn1|X0)H(Xn|X0) (利用平稳性)H(Xn|X0)H(Xn1|X0)
但如果序列具有周期性,且周期为 n,其中,x0,x1,,xn 是独立同分布的二元等概率序列,则
xk=xkn (kn)x0=xnH(Xn|X0)=0, H(Xn1|X0)=H(Xn1)=log2H(Xn|X0)<H(Xn1|X0)
(3) H(Xn|X1X2Xn1)n 的增函数;
错误
定理 3.3 对任意离散平稳信源,若 H1(X)< ,有:

H(XN|X1XN1) 不随 N 而增加;
HN(X)H(XN|X1XN1)
HN(X) 不随 N 而增加;
H(X) 存在,且 H(X)=limNH(XN|X1XN1)

(4) H(Xn|X1,,Xn1,Xn+1,,X2n)n 的非增函数。
正确
H(Xn+1|X1,,Xn,Xn+2,,X2n+2)=H(Xn|X0,,Xn1,Xn+1,,X2n+1)H(Xn|X1,,Xn1,Xn+1,,X2n)

习题 3.10
一个 2 状态马氏链的转移概率矩阵为
P=(3/41/41/43/4)
并假定初始状态概率矢量为 p(0)=(1 0);求
(1) Pnp(n), n=1,2,3
P1=P=(3/41/41/43/4),P2=P×P=(3/41/41/43/4)(3/41/41/43/4)=(5/83/83/85/8)
P3=P2×P=(5/83/83/85/8)(3/41/41/43/4)=(9/167/167/169/16)
所以,有
p(1)=p(0)P=(3/41/4);p(2)=p(0)P2=(5/83/8);p(3)=p(0)P3=(9/167/16)

(2) Pnp(n) 的一般形式。
由上面已经求得的结果可以类推得
Pn=(2n+12n+12n12n+12n12n+12n+12n+1),p(n)=p(0)Pn=(10)(2n+12n+12n12n+12n12n+12n+12n+1)=(2n+12n+12n12n+1)

习题 3.15
3.15 黑白气象传真图的消息只有黑色和白色两种,即信源 X=, ;设黑色出现的概率为 p()=0.3 ,白色的出现概率 p()=0.7

(1) 假设图上黑白消息出现前后没有关联,求熵 H(X)
假设黑白消息出现的前后没有关联,则等效于一个离散无记忆信源,概率空间为
(XP(x))=(0.30.7),i=12P(xi)=1
信源的熵为
H(X)=i=12P(xi)logP(xi)=0.7log20.70.3log20.3=0.881比特/符号
(2) 假设消息前后有关联,其依赖关系为 P(|)=0.9, P(|)=0.1, P(|)=0.2, P(|)=0.8,求此一阶马氏链的熵率 H2
假设黑白气象传真图的消息前后有关联,其状态集 E=A=, ,可以得到其状态转移矩阵为
P=(0.90.10.20.8)

此马尔可夫链状态转移矩阵有,则状态平稳分布存在。设状态的平稳分布为 πT=(π1,π2) ,有
πT=πTP,π1+π2=1,π1=2/3,π2=1/3.

此一阶马氏源的熵为
H2=23H(0.1,0.9)+13H(0.2,0.8)0.553比特/符号
(3) 分别求上述两种信源的剩余度,并比较 H(X)H2 的大小,并说明其物理意义。
前后没有关联情况下信源的剩余度:
γ1=10.881H(1/2)=10.8811=0.119

一阶马尔可夫信源的剩余度:
γ2=10.553H(1/2)=10.5531=0.447

H(X)>H2,γ1<γ2

物理意义:一阶马尔可夫信源考虑了消息前后的关联,使得符号熵减少,信源的冗余增加,同时编码的熵增加。

习题 3.21
给了一个一阶马氏链的状态转移图如图,符号集为 0,1,2

(1) 求状态平稳分布 (π0,π1,π2) 和马氏链熵率。
由状态图可得状态转移矩阵
P=(1pp/2p/2p/21pp/2p/2p/21p)
πTP=πT,即 π=(π0,π1,π2) 满足
(π0,π1,π2)(1pp/2p/2p/21pp/2p/2p/21p)=(π0,π1,π2)
以及 π0+π1+π2=1
解得
(π0,π1,π2)=(1/3,1/3,1/3)
信源的符号熵为
H=H2=π0H(1p,p2,p2)+π1H(p2,1p,p2)+π2H(p2,p2,1p)=(1p)log2(1p)plog2p+p 比特/符号

(2) 当 p 为何值时,信源熵率达到最大值?当 p=01 时,结果如何?
因为 H=(1p)log(1p)plogp+p,对 p 求一阶导数:
H(p)=log(1p)+1ln2logp1ln2+1=log2(1p)p

H(p)=0,得 log2(1p)p=0,所以 2(1p)p=1

所以 p=2/3 时,H 达到最大值:H 的最大值等于 log231.585 比特/符号;
p=0 时,H=0;当 p=1 时,H=(1p)log2(1p)plog2p+p=1 比特/符号。

(3) 如果将信源看成无记忆的且以平稳分布为概率分布,求信源的熵率。
(XP(ai))=(0121/31/31/3)


H(X)=i=13P(ai)logP(ai)=log231.585比特/符号

由此计算结果可知 H(X)=H

第四章 连续信息与连续信源

习题 4.1
(1) 指数概率密度 p(x)=λeλx,x0;
h(X)=0p(x)logp(x)dx=0λeλxlog(λeλx)dx=0λeλxlogλdx0λeλxlogeλxdx=logeλ
(2) 拉普拉斯概率密度 p(x)=(1/2)λeλ|x|
h(X)=p(x)logp(x)dx=12λeλ|x|log(12λeλ|x|)dx=log12λ12λeλ|x|dx+12λeλ|x|log(eλ|x|)dx=log2eλ

习题 4.13
给定两连续随机变量 XY ,其中 X 的概率密度是 p(x)=ex(0x<) ,条件概率密度是 p(y/x)=xexy(0y<) 。求 h(X),h(Y),h(XY),I(X;Y)
h(X)=0p(x)logp(x)dx=0exlogexdx=(loge)0xdex=logep(y)=0p(x)p(y|x)dx=0xex(1+y)dx=0x1+ydex(1+y)dx=1(1+y)2h(Y)=01(1+y)2log1(1+y)2dy=2011+ylog(1+y)d11+y=2logeh(XY)=00xex(1+y)log[xex(1+y)]dxdy=0xlogx(0ex(1+y)dy)dx+(loge)0x2(0(1+y)ex(1+y)dy)dx=0xlogx(exx)dx+(loge)0x(ex+0ex(1+y)dy)dx=0exlogxdx+(loge)x(ex+exx)dx=γloge+2loge=(γ+2)logeI(X;Y)=h(X)+h(Y)h(XY)=(1γ)loge

其中,γ=0exlnxdx 为欧拉常数,定义为调和级数与自然对数的差值,约等于 0.577 2。

习题 4.14
给定两连续随机变量 XY ,它们的联合概率密度是
p(xy)=12πσxσyexp{(xmx)22σx2(ymy)22σy2},<x,y<

(1) 求随机变量 U=X+YV=XY 的概率密度函数 p(u)p(v)
p(x)=p(xy)dy=12πσxexp[(xmx)22σx2]p(y)=p(xy)dx=12πσyexp[(ymy)22σy2]p(x)p(y)=p(xy)XYUmx+myσx2+σy2p(u)=12π(σx2+σy2)exp[(zmxmy)22(σx2+σy2)]Vmxmyσx2+σy2p(v)=12π(σx2+σy2)exp[(zmx+my)22(σx2+σy2)]

(2) 计算 h(U),h(V)I(U;V)

h(U)=h(V)=12log[2πe(σx2+σy2)]h(XY)=log(2πeσxσy)(uv)=(1111)(xy)|det(1111)|=2h(UV)=h(XY)+log2=log(4πeσxσy)I(U;V)=h(U)+h(V)h(UV)=log(2πe(σx2+σy2))log(4πeσxσy)=logσx2+σy22σxσy

习题 4.20
XN,YN 为定义在 RN 空间中的两个 N 维矢量,UN,VN 分别为 XN,YN 的可逆线性变换,即 u=Ax+α, v=By+β ,证明
I(UN;VN)=I(XN;YN)

(UNVN)=(A00B)(xy)+(αβ)

C=(A00B),z=(xy),γ=(αβ),(UNVN)=Cz+γ

I(UN;VN)=h(UN)+h(VN)h(UNVN)=h(Ax+α)+h(By+β)h(Cz+γ)=h(XN)+log|det(A)|+h(YN)+log|det(B)|h(XNYN)log|det(C)|=h(XN)+h(YN)h(XNYN)=I(XN;YN)

上面用到了行列式的性质:
det(C)=|A00B|=det(A)det(B)

第五章 无失真信源编码

习题 5.1
有一信源,它有 6 个可能的输出,其概率分布见下表,表中给出了对应的码 A,B,C,D,E 和 F。问:

(1) 这些码中哪些是即时码;
(2) 哪些是唯一可译码,并对所有唯一可译码,求出其平均码长 L

消息 P(ai) A B C D E F
a1 1/2 000 0 0 0 0 0
a2 1/4 001 01 10 10 10 100
a3 1/16 010 011 110 110 1100 101
a4 1/16 011 0111 1110 1110 1101 110
a5 1/16 100 01111 11110 1011 1110 111
a6 1/16 101 011111 111110 1101 1111 011

(1) 即时码就是没有一个码字是另一个码字的前缀。比如 B 里 0 是 01 的前缀,D 里 110 是 1101 的前缀,F 里 0 是 011 的前缀,所以只有 A, C, E 是即时码。

(2) 观察表中这些码组,A 是等长码,其中没有相同的码字,所以 A 是唯一可译码。其他码组都是变长码,可采用唯一可译变长码来判断:码组 B、C、E 是唯一可译码,码组 D、F 不是唯一可译码。唯一可译码的平均码长为
L=i=1qP(si)li

因此,LA=3 码符号/信源符号,LB=2.125 码符号/信源符号,LC=2.125 码符号/信源符号,LD=2 码符号/信源符号。

习题 5.5
是否存在码长分别为 1, 2, 2, 2, 2, 2, 3, 3, 3, 3 的唯一可译三元变长码?是否可以构造一个码长为 1, 2, 2, 2, 2, 2, 3, 3, 3 的即时码?存在多少这样的码?
(1) 因为 i2li=21+5×22+4×23=94>1,所以不存在码长分别为 1, 2, 2, 2, 2, 2, 3, 3, 3, 3 的唯一可译三元变长码。

(2) 可以构造码长为 1, 2, 2, 2, 2, 2, 3, 3, 3 的即时码,例如:0, 10, 11, 12, 20, 21, 220, 221, 222。

(3) 因为一个即时码码字集合与一棵码树的叶有一一对应的关系,不同的码树对应不同的码字集合。现分析这种码树生成过程:一个根节点延伸出 3 个 1 阶节点,其中有 2 个 1 阶节点要各自延伸出 3 个 2 阶节点,方法数是 3;6 个 2 阶节点中有 1 个延伸出 3 个 3 阶节点,方法数是 6,所以不同码树数目为 3 × 6 = 18,其中的一棵码树如图 5.5 所示。一个码字集合可以根据不同的对应关系分配给信源符号,如果信源符号和码字都是 9 个,则有 9! 种分配方式。所以,存在的编码方式数为 18 × 9!。

习题 5.7
设信源 S ,符号集为 s1,s2 ,其中,P(s1)=0.1

(1) 求信源的熵和信源剩余度。
H(S)=i=12p(si)logp(si)0.469 比特/符号
剩余度 γ=1H(S)log22=0.531=53.1
(2) 设码符号为 A=0,1 ,编出 S 的最优码,并求其平均码长。
码符号 X=0,1,对信源 S 编紧致码为:s00,s11。其平均码长为 1 码符号/信源符号。
(3) 把信源的 N 次扩展源 SN 编成最优码,求 N=2,3,4, 时的平均码长 LN/N
N=2 时,
S2=a1=s1s1,a2=s1s2,a3=s2s1,a4=s2s2
P(ai)=[0.01,0.09,0.09,0.81]

紧致码(即哈夫曼码)为:

码字 0 10 110 111
码长 1 2 3 3

平均码长 LNN=1Ni=14P(ai)li=0.645 码符号/信源符号

同理,当 N=3 时,平均码长 L33=0.533 码符号/信源符号

N=4 时,平均码长 L44=0.493 码符号/信源符号

N= 时,紧致码的平均码长为:
limNLNN=H(S)logr0.469 码符号/信源符号
(4) 计算当 N=1,2,3,4 时的编码效率和码剩余度。
编码效率 η=Hr(S)L=H(S)Lr=2),码剩余度 1η=1Hr(S)L=1H(S)L

所以 N=1 时,编码效率为 0.469,码剩余度为 0.531
N=2 时,编码效率为 0.727,码剩余度为 0.273
N=3 时,编码效率为 0.880,码剩余度为 0.12
N=4 时,编码效率为 0.951,码剩余度为 0.049

习题 5.8
某离散无记忆信源有8个信源符号 a0,a1,a2,a3,a4,a5,a6,a7 ,各符号的概率分别为 0.1, 0.1, 0.1, 0.1, 0.1, 0.4, 0.05, 0.05。

(1) 对信源进行码长方差最小的二元 Huffman 编码,求平均码长、码长的方差,以及码率和编码效率。

信源符号 码字
a0 010
a1 100
a2 101
a3 110
a4 111
a5 00
a6 0110
a7 0111

l=5×0.1×3+0.4×2+0.1×4=2.7 码元/信源符号

σ2=ipili2l2=5×0.1×32+0.4×22+2×0.05×422.72=0.41
H(X)=2.622 比特/信源符号,R=llog22=2.7 比特/信源符号

编码效率为:
η=H(X)R=2.6222.797.11

(2) 将信源符号编成香农码,求平均码长、码长的方差,以及码率和编码效率。
无参考答案,嘻嘻。

第六章 离散信道及其容量

习题 6.6
设二元对称信道的概率转移矩阵为
(3/41/41/43/4)

(1) 若 p(0)=1/3 ,求 I(x=0;y=1),I(x=1;Y),I(X;Y)
设输出概率为 q0,q1,有
(q0,q1)=(p0,p1)P=(1/32/3)(3/41/41/43/4)=(5/127/12)
I(x=0;y=1)=logp(y=1|x=0)p(y=1)=log21/47/12=1.22 bitI(x=1;Y)=jpijlogpijqj=14log21/45/12+34log23/47/12=0.088 比特/符号I(x=0;Y)=jpijlogpijqj=34log23/45/12+14log21/47/12=0.33 比特/符号I(X;Y)=ip(x)I(x;Y)=13I(x=0;Y)+23I(x=1;Y)=0.169 比特/符号

(2) 求该信道的容量及其达到容量时的输入概率分布。
该信道为二元对称信道,达到容量时输入等概,即 p0=p1=1/2,输出也等概。
C=1H(14)=10.811=0.189 比特/符号

习题 6.11
一个离散无记忆信道如图6.22所示:

(1) 写出该信道的转移概率矩阵;
信道的转移概率矩阵为
P=(1/21/41/40101/41/41/2)
(2) 该信道是否为对称信道?
该信道不是对称信道。
(3) 求该信道的信道容量;
Pβ=h,可得
(1/21/41/40101/41/41/2)(β0β1β2)=(12log12+14log14+14log14012log12+14log14+14log14)

解得 β0=2,β1=0,β2=2

C=log(22+20+22)=log232=0.585 比特/符号
(4) 求达到信道容量时的输出概率分布。
q0=q2=16,q1=23

p0=p2=29,p1=59

习题 6.12
一个 Z 信道的转移概率如图6.23所示:
(1) 求信道容量;
?考试考这个我似了算了
信道转移矩阵为
P=(10ϵ1ϵ),P1=(10ϵ1ϵ11ϵ)
h=(0H(ϵ)),β=P1h=(0H(ϵ)1ϵ),q0=11+2H(ϵ)1ϵ,q1=2H(ϵ)1ϵ1+2H(ϵ)1ϵ

(p0p1)=(q0q1)P1=(1ϵϵ1ϵ1+(1ϵ)ϵϵ1ϵϵϵ1ϵ1+(1ϵ)ϵϵ1ϵ)

输入概率均大于零,所以信道容量为
C=log(1+2H(ϵ)1ϵ)
(2) 若将两个同样的 Z 信道串接,求串接后信道的转移概率矩阵;
(3) 求 (2) 中串接信道的容量和达到容量时的输入的概率分布;
(4) 将 n 个同样的 Z 信道串接,求串接后信道的转移概率矩阵和信道容量。

习题 6.14
给定如图6.25所示的级联信道,求:
感觉这题答案有点问题
(1) XY 之间的信道容量 C1
由于 XY 之间的信道转移概率矩阵为
P1=(1ϵϵϵ1ϵ)

所以信道为弱对称信道,当输入等概时达到容量。

C1=log2H(ϵ)=log2+ϵlogϵ+(1ϵ)log(1ϵ)
(2) YZ 之间的信道容量 C2
YZ 之间的信道转移矩阵为
P2=(3/401/403/41/4)

由于信道为弱对称信道,当 p1=p2=1/2 时达到信道容量 C2,此时输出概率分布为 (3/8,3/8,1/4),所以信道容量
C2=H(3/8,3/8,1/4)H(1/4,3/4)=3/4 比特/符号
(3) XZ 之间的信道容量 C3 及达到容量时的输入概率分布。
XZ 之间的转移概率矩阵 P3=P1×P2,所以
P3=(1ϵϵϵ1ϵ)(3/41/4001/43/4)=(3(1ϵ)/41/4(3ϵ)/4(3ϵ)/41/43(1ϵ)/4)

该信道是弱对称信道,当输入等概率分布时达到信道容量 C3,此时
q0=q1=3/8,q2=1/4

所以
C3=H(38,38,14)H(3(1ϵ)4,14,3ϵ4)=34[log2H(ϵ)]

输入等概时达到容量。

第七章 有噪信道编码

习题 7.3
一信道输入符号集 A=0,1/2,1 ,输出符号集 B=0,1 ,信道的转移概率矩阵为
P=(101/21/201)

现有 4 个等概率消息通过此信道输出,若选择这样的信道编码:Cb:a1,a2,1/2,1/2, ai=01 (i=1,2) ,码长为 4,并选择如下译码规则:
f:(y1,y2,y3,y4)=(y1,y2,1/2,1/2)

(1) 编码后信息传输速率等于多少?
编码后信息传输速率:
R=log244=0.5 比特/符号
(2) 证明在此译码规则下,对于码字的译码错误率 PE(i)=0,i=1,2,3,4
设 4 个消息的编码分别为:(0, 0, 1/2, 1/2),(0, 1, 1/2, 1/2),(1, 0, 1/2, 1/2),(1, 1, 1/2, 1/2),通过信道传输后,每个码字的前两位无差错,所以不同的码字的传输得到不同的译码结果,无译码错误。

习题 7.4
一个二元对称信道的转移概率矩阵为
P=(1ppp1p) (p<1/2)

信道输入符号 0,1 的概率分别为 ω,1ω

(1) 求利用 MAP 准则的判决函数和平均错误率。
MAP 准则

联合概率矩阵:
((1p)ωpωp(1ω)(1p)(1ω))

当接收到 “0” 时:
(1p)ωp(1ω),则判断为 “0”;反之,则判断为 “1”。

所以
G(y=0)={0(ωp)1(ω<p)

同理,得
G(y=1)={0(ω1p)1(ω<1p)

因为 p<12,所以 1p>p。可总结如下:

  • MAP 判决函数:
    G(y)={0(ω1p)1(ω<p)
    G(y=0)=0,G(y=1)=1(pω<1p)

  • 平均错误率:
    pE={1ω(ω1p)ω(ω<p)p(pω<1p)


    pE=min(1ω,ω,p)

(2) 求利用 ML 准则的判决函数和平均错误率。
利用 ML 准则

当接收到 “0” 时,若 (1p)p,则判断为“0”;反之,则判断为 “1”,而根据题意有 1p>p

所以
G(y=0)=0

同理,得
G(y=1)=1

所以,有

  • ML 判决函数:
    G(y=0)=0,G(y=1)=1

  • 平均错误率:
    pE=p

(3) 什么时候上述两准则的判决结果相同?
p<ω1p 时,上述两准则的判决结果相同。

第八章 波形信道

这个例题是关于离散时间无记忆加性噪声信道的输入和输出的互信息、信道容量及其概率分布的求解。具体问题和解答如下:

习题 8.3

一个离散时间无记忆加性噪声信道的输入 X 限制在 [2,2],独立于 X 的噪声 Z[1,1] 区间均匀分布,熵为 h(Z)。信道输出 Y 的熵为 h(Y)

(1)写出信道输入 X 与输出 Y 的平均互信息 I(X;Y) 的表达式。
(2)求信道容量和达到容量时的输出概率分布。
(3)求达到容量时的输入概率分布。

(1)I(X;Y) 的表达式: I(X;Y)=h(Y)h(Z)

(2)因为 2X21Z1y=x+z,所以 y 的范围是:3<y<3,噪声熵 h(Z)=log(1+1)=log2,所以当 Y 有最大熵时,信道达到容量,此时 Y 应在 (3,3) 范围均匀分布,Y 的分布密度为
pY(y)={1/63<y<30其他

信道容量 maxI(X;Y)=log2(3+3)log22=log23=1.585 比特/自由度。

(3)因为 y=x+z,且 xz 相互独立,则 y 的概率密度可以由 xz 的概率密度卷积得到,设 pX(x)X(ω), pZ(z)Z(ω), pY(y)Y(ω),其中 表示傅里叶变换关系,有 X(ω)Z(ω)=Y(ω)

Y(ω)=33(16)ejωydy=sin(3ω)3ω

Z(ω)=11(12)ejωzdz=sinωω

X(ω)=Y(ω)Z(ω)=13(2cos2ω+1)

做反变换,得

X(ω)=Y(ω)Z(ω)=13(2cos2ω+1)13[δ(x2)+δ(x+2)+δ(x)]

所以达到容量时,X 的概率分布是

p(x=0)=p(x=2)=p(x=2)=13

习题 8.6

设离散时间连续信道的输入与输出分别为 XN=(X1,,XN)YN=(Y1,,YN),试证明:

(1)信源无记忆时,有 I(XN;YN)i=1NI(Xi;Yi),当且仅当信道无记忆时等式成立。
i=1NI(Xi;Yi)I(XN;YN)=i=1N[H(Xi)H(Xi|Yi)]i=1NH(Xi)+H(XN|YN)=H(XN|YN)i=1NH(Xi|Yi)=i=1NH(Xi|YNX1X2Xi1)i=1NH(Xi|Yi)0
故:
I(XN;YN)i=1NI(Xi;Yi)

(2)信道无记忆时,有 I(XN;YN)i=1NI(Xi;Yi),当且仅当信源无记忆时等式成立。
I(XN;YN)i=1NI(Xi;Yi)=H(YN)H(YN|XN)i=1N[H(Yi)H(Yi|Xi)]=H(YN)i=1NH(Yi)0
故:
I(XN;YN)i=1NI(Xi;Yi)

综上所述,当信源和信道无记忆时,这两个等式分别成立。

第九章 信息率失真函数

不用考,不学啦!

作者

未央

发布于

2024-06-21

更新于

2025-04-23

许可协议

评论

来发评论吧~
Powered By Valine
v1.4.16