声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1402|回复: 4

[滤波] 经FFT变换后,想要某一指定频段的频谱图

[复制链接]
发表于 2010-4-28 17:03 | 显示全部楼层 |阅读模式

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

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

x
一段振动信号,fs=200,N=24511,绘制时域图后,经过FFT变换后得其频谱图,如果我只想要其中25~27HZ频段的频谱图(其他频段均为0),同时反算出25~27HZ对应的时域图改怎么弄呢~
已写程序如下,望高手指点,
fs=200
N=24511
n=0:N-1
t=n/fs
load D4
s=D4
figure(2);
plot(t,s);
xlabel('t');
ylabel('y');
Y=fft(s)
Py=abs(Y)*2/N
f=(0:N/2-1)/N*fs
figure(3);
plot(f,Py(1:N/2))
回复
分享到:

使用道具 举报

发表于 2010-4-28 18:53 | 显示全部楼层
加个band pass filter试试
 楼主| 发表于 2010-4-28 21:04 | 显示全部楼层

具体该怎么弄呢,请指点。我自己弄了个不知道对不对,请指点
Y=fft(s)
y2=zeros(1,N)
Py=abs(Y)
f=(0:N/2-1)/N*fs
for ds=0:N-1
    if  ds/N*fs)<25
        ds/N*fs)>30
        y2(ds+1)=0
    else
        y2(ds+1)=Y(ds+1)
    end,
end
figure(2);
plot(f,Py(1:N/2))
figure(3);
plot(f,y2(1:N/2))
figure(4);
plot(t,real(ifft(y2)))
发表于 2010-4-28 22:13 | 显示全部楼层
help filter
help butter
...
发表于 2010-4-28 23:23 | 显示全部楼层
用FIR滤波器,频段范围设为25~27Hz即可
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 18:33 , Processed in 0.059063 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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