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

1、数字信号处理期末复习题一、单项选择题（在每个小题的四个备选答案中选出一个正确答案，并将正确答案的号码写在题干后面的括号内，每小题1分，共20分）1.要从抽样信号不失真恢复原连续信号，应满足下列条件的哪几条( )。()原信号为带限()抽样频率大于两倍信号谱的最高频率()抽样信号通过理想低通滤波器 .、 .、.、 .、2.在对连续信号均匀采样时，若采样角频率为s，信号最高截止频率为c，则折叠频率为(        )。 s         .c .c

2、/2        .s/23.若一线性移不变系统当输入为x(n)=(n)时输出为y(n)=R3(n)，则当输入为u(n)-u(n-2)时输出为( )。 .R3(n) .R2(n) .R3(n)+R3(n-1) .R2(n)-R2(n-1)4.已知序列Z变换的收敛域为z>1，则该序列为( )。 .有限长序列 .右边序列 .左边序列 .双边序列5.离散系统的差分方程为y(n)=x(n)+ay(n-1)，则系统的频率响应( )。 当a<1时，系统呈低通特性 .当a>1时，系统呈低通特性 .当0<

3、a<1时，系统呈低通特性 .当-1<a<0时，系统呈低通特性6.序列x(n)=R5(n)，其8点DFT记为X(k)，k=0,1,7，则X(0)为( )。 .2 .3 .4 .57.下列关于FFT的说法中错误的是( )。 .FFT是一种新的变换 .FFT是DFT的快速算法.FFT基本上可以分成时间抽取法和频率抽取法两类 .基2 FFT要求序列的点数为2L(其中L为整数)8.下列结构中不属于FIR滤波器基本结构的是( )。 .横截型 .级联型 .并联型 .频率抽样型9.已知某FIR滤波器单位抽样响应h(n)的长度为(M+1)，则在下列不同特性的单位抽样响应中可以用来设计线性相位滤

4、波器的是( )。 .hn=-hM-n .hn=hM+n .hn=-hM-n+1 .hn=hM-n+110.下列关于用冲激响应不变法设计IIR滤波器的说法中错误的是( )。 .数字频率与模拟频率之间呈线性关系 .能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器 .容易出现频率混叠效应 .可以用于设计高通和带阻滤波器11.利用矩形窗函数法设计FIR滤波器时，在理想特性的不连续点附近形成的过滤带的宽度近似等于( )。 .窗函数幅度函数的主瓣宽度 .窗函数幅度函数的主瓣宽度的一半 .窗函数幅度函数的第一个旁瓣宽度 .窗函数幅度函数的第一个旁瓣宽度的一半12.连续信号抽样序列在( 

5、60;      )上的Z变换等于其理想抽样信号的傅里叶变换。单位圆        .实轴.正虚轴        .负虚轴13.一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包含(        )。单位圆        .原点.

6、实轴       .虚轴14.已知某FIR滤波器单位抽样响应h(n)的长度为(M+1)，则在下列不同特性的单位抽样响应中可以用来设计线性相位滤波器的是( )。 .hn=-hM-n .hn=hM+n .hn=-hM-n+1 .hn=hM-n+115.序列x(n) = nR4(n)，则其能量等于( )。.5.10.14.2016.以下单位冲激响应所代表的线性移不变系统中因果稳定的是( )。.h(n) = u(n).h(n) = u(n +1).h(n) = R4(n).h(n) = R4(n +1)17.下列序列中z变换收敛域包

7、括z = 0的是( )。.u(n).-u(n).u(-n).u(n-1)18.实序列的傅里叶变换必是( )。.共轭对称函数.-.共轭反对称函数.线性函数.双线性函数19.已知序列x(n) =(n)，10点的DFTx(n) = X(k)（0 k 9），则X(5) =( 1 )。.10.1.0.-1020.欲借助FFT算法快速计算两有限长序列的线性卷积，则过程中要调用( )次FFT算法。.1.-.2.3.421.不考虑某些旋转因子的特殊性，一般一个基2 FFT算法的蝶形运算所需的复数乘法及复数加法次数分别为( )。.1和2.-.1和1.2和1.2和222.因果FIR滤波器的系统函数H(z)的全部极

