c语言调matlab fft,C语言及Matlab实现fft

C语言及Matlab实现fft

/*时间抽选基2FFT及IFFT算法C语言实现*/

学号3070503041 电信0702 李春阳

该程序有误输入{1，2，3，4，5，6，7，8}与MATLAB实现不一样，但可运行可以修改修改好请给我留言758677366@http://doc.docsou.com

估计数据排列子程序有误

#include

#include

#include

#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;i

scanf("%lf,%lf",&x[i].real,&x[i].img);

/*检验输入的数据*/

printf("我输入的数据\n");

for(i=0;i

printf("%4lf,+(%4lf)j\n",x[i].real,x[i].img);

initW();

printf("Use FFT(0) or IFFT(1)?\n");