c代码实现 ifft运算_2D-FFT及IFFT（C语言实现（转载）

前面编过2D-FFT的程序，现在把2D-IFFT的程序整合到一起，便于后面做图像变换反变换使用。

FFT与IFFT有如下关系：

相应的2D-FFT与2D-IFFT的关系如下：

所以可以利用一个FFT核心函数实现2D-FFT与2D-IFFT。代码如下：

#include

#include

#include

#define intsize sizeof(int)

#define complexsize sizeof(complex)

#define PI 3.1415926

int *a,*b;

int nLen,init_nLen,mLen,init_mLen,N,M;

FILE *dataFile;

typedef struct{

float real;

float image;

}complex;

complex *A,*A_In,*W;

complex Add(complex, complex);

complex Sub(complex, complex);

complex Mul(complex, complex);

int calculate_M(int);

void reverse(int,int);

void readData();

void fft(int,int);

void Ifft();

void printResult_fft();

void printResult_Ifft();

int main()

{

int i,j;

readData();

A = (complex *)malloc(complexsize*nLen);

reverse(nLen,N);

for(i=0; i

{

for(j=0; j

{

A[j].real = A_In[i*nLen+b[j]].real;

A[j].image = A_In[i*nLen+b[j]].image;

}

fft(nLen,N);

for(j=0; j

{

A_In[i*nLen+j].real = A[j].real;

A_In[i*nLen+j].image = A[j].image;

}

}

free(a);

free(b);

free(A);

A = (complex *)malloc(complexsize*mLen);

reverse(mLen,M);

for(i=0; i

{

for(j=0; j

{

A[j].real = A_In[b[j]*nLen+i].real;

A[j].image = A_In[b[j]*nLen+i].image;

}

fft(mLen,M);

for(j=0; j

{

A_In[j*nLen+i].real = A[j].real;

A_In[j*nLen+i].image = A[j].image;

}

}

free(A);

printResult_fft();

Ifft();

printResult_Ifft();

return 0;

}

void readData