基4—FFT程序代码

关于FFT程序，基2的可能能多，但基4—FFT程序代码相对较少，这个基4—FFT程序已经测试通过。

/*时间抽选基4 16点FFT及IFFT算法C语言实现*/
#include “stdio.h“
#include “math.h“
#define N 32

/*定义复数类型*/
typedef struct{
     double real;
     double img;
}complex;

complex x[N] *W; /*输入序列旋转因子*/
int size_x=64;     /*输入序列的大小，在本程序中仅限4的次幂*/
double PI;        /*圆周率*/

int main（）{
     int imethod;
     void fft（）;    /*快速傅里叶变换*/
     void ifft（）;
     void initW（）;  /*初始化变换核*/
     void change（）; /*变址*/
     void add（complex acomplex bcomplex *c）; /*复数加法*/
     void mul（complex acomplex bcomplex *c）; /*复数乘法*/
     void sub（complex acomplex bcomplex *c）; /*复数减法*/

     void output（）;                            /*输出结果*/
     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]:\n“）;
     for（i=0;i           scanf（“%lf%lf“&x.real&x.img）;
//     initW（）;
//     printf（“Use FFT（0） or IFFT（1）?\n“）;
//     scanf（“%d“&method）;
//     if（method==0）
//           fft（）;
//     else
//           ifft（）;
//     output（）;
//     return 0;
}

void add（complex acomplex bcomplex *c）{
     c->real=a.real+b.real;
     c->img=a.img+b.img;
}

void mul（complex acomplex bcomplex *c）{
     c->real=a.real*b.real - a.img*b.img;
     c->img=a.real*b.img + a.img*b.real;
}
void sub（complex acomplex bcomplex *c）{
     c->real=a.real-b.real;
     c->img=a.img-b.img;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1510  2010-05-23 22:17  fft-16-4.c

----------- ---------  ---------- -----  ----

                 1510                    1