《对称信道容量计算.docx》由会员分享，可在线阅读，更多相关《对称信道容量计算.docx（6页珍藏版）》请在三一文库上搜索。
1、二元对称信道容量计算m.tn: IBl =-(?)? *log2(s)-(l-a) * *1 Qi2(1-a);subpl nt (2, 2j 11plot U lriwN an J.ylabelCHf).fprimfBl=:*6. 3f* f Hl);lpll=p(yj/jrljpl2=pyl/K2)jpL3=v(y2/xl)tpl4=p( T2/K2)p 1=127/128, pl 2=1/1 28J plS=l/128A14=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=
2、(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)./ ().* (pld) + (a).:PFl= (-)* (l-p)+(p). *(- :)P?2= (pA, *(a)+(l-p ,. *(l-a):*Ut2924) subplnt (2, 2j 2),H2A-(pyl .*(p21.*loi2(p2l)- nyl).*(p22J
3、? *Lci2 (p22)-tpy2) t(p23). *1 QE( PZ3)-(py2) ? ?Cp24).plot (aAHSJAxlabelC aJ)T/labelC H2n).fprintf C 监:就卅胡 2):sm-船；subpl o+ (2r 2A3)plot (jj C) jilabelC aJ), ylabe 1CC*),050 0.5 1fprimf C 洁道吝里:?6.rC);二元对称信道容量随a变化时的变化曲线将程序中的a修改为01的某个确定值，即可用改程序计算出其相应的二元对称信道值，下例题所示,取a为0.5时，信源嫡H1 (即H(x)以及H(x/y)为定值，在图像上
4、显示为 定点，计算出信道容量为拆器？ FaMATBABR2014a uHftedaicSiQ 先Ca)-(：Untilsubp Lm 2A 1);plot (AHIAxlabeL pl1=127/128, pl2=l/： p21(a). ?11命常行口i：不施 MATLAB?Columns 79throuih 910.05750. 0575Columns 92through0.05750. 05750.942同理，当取a为0.3时，计算出信道容量为0.825a=0: 0. 01:1Hl=-(a). *10g 2 (a)-(l-a) * *lo (2(l-a)；subplot (2, 2j I
5、) zplcrt Hilij njrlabelC aJ)f ylabelCHl* );fpnntf H1=:S6. 3f* r Hl):Kpll-ptyl/xDj pl2=p(yl/x2iplSAp(y2/xl) f pl4A=p (y2/x2)P11=1/2, PHD P13=.l/2 f pl4=2/3:W 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,
6、 *(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 *(la):H2=- (pyl), *(p21). *lc2(p21)-(pyl)? *Cp22L *Lag2(p22)-(py2)? *(p23. *log p23J- (py2). *(p24)*lorI1I %pll=p(xl fyl , p p2L=. (a).10p22= 11*3) , * (pl 2 p2J=*a * n it f 1-31 . * (n 1 4命令行衰口-i 7M MATLAB?through0. ?9770, 8977