数字信号处理复习2

绪论：本章介绍数字信号处理课程的基本概念。

0.1信号、系统与信号处理

1.信号

信号是信息的载体。分类：模拟信号、量化信号、抽样信号和数字信号。 2.系统

系统定义为处理（或变换）信号的物理设备，或者说，凡是能将信号加以变换以达到人们要求的各种设备都称为系统。 3.信号处理

信号处理即是用系统对信号进行某种加工。包括：滤波、分析、变换、综合、压缩、估计、识别等等。

所谓“数字信号处理”，就是用数值计算的方法，完成对信号的处理。通过处理，往往可以达到两个目的：

0.2数字信号处理系统的基本组成

0.3数字信号处理的特点

（1）灵活性。

（2）高精度和高稳定性。 （3）便于大规模集成。

（4）对数字信号可以存储、运算、系统可以获得高性能指标。

0.4数字信号处理基本学科分支

数字信号处理（DSP）一般有两层含义，一层是广义的理解，为数字信号处理技术——Digital Signal Processing，另一层是狭义的理解，为数字信号处理器——Digital Signal Processor。

0.5课程内容

该课程在本科阶段主要介绍以傅里叶变换为基础的“经典”处理方法，包括： （1）离散傅里叶变换及其快速算法。

（2）滤波理论（线性时不变离散时间系统，用于分离相加性组合的信号，要求信号频谱占据不同的频段）。

在研究生阶段相应课程为“现代信号处理”（Advanced Signal Processing）。信号对象主要是随机信号，主要内容是自适应滤波（用于分离相加性组合的信号，但频谱占据同一频段）和现代谱估计。

1

一、重点与难点 1.信号及其分类；

2.数字信号处理系统的基本组成。 二、具体讲解 1．信号及其分类

信号是信息的载体，以某种函数的形式传递信息。这个函数可以是时间域、频率域或其它域，但最基础的域是时域。 分类：

周期信号/非周期信号 确定信号/随机信号 能量信号/功率信号

连续时间信号/离散时间信号/数字信号 按自变量与函数值的取值形式不同分类：

2．数字信号处理系统的基本组成

数字信号处理就是用数值计算的方法对信号进行变换和处理。不仅应用于数字化信号的处理，而且也可应用于模拟信号的处理。以下讨论模拟信号数字化处理系统框图。

（1）前置滤波器

将输入信号xa(t)中高于某一频率（称折叠频率，等于抽样频率的一半）的分量加以滤除。

（2）A/D变换器

在A/D变换器中每隔T秒（抽样周期）取出一次xa(t)的幅度，抽样后的信号称为离散信号。在A/D变换器中的保持电路中进一步变换为若干位码。 （3）数字信号处理器（DSP） （4）D/A变换器

2

按照预定要求，在处理器中将信号序列x(n)进行加工处理得到输出信号y(n)。由一个二进制码流产生一个阶梯波形，是形成模拟信号的第一步。 （5）模拟滤波器

把阶梯波形平滑成预期的模拟信号；以滤除掉不需要的高频分量，生成所需的模拟信号ya(t)。

1．周期信号的傅里叶级数

复指数函数集是一种常见的完备正交函数集，周期信号可以表示为复指数函数的线性组合。

一般是复函数。 2．傅里叶变换 （1）定义

电信号：时域 频域 f(t) F(Ω)

f(t)与F(Ω)之间的关系——傅里叶变换

正变换 积分因子： 反变换

积分因子：

3

（2）傅里叶变换性质

3．冲激函数δ(t) （1）定义

（2）性质

δ(t)＝0， 当t≠0时

抽样性： 搬移性： 4．卷积 卷积积分

卷积的傅里叶变换 时域卷积

4

频域卷积

第一章：本章概念较多，需要理解和识记的内容较多，学习时要注意。

1.2离散时间信号

1.离散时间信号的定义

离散时间信号是指一个实数或复数的数字序列，它是整数自变量n的函数，表示为

