dsp实现快速傅里叶的C语言程序,DSP-快速傅立叶变换(FFT)算法实验

《DSP-快速傅立叶变换(FFT)算法实验》由会员分享,可在线阅读,更多相关《DSP-快速傅立叶变换(FFT)算法实验(10页珍藏版)》请在人人文库网上搜索。

1、中 南 大 学DSP技术实验报告 实验名称:快速傅立叶变换(FFT)算法实验专业班级: 信息0602学生姓名: 张倩曦 (学号:24)指导老师: 陈宁完成日期: 2009年12月2日中南大学信息科学与工程学院快速傅立叶变换(FFT)算法实验一实验目的1掌握用窗函数法设计FFT 快速傅里叶的原理和方法;2熟悉FFT 快速傅里叶特性;3了解各种窗函数对快速傅里叶特性的影响。二实验设备PC 兼容机一台,操作系统为Windows2000(或Windows98,WindowsXP,以下默认为Windows2000),安装Code Composer Studio 2.0 软件。三实验原理1FFT 的原理和。

2、参数生成公式:公式(1)FFT 运算公式FFT 并不是一种新的变换,它是离散傅立叶变换(DFT)的一种快速算法。由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法;一次复数加法则需二次实数加法。每运算一个X(k)需要4N 次复数乘法及2N+2(N-1)=2(2N-1)次实数加法。所以整个DFT运算总共需要4N2 次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法。如此一来,计算时乘法次数和加法次数都是和N2 成正比的,当N 很大时,运算量是可观的,因而需要改进对DFT 的算法减少运算速度。根据傅立叶变换的对称性和周期性,我们可以将DFT 运算中有些项合并。我们先设序列长。

3、度为N=2L,L 为整数。将N=2L 的序列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)的操作和打开输出。

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 ( i=0;i0 ) 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;。

7、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 。

