一绘制二元熵函数曲线报告.doc

《一绘制二元熵函数曲线报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验一绘制二元熵函数曲线实验报告一、实验目的1.熟悉MATLAB工作环境及工具箱2.理解熵函数表达式及其性质二、实验内容用MATLAB软件编程绘制二元熵函数曲线三、实验过程1.复习二元熵函数，理解二元信源的熵H(w)=-wlogw-(1-w)log(1-w)表达式。2.熟悉MATLAB软件。1)MATLAB的操作界面MATLAB操作界面主要分为：任务栏、命令窗、命令历史窗、当前目录浏览器、工作空间浏览器及一个“启动按钮”。任务栏：位于软件的正上方。各个菜单分别为：文件、编辑、视窗、调试、桌面、窗体、帮助这几个窗口，点击每个窗口可以选择需要的操作。命令窗(

2、CommandWindow)：位于软件操作界面的右侧。在此窗口里，可以输入各种指令、函数、变量表达式并进行各种操作。该窗口用于输入命令并显示除图形以外的所有执行结果。窗口中的“>>”为命令提示符，直接在其后面输入命令并按下回车键后，会出现计算结果在命令后面。命令历史窗(CommandHistory)：位于软件操作界面的左下方。这个窗口记录了命令窗口已经运行过的所有命令（指令、函数等），允许用户对这些命令进行选择、复制。2)MATLAB的函数绘制二维图形最常用的就是plot函数，调用plot函数的三种形式：plot(x)、plot(x,y)、plot(x,

3、y,’r:x’)。还有就是如何添加横坐标和纵坐标标题的命令语句。3.实验程序。w=0.000001:0.0001:0.999999999%定义w的取值范围y=-w.*log2(w)-(1-w).*log2(1-w)%定义二元熵函数的表达式plot(w,y,'r')%画出二元熵函数的曲线图xlabel('w')%x轴的名称ylabel('H(w)')%y轴的名称gridon%给图形加上网格title('二元熵函数H(w)')%函数曲线的名称运行结果如下：一、实验结果分析从图中可以看出熵函数的一些性质，如果二元信源的输出概率是1或0（即二元信源的输出是确定的

4、），则该信源不提供任何信息。当二元信源符号等概率发生时，即w=0.5时，信源的熵达到最大值，等于1比特信息量，曲线关于w=0.5左右对称。二、实验总结对MATLAB掌握不够，还缺少很多的MATLAB知识，应加强学习MATLAB。实验二一般信道容量迭代算法实验报告一、实验目的1.熟悉MATLAB工作环境及工具箱2.掌握一般信道容量迭代算法原理二、实验内容用MATLAB软件编程实现一般信道容量迭代算法三、实验过程1.复习一般信道容量迭代算法，了解其基本思路。2.熟悉MATLAB的工作界面及所要用到的基本函数及语句，如：输入语句、循环语句、exp函数等。3.实

5、验程序。N=input('输入信源符号X的个数N=')M=input('输出信源符号Y的个数M=')p_yx=zeros(N,M)%程序设计需要信道矩阵初始化为零fprintf('输入信道矩阵概率')fori=1:Nforj=1:Mp_yx(i,j)=input('p_yx=');%输入信道矩阵概率ifp_yx(i)<0error('不符合概率分布')endendendfori=1:N%各行概率累加求和s(i)=0;forj=1:Ms(i)=s(i)+p_yx(i,j);endendfori=1:N%判断是否符合概率分布if(s(i)<=0.9999

6、99

7、

8、s(i)>=1.000001)error('不符合概率分布')endendb=input('输入迭代精度：')%输入迭代精度fori=1:Np(i)=1.0/N;%取初始概率为均匀分布endforj=1:M%计算q(j)q(j)=0;fori=1:Nq(j)=q(j)+p(i)*p_yx(i,j);endendfori=1:N%计算a(i)d(i)=0;forj=1:Mif(p_yx(i,j)==0)d(i)=d(i)+0;elsed(i)=d(i)+p_yx(i,j)*log(p_yx(i,j)/q(j));endenda(i)=exp(d(i

9、));endu=0;fori=1:N%计算uu=u+p(i)*a(i);endIL=log2(u)%计算ILIU=log2(max(a))%计算IUn=1while((IU-IL)>=b)%迭代计算fori=1:Np(i)=p(i)*a(i)/u;%重新赋值p(i)endforj=1:M%计算q(j)q(j)=0;fori=1:Nq(j)=q(j)+p(i)*p_yx(i,j);endendfori=1:N%计算a(i)d(i)=0;forj=1:Mif(p_yx(i,j)==0)d(i)=d(i)+0;elsed(i)=d(i)+p_yx(i,j)*l

10、og(p_yx(i,j)/q(j));endenda(i)=exp(d(i));