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

C语言及Matlab实现fft.

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

学号3070503041 电信0702 李春阳

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

估计数据排列子程序有误

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

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=( 1<

for(j=0;j

for(k=0;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]=