uestcer03 发表于 2010-12-16 22:51

请问我这样做小波包络谱对吗?(附代码)

clc
clear all
load 72HrStress_102_1
s=chanvals(:,1)';

s=(s-mean(s))/std(s,1);
fs=102400;

=wavedec(s,6,'db10');

for i=1:6
    d=wrcoef('d',c,l,'db10',i);
    y=abs(hilbert(d));
Sper=abs(fft(y));
Sper=Sper*2/length(Sper);
fn=(0:length(Sper)-1)*fs/length(Sper);
figure(1),
subplot(6,1,i),
plot(fn(1:length(fn)/2),Sper(1:length(Sper)/2))
axis();
end
页: [1]
查看完整版本: 请问我这样做小波包络谱对吗?(附代码)