数字信号处理fft

第一篇：数字信号处理fft

数字信号处理实验-FFT的实现

实

验

报

告

学生姓名：

学 号：

指导教师：

一、实验室名称：数字信号处理实验室

二、实验项目名称：FFT的实现

三、实验原理：

一.FFT算法思想：

1.DFT的定义：

对于有限长离散数字信号{x[n]}，0  n  N-1，其离散谱{x[k]}可以由离散付氏变换(DFT)求得。DFT的定义为：

N1X[k]通常令ej2Nx[n]en0j2Nnk，k=0,1,…N-1 WN，称为旋转因子。

2.直接计算DFT的问题及FFT的基本思想：

由DFT的定义可以看出，在x[n]为复数序列的情况下，完全直接运算N点DFT需要(N-1)2次复数乘法和N(N-1)次加法。因此，对于一些相当大的N值(如1024)来说，直接计算它的DFT所作的计算量是很大的。

FFT的基本思想在于，将原有的N点序列分成两个较短的序列，这些序列的DFT可以很简单的组合起来得到原序列的DFT。例如，若N为偶数，将原有的N

22点序列分成两个(N/2)点序列，那么计算N点DFT将只需要约[(N/2) ·2]=N/2次复数乘法。即比直接计算少作一半乘法。因子(N/2)2表示直接计算(N/2)点DFT所需要的乘法次数，而乘数2代表必须完成两个DFT。上述处理方法可以反复使用，即(N/2)点的DFT计算也可以化成两个(N/4)点的DFT(假定N/2为偶数)，从而又少作一半的乘法。这样一级一级的划分下去一直到最后就划分成两点的FFT运算的情况。

3.基2按时间抽取(DIT)的FFT算法思想：

设序列长度为N2L，L为整数(如果序列长度不满足此条件，通过在后面补零让其满足)。

将长度为N2L的序列x[n](n0,1,...,N1)，先按n的奇偶分成两组：

x[2r]x1[r]x[2r1]x2[r]，r=0,1,…,N/2-1 DFT化为：

N1N/21N/21X[k]DFT{x[n]}N/21n0x[n]WnkN2rkr0x[2r]W2rkNr0x[2r1]WN(2r1)kN/21r0N/21x1[r]Wx1[r]W2rkNWWkNr0N/21x2[r]WN

r0rkN/2kNr0x2[r]WN/22rkrk上式中利用了旋转因子的可约性，即：WNN/21NrkN/21rkWN/2。又令

rkX1[k]r0x[1r]W,/X2[k]2r0x[r]WN2，则上式可以写成： /2X[k]X1[k]WNX2[k](k=0,1,…,N/2-1)

k可以看出，X1[k],X2[k]分别为从X[k]中取出的N/2点偶数点和奇数点序列的N/2点DFT值，所以，一个N点序列的DFT可以用两个N/2点序列的DFT组合而成。但是，从上式可以看出，这样的组合仅表示出了X[k]前N/2点的DFT值，还需要继续利用X1[k],X2[k]表示X[k]的后半段本算法推导才完整。利用旋转因子的周期性，有：WN/2WN/2X1[N2N/21rkr(kN/2)，则后半段的DFT值表达式：

rkk]r0x1[r]W2N/2r(Nk)N/21r0x1[r]WN/2X1[k]，同样，X2[N2k]X2[k]

(k=0,1,…,N/2-1)，所以后半段(k=N/2,…,N-1)的DFT值可以用前半段k值表达式获得，中间还利用到WN(N2k)NWN2Wk得到后半段的X[k]值表达式W，

k为：X[k]X1[k]WNkX2[k](k=0,1,…,N/2-1)。

这样，通过计算两个N/2点序列x1[n],x2[n]的N/2点DFTX1[k],X2[k]，可以组合得到N点序列的DFT值X[k]，其组合过程如下图所示：

X1[k] X1[k]WNkX2[k]

X2[k] WNnk -1 X1[k]WNkX2[k]

比如，一个N = 8点的FFT运算按照这种方法来计算FFT可以用下面的流程图来表示：

x(0)W0x(1)W0x(2)W0x(3)W2W0W1W0x(5)W0x(6)W0x(7)W2X(7)W3X(6)W2X(5)X(3)X(2)X(1)X(0)x(4)X(4)

4.基2按频率抽取(DIF)的FFT算法思想：

设序列长度为N2L，L为整数(如果序列长度不满足此条件，通过在后面补零让其满足)。

在把X[k]按k的奇偶分组之前，把输入按n的顺序分成前后两半：

N1N/21nkNN1X[k]DFT{x[n]}N/21N/21x[n]Wn0(nn0N2)kx[n]WnkNnN/2x[n]WNnkn0N/21x[n]WnkNn0x[nNkN2]WNnk

