对称信道容量计算

二元对称信道容量计算m.tn：IBl=-(«)»*log2(s)-(l-a)**1qi2(1-a);subplnt(2,2j11plotU—lriwNanJ.ylabelCHf).fprimf「Bl=:*6.3f*fHl);lpll=p(yj/jrl>jpl2=p{yl/K2)jpL3=v(y2/xl)tpl4=p(t2/k2)p]1=127/128,pl2=1/128jplS=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)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)-<nyl).*(p22J•*Lci2(P22)-tPy2)t(p23).*1qe(pZ3)-(py2)•«Cp24).*Ut2924)subplnt(2,2j2),plot(a^HSJ^xlabelCaJ)T/labelCH2n).fprintfC监:辭.卅'胡2):sm-船；subplo+(2r2^3)plot(jjC)jilabelCaJ),ylabe1CC*),fprimfC洁道吝里:«6.rC);0 0.5 100.510 0.5 10 0.5 100.510 0.5 1二元对称信道容量随a变化时的变化曲线将程序中的a修改为0~1的某个确定值，即可用改程序计算出其相应的二元对称信道值，下例题所示，取a为0.5时，信源熵H1(即H(x))以及H(x/y)为定值，在图像上显示为定点，计算出信道容量为 0.942■'拆器•F:\MATLAB\R2014a\UntrtledSC.mUntil11asO>BHl=-(a|,#lQ£2Ca)-(：subpLm 2^1);plot■'拆器•F:\MATLAB\R2014a\UntrtledSC.mUntil11asO>BHl=-(a|,#lQ£2Ca)-(：subpLm 2^1);plot(^HI\^xlabeLpl1=127/128,pl2=l/：p21=(a).«<pll)./(dp22=(l-a).*(pl2k/p23=(ai.*1p13)./(I'ip24=(l-a).*(pl4k/pyl-(A1.*1l-p)*(p),py2=(p'丄乂ia)+tl-p)>命常行闔口 」i：'不施MATLAB?观看吃眩,Columns79throuih910.05750.0575Columns92through0.05750.0575同理，当取a为0.3时，计算出信道容量为 0.825二元非对称信道容量计算a=0:0.01:1Hl=-(a).*log2(a)-(l-a)**lo(2(l-a)；subplot(2,2jI)zplcrtHilijnjrlabelCaJ)fylabelCHl*);fpnntf『H1=:S6.3f*rHl):Kpll-ptyl/xDjpl2=p(yl/x2iplS^p(y2/xl)fpl4^=p(y2/x2)P11=1/2,PHDP13=.l/2fpl4=2/3:寫p1l=p(x1/yl)jp12=p(i2/yl》，pl3=p(xl/y2).plbp(x2/y2)|p21=(ah*(pll)+/((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(2f2r2):plot(a.HS),xlabel(JaJ)Tylabel(?H2Z):fprintf『H2:乐&3f?,H2):C=HJ-H2;subplot2*3)plot(pjpC)rxlabela*ylabelCC*);fprintf(:信道W1：K6.3f\C):

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

将程序中的a修改为0~1的某个确定值，即可用改程序计算出其相应的二元对称信道值，如下例题所示，取a为0.5时，H(x)及H(x/y)为定值，在图像中显示为一定点，同时信道容量计算出为0.102'謂釜-F:\MATLAB\I-1C□ntifledSO/na=0,5Hl=-(■).*log2(asubplot(2rZj1):plotIajH1i'謂釜-F:\MATLAB\I-1C□ntifledSO/na=0,5Hl=-(■).*log2(asubplot(2rZj1):plotIajH1iixlaIfifpviritfHl=：＞10%pll=p(xlfyl,pp2L=(a).p22=11*3),*(pl2p2J=*a''*n"'itf1-31.*(n14命令行衰口-i7^^MATLAB?through0,89770.®977C