小波分解信号遇到的问题

ytsfxyzyb · 发表于 2009-12-23 15:31

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本人想用小波分析进行模态参数辨识，对多模态实信号先进行小波分解为单模态信号，再对单模态信号小波脊提取进行模态参数辨识。各位前辈帮我看看这个方法理论上可行否？现在在小波分解信号上存在一些问题，程序及运行结果图如下：
clear;
load ScopeData;
t=ScopeData.time;
Delta=0.005;
Fs=1/Delta; % 采样频率
x=ScopeData.signals.values;

figure;
plot(t,x);
xlabel('time（s）');
ylabel('acc（m/s^2）');
figure;
nfft=512;
Y=fft(x,nfft);
f=Fs*(0:(nfft/2-1))/nfft;
fftY=abs(Y)/nfft;
plot(f,fftY(1:nfft/2)),grid %IMF分量的频域图
xlabel('Frequence');ylabel('Amp');

wpt=wpdec(x,3,'db4','shannon');

s130=wprcoef(wpt,[3,0]);
s131=wprcoef(wpt,[3,1]);
s132=wprcoef(wpt,[3,2]);
s133=wprcoef(wpt,[3,3]);
s134=wprcoef(wpt,[3,4]);
s135=wprcoef(wpt,[3,5]);
s136=wprcoef(wpt,[3,6]);
s137=wprcoef(wpt,[3,7]);
figure
subplot(4,2,1);plot(t,s130);
subplot(4,2,2);plot(t,s131);
subplot(4,2,3);plot(t,s132);
subplot(4,2,4);plot(t,s133);
subplot(4,2,5);plot(t,s134);
subplot(4,2,6);plot(t,s135);
subplot(4,2,7);plot(t,s136);
subplot(4,2,8);plot(t,s137);

figure;
nfft=512;
f=Fs*(0:(nfft/2-1))/nfft;
Y0=fft(s130,nfft);
fftY0=abs(Y0)/nfft;
Y1=fft(s131,nfft);
fftY1=abs(Y1)/nfft;
Y2=fft(s132,nfft);
fftY2=abs(Y2)/nfft;
Y3=fft(s133,nfft);
fftY3=abs(Y3)/nfft;
Y4=fft(s134,nfft);
fftY4=abs(Y4)/nfft;
Y5=fft(s135,nfft);
fftY5=abs(Y5)/nfft;
Y6=fft(s136,nfft);
fftY6=abs(Y6)/nfft;
Y7=fft(s137,nfft);
fftY7=abs(Y7)/nfft;

subplot(4,2,1);
plot(f,fftY0(1:nfft/2));
ylabel('P130');
subplot(4,2,2);
plot(f,fftY1(1:nfft/2));
ylabel('P131');
subplot(4,2,3);
plot(f,fftY2(1:nfft/2));
ylabel('P132');
subplot(4,2,4);
plot(f,fftY3(1:nfft/2));
ylabel('P133');
subplot(4,2,5);
plot(f,fftY4(1:nfft/2));
ylabel('P134');
subplot(4,2,6);
plot(f,fftY5(1:nfft/2));
ylabel('P135');
subplot(4,2,7);
plot(f,fftY6(1:nfft/2));
ylabel('P136');
subplot(4,2,8);
plot(f,fftY7(1:nfft/2));
ylabel('P137');

多模态信号图