Nn0[x[n]x[nN2NkN2]W2N]WN,k0,1,...,N1因为W2N1，则有WX[k](1)，所以：

kkN/21n0[x[n](1)x[nN2]]WN,k0,1,...,N1

nk按k的奇偶来讨论，k为偶数时：

N/21X[2r]n0[x[n]x[nN2]]WN,k0,1,...,N1 N22rnN/21k为奇数时：X[2r1]前面已经推导过WNN/21n0[x[n]x[n]]WN(2r1)n,k0,1,...,N1

2rkWN/2，所以上面的两个等式可以写为：

N2]]WN/2,r0,1,...,N/21 N2rnrkX[2r]n0[x[n]x[nN/21X[2r1]n0{[x[n]x[n]]WN}WN/2,r0,1,...,N/21

nnr通过上面的推导，X[k]的偶数点值X[2r]和奇数点值X[2r1]分别可以由组合而成的N/2点的序列来求得，其中偶数点值X[2r]为输入x[n]的前半段和后半段之和序列的N/2点DFT值，奇数点值X[2r1]为输入x[n]的前半段和后半段之差再与WN相乘序列的N/2点DFT值。

令x1[n]x[n]x[nN/21nN2]，x2[n][x[n]x[nN/21N2]]WN，则有：

nX[2r]n0x1[n]WrnN/2,X[2r1]n0x2[n]WrnN/2,r0,1,...,N21

这样，也可以用两个N/2点DFT来组合成一个N点DFT，组合过程如下图所示：

x[n] x[n]x[nN2]

x[nN2] -1 WNn [x[n]x[nN2]]WNn

二.在FFT计算中使用到的MATLAB命令：

函数fft(x)可以计算R点序列的R点DFT值;而fft(x,N)则计算R点序列的N点DFT，若R>N，则直接截取R点DFT的前N点，若R

四、实验目的：

离散傅氏变换(DFT)的目的是把信号由时域变换到频域，从而可以在频域分析处理信息，得到的结果再由逆DFT变换到时域。FFT是DFT的一种快速算法。在数字信号处理系统中，FFT作为一个非常重要的工具经常使用，甚至成为DSP运算能力的一个考核因素。

本实验通过直接计算DFT，利用FFT算法思想计算DFT，以及使用MATLAB函数中的FFT命令计算离散时间信号的频谱，以加深对离散信号的DFT变换及FFT算法的理解。

五、实验内容：

a) 计算实数序列x(n)cos516n,0n256的256点DFT。

b) 计算周期为1kHz的方波序列(占空比为50%，幅度取为+/-512，采样频率为25kHz，取256点长度) 256点DFT。

六、实验器材(设备、元器件)：

安装MATLAB软件的PC机一台，DSP实验演示系统一套。

七、实验步骤：

(1) 先利用DFT定义式，编程直接计算2个要求序列的DFT值。

(2) 利用MATLAB中提供的FFT函数，计算2个要求序列的DFT值。 (3) (拓展要求)不改变序列的点数，仅改变DFT计算点数(如变为计算1024点DFT值)，观察画出来的频谱与前面频谱的差别，并解释这种差别。通过这一步骤的分析，理解频谱分辨力的概念，解释如何提高频谱分辨力。

(4) 利用FFT的基本思想(基2-DIT或基2-DIF)，自己编写FFT计算函数，并用该函数计算要求序列的DFT值。并对前面3个结果进行对比。

(5) (拓展要求)尝试对其他快速傅立叶变换算法(如Goertzel算法)进行MATLAB编程实现，并用它来计算要求的序列的DFT值。并与前面的结果进行对比。

(6) (拓展要求)在提供的DSP实验板上演示要求的2种序列的FFT算法(基2-DIT)，用示波器观察实际计算出来的频谱结果，并与理论结果对比。

八、实验数据及结果分析：

程序： (1) 对要求的2种序列直接进行DFT计算的程序

(2) 对要求的2种序列进行基2-DIT和基2-DIF FFT算法程序 (3) 对要求的2种序列用MATLAB中提供的FFT函数进行计算的程序

结果：(1)对2种要求的序列直接进行DFT计算的频域波形

(2)对2种要求的序列进行基2-DIT和基2-DIF FFT算法频域波形 (3)对2种要求的序列用MATLAB中提供的FFT函数计算的频域波形。 (4)(拓展要求)分析利用上面的方法画出的信号频谱与理论计算出来的频谱之间的差异，并解释这种差异。

(5)(拓展要求)保持序列点数不变，改变DFT计算点数(变为1024点)，观察频谱的变化，并分析这种变化，由此讨论如何提高频谱分辨力的问题。

九、实验结论：

十、总结及心得体会：

十一、对本实验过程及方法、手段的改进建议：

第二篇：数字信号处理课程设计

目 录

摘

要........................................................................................................................................... 1 1 绪论 .............................................................................................................................................. 2

