数字信号处理第六章 习题及参考答案
第六章 习题及参考答案
一、习题
1、已知一个由下列差分方程表示的系统，x(n)、y(n)分别表示该系统的输
：
)1(2
1)()2(6
1)1(6
5)(-+
=-+
--
n x n x n y n y n y
（1）画出该系统的直接型结构； （2）画出该系统的级联型结构； （3）画出该系统的并联型结构。
2、已知某系统的系统函数为：
)
6.09.01)(5.01()
9.21)(1()(2
1
1
2
11------++-+-+=
z z
z
z
z z z H
请画出该系统的级联型结构。
3、已知FIR 滤波器的单位脉冲响应为)(8.0)(5n R n h n
=，
（1）求该滤波器的系统函数； （2）画出该滤波器的直接型结构。
4、已知滤波器的系统函数为：
3
2
1
3218.09.09.018.04.16.01)(-------+-+--=
z
z
z
z z z z H
请画出该滤波器的直接型结构。
5、已知滤波器的系统函数为：
)
8.027.11)(5.01()
44.11)(1(3)(2
1
1
2
11------+--+--=
z z z
z
z
z z H
请画出该滤波器的级联型结构和并联型结构。
6、已知某因果系统的信号流图如下图所示：
x(n)
y(n)
-2
5
-3
求该系统的系统函数和单位脉冲响应。 7、已知某系统的信号流图如下图所示：
x(n)
y(n)
求该系统的系统函数和极点。
8、已知IIR 滤波器的系统函数为：
4
.035.04.046.16.14)(2
3
2
3++++--=
z z z z z z z H
（1）画出级联型网络结构，要求利用MA TLAB 分解H(z); （2）用MA TLAB 验证所求的级联型结构是否正确。
9、已知IIR 滤波器的系统函数为：
3
2
1
3
2
1
4.03
5.04.01
6.141.158.12.5)(-------++-++=
z
z z
z
z
z z H
（1）画出该系统的并联型网络结构，要求用MA TLAB 分解； （2）用MA TLAB 验证（1）中所求的并联型结构是否正确。
10、已知FIR 滤波器的系统函数为：
4
3
2
1
12.011.072.03.01)(----++++=z
z
z
z
z H
请画出该系统的级联型结构，要求用MA TLAB 分解。
11、已知FIR 滤波器的系统函数为：
5
4
3
2
1
820520581)(-------+-+=z
z
z
z
z
z H
请画出该系统的直接型结构。
12、已知FIR 滤波器的系统函数为：
4
3
2
1
3531)(----++++=z
z
z
z
z H
请画出该系统的线性相位型结构。
13、已知FIR 滤波器的系统函数为：
3
2
1
72.16.18.02)(---+--=z
z
z
z H
要求用MA TLAB 求出该系统的全零点格型结构反射系数，并画出该全零点格型结构。
14、已知IIR 滤波器的系统函数为：
3
2
1
34711
)(----+-=
z
z
z
z H
要求用MA TLAB 求出该系统的全极点格型结构反射系数，并画出该全极点格型结构。
15、已知IIR 滤波器的系统函数为：
3
2
1
32125131321)(------++-+++=
z
z
z
z
z z z H
请画出该零-极点系统的格型结构。
二、答案
1、（1）2
1
1
6
16512
11)(---+
-
+=
z
z z z H
x(n)
y(n)
（2）两种级联型结构： a)
1
1
1
3111211211)(----
?
-
+
=
z
z
z z H
b)
1
1
1
3
112112111)(----
+?-
=
z
z z
z H
（3）1
1
31132116)(---
-+-
=
z
z
z H
x(n)
y(n)
2、 两种情况： a)
2
1
21
1
16.09.019.215.011)(------+++-?
-+=
z
z z
z
z
z
z H
x(n)
y(n)
b)
2
1
1
1
2
1
6.09.0115.019.21)(------+++?
-+-=
z
z z
z
z
z
z H
x(n)
y(n)
3、(1) 1
5
8.0132768.01)(----=
z
z z H
(2)
x(n)
y(n)
4、
x(n)
y(n)
5、(1)级联型：
x(n)
y(n)
(2)并联型：2
1
1
1
8.027.117549.06687.05.018313.35.7)(----+-+-+
--+
=z
z z z
z H
x(n)
y(n)
-0.8
6、1
11
513131112)(-----+
-
+
-=z
z z
z H
)
1()5(3)(5)()3
1()(2)(1
-?-++-=-n u n u n u n n h n n n δ
7、系统函数：)
2)(2(1
)(1
1
θ
θ
j j e
z
e
z
z H --=
---
极点：θ
θ
j j e
z e
z -==2,221
8、（1）整理得：3
2
1
321
4.03
5.04.014
6.16.14)(------++++--=z
z z
z
z
z z H
>>b=[4,-1.6,-1.6,4]; a=[1,0.4,0.35,0.4]; [sos,G]=tf2sos(b,a) 运行结果：
sos=1 1 0 1 -0.5 0
1 -1.4 1 1 0.9 0.8 G=4 所以，)
8.09.01)(5.01()
4.11)(1(4)(2
1
1
2
1
1
------++-+-+=
z z
z
z
z
z z H
级联结构有两种： a)
x(n)
y(n)
b)
x(n)
y(n)
（2）设输入信号x(n)为一个长度为6的单位脉冲序列，y1表示输入信号通过级联型滤波器的输出，y2表示输入信号通过直接型滤波器的输出，在命令窗口输入以下命令： >>format long; delta=impseq(0,0,5) y1=G*sosfilt(sos,delta) y2=filter(b,a,delta)
最后比较y1和y2，若两者相等，则验证了（1）中所求的级联型
结构是正确的。
9、（1）>>b=[5.2,1.58,1.41,-1.6]; a=[1,0.4,0.35,-0.4]; [C,B,A]=tf_par(b,a) 运行结果为： C=4 B=0.2 0 1 0.3 A=1 -0.5 0 1 0.9 0.8 根据运行结果，得：
2
1
1
1
8.09.013.015.012.04)(----++++
-+
=z
z
z z
z H
所以并联结构为：
x(n)y(n)
-0.8
（2）设输入x(n)是一个长度为6的单位脉冲序列，y1表示输入信号通过（1）中并联型滤波器的输出，y2表示输入x(n)通过直接型滤波器的输出，在MA TLAB命令窗口输入以下命令：
>>format long;
delta=impseq(0,0,5)
y1=par_filt(C,B,A,delta)
y2=filter(b,a,delta)
最后比较y1和y2，若两者相等，则验证了（1）中所求的并联型结构是正确的。
10、>>b=[1,0.3,0.72,0.11,0.12]
a=1
[sos,G]=tf2sos(b,a)
运行结果: sos=
1 0.
2 0.
3 1 0 0 1 0.1 0.
4 1 0 0 G=1
所以，)4.01.01)(3.02.01()(2121----++++=z z z z z H 所以，级联型结构为：
x(n)
y(n)
11、
x(n)
z -1
y(n)
z -1 z -1 z -1 z -1
12、2
3
1
4
5)(3)(1)(----++?++?=z
z
z
z
z z H
x(n)
z -1
y(n)
13、首项归一化：
3
2
1
______
86.08.04.012
)()(---+--==
z
z
z
z H z H
>>b=[1,-0.4,-0.8,0.86]; 运行程序，得： K= -1.4724 -1.7512 0.8600
所以，反射系数为：k1=-1.4724, k2= -1.7512, k3= 0.8600 所以全零点格型结构如下：
x(n)
z -1 z -1
z -1
14、>>b=[1,-7,4,-3] K=tf2latc(b)
运行程序，得：
K=
-0.2000
2.1250
-3.0000
所以反射系数为:
k1=-0.2,k2=2.125,k3=-3 所以全极点格型结构如下：x(n)
15、>>b=[1,2,3,1]
a=[1,-13,5,2]
[K,C]=tf2latc(b,a)
运行程序，得：
K =
-0.8214
-10.3333
2.0000
C =
61.1071
-125.6667
16.0000
1.0000
所以，该零-极点系统的格型结构如下：
x(n)
y(n)
数字信号处理习题集(附答案)
第一章数字信号处理概述 简答题： 1．在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器，它们分别起什么作用？ 答：在A/D变化之前为了限制信号的最高频率，使其满足当采样频率一定时，采样频率应大于等于信号最高频率2倍的条件。此滤波器亦称为“抗混叠”滤波器。 在D/A变换之后为了滤除高频延拓谱，以便把抽样保持的阶梯形输出波平滑化，故又称之为“平滑”滤波器。 判断说明题： 2．模拟信号也可以与数字信号一样在计算机上进行数字信号处理，自己要增加一道采样的工序就可以了。 （） 答：错。需要增加采样和量化两道工序。 3．一个模拟信号处理系统总可以转换成功能相同的数字系统，然后基于数字信号处理理论，对信号进行等效的数字处理。（） 答：受采样频率、有限字长效应的约束，与模拟信号处理系统完全等效的数字系统未必一定能找到。因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析，再考虑幅度量化及实现过程中有限字长所造成的影响。故离散时间信号和系统理论是数字信号处
理的理论基础。 第二章 离散时间信号与系统分析基础 一、连续时间信号取样与取样定理 计算题： 1．过滤限带的模拟数据时，常采用数字滤波器，如图所示，图中T 表示采样周期（假设T 足够小，足以防止混叠效应），把从)()(t y t x 到的整个系统等效为一个模拟滤波器。 （a ） 如果kHz T rad n h 101,8)(=π截止于，求整个系统的截止频 率。 （b ） 对于kHz T 201=，重复（a ）的计算。 采样（T） () n h () n x () t x () n y D/A 理想低通T c πω=() t y 解 （a ）因为当0)(8=≥ω πωj e H rad 时，在数 — 模变换中 )(1)(1)(T j X T j X T e Y a a j ωω=Ω= 所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为 8 π = ΩT c 因此 Hz T f c c 625161 2==Ω= π
数字信号处理复习题带答案
1.若一模拟信号为带限信号，且对其抽样满足奈奎斯特条件，则只要将抽样信号通过 _____A____即可完全不失真恢复原信号。 A 、理想低通滤波器 B 、理想高通滤波器 C 、理想带通滤波器 D 、理想带阻滤波器 2.下列哪一个单位抽样响应所表示的系统不是因果系统___D__? A 、.h(n)=δ(n)+δ(n -10) B 、h(n)=u(n) C 、h(n)=u(n)-u(n-1) D 、 h(n)=u(n)-u(n+1) 3.若序列的长度为M ，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N 需满足的条件是_____A_____。 A.N≥M B.N≤M C.N≤2M D.N≥2M 4.以下对双线性变换的描述中不正确的是__D_________。 A.双线性变换是一种非线性变换 B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s 平面的左半平面单值映射到z 平面的单位圆内 D.以上说法都不对 5、信号3(n)Acos(n )78 x ππ =-是否为周期信号，若是周期信号，周期为多少？ A 、周期N= 37 π B 、无法判断 C 、非周期信号 D 、周期N=14 6、用窗函数设计FIR 滤波器时，下列说法正确的是___a____。 A 、加大窗函数的长度不能改变主瓣与旁瓣的相对比例。 B 、加大窗函数的长度可以增加主瓣与旁瓣的比例。 C 、加大窗函数的长度可以减少主瓣与旁瓣的比例 。 D 、以上说法都不对。 7.令||()n x n a =，01,a n <<-∞≤≤∞，()[()]X Z Z x n =，则()X Z 的收敛域 为 __________。 A 、1||a z a -<< B 、1||a z a -<< C 、||a z < D 、1||z a -< 。
数字信号处理实验作业
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
数字信号处理经典例题解析
1：周期序列()()n n x 0cos ~ ω=， 0ω6 π =，()n x ~是由)(~ t x a ()t 0cos Ω=理想抽样而得。试求(1)()n x ~的周期; (2)()()[]n x F e X j ~ =ω (3) ()t x a ~=∑∞ -∞ =n nt j n 0 e Ωα；求n α (4) ()()[]t x F X a ~ =Ω 解：(1) 对于周期性序列()()n n x 0cos ~ ω= 因为 2ωπ = 6/2ππ =112=K N 所以序列周期12=N (2):由题意知()n x ~是由()t x a ~ 理想抽样所得，设抽样间隔为s T ，抽样输出为()t x a ?； 易得()()[]t x F X a ~ =Ω()[]t F 0cos Ω= ]2 [00t j t j e e F Ω-Ω+= =π()0Ω+Ωδ+π()0Ω-Ωδ 由采样序列()n x ~=()nt x a ?，由采样定理知： () ()[]n x F e X j ~=ω=()s T a X /?ω=ΩΩ =∑∞ ∞ --k s s s T k T X T )2( 1πω = ∑∞ ∞ --k s s T k X T )2(1 πω
=)]26()26([1s k s s T k T k T π πωπδππωπδ-++--∑∞∞- =)]26()26([ππ ωπδππωπδk k k -++--∑∞ ∞ - (3) 由)(~t x a ()t 0cos Ω== 2 00t j t j e e Ω-Ω+=∑∞ -∞ =n nt j n 0 e Ωα得: ?????=±==其他 n n n 0121 α (4)由(2)得:()ΩX =π()0Ω+Ωδ+π()0Ω-Ωδ 2：有限长序列()?? ? ??=n n x 6cos π ()n R 12求： (1))]([)(n R F e R n j n =ω (2) ()()[]n x F e X j =ω，用)(ωj N e R 表示; (3)求(2)中() ωj e X 的采样值??? ? ??k j e X 122 π 110≤≤k ; (4)()()[]n x DFT k X =; (5):求第(3)问中??? ? ??k j e X 122 π 的IDFT 变换; (6)：求() ()????????? ??=n R n F e X j 2416cos πω 的采样值??? ? ??k j e X 2421π 230≤≤k ; (7):求第(6)问中的采样序列()n x 1; (8):第(2)问中() ωj e X 的采样值??? ? ??k j e X 242 π 对应的采样序列。 .解:(1))]([)(n R F e R n j n =ω =∑-=1 )(N n n j N e n R ω
数字信号处理习题及答案1
数字信号处理习题及答案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. 7 3．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)
数字信号处理实验作业
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
《数字信号处理》复习题及答案
《数字信号处理》复习题 一、单项选择题（在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在题干的括号内。每小题2分） 1.在对连续信号均匀采样时，若采样角频率为Ωs，信号最高截止频率为Ωc，则折叠频率为( D)。 A. Ωs B. Ωc C. Ωc/2 D. Ωs/2 2. 若一线性移不变系统当输入为x(n)=δ(n)时输出为y(n)=R3(n)，则当输入为u(n)-u(n-2)时输出为( C)。 A. R3(n) B. R2(n) C. R3(n)+R3(n-1) D. R2(n)+R2(n-1) 3. 一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包含( A)。 A. 单位圆 B. 原点 C. 实轴 D. 虚轴 4. 已知x(n)=δ(n)，N点的DFT［x(n)］=X(k)，则X(5)=( B)。 A. N B. 1 C. 0 D. - N 5. 如图所示的运算流图符号是( D)基2 FFT算法的蝶形运算流图符号。 A. 按频率抽取 B. 按时间抽取 C. 两者都是 D. 两者都不是 6. 直接计算N点DFT所需的复数乘法次数与( B)成正比。 A. N B. N2 C. N3 D. Nlog2N 7. 下列各种滤波器的结构中哪种不是I I R滤波器的基本结构( D)。 A. 直接型 B. 级联型 C. 并联型 D. 频率抽样型 8. 以下对双线性变换的描述中正确的是( B)。 A. 双线性变换是一种线性变换 B. 双线性变换可以用来进行数字频率与模拟频率间的变换 C. 双线性变换是一种分段线性变换 D. 以上说法都不对 9. 已知序列Z变换的收敛域为｜z｜>1，则该序列为( B)。 A. 有限长序列 B. 右边序列 C. 左边序列 D. 双边序列 10. 序列x(n)=R5(n)，其8点DFT记为X(k)，k=0,1,…,7，则X(0)为( D)。 A. 2 B. 3
数字信号处理实验1认识实验
实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）；
数字信号处理习题及答案
三、计算题 1、已知10),()(<<=a n u a n x n ,求)(n x 的Z 变换及收敛域。 （10分） 解：∑∑∞ =-∞ -∞=-= = )()(n n n n n n z a z n u a z X 1 111 )(-∞=--== ∑ az z a n n ||||a z > 2、设)()(n u a n x n = )1()()(1--=-n u ab n u b n h n n 求 )()()(n h n x n y *=。（10分） 解：[]a z z n x z X -=? =)()(， ||||a z > []b z a z b z a b z z n h z H --=---= ?=)()(， ||||b z > b z z z H z X z Y -= =)()()( ， |||| b z > 其z 反变换为 [])()()()()(1n u b z Y n h n x n y n =?=*=- 3、写出图中流图的系统函数。(10分) 解：2 1)(--++=cz bz a z H 2 1124132)(----++= z z z z H ４、利用共轭对称性，可以用一次DFT 运算来计算两个实数序列的DFT ，因而可以减少计算量。设都是N 点实数序列，试用一次DFT 来计算它们各自的DFT ： [])()(11k X n x DFT = []) ()(22k X n x DFT =（10分）。 解:先利用这两个序列构成一个复序列，即 )()()(21n jx n x n w +=
即 [][])()()()(21n jx n x DFT k W n w DFT +== []()[]n x jDFT n x DFT 21)(+= )()(21k jX k X += 又[])(Re )(1n w n x = 得 [])(})({Re )(1k W n w DFT k X ep == [] )())(()(2 1*k R k N W k W N N -+= 同样 [])(1 })({Im )(2k W j n w DFT k X op == [] )())(()(21*k R k N W k W j N N --= 所以用DFT 求出)(k W 后，再按以上公式即可求得)(1k X 与)(2k X 。 5、已知滤波器的单位脉冲响应为)(9.0)(5n R n h n =求出系统函数，并画出其直接型 结构。（10分） 解： ｘ（ｎ） 1-z 1-z 1-z 1-z １ 9.0 2 9.0 3 9.0 4 9.0 ｙ（ｎ） ６、略。 ７、设模拟滤波器的系统函数为 31 11342)(2+-+=++=s s s s s H a 试利用冲激响应不变法，设计IIR 数字滤波器。（10分） 解 T T e z T e z T z H 31111)(-------=
数字信号处理习题集大题及答案
1设序列x(n)={4，3，2，1} ， 另一序列h(n) ={1，1，1，1}，n=0,1,2,3 （1）试求线性卷积 y(n)=x(n)*h(n) （2）试求6点圆周卷积。 （3）试求8点圆周卷积。 解：1．y(n)=x(n)*h(n)={4,7,9,10,6,3,1} 2．6点圆周卷积={5,7,9,10,6,3} 3．8点圆周卷积={4,7,9,10,6,3,1,0} 2二．数字序列 x(n)如图所示. 画出下列每个序列时域序列： (1) x(n-2); (2)x(3-n); (3)x[((n-1))6],(0≤n ≤5); (4)x[((-n-1))6],(0≤n ≤5); n 1 2 3 4 0.5 4 3210-1-2-3x(3-n) x[((n-1))6] n 5432104 3 2 1 0.5 n 1 2 3 4 0.5 5 43210x[((-n-1))6] 3．已知一稳定的LTI 系统的H(z)为) 21)(5.01() 1(2)(111------=z z z z H 试确定该系统H(z)的收敛域和脉冲响应h[n]。 解: 0.5 2Re Im 系统有两个极点，其收敛域可能有三种形式，|z|<0.5, 0.5<|z|<2, |z|>2 因为稳定，收敛域应包含单位圆，则系统收敛域为：0.5<|z|<2 1 1 111213 /25.013/4)21)(5.01()1(2)(--------=---=z z z z z z H )1(23 2 )()5.0(34)(--+= n u n u n h n n
4．设x(n)是一个10点的有限序列 x （n ）={ 2,3,1,4,-3,-1,1,1,0,6}，不计算DFT ，试确定下列表达式的值。 (1) X(0), (2) X(5), (3) ∑=9 0)(k k X ,（4） ∑=-9 5 /2)(k k j k X e π 解：（1） （2） （3） （4） 5． x(n)和h(n)是如下给定的有限序列 x(n)={5, 2, 4, -1, 2}， h(n)={-3, 2, -1 } (1) 计算x(n)和h(n)的线性卷积y(n)= x(n)* h(n)； (2) 计算x(n)和h(n)的6 点循环卷积y 1(n)= x(n)⑥h (n)； (3) 计算x(n)和h(n)的8 点循环卷积y 2(n)= x(n)⑧h (n)； 比较以上结果，有何结论？ 解：（1） 5 2 4 -1 2 -3 2 1 5 2 4 -1 210 4 8 -2 4-15 - 6 -12 3 -6 -15 4 -3 13 -4 3 2 14 ][]0[1 9 0===∑=n N n x X W 12 ][][]5[1 19 180510 -=-= ==???-=∑∑====奇 偶 奇数 偶数n n n n n n x n x X n n W 20 ]0[*10][] [101]0[9 9 ===∑∑==x k X k X x k k 0 ]8[*10][] [101]))210[((] []))[((2 )10/2(9 2 )10/2(9 10)/2(===-? --=-=-∑∑x k X e k X e x k X e m n x k j k k j k m N k j N πππ
数字信号处理第二章上机题作业
数字信号处理作业实验题报告 第一章16.(1) 实验目的： 求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。 实验要求： 运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。 源程序： B1=1;A1=[1, -0.6, 0.08]; ys=2; %设差分方程 xn=[1, zeros(1, 20)]; %xn=单位脉冲序列，长度N=31 xi=filtic(B1, A1, ys); hn1=filter(B1, A1, xn, xi); %求系统输出信号hn1 n=0:length(hn1)-1; subplot(2, 1, 1);stem(n, hn1, '.') title('单位脉冲响应'); xlabel('n');ylabel('h(n)') xn=ones(1, 20); sn1=filter(B1, A1, xn, xi); %求系统输出信号sn1 n=0:length(sn1)-1; Subplot(2, 1, 2); stem(n, sn1, '.') title('单位阶跃响应'); xlabel('n'); ylabel('s(n)')
运行结果： 实验分析： 单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。
第二章31题 实验目的： 用matlab判断系统是否稳定。 实验要求： 用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。 源程序： A=[2, -2.98, 0.17, 2.3418, -1.5147]; B=[0, 0, 1, 5, -50]; subplot(2,1,1); zplane(B,A); %求H(z)的极点 p=roots(A); %求H(z)的模 pm=abs(p); if max(pm)<1 disp('系统因果稳定'), else,disp('系统因果不稳定'),end un=ones(1,800); sn=filter(B, A, un); n=0:length(sn)-1; subplot(2, 1, 2);plot(n, sn) xlabel('n');ylabel('s(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)(T j X T j X T e Y a a j ωω=Ω= 所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为 8 π = ΩT c 因此 Hz T f c c 625161 2==Ω= π
(完整版)数字信号处理复习题-答案
、填空题 1．序列x(n) sin(3 n / 5)的周期为10 。2．线性时不变系统的性质有交换律律结合律分配律。 3．从奈奎斯特采样定理得出，要使实信号采样后能够不失真还原，采样频率 f 与信号最高频率fs 关系为：f>=2fs 4．若正弦序列x(n)=sin(30n π/120) 是周期的，则周期是N= 8 。 5．序列x(n) sin(3 n / 5)的周期为10 。 6．设LTI 系统输入为x(n) ，系统单位序列响应为h(n)，则系统零状态输出y(n)= 7．因果序列x(n) ，在Z→∞时，X(Z)= x(0) 。二、单项选择题 1．δ (n)的傅里叶变换是( A ) A. 1 B.δ (ω ) C.2πδ (ω) D.2π 2．序列x1(n)的长度为4，序列x2( n)的长度为3，则它们线性卷积的长度是( C ) A. 3 B. 4 C. 6 D. 7 3．LTI 系统，输入x(n)时，输出y( n)；输入为3x (n-2)，输出为( B ) A. y （n-2） B.3y （n-2） C.3y （n） D.y(n) 4．下面描述中最适合离散傅立叶变换DFT 的是(D ) A. 时域为离散序列，频域为连续信号 B. 时域为离散周期序列，频域也为离散周期序列 C. 时域为离散无限长序列，频域为连续周期信号 D. 时域为离散有限长序列，频域也为离散有限长序列 5．设系统的单位抽样响应为h(n)，则系统因果的充要条件为( C ) A．当n>0 时，h(n)=0 B．当n>0时，h(n) ≠0
6．下列哪一个系统是因果系统( 5．所谓采样，就是利用采样脉冲序列 p(t) 从连续时间信号 x a (t)中抽取一系列的离散样值。 ( 6．数字信号处理只有硬件方式实现。 ( × ) 7．对正弦信号进行采样得到的正弦序列一定是周期序列。 ( × ) 8．数字信号处理仅仅指的是数字处理器。 ( × ) 9．信号处理的两种基本方法：一是放大信号，二是变换信号。 ( × 10．在时域对连续信号进行抽样，在频域中，所得频谱是原信号频谱的周期延拓。 ( × ) 四、简答题 1．用 DFT 对连续信号进行谱分析的误差问题有哪些？ 答：混叠失真；截断效应(频谱泄漏) ；栅栏效应 2．画出模拟信号数字化处理框图，并简要说明框图中每一部分的功能作用。 1 2 3 部分：按照预制要 求对数字信号处理加工； 第 4部分：数字信号变为模拟信号； 第 5 部分：滤除高频部分， 平滑模拟信号。 A.N ≥M B.N ≤M C.N ≤ 2M D.N ≥ 2M 10 ．设因果稳定的 LTI 系统的单位抽样响应 h(n) ， 在 n<0 时， h(n)= ( A ) A.0 B.∞ C. - ∞ D.1 三、 判断题 1． 序列的傅立叶变换是频率ω的周期函数，周期是 2π。 ( √ ) 2 ． x(n)= sin (ω ( √ ) 0n) 所代表的序列不一定是周期 3． 卷积的计算过程包括翻转，移位，相乘，求和四个过程 ( √ ) 4． y(n)=cos[x(n)] 所代表的系统是非线性系统。 ( √ ) ) 则频域抽样点数 N 需满足的条件是 ( A C ．当 n<0 时， h(n)=0 D ．当 n<0 时， h(n) ≠0 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. x(n)= δ (n-3)的傅里叶变换为( A e 3jw B. e 3jw C.1 D.0 x(n) a n u(n),0 a 1 的傅里叶变换为 11 A. jw B. jw 1 ae 1-ae 8． C ) 1 C. -jw 1-ae 1 D.1 ae - jw 9．若序列的长度为 M ，要能够由频域抽样信号 X(k) 恢复原序列，而不发生时域混叠现象， √)
《数字信号处理》课程实验题目
计电学院《数字信号处理》课程实验 适用专业：电子通信工程专业；实验学时：9 学时 一、实验的性质、任务和基本要求 (一)本实验课的性质、任务 数字信号处理课程实验是数字信号处理课程的有效的补充部分，通过实验，使学生巩固和加深数字信号处理的理论知识的理解和掌握，在实验过程中了解简单但是完整的数字信号处理的工程实现方法和流程。通过实践进一步加强学生独立分析问题和解决问题的能力、实际动手能力、综合设计及创新能力的培养。 （二）基本要求 掌握数字信号处理基本理论知识和滤波器设计及应用。 （三）实验选项
二、实验教学内容 实验一 1、实验目的和要求 1）加深理解时域采样定理、体会使用MATLAB的离散FT函数fft( )来解决涉及模拟信号的问题； 2）加深理解对带通信号的采样特性，学会采用MATLAB解决该问题； 3）加深理解在频率采样法中，过渡点对所设计滤波器特性的影响。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第一组：张毅雷凌峰白法聪覃昱滔刘强何新文 第二组：邓志强林盛勇李日胜黎少锋梁聪杨晨 实验二 1、实验目的和要求 （1）加深理解采用数字信号处理方法对模拟信号处理的过程、掌握使用MATLAB处理的方法；对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线； （2）加深对截断效应的理解； （3）掌握使用MATLAB设计滤波器，并对语音信号处理的方法。对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第九组:汪涛张汉毅巫金敏张经中柳泽举 第六组：罗涛梁乐杰黄乃生 实验三 1、实验目的和要求 掌握采用MATLAB数字滤波器设计软件编制方法。软件要求在界面内有不同类型（高通低通带通带阻）滤波器的选择、或者只对低通滤波器采用不同方法设
数字信号处理复习题带答案
数字信号处理复习题带答案
1.若一模拟信号为带限信号，且对其抽样满足奈奎斯特条件，则只要将抽样信号通过_____A____即可完全不失真恢复原信号。 A、理想低通滤波器 B、理想高通滤波器 C、理想带通滤波器 D、理想带阻滤波器 2.下列哪一个单位抽样响应所表示的系统不是因果系统___D__? A、.h(n)=δ(n)+δ(n-10) B、h(n)=u(n) C、h(n)=u(n)-u(n-1) D、h(n)=u(n)-u(n+1) 3.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是_____A_____。 A.N≥M B.N≤M C.N≤2M D.N≥2M 4.以下对双线性变换的描述中不正确的是__D_________。 A.双线性变换是一种非线性变换 B.双线性变换可以用来进行数字频率与模拟频率间的变换 C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内 D.以上说法都不对 5、信号 3 (n)Acos(n) 78 x ππ =-是否为周期信号，若是周期信号，周期为多少？ A、周期N=3 7 π B、无法判断 C、非周期信号 D、周期N=14 6、用窗函数设计FIR滤波器时，下列说法正确的是___a____。
A 、加大窗函数的长度不能改变主瓣与旁瓣的相对比例。 B 、加大窗函数的长度可以增加主瓣与旁瓣的比例。 C 、加大窗函数的长度可以减少主瓣与旁瓣的比例 。 D 、以上说法都不对。 7.令||()n x n a =，01,a n <<-∞≤≤∞，()[()]X Z Z x n =，则()X Z 的收敛域 为 __________。 A 、1||a z a -<< B 、1||a z a -<< C 、||a z < D 、1||z a -< 。 8.N 点FFT 所需乘法（复数乘法）次数为 ____D___。 A 、2N log N B 、N C 、2N D 、 2 log 2 N N 9、δ(n)的z 变换是 A A. 1 B.δ(w) C. 2πδ(w) D. 2π 10、下列系统(其中y(n)是输出序列，x(n)是输入序列)中__ C___属于线性 系统。 A.y(n)=x 2(n) B.y(n)=4x(n)+6 C.y(n)=x(n-n 0) D.y(n)=e x(n) 11、在应用截止频率为Ωc 的归一化模拟滤波器的表格时，当实际Ωc ≠1时，代替表中的复变量s 的应为___B________。 A 、Ωc /s B 、s/Ωc C 、-Ωc /s D 、s/c Ω 12、用窗函数法设计FIR 数字滤波器时，在阶数相同的情况下，加矩形窗时所
数字信号处理习题解答1
第一章 第二章 11-=--m/2 m=-m -/2 12 m=--/2 -/21 2 m=-m=-()121.7DTFT[x(2n)]=(2n)e m=2n DTFT[x(2n)]=(m)e =[()(1) ()]e [()e e ()e ] [()()] j n n j m j m j m j m j m j j x x x m x m x m x m X e X e ωωωωπ ωωωπ∞ ∞∞ ∞∞ ∞∞ ∞ ∞ ∞-+-=+ =+∑∑ ∑∑∑，为偶数 求下列序列的傅里叶变换（）x(2n) 令，于是 -n 1 1 121 z (1) 2u(n)()2 ()2 1,|(2)|11(2),||n n n n n n X z u n z z z z z z z +∞ --=-∞+∞ --=-∞ --=== <-=>-∑∑14.求出下列序列的变换及收敛域 3.3(1).()cos(),781() 8 (2).()5.25n 640() (5)()x n A n A j n x n e x n y n e πππω=--==判断下面的序列是否周期的是常数 试判断系统是否为线性时不变的（）y(n)=x (n)(7) y(n)=x(n)sin() .试判断系统是否为因果稳定系统（）y(n)=x(n-n )
-1 -1-2 -1 -1112 1-317.X(z)=,2-5+2105< | z | < 2x(n)(2) | z | > 2x(n) 11 X(z)= -1-z 1-2z 05< | z | < 2(n)=2(-n-1)+()(n) | z | > 2(n)=()(n)-2(n)n n n n z z z u u u u 已知分别求：（）收敛域.对应的原序列收敛域对应的原序列解：收敛域.时： x 收敛域时： x -1-1 -1 -1-1 -1 21.(n)=0.9y(n-1)+x(n)+0.9x(n-1)(1)h(n)(2)H(e )1+0.9(1)H(z)=,|z|>0.91-0.91+0.9F(z)=H(z)z =z 1-0.9n 1z=0.9(n j n n z z z z h ω≥已知线性因果网络用下面差分方程表示： y 求网络的系统函数及单位脉冲响应写出网络频率响应函数的表达式，并定性画出其幅频特性曲线解： 令当时，有极点-1-1=0.9-112-1-1-1-1=0=0.9-1-1)=Res[F(z),0.9]1+0.9=z (z-0.9)|1-0.9=20.9(n)=0,n<0 n=0z =0,=0.9(n)=Res[F(z),0]+Res[F(z),0.9]1+0.91+0.9=z z|+z (z-0.9)|1-0.91-0.9=-1+2=1 h(n)=n z n z z z z z h z z z z ?∴因为系统是因果系统，所以有h 当时，有极点00000000=0n-m =0n -m =0 n n 20.9(n-1)+(n)+0.9 (2)H(e )=-0.9 (3)y(n)=h(n)*x(n) =(m)x(n-m) =(m)e =(m)e e =e H(e )+0.9=e -0.9 n j j j m j m j j m j j j j j u e e h h h e e ωω ω ωωωωωωωωδ∞ ∞ ∞ ?∑∑∑（ ）
数字信号处理实验十七
数字信号处理实验17 例题1： b=[8,-4,11,-2]; a=[1,-1.25,0.75,-0.125]; [sos,g]=tf2sos(b,a) [r,p,k]=residuez(b,a) 运行结果： sos = 1.0000 -0.1900 0 1.0000 -0.2500 0 1.0000 -0.3100 1.3161 1.0000 -1.0000 0.5000 g = 8 r =-8.0000 -12.0000i -8.0000 +12.0000i 8.0000 p = 0.5000 + 0.5000i 0.5000 - 0.5000i 0.2500 k = 16 例题2： sos=[1 0.9 0 1 -0.25 0;1 -3 2 1 1 0.5]; g=0.5; [b,a]=sos2tf(sos,g) [C,B,A]=dir2par(b,a) 子函数：dir2par（b，a）； function[C,B,A]=dir2par(num,den) M=length(num); N=length(den); [r1,p1,C]=residuez(num,den); p=cplxpair(p1,10000000*eps); I=cplxcomp(p1,p); r=r1(I); K=floor(N/2); B=zeros(K,2); A=zeros(K,3); if K*2==N; for i=1:2:N-2 Brow=r(i:1:i+1,:);
Arow=p(i:1:i+1,:); [Brow,Arow]=residuez(Brow,Arow,[]); B(fix((i+1)/2),:)=real(Brow); A(fix((i+1)/2),:)=real(Arow); end; [Brow,Arow]=residuez(r(N-1),p(N-1),[]); B(K,:)=[real(Brow),0]; A(K,:)=[real(Arow),0]; else for i=1:2:N-1 Brow=r(i:1:i+1,:); Arow=p(i:1:i+1,:); [Brow,Arow]=residuez(Brow,Arow,[]); B(fix(i+1)/2,:)=real(Brow); A(fix(i+1)/2,:)=real(Arow); end end 子函数：cplxcomp(p1,p2) function I=cplxcomp(p1,p2) I=[]; for j=1:length(p2) for i=1:length(p1) if(abs(p1(i)-p2(j))<0.0001) I=[I,i]; end; end end; I=I'; 运行结果： b = 0.5000 -1.0500 -0.3500 0.9000 a = 1.0000 0.7500 0.2500 -0.1250 C = -7.2000 B = 3.9846 1.6308 3.7154 0 A = 1.0000 1.0000 0.5000 1.0000 -0.2500 0 例题3： b=[8,-4,11,-2]; a=[1,-1.25,0.75,-0.125]; [K,C]=tf2latc(b,a) [b,a]=latc2tf(K,C) 运行结果： K = -0.7327 0.6032 -0.1250