请问如何进行FFT滤波

只看该作者 · 2017-9-6 15:22

FFT变换后，去掉不想要的频率，再IFFT,请问怎么能去掉不想要的频率？有固定公式吗？

2万 · 2017-9-6 15:55

“FFT变换后，去掉不想要的频率，再IFFT,请问怎么能去掉不想要的频率？”

既然频谱都有了，相应频率位置上将其“挖掉”便是。

“有固定公式吗？”

没有，看你要求。

只看该作者 · 2017-9-6 16:09

HWM 发表于 2017-9-6 15:55
“FFT变换后，去掉不想要的频率，再IFFT,请问怎么能去掉不想要的频率？”

既然频谱都有了，相应频率位置上 ...

谢谢，佩服高手，自己懂得太少了，另外看资料说加什么窗的，FIR和 IIR 滤波好像就得加窗，FFT也需要加窗吗？

2万 · 2017-9-6 17:30

"FFT也需要加窗吗？"

需要。

1万 · 2017-9-6 18:34

本帖最后由雪山飞狐D 于 2017-9-6 18:37 编辑

不是FFT以后随便挖掉再IFFT的，你查下卷积在时域与频域对应含义，时域上怎么做的滤波对应频域该怎么做，有对应的表达式

只看该作者 · 2017-9-6 20:02

雪山飞狐D 发表于 2017-9-6 18:34
不是FFT以后随便挖掉再IFFT的，你查下卷积在时域与频域对应含义，时域上怎么做的滤波对应频域该怎么做，有 ...

之前看过FFT后用“filter_resultR=dataR*(1/(1+pow(i/6,4)))”这个表达式，据说是二阶巴特沃斯滤波，请问表达式是不是不是固定的，都有哪几种能给几个例子吗？

只看该作者 · 2017-9-6 20:06

a3724395 发表于 2017-9-6 20:02
之前看过FFT后用“filter_resultR=dataR*(1/(1+pow(i/6,4)))”这个表达式，据说是二阶巴特沃斯滤波，请问 ...

之前看过FFT后用
for(i=0;i<SAMPLENUMBER;i++)
{
filter_resultR=dataR*(1/(1+pow(i/6,4)));//2阶的bw滤波器,wc=6Hz
filter_resultI=dataI*(1/(1+pow(i/6,4)));//2阶的bw滤波器,wc=6Hz
}
这个表达式，据说是截止频率6HZ的二阶巴特沃斯滤波，请问表达式是不是不是固定的，都有哪几种能给几个例子吗？

只看该作者 · 2017-9-6 20:13

HWM 发表于 2017-9-6 17:30
"FFT也需要加窗吗？"

需要。

比如说用个矩形窗，原始信号FFT后每个点的频率是{1HZ,2HZ,3HZ.......10HZ}，我想要3--6HZ的数据，那矩形窗的幅度和横坐标应该怎么选择才能互相对应上

2万 · 2017-9-6 20:21

a3724395 发表于 2017-9-6 20:13
比如说用个矩形窗，原始信号FFT后每个点的频率是{1HZ,2HZ,3HZ.......10HZ}，我想要3--6HZ的数据，那矩形 ...

矩形窗就是直接取窗内数据。

只看该作者 · 2017-9-6 20:24

HWM 发表于 2017-9-6 20:21
矩形窗就是直接取窗内数据。

“雪山飞狐D”说的卷积，您觉得对FFT滤波有关联吗

2万 · 2017-9-6 20:33

a3724395 发表于 2017-9-6 20:24
“雪山飞狐D”说的卷积，您觉得对FFT滤波有关联吗

你提到的“巴特沃斯滤波”是针对时域物理可实现的滤波器而言，既然你都说了频谱，那么直接在频谱上做**就可以了。

至于那位所说的，是时域和频域的关系。具体该如何，看你是想在时域内还是在频域内玩了。

只看该作者 · 2017-9-6 20:42

HWM 发表于 2017-9-6 20:33
你提到的“巴特沃斯滤波”是针对时域物理可实现的滤波器而言，既然你都说了频谱，那么直接在频谱上做** ...

个人感觉频域比时域好，想学习下频域的知识，我可以不用理频域和时域的具体关系吗，直接挖掉应该简单，重点是掌握FFT的变换？

只看该作者 · 2017-9-6 20:44

HWM 发表于 2017-9-6 20:33
你提到的“巴特沃斯滤波”是针对时域物理可实现的滤波器而言，既然你都说了频谱，那么直接在频谱上做** ...

另外矩形窗是不是应该有个表达式呢？这方面实在不太了解

1万 · 2017-9-6 23:33

本帖最后由雪山飞狐D 于 2017-9-6 23:53 编辑

a3724395 发表于 2017-9-6 20:02
之前看过FFT后用“filter_resultR=dataR*(1/(1+pow(i/6,4)))”这个表达式，据说是二阶巴特沃斯滤波，请问 ...

时域的滤波器是得到巴特沃斯滤波器的单位冲击响应表达式，用这个冲击响应跟时域波形进行卷积即可，FFT滤波就是滤波器单位冲击响应再做FFT，得到这个冲击响应的FFT再频域里面跟信号的FFT做乘法，最后再IFFT，里面的数据处理细节比较多，另外不能在频域里面“乱切”，要符合相关滤波器的响应式，不然的话，最终效果的平坦度还有相位等是会乱套，在论坛只能说大致原理和方向，具体的你可以找书看看，都有的。。。《数字信号处理（第三版）》第447页 [Understanding Digital Signal Processing, Third Edition ]