求一维序列的信息熵

《求一维序列的信息熵》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、求一维序列的信息爛(香浓爛)的matlab程序实例目录1124555556求一维序列的信息爛(香浓爛)的matlab程序实例弓丨言测试程序额外信息上代码实验一：计算离散信源的爛一、实验设备二、实验冃的三、实验内容四、实验报告要求习题引言对于一个二维信号，比如灰度图像，灰度值的范围是0-255,因此只要据像素灰度值(0-255)出现的概必就可以计算出信息爛。但是，对于一个一维信号，比如说心电信号，数据值的范围并不是确定的，不会是(0-255)这么确定，如果进行域值变换，使其转换到一个整数范用的话，就会丢失数据，请高手指

2、点，怎么计算。比如数字信号是x(n),n=l〜N(1)先用Hist函数对x(n)的赋值范围进行分块，比如赋值范围在0〜10的对应第—•块，10〜20的第二块，以此类推。这之前需要对x(n)做一些归一化处理(2)统计每一•块的数据个数，并求出相应的概率(3)用信息爛公式求解以上求解方法获得的虽然是近似的信息爛，但是一般认为，这么做是没有问题的求一维序列的信息爛的matlab程序代码如下：(已写成调用的函数形式)测试程序%%将下面程序拷贝保存成m文件,直接运行就行,不需要分开保存functioncal_entropy()

3、%主程序clearclccloseallfs=12000;N=12000;T=l/fs;t=(O:N・l产T;ff=104;sig=0・5'¥(1+sin(2忝pi器护f)).^sin(2^pi^3000^t)-¥rand(l,length(t));Hx=yyshang(sigf10)%%求一维离散序列信息爛matlab代码(子函数)functionHx=yyshang(y,duan)%不以原信号为参考的时间域的信号爛%输入:mxf:原信号的能量谱中能量最大的点%护待求信息爛的序列^oduan:待求信息爛的序列要被分

4、块的块数%Hx:y的信息爛%duan=10;%将序列按dug数等分,如果duan=lQ,就将序列分为10等份x_min=min(y);x_max=max(y);maxf(1)=^ibs(xjtnax-xjtnin);maxf(2)-xjnin;duanj—l.O/duan;jiatige=maxf(l)*duanjt;%fori=l:10%pnum(i)=length(find((y_p>-(i~l)^jiange)_p

5、ange));fori=2:duan-lpnum(i)=length(find((y>=maxf(2)+(i-1)^jiange)&(y-maxf(2)^duanA)^jiange));%sum(pnum)ppnum=pnum/sum(pnutn);%每段出现的概率%sum(ppnum)Hx=O;fori-l:duanifppnum(i)==OHi=O;elseHi—-ppnum(i)^log2(ppnum(i));en

6、dHx=Hx+Hi;end额外信息上代码%%为了方便以后查找,在这里把MATLAB中求信息爛的代码写下来,也方便大家查看。clear;clc;closeall;ImageData=iinread(9lena.jpg9);ifndims(ImageData)==3figure;imshow(ImageData);title「您选择的是RGB图,将转换为灰度图!•);ImageData=rgb2gray(ImageData);%如果是真彩色图,就将其装换为灰度图endp=itnhist(ImageData(:));p(p

7、~O)=[];p=pJnumel(ImageData);Hl=-sum(p.Hog2(p))%在这里Y是指的输出图像。H指的是信息爛。实验一：计算离散信源的爛一、实验设备1、计算机2、软件：Matlab二、实验目的1、熟悉离散信源的特点；2、学习仿真离散信源的方法3、学习离散信源平均信息量的计算方法4、熟悉Matlab编程；三、实验内容1、写出计算自信息虽的Matlab程序2、写出计算离散信源平均信息量的Matlab程序。3、掌握二元离散信源的最大信息量与概率的关系。4、将程序在计算机上仿真实现，验证程序的正确性并完

8、成习题。四、实验报告要求简要总结离散信源的特点及离散信源平均信息量的计算，写出习题的MATLAB实现语句。信息论棊础：自信息的计算公式Matlab实现：I=log2(l/p)或I=-log2(p)爛(平均白信息)的计算公式Matlab实现:HX=sum(-x.*log2(x)):或者h=h-x(i)*log2(x(i));习题1.甲地天气预报构