对称信道容量计算

m.tn： I

Bl =-(«)» *log2(s)-(l-a) * *1QI2(1-a); subpl nt (2, 2j 11

plot U—lriwN an J.ylabelCHf). fprimf「Bl=:*6. 3f*f Hl);

lpll=p(yj/jrl>jpl2=p{yl/K2)jpL3=v(y2/xl)tpl4=p(T2/K2) p] 1=127/128, pl 2=1/1 28J plS=l/128^14=127/128；

Kpl l=p (xl/yl), pl 2=p tx2/y 1) j. p 13=p (xl/y2), p 14=p i x2/y2) p2l= (n)!b«(pl 1). / ((s). * (pl 1) +(l-a) t * (pl2)); 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)-plot (jjC) jilabelC aJ), ylabe 1CC*),

}: ):

0 0.5 1 0 0.5 1 0 fprimf C 洁道吝里:«6.rC);

0.5 1 二元对称信道容量随 a变化时的变化曲线

将程序中的a修改为0~1的某个确定值，即可用改程序计算出其相应的二元对称信道值， 下例题所示，取 a为0.5时，信源熵H1 (即H(x))以及H(x/y)为定值，在图像上显示为 定点，计算出信道容量为

asO> B ■'拆器• F:\\MATLAB\\R2014a\\ Hl=-(a|, #lQ£2Ca)-(： UntrtledSC.m Until

subp Lm

2^ 1);

plot (^HI\\^xlabeL pl 1=127/128, pl2=l/： p21= (a). «

11 命常行闔口

」

i：'不施 MATLAB?观看吃眩,

Columns 79 throuih 91 0.0575

0. 0575

Columns 92 through 0.0575

0. 0575

0.942

同理，当取a为0.3时，计算出信道容量为

0.825

二元非对称信道容量计算

a=0: 0. 01:1

Hl=-(a). *log2 (a)-(l-a) * *lo(2(l-a)； subplot (2, 2j I ) z

plcrt Hilij njrlabelC aJ )f ylabelCHl*); fpnntf『H1=:S6. 3f* r Hl):

Kpll-ptyl/xDj pl2=p(yl/x2iplS^p(y2/xl)f pl4^=p (y2/x2) P11=1/2, PHD P13=.l/2f pl4=2/3:

寫p 1 l=p (x 1/yl)j p 12=p (i2/yl》，pl 3=p (x l/y2). pl bp (x2/y2)| p21= (ah *(pl l) + /((a). *(pll)+(l-a)* *(pl2)); p22= (1-a). *(pl2). /

* (pl2)+a, *(pll));

p23= (a). «(p)3)./((a). *(pl3)+(l-a)« *(pi4))； p24= < 1-a). *(pl4). / ((.L-a)・ * (p!4) + (a)a • (pl3)): pyl-(a). *(l-p) + (p)B *(Ha): py2= *(a) + (l-ph *(l~a):

H2=- (pyl), *(p21). *lc[2(p21)-(pyl)« *Cp22L *Lag2(p22)-(py2)« *(p23'. *log 'p23J- (py2). *(p24)*lo<2(p24) subplot (2f 2r 2): plot (a.HS), xlabel (J aJ ) T ylabel (? H2Z): fprintf『H2:乐& 3f?, H2): C=HJ-H2; subplot2* 3)

plot (pjp C) r xlabel a* ylabel CC*);

[

fprintf (:信道W1：K6. 3f\\C):

非对称信道容量随a变化时的变化曲线

0. ®977

将程序中的a修改为0~1的某个确定值，即可用改程序计算出其相应的二元对称信道 值，如下例题所示，取 a为0.5时，H(x)及H (x/y)为定值，在图像中显示为一定点，同时

'謂釜-F:\\MATLAB\\I -1C

□ntifledSO/n a=0, 5

Hl =-(■). *log2 (a subplot (2r ZJ 1): plot I aj H1 i i xla Ififpviritf Hl=：＞

%pll=p(xlfyl , p p2L= (a).

10

p22= 11*3), * (pl 2 p2J=*a'' * n \"'it f 1

-31 . * (n 1 4

命令行衰口

-i 7^^ MATLAB?

through

0, 8977 Columns S2 through

0. 8977

fit信1B苔重：0. L

信道容量计算出为

0.102 同理，取a为0.3时，计算出信道容量为

0.138

疑问：在代入a为一定值条件下，同 H( x)和0H(x/y)，计算出的信道容量也是一个 定值，那么在图像中也应为一定点，为何这两个程序输出结果图中显示为一直线？