数字信号处理期末试卷(含答案)全..

第一篇：数字信号处理期末试卷(含答案)全..

数字信号处理期末试卷(含答案)

一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。

1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过()即可完全不失真恢复原信号。

A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？()A.y(n)=x3(n)B.y(n)=x(n)x(n+2)C.y(n)=x(n)+2

D.y(n)=x(n2)3.．设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取()。A．M+N B.M+N-1

C.M+N+1

D.2(M+N)4.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是()。

A.N≥M B.N≤M C.N≤2M D.N≥2M 5.直接计算N点DFT所需的复数乘法次数与()成正比。A.N B.N2 C.N3 D.Nlog2N 6.下列各种滤波器的结构中哪种不是FIR滤波器的基本结构()。A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR滤波器的幅度响应H(w)特点()： A 关于w0、、2偶对称

B 关于w0、、2奇对称

C 关于w0、2偶对称 关于w奇对称

D关于w0、2奇对称 关于w偶对称 8.适合带阻滤波器设计的是：（）A h(n)h(N1n)N为偶数 B h(n)h(N1n)N为奇数 C h(n)h(N1n)N为偶数 D h(n)h(N1n)N为奇数

9.以下对双线性变换的描述中不正确的是()。A.双线性变换是一种非线性变换

B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内 D.以上说法都不对

10.关于窗函数设计法中错误的是：

A窗函数的截取长度增加，则主瓣宽度减小；

B窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的截取长度无关； C为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加； D窗函数法不能用于设计高通滤波器；

二、填空题(每空2分，共20分)1.用DFT近似分析连续信号频谱时, _________效应是指DFT只能计算一些离散点上的频谱。

2.有限长序列X(z)与X（k）的关系 X（k）与X(ejw)的关系 3.下图所示信号流图的系统函数为：

4.如果通用计算机的速度为平均每次复数乘需要４μs，每次复数加需要1μs，则在此计算机上计算210点的基2FFT需要__________级蝶形运算,总的运算时间是__________μs。

5.单位脉冲响应不变法优点 , 缺点____________,适合_______________________滤波器设计

6.已知FIR滤波器H(z)12z15z2az3z4具有线性相位，则a＝______,冲激响应h（2）＝___,相位(w)___ 3n)的周期__________________ 768.用频率采样法设计数字滤波器，对第二类型相位滤波器H(k)应具有的约束条件：幅值__________，相位_____________ 7.x(n)Acos(9．两序列h(n)=δ(n)+2δ(n-1)+3δ(n-2),x(n)=δ(n)+δ(n-1)，两者的线性卷积为y(n)，则y(2)_____ ________；若两者3点圆周卷积为y1(n)，则y1(0)＝__________________y1(2)＝__________________。三 计算题

1.有一个线性移不变的系统，其系统函数为：

3z112 H(z) z2

12(1z1)(12z1)21）用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

4．试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数：

H(s)=2其中抽样周期T=1s。

(s1)(s3)G

三、有一个线性移不变的因果系统，其系统函数为：

3z12 H(z) 11(1z)(12z1)21用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

七、用双线性变换设计一个三阶巴特沃思数字低通虑波器，采样频率为fs4kHz（即采样周期为T250s）,其3dB截止频率为fc1kHz。三阶模拟巴特沃思滤波器为：

Ha(s)

答案 12(sc1)2(sc)(s2c)3

一、选择题（10分，每题1分）

1.A 2.D 3.B 4.A 5.B 6.C 7.C 8.D 9.D 10.D

二、填空题（共25分 3、4、7、9每空2分；其余每空1分）1.栅栏效应 2.x(z)|z=wN-k x(k)＝X(ejw)|w＝2k 3.abz1cz2 4.8

N6144us 5.线性相位 频谱混迭、低通带通 6.2、5、－2w 7、14 9.HkHNk、k(11)10、5、4、5

N三计算题 1.（15分）

解1)H(z)1(1z1)(12z1)231z251z1z2231z2 ……………………………..2分

1时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分

3z1112………………………………..12分 H(z)1111(1z)(12z1)1z112z221h(n)()nu(n)2nu(n1)………………………………….15分

24.（10分）解： 当2zH(s)H(z)111………………1分

(1s)(s3)1ss3TT……………………3分

1eTZ1se3TZ10.318z1……………5分 10.418z10.018z22)H(z)H(s)|s21ZT1Z11221Z121Z1(1)(3)T1Z1T1Z1……8分

24z12z2…………………………… 10分 12152zz

三、（15）

1.解1)H(z)分 1(1z1)(12z1)231z251z1z2231z2 ……………………………..21时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分 2)当2zH(z)(111z)(12z1)231z211………………………………..12分 1112z11z21h(n)()nu(n)2nu(n1)………………………………….15分

2七、（12分）解：

wc2fcT0.5………………………………………3分

2wc2Ctan()………………………………………5分 T2THa(s)112(Ts)2(Ts)2(Ts)3……………………………8分 22221ZT1Z11H(z)Ha(s)|1211Z11Z1s2(1Z11Z1)(21Z11Z1)3

113z13z2z323z2

A

一、选择题（每题3分，共5题）

1、nj()36x(n)e，该序列是

。B.周期NnA.非周期序列 6 C.周期N6

D.周期N2

2、序列x(n)aA.3、对u(n1)，则X(Z)的收敛域为。

D.Za B.Za

C.Za Za

x(n)(0n7)和y(n)(0n19)分别作20点DFT，得X(k)和Y(k)，F(k)X(k)Y(k),k0,1,19，f(n)IDFT[F(k)],n0,1,19，n在 范围内时，f(n)是x(n)和y(n)的线性卷积。

A.0n7

B.7n19

C.12n19 D.0n19

4、x1(n)R10(n)，x2(n)R7(n)，用DFT计算二者的线性卷积，为使计算量尽可能的少，应使DFT的长度N满足。

A.N16 B.N16

C.N16

D.N16

5.已知序列Z变换的收敛域为｜z｜<1，则该序列为

。A.有限长序列

B.右边序列

C.左边序列

D.双边序列

二、填空题（每题3分，共5题）

1、对模拟信号（一维信号，是时间的函数）进行采样后，就是

信号，再进行幅度量化后就是

信号。

2、要想抽样后能够不失真的还原出原信号，则抽样频率必须

，这就是奈奎斯特抽样定理。

3、对两序列x(n)和y(n)，其线性相关定义为。

