對(duì)稱信道容量計(jì)算

1、二元對(duì)稱信道容量計(jì)算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 監(jiān):辭.卅'胡 2):sm-船；subpl o+ (2r 2A3)plot (jj C) jilabelC aJ), ylabe 1CC*),0 500.500 0.5 1fprimf C 潔道吝里:?6.rC);二元對(duì)稱信道容量隨 a變化時(shí)的變化曲線將程序中的a修改為01的某個(gè)確定值，即可用改程序計(jì)算出其相應(yīng)的二元對(duì)稱信道值

4、，下例題所示，取a為0.5時(shí)，信源熵H1 (即H(x)以及H(x/y)為定值，在圖像上顯示為定點(diǎn)，計(jì)算出信道容量為拆器？ 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同理，當(dāng)取a觀看吃眩,throuih 910. 0575through0. 0575為0.3時(shí)，0.942計(jì)算出信道容量為0.825二兀非對(duì)稱信道容量計(jì)算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:寫(xiě) 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:樂(lè)& 3f?, H2):C=HJ-H2;subplot2* 3)plot (pjp C) r xla

8、bel a* ylabel CC*);fprintf :信道 W1：K6. 3fC :非對(duì)稱信道容量隨a變化時(shí)的變化曲線將程序中的a修改為01的某個(gè)確定值，即可用改程序計(jì)算出其相應(yīng)的二元對(duì)稱信道值，如下例題所示，取a為0.5時(shí)，H(x)及H (x/y )為定值，在圖像中顯示為一定點(diǎn)，同時(shí)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