声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1969|回复: 4

[小波] 小波分解信号遇到的问题

[复制链接]
发表于 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');


多模态信号图

多模态信号图
  

多模态信号频谱图

多模态信号频谱图


以上为加速度信号及其频谱图,从频谱图上看只有三个主要的模态频率

分解成的单模态信号

分解成的单模态信号
  

分解成的单模态信号频谱图

分解成的单模态信号频谱图

以上为分解成的单模态信号及其频谱图,频谱图与想要的结果相差很多
理论上应该出现类似下列三幅频谱图,各位帮忙看看问题出在哪里,非常感谢!
101.jpg    201.jpg    301.jpg

ScopeData.mat

20.04 KB, 下载次数: 9

多模态信号

回复
分享到:

使用道具 举报

发表于 2009-12-24 19:55 | 显示全部楼层
分解单分量有问题。
发表于 2011-4-3 12:39 | 显示全部楼层
为什么没人回答呢,菜鸟初学路过~~!
发表于 2012-9-28 21:35 | 显示全部楼层
同问,不知道楼主解决了没?
发表于 2012-9-29 20:33 | 显示全部楼层
正在研究中
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 05:10 , Processed in 0.081118 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表