心胸决定格局,眼界决定境界...

梳状滤波器滤除谐波

假设原始信号采样率为4000kHz,由200Hz为基波的一组幅度为1的谐波构成,此处由正弦函数模拟。

由此得到200Hz,400Hz,600Hz...2000Hz,但是此处2000Hz的幅度几乎为0,后面详细说明。

 

200Hz幅度为1的时域波形。 FFT变换之换,对应频率的幅度依然为1,频率依然在200Hz,镜像频率对称,在频率(4000-2000)Hz直接忽略(奈奎斯特定律),这就是傅里叶变换的巧妙之处。400,600,…1800Hz同样如此。

200 

对于2000Hz,奇怪的是还有一些杂波。

2000_1

究其原因,感觉是pi的精度出现了问题,matlab默认为3.1416,于是强制pi=3.1415926;

2000_2

效果好了很多,但是幅度不是为1,不像其它的值,频率幅度不为1,接近0,sin(pi*n),所有的时域值都为0,所以频率对应的值接近0.

最终得到谐波图。

harmonic

应用梳妆滤波器滤波谐波。

comb

滤除后的频谱

filtered

实际和理论总有一些偏差

posted @ 2015-01-19 11:25  WELEN  阅读(1408)  评论(0编辑  收藏  举报