小波熵matlab程序

clc
clear all
x=importdata（‘C:\Documents and Settings\Administrator\桌面\参考文献\数据\N均匀磨损.txt‘）;
N=1024;
for k=1:20
y=x（:1）;
%用db3小波进行5层分解
[cl]=wavedec（y5‘db3‘）;

% %构造第1-5层低频系数 matlab388
% a5=wrcoef（‘a‘cl‘db3‘5）;
% a4=wrcoef（‘a‘cl‘db3‘4）;
% a3=wrcoef（‘a‘cl‘db3‘3）;
% a2=wrcoef（‘a‘cl‘db3‘2）;
% a1=wrcoef（‘a‘cl‘db3‘1）;


%构造第1-5层高频系数 matlab388
d5=wrcoef（‘d‘cl‘db3‘5）;
d4=