1、 天津工业大学信息论基础课程实验报告学 院 理学院 年 级 2009 级 班 级 信计 092 学 号 0910820217 姓 名 赵昊 20年 2 月 20 日1实验：信道容量的计算1. 实验目的熟悉和了解代码的应用2. 实验内容通过一些代码的编写，实现信道容量的计算3. 设计代码见附录。4. 说明体会通过本次试验，熟悉了通过编写代码实现信道容量的计算，增强了编写代码的能力。并对信道容量有了一定的认识，为学习增加了基础和理解附录代码部分：function C=xindao(PY_X,PX)e=0.01;S=size(PY_X);for k=1:1:S(1)I=0;PY=zeros(1,S(
2、2);for j=1:1:S(2)for i=1:1:S(1)PY(j)=PY(j)+PX(i)*PY_X(i,j);endI=I+PY_X(k,j)*log2(PY_X(k,j)/PY(j);endB(k)=expm(I);endC1=logm(PX*B)C2=logm(max(B)while (1) if abs(C1-C2)e2break;endfor i=1:1:S(1)PX(i)=PX(i)*B(i)/(PX*B);endfor k=1:1:S(1)I=0;PY=zeros(1,S(2);for j=1:1:S(2)for i=1:1:S(1)PY(j)=PY(j)+PX(i)*PY_X(i,j);endI=I+PY_X(k,j)*log2(PY_X(k,j)/PY(j);endB(k)=expm(I);endC1=logm(PX*B);C2=logm(max(B);endC=C1 控制窗口输入：PY_X=0.5 0.25 0.125 0.125;0.25 0.5 0.125 0.125;PX=0.5 0.5;C=xindao(PY_X,PX)结果为：C1 =0.0613C2 =0.0613C =0.0613C =0.0613