8、for (1) */for ( i=0;iSAMPLENUMBER/2;i+ ) wi=sqrt(dataRi*dataRi+dataIi*dataIi); /* END FFT */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;(。

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;数组mo:FFT变换输出数组,长度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算法的原理的理解,并且逐步开始应用于实际。而在个别程序和分析方面存在的不足还需要在今后的学习中慢慢进步。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/566585.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

windows rt c语言,有arm架构的windows平板吗??与x86的win有什么不同?arm写c语言怎样?

有的。现在微软已推出windows8,10月即将上市。这个系统是跨平台的,其中就包括了ARM架构的平板微软也推出了自己的surface平板电脑ARM处理器的三大特点是:耗电少功能强、16位/32位双指令集和合作伙伴众多。1、体积小、低功耗、低成本、高性能;…

c语言某小卖部6天卖出热茶,习题1—1第二课时教案

2.4《线性回归方程》教学设计教学目标:(1)收集现实问题中两个有关联变量的数据作散点图,利用散点图直观认识变量间的相关关系;(2)在两个变量具有线性相关关系时,在散点较长中作出线性直线,用线性回归方程进行预测;(3)理…

c语言开发unity,[Unity3D]U3D开发项目总结

从2月份到现在第一个U3D项目也基本收工,虽然项目结局不是太好,但总算也是成功卖掉并上线,总结将近10个月的时间大家从端游转到手游或从COCOS2转到U3D的整个开发过程。1.资源无疑这是整个项目我觉得做的最差的地方,也是前期最为忽略…

linux dd 进度条,Progress 进度条 – DDProgressHUD

DDProgressHUDProgress 进度条,UIActivityIndicatorView 小菊花,弹窗,状态显示,高度自定义DDProgressHUD的介绍提供了四种类型的展示:显示无限旋转的加载图(比如小菊花,可以自定义),显示文字信息…

android notification自动消失,Notification点击事件和点击消失

版权声明:本文为博主原创文章,未经博主允许不得转载。https://blog.csdn.net/u012691505/article/details/53004930用到Notification遇到一些问题记录一下1.Notification的点击事件(点击弹出对话框)网上有很多做法,我的做法是:(1)…

Android动态图标包制作教程,安卓手机ico图标制作美化图文教程

如何让手机更加与众不同?今天我们就来学习如何利用出色的ico图标制作软件——Axialis IconWorkshop制作出美化手机屏幕的个性图标!如今,每个人都有属于自己的手机,每天使用手机聊QQ、刷微博、玩游戏,可以说手机已经是很多人生活中不可缺少的…

android progressbar 水平动画,Android ProgressBar 自定义样式(三),动画模式

果:和之前的一样,在布局文件中:android:id"id/progressBar3"android:layout_width"wrap_content"android:layout_height"wrap_content"android:indeterminate"false"android:indeterminateDrawable&…

android工程师绩效,Deloitteandroid开发工程师工资待遇怎么样 - 德勤华永会计师事务所(特殊普通合伙) - 职友集...

职责描述:• Ability to design holistic and data-driven financial risk management solutions that strengthen client’s financial risk management capabilities, while enhancing shareholder value• Ability to build an in-depth knowledge of the client…

android gridview行分割线,Android使用GridView实现表格分割线效果

使用gridview实现表格分割线效果,网格布局表格布局也是可以实现的。效果如下:1.主函数代码:package com.example.qd.douyinwu;import android.app.activity;import android.content.context;import android.os.bundle;import android.support…

android代码打开数据库,android – 如何正确关闭并重新打开Room数据库

您好我有2个应用程序依赖于制作和恢复应用程序数据库的备份,只需将数据库文件复制到SD卡中,然后很难确定如何在关闭它之后重新打开Room Database单例以创建数据库’副本.构建数据库:Database(version 15, exportSchema true, entities [list of entities])abstra…

android 按钮控制线程,关于手机安卓-线程测试

如题,我测试时,只执行了一次,第二次卡住时间到后,没有调用超时判断,特请教一下,代码如下:/**********线程控制*******Dim 检测线程Thread.SetShareVar "是否卡主", false检测线程 Th…

android hook 实例,代码实例分析android中inline hook

以下内容通过1、实现目标注入程序,2、实现主程序,3、实现注入函数,4、thumb指令集实现等4个方面详细分析了android中inline hook的用法,以下是全部内容:最近终于沉下心来对着书把hook跟注入方面的代码敲了一遍&#xf…

dax 筛选 包含某个字_筛选状态(ALL与REMOVEFILTERS)

这一章比较绕,但是帮助我们理清切片器的筛选状态以及主要函数对于公式筛选上下文的改变。1.切片器的筛选状态切片器的全选与全不选的状态的区别,我们可以看下图:全不选全选多个点击筛选逐个点击至全部结论初始状态,没有选择任何元…

android数据共享 设计,水文数据共享平台移动端的设计与实现

摘要:水文观测数据是进行流域水文科学及相关科学研究的基础.2010年国家自然科学基金委员会启动重大研究计划"黑河流域生态—水文过程集成研究",提出流域数据工作是"黑河计划"的重要组成部分.本文以黑河流域为研究背景,通过对黑河流域已有各种观…

山东鲁能轨道智能巡检机器人_温湿度传感器在轨道巡检机器人中的应用

▲地下综合管廊智慧管廊建设,是智慧城市在地下的一个缩影,有助于缓解“大城市病”,实现精细化和动态管理。与此同时,各种传感器技术也将被运用到地下综合管廊运维的每一环,这其中,地下管廊巡检机器人便是其…

android显示网络图片控件,Android控件之ImageView(二)

前言在上一篇文章中,我们讲解了如何加载本地图片,那么在实际项目中 ImageView 大多数使用场景是加载网络图片,网络图片其实就是存储在服务器上的文件,我们需要从服务器获取到文件的二进制输入流 Inpustream ,然后将其转…

坯子库曲面推拉教程_一招曲面流动,搞定99%异形建模

曲面流动可以建什么模型?这样的▼这样的▼还有这样的▼那究竟如何使用曲面流动呢?本文告诉你!曲面流动是什么?曲面流动功能相当于Rhino(犀牛)中的“沿曲面流动”命令,可以使来源几何体群组或组件,根据基准平面为参照,…

android o 结构光流程,惊艳亮相!一分钟看懂OPPO Find X 3D结构光技术是什么鬼,太牛了...

法国时间6月19日,OPPO在巴黎卢浮宫正式举办未来旗舰Find X发布会。此次亮相的Find X新机,既有充满艺术感的3D玻璃机身设计,又有3D结构光、曲面全景屏、双轨潜望结构等多项黑科技。众多黑科技中,以3D结构光O-face最受数目。据悉&am…

jframe运行和预览大小不一样_同一款车型,为什么配的轮胎大小还会不一样?【内含福利】...

小编又来给大家送福利啦~阅读完文章之后 戳文章底部阅读原文 岁末巨献 国货好胎!免费体验价值1999元国产轮胎!相信有不少朋友,买车的时候一定会有类似的经历。选车、试驾的时候,看到试驾车配置相当漂亮,尤其搭配的轮圈…

关于鸿蒙系统传统文化的作文,关于生活中的传统文化作文(精选10篇)

关于生活中的传统文化作文(精选10篇)在学习、工作或生活中,大家都接触过作文吧,作文可分为小学作文、中学作文、大学作文(论文)。一篇什么样的作文才能称之为优秀作文呢?以下是小编帮大家整理的生活中的传统文化作文(精选10篇),仅…