离散傅里叶变换 - 快速计算方法及C实现 - 第二篇

最新推荐文章于 2022-09-26 22:03:44 发布

VIP文章 yu_xianguo

最新推荐文章于 2022-09-26 22:03:44 发布

阅读量1k

点赞数 2

分类专栏：编程文章标签： DFT FFT radix-2

本文链接：https://blog.csdn.net/YU_Xianguo/article/details/99591980

版权

DFT – Fast algorithms and C implementations - Part2

Radix-2 DFT

在我的上一篇博客里，已经介绍了基-2 离散傅里叶变换的基本原理（参见第4(4)部分），总结一下就是：

设N长序列 $\lbrace x_j \rbrace$ 的DFT为N长序列 $\lbrace X_k \rbrace$ ，若N被2整除，则：

$\forall ~k<N/2:~ \begin{cases} X_k = \sum_{j}{x_{2j}W_{N/2}^{jk}} + W_N^k \sum_{j}{x_{2j+1}W_{N/2}^{jk}}\\ X_{k+N/2} = \sum_{j}{x_{2j}W_{N/2}^{jk}} - W_N^k \sum_{j}{x_{2j+1}W_{N/2}^{jk}} \end{cases}$

radix-2算法将一个N长的DFT分解为两个N/2长的DFT的线性组合。如果N/2仍然可以被2整除，则可以继续对两个N/2长的子序列进行奇偶分解。如果序列长度为2的指数倍，即

$\exists m \in \lbrace 0,1,2,...\rbrace, s.t., N=2^m$

则可以应用radix-2算法将N长序列逐渐分解，直至序列长度变成1为止。由于标量的DFT等于其自身（见我的上一篇博客），因此我们实际上已经不需要执行任何一次DFT变换——我们只需要执行几次radix-2变换即可。

下面从360图书馆借用两张图来说明这个计算过程：

上图显示的是N=16的DFT计算过程：基于radix-2算法，首先将序列分成偶数列和奇数列，由于两个子序列长度为8，可以进一步对各自进行奇偶分解，经过4次分解以后，我们得到了16组长度为1的子序列，这些序列的DFT等于其自身。下面我们先根据radix-2算法求取8个N=2长序列[0,8], [4,2], [2,10], [6,14], [1,9], [5,13], [3,11], [7,15]的DFT，e.g.,

$\mathcal{F} \left( \begin{bmatrix} x_0 \\ x_8 \end{bmatrix} \right ) = \begin{bmatrix} x_0+W_2^0 x_8 \\ x_0-W_2^0 x_8 \end{bmatrix}$

然后再根据radix-2算法求取4个N=4长序列[0,4,8,12], [2,6,10,14], [1,5,9,13], [3,7,11,15]的DFT，e.g.,

最低0.47元/天解锁文章

yu_xianguo

关注

2
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
离散傅里叶变换 - 快速计算方法及C实现 - 第二篇

DFT – Fast algorithms and C implementations -Part2Radix-2DFT在我的上一篇博客里，已经介绍了基-2离散傅里叶变换的基本原理（参见第4(4)部分），总结一下就是：设N长序列的DFT为N长序列，若N被2整除，则：radix-2算法将一个N长的DFT分解为两个N/2长的DFT的线性组合。如果N/2仍然可以被2整除，则可以...
复制链接

扫一扫