图5.2 程序大致框图
4.FFT通过C++的实现
4.1基2的DIT蝶形算法源代码及注释框图
                              图6.1蝶形运算示意图
/************FFT***********/            //整个程序输入和输出利用同一个空间x[N]，节约空间
  #include   <stdio.h>
  #include   <math.h>
  #include   <stdlib.h>

  #define   N   1000                    //定义输入或者输出空间的最大长度
  typedef struct
  {
       double   real;
       double   img;
  }complex;                              //定义复数型变量的结构体
  void   fft();                            //快速傅里叶变换函数声明
  void   initW();                         //计算W(0)~W(size_x-1)的值函数声明
  void   change();                        //码元位置倒置函数函数声明
  void   add(complex,complex,complex *);   /*复数加法*/  
  void   mul(complex,complex,complex *);   /*复数乘法*/  
  void   sub(complex,complex,complex *);   /*复数减法*/  
  void   divi(complex,complex,complex *);   /*复数除法*/  
  void   output();                       /*输出结果*/

  complex   x[N],*W;                   /*输出序列的值*/
  int       size_x=0;                    /*输入序列的长度，只限2的N次方*/
  double    PI;                         //pi的值

  int   main()
  {
      int   i;

      system("cls");
      PI=atan(1)*4;
      printf("Please input the size of x:\n");                 
                                       /*输入序列的长度*/
      scanf("%d",&size_x);
      printf("Please input the data in x[N]:(such as:5 6)\n");    VC++的FFT快速傅里叶变换编程设计+流程图+源代码(13):http://www.751com.cn/jisuanji/lunwen_954.html