基于FPGA离散分数阶傅里叶变换算法的实现.pdfVIP

第2O卷 第2期 河南机电高等专科学校学报 Vo1．20No．2 2012年 3月 JournalofHenanMechanicalandElectricalEn~neefingCollege Mar．2012 基于 FPGA离散分数阶傅里叶变换算法的实现 史兆强 ，魏增辉 ，马 壮。 (1．河南机电高等专科学校 电子通信S-程系，河南 新乡453000； 2．黄河水利职业技术学院信息工程系，河南 开封475003；3．中国电子科技集团公司第二十七研究所，河南 郑州450047) 摘要：为在可编程逻辑门阵列(FPGA)上实现离散分数阶傅里叶变换(DFRFT)，对多种DFRFT算法进行比较分 析，并选用Ozaktas提出的离散快速算法做 了基于FPGA的实现。在对该算法进行理论分析的基础上，对快速算 法的计算过程做优化处理，并给出仿真结果比较。在保证精度要求下，提出详细的FPGA实现方案。工程实际应 用表明：该方法具有优 良性能。 关键词：离散分数阶傅里叶变换；现场可编程逻辑阵列；FFr核 中图分类号：TN911．72 文献标识码 ：A 文章编号：1008—2093(2012)02—0017—03 灵活高效的快速算法。近年来出现的比较成熟的快 1 引言 速离散算法可分三类 J： 分数阶傅里叶变换是对经典傅里叶变换的推广， 1)特征分解型，通过取离散傅里叶变换 (DFT)核 在分析和处理非平稳信号领域有着广泛的应用。矩阵的特征值和特征向量构造 DFT核矩阵的分数幂 1996年Ozaktas等提出了一种计算量与FFT相当的离 来计算DFRF．T。该类方法能保证分数阶傅里叶变换 散算法后，各种分数阶傅里叶变换的离散算法的相继 的大部分性质，但缺点是计算复杂度高，不利于实时 出现，进一步促进了分数阶傅里叶变换在光学信号处 处理。 理、时频分析在内的多个信号处理领域的发展。至今 2)线性组合型，通过对某些特定阶数的DFRFT 为止，有关分数阶傅里叶变换 已具有 比较成熟的快速 的线性组合得到任意阶数 DFRFT，缺点是计算精度与 离散算法，这就为分数阶傅里叶变换能进人数字信号 连续变换相比误差是最大的。 处理的工程实用阶段做了理论基础。近几年的研究 3)离散分解型，由Ozaktas等人提出，通过数学推 文献表明-】J，大多集中在对分数阶傅里叶变换的理论 导将连续FRFT计算转化为一个调制信号与chirp信 研究上，涉及具体工程实现方面甚少。本文主要在对 号相卷积形式，最后利用FFI运算得到分数阶傅里叶 现有 DFRFT算法的分析基础之上，选用 Ozaktas等 变换结果。 人 提出的快速算法，对其做了易于FPGA实现的仿 在 目前已有的各种 DFRFT算法中，离散分解型 真分析，给出可以直接应用在FPGA上的计算实现。 由于计算精度高、计算负担小、便于FPGA实现而成 2 离散分数阶傅里叶变换算法 为主流算法。 3 离散分解方法理论基础 对于连续的分数阶傅里叶变换定义式…如下： ()= []()=S_+-x() (t，“) ， (1) 分解法是根据 FRFI’表达式，将 FRFr分解为信 其 中 号的卷积形式，从而利用 FrI来计算FRFT。这种算 f、 -2 … ≠ 仃 法的机理在进行FRFvr计算之前须先对原始信号进行 ， (￡，)={(一“)， ：2nTr (2) 量纲归一化处理 J。理论上，任何信号和它的傅里叶 【(t+“)， =(2n+1)7r