【2018年最新整理】MATLAB应用——信息论论文.doc

青岛农业大学 信息与编码理论实验报告 姓名：韩祖良 班级：信计1201 学号指导老师：辛永训 2015年6月  信息与编码理论实验课：实验一 实验题目 Shannon码的编码 实验目的 掌握MATLAB基本操作；Shannon码的编码 实验地点及时间 信息楼214机房，周一上午1-2节；周三上午3-4节 实验内容 1. MATLAB简介及常用功能； 2．离散信源的MATLAB分析； 3．离散信道的MATLAB分析； 4．Shannon码的编码原理及软件实现。 实验习题 你对MATLAB软件的了解； 答：我参加过数学建模国赛和美赛，获得过省级二等奖和美赛成功参与奖，更期待今年10月份的国赛可以有所突破。我对MATLAB比较熟悉，在建模过程中经常用的是其强大的计算和编程绘图功能。其操作简单，不像c语言等需要抽象的编程语言，使用的大多直接是数学中的符号，而且快捷键很方便，比如输入“clc”命令就可以清屏、输入“plot(x,y)命令”就可以绘制二维图像。编程制图方面，可以直接打开代码本，不用将代码敲入主界面，建模过程中印象比较深刻的是去年国赛A题嫦娥二号着陆点的月球剖面图就可以用MATLAB绘制，其在代码中是一个1000*1000的矩阵，而矩阵数字的突变代表了月球剖面图的高低。同样，B题设计折叠椅的动图模型也是用MATLAB做出来的。下面较为详细的介绍一下MATLAB，MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，代表了当今国际科学计算软件的先进水平。? ?MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连??matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 通过具体实例说明离散信源的MATLAB分析：自信息、信息熵等； 答：注：本例题是本人从网上选择的例题，然后将原题修改而成，程序是参考网上相关问题自己编写。 自信息的计算公式：I=p237 Matlab实现：I=log2(1/p) 或I=-log2(p) 熵（平均自信息）的计算公式p237 Matlab实现：HX=sum(-x.*log2(x))；或者h=h-x(i)*log2(x(i)); 例题： 甲地天气预报构成的信源空间为： X 晴 云 大雨 小雨 乙地信源空间为： Y 晴 小雨 求此两个信源的熵。求各种天气的自信息量。 解： 运行程序： p1=[1/3,1/4,1/6,1/4];%p1代表甲信源对应的概率 p2=[5/8,3/8];%p2代表乙信源对应的概率 H1=0.0; H2=0.0; I=[]; J=[]; for i=1:4 H1=H1+p1(i)*log2(1/p1(i)); I(i)=log2(1/p1(i)); end disp(自信息量分别为：); I disp(H1信源熵为：); H1 for j=1:2 H2=H2+p2(j)*log2(1/p2(j)); J(j)=log2(1/p2(j)); end disp(自信息量分别为：); J disp(H2信源熵为：); H2 运行结果： 自信息量分别为： I = 1.5850 2.0000 2.5850 2.0000 H1信源熵为： H1 = 1.9591 自信息量分别为： J = 0.6781 1.4150 H2信源熵为： H2 = 0.9544 3．通过具体实例说明离散信道的MATLAB分析：信道容量、最佳输入概率分布等； 注:本例题摘自课本p238页。 答：例题：求下图所示的信道容量及其最佳的输入概率分布。 解：本题的MATLAB程序如下： clear all; %清除所有变量 a=0; p=[1/3 1/3 1/6 1/6;1/6 1/6 1/3 1/3]; %信道矩阵P for