4、快速傅里叶变换（FFT）算法基本可分为两大类，分别是：

。

5、无限长单位冲激响应滤波器的基本结构有直接Ⅰ型，______ 和

四种。

an

三、x(n)nb

四、求 B n0

求该序列的Z变换、收敛域、零点和极点。（10分）

n1X(Z)11z112z1，1z2 的反变换。（8分）

一、单项选择题（本大题12分，每小题3分）

1、x(n)cos(0.125n)的基本周期是

。（A）0.125（B）0.25（C）8（D）16。

2、一个序列x(n)的离散傅里叶变换的变换定义为

。（A）（C）N1X(e)X(z)jnx(n)znx(n)enjn

（B）

X(k)x(n)ej2nk/Nn0N1n0

（D）

X(zk)x(n)AnWkn。

3、对于M点的有限长序列，频域采样不失真恢复时域序列的条件是频域采样点数N

。（A）不小于M

（B）必须大于M

（C）只能等于M

（D）必须小于M。

4、有界输入一有界输出的系统称之为。

（A）因果系统

（B）稳定系统

（C）可逆系统

（D）线性系统。

三、填空题（本大题10分，每小题2分）

1、在对连续信号进行频谱分析时，频谱分析范围受

速率的限制。

2、(d。

3、对于一个系统而言，如果对于任意时刻n0，系统在该时刻的响应仅取决于在时刻及其以前的输入，则称该系统为

系统。

4、对一个LSI系统而言，系统的输出等于输入信号与系统单位采样响应的线性。

5、假设时域采样频率为32kHz，现对输入序列的32个点进行DFT运算。此时，DFT输出的各点频率间隔为

Hz。

七、综合题（本大题20分）已知连续时间信号xa(t)cos(16000t)，用T1/6000对其采样。

（1）求最小采样频率；（2）图示其频谱特性；（3）分析其频谱是否有混叠。

C

一、单项选择题(本大题共10小题，每小题2分，共20分)1.在对连续信号均匀采样时，要从离散采样值不失真恢复原信号，则采样角频率Ωs与信号最高截止频率Ωc应满足关系（）

A.Ωs>2Ωc

B.Ωs>Ωc

C.Ωs<Ωc

D.Ωs<2Ωc

2.下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？（）A.y(n)=y(n-1)x(n)B.y(n)=x(n)/x(n+1)C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1)3.已知某序列Z变换的收敛域为5>|z|>3，则该序列为（）A.有限长序列

B.右边序列

C.左边序列 D.双边序列 4.实偶序列傅里叶变换是（）

A.实偶序列

B.实奇序列

C.虚偶序列

D.虚奇序列 5.已知x(n)=δ(n)，其N点的DFT［x(n)］=X(k)，则X(N-1)=（）A.N-1

B.1

C.0 的点数至少应取（）

A.M+N

B.M+N-1

C.M+N+1 7.下面说法中正确的是（）A.连续非周期信号的频谱为周期连续函数 B.连续周期信号的频谱为周期连续函数 C.离散非周期信号的频谱为周期连续函数 D.离散周期信号的频谱为周期连续函数

8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构？（）A.直接型

B.级联型

C.频率抽样型 D.并联型 9.下列关于FIR滤波器的说法中正确的是（）A.FIR滤波器容易设计成线性相位特性 B.FIR滤波器的脉冲响应长度是无限的 C.FIR滤波器的脉冲响应长度是确定的

D.对于相同的幅频特性要求，用FIR滤波器实现要比用IIR滤波器实现阶数低 10.下列关于冲激响应不变法的说法中错误的是（）A.数字频率与模拟频率之间呈线性关系

B.能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器 C.具有频率混叠效应

D.可以用于设计低通、高通和带阻滤波器

三、填空题(本大题共5小题，每空2分，共20分)。16.线性移不变系统是因果系统的充分必要条件是________。17.傅里叶变换的四种形式________，________，________和________。

18.使用DFT分析模拟信号的频谱时，可能出现的问题有________、栅栏效应和________。19.下图所示信号流图的系统函数为________。

D.2(M+N)

D.-N+1 6.设两有限长序列的长度分别是M与N，欲通过计算两者的圆周卷积来得到两者的线性卷积，则圆周卷积

20.对于N点（N＝2L）的按时间抽取的基2FFT算法，共需要作________次复数乘和________次复数加。

四、计算题

1z423．(10分)考虑一个具有系统函数H(z)1611z4161）求系统的零点和极点，并作出图表示； 2）画出系统的级联型结构图。的稳定系统。

24.(10分)有一用于频谱分析的FFT处理器，其抽样点数必须是2的整数次幂，假定没有采用任何特殊的数据处理措施，已知条件为：1）频率分辨率小于10Hz；2）信号最高频率小于4kHz。试确定以下参量： 1）最小记录长度tp； 2）最大抽样间隔T； 3）在一个记录中的最少点数N。

25．(10分)将双线性变换应用于模拟巴特沃兹滤波器Ha(s)11s/c，设计一个3dB截止频率cΩc）

D 3的一阶数字滤波器。（注：式中模拟巴特沃兹滤波器的3dB截止频率为

一、单项选择题(每小题3分，共24分)

1、在对连续信号均匀采样时，要从离散采样值不失真恢复原信号，则采样周期Ts与信号最高截止频率fh应满足关系

A.Ts>2/fh

B.Ts>1/fh

C.Ts<1/fh

D.Ts<1/(2fh)

2、下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？（）A.y(n)=x3(n)

B.y(n)=x(n)x(n+2)

C.y(n)=x(n)+2

D.y(n)=x(n2)

3、已知某序列z变换的收敛域为|z|<1，则该序列为（）。

A．有限长序列 列

4、设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取（）。

A．M＋N

B.M＋N－1

C.M＋N＋1

D.2(M＋N)

B.右边序列

C.左边序列

D.双边序

5、计算N=2L（L为整数）点的按时间抽取基-2FFT需要（）级蝶形运算。

A．L B.L/2

C.N

D.N/2 6.、因果FIR滤波器的系统函数H(z)的全部极点都在（）处。

A.z = 0

B.z = 1

C.z = j

D.z =∞

7、下列对IIR滤波器特点的论述中错误的是（）。

A．系统的单位冲激响应h(n)是无限长的B.结构必是递归型的 C.系统函数H(z)在有限z平面（0<|z|<∞）上有极点

