请问我这样做小波包络谱对吗?(附代码)
clcclear 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]