实验七 快速傅立叶变换(FFT)算法的DSP实现(C 语言).docVIP

实验七 快速傅立叶变换（FFT）算法的DSP实现（C 语言） 一、实验目的 1.掌握FFT算法的基本思想。 2.掌握利用 CCS软件中的dsplib库进行fft算法的程序设计 二、实验环境 1.奔腾IV计算机 2.Code Composer Studio (CCS)软件 三、实验原理 1.FFT算法的基本思想 对于序列 x[n](0≤n≤N-1)，其频谱为: X(k)=DFT[x(n)]=(0≤n≤N-1)， 其中： 在x[n]为复数序列的情况下，完全可以直接运算N 点DFT 需要(N-1)2次复数乘法和N（N-1）次复数加法。因此，对于一些相当大的N 值来说，直接计算它的DFT 所需计算量很大。 FFT（快速傅立叶变换）的基本思想为：将原来的N点序列最终分成分成两点为一组序列，并将这些序列的DFT 通过蝶形运算（见下图)组合起来得到原序列的DFT。 N点FFT 仅需 次复数乘法和次复数加法。 图1 8点FFT运算流图 2. ccs软件中dsplib（信号处理库）简介 TMS320C54X系列函数库（DSPLIB）是对C语言编程可调用优化的DSP函数库，它含有50个通用目的的信号处理程序，全部由汇编语言编写，并可由C语言调用，方便C语言与汇编语言混合编程。这些程序用在计算强度大、执行速度重要的实时运算中。通过使用这些程序，可以取得较用C语言编写的相关程序快的多的运行速度，另外通过使用现成的程序可以使开发速度大大加快。DSPLIB可进行的运算有：FFT运算、滤波与卷积运算、自适应滤波运算、相关运算、数学函数运算、三角函数运算、矩阵运算等。 dsplib位于C:\ti\c5400\dsplib 其用户手册为:TMS320C54x dsp library programers reference.pdf 三、实验内容 1.信号x(n)是2个频率分别为1kHz和2kHz余弦信号的合成，取样频率为fs=16000Hz,采用FFT对该信号进行频域分析。 实验步骤： (1)将 fft.c,fft.cmd加入工程fft.pjt (2)将tms320.h,dsplib.h,c54xdsp.lib,rts.lib加入工程fft.pjt, 如下图： 其中： c54xdsp.lib在c:\ti\c5400\dsplib目录下； dsplib.h、tms320.h在c:\ti\c5400\dsplib\include目录下； rts.lib在出 c:\ti\c5400\cgtools\lib目录下。 (3)运行该程序，并进行以下设置： (4)画出结果波形： (5)观察波形中出现的峰值，峰值对应的横坐标分别为 9 、 17 。 它们的含义是： 频率点 ，与实际频率fk的关系是： fk=fs/N*k 。 2.参照上述程序： 产生一个双音多频信号”2”,该信号x(n)的低频和高频分别为697Hz和1336Hz,取样频率为fs=8000Hz,并采用FFT对其进行频谱分析,写出程序并画出频谱（链接命令文件与上述程序相同）。 程序：fftDTMF1.c #include math.h #include tms320.h #include dsplib.h #define NX 64 #define pi 3.1415926 #pragma DATA_SECTION (x,.input) int i; int scale = 1; int noscale = 0; DATA x[NX]; main() { for (i=0; iNX; i++) { x[i]=(int)0x2000*(sin(2*pi*i*697/8000)+sin(2*pi*i*1336/8000)); } // compute cbrev(x,x,NX/2); rfft(x,NX,scale); return; }