8、点都在( )处。.z = 0.z = 1.z = j.z =23.线性相位FIR滤波器主要有以下四类h(n)偶对称，长度N为奇数 .-h(n)偶对称，长度N为偶数h(n)奇对称，长度N为奇数 h(n)奇对称，长度N为偶数则其中不能用于设计高通滤波器的是( )。24、序列u (n)的Z变换及收敛域为（ ）,1<|z| ,1< |Z|< 1,0|z| 1, 0| z|<25、序列（）n u (n)的Z变换及收敛域为（ ）,<|z|< ,|z|< <|z|< |Z|<26、序列x(n)= ()|n| 的Z变换及收敛域为（ ） |z|<

9、 |z|< <|z|<2 <|z|<227、若X(z)=, |z|>|，则X(z)的Z反变换为（ ）x (n)=() n+1u(n+1)()n-1 u(n1) x(n)=()n-1u(n+1)()n+1 u(n1)x(n)=()n-1u(n1)()n+1u(n+1) x(n)=()n-1u(n1)()n+1u(n)28、序列x(m)，h(m)分别如图所示，y(n)=x(n)*h(n)，则y(4)为（ ）01231h(m)m01231/23/21x(m)m 3 5z-1z-1-3-2x(n)z-129、下面信号流图表示的系统函数为（ ）z-1z-1-1/42/

10、31/2x(n)y(n)H(z)= H(z)= H(z)= H(z)=30、下面信号流图表示的系统函数为（ ）y(n) H(z)=1+z-1+5 z -26 z -3 H(z)=1+3 z -1z -26 z -3 H(z)=13 z -1+5 z -26 z -3 H(z)=1z -15 z -2+6 z -331、若x(n)是长度为N的实序列，且DFTx(n) =X(k)，x(n)= x(N- n)，则有（ ）X(k)=X(Nk) X(k)= X(Nk) X(k)=X*( Nk) X(k)=X(N+k) 32、对实信号进行谱分析，若要求谱分辨率F50Hz，则最小记录时间Tpmin应为 （

11、）0.5S 0.05S 0.02S 0.2S33、对实信号进行谱分析，若信号最高频率为fc=10KHz，则最大采样间隔Tmax应为（ ） 0.1×10-3S 0.01×10-3S 0.5×10-3S 0.05×10-3S34、对于N=8点的基IFFT运算，在进行位倒序后，地址单元A(4)中存放的是输入序列x(n)中的哪一个值（ ） x(1) x(2) x(4) x(0)35、已知x(n)=(n)，N点的DFTx(n)=X(k)，则X(5)=(        )。.N 

12、0;  .1    .0     .- N40、已知DFTx(n)=X(k)，下面说法中正确的是(        )。.若x(n)为实数偶对称函数，则X(k)为虚数奇对称函数.若x(n)为实数奇对称函数，则X(k)为虚数奇对称函数.若x(n)为虚数偶对称函数，则X(k)为虚数奇对称函数.若x(n)为虚数奇对称函数，则X(k)为虚数奇对称函数36、如图所示的运算流图符号是(    

13、;    )基2 FFT算法的蝶形运算流图符号。 .按频率抽取       .按时间抽取 .两者都是        .两者都不是37、直接计算N点DFT所需的复数乘法次数与(        )成正比。.N          .N2.N3   

14、;      .Nlog2N38、下列各种滤波器的结构中哪种不是I I R滤波器的基本结构(        )。.直接型        .级联型.并联型       .频率抽样型39、以下对双线性变换的描述中正确的是(        )。A.双线性变

15、换是一种线性变换.双线性变换可以用来进行数字频率与模拟频率间的变换.双线性变换是一种分段线性变换.以上说法都不对40、若模拟滤波器Ha( s)= ，采样周期T=1S，则利用双线性变换法，将Ha(S)转换成数字滤波器H(z)应为（ ） 41、抽样频率确定时，DFT的频率分辨力取决于（ ）量化误差 信号带宽 抽样间隔 抽样点数42、如果一线性移不变系统的收敛域为一半径小于1的圆的外部，则该系统为（ ）因果稳定系统 因果非稳定系统 稳定非因果系统 因果稳定系统二、多项选择题（在每小题的五个备选答案中选出而二至五个正确答案，并将其号码分别写在题干后的括号内，未选全或有选错的，该题无分。每小题1分，共1

16、2分）1、下列系统中是因果的有（ ）Tx(n)=g(n)x(n) Tx(n)= Tx(n)=ex(n) Tx(n)=ax(n)+b2、下列说法中正确的有（ ） 因果序列的Z变换收敛域为Rx-<| z |,其中Rx-为收敛域最小半径 双边序列的Z变换收敛域为圆环 左边序列的Z变换收敛域为Rx-<| z |，其中Rx-为收敛域最小半径 右边序列的Z变换收敛域为0<|z|Rx+，其中Rx+为收敛域最大半径 有限长序列的Z变换总是收敛的3、两序列卷积运算包括的步骤有（ ）翻褶 取模 平移 相乘 相加4、下列说法正确的有（ ） FIR系统只能采用非递归结构的电路 IIR体统只能采用递归

17、结构的电路 FIR系统可以采用递归或非递归结构的电路 IIR系统可以采用递归或非递归结构的电路 FIR系统的H(z)在有限Z平面上无限点6、用DFT进行谱分析时，截断后序列的频谱Y (ej)与原序列频谱X(ej)的差别对谱分析的影响主要表现在（ ）频谱混叠 泄露 衰减 谱间干扰 栅栏效应7、实现FIR滤波器的基本网络结构主要有（ ）并联型 级联型 直接型 线性相位有限脉冲响应系统网络结构 频率采样型8、利用数字域频率变换设计数字高通滤波器的步骤有（ ） 将数字高通技术指标转换成模拟高通技术要求 将模拟高通技术指标转换成模拟低通技术要求 设计模拟低通滤波器 将模拟低通转换成数字低通 利用数字域频

18、率变换将数字低通滤波器转换成数字高通滤波器9、下列说法中不正确的有（ ） 在相同技术指标下，IIR滤波器可用比FIR滤波器较少的阶数 设计微分器或积分器等主要用IIR滤波器 FIR、IIR滤波器都可用快速傅立叶变换算法 FIR滤波器可以得到严格的线性相位 对图像信号处理，采用IIR滤波器较好10、由传输函数H(z)确定状态方程和输出方程的基本方法有（ ）直接法 巢式法 部分分式法 级联法 观察法三、说明题（认为正确的，在题干后的括号内打“”；认为错的打“×”，并说明理由，否则该题无分。）1、正弦序列x(n)=Asin(n0+ )为一周期序列 （ × ）2、实序列的傅立叶变换

19、的幅度是的奇函数 （ ）3、若信号持续时间无限长，则信号的频谱无限宽 （ ）4. 线性系统必然是移不变系统。 ( × )5. 非零周期序列的Z变换不存在。 ( )6. 按时间抽取的基2 FFT算法的运算量等于按频率抽取的基2 FFT算法。( )7. 通常FIR滤波器具有递归型结构。 ( × )8. 双线性变换法是非线性变换，所以用它设计IIR滤波器不能克服频率混叠效应。(× )9、H(z)唯一地对应一个系统冲激响应h(n) （ ）10、定点制比浮点制运算速度慢 （× ）11、移不变系统必然是线性系统。 (   ×&

20、#160;   )12、当输入序列不同时，线性移不变系统的单位抽样响应也不同。(  ×    )13、离散时间系统的滤波特性可以由其幅度频率特性直接看出。 (  ×    )14、因果稳定系统的系统函数的极点必然在单位圆内。 (        )15、与FIR滤波器相似，I I R滤波器的也可以方便地实现线性相位。(   × &

21、#160;   )16.非零周期序列的能量为无穷大。( × )17.序列的傅里叶变换就是序列z变换在单位圆上的取值。( × )18.离散傅里叶变换具有隐含周期性。( × )19.FIR滤波器必是稳定的。( × )20.当输入序列不同时，线性移不变系统的单位抽样响应也不同。( × )21.离散时间系统的滤波特性可以由其幅度频率特性直接看出。( × )22.用窗函数法设计FIR低通滤波器时，可以通过增加截取长度N来任意减小阻带衰减。( × )四、填空题1数字信号处理的主要对象是 数字信号 ，采用 数值运算 的

22、方法达到处理的目的；其实现方法主要有 硬件实现 和 软件实现 。2序列x(n)的能量定义为。3对正弦信号xa=sin314t进行采样，采样频率为fs=200Hz，则所得到的采样序列x(n)=4.我们可以从三个角度用三种表示方法描述一个线性时不变离散时间系统，它们是差分方程 、 系统函数和单位脉冲响应 。5线性时不变系统是因果系统的充分必要条件是。6线性时不变系统是稳定系统的充分必要条件是。 7线性常系数差分方程的求解方法有 经典法、递推法 和变域法 。8下图是模拟信号数字处理框图，填写出各个框内的处理方法：xa(t) ya(t)9设两个有限长序列的长度分别为N和M，则他们的线性卷积的结果序列的

23、长度为 N+M-1 。10一个长序列和一个短序列卷积时，有 重叠相加法 和 重叠保留法 两种分段卷积法。12数字频率只有相对意义，因为它是实际频率对 采样 频率的 归一化频率 。13数字频率是模拟频率对 采样频率 的 归一化 化，与数字频率相对应的模拟频率是 模拟折叠频率FS/2 。14正弦序列sin(n)不一定是周期序列，比如取为无理数 时就不是。15从满足采样定理的样值信号中可以不失真地恢复出原模拟信号。采用的方法，从时域角度看是： 时域采样定理 ；从频域角度看是： 频域采样定理 16.如果通用计算机的速度为平均每次复数乘需要s，每次复数加需要1s，则在此计算机上计算210点的基2FFT需

24、要 10级蝶形运算,总的运算时间是 s。17.在用DFT近似分析连续信号的频谱时,栅栏效应是指DFT只能计算一些离散点上的频谱。18.在FIR滤波器的窗函数设计法中，常用的窗函数有矩形和哈明等等。19.试写出（n）和（t）的两点区别：（1） （n）是离散的序列而 （t） 是连续的函数 ;（2） （t）在t=0时取值无穷大，而（n）在n=0时取值为1 。21. 判定某系统为因果系统的充要条件是：时域满足条件 n<0时，h(n)=0 ，等效于在频域满足条件 收敛域一定包含无穷点 。23. 研究一个周期序列的频域特性，应该用 DFS 变换。24. 脉冲响应不变法的基本思路是： 将s平面上的传递

25、函数Ha(s)转换为z平面上的系统函数H(z) 。25 要获得线性相位的FIR数字滤波器，其单位脉冲响应h(n)必须满足条件：（1） 奇对称 ；（2） 偶对称 。27， 借助模拟滤波器的H（s）设计一个IIR高通数字滤波器，如果没有强调特殊要求的话，宜选择采用 双线性变换法 法。28， 周期序列之所以不能进行Z变换，是因为 无限长序列 。30，某DFT的表达式是X（l）= ,则变换后数字频域上相邻两个频率样点之间的间隔是。35 FIR系统成为线性相位的充要条件是 。37 用窗口法设计出一个FIR低通滤波器后，发现它过度带太宽，这样情况下宜采用的修改措施是 。40.某线性移不变系统当输入x(n)

26、 =(n-1)时输出y(n) =(n -2) + (n -3)，则该系统的单位冲激响应h(n) =_(n -1) + (n -2)_。41.序列x(n) = cos (3n)的周期等于_2/3_。43.基2 FFT算法计算N = 2L（L为整数）点DFT需_L_级蝶形，每级由_N/2_个蝶形运算组成。44.下图所示信号流图的系统函数为H(z) =_。45.线性系统同时满足_和_两个性质。46.下图所示信号流图的系统函数为_。47.有限长单位冲激响应(FIR)滤波器的主要设计方法有_窗函数_，_等波纹逼近_两种。48.将模拟滤波器映射成数字滤波器主要有_脉冲响应不变法_及双线性变换法等。50.F

27、FT的基本运算单元称为_蝶形_运算。53.N点FFT的运算量大约是 N2 /2 。54.设计一个带阻滤波器，宜用第 第一类 类FIR数字滤波器。57.描述一个离散系统的方法，时域有 差分方程 Z域有 系统函数 欲求系统瞬态解，可用 差分方程 。58.某系统函数在单位圆外有极点但它却是稳定的，则该系统一定是 因果系统 。62.在 序列为无限长的 情况下，序列傅氏变换存在，但其DFT不存在。65.设计数字滤波器的方法之一是先设计模拟滤波器，然后通过模拟S域（拉氏变换域）到数字Z域的变换，将模拟滤波器变换成数字滤波器，其中常用的双线形变换的关系式是 。66.一个因果数字系统，如果系统的极点位于Z平面

28、的 范围，则该系统是稳定的。67.模拟信号是指其时域波形 连续 的信号： 数字信号是指其时域波形 离散 的信号， 时间离散和取值离散 的信号被称为离散时间信号。69.第二类线形相位FIR滤波器（h(n)偶对称，N为偶数）一定不能用做 高通和带阻 特性的滤波系统，第三类线形相位FIR滤波器（h(n)奇对称，N为奇数）一定不能用做 低通、高通和带阻 特性的滤波系统。四、分析、作图与计算题1.判断下列系统是否为：（1）稳定系统；（2）因果系统；（3）线性系统（4）移不变系统。(1) Tx(n)=ax(n)+b(2) Tx(n)=g(n) x(n)(3) Tx(n)=ex(n)(4) Tx(n)=x(

29、n)sin2.判断下列序列是否周期序列，试确定其周期(1) (2) x(n)=ej(n/6-)3. 设有一系统，其输入输出关系有以下差分方程确定设系统是因果系统（1）求系统的单位脉冲响应；（2）系统的系统函数；（3）系统的频率函数并画出系统的幅频特性曲线。4.求下列序列的Z变换(1) x(n)=a|n|, |a| < 1(2) x(n)=(3) x(n)=sin(0n)5.已知， |a|< 1,求其逆变换x(n).6.试求下列序列的N点DFT（闭和形式表达式）(1) x(n)=acos(0n) RN(n);(2) x(n)= (n-1)+(n-2)7、画出8点按时间抽取的基2 FF

30、T算法的运算流图。8.写出16点基2 FFT算法中位序颠倒的序列号。9.画出8点按频率抽取的基2 FFT算法的运算流图。10.某线性移不变系统的单位抽样响应为: h(n)=2(n)+(n-1)+(n-3)+2(n-4)求其系统函数，并画出该系统的横截型结构（要求用的乘法器个数最少），该滤波器是否具有线性相位特性，为什么？11.用双线性变换法设计无限长单位冲激响应（I I R）数字低通滤波器，要求通带截止频率c=0.5rad，通带衰减1不大于3dB，阻带截止频率st=0.75rad，阻带衰减2不小于20dB。以巴特沃思（Butterworth）模拟低通滤波器为原型，采样间隔T=2s,写出设计步骤

31、。12.某系统的差分方程为y(n)y(n1)+y(n2)= x(n)+x(n1)（1）该系统的系统函数；（2）试用典型范型（直接II型），一阶节的级联，一阶节的并联实现此方程。13.设滤波器的差分方程为y(n)= x(n)+ x(n-1)+y(n-1)+y(n-2) 求系统的频率响应； 系统函数； 试用典型范型（直接II型）及一阶节的级联，及一阶节的并联实现此方程。五、实验题（一）对模拟周期信号进行谱分析 （1）对于周期序列，如果周期不知道，如何用FFT进行谱分析？（2）如何选择FFT的变换区间？（包括非周期信号和周期信号）（3）当N=8时，和的幅频特性会相同吗？为什么？N=16 呢？实验程序

32、清单%第10章实验3程序exp3.m% 用FFT对信号作频谱分析clear all;close all%实验内容(1)=x1n=ones(1,4); %产生序列向量x1(n)=R4(n)M=8;xa=1:(M/2); xb=(M/2):-1:1; x2n=xa,xb; %产生长度为8的三角波序列x2(n)x3n=xb,xa;X1k8=fft(x1n,8); %计算x1n的8点DFTX1k16=fft(x1n,16); %计算x1n的16点DFTX2k8=fft(x2n,8); %计算x1n的8点DFTX2k16=fft(x2n,16); %计算x1n的16点DFTX3k8=fft(x3n,8)

33、; %计算x1n的8点DFTX3k16=fft(x3n,16); %计算x1n的16点DFT%以下绘制幅频特性曲线subplot(2,2,1);mstem(X1k8); %绘制8点DFT的幅频特性图title('(1a) 8点DFTx_1(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X1k8)subplot(2,2,3);mstem(X1k16); %绘制16点DFT的幅频特性图title('(1b)16点DFTx_1(n)');xlabel('/');

34、ylabel('幅度');axis(0,2,0,1.2*max(abs(X1k16)figure(2)subplot(2,2,1);mstem(X2k8); %绘制8点DFT的幅频特性图title('(2a) 8点DFTx_2(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X2k8)subplot(2,2,2);mstem(X2k16); %绘制16点DFT的幅频特性图title('(2b)16点DFTx_2(n)');xlabel('/'

35、);ylabel('幅度');axis(0,2,0,1.2*max(abs(X2k16)subplot(2,2,3);mstem(X3k8); %绘制8点DFT的幅频特性图title('(3a) 8点DFTx_3(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X3k8)subplot(2,2,4);mstem(X3k16); %绘制16点DFT的幅频特性图title('(3b)16点DFTx_3(n)');xlabel('/');ylabe

36、l('幅度');axis(0,2,0,1.2*max(abs(X3k16)%实验内容(2) 周期序列谱分析=N=8;n=0:N-1; %FFT的变换区间N=8x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k8=fft(x4n); %计算x4n的8点DFTX5k8=fft(x5n); %计算x5n的8点DFTN=16;n=0:N-1; %FFT的变换区间N=16x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k16=fft(x4n); %计算x4n的16点DFTX5k16=fft(x5n);

37、 %计算x5n的16点DFTfigure(3)subplot(2,2,1);mstem(X4k8); %绘制8点DFT的幅频特性图title('(4a) 8点DFTx_4(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X4k8)subplot(2,2,3);mstem(X4k16); %绘制16点DFT的幅频特性图title('(4b)16点DFTx_4(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.

38、2*max(abs(X4k16)subplot(2,2,2);mstem(X5k8); %绘制8点DFT的幅频特性图title('(5a) 8点DFTx_5(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X5k8)subplot(2,2,4);mstem(X5k16); %绘制16点DFT的幅频特性图title('(5b)16点DFTx_5(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(a

39、bs(X5k16)%实验内容(3) 模拟周期信号谱分析=figure(4)Fs=64;T=1/Fs;N=16;n=0:N-1; %FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)16点采样X6k16=fft(x6nT); %计算x6nT的16点DFTX6k16=fftshift(X6k16); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率Fk=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）subplot(3,1,1);stem(fk,abs

40、(X6k16),'.');box on %绘制8点DFT的幅频特性图title('(6a) 16点|DFTx_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k16)N=32;n=0:N-1; %FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)32点采样X6k32=fft(x6nT); %计算x6nT的32点DFTX6k32=fftshif

41、t(X6k32); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率Fk=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）subplot(3,1,2);stem(fk,abs(X6k32),'.');box on %绘制8点DFT的幅频特性图title('(6b) 32点|DFTx_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k32)N=64;n=0:N-1; %F

42、FT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)64点采样X6k64=fft(x6nT); %计算x6nT的64点DFTX6k64=fftshift(X6k64); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率Fk=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）subplot(3,1,3);stem(fk,abs(X6k64),'.'); box on%绘制8点DFT的幅频特性图title('(6a) 64点|DFTx

43、_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k64)实验程序运行结果实验3程序exp3.m运行结果如图10.3.1所示。图10.3.1（二）、IIR数字滤波器设计及软件实现（1）信号产生函数mstg清单function st=mstg%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg 返回三路调幅信号相加形成的混合信号，长度N=1600N=1600 %N为信号st的长度。Fs=10000;T=1/Fs;Tp=N*T; %采样频率F

44、s=10kHz，Tp为采样时间t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10;%第1路调幅信号的载波频率fc1=1000Hz,fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调幅信号的载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40;%第3路调幅信号的载波频率fc3=250Hz,fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号xt

45、2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号st=xt1+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);title('(a) s(t)的波形'

46、;)subplot(3,1,2)stem(f,abs(fxt)/max(abs(fxt),'.');grid;title('(b) s(t)的频谱')axis(0,Fs/5,0,1.2);xlabel('f/Hz');ylabel('幅度')(2) 思考题（1）请阅读信号产生函数mstg，确定三路调幅信号的载波频率和调制信号频率。（2）信号产生函数mstg中采样点数N=800，对st进行N点FFT可以得到6根理想谱线。如果取N=1000，可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg中采样点数N的值，观察频谱图验

47、证您的判断是否正确。（3）修改信号产生函数mstg，给每路调幅信号加入载波成分，产生调幅（AM）信号，重复本实验，观察AM信号与抑制载波调幅信号的时域波形及其频谱的差别。提示：AM信号表示式：。(3)滤波器参数及实验程序清单1、滤波器参数选取观察图10.4.1可知，三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz。带宽（也可以由信号产生函数mstg清单看出）分别为50Hz、100Hz、200Hz。所以，分离混合信号st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的指标参数选取如下：对载波频率为250Hz的条幅信号，可以用低通滤波器分离，其指标为带

48、截止频率Hz，通带最大衰减dB；阻带截止频率Hz，阻带最小衰减dB，对载波频率为500Hz的条幅信号，可以用带通滤波器分离，其指标为带截止频率Hz，Hz，通带最大衰减dB；阻带截止频率Hz，Hz，Hz，阻带最小衰减dB，对载波频率为1000Hz的条幅信号，可以用高通滤波器分离，其指标为带截止频率Hz，通带最大衰减dB；阻带截止频率Hz，阻带最小衰减dB，说明：（1）为了使滤波器阶数尽可能低，每个滤波器的边界频率选择原则是尽量使滤波器过渡带宽尽可能宽。（2）与信号产生函数mstg相同，采样频率Fs=10kHz。（3）为了滤波器阶数最低，选用椭圆滤波器。按照图10.4.2 所示的程序框图编写的实验

49、程序为exp4.m。2、实验程序清单%实验4程序exp4.m% IIR数字滤波器设计及软件实现clear all;close allFs=10000;T=1/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和通带截止频率wpB,A=ellip(N,rp,rs,wp); %调用e

50、llip计算椭圆带通DF系统函数系数向量B和Ay1t=filter(B,A,st); %滤波器软件实现% 低通滤波器设计与实现绘图部分figure(2);subplot(3,1,1);myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线yt='y_1(t)'subplot(3,1,2);tplot(y1t,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=0.1;rs=60; N,

51、wp=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wpB,A=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay2t=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和通带截止频率wpB,A=elli

52、p(N,rp,rs,wp,'high'); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay3t=filter(B,A,st); %滤波器软件实现% 高低通滤波器设计与实现绘图部分（省略）3 实验程序运行结果实验4程序exp4.m运行结果如图104.2所示。由图可见，三个分离滤波器指标参数选取正确，算耗函数曲线达到所给指标。分离出的三路信号y1(n)，y2(n)和y3(n)的波形是抑制载波的单频调幅波。(a) 低通滤波器损耗函数及其分离出的调幅信号y1(t)(b) 带通滤波器损耗函数及其分离出的调幅信号y2(t)(c)高通滤波器损耗函数及其分离出的调幅信号y3(t) 图

53、104. 实验4程序exp4.m运行结果4 简要回答思考题思考题（1）已经在10.4.2节解答。思考题（3）很简单，请读者按照该题的提示修改程序，运行观察。思考题（3） 因为信号st是周期序列，谱分析时要求观察时间为整数倍周期。所以，本题的一般解答方法是，先确定信号st的周期，在判断所给采样点数N对应的观察时间Tp=NT是否为st的整数个周期。但信号产生函数mstg产生的信号st共有6个频率成分，求其周期比较麻烦，故采用下面的方法解答。分析发现，st的每个频率成分都是25Hz的整数倍。采样频率Fs=10kHz=25×400Hz，即在25Hz的正弦波的1个周期中采样400点。所以，当N为400的整数倍时一定为st的整数个周期。因此，采