1.1 DSP系统特点和设计基本原则 ...................................................................................... 2 1.2 国内外研究动态 ............................................................................................................. 2 2系统设计........................................................................................................................................ 3 3硬件设计........................................................................................................................................ 5

3.1 硬件结构 ........................................................................................................................... 5 3.2 硬件电路设计 ................................................................................................................... 7

3.2.1 总输入电路 ........................................................................................................... 7 3.2.2 总输出电路 ........................................................................................................... 7 3.2.3 语音输入电路 ....................................................................................................... 9 3.2.4 语音输出电路 ....................................................................................................... 9

4 实验结果及分析 ......................................................................................................................... 10 4.1 实验结果 ......................................................................................................................... 10 4.2 实验分析 ......................................................................................................................... 12 5 总结与心得体会 ......................................................................................................................... 13 参考文献......................................................................................................................................... 14 致谢 ................................................................................................................................................ 15

摘

要

基于DSP的语音信号处理系统，该系统采用TMS320VC5509作为主处理器，TLV320AIC23B作为音频芯片，在此基础上完成系统硬件平台的搭建和软件设计，从而实现对语音信号的采集、滤波和回放功能，它可作为语音信号处理的通用平台。

语音是人类相互之间进行交流时使用最多、最自然、最基本也是最重要的信息载体。在高度信息化的今天，语音信号处理是信息高速公路、多媒体技术、办公自动化、现代通信及智能系统等新兴领域应用的核心技术之一。通常这些信号处理的过程要满足实时且快速高效的要求，随着DSP技术的发展，以DSP为内核的设备越来越多，为语音信号的处理提供了良好的平台。本文设计了一个基于TMS320VC5509定点的语音信号处理系统，实现对语音信号的采集、处理与回放等功能，为今后复杂的语音信号处理算法的研究和实时实现提供一个通用平台。

关键词：语音处理;DSP;TMS320VC5509;TLV320AIC23B

1 1 绪论

语音是人类相互间所进行的通信的最自然和最简洁方便的形式，语音通信是一种理想的人机通信方式。语音通信的研究涉及到人工智能、数字信号处理、微型计算机技术、语言声学、语言学等许多领域，所以说语音的通信是一个多学科的综合研究领域，其研究成果具有重要的学术价值。另外通过语音来传递信息是人类最重要的、最有效、最常用的交换信息的形式。语言是人类特有的功能，声音是人类常用的工具，是相互传递信息的主要手段。同时也是众构成思想交流和感情沟通的最主要的途径。

1.1 DSP系统特点和设计基本原则

DSP(digital signal processor)是一种独特的微处理器，是以数字信号来处理大量信息的器件。其工作原理是接收模拟信号，转换为0或1的数字信号。再对数字信号进行修改、删除、强化，并在其他系统芯片中把数字数据解译回模拟数据或实际环境格式。它不仅具有可编程性，而且其实时运行速度可达每秒数以千万条复杂指令程序，远远超过通用微处理器，是数字化电子世界中日益重要的电脑芯片。它的强大数据处理能力和高运行速度，是最值得称道的两大特色。

1.2 国内外研究动态

语音信号处理作为一个重要的研究领域，已经有很长的研究历史。但是它的快速发展可以说是从1940年前后Dudley的声码器和Potter等人的可见语音开始的;20世纪60年代中期形成的一系列数字信号处理的理念和技术基础;到了80年代，由于矢量量化、隐马尔可夫模型和人工神经网络等相继被应用于语音信号处理，并经过不断改进与完善，使得语音信号处理技术产生了突破性的进展。一方面，对声学语音学统计模型的研究逐渐深入，鲁棒的语音识别、基于语音段的建模方法及隐马尔可夫模型与人工神经网络的结合成为研究的热点。另一方面，为了语音识别实用化的需要，讲者自适应、听觉模型、快速搜索识别算法以及进一步的语言模型的研究等课题倍受关注。

在通信越来越发达的当今世界，尤其最近几十年，语音压缩编码技术在移动

2 通信、IP电话通信、保密通信、卫星通信以及语音存储等很多方面得到了广泛的应用。因此，语音编码一直是通信和信号处理的研究热点，并其取得了惊人的进展，目前在PC机上的语音编码已经趋于成熟，而如何在嵌入式系统中实时实现语音压缩编码则是近些年来语音信号处理领域的研究热点之一。

2系统设计

在实际生活中，当声源遇到物体时会发生反射，反射的声波和声源声波一起传输，听者会发现反射声波部分比声源声波慢一些，类似人们面对山体高声呼喊后可以在过一会儿听到回声的现象。声音遇到较远物体产生的反射会比遇到较近的反射波晚些到达声源位置，所以回声和原声的延迟随反射物体的距离大小改变。同时，反射声音的物体对声波的反射能力，决定了听到的回声的强弱和质量。另外，生活中的回声的成分比较复杂，有反射、漫反射、折射，还有回声的多次反射、折射效果。