x(n)。离散时间信号也常用图形描述。

2.几种基本离散时间信号

5

（1）单位采样序列

（2）单位阶跃序列

（3）矩形序列

（4）实指数序列 （5）正弦序列

ω是正弦序列数字域的频率，单位是弧度。

对连续信号中的正弦信号进行采样，可得正弦序列。设连续信号为它的采样值为

，因此

，

这个式子具有一般性，它反映了由连续信号采样得到的离散序列，其数字频率与模拟频率的一般关系。另外需要说明的是，ω的单位为弧度，Ω的单位为弧度/秒。本书中，我们一律以ω表示数字域频率，而以Ω及f表示模拟域频率。 （6）复指数序列

复指数序列是以余弦序列为实部、正弦序列为虚部所构成的一个复数序列。 3.信号运算

（1）加法：两个信号之和

由同序号的序列值逐点对应相加得到。 （2）乘法：两个信号之积

由同序号的序列值逐点对应相乘得到。

6

（3）移位： 当

，序列左移（称为超前）。

，序列右移（称为延时）；当

（4）翻转：

（5）尺度变换： 当为下采样。

时，序列

或，其中M和N都是正整数。

是通过取x(n)的每第M个采样形成，这种运算称

对于序列，定义如下

这种运算称为上采样。 4.信号分解

任一信号x(n)可表示成单位脉冲序列的移位加权和 简记为

一、重点与难点 1.几种常用的信号； 2.公式????的含义；

3.线性、时不变、因果和稳定系统的判别； 4.线性卷积的计算；

5.采样的框图、时域采样定理及信号内插恢复的过程。 二、具体讲解 1．线性卷积

线性卷积是一种非常重要的一种运算，对它的求解，一般我们采用作图法。线性卷积满足交换律，设两序列长度分别是N和M，线性卷积后序列的长度为N＋M－1。 卷积的计算过程包括翻转、移位、相乘、相加四个过程。

7

1）将和用和表示，画出和这两个序列； ；

2）选择一个序列 3）将 4）将

，并将其按时间翻转形成序列

；

移位n，得到和

相同m的序列值对应相乘后，再相加。

2．连续信号的采样

对连续信号进行理想采样，设采样脉冲，则采样输出

在讨论理想采样后，信号频谱发生的变化时，可遵循下面的思路： 1）由

；2）由

；

3）根据频域卷积定理，由 计算过程： 1）

计算出。

2）周期信号可以用傅里叶级数展开，因此

其中系数

所以

其傅里叶变换

8

3）

因此，采样后信号频谱产生周期延拓，周期为Ωs，同时幅度为原来的1/T倍。这是一个非常重要的性质，应熟练掌握。 例题

1.用单位脉冲序列及其加权和表示图所示序列

解：

2.判别系统y(n) =T[x(n)]=ax(n)+ b是否为线性系统，是否为时不变系统？

解：（1）线性 T[x1(n)]= ax1(n)+b T[x2(n)]= ax2(n)+b

而T[x1(n)+x2(n)]=a[x1(n)+x2(n)]+b≠ax1(n)+ b+ax2(n)+ b 故此系统不是线性系统。 （2）时不变性 T[x(n-n0)]=ax(n-n0)+b

y(n-n0) = ax (n-n0)+b= T[x(n-n0)]

故该系统是时不变系统。

3.判别系统y(n) =T[x(n)]=x(n)cos(ω0n+φ)的因果稳定性。

解：（1）因果性

因为y(n) =T[x(n)]=x(n)cos(ω0n+φ)只与x(n)的当前值有关，而与x(n+1)，x(n+2)……等未来值无关，故系统是因果的。 （2）稳定性

当| x(n)|

9

4.若LTI系统的输入x(n)和输出y(n)满足下列差分方程

y(n)=ay(n-1)+x(n)

