小波分析语音信号并计算各IMF包络的能量熵

小波分析语音信号并计算各IMF包络的能量熵

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）;%重构最后一层小波包系数
HE（i）=wentropy（x‘shannon‘i）;%求小波包能量熵
end
disp（‘小波包能量熵‘）;HE

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

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

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

                 2463                    1