基于DSP芯片的函数信号的FFT

信息与电气工程学院 课程设计任务书

2015 —2016 学年第 二 学期

专业： 通信工程 学号： 姓名： 课程设计名称： DSP课程设计

设计题目： 基于DSP芯片的函数信号的FFT 完成期限：自 2015 年 12 月 16 日至 2015 年 12 月 30 日共 2 周 设计依据、要求及主要内容（可另加附页）：

1、设计依据

1.1函数信号发生电路的构成，该电路的核心器件是ICL8038，产生高精度正弦波、三角波和方波。

2.2 FFT 的原理和参数生成公式。

2、设计要求

1． 熟悉FFT在实际中应用的效果； 2． 熟悉正弦、三角、方波信号的FFT变换；

3、主要内容：

利用函数信号发生器产生正弦波、三角波和方波。并且基于DSP芯片对其进行FFT变换。

1

指导教师（签字）： 批准日期： 年 月 日

2

摘 要

快速傅里叶变换（FFT）是将信号从时域变换到频域的一种方法，广泛应用于各种信号分析领域，文中介绍了FFT算法的基本原理和DSP中FFT算法的编程思想和设计原理及其硬件设计思想，基于TMS320C5502芯片用CCS仿真软件实现了FFT算法，文中以基2FFT为例，简要介绍了算法的实现，并画出了蝶形运行算图。然后编程实现算

法，输出波形。

关键字:DSP；CCS仿真软件；FFT

i

目 录

1.实验目的?????????????????4 2.实验设备?????????????????4 3.实验原理?????????????????4

3.1信号发生器原理 ? ??? ? ? ? ? ? ? ? ? ? ?? ? ? ?4 3.2 FFT 的原理和参数生成公式?????????????????5 4.实验步骤?????????????????7 5.运行并观察实验结果?????????????????8 6.实验结果?????????????????10 7.问题与思考?????????????????10 8.参考程序?????????????????11 9.总结?????????????????14 10.参考文献??????????????? ?? 14

ii

基于DSP芯片的函数信号的FFT

一、实验目的

1． 熟悉FFT在实际中应用的效果；

2． 熟悉正弦、三角、方波信号的FFT变换；

二、实验设备

１. 一台装有CCS软件的计算机； ２. DSP实验箱； ３. DSP硬件仿真器； ４. 网络交叉线。

三、实验原理

1.信号发生器原理

信号发生器模块可以产生20 Hz ~ 20 kHz 的 函数信号。信号波形有正弦波、三角波、方波；信号输出幅度 Vpp 在 0 ~ 3.3V 范围内可调；输出阻抗约600Ω左右。 信号发生器模块的主要函数信号发生电路及缓冲放大电路如下图所示。

首先简单介绍函数信号发生电路的构成，该电路的核心器件是ICL8038，它是单片多波形振荡集成电路，可以用很少的外围元件产生高精度正弦波、三角波和方波。其电路功能框图如下：

3

其工作波形如下图，

2.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。

4

我们称这样的RFFT 优化算法是包装算法：首先2N 点实数的连续输入称为“进包”。其次N 点的FFT 被连续运行。最后作为结果产生的N 点的合成输出是“打开”成为最初的与DFT 相符合的2N 点输入。

使用这一思想，我们可以划分FFT 的大小，它有一半花费在包装输入O（N）的操作和打开输出上。这样的RFFT 算法和一般的FFT 算法同样迅速，计算速度几乎都达到了两次DFT的连续输入。下列一部分将描述更多的在TMS320C55x 上算法和运行的细节。

2．FFT 的基本结构：

FFT信号流图如下：

X3(0)X1(0)X(0) x(0)N/4点X3(1)X1(1)DFTX(1) x(4)X4(0)X1(2)0 x(2)X(2)WN2N/4点X4(1)X1(3)DFT x(6)1X(3)WN2X2(0) x(1)0X(4)WNN/4点X2(1) 1DFTx(5)X(5)WN X2(2)20x(3)X(6)WNWN2N/4点 X2(3)DFT31WNx(7)X(7)WN2