求起始条件分别为h(n)=0，n<0和h(n)=0，n>0时的单位脉冲响应。 解：（1）令x(n)=δ(n)，根据起始条件可递推如下 y(0)=δ(0)=1，y(1)=ay(0)=a，……y(n)=ay(n-1)=a^-n 因此h(n)= y(n) =a^-n.u(n) （2）将差分方程改写成 y(n-1)=1/a[y(n)-x(n)]

n→n+1，则y(n)=1/a[y(n+1)-x(n+1)] 根据起始条件可递推如下

y(0)=1/a[y(1)-δ(1)]=0，y(-1)=1/a[y(0)-δ(0)]=-1/a，……y(n)=ay(n-1)=-a^-n 因此h(n)= y(n) =-a^-n.u(-n-1)

第二章：本章涉及信号及系统的频域分析方法，概念较多，但很基础，学习时要注意。

2.1序列的傅里叶变换的定义及性质

1.定义

若序列满足绝对可和条件

则其傅里叶变换（Discrete Time Fourier Transform-DTFT）定义为

10

如果对X(z)单位圆上进行等距离采样

现在要问，这样采样以后，信息有没有损失？或者说，采样后所获得的有限长序列xN(n)能不能代表原序列x(n)。

开始

为了弄清这个问题，我们从周期序列

由于

所以

也即

是原非周期序列x(n)的周期延拓序列，其时域周期为频域采样点数N。

在第一章我们看到，时域的采样造成频域的周期延拓，这里又对称的看到，频域采样同样造成时域的周期延拓。

因此，如果序列x(n)不是有限长的，则时域周期延拓时，必然造成混叠现象，因而一定会产生误差。

对于长度为M的有限长序列，只有当频域采样点数N大于或等于序列长度M时，才有

即可由频域采样值X(k)恢复出原序列x(n)，否则产生时域混叠现象，这就是所谓的频域采样定理。

2.用DFT进行谱分析的误差问题 （1）混叠现象

16

利用DFT逼近连续时间信号的傅里叶变换，为避免混叠失真，按照抽样定理的要求，采样频率至少是信号最高频率的两倍。

解决混叠问题的唯一方法是保证采样频率足够高。 （2）频谱泄露

任何带限信号都是非时限的，任何时限信号都是非带限的。实际问题中遇到的离散时间序列可能是非时限的、无限长序列，在对该序列利用DFT进行处理时，由于作DFT的点数总是有限的，因此就有一个必须将该序列截断的问题。序列截断的过程相当于给该序列乘上一个矩形窗口函数RN(n)。如果原来序列的频谱为数的频谱为

，则截断后有限长序列的频谱为

的频谱“泄露”

，矩形窗函

由于矩形窗函数频谱的引入，使卷积后的频谱被展宽了，即到其它频率处，称为频谱泄露。

在进行DFT时，由于取无限个数据是不可能的，所以序列的时域截断是必然的，泄露是难以避免的。为了尽量减少泄露的影响，截断时要根据具体的情况，选择适当形状的窗函数，如汉宁窗或汉明窗等。 （3）栅栏效应

由于DFT是有限长序列的频谱等间隔采样所得到的样本值，这就相当于透过一个栅栏去观察原来信号的频谱，因此必然有一些地方被栅栏所遮挡，这些被遮挡的部分就是未被采样到的部分，这种现象称为栅栏效应。由于栅栏效应总是存在的，因而可能会使信号频率中某些较大的频率分量由于被 “遮挡”而无法得到反映。此时，通常在有限长序列的尾部增补若干个零值，借以改变原序列的长度。这样对加长的序列作DFT时，由于点数增加就相当于调整了原来栅栏的间隙，可以使原来得不到反映的那些较大的频率分量落在采样点上而得到反映。 例题

1.设x(t)的最高频率fc不超过3Hz，现用fs＝10Hz，即T＝0.1s对其抽样，求

所得到的频率最大分辨率。如果信号x(n)由三个正弦组成，其频率分别f1＝2Hz，f2＝2.02Hz，f3＝2.07Hz，即

那么用DFT求其频谱时，能否分辨出三个频率分量？

