基于DSP的FFT的实现

1绪论    1
1.1课题来源    1
1.2课题研究的目的意义    1
1.3国内外现状及水平    1
2开发运行环境CCS    2
3系统方案设计    3
3.1设计原理    3
4软件设计    6
4.1程序流程图    6
寻隐者不遇教学设计
4.2源程序    8
4.3设计步骤    8
总 结    13
参考文献    14
致  谢    15
附录    16

1 绪论
1.1课题来源
傅立叶变换是一种将信号从时域变换到频域的变换形式,是声学,语音,电信和信号处理等领域中一种重要的分析工具。快速傅立叶变换(FFT)是快速计算DFT的一种高效方法,FFT的出现使DFT的运算大大简化,运算时间缩短一至两个数量级之多,DSP芯片的出现使FFT的实现变得更加方便。
1.2课题研究的目的意义
随着电子技术和集成电路技术的飞速发展,数字信号处理已经广泛地应用于通信、信号处理、生物医学以及自动控制等领域中。离散傅立叶变换(DFT)及其快速算法FFT作为数字信号处理的基本变换,有着广泛的应用。特别是近年来,基于FFTODFM技术的兴起,进一步推动了对高速FFT处理器的研究。FFT算法从出现到现在已有四十多年代历史,算法理论已经趋于成熟,但是其具体实现方法却值得研究。面向高速、大容量数据流的FFT实时处理,可以通过数据并行处理或者采用多级流水线结构来实现。特别是流水线结构使得FFT处理器在进行不同点数的FFT计算时可以通过对模板级数的控制很容易的实现。分析和比较了各种FFT算法后,选择基2和基4混合频域抽取算法作为FFT处理器的而实现算法,一种高速、处理点数可变的流水线结构FFT处理器的实现方法。
1.3国内外现状及水平
数字信号处理(Digital Signal Processing,简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科。DSP有两种含义:Digital Signal Processing(数字信号处理)、Digital Signal Processor (数字信号处理器)。我们常说的DSP指的是数字信号处理器。数字
信号处理器是一种适合完成数字信号处理运算的处理器。20世纪60年代以来,随着计算机和信息技术的飞速发展,数字信号处理技术应运而生并得到迅速的发展。在过去的二十多年时间里,数字信号处理已经在通信等领域得到极为广泛的应用。
数字信号处理是利用计算机或专用处理设备,以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理,以得到符合人们需要的信号形式。
数字信号处理是以众多学科为理论基础的,它所涉及的范围及其广泛。例如,在数学领域,微积分、概率统计、随机过程、数值分析等都是数字信号处理的基本工具,与网络理论、信号与系统、控制论、通信理论、故障诊断等也密切相关。近来新兴的一些学科,如人工智能、模式识别、神经网络等,都与数字信号处理密不可分。可以说,数字信号处理是把许多经典的理论体系作为自己的理论基础,同时又使自己成为一系列新兴学科的理论基础。
DSP主要应用在数字信号处理中,目的是为了能够满足实时信号处理的要求,因此需要将数字信号处理中的常用运算执行得尽可能快,这就决定了DSP的特点和关键技术。适合数字信号处理的关键技术:DSP包含乘法器、累加器、特殊地址产生器、领开销循环的等;
提高处理速度的关键技术:流水线技术、并行处理技术、超常指令(VLIW)、超标量技术、DMA等。从广义上讲,DSP、微处理器和微控制器(单片机)等都属于处理器,可以说DSP是一种CPUDSP和一般的CPU又不同,最大的区别在于:CPU是冯.诺伊曼结构的;DSP是数据和地址空间分开的哈佛结构。
第2章 系统开发平台与环境
TI Code Composer Studio (CCStudio) TI eXpressDSPTM 实时软件技术的重要组成部分 , 它可以使开发人员充分应用 DSP 的强大功能。随着 TI TMS 320C 5000 C5K TMS 320C 6000 C6K DSP 平台的应用范围不断扩大 , 已经由其应用于下载视频流的手持因特网接入产品扩展到蜂窝通信网络和光网络的通信基础设施 ,eXpressDSPTM 也便获得了越来越多软件工程师的青睐。 
    eXpressDSP 还包含了 DSP/BIOS 可伸缩内核 TMS320TMDSP 标准算法的应用互操作性和可重复使用性以及 400 多家第三方厂商支持。大部分厂商提供 eXpressDSP 兼容算法、即插式应用以及种类繁多的硬件配件和咨询服务。
   平衡电桥Code Composer Studio 3.1 能够使开发人员编制出更多面向高级 DSP 应用的、紧凑的
高性能代码。通过实时接入的 DSP 开发者之家网站 , 内置 Update Advisor 对最新的工具、驱动程序及其技术进行自动的流线式管理。只要确保代码和功能调用的正确输入 , 凭借编辑器程序中的 Dynamic CodeMaestro 技术即可快速生成 C C++ 编码。
TI eXpressDSP产品市场营销经理 Mike Trujillo 说:通过充分利用CCS的工具与功能,编程人员能够大大缩短应用开发的时间。使用 CCStudio 生成的高度优化代码,工程师能够最大限度地发挥高性能 DSP 的全部功能,或者,在其它情况下能够以成本更低的器件来满足其应用需求。
    Code Composer Studio v3.1 使开发人员能够无缝管理任何复杂程度的项目dif , 其项目管理器通过一个集成版本的控制接口与通用资源控制器连接 , 管理着成千上万的文件。同时支持外部文件制作功能 , 使项目能够在 PC UNIX 平台上交叉运行。他们可以通过采用一个改进的产品开发流程 , 就可实现同一组项目文件的共享。于是可以使他们的开发周期缩短数周 , 并获得时间上提前于竞争对手推向市场的优势。
    对于那些希望把业界领先的 C6000 TM DSP 平台的高性能与 C5000 TM DSP 平台的低功耗相结合的系统开发者来说 ,Code Composer Studio v3.1 为使其同时调试混合多处理器成为了可能 Code Composer Studio v3.1 ioa还增加了实时数据交换 (RTDX TM ) 仿真功能 ,
可支持来自任何地方的 2 50 C5000 C6000 DSP 器件同时运行。此外 , 支持 RTDX 的仿真器还实现了实时 DSP/BIOS TM 仿真调试 , 该高级调试功能可以使开发人员更深入地了解 DSP 代码在硬件或仿真状态中的运行情况。
3系统方案设计
3.1设计原理
3.1.1离散傅里叶变换DFT
对于长度为N的有限长序列x(n),它的离散傅里叶变换(DFT)为
X(k)= N-nk        1
式中WN=e-j*2π/N 称为旋转因子或蝶形因子
DFT的定义可以看出,在x(n)为复数序列的情况下,对某个k值,直接按(1)式计算X(k) 只需要N次复数乘法和(N-1)次复数加法。因此,对所有Nk值,共需要N2次复数乘法
N(N-1)次复数加法。对于一些相当大有N值(如1024点)来说,直接计算它的DFT所需要的计算量是很大的,因此DFT运算的应用受到了很大的限制。
3.1.2快速傅里叶变换FFT
旋转因子WN 有如下的特性。
对称性:WNk+N/2=-WNk                2
周期性:WNn(N-k)=WNk(N-n)=WN-nk              3
利用这些特性,既可以使DFT中有些项合并,减少了乘法积项,又可以将长序列的DFT分解成几个短序列的DFTFFT就是利用了旋转因子的对称性和周期性来减少运算量的。
FFT的算法是将长序列的DFT分解成短序列的DFT。例如:N为偶数时,先将N点的DFT分解为两个N/2点的DFT,使复数乘法减少一半:再将每个N/2点的DFT分解成N/4点的DFT,使复数乘又减少一半,继续进行分解可以大大减少计算量。最小变换的点数称为基数,对于基数为2FFT算法,它的最小变换是2DFT
一般而言,FFT算法分为按时间抽取的FFTDIT FFT和按频率抽取的FFT家有学子DIF FFT两大类。DIF FFT算法是在时域内将每一级输入序列依次按奇/偶分成2个短序列进行计算。而DIF FFT算法是在频域内将每一级输入序列依次奇/偶分成2个短序列进行计算。两者的区别是旋转因子出现的位置不同,得算法是一样的。在DIF FFT算法中,旋转因子WN出现在输入端,而在DIF FFT算法中它出现在输入端。
假定序列x(n)的点数N2的幂,按照DIF FFT算法可将其分为偶序列和奇序列。
偶序列:x(2r)=x1(r)
奇序列:x(2r+1)=x2(r)
其中:r=0,1,2,…,N/2-1x(n)DFT表示为
               
式中,X1 (k)X2(k)分别为X1(r)X2(r)N/2DFT
由于对称性,WNk+N/2=-WNk因此,NDFT可分为两部分:
前半部分:x(k)=x1(k)+WkNx2(k)                       4
后半部分:x(N/2+k)=x1(k)-WkNx2(k)  k=0,1,…,N/2-1  5
从式4和式5可以看出,只要求出0~N/2-1区间x时间是等人的1(k)x2(k)的值,就可求出0~N-1区间x(k)N点值。
以同样的方式进行抽取,可以求得N/4点的DFT,重复抽取过程,就可以使N点的DFT用上组2点的DFT来计算,这样就可以大减少运算量。
2 DIF FFT的蝶形运算如图3.1所示。设蝶形输入为x1(k)x2(k),输出为x(k)x(N/2+K),则有
x(k)=x1(k)+WkNx2(k)                          6
x(N/2+k)=x1(k)-WkNx2(k)                      7
在基数为2FFT中,设N=2M,共有M级运算,每级有N/22FFT蝶形运算,因此,NFFT总共有MN/2个蝶形运算。
3.1 2 DIF FFT的蝶形运算
例如:基数为2FFT,当N=8时,共需要3级,12个基2 DIT FFT的蝶形运算。其信号流程如图3.2所示。
x(0)                                                                              x(0)

本文发布于:2024-09-20 17:52:26,感谢您对本站的认可!

本文链接:https://www.17tex.com/xueshu/720991.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:信号处理   数字   算法
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议