快速傅立叶变换—算法及应用

来源 :湖南大学 | 被引量 : 3次 | 上传用户:afuren1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
离散傅立叶变换(DFT)广泛应用于几乎所有科学和工程领域。在信号处理领域中,傅立叶变换是最重要的信号处理工具。在通信领域,傅立叶变换大量应用于正交频分复用(OFDM)系统中。然而,直接计算傅立叶变换是一个计算强度很大的工作。所以,高效计算强度小的计算算法是必需的。快速傅立叶变换(FFT)就是一类高效计算强度小的傅立叶变换的计算算法。麻省理工的Johnson和Frigo教授提出的修正的分裂基FFT(MSRFFT)算法是2014之前计算复杂度最低的算法(2014我们提出2个计算复杂度更低的算法)。为了提高该算法的实用性,通过分析MSRFFT,我们发现MSRFFT有一组分解不是必需的。另外,通过适当的变换,可以减少查表访问旋转因子的次数。因此,我们对该算法进行了完善,将该算法的4个分解过程只保留了3个(但计算复杂度保持不变),将MSRFFT所需要的5/8N-2次因子访问次数减少到只要15/32N-2次。针对长度为N= 6m的DFT,我们提出了基-3/6FFT算法来对其进行计算分解,在计算分解的同时对带因子的子DFT进行旋转变换,使得旋转因子的下标不包括因子3,从而减少算法的计算复杂度。针对长度为q×2m的DFT,我们通过执行2m个长度为q的DFT将该DFT分解成q个长度为2m的子DFT。该方法通过旋转变换将带因子(wnN,wNN/4+n,...,WN(q-1)N/4+n的(共计q0项放在一列组成一个长度为q的子DFT,抽取该q个旋转因子的公共因子(wnN使得该子DFT变成带常数因子的DFT(SDFT)。SDFT将会使用我们提出的方法来高效实现,减少其计算复杂度。由于长度为q的FFT的计算精度比2的幂次方的FFT差,我们对以上算法进行改进,通过旋转变换将带因子为Wnq的(共计2m个)项放在一列组成一个长度为2m的子DFT,该SDFT同样可以高效实现,其计算复杂度和原算法一样,精确度提高了。针对长度为N= 2m的DFT,我们提出了四个算法。这四个算法,都是在MSRFFT算法的基础上提出来的。我们的算法克服了 MSRFFT算法的缺点,将MSRFFT算法的因子抽取方法扩展到复数。相比于分裂基FFT(SRFFT),我们所提四个计算2m长度FFT算法中的其中二个算法在计算复杂度、计算精度和旋转系数计算或查表次数三个方面的性能得到了提高。其它两个计算2m长度FFT算法是目前已经出版的算法中计算复杂度最低的算法。由于L-型的蝶的算法具有较高的实现复杂度,因此我们提出了一个新的实现方法,将一个L-型蝶分解成几个类似于基-2 FFT的蝶的执行单元(BU-2),将一个L-型蝶块分解成几个BU-2块。所有分裂基系列FFT算法都可以使用这个方法象Cooley-Tukey算法一样递推实现。采用这个方法,我们用递推的方法实现了第2章的简化的修正的MSRFFT算法和第6章所提的计算长度N=q× 2m的FFT算法。通过比较我们发现:(1)递推实现的MSFFT在DFT较小比其它算法速度快,但当DFT较大时,所需时间会会急剧增加,因为MSFFT需要的内存比其它算法要多。(2)递推实现的所提计算长度N = q × 2m的FFT算法,与其它算法相比所需要的时间相对要少;和FFTW相比,当DFT较小时,我们所提算法要比FFTW快,但当DFT较大时,我们的算法就比FFTW慢,因为我们的程序没有进行内存优化。剪切FFT就是针对输入信号数目和/或输出信号数目远小于DFT长度的一类优化算法。该类算法采用减少或消除标准FFT的多余操作的方式来达到提高效率的目的。考虑到剪切FFT的蝶不再对称的这一特点,我们提出四种基于共轭因子的方法来进行FFT剪切,减少剪切FFT算法在输入和/或输出阶段的计算复杂度。和其它剪切FFT算法相比,我们所提四种算法减少了计算复杂度。多序列比对(MSA)是生物信号学中最重要的分析工具之一。有一个叫MAFT的MSA程序采用FFT来进行序列相似性的计算,识别相似区域,减少序列比对的解空间,是序列比对精度最高的MSA程序中速度最快的程序。为了进一步完善MAFFT的计算效率,我们对MAFFT的序列相似性的计算方法进行了三个方面的修正。首先,基于复数的氨基酸和核苷酸表示被使用在修正的相关系数计算方案中,代替了原来基于实数的表示。由此,相比于原MAFFT方案,一半的内存和一多半的计算可以节省。其次,线性卷积代替了循环卷积来计算氨基酸和核苷酸队列的相关系数,在计算最优路径的过程中,更多的同态区域将会被发现。再其次,我们设计一个FFT算法来高效计算线性卷积。所设计FFT算法基于共轭FFT(CPFFT),不需要执行队列的排序。该FFT还是一个新颖的剪切FFT,他的输出只需要实数。模拟结果显示修正的方案计算队列的相关系数的速度是MAFFT的3倍。总之,我们提出的这几个快速算法是同类算法计算复杂度最低的算法。提出的执行方法可以将所有分裂基FFT算法进行递推执行。提出的修正的MAFFT算法,将MAFFT的速度提高了2~3倍。
其他文献
目的:研究青少年男子业余篮球运动员比赛训练过程中的血清肌酸激酶(CK)和血尿素的变化规律,探索用肌酸激酶和血尿素监控青少年男子业余篮球运动员的比赛训练。方法:10名青少年业余
目的:考评《初中学生生活质量评定量表》(Quality of Life Scale for Junior High School Students,JHSS-QOLS)的可行性、效度和信度。方法:通过随机抽取的278名初中学生的生活质
要想解决"中小学生课业负担重"的老问题,就要从国家、社会、学校三个层面对减负工作进行审视,分析过重课业负担产生的原因,才能对症下药。减负并非一刀切就能达到增效的目的,
音乐课程是初中阶段必不可少的教育课程,对学生审美情趣的培养具有很重要的作用。艺术的感染和熏陶不仅能够激发学生的创造能力,还能够提升学生的思维灵敏度。良好的初中音乐课
棉花是抗旱性较强的作物,但是干旱问题仍然是制约棉花种植发展的重要因素之一。为了应对日益严重的干旱灾害,本研究以17份棉花资源品种为试验材料,在花铃期进行控水胁迫试验,
目的比较空心钉与钢板内固定治疗老年后踝骨折的疗效。方法将87例老年后踝骨折患者按内固定方式分为空心钉组(39例)及钢板组(48例)。术后3周比较两组并发症情况。术后12个月
本文试图从异化与归化的角度来探讨《红楼梦》茶文化的英译策略。
目的:干扰素刺激基因15(Interferon Stimulating Gene 15,ISG15)蛋白,是一种由干扰素刺激基因所编码的小分子物质。研究表明ISG15具有广泛的生理学效应,参与炎症、细胞凋亡和
目的:探讨吸入麻醉药七氟醚对记忆巩固的影响及PICK1/ERK1/2/GSK3β信号通路在其中所起的作用。方法:动物实验应用C57BL/6 PICK1基因敲除小鼠进行七氟醚处理(2%,3h)。提取动