榆社网站建设,qq手机版在线登录入口,莱芜金点子信息港厂房出租,秦皇岛属于河北哪个市《DSP-快速傅立叶变换(FFT)算法实验》由会员分享#xff0c;可在线阅读#xff0c;更多相关《DSP-快速傅立叶变换(FFT)算法实验(10页珍藏版)》请在人人文库网上搜索。1、中 南 大 学DSP技术实验报告 实验名称#xff1a;快速傅立叶变换(FFT)算法实验专业班级#xff1a; 信…《DSP-快速傅立叶变换(FFT)算法实验》由会员分享可在线阅读更多相关《DSP-快速傅立叶变换(FFT)算法实验(10页珍藏版)》请在人人文库网上搜索。1、中 南 大 学DSP技术实验报告 实验名称快速傅立叶变换(FFT)算法实验专业班级 信息0602学生姓名 张倩曦 (学号24)指导老师 陈宁完成日期 2009年12月2日中南大学信息科学与工程学院快速傅立叶变换(FFT)算法实验一实验目的1掌握用窗函数法设计FFT 快速傅里叶的原理和方法2熟悉FFT 快速傅里叶特性3了解各种窗函数对快速傅里叶特性的影响。二实验设备PC 兼容机一台操作系统为Windows2000(或Windows98WindowsXP以下默认为Windows2000)安装Code Composer Studio 2.0 软件。三实验原理1FFT 的原理和。2、参数生成公式公式(1)FFT 运算公式FFT 并不是一种新的变换它是离散傅立叶变换(DFT)的一种快速算法。由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法一次复数加法则需二次实数加法。每运算一个X(k)需要4N 次复数乘法及2N2(N-1)2(2N-1)次实数加法。所以整个DFT运算总共需要4N2 次实数乘法和N*2(2N-1)2N(2N-1)次实数加法。如此一来计算时乘法次数和加法次数都是和N2 成正比的当N 很大时运算量是可观的因而需要改进对DFT 的算法减少运算速度。根据傅立叶变换的对称性和周期性我们可以将DFT 运算中有些项合并。我们先设序列长。3、度为N2LL 为整数。将N2L 的序列x(n)(n0,1,N-1)按N的奇偶分成两组也就是说我们将一个N 点的DFT 分解成两个N/2 点的DFT他们又重新组合成一个如下式所表达的N 点DFT一般来说输入被假定为连续的。当输入为纯粹的实数的时候我们就可以利用左右对称的特性更好的计算DFT。我们称这样的RFFT 优化算法是包装算法首先2N 点实数的连续输入称为“进包”。其次N 点的FFT 被连续运行。最后作为结果产生的N 点的合成输出是“打开”成为最初的与DFT 相符合的2N 点输入。使用这一思想我们可以划分FFT 的大小它有一半花费在包装输入O(N)的操作和打开输出。4、上。这样的RFFT 算法和一般的FFT 算法同样迅速计算速度几乎都达到了两次DFT的连续输入。下列一部分将描述更多的在TMS320C55x 上算法和运行的细节。5程序流程图四实验步骤(一)第一部分1实验准备-设置软件仿真模式。-启动CCS。2打开工程。浏览程序工程目录为C:ICETEK-VC5509-EDULabLab0503-FFTFFT.pjt。3编译并下载程序。4打开观察窗口选择菜单View-Graph-Time/Frequency5清除显示在以上打开的窗口中单击鼠标右键选择弹出式菜单中“Clear Display”功能。6设置断点在程序FFT.c 中有注释“break poi。5、nt”的语句上设置软件断点。7运行并观察结果。选择“Debug”菜单的“Animate”项或按F12 键运行程序。观察“Test Wave”窗口中时域图形在“Test Wave”窗口中点击右键选择属性更改图形显示为FFT。观察频域图形。观察“FFT”窗口中的由CCS 计算出的正弦波的FFT。8退出CCS。9. 实验结果10.源代码(C语言)#include myapp.h#include ICETEK-VC5509-EDU.h#include scancode.h#include #define PI 3.#define SAMPLENUMBER 128void InitForFFT(。6、);void MakeWave();int INPUTSAMPLENUMBER,DATASAMPLENUMBER;float fWaveRSAMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER;float sin_tabSAMPLENUMBER,cos_tabSAMPLENUMBER;main()int i;InitForFFT();MakeWave();for ( i0;i0 ) bb*2; i-; /* b 2(L-1) */for ( j0;j0 ) /* ppow(2,7-L)*j; */pp*2; i-;pp*j;for ( kj;。7、k128;kk2*b ) /* for (3) */TRdataRk; TIdataIk; tempdataRkb;dataRkdataRkdataRkb*cos_tabpdataIkb*sin_tabp;dataIkdataIk-dataRkb*sin_tabpdataIkb*cos_tabp;dataRkbTR-dataRkb*cos_tabp-dataIkb*sin_tabp;dataIkbTItemp*sin_tabp-dataIkb*cos_tabp; /* END for (3) */ /* END for (2) */ /* END 。8、for (1) */for ( i0;iSAMPLENUMBER/2;i ) wisqrt(dataRi*dataRidataIi*dataIi); /* END FFT */void InitForFFT()int i;for ( i0;iSAMPLENUMBER;i )sin_tabisin(PI*2*i/SAMPLENUMBER);cos_tabicos(PI*2*i/SAMPLENUMBER);void MakeWave()int i;for ( i0;iSAMPLENUMBER;i )INPUTisin(PI*2*i/SAMPLENUMBER*3)*1024;(。9、二)第二部分1.程序参数说明extern void InitC5402(void)extern void OpenMcBSP(void)extern void CloseMcBSP(void)extern void READAD50(void)extern void WRITEAD50(void)void kfft(pr,pi,n,k,fr,fi,l,il)基2快速傅立叶变换子程序n为变换点数应满足2的整数次幂k为幂次(正整数)数组x 输入信号数组数据存放于地址为3000H307FH存储 器中转为浮点型后生成x数组长度128数组moFFT变换输出数组长度128浮点型。10、整型后写入存储器中。2.子程序流程图3.启动CCS 2.0用Project/Open打开“ExpFFT01.pjt”工程文件双击“ExpFFT01.pjt” 及“Source”可查看各源程序加载“ExpFFT01.out”4.在主程序中k处设置断点5.单击“Run”运行程序或按F5运行程序程序将运行至断点处停止6.用View / Graph / Time/Frequency打开一个图形观察窗口7.单击“Animate”运行程序。或按F10运行调整观察窗口并观察变换结果。五.实验总结这次DSP技术的实验是在我们学习了数字信号处理、C语言程序设计的先修课程和DSP技术以及。11、熟悉了CCS软件之后完成的。实验一我们首先熟悉要用的各种软硬件环境接下来的实验二和三我们完成了有限冲击响应滤波器(FIR)算法实验和快速傅里叶变换(FFT)算法实验。这次实验我个人认为在以下几个方面收获最大1. 首先是实验前对自己过去几个学期已学理论知识的巩固加深、综合应用以及对实验相关资料的收集能力的提高。2. 实验中各个部分的要求要理解掌握认真对待。实验锻炼了我们认真的态度和严谨的精神。3. 当遇到问题和困难的时候一定要保持冷静慢慢检查自己多思考、多尝试才能获得最终的成功。4. 实验中遇到问题无法解决的时候要主动向同学请教同时也要乐于帮助同学。我认识到在实验的过程中要互相信任、互相帮助。通过实验我们加深了相互间的友情。总之这次实验让我得到了很多平时理论学习中不能获得的收获加深了我们对DSP特别是FIR滤波器和FFT算法的原理的理解并且逐步开始应用于实际。而在个别程序和分析方面存在的不足还需要在今后的学习中慢慢进步。
