- 1、本文档共51页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
基于位倒序寻址方式的FFT现 制作者 魏俊 基于位倒序寻址方式的FFT现 摘要:数字信号处理在许多领域中有重要的意义,数字信号处理目的是对数字信号进行处理和转换。它采用运算的方法来处理数字信号。由于傅里叶变换和z 变换在计算机上实现很不方便,所以引入了离散傅里叶变换,它在数字信号处理中有着重要的地位,面对大量的数据运算,为了提高运算速度,人们对离散傅里叶变换进行了改进,离散傅里叶运算效率大大提高,这种快速傅里叶算法的出现,在数字信号处理中发挥了巨大的作用。本文介绍了有关离散傅里叶变换和快速傅里叶变换性质、算法原理和应用举例。并用c语言编程实现快速傅里叶转换。 关键词:离散傅里叶转换、快速傅里叶转换、c语言编程 第一章 绪论 1.1数字信号处理的基本概念 1.2数字信号处理的特点 1.3数字信号处理的应用 1.1数字信号处理的基本概念 信号是信息的物理表现形式,信息是信号的具体内容。根据信息的载体不同,信号可以是电的、磁的、光的、声的、热的和机械的等不同种类的信号。 信号通常是一个或几个自变量的函数。如果只有一个自变量,则称为一维信号;如果是两个或以上自变量,则称为多维信号。 信号的自变量可以是时间、距离、电压、温度等不同形式。本书一般把信号看做时间的函数。 1.2数字信号处理的特点 1.灵活性高 2精度高 3易于大规模集成 4性能指标高 1.3数字信号处理的应用 1通用DSP:数字滤波、卷积、相关、FFT 2语音:语音通信、语言编码、识别、; 3图像图形:机器人视觉、图像传输和压缩、等; 4控制:磁盘控制器、机器人控制、激光打印机、电机控制; 5军事:雷达、保密通信、声纳、导航、传感器融合等; 6电讯:调制解调器、蜂窝电话、、视频会议等; 7汽车:自动驾驶控制、故障分析、导航、汽车音响等; 8消费:数字音响、数字电视、MP3播放器、数码相机等 第二章 离散傅里叶变换(DFT) 2.1. 傅里叶变换的几种可能形式及离散傅里叶变换定义 2.2DFT的性质 2.1傅里叶变换的几种可能的形式及离散傅里叶变换的定义 2连续时间·离散频率--傅里叶级数 3离散时间·连续频域-----序列的傅里叶变换 4离散时间·离散频率----离散傅里叶变换 由于 2.1.1DFS的性质 1线性 DFS = 其中 a,b为任意常数 所得到的频域序列也是周期序列,周期为N. 2序列的时域移位 4周期卷积 如果 , 则 5对偶性 2.2DFT的性质 2.2.1线性性质 如果和是两个有限长序列,长度分别 ,若 2.2.2圆周移位性质 2.2.3对偶性 2.3离散傅里叶的应用 2.3.1用DFT对信号进行谱分析 傅里叶变换理论 信号持续时间有限长,其频谱是无限宽。 信号的频谱有限长,在时域中,该信号的持续时间无限长。 上述两种情况,在时域或频域中进行采样,得到的序列都是无限长序列,不满足DFT的变换条件。 采用的处理方法:在频域中用滤波器滤除高于折叠频率的高频分量,在时域中则是截取有限点进行DFT。 第三章 快速傅里叶变换 3.1改进计算的方法 3.2按时间抽取的FFT算法 3.2.1算法原理 3.2.2按时间抽取的FFT算法的运算量与运算特点 3.3按频率抽取的FFT算法 3.3.1算法原理与运算特点 3.3.2按时间抽取与按频率抽取的同异 离散傅里叶的计算工作量 通常X(k), 都是复数,所以计算一个X(k)的值需要N次复数乘法运算,和N-1次复数加法运算,那么所有的X(k)就要NxN复数乘法运算,N(N-1)复数加法运算,当N很大时,计算量就相当惊人,如果当N=1024时,则要完成1048576次运算这样难易做到实时处理。 改进途径,利用 的周期性和对称性 3.2按时间抽取的FFT算法 3.2.1算法原理 一. 算法原理基于(2FFT) N/2点的DFT,先将x(n)按n的奇偶分成两组DFT,不足时补零 这样 就有: 由于 所以可表示为 X(k)的后一半的确定 蝶形运算 蝶形运算 前半部为X(0)~ X(3),后半部分为X(4)~X(7)整个过程如下图所示: 有图可知一个N点分解为两个N/2点DFT后,如果直接计算N/2点DFT,则每一个N/2点DFT只需要/4次复数乘法N/2(N/2-1)次复数加法。两个N/2点DFT共需/2次复数乘法和N(
文档评论(0)