D.肯定是稳定的8、线性相位FIR滤波器主要有以下四类(Ⅰ)h(n)偶对称，长度N为奇数

(Ⅱ)h(n)偶对称，长度N为偶数(Ⅲ)h(n)奇对称，长度N为奇数

(Ⅳ)h(n)奇对称，长度N为偶数 则其中不能用于设计高通滤波器的是（）。

A.Ⅰ、Ⅱ B.Ⅱ、Ⅲ

C.Ⅲ、Ⅳ

D.Ⅳ、Ⅰ

二、填空题（每题3分，共24分）

1、序列x(n)Asin(13n)的周期是。

32、序列R4(n)的Z变换为__

____，其收敛域为____

__。

3、对序列 力。

5、下图所示信号流图的系统函数为H(z)=_____

_____。

x(n)(nn0)，0n0N的N点的DFT为，0KN。

4、用DFT对连续信号进行频谱分析时，可能出现的问题有

、__、和DFT的分辨

6、有一模拟系统函数Ha(s)2，已知采样周期为T，采用脉冲响应不变法将其转换为数字系统函s3数H(z)是。

7、在利用窗函数法设计FIR滤波器时，一般希望窗函数能满足两项要求：①

；②

。但是，一般来说，以上两点很难同时满足。

8、IIR滤波器的有限字长效应与它的结构有关，结构的输出误差最小，结构输出误差其次，结构的输出误差最大。

五、用双线性变换法设计一个三阶巴特沃思数字低通滤波器，采样频率1.2kHz，截止频率为400Hz。要求⑴求该数字滤波器的系统函数，并画出其级联型结构；（归一化的三阶巴特沃思低通滤波器的模拟系统函数为

六、用矩形窗设计一线性相位低通FIR滤波器，设计要求：（1）若截止频率C、窗口长度N为已知，求该滤波器的单位抽样响应;（2）若C

E 1．序列1Ha(s)112s2s2s3）

（14分）

0.25，N＝33，x(n)anu(n)的Z变换为

，为/2时，信号的模拟角频率

x(n3)的Z变换是

。2．设采样频率fs1000Hz，则当和实际频率

f分别为

、。

3．N点序列x(n)的DFT表达式为

，其物理意义是

。4．序列x(n)和h(n)，长度分别为N和M（N>M），二者线性卷积的长度为

N点循环卷积中混叠的点有

个，循环卷积与线性卷积的关系是

5．全通系统的极零点分布特点是

三、分析计算题：（共 50分）

1．（15分）已知序列x(n){1,2,3,2,1}，n=0,1…,4(1)该序列是否可以作为线性相位FIR滤波器的单位脉冲响应？为什么？

(2)设序列x(n)的傅立叶变换用

X(ej)表示，不用求X(ej)，分别计算X(ej2j0)、X(ej)、X(ej)d、X(e)d。

(3)求x(n)与序列 y(n)R4(n)的线性卷积及7点循环卷积。

2．（15分）已知一因果系统的系统函数为

10.5z1H(z)3221z1z525试完成下列问题：

(1)系统是否稳定？为什么？

(2)求单位脉冲响应h(n)(3)写出差分方程；

(4)画出系统的极零图；(5)画出系统的所需存储器最少的实现结构。

sa(sa)2b2稳定，试用脉冲响应不变法将其转换成数字滤波器H(z)。3．（5分）已知模拟滤波器的传输函数Ha(s)：式中，a、b为常数，设Ha(s)因果

F

一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过（）即可完全不失真恢复原信号。

A.理想低通滤波器

B.理想高通滤波器

C.理想带通滤波器

D.理想带阻滤波器 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？（）A.y(n)=x3(n)

B.y(n)=x(n)x(n+2)

C.y(n)=x(n)+2

D.y(n)=x(n2)3.．设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取（）。A．M+N

B.M+N-1

C.M+N+1

D.2(M+N)4.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是（）。

A.N≥M

B.N≤M

C.N≤2M

D.N≥2M 5.直接计算N点DFT所需的复数乘法次数与（）成正比。

A.N

B.NC.N

3D.Nlog2N

6.下列各种滤波器的结构中哪种不是FIR滤波器的基本结构（）。

A.直接型

B.级联型

C.并联型

D.频率抽样型

7.第二种类型线性FIR滤波器的幅度响应H(w)特点()：

A 关于w0、、2偶对称

B 关于wC 关于w0、、2奇对称

奇对称

D关于w0、2奇对称 关于w偶对称 0、2偶对称 关于w8.适合带阻滤波器设计的是：（）A h(n)h(N1n)N为偶数

B h(n)h(N1n)N为奇数 C h(n)h(N1n)N为偶数

D h(n)h(N1n)N为奇数

9.以下对双线性变换的描述中不正确的是（）。A.双线性变换是一种非线性变换

B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内 D.以上说法都不对

10.关于窗函数设计法中错误的是：

A窗函数的截取长度增加，则主瓣宽度减小；

B窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的截取长度无关；

C为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加；

D窗函数法不能用于设计高通滤波器；

二、填空题(每空2分，共20分)

1.用DFT近似分析连续信号频谱时, _________效应是指DFT只能计算一些离散点上的频谱。

2.有限长序列X(z)与X（k）的关系

X（k）与X(ejw)的关系

3.下图所示信号流图的系统函数为：

4.如果通用计算机的速度为平均每次复数乘需要４μs，每次复数加需要1μs，则在此计算机上计算210点的基2FFT需要__________级蝶形运算,总的运算时间是__________μs。5.单位脉冲响应不变法优点, 缺点____________,适合_______________________滤波器设计 6.已知FIR滤波器H(z)12z15z2az3z4具有线性相位，则a＝______,冲激响应h（2）＝___,相位(w)___ 7.x(n)Acos(3n)的周期__________________ 768.用频率采样法设计数字滤波器，对第二类型相位滤波器H(k)应具有的约束条件：幅值__________，相位_____________

9．两序列h(n)=δ(n)+2δ(n-1)+3δ(n-2),x(n)=δ(n)+δ(n-1)，两者的线性卷积为y(n)，则y(2)_____ ________；若两者3点圆周卷积为y1(n)，则y1(0)＝__________________y1(2)＝__________________。三

计算题

1.有一个线性移不变的系统，其系统函数为：



H(z)(131z211z)(12z1)21 z2

21）用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

4．试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数：

H(s)=2其中抽样周期T=1s。

