调用DSP库函数实现FFT的运算(2)

#define SIGNAL1F 20 #define SAMPLEF 64 #define PI 3.1415926 #define SAMPLENUMBER 128

short INPUT[SAMPLENUMBER],x[SAMPLENUMBER]; float OUTPUT[SAMPLENUMBER]; void MakeWave();

void MakeWave() {

int i;

float sample_step=1.0/SAMPLEF; float j=0.0;

for ( i=0;i<SAMPLENUMBER/2;i++ ) {

INPUT[i]=sin(PI*2*j*SIGNAL1F)*1024; j=j+sample_step; } }

void main() { int i;

MakeWave();

for(i=0;i<SAMPLENUMBER;i=i+2) {

x[i]=INPUT[i/2]; }

for(i=1;i<SAMPLENUMBER;i=i+2)

x[i]=0.0; }

cbrev(x,x,SAMPLENUMBER/2);

cfft(x,SAMPLENUMBER/2, SCALE); //unpacki(x,SAMPLENUMBER/2);

//cbrev(x,x,SAMPLENUMBER/2);

//cifft(x,SAMPLENUMBER/2,SCALE);

for ( i=0;i<SAMPLENUMBER;i++ ) { OUTPUT[i]=x[i];}

while ( 1 ); // break point }

五、实验步骤

1.实验准备 设置软件仿真模式

2.启动CCS，打开工程，浏览程序

3.编译程序

4.导入.out文件并运行

5.分别设置窗口，并出图

六、实验结果

1、正弦输入波形（时域）

分析：由于采样频率为64HZ，相对于正弦函数频率他的采样频率较小，所以产生的时域图片会有失真。

输入波形（频域）

分析：有图可得峰值所在点为20π符合题目要求

经过FFT处理后的波

分析：出现峰值和谷值的地方为20和44且相对称，出现负值的原因是没有取模。另外在出图时设置grath的采样频率为2HZ，因为x[i]中有实数和虚数。

三个图放在一起比较

六、调试过程中遇到的问题和解决办法

1. 在程序运行的过程中会出现dsplib.h和tms320.h文件不存在的现象，

这是由于DSPLIB安装在固定的子目录上，而应用程序在桌面上，程序编译连接时，找不到DSPLIB中相应的程序。所以需要对工程的bulid options选项中的两个地方经行设置。首先是找不到dsplib.h、tms320.h文件时，在compiler标签下选中preprocessor选项，在Include Search Path栏中填入dsplib.h、tms320.h所在子目录，此处为下图所示：

当出现连接时找不到FFT运算相应的汇编程序，此时可在linker标签下选中basic

选项，在Library Search Path栏中填入55xdsp.lib库文件所在路径，本次试验示例如下：

此时完成设置，编译连接时可将应用程序及DSPLIB中的程序连接起来，程序可以

正常编译。

2． 出图的时候，通过view栏进行了出图窗口编辑，但是仍然出不了图，此时需要通过

debug栏中的Animate使程序运行出图。

七、心得体会

通过本次实验我理解了FFT快速傅里叶变化算法和CFFT算法，同时对

DSP课程我也有了更深入的了解，这一段时间的程序设计经历增加了我自主学习的能力。

虽然有参考案例，但是对于第一次接触CCS软件的我来说并不简单，

首先对于程序我需要看懂案例中的内容，然后在相应的地方进行改动，以符合我的实验案例。经过不断的摸索，现在终于能熟练的掌握CCS仿真软件的使用，同时也能就简单的dap问题进行编程运行模拟。

八、参考程序

参考程序：

G:\0f8d8f144f82e50a01dfc89984afbe02\c5500\dsplib\EXAMPLES\CFFT

设计成绩评定