DSP实验报告精选.docVIP

一．课程设计目的 1．掌握用窗函数法设计FFT 快速傅里叶的原理和方法； 2．熟悉FFT 快速傅里叶特性； 3．了解各种窗函数对快速傅里叶特性的影响。 二.课程设计要求 将所给的128点FFT算法的源程序改成八点或16点程序，并对其进行硬件仿真。 三．实验设备 PC 兼容机一台，操作系统为Windows2000(或Windows98，WindowsXP，以下默认为Windows2000)，安装Code Composer Studio 2.0 软件。ICETEK–VC5509-A 硬件评估板。 四．实验原理 FFT 的原理和参数生成公式 公式（1）FFT 运算公式FFT 并不是一种新的变换，它是离散傅立叶变换（DFT）的一种快速算法。由于我们在计算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=2^L，L 为整数。将N=2^L 的序列 x(n)(n=0,1,……，N-1)，按N 的奇偶分成两组，也就是说我们将一个N 点的DFT 分解成两个N/2 点的DFT，他们又从新组合成一个如下式所表达的N 点DFT：一般来说，输入被假定为连续的。当输入为纯粹的实数的时候，我们就可以利用左右对称的特性更好的计算DFT。 我们称这样的RFFT 优化算法是包装算法：首先2N 点实数的连续输入称为“进包”。其次N点的FFT 被连续被运行。最后作为结果产生的N 点的合成输出是“打开”成为最初的与DFT 相符合的2N 点输入。使用这战略，我们可以划分FFT 的大小，它有一半花费在包装输入O（N）的操作和打开输出上。这样的RFFT 算法和一般的FFT 算法同样迅速，计算速度几乎都达到了两次DFT 的连续输入。下列一部分将描述更多的在TMS320C54x 上算法和运行的细节。 五.实验步骤 1.设置硬件仿真 设置CCS 通过ICETEK-5100USB 仿真器连接ICETEK-VC5509-A 硬件环境进行软件调试和开发 (1)单击桌面上图标：进入CCS 设置窗口。 (2)在出现的窗口中按标号顺序进行如下设置： (3)接着在下面的窗口中按标号顺序进行如下选择： (4)在出现的窗口按标号顺序进行如下设置： (5)在出现的窗口按标号顺序进行如下设置： 以上设置完成后，CCS 已经被设置成Emulator 的方式(用仿真器连接硬件板卡的方式)，并且指定通过ICETEK-5100USB 仿真器连接ICETEK-VC5509-A 评估板。如果您需要一直使用这一方式就不需要重新进行以上设置操作了。 2．打开工程 浏览程序，工程目录为C:\ICETEK-VC5509-EDULab\Lab0503-FFT 3．编译并下载程序 4．打开观察窗口 选择菜单View-Graph-Time/Frequency…进行如下图所示设置。 5．清除显示 在以上打开的窗口中单击鼠标右键，选择弹出式菜单中“Clear Display”功能。 6．设置断点 在程序FFT.c 中有注释“break point”的语句上设置软件断点。 7.运行并观察结果 ⑴选择“Debug”菜单的“Animate”项，或按F12 键运行程序。 ⑵观察“IIR”窗口中时域图形；观察滤波效果。 六.实验代码 原始代码： #include myapp.h #include ICETEK-VC5509-EDU.h #include scancode.h #include math.h #define PI 3.1415926 #define SAMPLENUMBER 128 void InitForFFT(); void MakeWave(); int INPUT[SAMPLENUMBER],DATA[SAMPLENUMBER]; float fWaveR[SAMPLENUMBER],fWaveI[SAMPLENUMBER],w[SAMPLENUMBER]; float sin_tab[SAMPLENUMBER],cos_tab[SAMPLENUMBER]; main() { int i; InitForFFT(); MakeWave(); for ( i=0;iSAMPLENUMBER;i++ ) { fWaveR[i]=INPUT[i]; fWaveI[i]=0.0