当已知一个数字音源后，可以利用计算机的处理能力，用数字的方式通过计算模拟回声效应。简单的讲，可以在原声音流中叠加延迟一段时间后的声流，实现回声效果。当然通过复杂运算，可以计算各种效应的混响效果。如此产生的回声，我们称之为数字回声。

本次实验的程序流程图如下：

图2.1 程序流程图

3

本次实验的系统框图如下：

图2.2 系统框图

4

3硬件设计

3.1 硬件结构

图3.1是系统的硬件结构框图, 系统主要包括VC5509和A IC23 两个模块。

图3.1系统硬件结构框图

利用VC5509 的片上外设I2C( Inter - Integrated Circuit, 内部集成电路)模块配置AIC23 的内部寄存器; 通过VC5509 的McBSP (Multi channel Buffered Serial Ports, 多通道缓存串口)接收和发送采样的音频数据。控制通道只在配置AIC23 的内部寄存器时工作, 而当传输音频数据时则处于闲置状态。

AIC23通过麦克风输入或者立体声音频输入采集模拟信号, 并把模拟信号转化为数字信号, 存储到DSP的内部RAM中,以便DSP处理。

当DSP完成对音频数据的处理以后, AIC23再把数字信号转化为模拟信号, 这样就能够在立体声输出端或者耳机输出端听到声音。

5 AIC23能够实现与VC5509 DSP的McBSP端口的无缝连接, 使系统设计更加简单。接口的原理框图, 如下图所示。

图3.2 AIC23与VC5509接口原理图

系统中A IC23的主时钟12 MHz直接由外部的晶振提供。MODE接数字地, 表示利用I2 C控制接口对AIC23传输控制数据。CS接数字地, 定义了I2 C总线上AIC23的外设地址, 通过将CS接到高电平或低电平, 可以选择A IC23作为从设备在I2 C总线上的地址。SCLK和SDIN是AIC23控制端口的移位时钟和数据输入端,分别与VC5509的I2C模块端口SCL和SDA相连。

收发时钟信号CLKX1和CLKR1由A IC23的串行数据输入时钟BCLK提供, 并由A IC23的帧同步信号LRCIN、LRCOUT启动串口数据传输。DX1和DR1分别与A IC23 的D IN 和DOUT 相连, 从而完成VC5509与AIC23间的数字信号通信。

6 3.2 硬件电路设计

3.2.1 总输入电路

图3.3 总输入电路

从左到右各部分电路为：

话筒，开关，语音输入电路，UA741高增益放大电路，有源二阶带 通滤波器。

3.2.2 总输出电路

图3.4 总输出电路

从左到右各部分电路为：

LM386高频功率放大器及其外围器件连接电路，语音输出电路，开关，扬声器。

8 3.2.3 语音输入电路

图3.5语音输入电路

3.2.4 语音输出电路

图3.6 语音输出电路

语音信号通道包括模拟输入和模拟输出两个部分。模拟信号的输入输出电路如图所示。上图中MICBIAS 为提供的麦克风偏压，通常是3/4 AVDD，MICIN为麦克风输入，可以根据需要调整输入增益。下图中LLINEOUT 为左声道输出，RLINEOUT为右声道输出。用户可以根据电阻阻值调节增益的大小，使语音输入输出达到最佳效果。从而实现良好的模拟语音信号输入与模拟信号的输出。

9 4 实验结果及分析

4.1 实验结果

按“F5”键运行，注意观察窗口中的bEcho=0，表示数字回声功能没有激活。这时从耳机中能听到麦克风中的输入语音放送。将观察窗口中bEcho的取值改成非0值。这时可从耳机中听到带数字回声道语音放送。

分别调整uDelay和uEffect的取值，使他们保持在0-1023范围内，同时听听耳机中的输出有何变化。

当uDelay和uEffect的数值增大时，数字回声的效果就会越加的明显。

图4.1 修改前程序图

图4.2 修改前程序图

图4.3 频谱分析

图4.4 左声道及右声道波形

11 4.2 实验分析

所以，从本实验可知当已知一个数字音源后，可以利用计算机的处理能力，用数字的方式通过计算模拟回声效应。简单的讲，可以在原声音流中叠加延迟一段时间后的声流，实现回声效果。当然通过复杂运算，可以计算各种效应的混响效果。

声音放送可以加入数字回声，数字回声的强弱和与原声的延迟均可在程序中设定和调整。

12 5 总结与心得体会

