%振動信號的小波能譜熵
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