声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2078|回复: 7

[综合] 滤波后EMD分解

[复制链接]
发表于 2010-7-31 16:42 | 显示全部楼层 |阅读模式

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

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

x
我这有一个振动台试验数据,要进行模态识别,我对数据先进行了带通滤波处理,在进行EMD分解,但出现问题,滤波后得到的时域信号EMD分解不出来 IMF=emd(y)   得到的imf 和y一某一样,程序如下:
clear all
clc
fid=fopen('A7e1s1h.dat');        %读入数据,原始数据文件名在程序中输入
x=fscanf(fid,'%f',[1,inf]);
sf=255;                               %sf为采样频率
fmin=4;                             %根据傅立叶谱来判断各阶模态频率的范围,选取合适的滤波频率,单位:Hz
fmax=4.2;
status=fclose(fid);
n=length(x);
t=(0:1/sf:(n-1)/sf);                  %建立离散输出时间向量
nfft=2^nextpow2(n);                    % 建立FFT长度
ni=round(fmin*nfft/sf+1);           %四舍五入取整求最小截止频率对应数组元素的下标
na=round(fmax*nfft/sf+1);           %四舍五入取整求最大截止频率对应数组元素的下标
y=fft(x,nfft);                      %FFT变换,结果存于y
a=zeros(1,nfft);
a(ni:na)=y(ni:na);                   %将y的正频率带通内的元素赋值给a
a(nfft-na+1:nfft-ni+1)=y(nfft-na+1:nfft-ni+1);%将y的负频率带通内的元素赋值给a
y=ifft(a,nfft);                      %将滤波处理后数据恢复出时域信号
y=(real(y(1:n)));                    %取逆变换的实部n个元素为滤波结果列向量
subplot(2,1,1)                       %绘制滤波前的时程曲线图形
figure(1)
plot(t,x)
subplot(2,1,2)                       %绘制滤波后的时程曲线图形
plot(t,y)

%%%以上为输入数据和滤波过程%%%


imf=emd(y);            %调用EMD函数,得到滤波后的各阶IMF
emd_visu(y,1:length(y),imf)

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2010-7-31 16:44 | 显示全部楼层
请哪位高手给解答一下  谢谢啦

评分

1

查看全部评分

 楼主| 发表于 2010-7-31 16:46 | 显示全部楼层
A7e1s1h.dat (44.63 KB, 下载次数: 17)

评分

1

查看全部评分

发表于 2010-11-16 15:55 | 显示全部楼层
这个是 选择范围的问题吧,你试下改为
fmin=3.5;                             %根据傅立叶谱来判断各阶模态频率的范围,选取合适的滤波频率,单位:Hz
fmax=4.5;
发表于 2010-11-16 17:07 | 显示全部楼层
还有,楼主你用的是什么滤波方法啊?
发表于 2010-11-21 20:29 | 显示全部楼层
楼主用的好像是   频域滤波   
发表于 2010-11-21 20:32 | 显示全部楼层
再有,楼主之所以会出现这样的问题,主要是你的截止频率间隔设定的太小了,以至于EMD会认为所得到的滤波信号就是一个单分量信号,能否将其截止频率间隔设大点再试试呢
发表于 2012-11-5 20:10 | 显示全部楼层
我要分析一组振动信号,可是存在噪声,请问如何先进行滤波,在进行EMD?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 04:59 , Processed in 0.079342 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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