对称信道容量计算

1、二元对称信道容量计算m.t n: IBl =-(?)? *log2(s)-(l-a) * *1 Qi2(1-a);subpl nt (2, 2j 11plot U IriwN an J.ylabelCHf).fprimfBl=:*6. 3f* f Hl);lpll=p(yj/jrl>jpl2=pyl/K2)jpL3=v(y2/xl)tpl4=p( T2/K2)p仁 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)p

2、2l= (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)./ -a). * (pld) + (a).:PFl= ( ) (l-p)+(p). *(- ：)P?2= (pA, *(a)+(l-p ' *(l-a):H2A-(pyl .*(p21>.*loi2(p2l)- <n yl).*(p22J

3、? *Lci2 (p22)-tpy2) t(p23). *1 QE( PZ3)-(py2) ? ?Cp24). *Ut2924) subpl nt (2, 2j 2), plot (aAHSJAxlabelC aJ)T/labelC H2n).fprintf C 监:辭.卅'胡 2):sm-船；subpl o+ (2r 2A3)plot (jj C) jilabelC aJ), ylabe 1CC*),0 500.500 0.5 1fprimf C 洁道吝里:?6.rC);二元对称信道容量随 a变化时的变化曲线将程序中的a修改为01的某个确定值，即可用改程序计算出其相应的二元对称信道值

4、，下例题所示，取a为0.5时，信源熵H1 (即H(x)以及H(x/y)为定值，在图像上显示为定点，计算出信道容量为拆器？ FaSOATBABR2021aUkir=e(aS(C.#Q Ca)-(：U ntilsubp Lm2A 1);plot (AHIWxlabeL pl仁 127/128, pl2=l/ : p21(a). ?<pll)./(d p22=(l-a).*(pl2k/ p23= .* 1 p 13). / (I' i p24=(l-a).*(pl4k/ pyl-(A1.* 1l-p)*(p),py2=(p'丄 乂 ia)+tl-p)11命常行闔口i：'

5、不施 MATLAB?Columns 790.0575Colum ns 920.0575同理，当取a观看吃眩,throuih 910. 0575through0. 0575为0.3时，0.942计算出信道容量为0.825二兀非对称信道容量计算a=0: 0. 01:1Hl=-(a). *log 2 (a)-(l-a) * *lo (2(l-a)；subplot (2, 2j I ) zplcrt Hilij njrlabelC aJ)f ylabelCHl* );fpnntf H1=:S6. 3f* r Hl):Kpll-ptyl/xDj pl2=p(yl/x2ip£Ap(y2/xl)f

6、 pl4A=p (y2/x2)P1仁 1/2, PHD P13=.l/2 f 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

7、):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)*lo<2(p24) subplot (2 f 2r 2): plot (a.HS), xlabel (J aJ) t ylabel (?H2Z):fprintfH2:乐& 3f?, H2):C=HJ-H2;subplot2* 3)plot (pjp C) r xla

8、bel a* ylabel CC*);fprintf :信道 W1：K6. 3fC :非对称信道容量随a变化时的变化曲线将程序中的a修改为01的某个确定值，即可用改程序计算出其相应的二元对称信道值，如下例题所示，取a为0.5时，H(x)及H (x/y )为定值，在图像中显示为一定点，同时0. ?977'謂釜-F:MATLABI-1C10 n tifledSO/ na=0, 5Hl =-( ). *log2 (a subplot (2rZJ 1)： plot aj H1 i i xla IfifpviritfHl=:%pll=p(xl fyl , p p2L=(a).p22= 11*3) , * (pl 2p2J=*a'' * n "'it f 1-31 . * (n 1 4命令行衰口-i 7AA MATLAB?0, 89770. 8977fit信1B苔重：0. LthroughColu mn