整个过程共有log2N次，每次分组间隔为2^（L-1）----------------1=

第二次间隔为二，如第一个和第三个，第二个和第四个等

（2）基本运算单元以下面的蝶形运算为主：

计算公式如下：

rXm?1(p)?Xm?WNXm(q)Xm?1(p)?Xm?WXm(q)rN

（3）在FFT运算中，旋转因子WmN=cos(2πm/N)-jsin(2πm/N)，求正弦和余弦函数值的计算量是很大的。

W?erN?j2?RN?cos(2?R/N)?jsin(2?R/N)

（4）本程序采用的输入信号为：1024*sin(2*pi*3*t)，采样频率为1024

5

四、实验步骤

1． 连接好DSP开发系统，信号源模块的输出和AD50模块的输入，调节信号源使输入信号幅度小于2Vpp；

2． 测量AD50的5V电源，若电压太低，（少于5V比较多）更换LC01电感的值，或直接用导线取代LC01电感，或用0欧姆电阻也可以。 调节RPC03可变电阻，使Uc02运放的正输入端(3脚和5脚)，输入电平为2.5V把JC05跳到上面3.3V,AD50做Master；

3． 打开本实验工程文件（.\\singalsr\\signal.pjt），编译

，下载程序到DSP；

4． 在程序中这个位置上加上断点

：

5．选择view->graph->time/frequency? 。 设置对话框中的参数:

6

选择双时域，上端显示开始地址填入“input”，下端显示开始地址填入“output”，缓冲长度和显示长度输入NX的值，默认的NX是256。

6． 点击按F12运行程序形做实数FFT后的结果。

，观察上端波形和下端波形，验证是否下端波形为上端波

改变输入信号的频率，观察变化。 7．改变输入信号（三角或者方波），观察上下端波形的变化和相互关系。

8．改变tms320.h中NX的数值（64、128、256、512、1024）即FFT的点数，重复以上实验，观察变化。

五、运行并观察实验结果

⑴ 选择“Debug”菜单的“Animate”项，或按Alt+F5 键运行程序。 ⑵ 观察“Test Wave”窗口中时域图形；

7

图5

⑶ 在“Test Wave”窗口中点击右键，选择属性，更改图形显示为FFT。观察频域图形。

图6

⑷ 观察“FFT”窗口中的由CCS 计算出的正弦波的FFT。

图7

(5)改变输入函数 （方波）

INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*4+sin(PI*2*i/SAMPLENUMBER*8)))*1024;

8

图8

(5)改变输入函数（三角波）

INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*10)+sin(PI*2*i/SAMPLENUMBER*20)))*1024;

图9

六. 实验结果

通过观察频域和时域图，程序计算出了测试波形的功率谱，与CCS 计算的FFT 结果相近。

七.问题与思考

(1)观察图6和图7，可以看到二者波形相似，但横纵坐标均不相同，纵坐标大约是二倍的关系，横坐标大约为142倍。

9

(2)观察图8，因为两个频率比较相近，因此出现了前两个频谱交叠的现象。

八.参考程序

编程实现

int INPUT[SAMPLENUMBER],DATA[SAMPLENUMBER];

float fWaveR[SAMPLENUMBER],fWaveI[SAMPLENUMBER],w[SAMPLENUMBER]; float sin_tab[SAMPLENUMBER],cos_tab[SAMPLENUMBER]; /*说明：

分开计算cos(2*pi/N)及sin(2*pi/N)，合成蝶形运算的系数 */

void InitForFFT() {

int i;

for ( i=0;i

sin_tab[i]=sin(PI*2*i/SAMPLENUMBER); cos_tab[i]=cos(PI*2*i/SAMPLENUMBER); } }

/*说明：

输入信号，正弦函数，对1024*sin(2*pi*3*t的采样 */

void MakeWave() {

int i;

10

for ( i=0;i

INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024; } }

