62016701核心系统的FFT运算程序的设计与分析.pptVIP

利用6201/6701核心系统的FFT运算程序的设计与分析 1.前言 2.FFT运算的原理 3.FFT运算程序的设计 4.程序的运行分析 前言 快速傅里叶变换是一种高效实现离散傅里叶变换(DFT)的算法。FFT在Fourier变换中起着极其重要的作用。由于FFT的出现，是科学分析的许多方面有了完全的改观。在数字信号处理系统中,FFT作为一个非常重要的工具经常被使用，甚至成为DSP运算能力的一个考核因素。 FFT运算的原理 FFT 并不是一种新的变换，它是离散傅立叶变换（DFT）的一种快速算法。 FFT 运算公式： 由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法；一次复数加法则需二次实数加法。每运算一个X（k）需要4N 次复数乘法及2N+2（N-1）=2（2N-1）次实数加法。所以整个DFT 运算总共需要4N^2 次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法。如此一来，计算时乘法次数和加法次数都是和N^2 成正比的，当N 很大时，运算量是可观的，因而需要改进对DFT 的算法减少运算速度。 根据傅立叶变换的对称性和周期性，我们可以将DFT 运算中有些项合并。 通过上述方法N点DFT分解为两个N/2点的DFT N点DFT的一次时域抽取分解图 FFT运算程序的设计 程序流程图 程序的运行分析 设置断点 将光标移至源程序中ChangeOrder(x1r,x1i); /*调用倒序子程序进行排序*/ 的位置，然后点击按钮 。按F5全速运行程序，在断点处，程序将自动停止。 添加watch window观察变量 点击View-Watch Window，CCS界面右下方会弹出一个Watch Window窗口，选中窗口的Watch1 栏。 在Name栏单击表达式按钮 ，并将x1r作为要观察的变量名输入。然后将x1r栏展开，下面将会出现x1r数组中的所有数据。 单步运行程序 单步运行程序，同时观察watch window中x1r数组的排列顺序，排序结果与正确顺序相符。（如果不符，则可以及时找出修改，再编译运行，以得到正确的顺序。） 排序前的x1r ，下图为排序后的x1r 添加view graph比较输入输出波形 选择View—Graph—Time/Frequency，会弹出一个图像属性窗口，将Display设为FFT Magnitude，Graph Title设为x，Start Address设为x，Acquisition Buffer Size设为1024，FFT Framesize设为1024，其余默认(要设置的部分为下图中蓝色部分）。可得512点的FFT运算的输入数据的频谱图。同样可以得到512点FFT运算后的模的时域图。 * 512点FFT运算频谱图与实际运算后模的时域图 * * *