假设原始信号采样率为4000kHz，由200Hz为基波的一组幅度为1的谐波构成，此处由正弦函数模拟。
由此得到200Hz,400Hz,600Hz...2000Hz，但是此处2000Hz的幅度几乎为0，后面详细说明。
200Hz幅度为1的时域波形。 FFT变换之换，对应频率的幅度依然为1，频率依然在200Hz，镜像频率对称，在频率(4000-2000)Hz直接忽略(奈奎斯特定律)，这就是傅里叶变换的巧妙之处。400，600，…1800Hz同样如此。
对于2000Hz，奇怪的是还有一些杂波。
究其原因，感觉是pi的精度出现了问题，matlab默认为3.1416,于是强制pi=3.1415926;
效果好了很多，但是幅度不是为1，不像其它的值，频率幅度不为1，接近0，sin(pi*n)，所有的时域值都为0，所以频率对应的值接近0.
最终得到谐波图。
应用梳妆滤波器滤波谐波。
滤除后的频谱
实际和理论总有一些偏差