C# 实现 FFT 正反变换 和 频域滤波

要进行FFT运算首先要构造复数类，参考

http://blog.csdn.net/iamoyjj/archive/2009/05/15/4190089.aspx

 

下面的程序在依赖上述复数类的基础上实现了FFT正反变换算法和频域滤波算法，另外由于一般如果是对实数进行FFT的话，要将FFT得到的复数数组转为实数数组，下面类中的Cmp2Mdl方法的作用就是这个。这个FFT算法是基-2FFT算法，因此，如入的序列必须是2的n次方个点长。

频域滤波的基本原理是：

1、  对输入序列进行FFT

2、  得到的频谱乘以一个权函数（滤波器，系统的传递函数）

3、  得到的结果进行IFFT

4、  如果是实数运算的话用Cmp2Mdl方法转为实数

代码如下：

    /// <summary>

    /// 频率分析器

    /// </summary>

    public static class FreqAnalyzer

    {

 

        /// <summary>

        /// 求复数complex数组的模modul数组

        /// </summary>

        /// <param name="input">复数数组</param>

        /// <returns>模数组</returns>

        public static double[] Cmp2Mdl(Complex[] input)

        {

            ///有输入数组的长度确定输出数组的长度

            double[] output = new double[input.Length];

 

            ///对所有输入复数求模

            for (int i = 0; i < input.Length; i++)

            {