clear?all;
M=load（‘32.txt‘）;
s1=M;
subplot（921）?
plot（s1）?
title（‘原始信號‘）?
ylabel（‘S1‘）
wpt=wpdec（s13‘db1‘‘shannon‘）;
%plot（wpt）;?
s130=wprcoef（wpt[30]）;?
s131=wprcoef（wpt[31]）;?
s132=wprcoef（wpt[32]）;?
s133=wprcoef（wpt[33]）;?
s134=wprcoef（wpt[34]）;?
s135=wprcoef（wpt[35]）;?
s136=wprcoef（wpt[36]）;?
s137=wprcoef（wpt[37]）;?
s10=norm（s130）;?
s11=norm（s131）;?
s12=norm（s132）;?
s13=norm（s133）;?
s14=norm（s134）;?
s15=norm（s135）;?
s16=norm（s136）;?
s17=norm（s137）;
st10=std（s130）;?
st11=std（s131）;?
st12=std（s132）;?
st13=std（s133）;?
st14=std（s134）;?
st15=std（s135）;?
st16=std（s136）;?
st17=std（s137）;
disp（‘正常信號的特征向量‘）;?
snorm1=[s10s11s12s13s14s15s16s17]?
std1=[st10st11st12st13st14st15st16st17]
subplot（923）;plot（s130）;?
ylabel（‘S130‘）;?
subplot（925）;plot（s131）;?
ylabel（‘S131‘）;?
subplot（927）;plot（s132）;?
ylabel（‘S132‘）;?
subplot（929）;plot（s133）;?
ylabel（‘S133‘）;?
subplot（9211）;plot（s134）;?
ylabel（‘S134‘）;?
subplot（9213）;plot（s135）;?
ylabel（‘S135‘）;?
subplot（9215）;plot（s136）;?
ylabel（‘S136‘）;?
subplot（9217）;plot（s137）;?
ylabel（‘S137‘）;
%fft
y1=fft（s11024）;?
py1=y1.*conj（y1）/1024;?
y130=fft（s1301024）;?
py130=y130.*conj（y130）/1024;?
y131=fft（s1311024）;?
py131=y131.*conj（y131）/1024;?
y132=fft（s1321024）;?
py132=y132.*conj（y132）/1024;?
y133=fft（s1331024）;?
py133=y133.*conj（y133）/1024;?
y134=fft（s1341024）;?
py134=y134.*conj（y134）/1024;?
y135=fft（s1351024）;?
py135=y135.*conj（y135）/1024;?
y136=fft（s1361024）;?
py136=y136.*conj（y136）/1024;?
y137=fft（s1371024）;?
py137=y137.*conj（y137）/1024;
f=1000*（0:511）/1024;?
subplot（922）;?
plot（fpy1（1:512））;?
ylabel（‘P1‘）;?
title（‘原始信號的功率譜‘）?
subplot（924）;?
plot（fpy130（1:512））;?
ylabel（‘P130‘）;?
subplot（926）;?
plot（fpy131（1:512））;?
ylabel（‘P131‘）;?
subplot（928）;?
plot（fpy132（1:512））;?
ylabel（‘P132‘）;?
subplot（9210）;?
plot（fpy133（1:512））;?
ylabel（‘P133‘）;?
subplot（9212）;?
plot（fpy134（1:512））;?
ylabel（‘P134‘）;?
subplot（9214）;?
plot（fpy135（1:512））;?
ylabel（‘P135‘）;?
subplot（9216）;?
plot（fpy136（1:512））;?
ylabel（‘P136‘）;?
subplot（9218）;?
plot（fpy137（1:512））;?
ylabel（‘P137‘）;?
%?plottree（wpt）

%%%%%%%%%%%求能量熵%%%%%%%%%%%%%%%%%%%
[T]=wpdec（s13‘db1‘‘shannon‘）;%小波包分解
for?i=1:2^3
x=wprcoef（Ti+2^3-2）;%重構(gòu)最后一層小波包系數(shù)
HE（i）=wentropy（x‘shannon‘i）;%求小波包能量熵
end
disp（‘小波包能量熵‘）;HE

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2463??2014-07-16?16:37??wavenengliang.m

-----------?---------??----------?-----??----

?????????????????2463????????????????????1