SE.rar_SE熵_matlab小波熵_信号 熵 特征_小波熵_小波特征

clc clear all x=importdata('D:\Program Files\MATLAB\R2013a\bin\shu.mat'); N=1024; for k=1:20 y=x(:,1); %用db3小波进行5层分解 [c,l]=wavedec(y,5,'db3'); % %构造第1-5层低频系数 matlab388 % a5=wrcoef('a',c,l,'db3',5); % a4=wrcoef('a',c,l,'db3',4); % a3=wrcoef('a',c,l,'db3',3); % a2=wrcoef('a',c,l,'db3',2); % a1=wrcoef('a',c,l,'db3',1); %构造第1-5层高频系数 matlab388 d5=wrcoef('d',c,l,'db3',5); d4=wrcoef('d',c,l,'db3',4); d3=wrcoef('d',c,l,'db3',3); d2=wrcoef('d',c,l,'db3',2); d1=wrcoef('d',c,l,'db3',1) %构造第5层高频细节系数 d5=wrcoef('d',c,l,'db3',5); % %合并为一个矩阵 % A=[a1 a2 a3 a4 a5 d5]; A=[ d1 d2 d3 d4 d5 a5 ]; for t=1:6 for s=1:N T(s)=A(s,t)^2 ; end e(t)= sum(T); end %小波熵 E=sum(e); p=e/E; for i=1:6; h(i)=p(i)*log(p(i)); end H(k)=-sum(h); end