基于DSPLib的FFT程序设计.docVIP

基于DSPLib的FFT程序设计实验 一、 实验目的 1、了解FFT的原理； 2、了解在DSP中FFT 的设计及编程方法； 3、熟悉对DSPLIB的调用方法； 二 、实验内容(要求) 信号x(n)是2个频率分别为1kHz和2kHz余弦信号的合成，取样频率为fs=16000Hz,采用FFT对该信号进行频域分析。 三 、实验原理（背景知识） 时间抽选基2FFT算法的基本原理。 对于序列 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运算流图 DSPLib的应用。CCS提供DSPLib，其中包含了基本的数学计算和数字信 号处理常用算法模块的函数，可直接调用这些函数实现一些数字信号处理算法。本实验主要使用了DSPLib的cbrev，rfft两个函数实现数字序列的逆序和DFT的计算。dsplib位于C:\ti\c5400\dsplib 其用户手册为:TMS320C54x dsp library programers reference.pdf 3、CCS中的graph详细使用说明。 CCS提供显示时频图（Time/Frequency），运算结果也可以通过CCS提供的图形功能经过一定处理显示出来，CCS提供的图形显示包括时频分析、星座图、眼图和图像显示。如表1所示，用户准备好需要显示的数据后，选择命令View→Graph，设置相应的参数，即可按所选图形类型显示数据。 程序框图如下： 四 、实验步骤 1、打开CCS，并设置好相对应的参数，采用软件仿真。 2、编写主函数，包含必要的头文件，并设置好相对应的路径。完毕将支持的库函数加入工程中。 3、接着编写中断向量表，中断向量表只要和主函数main中断一一对应就行。 4、编译，看一下会产生哪些链接文件，然后编写链接文件。在这里注意字对齐的格式要求。 5、编译通过后，用软件仿真，选择命令View→Graph→Time/Frequency对话框，验证是否有逻辑错误，然后进行修改。 五、程序源代码 1、主函数： #include math.h #include tms320.h #include dsplib.h #define Nx 1024 #define pi 3.1415924 #define PMST (unsigned int*)0x1D int i,j; int scale=1; int noscale=0; DATA x[Nx],x_tempt[Nx]; DATA y[Nx/2]; #pragma DATA_SECTION(x,.input) #pragma DATA_SECTION(x_tempt,.input_tempt) #pragma DATA_SECTION(y,.output) void dataio(); void main() { *PMST=0x00a0; dataio(); for(i=0;iNx;i++) { x_tempt[i]=x[i]; } cbrev(x,x,Nx/2); rfft(x,Nx,scale); y[0]=x[0]; y[Nx/2-1]=abs(x[1]); j=1; for(i=2;iNx;i=i+2) { y[j]=sqrt(pow(x[i],2)+pow(x[i+1],2)); j++; } return; } void dataio() { //get data for io } 中断函数： .ref _c_int00 .sect .vectors rs: BD _c_int00 ;Rest nop nop nmi: .space 4*16 ;NMI, SINT16 sint17:.space 4*16 ;SINT17 sint18: .space 4*16 ;SINT18 sint19: .space 4*16 ;SINT19 sint20: .space 4*16 ;SINT20 sint21: .space 4*16 ;SINT21 sint22: .sp