main() {

int i;

InitForFFT(); MakeWave();

for ( i=0;i

fWaveR[i]=INPUT[i]; fWaveI[i]=0.0f; w[i]=0.0f; }

FFT(fWaveR,fWaveI);

for ( i=0;i

DATA[i]=w[i]; }

while ( 1 ); // break point }

void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER]) {

int x0,x1,x2,x3,x4,x5,x6,xx; int i,j,k,b,p,L; float TR,TI,temp;

/********** following code invert sequence ************/ /*说明：

实现比特反转，改变输入信号的顺序，以方便使输出信号按自然顺序输出 */

for ( i=0;i

x0=x1=x2=x3=x4=x5=x6=0;

x0=i&0x01; x1=(i/2)&0x01; x2=(i/4)&0x01; x3=(i/8)&0x01;x4=(i/16)&0x01; x5=(i/32)&0x01; x6=(i/64)&0x01;

xx=x0*64+x1*32+x2*16+x3*8+x4*4+x5*2+x6;

11

dataI[xx]=dataR[i]; }

for ( i=0;i

dataR[i]=dataI[i]; dataI[i]=0; }

/************** following code FFT *******************/ /*说明：

基2fft算法，蝶形运算为核心 */

for ( L=1;L<=7;L++ ) { /* for(1) */ b=1; i=L-1; while ( i>0 ) {

b=b*2; i--; } /* b= 2^(L-1) */

for ( j=0;j<=b-1;j++ ) /* for (2) */ {

p=1; i=7-L;

while ( i>0 ) /* p=pow(2,7-L)*j; */ {

p=p*2; i--; }

p=p*j;

for ( k=j;k<128;k=k+2*b ) /* for (3) */ {

TR=dataR[k]; TI=dataI[k]; temp=dataR[k+b];

dataR[k]=dataR[k]+dataR[k+b]*cos_tab[p]+dataI[k+b]*sin_tab[p]; dataI[k]=dataI[k]-dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p]; dataR[k+b]=TR-dataR[k+b]*cos_tab[p]-dataI[k+b]*sin_tab[p]; dataI[k+b]=TI+temp*sin_tab[p]-dataI[k+b]*cos_tab[p]; } /* END for (3) */ } /* END for (2) */ } /* END for (1) */

//频域振幅平方和计算

for ( i=0;i

w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]); }

12

} /* END FFT */

九.总结

在毕业设计期间，有不少老师与同学给予我支持和帮助，在此表示衷心感谢！特别

感谢孙炎增老师对我的精心指导，孙老师有丰富多年的教学经验，给我们进行了深入而系统的指导。老师认真负责的工作态度，严谨的治学精神和深厚的理论水平都使我受益匪浅。由于我们经验不足和时间紧迫，我们在设计过程遇到了好多的疑难问题，这个过程中孙老师给我们提供了好多的指导和帮助，用他自己的丰富的设计经验帮我们解决了问题，同时我们也从他那里学习到了好多的知识，为今后我们走上社会工作岗位奠定了坚实的基础。还要感谢院领导和计算机房的老师给我们创造了良好的条件，以及其他帮助我的老师和同学。再次对他们表示衷心感谢！

随着毕业设计的即将结束，我将踏入社会，走向自己的工作岗位，拥有新的生活和

工作环境，大学四年的时间里我学会了扎实的专业知识，学会了如何独立生活和如何与人相处，在此，向母校的老师和帮助我的同学表示衷心感谢！

十.参考文献

【1】邹颜. DSP原理及应用.北京：电子工业出版社，2005.1

【2】纪震. 钟锦春， 强乐. DSP系统入门与实践.北京： 电子工业出版社2006.3 【3】汪安民. TMS320C54XX .北京： 清华大学出版社， 2002

【4】吴冬梅. 张玉洁. DSP技术及应用.北京： 北京大学出版社， 2006.8 【5】赵红怡. DSP技术与应用实例.北京： 电子工业出版社， 2003.6

【6】支长义. 程志平， 陈书立. DSP原理及开发应用 .北京：北京航空航天大学出版社，

2006.8

【7】谭浩强. C程序设计. 第二版.北京：清华大学出版社，1999.12

13

本文来源：https://www.bwwdw.com/article/2xtw.html