VC++的FFT快速傅里叶变换编程设计+流程图+源代码(13)

菜单

#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");