作者:cc_vb8 | 来源:互联网 | 2023-05-17 14:14
原理自己查书,代码如下:****************************************************************功能:一维快
原理自己查书,代码如下:
/****************************************************************
功能: 一维快速傅里叶反变换
参数: complex *TD: 指向时域数组的指针
complex *FD: 指向频域数组的指针
int r: 2的幂数,即迭代次数
返回值: 无
***************************************************************/
void Ctry::IFFT(complex *FD, complex *TD, int r)
{
LONG count; //傅里叶变换点数
int i; //循环变量
complex *x;
count = 1 <[count]; //分配运算所需要的存储器
//将频域点写入x
memcpy(x, FD, sizeof(complex)* count);
//求共轭
for (i = 0; i (x[i].real, -x[i].imag());
}
FFT(x, TD, r); //调用快速傅里叶变换
//求时域点共轭
for (i = 0; i (TD[i].real/count, -TD[i].imag()/count);
}
delete x; //释放内存
}
FFT函数见上一篇博客