基于Android系统的傅立叶变换.docVIP

基于Android系统的傅立叶变换 　　摘要：傅立叶变换在通信、光学、天文方面都有着极其广泛的应用。基于PC平台的研究已经相当成熟，具有较多的工具和库可以使用（如，MATLAB和FFTW开源库）。但是在嵌入式系统上却很难找到相关的库和分析工具，在Android系统上更是如此。然而音频分析和游戏开发等领域傅立叶变换有着不可或缺的作用，所以在嵌入式系统上实现快速高效的傅立叶变换有着实际的应用意义。 　　关键词：傅立叶变换；Android；音频分析 　　中图分类号：TP301 文献标识码：A文章编号：1007-9599 (2011) 09-0000-03 　　Fourier Transform Based on Android System 　　Gu Qiwei 　　(School of Computer Engineering and Science Shanghai University,Shanghai200072,China) 　　Abstract:Fourier transform in the communications,optics,astronomy have a very wide range of applications.PC-based research is already quite mature,with more tools and libraries can be used(eg,MATLAB and open source library FFTW).But in the embedded system is difficult to find the relevant libraries and analysis tools,especially in the Android system.However,audio analysis and Fourier transform the field of game development has an indispensable role,so in the embedded systems to achieve fast and efficient application Fourier transform has a real meaning. 　　Keywords:Fourier transform;Android;Audio analysis 　　傅立叶变换在科学计算中有着极其重要的意义，因为Android系统开放的SDK是基于JAVA语言实现的，如果想通过SDK实现其算法在性能上由于JAVA语言自身的特性很难达到较高的效率。如何解决在Android系统上实现高性能的傅立叶变换是本文要着手解决的问题。 　　一、离散傅立叶分析 　　随着计算机技术的发展与完善。科学研究和工程计算中的很多问题跟计算机已经密不可分了。计算机的一个典型特征是离散化，而傅立叶变换的本质是一个积分计算，体现为连续化特征[2]。同时在实际应用中信号都是离散化采样得到的（例如来自CD的信号），为了能够有效利用计算机实现傅立叶变换的需求，需要对傅立叶级数实现高效、高精度的离散化，为此提出了离散化傅立叶变换（DFT）的概念。由1.3式可得傅立叶变换的离散形式： 　　 ……（2.1） 　　令 可得2.1式的矩阵形式： 　　 ……（2.2） 　　其中， ， ， 　　 　　根据2.2式，程序不难实现，下面给出代码： 　　DFT（input） 　　for k-0 to length（input） 　　for j-0 to length（input） 　　wjk=cos（2*PI/n*k*j）+sin（2*PI/n*k*j）i； 　　Output+=input[j]*wjk 　　将2.1式分为奇、偶两部分则可得到下式： 　　 ，其中 　　令， ， ，，将其代入上式可得： 　　 ……（2.3） 　　上式用矩阵的形式可写成 ， 和 分别表示偶数序列和奇数序列的傅立叶变换。根据以下性质： 　　（1） 和 是以N为周期的序列。 　　（2） 。 　　可将2.3式表示成： 　　 ……（2.4） 　　由2.2式可知，离散傅立叶变换需要与矩阵 相乘，所以2.2式的时间复杂度为 。而2.4式中 和 的时间复杂度均为 ，所以2.4式的时间复杂度为 。由此可见2.4式的算法要比2.2式计算次数少了一半。如果对 和 做同样的分解，便可将计算量继续减半。如此迭代便得到了快速傅立叶算法[3]，最终的时间复杂度为 。 　　下面以序列 为例推导整个快速傅立叶变换过程。 　　首先根据2.4式将2.2式的傅立叶变换序列分成奇、偶两部分。