典型的关于信道容量的计算

p_yx= input('信道矩阵P_YX='); %输入信道矩阵 [r,s]=size(p_yx); b=0.01; %迭代精度 for i=1:r p(i)=1.0/r; %初始概率为均匀分布 end C0=-inf; C=1; n=0; while(abs(C-C0)/C>=b) for j=1:s q(j)=0; for i=1:r q(j)=q(j)+p(i)*p_yx(i,j); end end for i=1:r d(i)=0; for j=1:s if(p_yx(i,j)==0) d(i)=d(i)+0; else d(i)=d(i)+p_yx(i,j)*log(p_yx(i,j)/q(j)); end end a(i)=exp(d(i)); end u=0; for i=1:r u=u+p(i)*a(i); end for i=1:r p(i)=p(i)*a(i)/u; end C0=C; C=log2(u); n=n+1; end disp('信道容量'); C disp('输入分布：') p