matlab基2时间抽选法,按时间抽取的基2FFT算法分析及MATLAB实现

电子技术研发ＥｌｅｃｔｒｏｎｉｃｓＲ＆Ｄ

电一子一技一术～

按时问抽取的基２ＦＦＴ算法分析及ＭＡＴＬＡＢ实现

张登奇李宏民李丹

(湖南理工学院信息与通信工程学院)

摘要：ＤＦＴ是一种应用广泛的数学变换工具，ＭＡＴＬＡＢ是～款功能强大的科学计算语言。ＭＡＴＬＡＢ提供的ＦＦＴ函数解决了ＤＦＴ的快速计算问题，但由于它是内建函数而不能了解到软件实现的过程。文章以按时间抽取的基２ＦＦＴ算法为例，根据快速傅里叶变换的原理和规律，绘出了算法实现的程序框图，列出ＴＩＶＩＡＴＬ＆Ｂ环境下软件实现的程序，建立了从算法理论到程序实现的完整概念。

关键词：数字信号处理；离散傅里叶变换；快速傅里叶变换：按时间抽取；ＭＡＴＬＰｄ３ＡｎａｌｙｓｉｓｏｆＤｅｃｉｍａｔｉｏｎ．ｉｎ．ｔｉｍｅＲａｄｉｘ２ＦＦＴＡｌｇｏｒｉｔｈｍａｎｄ

ＺｈａｎｇＤｅｎｇｑｉＬｉＨｏｎｇｒｎｉｎＬｉＤａｎ

ＭＡＴＬＡＢ

ＢａｓｅｄＲｅａｌｉｚａｔｉｏｎ

(Ｃｏｌｌｅｇｅ

ｏｆＩｎｆｏｒｍａｔｉｏｎ＆ＣｏｍｍｕｎｉｃａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＨｕｎａｎＩｎｓｔｉｔｕｔｅｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ)

Ａｂｓｔｒａｃｔ：ＤＦＴｉＳａｗｉｄｅｌｙｕｓｅｄｍａｔｈｅｍａｔｉｃａＩｔｒａｎｓｆｏｌ＇ｎｌｔｏｏｌ。ａｎｄＭＡＴＬＡＢｉＳａｐｏｗｅｒｆｕｌｓｃｉｅｎｔｉｆｉｃｃｏｍｐｕｔｉｎｇｌａｎｇｕａｇｅ．ＴｈｅＦＦＴｆｕｎｃｔｉｏｎｏｆＭＡＴＬＡＢｃａｎｓｏｌｖｅｔｈｅｐｒｏｂｌｅｍｓｏｆｆａｓｔｃａｌｃｕｌａｔｉｏｎｏｆＤＦＴ，ｂｕｔｉｔｃａｎ’ｔｕｎｄｅｒｓｔａｎｄｔｈｅｐｒｏｃｅｓｓｏｆｓｏｆｔｗａｒｅｒｅａｌｉｚａｔｉｏｎｂｅｃａｕｓｅｉｔｉｓａｂｕｉｌｔ－ｉｎｆｕｎｃｔｉｏｎ．Ｔａｋｉｎｇｔｈｅａｌｇｏｒｉｔｈｍｏｆｔｈｅｒａｄｉｘ一２ＤＩＴＦＦＴａｓｔｈｅｅｘａｍｐｌｅａｎｄｂａｓｅｄｏｎｔｈｅｐｒｉｎｃｉｐｌｅａｎｄｒｕｌｅｓｏｆＦＦＴ，ｔｈｉｓｐａｐｅｒｄｒａｗｓｔｈｅｂｌｏｃｋｄｉａｇｒａｍｓｏｆｔｈｅａｌｇｏｒｉｔｈｍｒｅａｌｉｚａｔｉｏｎ，ｌｉｓｔｓｔｈｅｐｒｏｃｅｄｕｒｅｓｏｆｔｈｅｓｏｆｔｗａｒｅｒｅａｌｉｚａｔｉｏｎｕｎｄｅｒＭＡＴＬＡＢｅｎｖｉｒｏｎｍｅｎｔ．ａｎｄａｌｓｏｅｓｔａｂｌｉｓｈｅｓｔｈｅｃｏｍｐｌｅｔｅｃｏｎｃｅｐｔｆｒｏｍａｌｇｏｒｉｔｈｍｔｈｅｏｒｉｅｓｔｏｐｒｏｇｒａｍｒｅａｌｉｚａｔｉｏｎ．

