参考 《C》C语言实现FFT算法_杨贵安的博客-优快云博客_c语言 fft 实现了 fft 运算;
如果要实现 ifft 运算,只需要在输入和输出时取共轭就好了,故 ifft 的 C 语言实现如下(已经通过matlab 验证)。
// n 为信号长度, n=2^k
// pr 为输入信号实部,pi 为输入信号虚部
// fr 为输出信号实部,fi 为输出信号虚部
// 函数执行后,pr 变为输出信号的模
void ifft(int n, int k,
long double pr[], long double pi[], long double fr[], long double fi[])
{
int it, m, is, i, j, nv, l0;
long double p, q, s, vr, vi, poddr, poddi;
for (i = 0; i < n; i++) // 取共轭
{
pi[i] = -pi[i];
}
for (it = 0; it <= n - 1; it++) //将pr[0]和pi[0]循环赋值给fr[]和fi[]
{
m = it;
is = 0;
for (i = 0; i <= k - 1;