数字信号处理2009试题-电信-A-答案
哈工大威海数字信号处理2009试题电信a答案
哈尔滨工业大学(威海) 2009 / 2010 学年 秋 季学期数字信号处理 试题卷( A )答案一、填空题(1~5题每空2分,其余每空1分,共20分)1、 设)2(3)1(6)()(−+−+=n n n n h δδδ为一个LSI 系统的单位抽样响应,则该系统的频率响应=)(ωj eH ωω2361j j e e −−++。
2、 0()cos()x n n ω=中仅包含频率为0ω的信号,)4/cos()()(πn n x n y =中包含的频率为40πω±。
3、 一个长度1001=N 点的序列与长度为)(n x 642=N 点的序列用N=128点的DFT计算循环卷积时,循环卷积等于线性卷积的n 的取值范围为:)(n h 12735≤≤n 。
4、 是序列[], 07X k k ≤≤[]{ -1, 2, -3, 2, 0, -4, 6, 2}x n =的8点DFT 。
则7(/4)[]j k k eX k π−==∑ 16 。
5、 的16点DFT为,)(n x )(k X )3162cos()(k k X π=,则 =)(n x 2/)]13()3([−+−n n δδ。
6、 在利用窗函数法设计FIR 滤波器时,由于窗函数截短造成滤波器通带和阻带内的波动,这种现象称为 吉布斯效应 。
7、 无限长单位冲激响应(IIR)滤波器的结构上有反馈,因此是__递归 ___型的。
8、 已知因果序列的Z 变换为,则= 1 )(n x zeZ X /1)(=)0(x 。
9、 对长度为N 的序列向右循环移位m 位得到的序列用表示,其数学表达式为 =)(n x )(n x m )(n x m )())((n R m n x N N −。
10、 的零、极点分布关于单位圆 镜像对称 )()(1−Z H Z H 。
11、 如果通用计算机的速度为平均每次复数乘需要5μs ,每次复数加需要1μs ,则在此计算机上计算210点的基2FFT 需要 10 级蝶形运算,总的运算时间是 35840 μs 。
数字信号处理期末A卷试卷答案及评分标准-09电子
4. (8 分) 以 20kHz 的采样率对最高频率为 10kHz 的带限信号������������ (������), 然后计算������(������)的 N=1000 个采样点的 DFT,即 ������(������) = � ������(������)������ −������ ������ ������������ ,������ = 1000
2������ ������������� ������
������������ (������) =
1 − ������
1 − ������
�������������0 −
�������������0 −
2������ ������������� ������ 2������ ������� ������
������=−∞
∞
8
。 四、计算题(共计 60 分) 级联型 和 1. (10 分)求下列序列������(������)的离散傅里叶变换(DFT)������(������ )。 解: ������(������) = ������ ������������0 ������ ������������ (������)
������−1 ������=0
0
。
4. 无限长单位冲激响应(IIR)滤波器的基本结构有直接 I 型、直接 II 型、
5. DFT 与 DFS 有密切关系,因为有限长序列可以看成周期序列的 主值序列 ,而周期序列 可以看成有限长序列的 周期延拓 。 6. 无限长单位冲激响应(IIR)滤波器的结构上有反馈,因此是 递归 型的。 5 。 ������(������) = �
3������ 3������ + ������������������������ � � ������ −������ 2 ������ = �������1 − ������ −������������������ � 2
(完整word版)数字信号处理期末试卷(含答案)全..(word文档良心出品)
数字信号处理期末试卷(含答案)一、单项选择题(在每小题的四个备选答案中,选出一个正确答案,并将正确答案的序号填在括号内。
1.若一模拟信号为带限,且对其抽样满足奈奎斯特采样定理,则只要将抽样信号通过( )即可完全不失真恢复原信号。
A.理想低通滤波器B.理想高通滤波器C.理想带通滤波器D.理想带阻滤波器 2.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?( )A.y(n)=x 3(n)B.y(n)=x(n)x(n+2)C.y(n)=x(n)+2D.y(n)=x(n 2)3..设两有限长序列的长度分别是M 与N ,欲用圆周卷积计算两者的线性卷积,则圆周卷积的长度至少应取( )。
A .M+NB.M+N-1C.M+N+1D.2(M+N)4.若序列的长度为M ,要能够由频域抽样信号X(k)恢复原序列,而不发生时域混叠现象,则频域抽样点数N 需满足的条件是( )。
A.N ≥MB.N ≤MC.N ≤2MD.N ≥2M 5.直接计算N 点DFT 所需的复数乘法次数与( )成正比。
A.N B.N 2 C.N 3 D.Nlog 2N6.下列各种滤波器的结构中哪种不是FIR 滤波器的基本结构( )。
A.直接型 B.级联型 C.并联型 D.频率抽样型7.第二种类型线性FIR 滤波器的幅度响应H(w)特点( ): A 关于0=w 、π、π2偶对称 B 关于0=w 、π、π2奇对称C 关于0=w 、π2偶对称 关于=w π奇对称D 关于0=w 、π2奇对称 关于=w π偶对称 8.适合带阻滤波器设计的是: ( ) A )n N (h )n (h ---=1 N 为偶数 B )n N (h )n (h ---=1 N 为奇数 C )n N (h )n (h --=1 N 为偶数D )n N (h )n (h --=1 N 为奇数9.以下对双线性变换的描述中不正确的是( )。
A.双线性变换是一种非线性变换B.双线性变换可以用来进行数字频率与模拟频率间的变换C.双线性变换把s 平面的左半平面单值映射到z 平面的单位圆内D.以上说法都不对10.关于窗函数设计法中错误的是:A 窗函数的截取长度增加,则主瓣宽度减小;B 窗函数的旁瓣相对幅度取决于窗函数的形状,与窗函数的截取长度无关;C 为减小旁瓣相对幅度而改变窗函数的形状,通常主瓣的宽度会增加;D 窗函数法不能用于设计高通滤波器; 二、填空题(每空2分,共20分)1. 用DFT 近似分析连续信号频谱时, _________效应是指DFT 只能计算一些离散点上的频谱。
数字信号处理习题及答案(精编文档).doc
【最新整理,下载后即可编辑】==============================绪论============================== 1. A/D 8bit 5V 00000000 0V 00000001 20mV 00000010 40mV 00011101 29mV==================第一章 时域离散时间信号与系统================== 1.①写出图示序列的表达式答:3)1.5δ(n 2)2δ(n 1)δ(n 2δ(n)1)δ(n x(n)-+---+++= ②用δ(n) 表示y (n )={2,7,19,28,29,15}2. ①求下列周期)54sin()8sin()4()51cos()3()54sin()2()8sin()1(n n n n n ππππ-②判断下面的序列是否是周期的; 若是周期的, 确定其周期。
(1)A是常数 8ππn 73Acos x(n)⎪⎪⎭⎫⎝⎛-= (2))81(j e)(π-=n n x解: (1) 因为ω=73π, 所以314π2=ω, 这是有理数, 因此是周期序列, 周期T =14。
(2) 因为ω=81, 所以ωπ2=16π, 这是无理数, 因此是非周期序列。
③序列)Acos(nw x(n)0ϕ+=是周期序列的条件是是有理数2π/w 0。
3.加法乘法序列{2,3,2,1}与序列{2,3,5,2,1}相加为__{4,6,7,3,1}__,相乘为___{4,9,10,2} 。
移位翻转:①已知x(n)波形,画出x(-n)的波形图。
②尺度变换:已知x(n)波形,画出x(2n)及x(n/2)波形图。
卷积和:①h(n)*求x(n),其他02n 0n 3,h(n)其他03n 0n/2设x(n) 例、⎩⎨⎧≤≤-=⎩⎨⎧≤≤=}23,4,7,4,23{0,h(n)*答案:x(n)=②已知x (n )={1,2,4,3},h (n )={2,3,5}, 求y (n )=x (n )*h (n )x (m )={1,2,4,3},h (m )={2,3,5},则h (-m )={5,3,2}(Step1:翻转)解得y (n )={2,7,19,28,29,15}③(n)x *(n)x 3),求x(n)u(n u(n)x 2),2δ(n 1)3δ(n δ(n)2、已知x 2121=--=-+-+=}{1,4,6,5,2答案:x(n)=4. 如果输入信号为,求下述系统的输出信号。
数字信号处理期末考试试题以及参考答案
数字信号处理期末考试试题以及参考答案1. 说明数字信号处理的基本概念和应用领域。
数字信号处理(Digital Signal Processing,简称DSP)是利用计算机和数字技术对信号进行处理的一种方法。
与传统的模拟信号处理相比,数字信号处理具有精度高、灵活度大以及易于集成等优势。
它广泛应用于通信、音频处理、图像处理、雷达信号处理等领域。
2. 解释采样定理的原理,并举例说明其应用。
采样定理是数字信号处理的基础理论,它规定了采样频率必须满足一定条件,以保证从连续信号中恢复出完整的原始信息。
根据采样定理,采样频率必须大于信号最高频率的两倍,即Nyquist采样频率。
例如,对于音频信号处理,人耳可以接受的最高频率为20kHz,因此需要以至少40kHz的采样频率进行采样,才能保证恢复出高质量的音频信号。
3. 描述离散时间信号和离散序列的特点,并给出示例。
离散时间信号是在离散时间点上获取的信号,相邻时间点之间存在离散性。
离散时间信号可以用离散序列来表示,离散序列是按照离散时间点取样的数字信号。
例如,某地区每天的气温是一个离散时间信号,每天不同的时间点测量一次气温,将其离散化后可以得到一个离散序列,表示该地区每天的气温变化。
4. 详述时域和频域分析在数字信号处理中的作用。
时域分析是对信号在时间上进行分析,通过观察信号的波形和幅度变化,可以了解信号的时序特性、周期性以及脉冲等特征。
频域分析是将信号变换到频率域进行分析,通过观察信号的频谱和频率特征,可以了解信号的频率分布、频率成分以及谐波情况等。
在数字信号处理中,时域分析和频域分析是互补的工具。
通过时域分析可以了解信号的时间特性,而频域分析则更适合对信号的频率特性进行研究,两者结合可以全面分析信号的性质和特点。
5. 介绍常见的数字滤波器类型,并分别阐述其特点和应用场景。
常见的数字滤波器类型有低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
- 低通滤波器:可以通过滤除高频噪声、保留低频信号来平滑信号。
数字信号处理试题和答案优选资料
一. 填空题1、一线性时不变系统,输入为x(n)时,输出为y(n);则输入为2x(n)时,输出为2y(n) ;输入为x(n-3)时,输出为y(n-3) 。
2、从奈奎斯特采样定理得出,要使实信号采样后能够不失真还原,采样频率fs与信号最高频率f max关系为:fs>=2f max。
3、已知一个长度为N的序列x(n),它的离散时间傅立叶变换为X(e jw),它的N点离散傅立叶变换X(K)是关于X(e jw)的N 点等间隔采样。
4、有限长序列x(n)的8点DFT为X(K),则X(K)= 。
5、用脉冲响应不变法进行IIR数字滤波器的设计,它的主要缺点是频谱的交叠所产生的现象。
6.若数字滤波器的单位脉冲响应h(n)是奇对称的,长度为N,则它的对称中心是(N-1)/2 。
7、用窗函数法设计FIR数字滤波器时,加矩形窗比加三角窗时,所设计出的滤波器的过渡带比较窄,阻带衰减比较小。
8、无限长单位冲激响应(IIR)滤波器的结构上有反馈环路,因此是递归型结构。
9、若正弦序列x(n)=sin(30nπ/120)是周期的,则周期是N= 8 。
10、用窗函数法设计FIR数字滤波器时,过渡带的宽度不但与窗的类型有关,还与窗的采样点数有关11.DFT与DFS有密切关系,因为有限长序列可以看成周期序列的主值区间截断,而周期序列可以看成有限长序列的周期延拓。
12.对长度为N的序列x(n)圆周移位m位得到的序列用xm(n)表示,其数学表达式为xm (n)= x((n-m))NRN(n)。
13.对按时间抽取的基2-FFT流图进行转置,并将输入变输出,输出变输入即可得到按频率抽取的基2-FFT流图。
14.线性移不变系统的性质有 交换率 、 结合率 和分配律。
15.用DFT 近似分析模拟信号的频谱时,可能出现的问题有混叠失真、 泄漏 、 栅栏效应 和频率分辨率。
16.无限长单位冲激响应滤波器的基本结构有直接Ⅰ型,直接Ⅱ型, 串联型 和 并联型 四种。
(完整)数字信号处理试卷及答案,推荐文档
数字信号处理试卷及答案1一、填空题(每空1分, 共10分)1.序列()sin(3/5)x n n π=的周期为 。
2.线性时不变系统的性质有 律、 律、 律。
3.对4()()x n R n =的Z 变换为 ,其收敛域为 。
4.抽样序列的Z 变换与离散傅里叶变换DFT 的关系为 。
5.序列x(n)=(1,-2,0,3;n=0,1,2,3), 圆周左移2位得到的序列为 。
6.设LTI 系统输入为x(n) ,系统单位序列响应为h(n),则系统零状态输出y(n)= 。
7.因果序列x(n),在Z →∞时,X(Z)= 。
二、单项选择题(每题2分, 共20分)1.δ(n)的Z 变换是 ( )A.1 B.δ(ω) C.2πδ(ω) D.2π 2.序列x 1(n )的长度为4,序列x 2(n )的长度为3,则它们线性卷积的长度是 ( )A. 3 B. 4 C. 6 D. 73.LTI 系统,输入x (n )时,输出y (n );输入为3x (n-2),输出为 ( ) A. y (n-2) B.3y (n-2) C.3y (n ) D.y (n )4.下面描述中最适合离散傅立叶变换DFT 的是 ( )A.时域为离散序列,频域为连续信号B.时域为离散周期序列,频域也为离散周期序列C.时域为离散无限长序列,频域为连续周期信号D.时域为离散有限长序列,频域也为离散有限长序列5.若一模拟信号为带限,且对其抽样满足奈奎斯特条件,理想条件下将抽样信号通过 即可完全不失真恢复原信号 A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器6.下列哪一个系统是因果系统 ( )A.y(n)=x (n+2) B. y(n)= cos(n+1)x (n) C. y(n)=x (2n) D.y(n)=x (- n)7.一个线性时不变离散系统稳定的充要条件是其系统函数的收敛域包括 ( )A. 实轴B.原点C.单位圆D.虚轴8.已知序列Z 变换的收敛域为|z |>2,则该序列为A.有限长序列B.无限长序列C.反因果序列D.因果序列 9.若序列的长度为M ,要能够由频域抽样信号X(k)恢复原序列,而不发生时域混叠现象,则频域抽样点数N 需满足的条件是 A.N≥M B.N≤M C.N≤2M D.N≥2M 10.设因果稳定的LTI 系统的单位抽样响应h(n),在n<0时,h(n)= ( )A.0 B .∞ C. -∞ D.1 三、判断题(每题1分, 共10分)1.序列的傅立叶变换是频率ω的周期函数,周期是2π。
数字信号处理试卷及详细答案(三套)
数字信号处理试卷答案完整版一、填空题:(每空1分,共18分)1、 数字频率ω是模拟频率Ω对采样频率s f 的归一化,其值是连续(连续还是离散?)。
2、 双边序列z 变换的收敛域形状为圆环或空集。
3、 某序列的DFT 表达式为∑-==1)()(N n knMWn x k X ,由此可以看出,该序列时域的长度为N ,变换后数字频域上相邻两个频率样点之间的间隔是Mπ2。
4、 线性时不变系统离散时间因果系统的系统函数为252)1(8)(22++--=z z z z z H ,则系统的极点为2,2121-=-=z z ;系统的稳定性为不稳定。
系统单位冲激响应)(n h 的初值4)0(=h ;终值)(∞h 不存在。
5、 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点的有限长序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为64+128-1=191点点的序列,如果采用基FFT2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为256点。
6、 用冲激响应不变法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的映射变换关系为Tω=Ω。
用双线性变换法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的映射变换关系为)2tan(2ωT =Ω或)2arctan(2T Ω=ω。
7、当线性相位FIR 数字滤波器满足偶对称条件时,其单位冲激响应)(n h 满足的条件为)1()(n N h n h --=,此时对应系统的频率响应)()()(ωϕωωj j e H eH =,则其对应的相位函数为ωωϕ21)(--=N 。
8、请写出三种常用低通原型模拟滤波器巴特沃什滤波器、切比雪夫滤波器、 椭圆滤波器。
二、判断题(每题2分,共10分)1、 模拟信号也可以与数字信号一样在计算机上进行数字信号处理,只要加一道采样的工序就可以了。
数字信号处理09A卷答案
NOA 卷课程名称:数字信号处理 使用班级:电信 考试形式:闭卷班级: 姓名:学号:一、单项选择题(每小题1.5分,共15分) (1) 序列)3c os()(nA n x =的周期为[ D ]。
(2) 系统的单位抽样响应h(n)= 0.3n u(n),则该系统是[ A ](3) 利用模拟滤波器设计法设计IIR 数字滤波器的方法是先设计满足相应指标的模拟滤波器,再按某种方法将模拟滤波器转换成数字滤波器。
脉冲响应不变法是一种时域上的转换方法,即它使[ D ](4) 连续时间信号等时间间隔理想抽样后,其频域信号是[ C ]。
(5) 线性移不变系统是稳定系统的充分必要条件是 [ B ](6)两个有限长序列x1(n )和x2(n ),长度分别为N1和N2,若x1(n )与x2(n )线性卷积后的结果序列为x (n ),则x (n )的长度为:[ A ](7) 用DFT 对一个32点的离散信号进行谱分析,其谱分辨率决定于谱采样的点数N ,即[ C ],分辨率越高。
(8) 一离散序列x(n),其定义域为-5≤n<∞,若其Z 变换存在,则其Z 变换X(z)的收敛域为: [ C ](9) 对连续时间非周期信号的DFT 逼近过程的步骤为:[ A ](10) 数字滤波器可分为低通、高通、带通、带阻和全通几种,其具有[ A ] 二、填空题(每小题1.5分,共15分)1) 有限长序列的Z 变换得到的X (Z ),在Z 平面的最小收敛域为 0<Z<∞ 。
2) 单位阶跃序列u (n)的Z 变换为X (Z )= Z/(Z-1) 。
全通网络是指: 对信号的任意频率分量具有相同的幅度衰减的系统 。
4) 线性移不变系统是因果系统的充分必要条件是单位抽样响应 在n<0时为0 。
5) 如离散时间系统h(n)由两子系统h 1(n)和h 2(n)并联构成,则h(n)= h 1(n)+h 2(n) 。
6) 在数字信号处理中,FIR 系统最主要优点是: 容易实现线性相位 。
数字信号处理期末考试及习题答案
北京邮电大学电信工程学院<<数字信号处理>>期末考试试题范本
(2)
z >2
x(n)
=
[(1)n 3
−
2n
]u(n);
1 3
<
z
<2
x(n)=(
1)n 3
u
(n)
+
2n
u(−n
−
1);
z
<
1 3
x(n)=[2n
−
(1 )n 3
]u(−n
−1)
三、线性卷积 (12 分)
设信号 x(n) = [1,1,1,1,3,3,3,3,1] 通过 LTI 离散系统 h(n) = [1,-1,1],分别按 下列方法计算此离散系统的输出 y(n)。 (1) 采用时域线性卷积 (2) 采用 N = 6 的重叠保留(舍去)法
北京邮电大学电信工程学院<<数字信号处理>>期末考试试题范本
<<数字信号处理>>
班级 ______ 姓名 _______ 学号 ________ 成绩 ______
一、填空 (18分, 请直接写在此试题纸的空格处)
(1) 设 h(n)是一个线性非移变系统的单位取样响应。若该系统又是因果的,则 h(n)应满足
ϖ1
∫ ∫ =
1 2π
ϖ2
(e jϖ (α −n)
ϖ1
+ e− jϖ (α −n) )dϖ
=
π
(
1 n−
α
)
ϖ2 ϖ1
cos
ϖ
(
n
−
α
)dϖ
(n
−
α
数字信号处理习题集(附答案)
第一章数字信号处理概述简答题:1.在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器,它们分别起什么作用?答:在A/D变化之前为了限制信号的最高频率,使其满足当采样频率一定时,采样频率应大于等于信号最高频率2倍的条件。
此滤波器亦称为“抗混叠”滤波器。
在D/A变换之后为了滤除高频延拓谱,以便把抽样保持的阶梯形输出波平滑化,故又称之为“平滑”滤波器。
判断说明题:2.模拟信号也可以与数字信号一样在计算机上进行数字信号处理,自己要增加一道采样的工序就可以了。
()答:错。
需要增加采样和量化两道工序。
3.一个模拟信号处理系统总可以转换成功能相同的数字系统,然后基于数字信号处理理论,对信号进行等效的数字处理。
()答:受采样频率、有限字长效应的约束,与模拟信号处理系统完全等效的数字系统未必一定能找到。
因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析,再考虑幅度量化及实现过程中有限字长所造成的影响。
故离散时间信号和系统理论是数字信号处理的理论基础。
第二章 离散时间信号与系统分析基础一、连续时间信号取样与取样定理计算题:1.过滤限带的模拟数据时,常采用数字滤波器,如图所示,图中T 表示采样周期(假设T 足够小,足以防止混叠效应),把从)()(t y t x 到的整个系统等效为一个模拟滤波器。
(a ) 如果kHz rad n h 101,8)(=π截止于,求整个系统的截止频率。
(b ) 对于kHz T 201=,重复(a )的计算。
解 (a )因为当0)(8=≥ωπωj e H rad 时,在数 — 模变换中)(1)(1)(Tj X Tj X Te Y a a j ωω=Ω=所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为8π=ΩT c因此 Hz Tf c c 6251612==Ω=π 由于最后一级的低通滤波器的截止频率为Tπ,因此对T8π没有影响,故整个系统的截止频率由)(ωj e H 决定,是625Hz 。
数字信号处理真题
数字信号处理试题
考生注意:答案写在答题纸上(包括填空题等),保持卷满面整洁。
一.填空题(每空2分,共20分)
1. 线性时不变离散因果系统的差分方程为y(n)= ―2x(n)+5x(n-1)-x(n-4),则该系统的单位脉冲响应为_______________。
五.分析计算题(共40分)
1.(8分)输入信号x(t)= + 经过一个采样频率为 =6 的理想采样系统后,又经理想低通滤波器H(j )还原,H(j )=
1/2, | |<3π 0, | |3π
求低通滤波器H(j )的输出信号y(t)。
2.(8分)已知 (n)={1,0,1}, (n)={1,1,1,1,1,},
六.设计题(共40分)
1.(10分)FFT的应用之一是快速计算线性卷积,假如一个信号序列x(n)通过一个M阶的、单位脉冲响应为h(n)的FIR滤波器,那么可以用FFT运算来快速计算滤波器的输出序列y(n),试设计一个快速求解输出序列y(n)的实现步骤,其中序列x(n)的长度设为N,
2. (10分)用脉冲响应不变法设计一个低通数字滤波器,已知模拟低通原型滤波器的传递函数为 (s)= ,系统采样频率为 ,设计该低通数字滤波器的系统函数H(z)。
3. (12分)用双线性变换法设计一个三阶巴特沃兹(Butterworth)低通数字滤波器,采样频率为 =8kHz,3dB截止频率为2kHz,已知三阶巴特沃兹滤波器的归一化低通原型为H(s)= 要求:
(1)设计该低通滤波器的系统函数H(z);
(2)画出该滤波器的直接II型(正准型)实现结构。
2.已知系统的输入输出关系为y(n)= +5,则该系统为( )
数字信号处理习题答案作者杨毅明习题解答
(11.18)
如果原来的单位脉冲响应乘上一个绝对值小于 1 的指数序列,则新的系统可以成为稳定 系统。例如:
h(n)
=
0.7 n
π sin(
n)u(n)
2
(11.12)
5
它的绝对值小于等于 0.7n,根据等比数列前 N 项之和=a1(1-rN)/(1-r),a1 是数列的首项,r 是公比,N 是数列的项的个数,单位脉冲响应(11.12)的绝对值之和
h(n) = 2(−0.5)n u(n)
(11.16)
17. (1)语言法
观察公式(2.104),h(n)是指数序列。根据因果性,y(n)的第 1 个非 0 值是输入 3x(n)产生 的,因为 h(n)是指数序列,所以后面的 y(n)是先前的 y(n-1)的 0.6 倍产生的,分析得到该 系统的差分方程是
5. 因为语文成绩的等级是离散的自变量,计算比例和表示比例时都是使用有限长的数字, 所以统计是数字信号处理。
6. 环境的温度变化是非常缓慢的,观察这种变化时没必要连续进行,记录这种温度没必要 也不可能百分之百准确,还是用数字信号处理的方法好。
7. 放大电信号的幅度是最简单的信号处理,选择功率放大器就可以增大声音的传播距离。 功率放大器是模拟电路,它是最简单的模拟信号处理电路。
天数=要求的单词量÷平均每天记忆的单词量=1 万÷平均每天记忆的单词量,达到阅读 英语书籍需要的年数=需要的学习天数÷365 天。 12. 选择数字方式较好,理由是:用数字方式记录的信号容易反复地查找、观察、分析和比 较,就像使用 CD 机查找一首歌或反复播放一首歌;地震信号的观测是每天 24 小时连续 不停的,模拟方式使用热敏纸记录或磁带记录,消耗材料较多。
信号处理习题答案
数字信号处理习题解答 第二章 数据采集技术基础有一个理想采样系统,其采样角频率Ωs =6π,采样后经理想低通滤波器H a j Ω还原,其中⎪⎩⎪⎨⎧≥Ω<Ω=Ωππ30321)(,,j H a 现有两个输入,x 1t =cos2πt ,x 2t =cos5πt ;试问输出信号y 1t ,y 2t 有无失真为什么分析:要想时域采样后能不失真地还原出原信号,则采样角频率Ωs 必须大于等于信号谱最高角频率Ωh 的2倍,即满足Ωs ≥2Ωh ; 解:已知采样角频率Ωs =6π,则由香农采样定理,可得 因为x 1t =cos2πt ,而频谱中最高角频率πππ32621=<=Ωh ,所以y 1t 无失真;因为x 2t =cos5πt ,而频谱中最高角频率πππ32652=>=Ωh ,所以y 2t 失真;设模拟信号xt =3cos2000πt +5sin6000πt +10cos12000πt ,求: (1) 该信号的最小采样频率;(2) 若采样频率f s =5000Hz,其采样后的输出信号; 分析:利用信号的采样定理及采样公式来求解;错误!采样定理采样后信号不失真的条件为:信号的采样频率f s 不小于其最高频率f m 的两倍,即f s ≥2f m错误!采样公式)()()(s nT t nT x t x n x s===解:1在模拟信号中含有的频率成分是f 1=1000Hz,f 2=3000Hz,f 3=6000Hz∴信号的最高频率f m =6000Hz由采样定理f s ≥2f m ,得信号的最小采样频率f s =2f m =12kHz 2由于采样频率f s =5kHz,则采样后的输出信号⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛====n n n n n n n n n n n f n x nT x t x n x s s nTt s522sin 5512cos 13512cos 10522sin 5512cos 35112cos 105212sin 5512cos 3562cos 10532sin 5512cos 3)()()(πππππππππππ 说明:由上式可见,采样后的信号中只出现1kHz 和2kHz 的频率成分,即kHzf f f kHzf f f ss 25000200052150001000512211======,,若由理想内插函数将此采样信号恢复成模拟信号,则恢复后的模拟信号()()t t t f t f t y ππππ4000sin 52000cos 132sin 52cos 13)(21-=-=可见,恢复后的模拟信号yt 不同于原模拟信号xt ,存在失真,这是由于采样频率不满足采样定理的要求,而产生混叠的结果;第三章 傅里叶分析I. 傅里叶变换概述习题设序列xn=δn-m ,求其频谱Xe j ω,并讨论其幅频和相频响应 分析:求解序列的频谱有两种方法:错误!先求序列的z 变换Xz ,再求频谱ωωj e z j z X eX ==)()(,即Xe j ω为单位圆上的z 变换;错误!直接求序列的傅里叶变换∑∞-∞=-=n nj j en x e X ωω)()(解:对序列xn 先进行z 变换,再求频谱,得m z m n ZT n x ZT z X -=-==)]([)]([)(δ则ωωωjm e z j e z X e X j -===)()(若系统的单位采样响应hn=xn ,则系统的频率响应)}(exp{)(1)()(ωϕωωωωωj e H e e e X e H j jm jm j j ====--•故其幅频和相频响应如图分别为幅频响应 1)(=ωj e H 相频响应 ωωϕm -=)(由图可见,该系统的频率响应具有单位幅值以及线性相位的特点; 设xn 的傅里叶变换为Xe j ω,试利用Xe j ω表示下列序列的傅里叶变换:(1))1()1()(1n x n x n x --+-= (2) )]()([21)(2n x n x n x -+=*分析:利用序列翻褶后的时移性质和线性性质来求解,即)()(ωj e X n x ⇔,)()(ωj e X n x -⇔-)()(ωωj m j e X e n m x --⇔- 解:1由于)()]([ωj e X n x DTFT =,)()]([ωj e X n x DTFT -=-,则)()]1([ωωj j e X e n x DTFT --=- )()]1([ωωj j e X e n x DTFT -=--故ωωωωωcos )(2])[()]([1j j j j e X e e e X n x DTFT ---=+= 2由于)()]([ωj e X n x DTFT **=-故)](Re[2)()()]([2ωωωj j j e X e X e X n x DTFT =+=* 设Xe j ω是如图所示的信号xn 的傅里叶变换,不必求出Xe j ω,试完成下列计算:(1))(0j e X(2) ⎰-ππωωd e X j )((3) ωππωd e X j ⎰-2)(分析:利用序列傅里叶变换的定义以及帕塞瓦定理来求解; (1) 序列的傅里叶变换公式为:正变换 ∑∞-∞=-=n nj j en x e X ωω)()(反变换 ⎰-=ππωωωπd e e X n x n j j )(21)((2) 帕塞瓦定理⎰∑-∞-∞==ππωωπd e X n x j n 22)(21)(解:1由傅里叶正变换公式可知ω=0,则6)()()(00===∑∑∞-∞=∞-∞=⋅-n n nj j n x en x e X2由于e j0=1,则由傅里叶反变换公式可知n=0,故πππωωππωππω422)(2)()(00====⋅=--⎰⎰n j j j n x d e e X d e X(3) 由帕塞瓦定理,得ππωππω28)(2)(22==∑⎰∞-∞=-n j n x d e XII. 周期序列的离散傅里叶级数DFS如图所示,序列xn 是周期为6的周期性序列,试求其傅里叶级数的系数;分析:利用DFS 的定义求解,即∑-===1)(~)](~[)(~N n kn N W n x n x DFS k X ,其中k = 0 ~ N-1解:已知N = 6,则由DFS 的定义得k jk jk j k j k j n nk j n kn eeee e en x W n x k X 5624623622626250625061068101214)(~)(~)(~ππππππ-----=-=+++++===∑∑对上式依次取k = 0 ~ 5,计算求得339)5(~33)4(~0)3(~33)2(~339)1(~60)0(~j X j X X j X j X X +=-==+=-==,,,, 设⎩⎨⎧≤≤+=n n n n x 其他,,0401)(,)2()(4-=n R n h令6))(()(~n x n x =,6))(()(~n h n h =,试求)(~n x 与)(~n h 的周期卷积;分析:可以利用列表法求解,直观方便;由于)(~)(~n x n y =错误!∑-=-=1)(~)(~)(~N m m n h m x n h 只要将列表中对应于某个n 的一行中的)(~m n h -值和第一行中与之对应的)(~m x 值相乘,然后再将所有乘积结果相加,就得到此n 的)(~n y 值 解:注意:本题需要利用下一节中有限序列与周期序列的关系以及序.........................列循环移位的概念........;.在一个周期N =6内的计算卷积值)(~)(~n x n y =错误!∑-=-=1)(~)(~)(~N m m n h m x n h 则)(~n x 与)(~n h 的周期卷积)(~n y 值n =0~5如下表所示:III. 离散傅里叶变换DFT已知xn 如图所示,为{1,1,3,2},试画出序列x-n 5,x-n 6 R 6n,xn 3R 3n,xn 6, xn-35R 5n 和xn 7 R 7n 的略图;分析:此题需注意周期延拓的数值,也就是xn N 中N 的数值;如果N 比序列的点数多,则需补零;如果N 比序列的点数少,则需将序列按N 为周期进行周期延拓,造成混叠相加形成新的序列; 解:各序列的略图如图所示;试求下列有限长序列的N 点离散傅里叶变换闭合形式表达式:(1))()(n R a n x N n = (2)N n n n n x <<-=000)()(,δ (3))()(n nR n x N = (4))()(2n R n n x N = 分析:利用有限长序列的DFT 的定义,即10)()(10-≤≤=∑-=N k W n x k X N n knN ,解:1因为)()(n R a n x N n =,所以k Nj N N n nk Njn N n knNn aea ea Wa k X ππ2121011)(--=--=--===∑∑2因为N n n n n x <<-=000)()(,δ,所以k n Njn n knNN n knNeW W n n k X 002100)()(πδ-=-===-=∑3由)()(n nR n x N =,得∑-==10)(N n knN nW k X注意:为了便于求解......,.必须利用代数简化法消除...........掉上式中的变量.......n .;.∑-=+=10)1()(N n n k N kNnW k X W NW W N WN W N W N W W W N W W W nW nWW k X kNk N N n knNkNN N k N k N k N N k N k N k N k N N n n k N N n kn Nk N-=--+--=+--=-+-+++--++++=-=-∑∑∑-=---=+-=11)1()1(])1()2(2[])1(32[)1)((11)1(32)1(321)1(1则所以kNW Nk X --=1)( 4注意:本题可利用上题的结论来进行化简...............;.由)()(2n R n n x N =,则∑-==102)(N n knN W n k X根据第3小题的结论:若)()(1n nR n x N = 则kNN n kn N W NnW k X --==∑-=1)(11 与上题同理,得kNN n knNN n knNkNN N k N k N k N N k N k N k N k N N n n k N N n kn Nk NW NN N k X N N nW N N W n N W N W N W W W N W W W W n Wn W k X ----=+--=+--=-+--=-+-+++--++++=-=-∑∑∑∑-=-=---=+-=12)2()(2)2(2)2()12()1(])1()2(4[])1(94[)1)((1111122)1(232)1(23210)1(2102所以10)1()2()(22-≤≤---=N k W N W N N k X k N kN , 试画出图示的两个有限长序列的六点循环卷积;分析:本题可以直接利用循环卷积的公式求解,也可以利用循环移位的概念来求解,即:有限长序列xn 左移mm 为正整数位的循环移位定义为)())(()(n R m n x n x N N m +=且移位时,在主值区间n =0~N-1内,当某序列值从区间的一端移出时,与它同值的序列值又从区间的另一端移入; 解:由循环卷积的定义,可知)()(1n x n y =错误!612))(([)(n x n x =错误!)(]))((662n R n x61))(([n x =错误!)(]))3((366n R n -δ )())3((3661n R n x -=则根据循环移位的概念,将序列x 1n 循环右移3个单位后乘以3并取其主值序列n =0~5即可,其结果如图所示;如图所示的5点序列xn ,试画出: (1) x nxn (2) x n 错误!xn (3) x n 错误!xn分析:本题可由图解法来计算循环卷积,并利用循环卷积来求解线性卷积;同时应注意循环卷积代替线性卷积的条件:设两个有限长序列xn 、hn 的点数分别为N 和M ,其循环卷积的长度为L ,则要用循环卷积代替线性卷积的条件是:循环卷积的长度L 必须不小于线性卷积的长度N +M-1,即L ≥N +M-1否则,在循环卷积周期延拓时会产生混叠;解:由于xn 是5点序列,所以xn xn 是5+5-1=9点序列,因此,xn 错误!xn 的前9个点n =0,1,…,8就是xn xn 值,后一个点n =9为零,因为L点循环卷积等于线性卷积结果的L 点周期延拓、混叠相加后的主值区间内的序列L 可以是任意整数值;其运算结果分别如图a 、b 、c 所示;已知两个有限长序列为⎩⎨⎧≤≤≤≤-=⎩⎨⎧≤≤≤≤+=651401)(640301)(n n n y n n n n x ,,,,试作图表示xn ,yn 以及fn =xn 错误!yn ; 分析:直接利用循环卷积公式或图解法求解; 解:其结果如图所示;习题已知xn 是N 点有限长序列,且Xk = DFT xn ;现将它补零扩展成长度为rN 点的有限长序列yn ,即⎩⎨⎧-≤≤-≤≤=1010)()(rN n N N n n x n y ,, 试求rN 点DFT yn 与Xk 的关系; 分析:利用DFT 定义求解;yn 是rN 点序列,因而结果相当于在频域序列进行插值; 解:由10)()]([)(102-≤≤==∑-=-N k en x n x DFT k X N n nk Nj ,π可得110)()()()]([)(10211-==⎪⎭⎫⎝⎛=====∑∑∑-=--=-=N l lr k r k X en x W n x Wn y n y DFT k Y N n rk n N jnkrN N n rN n nk rN,,,,, π所以在一个周期内,Yk 的采样点数是Xk 的r 倍Yk 的周期为rN ,相当于在Xk 的每两个值之间插入r-1个其它的数值不一定为零,而当k 为r 的整数l 倍时,Yk 与⎪⎭⎫⎝⎛r k X 相等; 习题频谱分析的模拟信号以8kHz 被采样,计算了512个采样点的DFT,试确定频谱采样之间的间隔,并证明你的回答; 分析:利用频域采样间隔F 0和时域采样频率f s 以及采样点数N 的关系f s =N F 0;证:由ππ2200Ω=Ω=F f s s , 得0ΩΩ=ss F f 其中Ωs 是以角频率为变量的频谱周期,Ω0是频谱采样之间的频谱间隔;又N F f ss =ΩΩ=00 则Nf F s=0 对于本题有f s =8kHz,N =512 所以 Hz F 625.1551280000==习题设有一个频谱分析用的信号处理器,采样点数必须为2的整数幂,假定没有采用任何特殊数据处理措施,要求频率分辨力≤10Hz,如果采用的采样时间间隔为,试确定: (1) 最小记录长度;(2) 所允许处理信号的最高频率; (3) 在一个记录中的最小点数; 分析:采样间隔T 和采样频率f s 满足f s =1/T ,记录长度T 0和频域分辨力F 0的关系为T 0=1/ F 0,采样定理为f s ≥2f h f h 为信号最高频率分量,一个记录中最少的采样总数N 满足002F f F f T T N hs ≥==解:1因为T 0=1/ F 0,而F 0≤10Hz,所以s T 1010≥即最小记录长度为; 2因为kHz T f s 10101.0113=⨯==,而f s ≥2f h 所以kHz f f s h 521=≤即允许处理信号的最高频率为5kHz; 31000101.01.030=⨯≥=T T N 又因N 必须为2的整数幂所以一个记录中的最少点数为N =210=1024; IV. 快速傅里叶变换FFT如果一台通用计算机的速度为平均每次复乘5μs,每次复加μs,用它来计算512点的DFT xn ,问直接计算需要多少时间,用FFT 运算需要多少时间 分析:错误!直接利用DFT 计算:复乘次数为N 2,复加次数为NN-1;错误!利用FFT 计算:复乘次数为N N2log 2,复加次数为N N 2log ; 解: 1直接计算复乘所需时间s N T 31072.151210*********=⨯⨯=⨯⨯=-- 复加所需时间s N N T 130816.0)1512(512105.0)1(105.0662=-⨯⨯⨯=-⨯⨯⨯=--所以s T T T 441536.121=+=2用FFT 计算复乘所需时间s N N T 01152.0512log 2512105log 210526261=⨯⨯=⨯⨯=-- 复加所需时间s N N T 002304.0512log 512105.0log 105.026262=⨯⨯=⨯⨯=-- 所以s T T T 013824.021=+=已知Xk ,Yk 是两个N 点实序列xn ,yn 的DFT 值,今需要从Xk ,Yk 求xn ,yn 的值,为了提高运算效率,试用一个N 点IFFT 运算一次完成;分析:我们来组成一个新的序列Xk +j Yk 序列,则有)()()]([)]([)]()([n jy n x k Y jIDFT k X IDFT k jY k X IDFT +=+=+它的实部即为实序列xn ,虚部即为实序列yn ; 解:依据题意,可知)()()()(k Y n y k X n x ⇔⇔,取序列)()()(k jY k X k Z +=对Zk 作N 点IFFT 可得序列zn ;又根据DFT 线性性质)()()]([)]([)]()([n jy n x k Y jIDFT k X IDFT k jY k X IDFT +=+=+由原题意可知,xn ,yn 都是实序列; 再根据zn = xn +j yn ,可得)](Im[)()](Re[)(n z n y n z n x ==习题, N =16时,画出基-2按时间抽取法DIT 及按频率抽取法DIF 的FFT 流图时间抽取采用输入倒位序,输出自然数顺序,频率抽取采用输入自然顺序,输出倒位序; 分析:错误!DIF法与DIT 法的异同:不同点:DIT 与DIF 的基本蝶形图不同,DIF 的复数乘法出现在减法之后,DIT 的复数乘法出现在减法之前;相同点:DIT 与DIF 的运算量是相同的;错误!DIF法与DIT 法的关系:它们的基本蝶形是互为转置的;解:1按时间抽取DIT 如图所示2按频率抽取DIF 如图所示课堂思考题若)(),(21n x n x 是因果稳定序列,求证:⎰⎰⎰---=ππωππωππωωωπωπωπ})(21}{)(21{)()(212121d e X d e X d e X e X j j j j证:设)()()(21n x n x n y *= 则由时域卷积定理,得)()()(21ωωωj j j e X e X e Y =即⎰⎰--===*ππωωωππωωωπωπd e e X e X d e e Y n y n x n x n j j j n j j )()(21)(21)()()(2121令上式的左右两边n=0,得)0()0()()()()()()(2121002102121x x k n x k x n x n x d e X e Xn n k n j j ⋅=⎥⎦⎤⎢⎣⎡-=*====-∑⎰ππωωωπ又傅里叶反变换公式,得⎰-=ππωωωπd e e X n x n j j )(21)(11,⎰-=ππωωωπd e e X n x n j j )(21)(22则⎰-=ππωωπd e X x j )(21)0(11,⎰-=ππωωπd e Xx j )(21)0(22所以⎰⎰⎰---=ππωππωππωωωπωπωπ})(21}{)(21{)()(212121d e X d e X d e X e X j j j j课堂思考题在N =16时按时间抽取的基-2FFT 算法中,若输入序列xn 采用倒位序,输出序列Xk 采用自然数顺序,试写出输入序列xn 的排列顺序,并简述理由;答:N =16的基-2FFT 算法中,输入序列xn 倒位序排列顺序为x 0、x 8、x 4、x 12、x 2、x 10、x 6、x 14、x 1、x 9、x 5、x 13、x 3、x 11、x 7、x 15;其倒位序排序规则如表所示:第五章 时域分析随机相位正弦波)sin()(0ϕω+=t x t x式中,x 0,ω均为常数,φ在0~2π内随机取值,试求其自相关函数并作图; 分析:利用自相关函数的定义求解,即⎰+=∞→TT xx dt t x t x T R 0)()(1lim)(ττ 解:由自相关函数的定义式,得[]()ωταωτααωταπτπωαωαϕωϕτωϕωττϕπϕπcos 2sin cos sin cos sin 2lim )(21)(sin )sin(1lim )()(1lim )(202202/2/200x d x R T d dt t dt t t x T dtt x t x T R T xx T T T TT xx =+====++++=+=⎰⎰⎰++-∞→-∞→∞→故且则令,可见,该随机相位正弦波的自相关函数只与角频率ω有关,而不含相...位信息...,这表明:正弦函数的自相关函数为失去了相位信息的同频率余弦函数;其自相关函数图形如图所示;两个随机相位正弦波)sin()()sin()(00ϕθωθω-+=+=t B t y t A t x式中,A 0, B 0,ω, φ均为常数,θ在0~2π内的取值概率相同,即满足R xx ττx 02/2⎪⎩⎪⎨⎧≤≤=其它,,02021)(πθπθp 试求其互相关函数并作图; 分析:利用互相关函数的定义求解,即⎰+=∞→TT xy dt t y t x T R 0)()(1lim)(ττ 解:由互相关函数的定义式,得[])cos(21)(sin )sin(21)()(1lim)(0020000ϕωτθϕθτωθωπττπ-=-+++=+=⎰⎰∞→B A d t t B A dt t y t x T R TT xy 可见,两个正弦函数的互相关函数仍为同频率的余弦函数,其最大峰值出现在τ=φ/ω处;其互相关函数图形如图所示;第六章 数字滤波器设计已知模拟滤波器的模方函数R xyττφ/ω)16)(9()4(20)(22222Ω+Ω+Ω-=Ωj H 求模拟滤波器的传递函数;分析:利用模拟滤波器的模方函数|Hj Ω|2与其传递函数H s 之间的关系式求解,即Ω=Ω=-==Ωj s j s s H s H s H j H )()()()(22解:将s=j Ω,即Ω2 = s 2代入|Hj Ω|2,得())4)(4)(3)(3()2()2(52)16)(9()4(20)()()(22222222-+-+-+=--+=-=s s s s j s j s s s s s H s H s H 可见,系统有四个极点s 1, 2=±3,s 3, 4=±4和两对零点z 1, 2=±j2;为了得到一个稳定的滤波系统,则将左半平面的极点分配给H s ;并取虚轴上的一对共轭零点作为H s 的零点,以保证H s 收敛,故模拟滤波器的传递函数为)4)(3()2)(2(52)(++-+=s s j s j s s H试设计一个巴特沃思BW 低通模拟滤波器,使滤波器的幅度响应在通带截止频率105rad/s 处的衰减不大于3 dB,在阻带截止频率4×105 rad/s 处的衰减不小于35 dB;分析:按照§中所述的巴特沃思低通滤波器的设计过程来实现;错误!先确定滤波器的阶数N由于公式1()()()()()()λλααγααγλγ,求解令令⇒⎪⎭⎪⎬⎫+ΩΩ+=Ω=+ΩΩ+=Ω=ΩΩ=ΩΩ=22221lg 10]1lg[10)(1lg 10]1lg[10)(NNcscp Nc s s s N c p p p 则滤波器的阶数公式2()()s p N ΩΩ≥lg /lg λγ注意:N 为正整数 且截止频率公式3s N c pN c Ω=ΩΩ=Ω--/1/1λγ或错误!求解位于左半S 平面上的极点公式4()N k es NN k j c k 2,,2,1212 =Ω=-+,π错误!确定N 阶巴特沃斯低通滤波器的传递函数公式5()()()()N N cNk kN c s s s s s s s s s H ---Ω=-Ω=∏= 211)( 解:错误!先确定滤波器的阶数N由题意可知,Ωp =105rad/s 时,通带最大衰减αp =3 dBΩs =4×105rad/s 时,阻带最小衰减αs =35 dB则代入公式1,求得参数γ和λ()()()()⎩⎨⎧==⇒⎪⎩⎪⎨⎧+=+=⇒⎪⎩⎪⎨⎧+=Ω=+=Ω=2.5611lg 10351lg 1031lg 10)(1lg 10)(2222λγλγλααγααs s p p 将参数γ、λ、Ωp 和Ωs 代入公式2,则滤波器的阶数()()39.2lg /lg =⇒=ΩΩ≥N N s p 取λγ将参数N 、γ和Ωp 代入公式3,可得截止频率s rad p p N c /105/1=Ω=Ω=Ω-γ错误!求解位于左半S 平面上的极点将参数Ωc 和N 代入公式4,得极点()3,2,13/)1(212=Ω=Ω=+-+k ee s k j c NN k j c k ,ππ即2/)31(2/)31(3/4323/21c j c cj c c j c j e s e s j e s Ω--=Ω=Ω-=Ω=Ω+-=Ω=πππ错误!确定巴特沃斯低通滤波器的传递函数H s将参数N 、Ωc 和s k 代入公式5,得巴特沃斯低通滤波器的传递函数式中Ωc =105rad/s()()()()322333213122)(cc c cc Nk kN c s s s s s s s s s s s s H Ω+Ω+Ω+Ω=---Ω=-Ω=∏= 试导出二阶巴特沃斯低通滤波器的系统函数,设Ωc =3 rad/s; 分析:本题利用模方函数求出其左半S 平面极点,而求得系统函数;N 阶巴特沃斯低通滤波器的模方函数定义为()Nc j j j H 2211)(ΩΩ+=Ω在上式中代入j Ω= s,可得()Nc j s s H s H 211)()(Ω+=- 而HsH-s 在左半S 平面的极点即为Hs 的极点,因此()∏=-=Nk ks s k s H 1)(其中()N k es NN k j c k ,,2,1212 =Ω=-+,π,k 0由1)(0==s s H 来确定;注意:可以证明,系数k 0=Ωc N ;解:对于二阶N =2巴特沃斯低通滤波器,其模方函数为()()4221111)(c Nc j j j j j H ΩΩ+=ΩΩ+=Ω令j Ω= s,则有()411)()(c j s s H s H Ω+=- 各极点满足()()4,3,2,1412212=Ω=Ω=+-+k ees k j c NN k j c k ,ππ则k =1, 2时,所得的s k 位于左半S 平面,即为Hs 的极点223223223223452431jes j e s j c j c --=Ω=+-=Ω=ππ 由以上两个极点构成的系统函数为()()923)(20210++=--=s s k s s s s k s H代入条件1)(0==s s H ,可得k 0 =9 注:k 0 =Ωc 2,故二阶巴特沃斯低通滤波器的系统函数9239)(2++=s s s H试导出三阶巴特沃斯低通滤波器的系统函数,设Ωc =2 rad/s; 分析:与习题6. 3同理,利用模方函数求出其左半S 平面极点,而求得系统函数;解:对于三阶N =3巴特沃斯低通滤波器,其模方函数为()()6221111)(c Nc j j j j j H ΩΩ+=ΩΩ+=Ω 令j Ω= s,则有()611)()(c j s s H s H Ω+=-各极点满足()()6,,2,1231212 ==Ω=+-+k ees k j NN k j c k ,ππ不难得知,当k =1, 2, 3时,相应的极点s k 均位于左半S 平面;则滤波器的系统函数Hs 的极点312223123432321j es e s j es j j j --==-==+-==πππ因此,三阶巴特沃斯低通滤波器的系统函数为()()()8848)(233213+++=---Ω=s s s s s s s s s s H c设模拟滤波器的系统函数为22)()(ba s as s H +++=试利用冲激响应不变法,设计IIR 数字低通滤波器;示解:将H s 展开成部分分式,得jba s jb a s b a s a s s H -++++=+++=2/12/1)()(22对H s 取拉氏反变换,得tjb a t jb a ee t h )()(2121)(--+-+=对ht 作周期为T 的等间隔采样,得[]nT jb a nTjb a nT t e e t h n h )()(21)()(--+-=+== 对hn 取Z 变换,得IIR 数字低通滤波器的系统函数为22111)(1)(0)cos 2(1)cos (1111121)()(----------+-∞=-+--=⎥⎦⎤⎢⎣⎡-+-==∑z e z bT e zbT e z e z e z n h z H aT aT aTT jb a T jb a n n设有一模拟滤波器()11)(2++=s s s H采样周期T =2,试用双线性变换法将它转变为数字系统函数Hz ; 分析:双线性变换法是模拟系统函数的S 平面和数字系统函数的Z 平面之间是一一对应的关系,消除了频谱的混叠现象,其变换关系为11112--+-=z z T s 解:将T =2代入变换公式,可得1111--+-=zz s 则数字系统函数()121112111131111111)()(11------+-=++=+⎪⎪⎭⎫ ⎝⎛+-+⎪⎪⎭⎫ ⎝⎛+-==--z z z z z z s H z H z z s用双线性变换法设计一个数字低通滤波器,采样频率f s = ,截止频率f c = 400Hz;分析:按照§中所述的采用双线性变换法的设计过程来实现;错误!利用关系式ω=T Ω将给定的模拟域频率指标转化为数字域频率指标错误!利用如下的预畸变补偿公式将数字域频率指标变换为补偿..后的..模拟域频率指标 ⎪⎭⎫⎝⎛=Ω'2tan 2ωT 错误!按补偿后的....模拟域频率指标设计三阶巴特沃斯模拟滤波器Hs 参见例6.2.4错误!利用双线性变换公式,将模拟滤波器Hs 变换为数字滤波器Hz ,即11112)()(--+-==z z T s s H z H T .为采样周期.....解:此数字滤波器的截止频率3212001400212πππω=⨯⨯==Ω=s cc c f f T 由预畸变补偿,得相应的模拟滤波器的截止频率s s c c f f T 323tan 22tan 2==⎪⎭⎫⎝⎛=Ω'πω 由习题6. 4可知,三阶巴特沃斯模拟滤波器的系统函数()()()3213)(s s s s s s s H c---Ω'=其中,滤波器的系统函数Hs 的极点πππ3432321j c j c j c es e s es Ω'=Ω'=Ω'=,,故有3223322)(cc c cs s s s H Ω'+Ω'+Ω'+Ω'= 将双线性变换公式和s c f 32=Ω'代入,可得三阶巴特沃斯数字滤波器的系统函数312111213131313211221212313313112112)1(33)1)(1(6)1()1(32)1()1(33)1()1)(1()2(2)1()1()2(2)1()2()1()()()(1111--------------+-=+-=+++-++-+-+=+Ω'++-Ω'++-Ω'+-+Ω'===----z z z z z z z z z z f z z f z f z s H s H z H c s c s c s c z z fs z z T s s请选择合适的窗函数及窗宽N 来设计一个线性相位低通滤波器⎩⎨⎧≤≤<≤=-πωωωωωαωc cj j d e e H ,,00)( 要求其阻带最小衰减为 45 dB,过渡带宽为8π/51,试求出hn 设截止频率ωc =π;分析:本题是真正实用的设计题,从中可以看到阻带衰减影响窗形状的选择当然用凯塞窗则可改变β来满足阻带衰减的要求,而窗宽N 的选择则影响过渡带宽;按照§中所述的线性相位FIR 数字低通滤波器的设计步骤来实现;错误!给定所要求的理想频率响应函数H d e j ω;错误!确定相应的理想脉冲响应序列[]⎰-==ππωωωωπd e e H e H IDTFT n h n j j d j d d )(21)()(错误!由阻带最小衰减及过渡带宽的要求,利用表6.4.1参见教材P156表,确定窗函数wn 的形状及其宽度N ;表6.4.1 常用窗函数及加窗后FIR 滤波器的特性结论:过渡带的宽度随窗宽.........N .的增加而减小......,.而阻带最小衰减则仅由..........窗的形状决定......,.不受..N .的影响...;.错误!求得所设计的FIR 滤波器的单位脉冲响应1,,1,0)()()(-==N n n w n h n h d ,错误!计算FIR 滤波器的频率响应[])()(21)()()(10ωωωωπj j d N n nj j e W e H e n h n h DTFT e H *===∑-=-检验是否满足设计要求,如不满足,则需重新设计;解:根据题目所给的低通滤波器频响的表达式H d e j ω,可得其脉冲响应序列[]()[]()αωαωπωωπωπωωωωαππωωω--====•⎰⎰---n n d e e d e e He H IDTFT n h c c c n j j n j j dj d d c csin 21)(21)()(因为题目要求设计的低通滤波器的阻带最小衰减为 45 dB,对照教材P156表可知,矩形窗、三角形窗、汉宁窗都不符合条件,所以应该选择哈明窗;由于加窗后滤波器的过渡带宽见表应小于所需的过渡带宽,即075.425186.6>∴<N N ππ 这里选窗宽N = 43,以满足要求;由于哈明窗函数)(12cos 46.054.0)(n w N n n w R ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛--=π则所设计的FIR 滤波器的单位脉冲响应()[]()⎪⎩⎪⎨⎧=--⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-==•其它,,042,,1,021215.0sin 21cos 46.054.0)()()( n n n n n w n h n h d πππ 用矩形窗设计一个FIR 线性相位低通数字滤波器;已知ωc =π,N = 21,试求出hn ;分析:本题给定的是理想线性低通滤波器,故⎩⎨⎧-<<≤<≤≤-=-c c cc j jde e H ωωππωωωωωωαω,,,0)(解:因为理想线性低通滤波器的脉冲响应序列[]()[]()αωαωπωωπωπωωωωαππωωω--====•⎰⎰---n n d e e d e e H e H IDTFT n h c c c n j j n j j d j d d ccsin 21)(21)()(由于矩形窗函数⎩⎨⎧-≤≤=其它,,0101)(N n n w R 则所设计的FIR 滤波器的单位脉冲响应()[]()()()⎪⎩⎪⎨⎧=--=--==其它,,020,,1,0102sin 10210sin )()()( n n n n n n w n h n h R d ππππ。
数字信号处理试题及参考答案
数字信号处理期末复习题一、单项选择题(在每个小题的四个备选答案中选出一个正确答案,并将正确答案的号码写在题干后面的括号内,每小题1分,共20分)1.要从抽样信号不失真恢复原连续信号,应满足下列条件的哪几条( ① )。
(Ⅰ)原信号为带限(Ⅱ)抽样频率大于两倍信号谱的最高频率(Ⅲ)抽样信号通过理想低通滤波器①.Ⅰ、Ⅱ②.Ⅱ、Ⅲ③.Ⅰ、Ⅲ④.Ⅰ、Ⅱ、Ⅲ2.在对连续信号均匀采样时,若采样角频率为Ωs,信号最高截止频率为Ωc,则折叠频率为( ④ )。
①Ωs ②.Ωc③.Ωc/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<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),则在下列不同特性的单位抽样响应中可以用来设计线性相位滤波器的是( ④ )。
数字信号处理试题A卷(电信用)
数字信号处理试卷A 卷<电信2002用)姓名: 班级: 学号: 成绩:一、选择. 共10题,每空2分(注:Z 指Z 变换>1.______________________。
A. B. C. D.2.如果y(n>=x((n+5>>8R8(n>;则< ) A., B.,3.若N1 点序列x(n>和N2点序列y(n>的线性卷积和L 点圆周卷积相等,则他们相等的条件是< )。
b5E2RGbCAP A. L=N1 B. L=N2 C. L=N1+N2-1 D. L N1+N2-1 p1EanqFDPw 4.一个因果稳定系统的系统函数H(z>的全部极点在< )。
A. 单位圆外B. 单位圆内C. 整个Z 平面D.左半平面5.如果下面的哪一个序列能够使它的成为实数< )n 03-3 -4 -5-6 -2 n6.对频域的抽样会造成时域的< ),如果x(n>是有限长序列,点数为M, 频域抽样N点,为避免混叠现象,应满足< )。
DXDiTa9E3dA.离散化B. 周期延拓C. N>M D . N<M7. 若要使有限长单位冲激响应<FIR)滤波器具有严格的第一线性相位,则其单位冲激响应h(n>应满足< )。
RTCrpUDGiTA. h(n>=-h(N-1-n>B.h(n>=+h(N-1-n>C.h(n>序列以n=0为偶对称中心D. h(n>序列以n=0为奇对称中心8.M点DFT离散谱的周期为__________。
A. M点B. M+1点C. 少于M点 D . 多于M点9.若已知模拟滤波器的幅度平方响应函数的全部零极点分布如图所示,则属于其系统函数Ha<S)的有< )5PCzVD7HxA-2jLBHrnAILg10、用单位抽样序列和及其加权和表示图所示的序列。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
,N=1000, ωk
= −0.4π
,
对应的模拟频率为 Ωk = −0.4π × 20000 = −8000π ,或者 f = −4kHz
频谱采样点之间的间隔为: Δf = fs = 20000 = 20Hz N 1000
5 (8 分)已知模拟滤波器的系统函数为 H(s) = (s + a)/((s + a)2 + b2 ) ,采样周期为 T,
由 傅 里 叶 变 换 性 质 , 实 序 列 的 DFT 必 然 为 共 轭 对 称 函 数 , 即
X (k) = X * (N − k) = X (N − k) ,所以 X (k ) 为实偶对称。
3 (10 分)如何用一个 N 点 DFT 变换计算两个实序列 x1(n) 和 x2 (n) 的 N 点 DFT 变换?
12、 实序列 x(n) 的 8 点 DFT 表示为 X (k ) ,已知 X (1) = 1 − 2 j ,则 X (7) = 1+2j 。
13、 不考虑某些旋转因子的特殊性,一般一个基 2 FFT 算法的蝶形运算所需 复数乘法及 2 次复数加法数。
1次
二、选择题(每题 2 分,共 30 分)
1 采样频率确定时,DFT 的频率分辨率取决于( A )。
两个实序列的 N 点 DFT 可由一个 N 点复序列的 DFT 求得。 首先构造一个 N 点复序列:x(n) = x1(n) + jx2 (n) ,求得 x(n) 的 N 点 DFT X (k) 后,
利用 DFT 的对称性从 X (k ) 中提取 X1(k) 和 X 2 (k) 。
X 1 (k )
=
4、 X [k], 0 ≤ k ≤ 7 是 序 列 x[n] = { -1, 2, - 3, 2, 0, - 4, 6, 2} 的 8 点 DFT 。 则
7
∑ e− j(π k / 4) X [k ] = 16 。
k =0
5、 x(n) 的 16 点 DFT 为 X (k ) , X (k) = cos( 2π 3k) , 则 x(n) = 16
哈尔滨工业大学(威海) 2009 / 2010 学年 秋 季学期
数字信号处理 试题卷( A )答案
一、填空题(1~5 题每空 2 分,其余每空 1 分,共 20 分)
1、 设 h(n) = δ (n) + 6δ (n −1) + 3δ (n − 2) 为一个 LSI 系统的单位抽样响应,则该系统的
频率响应 H (e jω ) = 1 + 6e − jω + 3e − j2ω 。
解:采样率 Ωs = 2π / Ts = 2π (20 ×103 ) ,离散频率ω 与模拟频率 Ω 的关系是:
ω
=
ΩTs
或ω
=
Ω 20000
,N
点
DFT
是对
DTFT
在
N
个频率点上的采样:
ωk
=
2π N
k
k = 0,1,L, N −1,所以, X (k ) 对应的模拟频率为:
Ωk
= 20000ωk
=
2π N
[δ (n − 3) + δ (n − 13)] / 2 。
6、 在利用窗函数法设计 FIR 滤波器时,由于窗函数截短造成滤波器通带和阻带内的波
动,这种现象称为 吉布斯效应
。
7、 无限长单位冲激响应(IIR)滤波器的结构上有反馈,因此是__递归 ___型的。
8、 已知因果序列 x(n) 的 Z 变换为 X (Z ) = e1/ z ,则 x(0) = 1
z z
−1 −1
三、计算题(共 50 分)
1 (5 分)已知 xa (t) 的傅里叶变换如下图所示,对 xa (t) 进行等间隔抽样而得 x(n) ,抽样
周期为 0.25ms,试画出 x(n) 的傅里叶变换 X (e jω ) 的图形。
X ( jΩ)
a
f -1kHz 0 1kHz
∑ X (e jω )
A 抽样点数
B. 抽样间隔
C. 信号带宽
D. 量化误差
2 如图所示的运算流图符号是( D )基 2 FFT 算法的蝶形运算流图符号。
A.按频率抽取 B.按时间抽取
WNk
−1
C.两者都是
D.两者都不是
3 已知序列 Z 变换的收敛域为|z|>0,则该序列为( B )。
A.有限长序列
B.(n)] =
DFT[Re[x(n)]] =
1 [X (k) + 2
X * ((N
−
k )) N
RN (k )]
,这是 X (k) 的共轭对称部分。
jX 2 (k)
=
DFT[ jx2 (n)] =
DFT[ j Im[x(n)]] =
1 [X (k) − 2
X * ((N
− k )) N
RN (k )]
) z −1 + e−2
aT
z
−2
,这是一个二阶递归系统。
直接型结构:
x(n) z −1
y (n) z −1
− e−aT cos bT 2e−aT cos bT z −1
− e −2aT
6 (12 分 ) 一 个 模 拟 信 号 xa (t) 通 过 一 个 模 拟 低 通 滤 波 器 , 滤 波 器 的 截 止 频 率 为
C 窗函数法和频率抽样法设计 FIR 滤波器的不同在于前者在时域后者在频域中进行。
D 用窗函数法设计 FIR 数字滤波器时,改变窗函数的种类可以改变阻带衰减。
9 下列关于用冲激响应不变法设计 IIR 滤波器的说法中错误的是( D )。
A.数字频率与模拟频率之间呈线性关系
B.能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器
fc = 2kHz ,过渡带宽度 Δf = 600Hz ,阻带衰减 50dB,该滤波器以数字方式实现,
如下图所示,用窗函数法设计一个刚好满足模拟滤波器技术指标的线性相位数字滤波
器,采样频率 f s = 10kHz 。
xa (t)
Ts
A/D
x(n) H (e jω ) y(n)
ya (t)
1
1
H (z)
=
1−
2
e−aT e− jbT
z −1
+
2
1 − e−aT e jbT
z −1
=
(1 −
1− (e−aT cos bT )z−1 e−aT e− jbT z−1)(1− e−aT e jbT
z −1 )
第 6 页 (共 8 页)
=
1
−
1− (e−aT cos bT (2e−aT cos bT )z−1
用冲激响应不变法求相应的数字滤波器的系统函数 H(z),并画出直接型实现结构。
将 H(s) 展开成部分分式, H (s) = s + a = 1/ 2 + 1/ 2 ,极点: (s + a)2 + b2 s + a + jb s + a − jb
s1,2 = −(a ± jb) ,系数 A1,2 = 1/ 2 ,所以 p1,2 = e s1,2t = e −(a± jb)T ,因此,
D.频率抽样型
12 模拟-数字转换法设计数字滤波器时,S 平面左半平面必须映射到 Z 平面的( D )。
A 实轴上
B.单位圆上 C. 单位圆外部 D. 单位圆内部
第 3 页 (共 8 页)
13 设系统的单位抽样响应为 h(n) = δ (n -1) + δ (n + 1) ,其频率响应为( A )
2、 x(n) = cos(ω0n) 中仅包含频率为ω0 的信号,y(n) = x(n) cos(nπ / 4) 中包含的频率为
ω0
±
π 4
。
3、 一个长度 N1 = 100 点的序列 x(n) 与长度为 N2 = 64 点的序列 h(n) 用 N=128 点的 DFT
计算循环卷积时,循环卷积等于线性卷积的 n 的取值范围为: 35 ≤ n ≤ 127 。
D.双边序列
4 对连续信号均匀采样,采样角频率为Ωs,信号截止频率为Ωc,则折叠频率为( D )。
A.Ωs
B.Ωc
C.Ωc/2
D.Ωs/2
5 借助 FFT 算法快速计算两有限长序列线性卷积,要调用( 3 )次 FFT 算法。
A.1
B.2
C.3
D.4
第 2 页 (共 8 页)
6 已知某 FIR 滤波器单位抽样响应 h(n)的长度为(M+1),则在下列不同特性的单位抽样响
C.容易出现频率混叠效应
D.可以用于设计高通和带阻滤波器
10 下列滤波器中,幅度响应在通带内单调下降,阻带内具有等纹波特性的是( C )。
A 巴特沃思滤波器 B 切比雪夫滤波器 C 反切比雪夫滤波器 D 椭圆滤波器
11 下列结构中不属于 IIR 滤波器基本结构的是( D )
A.直接型
B.级联型
C.并联型
=
1 T
[
r
∞ =−∞
X
a
(
jΩ
+
j
2π T
r) ,
模拟信号 xa (t) 等间隔抽样而得到的序 列 x(n) 的频谱,师 xa (t) 的频谱以 2π / T 为周期的周期性开拓,并乘以因子1/ T 。
其中数字角频率的坐标值是模拟角频率 的 T 倍。
X (e jω )
a/T
− 2π − π / 2 0 π / 2
20000k
或
fk
=
20000k N
, N=1000