一种直接计算实序列DFT的基2DIT-FFT算法.docVIP

中围电子学会第十四届信息论学术年会，２００７年１２月１５．１８日，广州 一种直接计算实序列ＤＦＴ的基２ ＤＩＴ－ＦＦＴ算法 李海军严朝军彭文标 三峡大学电气信息学院 宜昌４４ ３００２ 摘要：本文钊‘对Ｎ＝２。点实序列的ＤＦＴ提出了一种直接按ＦＦＴ运算流图进行计算的算法。该算法由Ｎ个存 储单元存储Ｎ点实序列ＤＦＴ和Ｉ蝶形运算网络算法两部分组成。利用本文提出的算法计算Ｎ点实序列ＤＦＴ时， 只需要Ｎ个存储单元，计算总量不到基２ ＦＦＴ算法计算量的一半。 关键词：实序列ＤＦＴ基２ ＦＦＴ；算法 Ａｎ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｃｏｍｐｕｔｉｎｇ ｄｉｒｅｃｔｌｙ ＤＦＴ ｏｆ ｒｅａｌ ｓｅｑｕｅｎｃｅ ｗｉｔｈ ｒａｄｉｘ．２ ＤＩＴ．ＦＦＴ Ａｂｓｔｒａｃｔ：Ａｎ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｃｏｍｐｕｔｉｎｇ ｄｉｒｅｃｔｌｙ Ｎ＝２‘”一ｐｏｉｎｔ ＤＦＴ ｏｆ ｒｅａｌ—ｖａｌｕｅｄ ｓｅｑｕｅｎｃｅ ｗｉｔｈ ｔｈｅ ｆｌｏｗ．ｇｒａｐｈ ｏｆ ｒａｄｉｘ一２ ＤＩＴ－ＦＦＴ ａｌｇｏｒｉｔｈｍ ｉＳ ｄｅｖｅｌｏｐｅｄ．Ｔｈｅ ａｌｇｏｒｉｔｈｍ ｉＳ ｃｏｍｐｏｓｅｄ ｏｆ ｔｈｅ ｐｒｉｎｃｉｐｌｅ ｏｆ Ｎ ｍｅｍｏｒｙ ｌｏｃａｔｉｏｎｓ ｆｏｒ ｓｔｏｒｉｎｇ Ｎ—ｐｏｉｎｔ ＤＦＴ ａｎｄ ａｎ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｃｏｍｐｕｔｉｎｇ ｂｕｔｔｅｒｆｌｙ ｎｅｔｗｏｒｋ．Ｉｔ ｅｍｐｌｏｙｓ Ｎ ｍｅｍｏｒｙ ｌｏｃａｔｉｏｎｓ ａｎｄ ｌｅｓｓ ｔｈａｎ ｈａｌｆ ｃｏｍｐｕｔａｔｉｏｎ ａｍｏｕｎｔ ｏｆ ｒａｄｉｘ－２ ＦＦ Ｔ． Ｋｅｙｗｏｒｄ：Ｒｅａｌ Ｓｅｑｕｅｎｃｅ ＤＦＴ Ｒａｄｉｘ一２ ＦＦＴ ａｌｇｏｒｉｔｈｍ  一、引言 ＤＦＴ是信号分析与处理中的一种重要变换。由于直接计算ＤＦＴ的运算量很大，影响了其 实际应用，因此研究ＤＦＴ的快速算法是有意义的。自从Ｃｏｏｌｅｙ－Ｔｕｋｅｙ的傅里叶变换快速算 法提出后，新的算法不停涌现，有影响的算法主要有基２、基４、实因子、分裂基、素因子 等算法。这些算法都是针对复数形式的输入数据提出的，而输入数据为实序列时，很多文献， 如ｆ１１『２１ｆ３１，认为有三种处理方法。一是用一个Ｎ点ＦＦＴ计算两个Ｎ点实序列的ＤＦＴ，二是 用Ｎ／２点ＦＦＴ计算一个Ｎ点实序列的ＤＦＴ，三是用离散哈特莱变换（ＤＨＴ）的快速算法（ＦＨＴ） 得到实信号的ＤＨＴ谱。这三种方法不能直接得出Ｎ点ＤＦＴ的结果，要经过转换才能得到。 近些年来，有不少论文对实序列ＤＦＴ的快速算法进行了探讨，如文献【４１讨论了实序列ＤＦＴ 的基２ ＦＦＴ递推算法，该算法在减少乘法次数的同时增加了加法次数；文献ｆ５１提出了实序列 傅里Ｉ｝｜．变换定义及其分裂基快速算法，但其规律性并不是很强。本文提出的算法，以规律性 最强、最易于理解的基２ ＦＦＴ运算流图为基础，直接计算实序列ＤＦＴ，实现该算法的程序简 洁明了，所需存储单元个数与ＤＦＴ点数相同，实乘与实加次数总和不到基２ ＦＦＴ运算量总和 的一半。 二、利用Ｎ个存储单元存储Ｎ点实序列ＤＦＴ的原理 众所周知，实序列ｘ（ｎ）的Ｎ点ＤＦＴ具有共轭对称性，用公式表示为ｘ（ｋ）＝ｘ＋（Ｎ—ｋ），Ｏ≤ｋ ≤Ｎ一１。根据共轭对称的特点，可以得到Ｒｅ［Ｘ（ｋ）】＿Ｒｅ［Ｘ（Ｎ—ｋ）］，Ｉｍ［Ｘ（ｋ）】＝Ｉｍ［Ｘ（Ｎ—ｋ）］。由此我 们可以推出：实序列ｘ（ｎ）的Ｎ点ＤＦＴ应该可以用Ｎ个存储单元来存储，其中Ｎ／２个用来存储 实部，Ｎ／２个用来存储虚部。下面以Ｎ＝８点为例对Ｎ个存储单元存储Ｎ点实序列ＤＦＴ的原理 进行推导。 Ｎ＝８点实序列的ＤＦＴ公式为：（ｔ）：童工（打）时，七：ｏ Ｊ，…，７ （１） 为方便推导，我们将式（１）改用图１所示的矩阵（２）来描述。下面我们分三步对该原理进行 推导。 中国电子学会第十四届信息论学术年会，２００７年１２月１５．１８日，广州  Ｏ  Ｄ  ０ ： ３ ｛ 。Ｗ １ 坩蚶Ｗ ４ ６ 。ⅥＪ ｌ』１ ２ Ｊ ０ １：４ ５ ６ ８ 旺 ９ Ｕ 比 。Ⅵ％Ｗ ｏ Ｗ ｗ。ｍ  。吣Ｗ％Ｊ Ｗ Ｊ １ Ｊ ｄ ｌ ０ ３ Ｗ吖吖Ｗ  ６： 比 ¨ Ｈ 珀  ㈣…㈦㈣㈨㈣㈣㈣  闵阻阿阻瞰瞰冈瞰 ∞”∞”曲印回” Ｗ Ｊ Ｋ Ｗ ｗ ｗ １ Ｊ  Ｊ％ｗ№ 之 ３  Ｋ。ｗ—ｍ ｗ盯Ｊ％％Ｗ ２ ｌ  怙ｊ＂忆旷＝兰ｊ”忆旷‰肚  ㈣㈩㈣㈥㈥㈣㈦㈣ 图１矩阵（２） 图２矩阵（３） 第一步：对矩阵（２）进行化简。根据旋转因子ＷＮ＝ｅ－Ｊ２删的特性，矩