快速傅氏变换之旅(五) 嵌入式中的FFT（最好选择DSP或FPGA）

欢迎大家提意见讨论。

转载请标明是引用于 http://blog.csdn.net/chenyujing1234

例子代码:(编译工具：VS2005)

==================理论==================================================================================

傅立叶变换的重要性不用我说，想必大家也很清楚，有了傅立叶变换，我们就可以从信号的频域特征去分析信号。尤其在无线通信系统中，傅里叶变换的重要性就更加明显了，无论是设计者还是测试工程师，在工作中都会和傅立叶变换打交道。在以下的文章中，我给出一种傅里叶变换的C语言实现方法（参考了C常用算法集），可以用于在嵌入式系统中实现傅立叶变换。

常规的傅立叶变换算法并不适用于嵌入式控制系统，原因是运算量太大（涉及到复数运算），比如离散的傅立叶变换等同于用序列Y(n×1列矢量)乘以n×n 矩阵Fn，需要n×n次乘法。若n=1024,则是104，8576次乘法运算。哇，这么多呀！什么概念呢？如果你选用的CPU单周期指令为25ns, 单周期也可以完成一次乘法运算，那么要计算1024点的傅立叶变换则需要26.2144ms,这还不包括加法或其它运算，对于大多数实时系统，这个处理时间实在太长。于是寻找一个快速的傅立叶变换算法是人们所期望的。

本来我想把FFT的整个数学推导过程列完出来，但当自己硬着头皮看完后，发现对我没有任何用处，我又不是专门研究数学算法的，哪有那么多时间跟着书本的公式去慢慢推导。我想，这些推导问题还是让数学家想去吧。我需要的不过是理解它，然后学会应用它就行。有兴趣的读者可以参考相关的资料，这方面的资料实在太多了。

虽然FFT大幅度地降低了常规傅立叶变换的运算量，但对于一般的单片机或嵌入式设备（如ARM MIPS Linux）而言，处理FFT运算还是力不从心。主要原因是FFT计算过程中的蝶形运算是复数运算，要分开实部和虚部分别计算，想想这是多么繁琐的事情。可能会有些初学者认为，有这么复杂吗？我在PC上使用C++一样可以对复数直接进行加、减、乘、除运算。你说得不错，可以这么做，但那是C++封装了对复数处理的类，直接调用就行。在PC上运算这种类型的算法一般不考虑时间和空间，多一两秒的运行时间不会有什么灾难性的结果。

所以我们要衡量一个处理器有没有足够的能力来运行FFT算法，根据以上的简单介绍可以得出以下两点：

处理器要在一个指令周期能完成乘和累加的工作，因为复数运算要多次查表相乘才能实现。
间接寻址，可以实现增/减1个变址量，方便各种查表方法。FFT要对原始序列进行反序排列，处理器要有反序间接寻址的能力。

所以，在数字信号的分析处理应用中，DSP比其它的处理器有绝对的优势，因为DSP完全具备以上条件。这就是单片机（51系列，AVR，PIC等等）或ARM处理器很少用来进行数字信号分析的原因。

==================实践==================================================================================

在实践的偿试中比如从AD得到的数字信号的采样率为30K,

这么多的数据在800M的机器里做FFT得用700ms.而AD又一直产生，对于实时要求如此之高的应用这显示是不满足要求的.

当然不同的FFT算法可能得到的结果不一样

（算法可以参照我的两篇文章http://blog.csdn.net/chenyujing1234/article/details/7419863 http://blog.csdn.net/chenyujing1234/article/details/7416775)

如果在类似于示波器的模拟信号实时采测系统中，用到的FFT处理速度相当关键，系统所耗的时间为两部分：前端数据采集

（请参考http://blog.csdn.net/chenyujing1234/article/details/7545658），另一部分为FFT运算。

两者加起来时间如果过来就会造成前端采集数据丢失，有些网友讲把数据先存在RAM时，这必给RAM带来压力。

当然在非DSP上做FFT有一个好处就是能进行画像显示。但一屏数据在OnPaint时做一次也要4ms。（第一次初始化时是500ms）.

这只能开发时用到。

综上所述，做FFT运算最好的平台为DSP或FPGA

快速傅氏变换之旅(五) 嵌入式中的FFT（最好选择DSP或FPGA）相关推荐

快速傅氏变换之旅(二) 七种FFT算法速度比较(含代码)
转载请标明是引用于 http://blog.csdn.net/chenyujing1234 例子代码:(编译工具:VS2005) http://www.rayfile.com/zh-cn/files/ ...
数字信号处理学习笔记[0] 连续信号的频谱和傅氏变换
文章目录绪论 1 连续信号的频谱和傅氏变换 1.1 有限区间上连续信号的傅氏级数和离散频谱 1.2 傅氏变换,连续信号与频谱 1.2.3 频谱的基本性质实际应用举例习题绪论 Q: 举例说明&q ...
FFT快速傅式变换算法halcon算子，持续更新
目录 convol_fft convol_gabor correlation_fft deserialize_fft_optimization_data energy_gabor fft_generi ...
【数理方程】傅氏变换拉氏变换
卷积定义: f1(x)f_1(x)f1(x)和f_2(x)都可进行Fourier变换: f1(x)∗f2(x)=∫−∞+∞f1(x−ξ)f2(ξ)dξ=∫−∞+∞f1(ξ)f2(x−ξ)dξf_1 ...
【地震波滤波】保边滤波、傅氏变换干扰波去噪滤波、基于小波分解和重建的干扰波去噪、基于维纳滤波的去噪、中值滤波、视速度滤波
1.保边滤波保护边缘滤波器,通常有四种类型,其中性能较为优良的是双边滤波器,其主要原理为: 双边滤波方法(Bilateral filtering)是基于Gsuss滤波方法提出的,主要是针对Gauss ...
matlab z变换离散化_大学学的傅里叶变换、拉氏变换、z变换，这些还能搞得懂不？...
1.关于傅里叶变换变换? 答:fourier变换是将连续的时间域信号转变到频率域:它可以说是laplace变换的特例,laplace变换是fourier变换的推广,存在条件比fourier变换要宽,是 ...
matlab拉氏逆变换,偏微分方程的拉氏变换解法.ppt
偏微分方程的拉氏变换解法第15章积分变换的MATLAB求解编者 Outline 15.1 傅里叶变换 15.2 拉普拉斯变换 15.3 Z 变换 15.1 傅里叶变换 1.傅里叶变换的概念若函 ...
离散哈特莱变换(DHT)及快速哈特莱变换(FHT)学习
离散哈特莱变换(DHT)及快速哈特莱变换(FHT)学习说在前边最近复习\(DSP\)的时候,发现了一个号称专门针对离散实序列的变换,经分析总运算量为普通\(FFT\)的几乎一半,而且完全没有复数. ...
数学实验8 用matlab软件求级数的和_函数的泰勒级数和傅氏级数,数学实验8用matlab软件求级数的和函数的泰勒级数和傅氏级数.pdf-资源下载蚂蚁文库...
数学实验8 用matlab软件求级数的和函数的泰勒级数和傅氏级数.pdf 数学实验八用Matlab软件求级数的和.函数的泰勒级数和傅氏级数一.求级数的和在Matlab中,可用symsum函数求数列或级 ...

快速傅氏变换之旅(五) 嵌入式中的FFT（最好选择DSP或FPGA）

快速傅氏变换之旅(五) 嵌入式中的FFT（最好选择DSP或FPGA）相关推荐

最新文章

热门文章