解：x(t)的最高频率fc不超过3Hz，现用fs＝10Hz对其抽样，由抽样定理可知，不会发生混叠问题。Tp＝25.6s，对x(n)做DFT时，所得到的频率最大分辨率

17

如果信号x(n)由三个正弦组成，其频率分别f1＝2Hz，f2＝2.02Hz，f3＝2.07Hz，那么用DFT求其频谱时，由于又由于

，所以不能分辨出由f2产生的正弦分量；

，所以能分辨由f3产生的正弦分量。

2.对有限长序列

，即

，

的Z变换在单位圆上进行5等份采

样，得到采样值

试根据频率采样定理求的逆离散傅里叶变换。

解：

第四章：快速傅里叶变换并不是一种新的变换，它是离散傅里叶变换的一种快速算法。

4.2直接计算DFT的问题及改进的途径

直接计算DFT，需要 （1）

的对称性

次复数乘法，

次复数加法。

18

（2）

的周期性

利用

的对称性和周期性，将大点数的DFT分解成若干个小点数的DFT，FFT

正是基于这个基本思路发展起来的。

分类：按时间抽取（DIT）算法和按频率抽取（DIF）算法。 一、重点与难点

1.DFT提高运算量的途径；

2.基2FFT的算法原理和FFT运算特点； 3.实序列FFT算法思路； 二、具体讲解

1.DFT提高运算量的途径

直接计算离散傅里叶变换，由于计算量近似正比于N2，显然对于很大的N值，直接计算离散傅里叶变换要求的算术运算量非常大。我们可以利用系数WNnk的特性来改善离散傅里叶变换的计算效率。 （1） （2） 利用

的周期性

的对称性和周期性，可以将大点数的DFT分解成若干个小点数的DFT，的对称性

快速傅里叶变换正是基于这个基本思路发展起来的。FFT算法基本上可分为两大类，即按时间抽取（DIT）算法和按频率抽取（DIF）算法。 2.基2FFT的算法原理和FFT运算特点 对于算法原理，要求能够看懂分解流图。 例题

若已知有限长序列

（

和

19

，画出其按时间抽取的基2FFT流图，）的值。

并按FFT运算流程计算

解：（1）基2FFT流图

（2）

，

，

解：

2.画一个按时间抽取4点序列的基2FFT流图。在图上标明时域、频域各输入、

输出项的排列顺序，并标出由第4根水平线（从上往下数）发出的所有支路的系数。

第五章：本章讲授了设计IIR滤波器常用的两种设计方法——脉冲响应不变法和双线性变换法。

20

5.1引言

1.数字滤波器的分类

（1）IIR和FIR数字滤波器

这是根据滤波器的单位脉冲响应h(n)的长度是否有限来划分的。若h(n)是一个长度为M+1的有限长序列，通常将此时的系统称为有限长单位脉冲响应（FIR，Finite Impulse Response）系统。

如果系统函数的分母中除a0外，还有其它的ak不为零，则相应的h(n)将是无限长序列，称这种系统为无限长单位脉冲响应（IIR，Infinite Impulse Response）系统。 （2）低通、高通、带通、带阻滤波器

21

注意：数字滤波器（DF）与模拟滤波器（AF）的区别 数字滤波器的频率响应都是以2π为周期的，滤波器的低通频带处于2π的整数倍处，而高频频带处于π的奇数倍附近。 2.设计指标描述

滤波器的指标通常在频域给出。数字滤波器的频率响应常表示为 其中，

称为幅频响应，

称为相频响应。对IIR数字滤波器，通常用幅

一般为复函数，通

频响应来描述设计指标，而对于线性相位特性的滤波器，一般用FIR滤波器设计实现。 IIR低通滤波器指标描述：

——通带截止频率，

——阻带截止频率，

——通带最大衰减，

——阻

带最小衰减，——3dB通带截止频率

3.设计方法 三步：（1）按照实际需要确定滤波器的性能要求。

