C语言及Matlab实现fft.docVIP

/*时间抽选基2FFT及IFFT算法C语言实现*/ 学号3070503041 电信0702 李春阳 该程序有误输入{1，2，3，4，5，6，7，8}与MATLAB实现不一样，但可运行可以修改修改好请给我留言758677366@ 估计数据排列子程序有误 #include stdio.h #include math.h #include stdlib.h #define N 1000 /*定义复数类型*/ typedef struct{ double real; double img; }complex; complex x[N], *W; /*输入序列,变换核*/ int size_x=0; /*输入序列的大小，在本程序中仅限2的次幂*/ double PI; /*圆周率*/ int main(){ int i,method; void fft(); /*快速傅里叶变换*/ void ifft(); void initW(); /*初始化变换核*/ void change(); /*变址*/ void add(complex a,complex b,complex *c); /*复数加法*/ void mul(complex a,complex b,complex *c); /*复数乘法*/ void sub(complex a,complex b,complex *c); /*复数减法*/ void divi(complex a,complex b,complex *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;isize_x;i++) scanf(%lf,%lf,x[i].real,x[i].img); /*检验输入的数据*/ printf(我输入的数据\n); for(i=0;isize_x;i++) printf(%4lf,+(%4lf)j\n,x[i].real,x[i].img); initW(); printf(Use FFT(0) or IFFT(1)?\n); scanf(%d,method); if(method==0) fft(); else ifft(); output(); return 0; } /*快速傅里叶变换*/ void fft(){ int i=0,j=0,k=0,l=0; complex up,down; //蝶形输出端的上和下 change(); //调整输入的序列，即蝶形的输入端的相应数据排列 for(i=0;i (int)( log(size_x)/log(2) );i++) { l=( 1i ); for(j=0;jsize_x;j+= (1l) ){ /*一组蝶形运算*/ for(k=0;kl;k++){ mul(x[j+k+l],W[size_x*k/2/l],up); add(x[j+k],up,up); mul(x[j+k+l],W[size_x*k/2/l],down); sub(x[j+k],down,down); x[j+k]=up; x[j+k+l]=down