傅里叶变换(FFT)的主要思想与算法

简 介: 利用fft算法实现快速傅里叶变换, 在理论、工程中具有非常广泛的应用。除了能够在合适的计算平台完成fft算法,同时还需要注意到它在频谱分析中可能带来的频率混叠以及频率泄露等问题。
关键词: fft,算法实现
01 python算法
  今天下午的信号与系统, 给同学们介绍了离散傅里叶变换的基本应用, 并且介绍了快速傅里叶变换(fft)的主要思想与算法。fft算法因其优异的性能和广泛的应用, 堪称信息处理领域的原子武器。实现fft编程语言很多, 比较来比较去, 利用python语音所描述的该算法最为简明和优雅。
1.1 fft算法代码
  下面的代码是在  the fast fourier transform (fft): most ingenious algorithm ever?  视频中给出的 fft 递归算法形式, 最大精度反映了fft算法核心。
  这个代码实现了dif(时域抽取快速傅里叶变换), 利用递归定义,将fft核心算法中的分而治之体现的淋漓尽致, 突出了递归核心中的核心思想。
def fft(p):    n = len(p)    if n * 1: return p    ye = fft(p[0::2])    yo = fft(p[1::2])    y = [0]*n    w = exp(-1j*2*pi/n)    for j in range(n//2):        yow = w**j * array(yo)        y[j] = ye[j] + yow[j]        y[j+n//2] = ye[j] - yow[j]    return y  
  利用python语音中对于数组切片操作语法, 还可以将上面fft算法中的循环部分都替换成关于数组的操作, 使得实际运算速度得到提高。
def fft1(p):    n = len(p)    if n * 1: return p    ye = fft(p[0::2])    yo = fft(p[1::2])    w = exp(-1j*2*pi/n)**array(list(range(n//2)))    yow = w*yo    y = [0]*n    y[:n//2] = ye + yow    y[n//2:] = ye - yow    return y  
1.2 fft 算法测试
  为了测试算法的有效性, 下面对于一个方波信号计算对应的fft结果。
  测试算法代码如下:
len = 1024onelen=10p1 = [1]*onelen+[0]*(len-onelen)y = fft(p1)plt.plot(abs(array(y)), label='abs(fft)')plt.plot(p1, label='data')plt.xlabel(y)plt.ylabel(abs(fft(y)))plt.grid(true)plt.legend(loc='upper right')plt.tight_layout()plt.show()  
  下面是测试利用python语言实现的fft算法计算结果。
▲ 图1.2.1 利用python语音实现的fft算法测试结果
02 其它语言fft
  fft算法贵在计算效率,前面使用python实现fft,虽然形式上优雅,但实际执行效率不高。提高执行效率,还是需要使用编译语言。
2.1 fortran fft算法
  在我上大学期间所学的编程语言为fortran, 估计现在没有多少同学学习这个算法语言。下面给出了利用fortran语言实现的fft算法程序。
  算法整体上包括有两个阶段:
第一个阶段实现了对输入数据进行倒读顺序排列;
第二阶段利用三重循环实现了分组蝶形运算。
  当然了,时过三十年再看fortran感觉十分酸爽, 但它简练语言和执行高效还是让我们回忆起当年编程时所感觉到的快乐。
▲ 图 fortran 语言实现的fft算法
2.2 c语言fft算法
  下面是在网络上博文 c++ program to compute discrete fourier transform using fast fourier transform approach[1]  给出的fft算法, 没有对其功能进行测试。相比于前面利用python,fortran来看, c语言实现fft就显得非常啰嗦了。
#include #include #include #include using namespace std;unsigned int bitreverse(unsigned int x, int log2n) {    int n = 0;    int mask = 0x1;    for (int i = 0; i < log2n; i++) {        n = 1;    }    return n;}const double pi = 3.1415926536;templatevoid fft(iter_t a, iter_t b, int log2n) {    typedef typename iterator_traits::value_type complex;    const complex j(0, 1);    int n = 1 << log2n;    for (unsigned int i = 0; i < n; ++i) {        b[bitreverse(i, log2n)] = a[i];    }     for (int s = 1; s <= log2n; ++s)  {        int m = 1  1;        complex w(1, 0);        complex wm = exp(-j * (pi / m2));        for (int j = 0; j < m2; ++j) {            for (int k = j; k < n; k += m) {                complex t = w * b[k + m2];                complex u = b[k];                b[k] = u + t;                b[k + m2] = u - t;            }            w *= wm;        }    }}  
※ 总 结 ※
  利用fft算法实现快速傅里叶变换, 在理论、工程中具有非常广泛的应用。除了能够在合适的计算平台完成fft算法,同时还需要注意到它在频谱分析中可能带来的频率混叠以及频率泄露等问题。
原文标题:优雅的fft算法
文章出处:【微信公众号:硬件攻城狮】欢迎添加关注!文章转载请注明出处。


超级电容出现这些问题不能使用
SiO2在氢氟酸中的刻蚀机理
无线防水智能话筒的结构设计是怎样的
导致汽车车门闭锁器无法正常工作存在哪些原因
岂止一箭多星 印度极轨卫星运载火箭要玩一箭103星
傅里叶变换(FFT)的主要思想与算法
vivoX23炫彩版发布 11月23日开卖售价2798元
央行数字货币大概什么时候推出
光程研创将推出全球首款基于锗硅技术的宽谱3D飞行时间传感技术
3D打印为汽车行业带去无限可能 或将实现大规模量产
电动汽车是不是比燃油车更容易自燃
外骨骼机器人发展前景之家政搬家应用
有哪些方法可以解决压力传感器的故障?压力传感器的正确安装方法是什么?
中国智能家居产品与人工智能、物联网等新兴技术结合
联想集团主办的全生态研发中心启动仪式在沈阳举行
AI芯片浪潮:创新才是关键
Silicon Labs利用基于标准的Wi-SUN技术扩展物联网无线产品组合
跨界合作与创新:推动汽车半导体技术快速发展的关键
Nanoleaf引领一个以完全个性化为中心概念的智能家居新时代
水雨情监测系统-实时监测水雨情设备