中南大学信息论编码.doc

《中南大学信息论编码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、中南大学信息论编码实验报告题目信息论编码学生姓名指导教师学院学号专业班级完成时间21实验一关于信源熵的实验一、实验目的1.掌握离散信源熵的原理和计算方法。2.熟悉matlab软件的基本操作，练习使用matlab求解信源的信息熵。3.自学图像熵的相关概念，并应用所学知识，使用matlab或其他开发工具求解图像熵。4.掌握Excel的绘图功能，使用Excel绘制散点图、直方图。二、实验原理1.离散信源相关的基本概念、原理和计算公式产生离散信息的信源称为离散信源。离散信源只能产生有限种符号。随机事件的自信息量I（xi）为其对应的随机变量

2、xi出现概率对数的负值。即：I（xi）=-log2p(xi)随机事件X的平均不确定度（信源熵）H（X）为离散随机变量xi出现概率的数学期望，即：H(X)=-∑p(x)I(x)=-∑p(x)logp(x)2.二元信源的信息熵设信源符号集X={0，1}，每个符号发生的概率分别为p(0)=p，p(1)=q，p+q=1，即信源的概率空间为则该二元信源的信源熵为：H(X)=-plogp–qlogq=-plogp–(1-p)log(1-p)即：H(p)=-plogp–(1-p)log(1-p)其中0≤p≤13.MATLAB二维绘图用matla

3、b中的命令plot(x,y)就可以自动绘制出二维图来。例1-2，在matlab上绘制余弦曲线图，y=cosx，其中0≤x≤2p。>>x=0:0.1:2*pi；%生成横坐标向量，使其为0，0.1，0.2，…，6.2>>y=cos(x)；%计算余弦向量>>plot(x,y)%绘制图形4.MATLAB求解离散信源熵求解信息熵过程：1)输入一个离散信源，并检查该信源是否是完备集。212)去除信源中符号分布概率为零的元素。3)根据平均信息量公式，求出离散信源的熵。5.图像熵的相关知识图像熵是一种特征的统计形式，它反映了图像中平均信息量的多少

4、。图像的一维熵表示图像中灰度分布的聚集特征所包含的信息量，令Pi表示图像中灰度值为i的像素所占的比例，则定义灰度图像的一元灰度熵为：图像熵计算过程：1)输入一幅图像，并将其转换成灰度图像。2)统计出图像中每个灰度阶象素概率。3)计算出一幅图像的一维熵。6.Excel的绘图功能比如：用Excel或制作二元熵函数曲线。具体步骤如下：1）启动Excel应用程序。2）准备一组数据p。在Excel的一个工作表的A列（或其它列）输入一组p，取步长为0.01，从0至100产生101个p（利用Excel填充功能）。3）使用Excel的计算功能，在

5、B列中用二元熵函数计算公式，求得A列中各数值对应的二元熵值。比如：在单元格B2中输入公式：=-A2*LOG(A2,2)-(1-A2)*LOG(1-A2,2)。4）使用Excel的图表向导，图表类型选“XY散点图”，子图表类型选“无数据点平滑散点图”，绘制二元熵函数散点图。三、实验内容1、使用matlab软件绘制二元信源熵函数曲线，并说明其物理意义。其程序源代码为：p=0.00000:0.001:1;h=-p.*log2(p)-(1-p).*log2(1-p);plot(p,h);title(‘二元信息熵函数曲线’)；结果为：21因

6、为信源熵表征信源的平均不确定度，而由图表可以看出，当二元信源的符号的发生概率越高或越低时，信源的不确定性越低，反之，则不确定性越高。2、使用matlab软件求解离散单符号信源熵，请自己构造两个信源空间，根据求解结果说明其物理意义。（1）设事件为今天的天气状况为X，有三种情况，分别为晴天，阴天，雨天其概率分别为1/2，1/31/6，则：其源代码为：p1=[1/21/31/6]h=sum(-p1.*log2(p1))结果为：h=1.4591结果表示其事件的信源熵。（2）设学校有四个食堂，去四个食堂的概率分别为1/4，1/8，1/8，1

7、/8，则有：p2=[1/41/81/21/8]h=sum(-p2.*log2(p2))结果为：H=1.75003、使用matlab软件计算图像一维图像熵，请自己选择任意两幅图像，根据求解结果说明其物理意义。其源代码程序为：21lename='mao3.jpg';imfinfo(filename)¡¡¡¡imgRgb=imread(filename);¡¡¡¡imshow(imgRgb);¡¡¡¡A=rgb2gray(imgRgb);¡¡¡¡figure¡¡imshow(A);¡¡¡¡imwrite(A,'mao4.jpg');[M

8、,N]=size(A);fip=zeros(256,1);form=1:M;forn=1:N;ifA(m,n)==0;i=1;elsei=A(m,n);endp(i)=p(i)+1;endendp=p./(M*N)result=0;fori=1:l