通过本次课程设计，我明白了细节决定成败这句话的道理，在实验中，有很多注意的地方，都被忽视了，导致再花费更多的时间去修改，这严重影响了试验的进度。同时，在本次实验中我了解了ICETEK – VC5509 – A板上语音codec芯片TLV320AIC23的设计和程序控制原理，并进一步掌握了数字回声产生原理、编程及其参数选择、控制，以及了解了VC5509DSP扩展存储器的编程使用方法。

这一学期的理论知识学习加上这次课程设计，使我对DSP有了更加深刻的了解，对数字信号的处理功能，软硬件相结合，语音信号的采集与放送等等方面都有了很深的了解，相信本次课程设计，无论是对我以后的学习，还是工作等方面都有一个很大的帮助。因此，本次课程设计让我受益匪浅。

参考文献

[1]李利.DSP原理及应用[M].北京:中国水利水电出版社,2004. [2]王安民,陈明欣,朱明.TMS320C54xxDSP实用技术[M].北京:清华大学出版社,2007 [3]彭启琮,李玉柏.DSP技术[M].成都:电子科技大学出版社,1997 [4]李宏伟,等.基于帧间重叠谱减法的语音增强方法[J].解放军理工大学学报,2001(1):41~44 [5]TexasInstrumentsIncorporated.TMS320C54x系列DSP的CPU与外设[M].梁晓雯,裴小平,李玉虎,译.北京:清华大学出版社,2006 [6]赵力.语音信号处理[M].北京:机械工业出版社,2003比较图4和图5,可以看到1200Hz以上的频谱明显得到了抑制。

[7]江涛,朱光喜.基于TMS320VC5402的音频信号采集与系统处理[J].电子技术用,2002,28(7):70~72[8]TexasInstrumentsIncorporated:TMS320VC5402Datasheet,2001

致谢

在本次课程设计的即将完成之际，笔者的心情无法平静，本文的完成既是笔者孜孜不倦努力的结果，更是指导老师樊洪斌老师亲切关怀和悉心指导的结果。在整个课程设计的选题、研究和撰写过程中，老师都给了我精心的指导、热忱的鼓励和支持，他的精心点拨为我开拓了研究视野，修正了写作思路，对课程设计的完善和质量的提高起到了关键性的作用。另外，导师严谨求实的治学态度、一丝不苟的工作作风和高尚的人格魅力，都给了学生很大感触，使学生终生受益。在此，学生谨向老师致以最真挚的感激和最崇高的敬佩之情。

另外，还要感谢这段时间来陪我一起努力同学，感谢我们这个小团队，感谢每一个在学习和生活中所有给予我关心、支持和帮助的老师和同学们，几年来我们一起学习、一起玩耍，共同度过了太多的美好时光。我们始终是一个团结、友爱、积极向上的集体。

第三篇：数字信号处理课程总结

以下图为线索连接本门课程的内容：

xa(t)数字信号前置滤波器A/D变换器处理器D/A变换器AF(滤去高频成分)ya(t)x(n)

一、 时域分析

1. 信号

 信号：模拟信号、离散信号、数字信号(各种信号的表示及关系)  序列运算：加、减、乘、除、反褶、卷积  序列的周期性：抓定义

njwna、e(n)(可表征任何序列)cos(wn) u(n)、 典型序列：、、RN(n)、x(n)x(m)(nm)

m特殊序列：h(n) 2. 系统

 系统的表示符号h(n)  系统的分类：y(n)T[x(n)]

线性：T[ax1(n)bx2(n)]aT[x1(n)]bT[x2(n)] 移不变：若y(n)T[x(n)]，则y(nm)T[x(nm)] 因果：y(n)与什么时刻的输入有关 稳定：有界输入产生有界输出

 常用系统：线性移不变因果稳定系统  判断系统的因果性、稳定性方法  线性移不变系统的表征方法：

线性卷积：y(n)x(n)*h(n)

NMk差分方程： y(n)ak1y(nk)bk0kx(nk) 3. 序列信号如何得来?

xa(t)x(n)抽样

 抽样定理：让x(n)能代表xa(t)  抽样后频谱发生的变化?  如何由x(n)恢复xa(t)?

sin[xa(mT)T(tmT)]

xa(t)=mT

(tmT)

二、 复频域分析(Z变换)

时域分析信号和系统都比较复杂，频域可以将差分方程变换为代数方程而使分析简化。 A. 信号 1.求z变换

定义：x(n)X(z)x(n)znn

收敛域：X(z)是z的函数，z是复变量，有模和幅角。要其解析，则z不能取让X(z)无穷大的值，因此z的取值有限制，它与x(n)的种类一一对应。

 x(n)为有限长序列，则X(z)是z的多项式，所以X(z)在z=0或∞时可能会有∞，所以z的取值为：0z;

 x(n)为左边序列，0zRx，z能否取0看具体情况;

 x(n)为右边序列，Rxz，z能否取∞看具体情况(因果序列);  x(n)为双边序列，RxzRx 2.求z反变换：已知X(z)求x(n)

 留数法

 部分分式法(常用)：记住常用序列的X(z)，注意左右序列区别。  长除法：注意左右序列 3.z变换的性质：

 由x(n)得到X(z)，则由x(nm)zmX(z)，移位性;  初值终值定理：求x(0)和x();

 时域卷积和定理：y(n)x(n)*h(n)Y(z)X(z)H(z);  复卷积定理：时域的乘积对应复频域的卷积;  帕塞瓦定理：能量守恒



