C语言实现 IFFT 运算

参考《C》C语言实现FFT算法_杨贵安的博客-CSDN博客_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; i++){j = m / 2;is = 2 * is + (m - 2 * j);m = j;}fr[it] = pr[is];fi[it] = pi[is];}pr[0] = 1.0;pi[0] = 0.0;p = 6.283185306 / (1.0 * n);pr[1] = cos(p); //将w=e^-j2pi/n用欧拉公式表示pi[1] = -sin(p);for (i = 2; i <= n - 1; i++)  //计算pr[]{p = pr[i - 1] * pr[1];q = pi[i - 1] * pi[1];s = (pr[i - 1] + pi[i - 1]) * (pr[1] + pi[1]);pr[i] = p - q; pi[i] = s - p - q;}for (it = 0; it <= n - 2; it = it + 2){vr = fr[it];vi = fi[it];fr[it] = vr + fr[it + 1];fi[it] = vi + fi[it + 1];fr[it + 1] = vr - fr[it + 1];fi[it + 1] = vi - fi[it + 1];}m = n / 2;nv = 2;for (l0 = k - 2; l0 >= 0; l0--) //蝴蝶操作{m = m / 2;nv = 2 * nv;for (it = 0; it <= (m - 1) * nv; it = it + nv)for (j = 0; j <= (nv / 2) - 1; j++){p = pr[m * j] * fr[it + j + nv / 2];q = pi[m * j] * fi[it + j + nv / 2];s = pr[m * j] + pi[m * j];s = s * (fr[it + j + nv / 2] + fi[it + j + nv / 2]);poddr = p - q;poddi = s - p - q;fr[it + j + nv / 2] = fr[it + j] - poddr;fi[it + j + nv / 2] = fi[it + j] - poddi;fr[it + j] = fr[it + j] + poddr;fi[it + j] = fi[it + j] + poddi;}}for (i = 0; i < n; i++)  // 取共轭{fr[i] = fr[i] / n;fi[i] = -fi[i] / n;}for (i = 0; i <= n - 1; i++){pr[i] = sqrt(fr[i] * fr[i] + fi[i] * fi[i]);  //幅值计算}return;
}

C语言实现 IFFT 运算相关推荐

c代码实现 ifft运算_X^n+1=0上的FFT和IFFT（基2）——C语言实现
我们一般意义上学习的FFT都是基于的,即FFT中的单位根我们取的是 ,但是在某些情况下我们需要上的FFT和IFFT变换. 1.直接想到的思路是把的根替换成的根. 解法: 的根可以使用的2n个 ...
c代码实现 ifft运算_fft算法c语言_matlab fft算法_ifft c语言
FFT快速算法C程序_工学_高等教育_教育专区.电子信息工程综合课程设计报告书 DSP 课程设计报告题学目: 院: FFT 快速算法 C 程序计算机与信息工程学院 09 ... fft算法代码 ...
c语言指针的相关运算,C语言指针的运算
本文讨论使用指针进行的运算,最重要的运算是获取指针所引用的对象或函数.也可以比较指针,使用指针来遍历一个内存区域. 使用指针读取与修改对象间接运算符 * 生成内存中的位置,该位置的地址存储在一个指针 ...
c语言枚举,位运算
c语言枚举,位运算枚举枚举类型的定义和枚举变量的说明一.枚举的定义枚举类型定义的一般形式为二.枚举变量的说明枚举类型变量的赋值和使用位运算位域一.位域的定义和位域变量的说明二.位域 ...
单片机c语言位运算写法,单片机与嵌入式系统中C语言的位运算小记
原标题:单片机与嵌入式系统中C语言的位运算小记编了个PIC的项目,对里的还是蛮有感悟的,特此记录一下. 譬如说,在程序中定义了一个char类型的变量, purge_short_enable_flag ...
C语言实现乘方运算---m的n次方(附完整源码)
C语言实现乘方运算---m的n次方 int power(int m, int n) {int i = 1;int result = 1;if (n == 0)
c语言集合运算有菜单,怎样用C语言实现集合运算
[转]STL算法-合并(集合) hechao322514692017-02-24 R:集合运算 thoixy29792014-10-27 C++ STL 集合运算 weixin_42993054605 ...
c语言程序设计施莹答案,C语言课件-位运算.ppt
C语言课件-位运算.ppt C语言程序设计,信息工程学院施莹 2019年11月,第12章位运算,C语言是为描述系统而设计的,其功能应当具有汇编的功能故C语言具有高级语言的的功能,也肯定低级语言的 ...
谭浩强c语言不讲位运算呢,谭浩强C语言_CHAR12位运算.DOC
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp计算机&nbsp>&nbspC/C++资料谭浩强C语言_CHAR12位运算.DOC11页 ...

C语言实现 IFFT 运算

C语言实现 IFFT 运算相关推荐

最新文章

热门文章