(s1)(s3)G

三、有一个线性移不变的因果系统，其系统函数为：



H(z)1(1z1)(12z1)231z2

1用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

七、用双线性变换设计一个三阶巴特沃思数字低通虑波器，采样频率为

fs4kHz（即采样周期为T250s）,其3dB截止频率为fc1kHz。三阶模拟巴特沃思滤波器为：

1Ha(s)

23sss12()2()()ccc

答案

二、选择题（10分，每题1分）

1.A 2.D 3.B 4.A 5.B 6.C 7.C 8.D 9.D 10.D

二、填空题（共25分3、4、7、9每空2分；其余每空1分）1.栅栏效应

2.x(z)|z=wN-k

x(k)＝X(ejw)|w＝

2k

3.aNbz1cz2 4.8

6144us

5.线性相位

频谱混迭、低通带通

6.2、5、－2w 7、14

9.1)10、5、、5 HkHNk、k(1N三计算题 1.（15分）

解1)H(z)1(1z1)(12z1)231z251z1z2231z2 ……………………………..2分

当2z1时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分

H(z)1(1z1)(12z1)21h(n)()nu(n)2nu(n1)………………………………….15分

24.（10分）解： 31z2111112z11z2………………………………..12分

111………………1分(1s)(s3)1ss3TTH(z)T13T1……………………3分 1eZseZ0.318z1……………5分 10.418z10.018z2H(s)2)H(z)H(s)|s21ZT1Z11221Z121Z1(1)(3)T1Z1T1Z1……8分

24z12z2……………………………

10分 152z1z

2三、（15）

1.解1)H(z)1(1z1)(12z1)231z251z1z2231z2 ……………………………..2分

2)当2z1时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分

H(z)1(1z1)(12z1)21h(n)()nu(n)2nu(n1)………………………………….15分

231z2111112z11z2………………………………..12分

七、（12分）解：

wc2fcT0.5………………………………………3分

………………………………………5分 C2w2tan(c)T2THa(s)112(Ts)2(Ts)2(Ts)3222s21ZT1Z11……………………………8分

H(z)Ha(s)|1211Z12(1Z1)(21Z1)31Z11Z1

113z13z2z323z2

1Z1

第二篇：数字信号处理期末试卷(含答案)2

数字信号处理期末试卷(含答案)

一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。

1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过()即可完全不失真恢复原信号。

A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？()A.y(n)=x3(n)B.y(n)=x(n)x(n+2)C.y(n)=x(n)+2

D.y(n)=x(n2)3.．设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取()。A．M+N B.M+N-1

C.M+N+1

D.2(M+N)4.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是()。

A.N≥M B.N≤M C.N≤2M D.N≥2M 5.直接计算N点DFT所需的复数乘法次数与()成正比。A.N B.N2 C.N3 D.Nlog2N 6.下列各种滤波器的结构中哪种不是FIR滤波器的基本结构()。A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR滤波器的幅度响应H(w)特点()： A 关于w0、、2偶对称

B 关于w0、、2奇对称

C 关于w0、2偶对称 关于w奇对称

D关于w0、2奇对称 关于w偶对称 8.适合带阻滤波器设计的是：（）A h(n)h(N1n)N为偶数 B h(n)h(N1n)N为奇数 C h(n)h(N1n)N为偶数 D h(n)h(N1n)N为奇数

9.以下对双线性变换的描述中不正确的是()。A.双线性变换是一种非线性变换

B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内 D.以上说法都不对

10.关于窗函数设计法中错误的是：

A窗函数的截取长度增加，则主瓣宽度减小；

B窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的截取长度无关； C为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加； D窗函数法不能用于设计高通滤波器；

二、填空题(每空2分，共20分)1.用DFT近似分析连续信号频谱时, _________效应是指DFT只能计算一些离散点上的频谱。

2.有限长序列X(z)与X（k）的关系 X（k）与X(ejw)的关系 3.下图所示信号流图的系统函数为：

4.如果通用计算机的速度为平均每次复数乘需要４μs，每次复数加需要1μs，则在此计算机上计算210点的基2FFT需要__________级蝶形运算,总的运算时间是__________μs。

5.单位脉冲响应不变法优点 , 缺点____________,适合_______________________滤波器设计

6.已知FIR滤波器H(z)12z15z2az3z4具有线性相位，则a＝______,冲激响应h（2）＝___,相位(w)___ 37.x(n)Acos(n)的周期__________________ 768.用频率采样法设计数字滤波器，对第二类型相位滤波器H(k)应具有的约束条件：幅值__________，相位_____________ 9．两序列h(n)=δ(n)+2δ(n-1)+3δ(n-2),x(n)=δ(n)+δ(n-1)，两者的线性卷积为y(n)，则y(2)_____ ________；若两者3点圆周卷积为y1(n)，则y1(0)＝__________________y1(2)＝__________________。三 计算题

1.有一个线性移不变的系统，其系统函数为：

3z112 H(z) z2

12(1z1)(12z1)21）用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

答案

一、选择题（10分，每题1分）

1.A 2.D 3.B 4.A 5.B 6.C 7.C 8.D 9.D 10.D

二、填空题（共25分 3、4、7、9每空2分；其余每空1分）

12k 3.abzcz 4.8 1.栅栏效应 2.x(z)|z=wN-k x(k)＝X(ejw)|w＝2N6144us 5.线性相位 频谱混迭、低通带通 6.2、5、－2w 7、14 9.HkHNk、k(11)10、5、4、5

N三计算题 1.（15分）

解1)H(z)1(1z1)(12z1)231z251z1z2231z2 ……………………………..2分

1时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分

3z1112………………………………..12分 H(z)111112z(1z)(12z1)1z1221h(n)()nu(n)2nu(n1)………………………………….15分 当2z

第三篇：数字信号处理期末试卷(含答案)1

数字信号处理期末试卷(含答案)

一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。

1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过(a)即可完全不失真恢复原信号。

A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？(d)A.y(n)=x3(n)B.y(n)=x(n)x(n+2)C.y(n)=x(n)+2

D.y(n)=x(n2)3.．设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取(b)。A．M+N B.M+N-1

C.M+N+1

D.2(M+N)4.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是(a)。

A.N≥M B.N≤M C.N≤2M D.N≥2M 5.直接计算N点DFT所需的复数乘法次数与(b)成正比。A.N B.N2 C.N3 D.Nlog2N 6.下列各种滤波器的结构中哪种不是FIR滤波器的基本结构(c)。A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR滤波器的幅度响应H(w)特点(c)： A 关于w0、、2偶对称

B 关于w0、、2奇对称

C 关于w0、2偶对称 关于w奇对称

D关于w0、2奇对称 关于w偶对称 8.适合带阻滤波器设计的是：（d）A h(n)h(N1n)N为偶数 B h(n)h(N1n)N为奇数 C h(n)h(N1n)N为偶数

D h(n)h(N1n)N为奇数

9.以下对双线性变换的描述中不正确的是(d)。A.双线性变换是一种非线性变换

B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内 D.以上说法都不对

10.关于窗函数设计法中错误的是：d A窗函数的截取长度增加，则主瓣宽度减小；

B窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的截取长度无关； C为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加； D窗函数法不能用于设计高通滤波器；

二、填空题(每空2分，共20分)1.用DFT近似分析连续信号频谱时, _________效应是指DFT只能计算一些离散点上的频谱。

2.有限长序列X(z)与X（k）的关系 X（k）与X(ejw)的关系 3.下图所示信号流图的系统函数为：

4.如果通用计算机的速度为平均每次复数乘需要４μs，每次复数加需要1μs，则在此计算机上计算210点的基2FFT需要__________级蝶形运算,总的运算时间是__________μs。

5.单位脉冲响应不变法优点 , 缺点____________,适合_______________________滤波器设计

6.已知FIR滤波器H(z)12z15z2az3z4具有线性相位，则a＝______,冲激响应h（2）＝___,相位(w)___ 3n)的周期__________________ 768.用频率采样法设计数字滤波器，对第二类型相位滤波器H(k)应具有的约束条件：幅值__________，相位_____________ 7.x(n)Acos（9．两序列h(n)=δ(n)+2δ(n-1)+3δ(n-2),x(n)=δ(n)+δ(n-1)，两者的线性卷积为y(n)，则y(2)_____ ________；若两者3点圆周卷积为y1(n)，则y1(0)＝__________________y1(2)＝__________________。三 计算题

1.有一个线性移不变的系统，其系统函数为：

3z112 H(z) z2

12(1z1)(12z1)21）用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

4．试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数：

H(s)=2其中抽样周期T=1s。

(s1)(s3)G

三、有一个线性移不变的因果系统，其系统函数为：

3z12 H(z) 11(1z)(12z1)21用直接型结构实现该系统

2)讨论系统稳定性，并求出相应的单位脉冲响应h(n)

七、用双线性变换设计一个三阶巴特沃思数字低通虑波器，采样频率为fs4kHz（即采样周期为T250s）,其3dB截止频率为fc1kHz。三阶模拟巴特沃思滤波器为：

Ha(s)

答案 12(sc1)2(sc)(s2c)3

一、选择题（10分，每题1分）

1.A 2.D 3.B 4.A 5.B 6.C 7.C 8.D 9.D 10.D

二、填空题（共25分 3、4、7、9每空2分；其余每空1分）

1.栅栏效应 2.x(z)|z=wN-k x(k)＝X(ejw)|w＝2k 3.abz1cz2 4.8

N6144us 5.线性相位 频谱混迭、低通带通 6.2、5、－2w 7、14 9.HkHNk、k(11)10、5、4、5

N三计算题 1.（15分）

解1)H(z)1(1z1)(12z1)231z251z1z2231z2 ……………………………..2分

1时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分

3z1112………………………………..12分 H(z)1111(1z)(12z1)1z112z221h(n)()nu(n)2nu(n1)………………………………….15分

24.（10分）解： 当2zH(s)H(z)111………………1分

(1s)(s3)1ss3TT……………………3分

1eTZ1se3TZ10.318z1……………5分 10.418z10.018z22)H(z)H(s)|s21ZT1Z11221Z121Z1(1)(3)T1Z1T1Z1……8分

24z12z2…………………………… 10分 12152zz

三、（15）

1.解1)H(z)分 1(1z1)(12z1)231z251z1z2231z2 ……………………………..21时： 2收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分 2)当2z

H(z)(111z)(12z1)231z211………………………………..12分 1112z11z21h(n)()nu(n)2nu(n1)………………………………….15分

2七、（12分）解：

wc2fcT0.5………………………………………3分

2wc2Ctan()………………………………………5分 T2THa(s)112(Ts)2(Ts)2(Ts)3……………………………8分 22221ZT1Z11H(z)Ha(s)|1211Z11Z1s2(1Z11Z1)(21Z11Z1)3

113z13z2z3 23z2

第四篇：数字信号处理课程总结(全)

数字信号处理课程总结

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

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

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

课程设计报告

课程名称： 数字信号处理 课题名称： 语音信号的处理与滤波

姓 名： 学 号： 院 系： 专业班级： 指导教师： 完成日期： 2013年7月2日

目录

第1部分 课程设计报告………………………………………3 一．设计目的……………………………………………3 二．设计内容……………………………………………3 三．设计原理……………………………………………3 四．具体实现……………………………………………5 1.录制一段声音…………………………………5 2.巴特沃斯滤波器的设计………………………8 3.将声音信号送入滤波器滤波…………………13 4.语音信号的回放………………………………19 5.男女语音信号的频谱分析……………………19 6.噪声的叠加和滤除……………………………22 五． 结果分析……………………………………………27 第2部分 课程设计总结………………………………28 一． 参考文献……………………………………………28

第1部分 课程设计报告

一．设计目的

综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。

二．设计内容

录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法和双线性变换法设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号；换一个与你性别相异的人录制同样一段语音内容，分析两段内容相同的语音信号频谱之间有什么特点；再录制一段同样长时间的背景噪声叠加到你的语音信号中，分析叠加前后信号频谱的变化，设计一个合适的滤波器，能够把该噪声滤除；

三．设计原理

1.在Matlab软件平台下，利用函数wavrecord()，wavwrite(),wavread(),wavplay()对语音信号进行录制，存储，读取，回放。

2.用y=fft(x)对采集的信号做快速傅立叶变换，并用[h1,w]=freqz(h)进行DTFT变换。

3.掌握FIR DF线性相位的概念，即线性相位对h(n)、H()及零点的约束，了解四种FIR DF的频响特点。

4.在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波。

5.抽样定理

连续信号经理想抽样后时域、频域发生的变化（理想抽样信号与连续信号频谱之间的关系）

理想抽样信号能否代表原始信号、如何不失真地还原信号即由离散信号恢复连续信号的条件（抽样定理）

理想采样过程描述： 时域描述：

ˆa(t)xa(t)T(t)xa(t)(tnT)xa(nT)(tnT)xnnT(t)频域描述：利用傅氏变换的性质，时域相乘频域卷积，若

n(tnT)ˆa(t)Xa(j)xXa(j)xa(t)T(j)T(t)

则有

ˆ(j)1X(j)(j)XaaT2121ˆXa(j)Xa(jjk)Xa(jjks)TkTTkˆ(j)与X(j)的关系：理想抽样信号的频谱是连续信号频谱的Xaa

周期延拓,重复周期为s(采样角频率)。如果：

X(j)Xa(j)a0s/2s/2即连续信号是带限的，且信号最高频率不超过抽样频率的二分之一，则可不失真恢复。

奈奎斯特采样定理：要使实信号采样后能够不失真还原，采样频率必须大于信号最高频率的两倍：s2h 或 fs2fh

四．具体实现

1.录制一段声音

1.1录制并分析

在MATLAB中用wavrecord、wavread、wavplay、wavwrite对声音进行录制、读取、回放、存储。

程序如下：

Fs=8000;%抽样频率 time=3;%录音时间 fprintf('按Enter键录音%ds',time);%文字提示 pause;%暂停命令 fprintf('录音中......');x=wavrecord(time*Fs,Fs,'double');%录制语音信号 fprintf('录音结束');%文字提示 fprintf('按Enter键回放录音');pause;%暂停命令

wavplay(x,Fs);%按任意键播放语音信号

wavwrite(x,Fs,'C:UsersacerDesktop数字信号sound.wav');%存储语音信号

N=length(x);%返回采样点数 df=fs/N;%采样间隔 n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 figure(2);subplot(2,1,1);plot(x);%录制信号的时域波形 title('原始信号的时域波形');%加标题 ylabel('幅值/A');%显示纵坐标的表示意义 grid;%加网格

y0=fft(x);%快速傅立叶变换 figure(2);subplot(2,1,2);plot(f,abs(y0(n1)));%原始信号的频谱图 title('原始信号的频谱图');%加标题 xlabel('频率w/pi');%显示横坐标表示的意义 ylabel('幅值 ');%显示纵坐标表示的意义 title('原始信号的频谱图');%加标题

grid;%加网格

图1.1 原始信号的时域与频谱图

1.2滤除无效点

针对实际发出声音落后录制动作半拍的现象，如何拔除对无效点的采样的问题: 出现这种现象的原因主要是录音开始时，人的反应慢了半拍，导致出现了一些无效点，而后而出现的无效的点，主要是已经没有声音的动作，先读取声音出来，将原始语音信号时域波形图画出来，根据己得到的信号，可以在第二次读取声音的后面设定采样点，取好有效点，画出滤除无效点后的语音信号时域波形图，对比可以看出。这样就可以解决这个问题。

x=wavread('C:UsersacerDesktop数字信号sound.wav', 7

[4000,24000]);%从4000点截取到24000结束 plot(x);%画出截取后的时域图形 title('截取后的声音时域图形');%标题 xlabel('频率');ylabel('振幅');grid;%画网格

图1.2 去除无效点

2.巴特沃斯滤波器的设计

2.1设计巴特沃思低通滤波器

MATLAB程序如下。滤波器图如图3.3所示。

%低通滤波

fp=1000;fs=1200;Fs=22050;rp=1;rs=100;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);[N,wc]=buttord(wap,was,rp,rs,'s');[B,A]=butter(N,wc,'s');[Bz,Az]=bilinear(B,A,Fs1);figure(1);[h,w]=freqz(Bz,Az,512,Fs1*22050);plot(w,abs(h));title('巴特沃斯低通滤波器');xlabel('频率（HZ）');ylabel('耗损（dB）');gridon;9

图2.1 巴特沃思低通滤波器

2.2设计巴特沃思高通滤波器

MATLAB程序如下。滤波器图如图3.5所示。%高通滤波

fp=4800;fs=5000;Fs=22050;rp=1;rs=100;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;T=1;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);10

[N,wc]=buttord(wap,was,rp,rs,'s');[B,A]=butter(N,wc,'high','s');[Bz,Az]=bilinear(B,A,Fs1);figure(1);[h,w]=freqz(Bz,Az,512,Fs1*22050);plot(w,abs(h));title('巴特沃斯高通滤波器');xlabel('频率（HZ）');ylabel('耗损（dB）');grid on;

图2.2巴特沃思高通滤波器

2.3设计巴特沃思带通滤波器

MATLAB程序如下。滤波器图如图3.7所示。%带通滤波

fp=[1200,3000];fs=[1000,3200];Fs=8000;rp=1;rs=100;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;T=1;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);[N,wc]=buttord(wap,was,rp,rs,'s');[B,A]=butter(N,wc,'s');[Bz,Az]=bilinear(B,A,Fs1);figure(4);[h,w]=freqz(Bz,Az,512,Fs1*1000);plot(w,abs(h));title('巴特沃斯带通滤波器');xlabel('频率（HZ）');ylabel('耗损（dB）');grid on;12

图2.3巴特沃思带通滤波器

3．将声音信号送入滤波器滤波

x=wavread('C:UsersacerDesktop数字信号sound.wav');%播放原始信号

wavplay(x,fs);%播放原始信号 N=length(x);%返回采样点数 df=fs/N;%采样间隔 n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 figure(4);subplot(4,2,1);plot(x);%录制信号的时域波形

title('原始信号的时域波形');%加标题 ylabel('幅值/A');%显示纵坐标的表示意义 grid;%加网格

y0=fft(x);%快速傅立叶变换 subplot(4,2,3);plot(f,abs(y0(n1)));%原始信号的频谱图 title('原始信号的频谱图');%加标题 xlabel('频率w/pi');%显示横坐标表示的意义 ylabel('幅值 ');%显示纵坐标表示的意义 title('原始信号的频谱图');%加标题 grid;%加网格

3．1低通滤波器滤波 fs=8000;beta=10.056;wc=2*pi*1000/fs;ws=2*pi*1200/fs;width=ws-wc;wn=(ws+wc)/2;n=ceil(12.8*pi /width);h=fir1(n,wn/pi,'band',kaiser(n+1,beta));[h1,w]=freqz(h);