（2）用一个因果稳定的系统函数去逼近这个性能要求。 （3）用一个有限精度的算法去实现这个系统函数。

IIR滤波器常借助模拟滤波器理论来设计数字滤波器，设计步骤为：先根据所给的滤波器性能指标设计出相应的模拟滤波器传递函数Ha(s)，然后由ha(s)经变换得到所需的数字滤波器的系统函数H(z)。在变换中，一般要求所得到的数字滤波器频率响应应保留原模拟滤波器频率响应的主要特性。为此要求：

（1）因果稳定的模拟滤波器必须变成因果稳定的数字滤波器； （2）数字滤波器的频响应模仿模拟滤波器的频响。 工程上常用的方法：脉冲响应不变法和双线性变换法。 一、重点与难点

1.数字滤波器的分类，理想滤波器的幅频特性，性能指标； 2.模拟滤波器的分类、特点和设计方法；

3.脉冲响应不变法和双线性变换法的变换原理，s域和z域的映射关系，优缺点； 4.数字滤波器的频率变换。

22

二、具体讲解 1.预畸变

对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但

是各个分段边缘的临界频率点产生了畸变。比如说，经过转换后可能等于

3。这样，如果设计模拟滤波器时按照波器原型经双线性变换后得到的数字滤波器其

来计算，则这样设计出来的模拟滤必然不是2，而可能是（比如说）

1.5，于是数字滤波器的性能与模拟原型的性能就有了失真，这种失真称为频率失真。可以看出，这种失真是双线性变换所固有的，这是双线性变换法的一个很大的缺点。 怎样弥补这个缺点呢？一般的做法是，在给定数字滤波器的通带截止频率带截止频率

和阻

后，我们并不直接按照这个给定的数据去设计原型模拟滤波器，而是先

根据的关系求出相应的和，然后根据和设计模拟原型，求

出其传输函数，最后再利用代入求出数字滤波器的系统函数

。这种方法称为预畸变。

23

因此，预畸变就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数字频率上，此时数字滤波器将满足要求的技术指标。 2.数字滤波器的频率变换

设计IIR数字滤波器时我们常常借助于模拟滤波器，即先将所需要的数字滤波器技术要求转换为一个低通模拟滤波器的技术要求，然后设计这个原型低通模拟滤波器。在得到低通模拟滤波器的传输函数H(p)（或H(s)）后，再变换为所需要的数字滤波器的系统函数H(z)。

将H(p)（或H(s)）变换为H(z)的方法有两种。一种是先将设计出来的模拟低通原型通过频率变换变换成所需要的模拟高通、带通或带阻滤波器，然后再利用脉冲响应不变法或双线性变换法将其变换为相应的数字滤波器，这种方法的频率变换是在模拟滤波器之间进行的。

另一种方法是先将设计出来的模拟低通原型通过脉冲响应不变法或双线性变换法转换为归一化数字低通，最后通过频率变换把数字低通变换成所需要的数字高通、带通或带阻滤波器，这种方法的频率变换是在数字滤波器之间进行的。

24

例题

1.图示是由RC组成的模拟滤波器

（1）写出传输函数；

转换成数字滤波器

，设采样间隔为T；

（2）选用一种合适的转换方法将

（3）比较脉冲响应不变法和双线性变换法的优缺点。 解：（1）由图可知该滤波器为模拟高通滤波器

（2）采用双线性变换法，因为选用脉冲响应不变法，会在高频处发生频率混叠现象。

（3）脉冲响应不变法： 优点：时域逼近良好，

；

缺点：容易产生混叠失真，只适用于带限滤波器； 双线性变换法：

优点：设计运算简单；避免了频谱的混叠效应，适合各种类型滤波器；

缺点：

，会产生非线性频率失真。

2.如图所示系统

（1）写出该系统的系统函数通带滤波器？

25

，画出系统的幅频响应，并问这一系统是哪一种

本文来源：https://www.bwwdw.com/article/uewh.html