Ｋｅｙｗｏｒｄｓ：ｄｉｇｉｔａｌｓｉｇｎａｌ

ｐｒｏｃｅｓｓｉｎｇ(ＤＳＰ)；ｄｉｓｃｒｅｔｅ

Ｆｏｕｒｉｅｒ

ｔｒａｎｓｆｏｒｍ(ＤＦＴ)；ｆａｓｔＦｏｕｒｉｅｒｔｒａｎｓｆｏｒｍ(ＦＦＴ)；ｄｅｃｉｍａｔｉｏｎｉｎ

ｔｉｍｅ(ＤＩＴ)：Ⅳ队ＴＬＡＢ

０引言

离散傅里变换(ＤＦＴ)是数字信号处理中重要的数学变换工具，自１９６５年库利和图基提出基２ＦＦＴ算法［１］以来，现已有多种快速算法，且还在不断地进行研究和探索。本文以按时间抽取的基２ＦＦＴ(ＤＩＴ－ＦＦＴ)算法为例，介绍ＦＦＴ算法的基本原理，归纳ＦＦＴ算法的运算规律，总结ＦＦＴ算法的编程框图，列出ＭＡＴＬＡＢ环境下软件实现的程序，从而建立起一种从算法理论到程序实现的全过程的完整概念。

１

主：珲嚣裟嚣篙∞

图１

ＤＩＴ蝶形运算流图符号

２

ＤＩＴ－ＦＦＴ算法的运算规律及编程思想椭聃蛹椭聃鹏埘册

图２

舳枷均珊聃栅珊册

８点ＤＩＴ－ＦＦｒ运算流程

ＤＩＴ－ＦＦＴ算法的基本原理。

有限长序ｙｌＪｘ(ｎ)的Ｎ，点ＤＦＴ定义为：ｚ(｝)‘薹。(”)瞄”，式

中％：ｅ１等，其整数次幂简称为旋转因子。直接进行ＤＦＴ运算大约需要２Ｎ２次三角函数计算，４Ⅳ２次实数乘法计算和２Ｎ(２Ｎ—１)次实数加法计算，且需许多索引和寻址操作［２Ｉ。文献［３】列出了直接ＤＦＴ的ＭＡＴＬＡＢ程序，这种直接ＤＦＴ运算概念清楚、编程简单，但占用内存大、运算速度低，在实际工作中并不实用。基２ＦＦＴ算法的基本思想是把原始的Ⅳ点序列依次分解成一系列短序列，充分利用旋转因子的周期性和对称性，分别求出这些短序列对应的ＤＦＴ，再进行适当的组合，得到原ＪｊＶ点序列的ＤＦＴ，最终达到减少运算次数，提高运算速度的目的。按时间抽取的基２ＦＦＴ算法，先是将Ⅳ点输入序列ｘ(胛)在时域按奇偶次序分解成２个Ｎ／２点序列Ｘ，仰)和Ｘ，∽)，再分别进行ＤＦＴ运算，求出与之对应的ｘｌ(ｋ)和Ｘ２(ｋ)，然后利用图ｌ所示的运算流程进行蝶形运算，得到原Ｎ点序列的ＤＦＴ。只要Ⅳ是２的整数次幂，这种分解就可一直进行下去，直到其ＤＦＴ就是本身的ｌ点时域序列。一个完整的８点ＤＩＴ－ＦＦＴ运算流程如图２所示［４］。图中的输入序列不再是顺序排列但有规律可循，数组Ａ(存储地址)用于存放输入数据和每级运算的结果。

为了编写ＤＩＴ－ＦＦＴ算法的运算程序，首先要分析其运算规律，总结编程思想并绘出程序框图。由图２可知，ＤＩＴ－ＦＦＴ算法的运算过程很有规律。

２．１原位计算

对Ｎ：２ｕ点的ＦＦＴ共进行肘级运算，每级由Ⅳ／２个蝶形运算组成。在同一级中，每个蝶的输入数据只对本蝶有用，且输出节点与输入节点在同一水平线上，这就意味着每算完一个蝶后，所得数据可立即存入原输入数据所占用的数组元素(存储单元)，这种原位(址)计算的方法可节省大量内存。

２．２蝶形运算

实现ＦＦＴ运算的核心是蝶形运算，找出蝶形运算的规律是编程的基础。蝶形运算是分级进行的；每级的蝶形运算可以按旋转因子的指数大小排序进行；如果指数

万方数据

７５