快速傅立叶变换(FFT)C#.docVIP

快速傅立叶变换（FFT）的C#代码 这个代码是从《快速傅立叶变换（FFT）的C++实现与Matlab实验》这篇文章里的源代码转换而来，请注意查看原文。 在这里自己转换成了C#代码，并作了一些改动，主要是对N值的确定，原文的N值为常量1024，自己通过对输入的数组的长度来确定N值，N值的确定符合2的n次方，函数返回N值。通过作者提供的测试变量进行测试，能得到相应的结果。代码如下： FFT代码： using System;using System.Collections.Generic;using System.Text; namespace ConsoleApplication1{??? /// summary??? /// 快速傅立叶变换(Fast Fourier Transform)。???? /// /summary??? public class TWFFT??? {??????? private TWFFT()??????? {??????? } ??????? private static void bitrp(float[] xreal, float[] ximag, int n)??????? {??????????? // 位反转置换 Bit-reversal Permutation??????????? int i, j, a, b, p; ??????????? for (i = 1, p = 0; i n; i *= 2)??????????? {??????????????? p++;??????????? }??????????? for (i = 0; i n; i++)??????????? {??????????????? a = i;??????????????? b = 0;??????????????? for (j = 0; j p; j++)??????????????? {??????????????????? b = b * 2 + a % 2;??????????????????? a = a / 2;??????????????? }??????????????? if (b i)??????????????? {??????????????????? float t = xreal[i];??????????????????? xreal[i] = xreal[b];??????????????????? xreal[b] = t; ??????????????????? t = ximag[i];??????????????????? ximag[i] = ximag[b];??????????????????? ximag[b] = t;??????????????? }??????????? }??????? } ??????? public static int FFT(float[] xreal, float[] ximag)??????? {??????????? //n值为2的N次方??????????? int n = 2;??????????? while (n = xreal.Length)??????????? {??????????????? n *= 2;??????????? }??????????? n /= 2; ??????????? // 快速傅立叶变换，将复数 x 变换后仍保存在 x 中，xreal, ximag 分别是 x 的实部和虚部??????????? float[] wreal = new float[n / 2];??????????? float[] wimag = new float[n / 2];??????????? float treal, timag, ureal, uimag, arg;??????????? int m, k, j, t, index1, index2; ??????????? bitrp(xreal, ximag, n); ??????????? // 计算 1 的前 n / 2 个 n 次方根的共轭复数 Wj = wreal [j] + i * wimag [j] , j = 0, 1, ... , n / 2 - 1??????????? arg = (float)(-2 * Math.PI / n);??????????? treal = (float)Math.Cos(arg);??????????? timag = (float)Math.Sin(arg);????????