对称信道容量计算.pdf
二元对称信道容量计算
m.tn:I
Bl=-(«)»*log2(s)-(l-a)**1QI2(1-a);
subplnt(2,2j11
plotU—lriwNanJ.ylabelCHf).
fprimf「Bl=:*6.3f*Hl);
f
lpll=p(yj/jrljpl2=p{yl/K2)jpL3=v(y2/xl)pl4=p(2/2)
TK
t
p]1=127/128,pl2=1/128plS=l/128^14=127/128;
Kpll=p(xl/yl),pl2=ptx2/y1)j.p13=p(xl/y2),p14=pix2/y2)
p2l=(n!b«(pl1)./((s).*(pl1)+(l-a)*(pl2));
)t
p22=(]-a),*(pl2)./1]-a).*pl2)+a.):
p23=(a)-*(pl3)./((a).*(pl3)+(Ha).*(pl4));
p24-(1-a).■(pit)./(1-a).*(pld)+(a).}:
PFl=(■)-*(l-p)+(p).*(]-■):
|
P?2=(p^,*(a)+(l-p.*(l-a:
)
H2^-(pyl}.*(p21.*loi2(p2l)-nyl).*(p22J•*Lci2(22)-ty2)t(p23).*1(Z3)-(py2)•«Cp24).*U2924)subplnt(2,2j2),
PPQEPt
n
plot(a^HSJ^xlabelCa)/labelCH2).
T
fprintfC监:辭.卅胡2):
sm-船;
subplo+(22^3)
r
plot(jjC)jilabelCa),ylabe1CC*),
00.5100.51
00.51
fprimfC洁道吝里:«6.C);
r
二元对称信道容量随a变化时的变化曲线
将程序中的a修改为0~1的某个确定值,即可用改程序计算出其相应的二元对称信道值,下例题所示,
取a为0.5时,信源熵H1(即H(x))以及H(x/y)为定值,在图像上显示为定点,计算出信道容量为