快速傅里叶反变换（一维）

来源：互联网发布：淘宝精选是怎么报名编辑：程序博客网时间：2024/05/16 18:20

原理自己查书，代码如下：

/****************************************************************功能：       一维快速傅里叶反变换参数：       complex<double> *TD： 指向时域数组的指针             complex<double> *FD：  指向频域数组的指针             int r：    2的幂数，即迭代次数返回值：   无***************************************************************/void  Ctry::IFFT(complex<double> *FD, complex<double> *TD, int r){LONG count;     //傅里叶变换点数int       i;            //循环变量complex<double> *x;count = 1 << r;      //计算傅里叶变换点数x = new complex<double>[count];      //分配运算所需要的存储器//将频域点写入xmemcpy(x, FD, sizeof(complex<double>)* count);//求共轭for (i = 0; i < count; i++){x[i] = complex<double>(x[i].real, -x[i].imag());}FFT(x, TD, r);     //调用快速傅里叶变换//求时域点共轭for (i = 0; i < count; i++){TD[i] = complex<double>(TD[i].real/count, -TD[i].imag()/count);}delete x;   //释放内存}

FFT函数见上一篇博客

0 0