快速傅里叶变换FFT实验.pdfVIP

实验八、快速傅立叶变换（FFT）的实现 一、实验目的 在数字信号处理系统中，FFT 作为一个非常重要的工具经常使用，甚至成为 DSP 运算能力的一个考核因素。FFT 是一种高效实现离散付氏变换的算法。离散付氏变换 的目的是把信号由时域变换到频域，从而可以在频域分析处理信息，得到的结果再由 付氏逆变换到时域。 本实验的目的在于学习FFT 算法，及其在TMS320C54X 上的实现，并通过编程掌 握 C54X 的存储器管理、辅助寄存器的使用、位倒序寻址方式等技巧，同时练习使用 CCS 的探针和图形工具。另外在BIOS 子目录下是一个使用DSP/BIOS 工具实现FFT 的程序。通过该程序，你可以使用DSP/BIOS 提供的分析工具评估FFT 代码执行情况。 二、实验原理 ㈠ 基—2 按时间抽取FFT 算法 对于有限长离散数字信号{x[n]}，0 ≤ n ≤ N-1 ，其离散谱{x[k]}可以由离散付氏变 N −1 2π −j ( )nk ( ) N X k ∑x[n]e k 0,1,..., N −1 n 0 换（DFT ）求得。DFT 的定义为 可以方便的把它改写为如下形式： 不难看出，WN 是周期性的，且周期为N ，即 WN 的周期性是DFT 的关键性质之一。为了强调起见，常用表达式WN 取代W 以 N −1 ( ) nk X k ∑x[n]W k 0,1,..., N −1 N n 0 W (n+mN )(k +lN ) W nk m,l 0,±1,±2... N N 便明确其周期是N 。 由DFT 的定义可以看出，在x[n]为复数序列的情况下，完全直接运算N 点DFT 2 次复数乘法和N （N-1 ）次加法。因此，对于一些相当大的N 值（如1024） 需要（N-1 ） 来说，直接计算它的DFT 所作的计算量是很大的。FFT 的基本思想在于，将原有的N 150 点序列序列分成两个较短的序列，这些序列的DFT 可以很简单的组合起来得到原序列 的DFT 。例如，若N 为偶数，将原有的N 点序列分成两个（N/2 ）点序列，那么计算 2 2 N 点DFT 将只需要约[(N/2) ·2]=N /2 次复数乘法。即比直接计算少作一半乘法。因 2 表示直接计算（N/2 ）点DFT 所需要的乘法次数，而乘数2 代表必须完成 子（N/2 ） 两个DFT 。上述处理方法可以反复使用，即（N/2 ）点的DFT 计算也可以化成两个（N/4 ） 点的DFT （假定N/2 为偶数），从而又少作一半的乘法。这样一级一级的划分下去一直 到最后就划分成两点的FFT 运算的情况。比如，一个N = 8 点的FFT 运算按照这种方 法来计算FFT 可以用下面的流程图来表示： x(0) X(0) W0 x(1) X(1) W