matlab-信道容量的迭代計(jì)算

1、信息論基礎(chǔ)課程實(shí)驗(yàn)報(bào)告（四） 班級： 計(jì)算072班 姓名： 白梅 學(xué)號： 實(shí)驗(yàn)內(nèi)容：信道容量的迭代計(jì)算。實(shí)驗(yàn)要求：1. 輸入一個(gè)信道矩陣；2. 輸出下列內(nèi)容：（1）最佳輸入概率分布；（2）信道容量；（3）迭代次數(shù)；（4）輸入符號的個(gè)數(shù)；（5）輸出符號的個(gè)數(shù)；（6）迭代計(jì)算精度；3. 實(shí)驗(yàn)報(bào)告中要包括源代碼和輸入、輸出結(jié)果。實(shí)驗(yàn)源程序如下：functionCC,Paa=ChannelCap(P,k)if(length(find(P10e-10) error(Not a prob.vector,component do not add up to 1。不是概率分布);endr,s=size(P);

2、Pa=(1/(r+eps)*ones(1,r);sumrow=zeros(1,r);Pba=P;n=0;C=0;CC=1;while abs(CC-C)=k n=n+1; Pb=zeros(1,s); for j=1:s for i=1:r Pb(j)=Pb(j)+Pa(i)*Pba(i,j); end end suma=zeros(1,s); for j=1:s for i=1:r Pab(j,i)=Pa(i)*Pba(i,j)/(Pb(j)+eps); suma(j)=suma(j)+Pa(i)*Pba(i,j)*log2(Pab(j,i)+eps)/(Pa(i)+eps); end en

3、d C=sum(suma); L=zeros(1,r); sumaa=0; for i=1:r for j=1:s L(i)=L(i)+Pba(i,j)*log(Pab(j,i)+eps); end a(i)=exp(L(i); end sumaa=sum(a); for i=1:r Paa(i)=a(i)/(sumaa+eps); end CC=log2(sumaa); Pa=Paa; s0=很好！輸入正確，迭代結(jié)果如下: ; s1=最佳輸入概率分布Pa:; s2=信道容量C:; s3=迭代次數(shù)n:; s4=輸入符號數(shù)r:; s5=輸出符號數(shù)s:; s6=迭代計(jì)算精度k:; for i=1:

4、r Bi=i; endenddisp(s0);disp(s1),disp(B),disp(Paa);disp(s4),disp(r);disp(s5),disp(s);disp(s2),disp(CC);disp(s6),disp(k);disp(s3),disp(n);設(shè)置k=0.00001 k=0.00001k = 1.0000e-005輸入信道矩陣P： P=1/2 1/2 0;0.2 0.3 0.5P = 0.5000 0.5000 00.2000 0.3000 0.5000得到輸出結(jié)果： ChannelCap(P,k)很好！輸入正確，迭代結(jié)果如下: 最佳輸入概率分布Pa: 1 2 0.5949 0.4051輸入符號數(shù)r: