数字信号处理试题及参考答案

数字信号处理期末复习题

一、单项选择题(在每个小题的四个备选答案中选出一个正确答案，并将正确答案的号码写在题干后面

的括号内，每题1分，共20分)

1.要从抽样信号不失真恢复原连续信号，应满足以下条件的哪几条(①)。

(I)原信号为带限

(II)抽样频率大于两倍信号谱的最高频率

(III)抽样信号通过理想低通滤波器

①.I、II②.II、III

③.I、III④.I、II、III

2.在对连续信号均匀采样时，假设采样角频率为Qs,信号最高截止频率为Qc,则折叠频率为

(④)。

①Qs②.Qc

③.Qc/2④.Qs/2

3.假设一线性移不变系统当输入为x(n)=5(n)时输出为y(n)=R3(n),则当输入为u(n)-u(n-2)时输出

为(②)o

①.Rs(n)②.R2(n)

®.Ra(n)+Rs(n—1)④.Rz(n)—Rs(n—1)

4.序列Z变换的收敛域为IzI〉1,则该序列为(②)。

①.有限长序列②.右边序列

③.左边序列④.双边序列

5.离散系统的差分方程为y(n)=x(n)+ay(nT),则系统的频率响应(③)。

①当IaI〈1时，系统呈低通特性

②.当IaI>1时，系统呈低通特性

③.当0<a〈l时，系统呈低通特性

④.当-l<a<0时，系统呈低通特性

6.序歹!jx(n)=R5(n),其8点DFT记为X(k),k=0,1,…，7,则*(0)为(④)。

①.2②.3

③.45

7.以下关于FFT的说法中错误的选项是（①）。

①.FFT是一种新的变换

②.FFT是DFT的快速算法

③.FFT根本上可以分成时间抽取法和频率抽取法两类

④.基2FFT要求序列的点数为2,（其中L为整数）

8.以下结构中不属于FIR滤波器根本结构的是（③）。

①.横截型②.级联型

③.并联型④.频率抽样型

9.某FIR滤波器单位抽样响应h（n）的长度为（M+1）,则在以下不同特性的单位抽样响应中可以用来设计

线性相位滤波器的是（④）。

①.h[n]=一h[M—n]

②.hEn]=h[M+n]

③.h[n]=-h[M-n+1]

④.h[n]=h[M-n+1]

10.以下关于用冲激响应不变法设计HR滤波器的说法中错误的选项是（④）。

①.数字频率与模拟频率之间呈线性关系

②.能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器

③・容易出现频率混叠效应

④.可以用于设计高通和带阻滤波器

11.利用矩形窗函数法设计FIR滤波器时，在理想特性的不连续点附近形成的过滤带的宽度近似等于

（①兀

①.窗函数幅度函数的主瓣宽度

②.窗函数幅度函数的主瓣宽度的一半

③.窗函数幅度函数的第一个旁瓣宽度

窗函数幅度函数的第一个旁瓣宽度的一半

12.连续信号抽样序列在(①)上的Z变换等于其理想抽样信号的傅里叶变换。

①单位圆②.实轴

③.正虚轴④.负虚轴

13.一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包含(①)o

①单位圆②.原点

③.实轴虚轴

14.某FIR滤波器单位抽样响应h(n)的长度为(M+1),则在以下不同特性的单位抽样响应中可以用来设计

线性相位滤波器的是(①)。

①.h[n]=一h[M一n]

②.h[n]=h[M+n]

③.h[n]=-h[M-n+1]

h[n]=h[M—n+1]

15.序列x(n)=nR4(n),则其能量等于(③)。

①.5②.10

③.14④.20

16.以下单位冲激响应所代表的线性移不变系统中因果稳定的是(③)。

①.h(n)=u(n)②.h(n)=u(n+1)

③.h(n)=R4(n)h(n)=R4(n+1)

17.以下序列中z变换收敛域包括z=0的是(③)。

①.u(n)②.-u(n)

③.u(-n)u(n-l)

18.实序列的傅里叶变换必是(①)。

①.共轨对称函数②.二共软反对称函数

③.线性函数④.双线性函数

19.序列x(n)=6(n),10点的DFT[x(n)]=X(k)〔0WkW9),则X(5)=(1)。

0.10②.1③.0-10

20.欲借助FFT算法快速计算两有限长序列的线性卷积，则过程中要调用(③)次FFT算法。

①.1②.2③.3@.4

21.不考虑某些旋转因子的特殊性，一般一个基2FFT算法的蝶形运算所需的复数乘法及复数加法次数

分别为(①)。

①.1和2②.1和1

③.2和1④.2和2

22.因果FIR滤波器的系统函数H(z)的全部极点都在(①)处。

①.z=0②.z=1

③.z=j④.z=8

23.线性相位FIR滤波器主要有以下四类

①h(n)偶对称，长度N为奇数②.-h(n)偶对称，长度N为偶数

③h(n)奇对称，长度N为奇数④h(n)奇对称，长度N为偶数

则其中不能用于设计高通滤波器的是(③)。

24、序列u(〃)的Z变换及收敛域为(①

zz

①----，l<|z|W8②----,1<|z|<°°③l,0W|z|W8④1,0W|z|<8

z-1z-1

25、序列的Z变换及收敛域为1①

①告,"Z1〈8②告Jzl〈/

z—z—

44

③二j—41Z|〈8|Z|<|

1--Z1—z

44

26、序列x(〃)=(-)n的Z变换及收敛域为［③

2

z”1)

■〈J

②-----------------——

①二11

z”1)

-<|z|<2④----------------^―-<1z1<2

2(l-2)(Z—；)2

l-2z^1

27、假设X(z)=「一，—I，则才(z)的Z反变换为〔④)

z-1-22

①x(〃)=(1)"7(加1)—(L)zu(n—l)②x5)=(L)k'u(加1)—(工)用u(nT)

2222

③X5)=(L)，〃(〃一1)—(1)"%(加1)④x(〃)=(1)"%(〃—1)—(1)"明(〃)

2222

28、序列x®,力(血分别如下图，y(ri)=x(玲*h(ri),则y(4)为(③)

①？3②5士③3@5

22

1211：

1+-Z-l+Z--Z-

③〃(2)=——----④H(分=——-----—

1+二」2-21+-Z-'

342

30、下面信号流图表示的系统函数为(④)>(力

①〃(z)=-1+2-'+5z2—6z3②〃(z)=l+3z'—z2—6z3

③〃(z)=l—3z'+5zJ—6④〃(z)=l—z'—5z"+6z'1

31、假设x(〃)是长度为"的实序列，且DFT[x(〃)]=X(/,x(n)=x(炉〃)，则有(②)

①X5)=-XlN—吩②X(%=X〈N—b

③T(A)=-法(N-%④％(4)=-才(M")

32、对实信号进行谱分析，假设要求谱分辨率尸W50Hz,则最小记录时间。点应为(③)

①0.5S②0.05S③0.02S④0.2S

33、对实信号进行谱分析，假设信号最高频率为£=10KHz,则最大采样间隔备.应为(③)

①0.1X10飞②0.01*10飞③0.5X1(TS④0.05X10飞

34、对于淤8点的基IFFT运算，在进行位倒序后，地址单元4(4)中存放的是输入序列x(④中的哪一个

值(①〕

①x⑴②x⑵③*⑷④x(0)

35、x(n)=6(n),N点的DFT[x(n)]=X(k),则X⑸=()o

①.N②.1③.0-N

40、DFT[x(n)]=X(k),下面说法中正确的选项是(②)«

①.假设x(n)为实数偶对称函数,则X(k)为虚数奇对称函数

②.假设x(n)为实数奇对称函数,则X(k)为虚数奇对称函数

③.假设x(n)为虚数偶对称函数,则X(k)为虚数奇对称函数

④.假设x(n)为虚数奇对称函数,则X(k)为虚数奇对称函数

36、如下图的运算流图符号是()基2FFT算法的蝶形运算流图符号。

①.按频率抽取②.按时间抽取

③.两者都是④.两者都不是

37、直接计算N点DFT所需的复数乘法次数与(②)成正比。

①.N②.N2

③.N3④.Nlog2N

38、以下各种滤波器的结构中哪种不是IIR滤波器的根本结构(④)。

①.直接型②.级联型

③.并联型④.频率抽样型

39、以下对双线性变换的描述中正确的选项是(②)。

A.双线性变换是一种线性变换

②.双线性变换可以用来进行数字频率与模拟频率间的变换

③.双线性变换是一种分段线性变换

④.以上说法都不对

40、假设模拟滤波器〃(s)=Y-,采样周期Z1S，则利用双线性变换法，将〃(S转换成数字滤

s2+2s

波器〃(力应为[②)

1z~+2z+11z2+2z+11z"+2z+11z~+2z+1

①-----5------②-----------③-----------④-----------

4z2-z42z~-2z42z~+2z4z~+z

41、抽样频率确定时，DFT的频率分辨力取决于(④)

①量化误差②信号带宽③抽样间隔④抽样点数

42、如果一线性移不变系统的收敛域为一半径小于1的圆的外部，则该系统为〔②）

①因果稳定系统②因果非稳定系统③稳定非因果系统④因果稳定系统

二、多项选择题（在每题的五个备选答案中选出而二至五个正确答案，并将其号码分别写在题干后的括

号内，未选全或有选错的，该题无分。每题1分，共12分）

1、以下系统中是因果的有〔①②）

①7[x（A）]=g（A）X（〃）②T[x（A）]=左）

k=%

③T[x（77）]=e"*④T\_x（n）]-ax{n}+b

2、以下说法中正确的有1①②）

①因果序列的z变换收敛域为办〈1z8,其中a为收敛域最小半径

②双边序列的Z变换收敛域为圆环

③左边序列的Z变换收敛域为此z1<8,其中醇为收敛域最小半径

④右边序列的Z变换收敛域为0〈|z|<R：,其中痣为收敛域最大半径

⑤有限长序列的Z变换总是收敛的

3、两序列卷积运算包括的步骤有（①③④⑤）

①翻褶②取模③平移④相乘⑤相加

4、以下说法正确的有（②③）

①FIR系统只能采用非递归结构的电路

②IIR体统只能采用递归结构的电路

③FIR系统可以采用递归或非递归结构的电路

④IIR系统可以采用递归或非递归结构的电路

⑤FIR系统的〃（z）在有限Z平面上无限点

6、用DFT进行谱分析时，截断后序列的频谱F3”）与原序列频谱才（4"）的差异对谱分析的影响主要表

现在（①②④⑤）

①频谱混叠②泄露③衰减④谱间干扰⑤栅栏效应

7、实现FIR滤波器的根本网络结构主要有1②③④⑤）

①并联型②级联型③直接型

④线性相位有限脉冲响应系统网络结构⑤频率采样型

8、利用数字域频率变换设计数字高通滤波器的步骤有1①②③⑤）

①将数字高通技术指标转换成模拟高通技术要求

②将模拟高通技术指标转换成模拟低通技术要求

③设计模拟低通滤波器

④将模拟低通转换成数字低通

⑤利用数字域频率变换将数字低通滤波器转换成数字高通滤波器

9、以下说法中不正确的有〔①④）

①在相同技术指标下，HR滤波器可用比FIR滤波器较少的阶数

②设计微分器或积分器等主要用HR滤波器

③FIR、IIR滤波器都可用快速傅立叶变换算法

④FIR滤波器可以得到严格的线性相位

⑤对图像信号处理，采用HR滤波器较好

10、由传输函数〃（力确定状态方程和输出方程的根本方法有（）

①直接法②巢式法③局部分式法④级联法⑤观察法

三、说明题（认为正确的，在题干后的括号内打“;认为错的打“X",并说明理由，否则该题无

分。）

1、正弦序列x（A）=/sin（A0o+。）为一周期序列（X）

2、实序列的傅立叶变换的幅度是。的奇函数[V）

3、假设信号持续时间无限长，则信号的频谱无限宽[V）

4.线性系统必然是移不变系统。（X）

5.非零周期序列的Z变换不存在。（V）

6.按时间抽取的基2FFT算法的运算量等于按频率抽取的基2FFT算法。（V）

7.通常FIR滤波器具有递归型结构。（X）

8.双线性变换法是非线性变换，所以用它设计HR滤波器不能克服频率混叠效应。（X）

9、H（z）唯一地对应一个系统冲激响应h（n）（V）

10、定点制比浮点制运算速度慢(X)

11、移不变系统必然是线性系统。X)

12、当输入序列不同时，线性移不变系统的单位抽样响应也不同。（X）

13、离散时间系统的滤波特性可以由其幅度频率特性直接看出。（X）

14、因果稳定系统的系统函数的极点必然在单位圆内。（V）

15、与FIR滤波器相似，IIR滤波器的也可以方便地实现线性相位。（X）

16.非零周期序列的能量为无穷大。（X）

17.序列的傅里叶变换就是序列z变换在单位圆上的取值。（X）

18.离散傅里叶变换具有隐含周期性。（X）

19.FIR滤波器必是稳定的。（X）

20.当输入序列不同时，线性移不变系统的单位抽样响应也不同。（X）

21.离散时间系统的滤波特性可以由其幅度频率特性直接看出。（X）

22.用窗函数法设计FIR低通滤波器时，可以通过增加截取长度N来任意减小阻带衰减。（X）

四、填空题

1.数字信号处理的主要对象是数字信号,采用数值运算的方法到达处理的目的;

其实现方法主要有硬件实现和软件实现O

2

=00

2.序列x（n）的能量定义为。

n=-<x>

3.对正弦信号xa=sin314t进行采样，采样频率为fs=200Hz,则所得到的采样序列x（n）=sinL〃

2

4.我们可以从三个角度用三种表示方法描述一个线性时不变离散时间系统，它们是差分方程、

系统函数和单位脉冲响应。

5.线性时不变系统是因果系统的充分必要条件是n<0,h（n）=0。

6.线性时不变系统是稳定系统的充分必要条件是£口（〃）|<8。

_00

7.线性常系数差分方程的求解方法有经典法、递推法和变域法

8.以下图是模拟信号数字处理框图，填写出各个框内的处理方法:

9.设两个有限长序列的长度分别为N和M,则他们的线性卷积的结果序列的长度为N+M-1。

10.一个长序列和一个短序列卷积时，有重叠相加法和重叠保存法两种分段卷积法。

12数字频率只有相对意义，因为它是实际频率对采样频率的归一化频率。

13.数字频率3是模拟频率Q对采样频率的归一化化,与数字频率”相对应的模拟频

率是模拟折叠频率Fs/2o

14.正弦序列sin(n3)不一定是周期序列，比方3取2不/。为无理数时就不是。

15.从满足采样定理的样值信号中可以不失真地恢复出原模拟信号。采用的方法，从时域角度看是:_时

域采样定理；从频域角度看是：频域采样定理

16.如果通用计算机的速度为平均每次复数乘需要4us,每次复数加需要1NS,则在此计算机上计算

2,点的基2FFT需要」^级蝶形运算，总的运算时间是」0x(N/2)x(N/2x4+Nxl)口s。

17.在用DFT近似分析连续信号的频谱时，栅栏效应是指DFT只能计算一些离散点上的频谱。

18.在FIR滤波器的窗函数设计法中，常用的窗函数有矩形和哈明等等。

19.试写出6(n)和6(t)的两点区别：

⑴6(n)是离散的序列而6⑴是连续的函数;

⑵6(t)在t=0时取值无穷大，而6[n)在n=0时取值为1。

21.判定某系统为因果系统的充要条件是：时域满足条件n〈0时，h(n)=0,等效于在频域满足条

件收敛域一定包含无穷点o

23.研究一个周期序列的频域特性，应该用DFS变换。

24.脉冲响应不变法的根本思路是：将s平面上的传递函数E(s)转换为z平面上的系统函数H(z)。

25.要获得线性相位的FIR数字滤波器，其单位脉冲响应h(n)必须满足条件：[1)奇对称；[2)

偶对称。

27,借助模拟滤波器的H(s)设计一个HR高通数字滤波器，如果没有强调特殊要求的话，宜选择采

用双线性变换法法。

28,周期序列之所以不能进行Z变换，是因为无限长序列。

N-1

30,某DFT的表达式是X[1)=\x(公阅'，则变换后数字频域上相邻两个频率样点之间的间隔是

G=0

2出N。

35.FIR系统成为线性相位的充要条件是。

37.用窗口法设计出一个FIR低通滤波器后，发现它过度带太宽，这样情况下宜采用的修改措施是

40.某线性移不变系统当输入x(n)=5(n-1)时输出y(n)=6(n-2)+8(n-3),则该系统的单位冲激

响应h(n)=3(nT)+6(n-2)。

41.序列x(n)=cos(3nn)的周期等于2/3□

43.基2FFT算法计算N=2L(L为整数)点DFT需L级蝶形,每级由N/2个蝶形运

算组成。

44.以下图所示信号流图的系统函数为H(z)=o

45.线性系统同时满足和两个性质。

46.以下图所示信号流图的系统函数为。

y(")

47.有限长单位冲激响应(FIR)滤波器的主要设计方7

48.将模拟滤波器映射成数字滤波器主要有__典

50.FFT的根本运算单元称为蝶形运算。

53.N点FFT的运算量大约是N/2。

54.设计一个带阻滤波器，宜用第第一类类FIR数字滤波器。

57.描述一个离散系统的方法，时域有差分方程Z域有系统函数欲求系统瞬态

解，可用差分方程。

58.某系统函数在单位圆外有极点但它却是稳定的，则该系统一定是因果系统。

62.在序列为无限长的情况下,序列傅氏变换存在，但其DFT不存在。

65.设计数字滤波器的方法之一是先设计模拟滤波器，然后通过模拟S域(拉氏变换域)到数字Z域的

变换，将模拟滤波器变换成数字滤波器，其中常用的双线形变换的关系式是—o

66.一个因果数字系统，如果系统的极点位于Z平面的范围，则该系统是稳定的。

67.模拟信号是指其时域波形连续的信号：

数字信号是指其时域波形离散的信号,时间离散和取值离散的信号被称为离散时间信号。

69.第二类线形相位FIR滤波器［h(n)偶对称，N为偶数)一定不能用做高通和带阻特性的滤波系统,

第三类线形相位FIR滤波器(h(n)奇对称，N为奇数)一定不能用做低通、高通和带阻特性的滤波

系统。

四、分析、作图与计算题

L判断以下系统是否为：门)稳定系统；(2)因果系统；[3)线性系统[4)移不变系统。

(1)T[x(n)]=ax(n)+b

(2)T[x(n)]=g(n)x(n)

(3)T[x(n)]=ex(n)

(4)T[x(n)]=x(n)sin(^-n+—)

97

2.判断以下序列是否周期序列，试确定其周期

(1)x(n)=Acosf^n-^)

(2)x(n)=e

3.设有一系统，其输入输出关系有以下差分方程确定

设系统是因果系统

m求系统的单位脉冲响应；

[2)系统的系统函数；

13)系统的频率函数并画出系统的幅频特性曲线。

4.求以下序列的Z变换

(1)x(n)=a|n1,|a|<1

(2)x(n)=(^-)Hw(n-l)

(3)x(n)=sin(won)

l-2

5.X(z)=———a——-,|a|<1,求其逆变换x(

(l-az^-az1)

6.试求以下序列的N点DFT〔闭和形式表达式)

(1)x(n)=acos(won)RN(n);

(2)x(n)=6(n-1)+8(n-2)

7、画出8点按时间抽取的基2FFT算法的运算流图。

8.写出16点基2FFT算法中位序颠倒的序列号。

9.画出8点按频率抽取的基2FFT算法的运算流图。

10.某线性移不变系统的单位抽样响应为：

h(n)=26(n)+6(n-l)+6(n-3)+26(n-4)

求其系统函数，并画出该系统的横截型结构〔要求用的乘法器个数最少)，该滤波器是否具有线性相位

特性，为什么？

11.用双线性变换法设计无限长单位冲激响应(IIR)数字低通滤波器，要求通带截止频率3c=0.5Jrrad,

通带衰减61不大于3dB,阻带截止频率3st=0.75mrad,阻带衰减82不小于20dB。以巴特沃思

(Butterworth)模拟低通滤波器为原型，采样间隔T=2s,写出设计步骤。

311

12.某系统的差分方程为y(n)-1)+—y(〃―2)=x(〃)+—x(〃—1)

483

[1)该系统的系统函数；

12)试用典型范型(直接H型)，一阶节的级联，一阶节的并联实现此方程。

13.设滤波器的差分方程为y(〃)=x(n)+x(hl)+』y(kl)+』y(h2)

34

①求系统的频率响应；

②系统函数；

③试用典型范型〔直接n型)及一阶节的级联，及一阶节的并联实现此方程。

五、实验题

(一)对模拟周期信号进行谱分析

11)对于周期序列，如果周期不知道，如何用FFT进行谱分析？

[2)如何选择FFT的变换区间？(包括非周期信号和周期信号〕

[3)当N=8时，/(〃)和/(〃)的幅频特性会相同吗？为什么？N=16呢？

实验程序清单

%第10章实验3程序exp3.m

%用FFT对信号作频谱分析

clearall;closeall

%实验内容(1)===================================================

xln=[ones(l,4)]；%产生序列向量xl(n)=R4(n)

M=8;xa=l:(M/2);xb=(M/2):-l:l;x2n=[xa,xb];%产生长度为8的三角波序列x2(n)

x3n=[xb,xa];

Xlk8=fft(xln,8);%计算xln的8点DFT

Xlkl6=fft(xln,16);%计算xln的16点DFT

X2k8=fft(x2n,8);%计算xln的8点DFT

X2kl6=fft(x2n,16);%计算xln的16点DFT

X3k8=fft(x3n,8);%计算xln的8点DFT

X3kl6=fft(x3n,16);%计算xln的16点DFT

%以下绘制幅频特性曲线

subplot(2,2,l);mstem(Xlk8);%绘制8点DFT的幅频特性图

title('(la)8点DFT[x_l(n)]');xlabel('6>/n);ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X1k8))])

subplot(2,2,3)；mstem(X1k16);%绘制16点DFT的幅频特性图

title('(lb)16点DFT[x_l(n)T);xlabel('3/n)ylabelC幅度为

axis([0,2,0,1.2*max(abs(X1k16))])

figure(2)

subplot(2,2,1);mstem(X2k8);%绘制8点DFT的幅频特性图

title('(2a)8点DFT[x_2(n)]');xlabel('«/n);ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X2k8))])

subplot(2,2,2);mstem(X2kl6);%绘制16点DFT的幅频特性图

title('(2b)16点DFT[x_2(n)]');xlabel('/u');ylabel('幅度');

axis([0,2,0,1.2*max(abs(X2kl6))])

subplot(2,2,3)；mstem(X3k8);%绘制8点DFT的幅频特性图

title('(3a)8点DFT[x_3(n)]');xlabel(,<o/n);ylabel('幅度');

axis([0,2,0,1.2*max(abs(X3k8))])

subplot(2,2,4)；mstem(X3kl6);%绘制16点DFT的幅频特性图

title('(3b)16点DFT[x_3(n)]');xlabel('3/n)ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X3kl6))])

%实验内容⑵周期序列谱分析==================================

N=8;n=0:N-l;%FFT的变换区间N=8

x4n=cos(pi*n/4);

x5n=cos(pi*n/4)+cos(pi*n/8);

X4k8=fft(x4n);%计算x4n的8点DFT

X5k8=fft(x5n);%计算x5n的8点DFT

N=16;n=0:N-l;%FFT的变换区间N=16

x4n=cos(pi*n/4);

x5n=cos(pi*n/4)+cos(pi*n/8);

X4kl6=fft(x4n);%计算x4n的16点DFT

X5kl6=ffi(x5n);%计算x5n的16点DFT

figure(3)

subplot(2,2,l);mstem(X4k8);%绘制8点DFT的幅频特性图

title('(4a)8点DFT[x_4(n)]');xlabel('o/n);ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X4k8))])

subplot(2,2,3)；mstem(X4k16);%绘制16点DFT的幅频特性图

title('(4b)16点DFT[x_4(n)]');xlabel('w/兀');ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X4kl6))])

subplot(2,2,2);mstem(X5k8);%绘制8点DFT的幅频特性图

title('(5a)8点DFT[x_5(n)]');xlabel('s/JT);ylabel('幅度')；

axis([0,2,0,1.2*max(abs(X5k8))])

subplot(2,2,4);mstem(X5k16);%绘制16点DFT的幅频特性图

title('(5b)16点DFT[x_5(n)T);xlabel('3/n)；ylabel('幅度');

axis([0,2,0,1.2*max(abs(X5kl6))])

%实验内容⑶模拟周期信号谱分析===============================

figure⑷

Fs=64;T=l/Fs;

N=16;n=0:N-l;%FFT的变换区间N=16

x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);%对x6(t)16点采样

X6kl6=fft(x6nT);%计算x6nT的16点DFT

X6kl6=fftshift(X6kl6);%将零频率移到频谱中心

Tp=N*T;F=l/Tp;%频率分辨率F

k=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率(以零频率为中心)

subplot(3,l,l);stem(fk,abs(X6kl6)；.');boxon%绘制8点DFT的幅频特性图

title(〈6a)16点|DFT[x_6(nT)V);xlabelCf(Hz));ylabelC幅度)

axis([-N*F/2-1,N*F/2-1,0,1,2*max(abs(X6kl6))])

N=32;n=0:N-1;%FFT的变换区间N=16

x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);%对x6(t)32点采样

X6k32=fft(x6nT);%计算x6nT的32点DFT

X6k32=fftshift(X6k32);%将零频率移到频谱中心

Tp=N*T;F=l/Tp;%频率分辨率F

k=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率(以零频率为中心)

subplot(3,l,2);stem(fk,abs(X6k32),'.');boxon%绘制8点DFT的幅频特性图

title('(6b)32点|DFT[x_6(nT)][');xlabel('f(Hz)');ylabel('幅度')；

axis([-N*F/2-1,N*F/2-l,0,1,2*max(abs(X6k32))])

N=64;n=0:N-l;%FFT的变换区间N=16

x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);%对x6(t)64点采样

X6k64=fft(x6nT);%计算x6nT的64点DFT

X6k64=fftshift(X6k64);%将零频率移到频谱中心

Tp=N*T;F=l/Tp;%频率分辨率F

k=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率(以零频率为中心)

subplot(3,l,3)；stem(fk,abs(X6k64)；.');boxon%绘制8点DFT的幅频特性图

title('(6a)64点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度')；

axis([-N*F/2-1,N*F/2-l,0,1,2*max(abs(X6k64))])

实验程序运行结果

实验3程序exp3.m运行结果如图10.3.1所示。

(二)、IIR数字滤波器设计及软件实现

(1)信号产生函数mstg清单

functionst=mstg

%产生信号序列向量St,并显示st的时域波形和频谱

%st=mstg返回三路调幅信号相加形成的混合信号，长度N=1600

N=1600%N为信号st的长度。

Fs=10000;T=l/Fs;Tp=N*T;%采样频率Fs=10kHz,Tp为采样时间

t=O:T:(N-l)*T;k=O:N-l;f=k/Tp;

fcl=Fs/10;%第1路调幅信号的载波频率fcl=1000Hz,

fml=fcl/10;%第1路调幅信号的调制信号频率fml=100Hz

fc2=Fs/20;%第2路调幅信号的载波频率fc2=500Hz

fm2=fc2/10;%第2路调幅信号的调制信号频率fm2=50Hz

fc3=Fs/40;%第3路调幅信号的载波频率fc3=250Hz,

fm3=fc3/10;%第3路调幅信号的调制信号频率fm3=25Hz

xtl=cos(2*pi*fml*t).*cos(2*pi*fcl*t);%产生第1路调幅信号

xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t);%产生第2路调幅信号

xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t);%产生第3路调幅信号

st=xtl+xt2+xt3;%三路调幅信号相加

fxt=fft(st,N);%计算信号st的频谱

%====以下为绘图局部，绘制st的时域波形和幅频特性曲线=================

subplot(3,1,1)

plot(t,st);grid;xlabel('t/s');ylabel('s(t)');

axis([0,Tp/8,min(st),max(st)]);titleC(a)s(t)的波形')

subplot⑶1,2)

stem(f,abs(fxt)/max(abs(fxt)grid;title('(b)s(t)的频谱')

axis([0,Fs/5,0,1.2]);

xlabel(*f/Hz,);ylabelC幅度')

(2)思考题

11)请阅读信号产生函数mstg,确定三路调幅信号的载波频率和调制信号频率。

12)信号产生函数mstg中采样点数N=800,对st进行N点FFT可以得到6根理想谱线。如果取

N=1000,可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg中采样点数N的值，观察频谱

图验证您的判断是否正确。

13)修改信号产生函数mstg,给每路调幅信号参加载波成分，产生调幅[AM)信号，重复本实验,

观察AM信号与抑制载波调幅信号的时域波形及其频谱的差异。

提示：AM信号表示式：s(t)=[1+cos(2〃4/)]cos(2万力/)。

(3)滤波器参数及实验程序清单

1、滤波器参数选取

观察图10.4.1可知，三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz。带宽〔也可以由信

号产生函数mstg清单看出)分别为50Hz、100Hz、200Hz»所以，别离混合信号st中三路抑制载波单频

调幅信号的三个滤波器(低通滤波器、带通滤波器、高通滤波器)的指标参数选取如下：

对载波频率为250Hz的条幅信号，可以用低通滤波器别离，其指标为

带截止频率/=280HZ,通带最大衰减=O.ldBdB；

阻带截止频率fs-450HZ,阻带最小衰减々s=60dBdB,

对载波频率为500Hz的条幅信号，可以用带通滤波器别离，其指标为

带截止频率力j=440Hz,于—560HZ,通带最大衰减戊夕—0.IdBJB；

阻带截止频率//=275Hz,fsu-900Hz,Hz,阻带最小衰减=60dBdB,

对载波频率为1000Hz的条幅信号，可以用高通滤波器别离，其指标为

带截止频率力=890HZ,通带最大衰减tZp=0.1dBdB；

阻带截止频率£=550Hz,阻带最小衰减。$=60dBdB,

说明：[1)为了使滤波器阶数尽可能低，每个滤波器的边界频率选择原则是尽量使滤波器过渡带宽

尽可能宽。

[2)与信号产生函数mstg相同，采样频率Fs=10kHz。

(3)为了滤波器阶数最低，选用椭圆滤波器。

按照图10.4.2所示的程序框图编写的实验程序为exp4.mo

2、实验程序清单

%实验4程序exp4.m

%IIR数字滤波器设计及软件实现

clearall;closeall

Fs=10000;T=l/Fs;%采样频率

%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st

st=mstg;

%低通滤波器设计与实现=====================

fp=280;fs-450;

wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;%DF指标〔低通滤波器的通、阻带边界频)

[N,wp]=ellipord(wp,ws,rp,rs);%调用ellipord计算椭圆DF阶数N和通带截止频率wp

[B,A]=ellip(N,rp,rs,wp);%调用ellip计算椭圆带通DF系统函数系数向量B和A

ylt=fliter(B,A,st);%滤波器软件实现

%低通滤波器设计与实现绘图局部

figure(2);subplot(3,1,1);

myplot(B,A);%调用绘图函数myplot绘制损耗函数曲线

yt=,y」(t)';

subplot(3,1,2);tplot(ylt,T,yt);%调用绘图函数tplot绘制滤波器输出波形

%带通滤波器设计与实现==========================

fpl=440;fpu=560;fsl=275;fsu=900;

wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=O.l;rs=60;

[N,wp]=ellipord(wp,ws,rp,rs);%调用ellipord计算椭圆DF阶数N和通带截止频率wp

[B,A]=ellip(N,rp,rs,wp);%调用ellip计算椭圆带通DF系统函数系数向量B和A

y2t=filter(B,A,st);%滤波器软件实现

%带通滤波器设计与实现绘图局部(省略)

%高通滤波器设计与实现=========================

fp=890;fs=600;

wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;%DF指标(低通滤波器的通、阻带边界频)

[N,wp]=ellipord(wp,ws,rp,rs);%调用ellipord计算椭圆DF阶数N和通带截止频率wp

[B,A]=ellip(N,rp,rs,wp,5high));%调用ellip计算椭圆带通DF系统函数系数向量B和A

y3t=filter(B,A,st);%滤波器软件实现

%上下通滤波器设计与实现绘图局部(省略)

3实验程序运行结果

实验4程序exp4.m运行结果如图104.2所示。由图可见，三个别离滤波器指标参数选取正确，算耗

函数曲线到达所给指标。别离出的三路信号yl(n),y2(n)和y3(n)的波形是抑制载波的单频调幅波。

(a)低通滤波器损耗函数及其别离出的调幅信号y“t)

(b)带通滤波器损耗函数及其别离出的调幅信号y2(t)

(c)高通滤波器损耗函数及其别离出的调幅信号ys(t)

图104.实验4程序exp4.m运行结果

4简要答复思考题

思考题(1)已经在10.4.2节解答。思考题(3)很简单，请读者按照该