ifft的c语言编程

#include 
#include 


void IFFT（double x[]double y[]int n）//数组x存储时域序列的实部，数组y存储时域序列的虚部n代表N点FFT.
{
int ijklmn1n2;
double cc1ess1ttrti;


for（j = 1i = 1; i<16; i++）//计算 i = log2N；变量j进行每次自乘2，当j=n，即算出i的值.
{
     m = i;
     j = 2*j;
     if（（j == n））//此时判断j是否与n相等时，j与2乘，所以n输入时应该选为2的次方。
         break;
}


n1 = n - 1; //计算蝶形图的输入下标（码位倒读）
for（j=0i=0; i{
    if（i    {
         tr = x[j];
         ti = y[j];
         x[j] = x[i];
         y[j] = y[i]; 
         x[i] = tr;
         y[i] = ti;                 
    }
    k = n/2;
    while（k<（j+1））
    {
        j = j - k;
        k = k/2;              
    }
    j = j + k;
}


n1 = 1;  //计算每一级的输出，l为某一级，i为同一级的不同群，使用同一内存（即位运算）
for（l=1; l<=m; l++）
{
    n1 = 2*n1;
    n2 = n1/2;
    e = 3.1415926/n2;
    c = 1.0;
    s = 0.0;
    c1 = cos（e）;
    s1 = sin（e）;
    for（j=0; j    {
        for（i=j; i