nx(n)212X(ejw)dw2

4.序列的傅里叶变换

公式：X(ejw)x(n)enjwn

x(n)12X(ej)ejnd

注意：X(ejw)的特点：连续、周期性;X(ejw)与X(z)的关系 B. 系统

由h(n)H(z)，系统函数，可以用来表征系统。

 H(z)的求法：h(n)H(z);H(z)=Y(z)/X(z);  利用H(z)判断线性移不变系统的因果性和稳定性  利用差分方程列出对应的代数方程

MNMy(n)ak1y(nk)kbk0x(nk)kY(z)X(z)bk0Nkzk

k1ak1zk 系统频率响应H(ejw)：以2为周期的的连续函数



H(e)jwh(n)enjwn

H(ejw)h(n)enjwn，当h(n)为实序列时，则有H(ejw)=H*(ejw)

三、 频域分析

根据时间域和频域自变量的特征，有几种不同的傅里叶变换对

 时间连续，非周期频域连续(由时域的非周期造成)，非周期(由时域的连续造成); X(j)x(t)ejtdt

x(t)12X(j)ejtd

 时间连续，周期频域离散，非周期

X(jk0)1T0T0/2x(t)ejk0tdt

T0/2x(t)X(jk0)ejk0t

 时间离散，非周期频域连续，周期



X(e)jwx(n)enjwn

x(n)12X(ej)ejnd，wT(数字频率与模拟频率的关系式)

 时间离散，周期频域离散，周期

~X(k)N1n0~x(n)ej2Nkn~x(n)W

knNn0N11~x(n)NN1n0~X(k)ej2Nkn1NN1n0~knX(k)WN

 本章重点是第四种傅里叶变换-----DFS  注意：

x(n)和X(k)都是以N为周期的周期序列; 1)~x(n)和X(k)的定义域都为(,)

2)尽管只是对有限项进行求和，但~;

~~~例如：k0时，X(0)N1x(n)

n0~~k1时，X(1)N1n0~x(n)ej2Nn

2NNnN1~kN时，X(N)N1n0j~x(n)en02N~~x(n)=X(0)

~kN1时，X(N1)N1n0~x(n)ej(N1)n~X(1)

x(n)也有类似的结果。x(n)和X(k)一

同理也可看到~可见在一个周期内，~~一对应。

 比较X(e)jwx(n)enjwn~和X(k)N1n0~x(n)ej2Nkn~x(n)W，当x(n)knNn0N1x(n)的一个周期内有定义时，即x(n)=~x(n)，0nN1，则在只在~N12Nj2Nk时，X(ejw)X(k)。

1,kr 0,kr~ en0(kr)nx(n)和X(k)的每个周期值都只是其主值区间的周期延拓，所以求和 因为~~在任一个周期内结果都一样。

 DFT：有限长序列x(n)只有有限个值，若也想用频域方法分析，它只属于序列的傅里叶变换，但序列的傅氏变换为连续函数，所以为方便计算机处理，也希望能像DFS一样，两个域都离散。将x(n)想象成一个周期x(n)的一个周期，然后做DFS，即 序列~

~X(k)N1n0~x(n)ej2NknN1n0x(n)ej2Nkn

x(n)只有x(n)，不是真正的周期序列，但因为求和只需N注意：实际上~个独立的值，所以可以用这个公式。同时，尽管x(n)只有N个值，但依上式求出的X(k)还是以N为周期的周期序列，其中也只有N个值独立，这样将~X(k)规定在一个周期内取值，成为一个有限长序列，则会引出

N1j2Nkn~DFT X(k)x(n)en0RN(k)

x(n)1NN1n0X(k)ej2NknRN(n)

比较：三种移位：线性移位、周期移位、圆周移位

三种卷积和：线性卷积、周期卷积、圆周卷积

重点：1)DFT的理论意义，在什么情况下线性卷积=圆周卷积 2)频域采样定理：掌握内容，了解恢复

3)用DFT计算模拟信号时可能出现的几个问题，各种问题怎样引起?

混叠失真、频谱泄漏、栅栏效应

 FFT：为提高计算速度的一种算法

1) 常用两种方法：按时间抽取基2算法和按频率抽取基2算法，各自的原理、特点是什么，能自行推导出N小于等于8的运算流图。 2) 比较FFT和DFT的运算量; 3) 比较DIT和DIF的区别。

