快速傅里叶变换C语言程序.docxVIP

快速傅里叶变换C语言程序

代码解析

2.位反转函数：`bit_reverse`函数负责将输入序列按照位反转的规则重新排序。这是DIT-FFT算法的要求，确保后续的蝶形运算能够正确进行。其核心思想是通过寻找每个序号的“位反转序号”并进行交换。

3.FFT主函数：`fft`函数是实现的核心。

*首先调用`bit_reverse`对输入序列进行重排。

*然后通过三重循环实现各级的蝶形运算：

*外层循环（stage）：对应FFT的级数，从1到m（m=log2(n)）。每一级处理不同大小的蝶形单元。

*中层循环（group）：遍历当前级中的所有蝶形运算组。每组包含`LE`个点（`LE=2^stage`）。

*内层循环（butterfly）：遍历一个蝶形组内的所有蝶形单元。每个蝶形单元包含两个输入点，它们之间的距离为`LE_half=LE/2`。

*旋转因子（twiddlefactor）：在每个蝶形单元中，需要计算并应用旋转因子`W`。基本旋转因子`W_m`根据当前级的`LE`计算得出，然后通过不断自乘来获得该组内所有蝶形单元所需的旋转因子。

4.主函数测试：主函数中提供了一个简单的测试用例。我们对一个8点的实序列`[1,1,1,1,0,0,0,0]`进行FFT。根据理论分析，该序列的FFT结果在0Hz、fs/8、2fs/8、3fs