上傳人：飛***

認證信息

認證類型：個人認證

認證主體：常**（實名認證）

IP屬地：河北

IP屬地：河北 上傳時間：2020-05-18 格式：DOC 頁數(shù)：4 大小：17.50KB 積分：15 舉報 版權(quán)申訴

求一維序列的信息熵（香濃熵）的matlab程序?qū)嵗龑τ谝粋€二維信號，比如灰度圖像，灰度值的范圍是0-255，因此只要根據(jù)像素灰度值(0-255)出現(xiàn)的概率，就可以計算出信息熵。但是，對于一個一維信號，比如說心電信號，數(shù)據(jù)值的范圍并不是確定的，不會是(0-255)這么確定，如果進行域值變換，使其轉(zhuǎn)換到一個整數(shù)范圍的話，就會丟失數(shù)據(jù)，請高手指點，怎么計算。比如數(shù)字信號是x(n),n=1N(1)先用Hist函數(shù)對x(n)的賦值范圍進行分塊，比如賦值范圍在010的對應(yīng)第一塊，1020的第二塊，以此類推。這之前需要對x(n)做一些歸一化處理(2)統(tǒng)計每一塊的數(shù)據(jù)個數(shù)，并求出相應(yīng)的概率(3)用信息熵公式求解以上求解方法獲得的雖然是近似的信息熵，但是一般認為，這么做是沒有問題的求一維序列的信息熵的matlab程序代碼如下：（已寫成調(diào)用的函數(shù)形式）測試程序：fs=12000;N=12000;T=1/fs;t=(0:N-1)*T;ff=104;sig=0.5*(1+sin(2*pi*ff*t).*sin(2*pi*3000*t)+rand(1,length(t);Hx=yyshang(sig,10)%求一維離散序列信息熵matlab代碼function Hx=yyshang(y,duan)%不以原信號為參考的時間域的信號熵%輸入：maxf:原信號的能量譜中能量最大的點%y:待求信息熵的序列%duan:待求信息熵的序列要被分塊的塊數(shù)%Hx:y的信息熵%duan=10;%將序列按duan數(shù)等分，如果duan=10,就將序列分為10等份x_min=min(y);x_max=max(y);maxf(1)=abs(x_max-x_min);maxf(2)=x_min;duan_t=1.0/duan;jiange=maxf(1)*duan_t;% for i=1:10% pnum(i)=length(find(y_p=(i-1)*jiange)&(y_pi*jiange);% endpnum(1)=length(find(y=maxf(2)+(i-1)*jiange)&(y=maxf(2)+(duan-1)*jiange);%sum(pnum)ppnum=pnum/sum(pnum);%每段出現(xiàn)的概率%sum(ppnum)Hx=0;for i=1:duanif ppnum(i)=0Hi=0;elseHi=-ppnum(i)*log2(ppnum(i);endHx=Hx+Hi;endend%-擴展閱讀：實驗一：計算離散信源的熵一、實驗設(shè)備:1、計算機2、軟件：Matlab二、實驗?zāi)康?1、熟悉離散信源的特點；2、學(xué)習(xí)仿真離散信源的方法3、學(xué)習(xí)離散信源平均信息量的計算方法4、熟悉Matlab編程；三、實驗內(nèi)容:1、寫出計算自信息量的Matlab程序2、寫出計算離散信源平均信息量的Matlab程序。3、掌握二元離散信源的最大信息量與概率的關(guān)系。4、將程序在計算機上仿真實現(xiàn)，驗證程序的正確性并完成習(xí)題。四、實驗報告要求簡要總結(jié)離散信源的特點及離散信源平均信息量的計算,寫出習(xí)題的MATLAB實現(xiàn)語句。信息論基礎(chǔ)：自信息的計算公式Matlab實現(xiàn)：I=log2(1/p)或I=-log2(p)熵（平均自信息）的計算公式Matlab實現(xiàn)：HX=sum(-x.*log2(x)；或者h=h-x(i)*log2(x(i);習(xí)題：1.甲地天氣預(yù)報構(gòu)成的信源空間為：乙地信源空間為：求此兩個信源的熵。求各種天氣的自信息量。案：運行程序：p1=1/2,1/4,1/8,1/8;%p1代表甲信源對應(yīng)的概率p2=7/8,1/8;%p2代表乙信源對應(yīng)的概率H1=0.0;H2=0.0;I=;J=;for i=1:4H1=H1+p1(i)*log2(1/p1(i);I(i)=log2(1/p1(i);enddisp(自信息量分別為：);Idisp(H1信源熵為：);H1fo