四、 数字滤波器(DF)

一个离散时间系统可以用h(n)、H(z)、差分方程和H(ejw)来表征。 问题：

1、各种DF的结构

2、如何设计满足要求指标的DF?

3、如何实现设计的DF?

A. 设计IIR DF，借助AF来设计，然后经S---Z的变换即可得到。

1) 脉冲响应不变法：思路、特点 2) 双线性变换法：思路、特点、预畸变 3) 模拟滤波器的幅度函数的设计 B. 设计FIR DF 1) 线性相位如何得到?条件是什么?各种情况下的特点。 2) 窗函数设计法：步骤、特点 3) 频率抽样法：步骤、特点 C. 实现DF

Ma

标准形式：H(z)k0Nkzk

bkzk1k1

第四篇：数字信号处理学习心得

通信工程 0801 赖立根

《数字信号处理》是我们通信工程和电子类专业的一门重要的专业基础课程，主要任务是研究数字信号处理理论的基本概念和基本分析方法，通过建立数学模型和适当的数学分析处理，来展示这些理论和方法的实际应用。

数字信号处理技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连;它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。信息科学是研究信息的获取、传输、处理和利用的一门科学，信息要用一定形式的信号来表示，才能被传输、处理、存储、显示和利用，可以说，信号是信息的表现形式，而信息则是信号所含有的具体内容。

一单元的课程我们深刻理解到时域离散信号和时域离散系统性质和特点;时域离散信号和时域离散系统时域分析方法;模拟信号的数字处理方法。

二单元的课程我们理解了时域离散信号(序列)的傅立叶变换，时域离散信号Z变换，时域离散系统的频域分析。

三单元的课程我们学习了离散傅立叶变换定义和性质，离散傅立叶变换应用——快速卷积，频谱分析。

四单元的课程我们重点理解基2 FFT算法——时域抽取法﹑频域抽取法，FFT的编程方法，分裂基FFT算法。

五单元的课程我们学了网络结构的表示方法——信号流图，无限脉冲响应基本网络结构，有限脉冲响应基本网络结构，时域离散系统状态变量分析法。

六单元的课程我们理解数字滤波器的基本概念，模拟滤波器的设计,巴特沃斯滤波器的设计,切比雪夫滤波器的设计，脉冲响应不变法设计无限脉冲响应字数字滤波器，双线性变换法设计无限脉冲响应字数字滤波器，数字高通﹑带通﹑带阻滤波器的设计。

七单元的课程我们学习了线性相位有限脉冲响应(FIR)数字滤波器，窗函数法设计有限脉冲响应(FIR)数字滤波器，频率采样法设计有限脉冲响应(FIR)数字滤波器

通信工程是一门工程学科，主要是在掌握通信基本理论的基础上，运用各种工程方法对通信中的一些实际问题进行处理。通过该专业的学习，可以掌握电话网、广播电视网、互联网等各种通信系统的原理，研究提高信息传送速度的技术，根据实际需要设计新的通信系统，开发可迅速准确地传送各种信息的通信工具等。

对于我们通信专业，我觉得是个很好的专业，现在这个专业很热门，这个专业以后就业的方向也很多，就业面很广。我们毕业以后工作，可以进入设备制造商、运营商、专有服务提供商以及银行等领域工作。当然,就业形势每年都会变化，所以关键还是要看自己。可以从事硬件方面，比如说PCB，别小看这门技术，平时我们在试验时制作的简单，这一技术难点就在于板的层数越多，要做的越稳定就越难，这可是非常有难度的，如果学好了学精了，也是非常好找工作的。也可以从事软件方面，这实际上要我们具备比较好的模电和数电的基础知识。

我选择了这个专业，在这里读了三年关于通信知识的书，我还是想以后毕业能够从事这个方面的工作，现在学了通信原理、数字信号处理这些很有用的专业课，所以，我在以后的学习中，我会把这些方面的知识学扎实，从事技术这一块要能吃苦，我也做好了准备，现在还很年轻，年轻的时候多吃点苦没什么，为了我自己美好的将来，我会努力学好这个专业的。

数字信号处理课程属于专业基础课，所涵盖的内容主要有：离散时间信号与系统的基本概念及描述方法，离散傅立叶变换及快速傅立叶变换，数字滤波器结构及设计等。对于电气信息类专业的学生来说，这些内容是学习后续专业课程的重要基础，也是实际工作中必不可少的专业基础知识。目前几乎所有的高等院校都在电子工程类、信息工程类、通信工程类、电子技术类、自动控制类、电气工程类、机电工程类、计算机科学类等工科电类及其他相关专业的本科生中开设了该门课程。随着计算机技术、微电子技术、数字信号处理理论和方法的发展，半个世纪以来，尤其是最近的三十来年里，数字信号处理的方法和应用得到了飞跃式的发展，数字信号处理的地位和作用变得越来越重要。因此，加强该课程的建设具有重要的意义。。