ys=fftfilt(h,x);%信号送入滤波器滤波，ys为输出 fftwave=fft(ys);%将滤波后的语音信号进行快速傅立叶变换 figure(4);subplot(4,2,2);%在四行两列的第二个窗口显示图形 plot(ys);%信号的时域波形

title('低通滤波后信号的时域波形');%加标题 xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义 grid;%网格

subplot(4,2,4);%在四行两列的第四个窗口显示图形 plot(f, abs(fftwave(n1)));%绘制模值 xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义

title('低通滤波器滤波后信号的频谱图');%标题 grid;%加网格

wavplay(ys,8000);%播放滤波后信号

3．2高通滤波器滤波 fs=8000;beta=10.056;ws=2*5000/fs;wc=2*4800/fs;

width=ws-wc;wn=(ws+wc)/2;n=ceil(12.8*pi/width);h=fir1(n,wn/pi, 'high',kaiser(n+2,beta));[h1,w]=freqz(h);ys=fftfilt(h,x);%将信号送入高通滤波器滤波 subplot(4,2,5);%在四行两列的第五个窗口显示图形 plot(ys);%信号的时域波形 xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义 title('高通滤波后信号的时域波形');%标题 ylabel('幅值/A');%显示纵坐标的表示意义 grid;%网格

fftwave=fft(ys);%将滤波后的语音信号进行快速傅立叶变换 subplot(4,2,7);%在四行两列的第七个窗口显示图形 plot(f,abs(fftwave(n1)));%绘制模值 axis([0 1 0 50]);xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义

title('高通滤波器滤波后信号的频谱图');%标题 grid;%加网格

wavplay(ys,8000);%播放滤波后信号

3．3带通滤波器 fs=8000;beta=10.056;wc1=2*pi*1000/fs;wc2=2*pi*3200/fs;ws1=2*pi*1200/fs;ws2=2*pi*3000/fs;width=ws1-wc1;wn1=(ws1+wc1)/2;wn2=(ws2+wc2)/2;wn=[wn1 wn2];n=ceil(12.8/width*pi);h=fir1(n,wn/pi,'band',kaiser(n+1,beta));[h1,w]=freqz(h);ys1= fftfilt(h,x);%将信号送入高通滤波器滤波 figure(4);subplot(4,2,6);%在四行两列的第六个窗口显示图形 plot(ys1);%绘制后信号的时域的图形 title('带通滤波后信号的时域波形');%加标题 xlabel('频率w/pi');ylabel('幅值/A');%显示纵坐标表示的意义 grid;%网格

fftwave=fft(ys1);%对滤波后的信号进行快速傅立叶变换 subplot(4,2,8);%在四行两列的第八个窗口显示图形

plot(f, abs(fftwave(n1)));%绘制模值 axis([0 1 0 50]);xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义 title('带通滤波器滤波后信号的频谱图');%加标题 grid;%网格

wavplay(ys1,8000);%播放滤波后信号 图形如下：

原始信号的时域波形幅值/A0-1012x 10原始信号的频谱图34幅值/A1低通滤波后信号的时域波形0.50-0.5012频率w/pi3400.51频率w/pi高通滤波后信号的时域波形幅值/A0幅值/A0幅值/Ax 10高通滤波器滤波后信号的频谱图5012频率w/pi34幅值/A0.20-0.2幅值/A2001000x 10低通滤波器滤波后信号的频谱图200100000.51频率w/pi带通滤波后信号的时域波形0.50-0.501234频率w/pix 10带通滤波器滤波后信号的频谱图50幅值 00.5频率w/pi1000.5频率w/pi1

分析：三个滤波器滤波后的声音与原来的声音都发生了变化。其中低

通的滤波后与原来声音没有很大的变化，其它两个都又明显的变化

4.语音信号的回放

sound(xlow,Fs,bits);%在Matlab中，函数sound可以对声音进行回放，其调用格式: sound(xhigh, Fs,bits);%sound(x, Fs, bits);sound(xdaitong, Fs,bits);5.男女语音信号的频谱分析

5.1 录制一段异性的声音进行频谱分析

Fs=8000;%抽样频率 time=3;%录音时间 fprintf('按Enter键录音%ds',time);%文字提示 pause;%暂停命令 fprintf('录音中......');x=wavrecord(time*Fs,Fs,'double');%录制语音信号 fprintf('录音结束');%文字提示 fprintf('按Enter键回放录音');pause;%暂停命令 wavplay(x,Fs);%按任意键播放语音信号

wavwrite(x,Fs,'C:UsersacerDesktop数字信号sound2.wav');%存储语音信号

5.2 分析男女声音的频谱

x=wavread(' C:UsersacerDesktop数字信号sound2.wav ');%播放原始信号,解决落后半拍

wavplay(x,fs);%播放原始信号 N=length(x);%返回采样点数 df=fs/N;%采样间隔 n1=1:N/2;

f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 figure(1);subplot(2,2,1);plot(x);%录制信号的时域波形

title('原始女生信号的时域波形');%加标题 ylabel('幅值/A');%显示纵坐标的表示意义 grid;%加网格

y0=fft(x);%快速傅立叶变换 subplot(2,2,2);plot(f,abs(y0(n1)));%原始信号的频谱图 title('原始女生信号的频谱图');%加标题 xlabel('频率w/pi');%显示横坐标表示的意义 ylabel('幅值 ');%显示纵坐标表示的意义 grid;%加网格

[y,fs,bits]=wavread(' C:UsersacerDesktop数字信号sound.wav ');% 对语音信号进行采样

wavplay(y,fs);%播放原始信号 N=length(y);%返回采样点数 df=fs/N;%采样间隔 n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 subplot(2,2,3);plot(y);%录制信号的时域波形

title('原始男生信号的时域波形');%加标题 ylabel('幅值/A');%显示纵坐标的表示意义 grid;%加网格

y0=fft(y);%快速傅立叶变换

subplot(2,2,4);%在四行两列的第三个窗口显示图形 plot(f,abs(y0(n1)));%原始信号的频谱图 title('原始男生信号的频谱图');%加标题 xlabel('频率w/pi');%显示横坐标表示的意义 ylabel('幅值 ');%显示纵坐标表示的意义 grid;%加网格

5.3男女声音的频谱图

原始女生信号的时域波形0.50-0.5-1150100原始女生信号的频谱图幅值/A幅值 012345000x 10原始男生信号的时域波形0.50.5频率w/pi原始男生信号的频谱图1300200幅值/A0幅值 012x 1034100-0.5000.5频率w/pi1

