小波熵程序

matlab中的小波熵程序

%振动信号的小波能谱熵
clear;
clc;
close all;
%
% 
load zc-nlh-10.mat;
y0=Data1_AI_4（1:160000）; %内圈裂纹故障有160000个数据，分成40组，每组4000个数据
y1=reshape（y0400040）;  %分成40组，每组4000个数据
y1=y1‘;
y=y1（1:）;
fs=10000;
t=（（1:length（y））/fs）‘; 
figure;
h = gca; % 获取当前绘图坐标的指针
set（h‘FontSize‘12）; % 设置文字大小，同时影响坐标轴标注、图例、标题等。
plot（ty）;
xlabel（‘\sl{t/s}‘）;
ylabel（‘\sl{a（m/s^2）}‘）;
title（‘原始信号‘）;
set（gca‘ylim‘[-150150]）;


%计算40组样本各自的小波能谱熵得到40x7的矩阵WEE
for i=1:40
   WEE（i:）=wentropy0（y1（i:））;
end


%组合样本160x7
x1=load（‘zc.mat‘）;
x2=load（‘nlh.mat‘）;
x3=load（‘sd.mat‘）;
x4=load（‘dz.mat‘）;
xx=[x1.WEE;x2.WEE;x3.WEE;x4.WEE];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        730  2012-11-17 15:25  小波熵程序\Waveletentropy.m

     文件       1689  2012-11-25 15:09  小波熵程序\wentropy0.m

     目录          0  2012-12-26 09:44  小波熵程序

----------- ---------  ---------- -----  ----

                 2419                    3