我们的数字信号处理课是罗老师教的，罗老师有过实际工作的经验，对于这门课的实际用途很了解，罗老师对于这门课采用多种教学方法，丰富教学内容，吸引学生对课程的关注。利用实验课使学生亲自编程，体会信号处理课程的乐趣，这样子激发了学生的兴趣、提高了教学的效果。因此，我们班的同学在这一个学期的学习中，这门课都学的比较好。

数字信号处理课程的特点是课程本身理论性强、公式推导较多、概念比较抽象，学生常有枯燥难学之感。近年来，国外及国内有些学校对一般电类专业该课程的教学主要强调应用性学习，主要介绍数字信号处理的用途和用法，而对其深奥的理论推导仅做一般介绍，并给学生提供进行实验的机会，以激发学生对该课程的兴趣和学习主动性。

对该课程的改革思想主要是课程内容要适应数字信号处理技术的发展现状，淡化枯燥的数学推导，辅助以现代化教学手段，并开设相应的实验课。结合专业现状，将课堂教学一部分变为多媒体教学，尽量将一些理论分析用图形手段展示出来，以增强学生的感性认识。实验课主要是以MATLAB为平台，充分利用MATLAB的数字信号处理工具箱提供的各种功能让学生亲自动手将课堂所学进行仿真实现。实验课还可以通过用DSP试验箱实现数字信号处理的功能向学生进行演示。

第五篇：数字信号处理实验讲稿

邯 郸 学 院

讲 稿

2010 ～2011 学年 第 一 学期

分院(系、部)： 信息工程学院 教 研 室： 电子信息工程 课 程 名 称： 数字信号处理

授 课 班 级： 07级电子信息工程

主 讲 教 师： 王苗苗 职

称：

助教(研究生)

使 用 教 材： 《数字信号处理》

制 作 系 统：

Word2003

邯郸学院制

实验一.. Matlab仿真软件介绍

一、实验目的

熟悉Matlab仿真软件

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、学习Matlab仿真软件的安装

2、熟悉Matlab仿真软件的操作环境

3、直接在Matlab仿真软件的命令窗口实现数值计算

4、编写M文件

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

1、熟悉Matlab仿真软件

2、参阅Matlab及在电子信息类课程中的应用(第2版)唐向宏 电子工业出版社

实验二 离散信号和系统分析的Matlab实现

一、实验目的

1、Matlab实现离散信号和系统分析

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、利用Matlab产生离散信号

2、利用Matlab计算离散卷积

3、利用Matlab求解离散LTI系统响应

4、利用Matlab计算DTFT

5、利用Matlab实现部分分式法

6、利用Matlab计算系统的零极点

7、利用Matlab进行简单数字滤波器设计

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验三 利用Matlab实现信号DFT的计算

一、实验目的

1、Matlab实现信号DFT的计算

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、利用Matlab计算信号的DFT

2、利用Matlab实现由DFT计算线性卷积

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验四 利用Matlab实现滤波器设计

一、实验目的

1、Matlab实现实现滤波器设计

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、 利用Matlab实现模拟低通滤波器的设计

2、 利用Matlab实现模拟域频率变换

3、 利用Matlab实现脉冲响应不变法

4、 利用Matlab实现双线性变换法

5、 利用Matlab实现数字滤波器设计

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验五 利用Matlab实现FIR滤波器设计

一、实验目的

1、Matlab实现实现滤波器设计

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、利用Matlab实现窗函数法

2、利用Matlab实现频率取样法

3、利用Matlab实现优化设计

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验六.. 随机信号功率谱估计的Matlab实现

一、实验目的

1、Matlab实现实现滤波器设计

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、利用Matlab实现随机序列

2、利用Matlab计算相关函数的估计

3、利用Matlab进行非参数功率谱估计

4、利用Matlab进行AR模型功率谱估计

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验七.. 数字滤波器结构的Matlab实现

一、实验目的

1、Matlab实现实现滤波器设计

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、利用Matlab实现数字滤波器直接型设计

2、利用Matlab实现数字滤波器级联设计

3、利用Matlab实现数字滤波器并联型设计

4、利用Matlab实现数字滤波器格型设计

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容

实验八.... 利用Matlab实现信号小波分析

一、实验目的

1、Matlab实现实现滤波器设计

2、进一步熟悉Matlab软件操作

二、实验设备和元器件

含Matlab仿真软件的计算机

三、实验内容和步骤

1、小波测试信号

2、分解与重构滤波器组

3、离散小波变换

4、离散小波反变换

5、基于小波的信号去噪

6、基于小波的信号压缩

四、实验报告要求

按照《Matlab程序设计》模板提交实验报告

五、预习要求

预习课本上的相关内容