图5.3男女声音信号波形与频谱对比

分析：就时域图看，男生的时域图中振幅比女生的高，对于频谱图女生的高频成分比较多

6.噪声的叠加和滤除

6.1录制一段背景噪声

Fs=8000;%抽样频率 time=3;%录音时间 fprintf('按Enter键录音%ds',time);%文字提示 pause;%暂停命令 fprintf('录音中......');x=wavrecord(time*Fs,Fs,'double');%录制语音信号

fprintf('录音结束');%文字提示 fprintf('按Enter键回放录音');pause;%暂停命令 wavplay(x,Fs);%按任意键播放语音信号 wavwrite(x,Fs,'C:UsersacerDesktop数字信号噪音.wav');%存储语音信号

6.2 对噪声进行频谱的分析

[x1,fs,bits]=wavread(' C:UsersacerDesktop数字信号噪音.wav ');%对语音信号进行采样

wavplay(x1,fs);%播放噪声信号 N=length(x1);%返回采样点数 df=fs/N;%采样间隔

n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 figure(5);subplot(3,2,1);plot(x1);%信号的时域波形 title('噪声信号的时域波形');grid;ylabel('幅值/A');y0=fft(x1);%快速傅立叶变换

subplot(3,2,2);plot(f,abs(y0(n1)));%噪声信号的频谱图 ylabel('幅值');title('噪声信号的频谱图');

6.3原始信号与噪音的叠加

fs=8000;[x,fs,bits]=wavread(' C:UsersacerDesktop数字信号sound.wav ');%对录入信号进行采样

[x1,fs,bits]=wavread(' C:UsersacerDesktop数字信号噪音.wav ');%对噪声信号进行采样

yy=x+x1;%将两个声音叠加

6.4叠加信号的频谱分析：

wavplay(yy,fs);%播放叠加后信号 N=length(yy);%返回采样点数 df=fs/N;%采样间隔 n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;%频带宽度 figure(5);subplot(3,2,3);plot(yy,'LineWidth',2);%信号的时域波形

title('叠加信号的时域波形');xlabel('时间/t');ylabel('幅值/A');grid;y0=fft(yy);%快速傅立叶变换 subplot(3,2,4);plot(f,abs(y0(n1)));%叠加信号的频谱图 title('叠加信号的频谱图');xlabel('频率w/pi');ylabel('幅值/db');grid;

6.5 设计一个合适的滤波器将噪声滤除 fs=18000;%采样频率 Wp=2*1000/fs;%通带截至频率 Ws=2*2000/fs;%阻带截至频率 Rp=1;%最大衰减 Rs=100;%最小衰减

[N,Wn]=buttord(Wp,Ws,Rp,Rs);%buttord函数（n为阶数，Wn为截至频率）

[num,den]=butter(N,Wn);%butter函数（num为分子系数den为分母系数）

[h,w]=freqz(num,den);%DTFT变换

ys=filter(num,den,yy);%信号送入滤波器滤波，ys为输出 fftwave=fft(ys);%将滤波后的语音信号进行快速傅立叶变换 figure(5);subplot(3,2,5);plot(ys);%信号的时域波形

title('低通滤波后信号的时域波形');%加标题 ylabel('幅值/A');%显示标表示的意义 grid;%网格 subplot(3,2,6);plot(f, abs(fftwave(n1)));%绘制模值 title('低通滤波器滤波后信号的频谱图');%标题 xlabel('频率w/pi');ylabel('幅值/A');%显示标表示的意义 grid;%加网格

wavplay(ys,8000);%播放滤波后信号 grid;图形如下：

噪声信号的时域波形1100噪声信号的频谱图幅值/A0-1幅值0123450000.5叠加信号的频谱图1x 10叠加信号的时域波形10-101时间/t2200幅值/db34幅值/A100000.5频率w/pi1x 10低通滤波后信号的时域波形0.5低通滤波器滤波后信号的频谱图200幅值/A0-0.5幅值/A012x 1034100000.5频率w/pi1

图6.1噪音的叠加与滤除前后频谱对比

7.结果分析

1.录制刚开始时，常会出现实际发出声音落后录制动作半拍，可在[x,fs,bits]=wavread('d:matlavworkwomamaaiwo.wav')加 窗[x,fs,bits]=wavread('d:matlavworkwomamaaiwo.wav'，[100 10000])，窗的长度可根据需要定义。

2.语音信号通过低通滤波器后，把高频滤除，声音变得比较低沉。当通过高通滤波器后，把低频滤除，声音变得比较就尖锐。通过带通滤波器后，声音比较适中。

3.通过观察男生和女生图像知：时域图的振幅大小与性别无关，只与说话人音量大小有关，音量越大，振幅越大。频率图中，女生高 27

频成分较多。

4.叠加噪声后，噪声与原信号明显区分,但通过低通滤波器后，噪声没有滤除，信号产生失真。原因可能为噪声与信号频率相近无法滤除。

第2部分 课程设计总结

通过本次课程设计，使我们对数字信号处理相关知识有了更深刻的理解，尤其是对各种滤波器的设计。在设计的过程中遇到了很多问题，刚刚开始时曾天真的认为只要把以前的程序改了参数就可以用了，可是问题没有我想象中的那么简单，单纯的搬程序是不能解决问题的。通过查阅资料和请教同学收获了很多以前不懂的理论知识。再利用所学的操作，发现所写的程序还是没有能够运行，通过不断地调试，运行，最终得出了需要的结果。整个过程中学到了很多新的知识，特别是对Matlab的使用终于有些了解。在以后的学习中还需要深入了解这方面的内容。在这次的课程设计中让我体会最深的是：知识来不得半点的马虎。也认识到自己的不足，以后要进一步学习。

八．参考文献

[1]数字信号处理教程(第三版)程佩青 清华大学出版社 [2]MATLAB信号处理 刘波 文忠 电子工业出版社 [3]MATLAB7.1及其在信号处理中的应用 王宏 清华大学出版社

[4]MATLAB基础与编程入门 张威 西安电子科技大学出版社

[5] 数字信号处理及其MATLAB实验 赵红怡 张常 化学工业出版社

[6]MATLAB信号处理详解 陈亚勇等 人民邮电出版社 [7] 数字信号处理

钱同惠 机械工业出版社 29