快速傅里叶变换FFT

快速傅里叶变换FFT

DFT是信号分析与处理中的一种重要变换。但直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大时，计算量太大，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。
1.直接计算DFT
长度为N的有限长序列x(n)的DFT为：

2.减少运算量的思路和方法
思路：N点DFT的复乘次数等于N2。把N点DFT分解为几个较短的DFT，可使乘法次数大大减少。另外，旋转因子WmN具有周期性和对称性。
(考虑x(n)为复数序列的一般情况，对某一个k值，直接按上式计算X(k)值需要N次复数乘法、(N-1)次复数加法.)
方法：
分解N为较小值：把序列分解为几个较短的序列，分别计算其DFT值，可使乘法次数大大减少；
利用旋转因子WNk的周期性、对称性进行合并、归类处理，以减少DFT的运算次数。
周期性：
对称性：
3.FFT算法思想
不断地把长序列的DFT分解成几个短序列的DFT，并利用旋转因子的周期性和对称性来减少DFT的运算次数。








再次分解，对N=8点，可分解三次。


c语言程序:

// FFT.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<windows.h>
#define N 1000
typedef struct{
    double real;
    double img;
 }