DSP的FFT实现设计报告

湖南科技大学信息与电气工程学院课程设计报告题 目： 基于DSP芯片的函数信号的FFT 专 业： 通信工程 班 级： 01班 姓 名： 学 号： 指导教师： 2015年 12 月 24 日信息与电气工程学院课程设计任务书20 20 学年第 学期专业： 通信工程 学号： 姓名： 课程设计名称： DSP原理及应用 设计题目： 基于DSP芯片的函数信号的FFT 完成期限：自 2015 年 12 月 21 日至 2015 年 12 月 24 日共 1 周设计依据、要求及主要内容（可另加附页）：DSP原理及应用课程设计是配合DSP原理及应用理论教学，为通信工程和电子信息工程专业开设的专业技术技能设计，课程设计对通信工程和电子信息工程专业的学生是一个非常重要的实践教学环节。 通过课程设计，使学生综合运用DSP技术课程和其他有关先修课程的理论和生产实际知识去分析和解决具体问题的能力得到提高，并使其所学知识得到进一步巩固、深化和发展。通过课程设计初步培养学生对工程设计的独立工作能力，学习设计的一般方法。通过课程设计树立正确的设计思想，提高学生分析问题、解决问题的能力。通过课程设计训练学生的设计基本技能，如计算、绘图、查阅设计资料和手册，熟悉标准与规范等。要求：（1）按时参加课程设计指导，定期汇报课程设计进展情况。（2）广泛收集相关技术资料。（3）独立思考，刻苦钻研，严禁抄袭。（4）按时完成课程设计任务，认真、正确地书写课程设计报告。（5）培养实事求是、严谨的工作态度和认真的工作作风。主要内容：（1）熟悉DSP处理器及其结构性能，掌握DSP芯片配套开发工具的使用方法。（2）按要求设计出硬件电路。（3）画出硬件连接原理图，并对硬件工作原理进行说明。（4）给出软件流程图及编写程序，每一条指令的后面附上相应的注释。（5）进行软、硬件调试，检查是否达到相关的功能。（6）写出调试方法。（7）设计报告结尾附上心得体会。指导教师（签字）： 批准日期： 年 月 日目 录摘要4一、实验目的4二、实验设备4三、实验内容 4四、实验原理4五、FFT算法的DSP实现过程5六、设计步骤7七、编译程序9八、实验过程10九、实验结果13十、问题与思考13十一、心得体会14十二、参考文献14摘要随着计算机和微电子技术的飞速发展，基于数字信号处理的频谱分析已经应用到各个领域并且发挥着重要作用。信号处理方法是当前机械设备故障诊断中重要的技术基础之一，分析结果的精确程度是诊断成功与否的关键因素。研究频谱分析是当前主要的发展方向之一。数字信号处理基本上从两个方面来解决信号的处理问题：一个是时域方法，即数字滤波；另一个是频域方法，即频谱分析.本文主要介绍了离散傅里叶变换以及快速傅里叶变换，通过对DFT以及FFT算法进行研究，从基础深入研究和学习，掌握FFT算法的关键。通过对DSP芯片工作原理以及开发环境的学习，掌握CCS的简单调试和软件仿真，在DSP芯片上实现对信号的实时频谱分析。关键字：DSP FFT算法 频谱分析一、实验目的1. 加深对DFT算法原理和基本性质的理解；2. 熟悉FFT的算法原理和FFT子程序的算法流程和应用；3. 学习DSP中FFT的设计和编程思想；4. 学习使用CCS的波形观察器观察波形和频谱情况；5. 熟悉FFT在实际中应用的效果；6. 熟悉正弦、三角、方波信号的FFT变换；二、实验设备.一台装有CCS软件的计算机；.DSP实验箱；.DSP硬件仿真器；.网络交叉线。三、实验内容 用DSP汇编语言及C语言进行编程，实现FFT运算、对输入信号进行频谱分析。四、实验原理信号发生器模块可以产生20 Hz 20 kHz 的 函数信号。信号波形有正弦波、三角波、方波；信号输出幅度 Vpp 在 0 3.3V 范围内可调；输出阻抗约600左右。信号发生器模块的主要函数信号发生电路及缓冲放大电路如下图所示。首先简单介绍函数信号发生电路的构成，该电路的核心器件是ICL8038，它是单片多波形振荡集成电路，可以用很少的外围元件产生高精度正弦波、三角波和方波。其电路功能框图如下：其工作波形如下图，五、FFT算法的DSP实现过程：DSP芯片的出现使FFT的实现方法变得更为方便。由于大多数DSP芯片都具有在单指令周期内完成乘法累加操作，并且提供了专门的FFT指令，使得FFT算法在DSP芯片实现的速度更快。FFT算法可以分为按时间抽取FFT和按频率抽取FFT两大类，输入也有实数和复数之分，一般情况下，都假定输入序列为复数。（一）FFT运算序列的存储分配FFT运算时间是衡量DSP芯片性能的一个重要指标，因此提高FFT的运算速度是非常重要的。在用DSP芯片实现FFT算法时，应允许利用DSP芯片所提供的各种软、硬件资源。如何利用DSP芯片的有限资源，合理地安排好所使用的存储空间是十分重要的。(二)FFT运算的实现用汇编程序实现FFT算法主要分为四步：1.实现输入数据的比特反转输入数据的比特反转实际上就是将输入数据进行码位倒置，以便在整个运算后的输出序列是一个自然序列。在用汇编指令进行码位倒置时，使用码位倒置可以大大提高程序执行速度和使用存储器的效率。在这种寻址方式下，AR0存放的整数N是FFT点的一半，一个辅助寄存器指向一个数据存放的单元。当使用位码倒置寻址将AR0加到辅助寄存器时，地址将以位码倒置的方式产生。2.实现N点复数FFTN点复数FFT算法的实现可分为三个功能块，即第一级蝶形运算、第二级蝶形运算、第三级至log2N级蝶形运算。对于任何一个2的整数幂，总可以通过M次分解最后成为2点的DFT计算。通过这样的M次分解，可构成M（即log2N）级迭代计算，每级由N/2个蝶形运算组成。3.功率谱的计算用FFT计算想x(n)的频谱，即计算X（k）=X(k)一般是由实部XR(k)和虚部Xl(k)组成的复数，即X（k）=XR(k)+jXl(k)因此，计算功率谱时只需将FFT变换好的数据，按照实部实部XR(k)和虚部Xl(k)求它们的平方和，然后对平方和进行开平方运算。但是考虑到编程的难度，对于求FFT变换后数据的最大值，不开平方也可以找到最大值，并对功率谱的结果没有影响，所以在实际的DSP编程中省去了开方运算。4.输出FFT结果(三)汇编语言程序程序主体由rfft-task、bit-rev、fft和power四个子程序组成。rfft-task：主调用子程序，用来调用其他子程序，实现统一的接口。bit-rev：位码倒置子程序，用来实现输入数据的比特反转。fft：FFT算法子程序，用来完成N点FFT运算。在运算过程中，为避免运算结果的溢出，对每个蝶形的运算结果右移一位。fft子程序分为三个功能块：第一级蝶形运算、第二级蝶形运算、第三级至至log2N级蝶形运算。(四）正弦系数表和余弦系数表：正弦系数表和余弦系数表可以由数据文件coeff.inc给出，主程序通过.copy汇编命令将正弦和余弦系数表与程序代码汇编在一起。在本例中，数据文件coeff.inc给出1024复数点FFT的正弦、余弦系数各512个。利用此系数表可完成81024点FFT的运算。（五）FFT算法的模拟信号输入：FFT算法的模拟信号输入可以采用C语言编程来生成一个文本文件sindata，然后在rfft-task汇编程序中，通过.copy汇编命令将生成的数据文件复制到数据存储器中，作为FFT算法的输入数据参与FFT运算。这种方法的优点是程序的可读性强，缺点是当输入数据修改后，必须重新编译、汇编和链接。六、设计步骤：1.启动CCS，在CCS中建立一个C源文件和一个命令文件，并将这两个文件添加到工程，再编译并装载程序：阅读Dsp原理及应用中fft 用dsp实现的有关程序。2.双击，启动CCS的仿真平台的配着选项。选择C5502 Simulator。3.启动ccs2后建立工程文件FFT.pjt4.建立源文件FFT.c与链接文件FFT.cmd5.将这两个文件加到FFT.pjt这个工程中。7. 创建out文件8. 加载out文件七、编译程序int INPUTSAMPLENUMBER,DATASAMPLENUMBER;float fWaveRSAMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER;float sin_tabSAMPLENUMBER,cos_tabSAMPLENUMBER;void InitForFFT()int i;for ( i=0;iSAMPLENUMBER;i+ )sin_tabi=sin(PI*2*i/SAMPLENUMBER);cos_tabi=cos(PI*2*i/SAMPLENUMBER);void MakeWave()int i;for ( i=0;iSAMPLENUMBER;i+ )INPUTi=sin(PI*2*i/SAMPLENUMBER*3)*1024;main()int i;InitForFFT();MakeWave();for ( i=0;iSAMPLENUMBER;i+ )fWaveRi=INPUTi;fWaveIi=0.0f;wi=0.0f;FFT(fWaveR,fWaveI);for ( i=0;iSAMPLENUMBER;i+ )DATAi=wi;while ( 1 );/ break pointvoid FFT(float dataRSAMPLENUMBER,float dataISAMPLENUMBER)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;iSAMPLENUMBER;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;dataIxx=dataRi;for ( i=0;iSAMPLENUMBER;i+ )dataRi=dataIi; dataIi=0; for ( L=1;L0 ) b=b*2; i-; /* b= 2(L-1) */for ( j=0;j0 ) /* p=pow(2,7-L)*j; */p=p*2; i-;p=p*j;for ( k=j;k128;k=k+2*b ) /* for (3) */TR=dataRk; TI=dataIk; temp=dataRk+b;dataRk=dataRk+dataRk+b*cos_tabp+dataIk+b*sin_tabp;dataIk=dataIk-dataRk+b*sin_tabp+dataIk+b*cos_tabp;dataRk+b=TR-dataRk+b*cos_tabp-dataIk+b*sin_tabp;dataIk+b=TI+temp*sin_tabp-dataIk+b*cos_tabp; /* END for (3) */ /* END for (2) */ /* END for (1) */for ( i=0;iGraph-Time/Frequency进行如下图所示设置。6设置断点：在程序FFT.c 中有注释“break point”的语句上设置软件断点。7运行并观察结果。 选择“Debug”菜单的“Animate”项，或按Alt+F5 键运行程序。 观察“Test Wave”窗口中时域图形； 在“Test Wave”窗口中点击右键，选择属性，更改图形显示为FFT。观察频域图形。 观察“FFT”窗口中的由CCS 计算出的正弦波的FFT。 (5) 改变输入信号（三角或者方波），观察上下端波形的变化和相互关系以及改变tms320.h中NX的数值（64、128、256、512、1024）即FFT的点数，重复以上实验，观察变化。INPUTi=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*4+sin(PI*2*i/SAMPLENUMBER*8)*1024; (5)改变输入函数INPUTi=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*10)+sin(PI*2*i/SAMPLENUMBER*20)*1024;九、实验结果通过观察频域和时域图，程序计