信息论实验MATLAB求信道容量.docx

《信息论实验MATLAB求信道容量.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、formatrat%所有的数全部用分数显示sum1=0;n=input('请输入一个数字')A=ones(n)*2^-50;%定义一个极小值矩阵D=input('请输入一个矩阵')B=A+D;%用极小值矩阵加上矩阵B避免出现log2（0）的情况fori=1:nforj=1:nP(i,j)=B(i,j)*log2(B(i,j));endendF=sum(P,2);%每行求和Σr=DF;%计算公式中的β值记为rfork=1:nsum1=sum1+2^r(k,1);endC=log2(sum1)%信道容量%以下代码段是求信道容量的最佳输入分布GQ=z

2、eros(1,n);fort=1:na=2^(r(t,1)-C);Q(t)=a;endG=Q*inv(D)%因为G=*D^(-1)(D转置矩阵)