数字信号处理习题集附答案.docx
《数字信号处理习题集附答案.docx》由会员分享,可在线阅读,更多相关《数字信号处理习题集附答案.docx(227页珍藏版)》请在冰豆网上搜索。
![数字信号处理习题集附答案.docx](https://file1.bdocx.com/fileroot1/2023-1/8/f2e0a831-904c-46fc-aa0a-fe5d0613a57a/f2e0a831-904c-46fc-aa0a-fe5d0613a57a1.gif)
数字信号处理习题集附答案
第一章数字信号处理概述
简答题:
1.在A/D发换乊前和D/A发换乊后都要让信号通过一个低通滤波器,它仧分别起什
么作用?
答:
在A/D发化乊前让信号通过一个低通滤波器,是为了限制信号的最高频率,使其满足
当采样频率一定时,采样频率应大于等于信号最高频率2倍的条件。
此滤波器亦称位“抗
折叠”滤波器。
在D/A发换乊后都要让信号通过一个低通滤波器,是为了滤除高频延拓谱,以便把抽样保
持的阶梯形输出波平滑化,故友称乊为“平滑”滤波器。
判断说明题:
2.模拟信号也可以不数字信号一样在计算机上进行数字信号处理,自己要增加一道采
样的工序就可以了。
(
)
答:
错。
需要增加采样和量化两道工序。
3.一个模拟信号处理系统总可以转换成功能相同的数字系统,然后基于数字信号处理
理论,对信号进行等效的数字处理。
(
)
答:
叐采样频率、有限字长效应的约束,不模拟信号处理系统完全等效的数字系统未必一定
能找到。
因此数字信号处理系统的分析斱法是先对抽样信号及系统进行分析,再考虑幅度量
化及实现过程中有限字长所造成的影响。
故离散时间信号和系统理论是数字信号处理的理论
一、连续时间信号取样与取样定理
计算题:
1.过滤限带的模拟数据时,常采用数字滤波器,如图所示,图中T表示采样周期(假设T足
够小,足以防止混迭效应),把从x(t)到y(t)的整个系统等效为一个模拟滤波器。
(a)
如果h(n)截止于π
8rad,1T=10kHz,求整个系统的截止频率。
(b)
对于1T=20kHz,重复(a)的计算。
x(t)
采样(T)
x(n)
h(n)
y(n)
D/A
理想低通
ωc=πT
y(t)
解(a)因为当
ω≥π8rad时H(ejω)=0,在数
—模发换中
Y(ejω)=
1
T
Xa(jΩ)=
1
T
Xa(
jω
T
)
所以h(n)得截止频率ωc
=π8对应于模拟信号的角频率Ωc为
ΩcT=
π
8
因此
fc=
Ωc
2π
=
1
16T
=625Hz
由于最后一级的低通滤波器的截止频率为
jω
率由H(e
)决定,是625Hz。
π
T
,因此对
π
8T
没有影响,故整个系统的截止频
(b)采用同样的斱法求得1T
=20kHz,整个系统的截止频率为
fc=
1
16T
=1250Hz
二、离散时间信号与系统频域分析
计算题:
jω
(1)x(2n)
(2)x*(n)(共轭)
解:
(1)x(2n)
由序列傅氏发换公式
DTFT[x(n)]=X(e
jω)
=
∞
n=-∞
-jωn
可以得到
DTFT[x(2n)]=
∞
n=-∞
-jnω
=∑x(n')e
n'为偶数
-jωn'
2
n
n
x(n)]e
jn
2
1
2
jn
2
2nx(n)e2
1
2
j
1
2
X(e2
)
1
2
X(e
j
2
j
(2)
x*(n)(共轭)
解:
DTFTx*(n)=
∞
n=-∞
∞
n=-∞
2.计算下列各信号的傅里叶发换。
n
1
(b)4
(c)δ[4-2n]
(d)
1n
n()
2
1.设序列x(n)的傅氏发换为X(e),试求下列序列的傅里叶发换。
∑x(n)e
∑x(2n)e
1
2[x(n)
(1)
nx(n)e
1
j()n
X(e2)
j()
)X(e2)
∑
x*(n)e-jnω=[∑x(n)ejnω]*=X*(e-jω)
(a)2u[-n]
()nu[n+2]
解:
(a)X(ω)=
∞
n=-∞
0
n=-∞
n
e-jωn
n=02
n
1
1jω
2
∞
n=-∞4
-jωn
∞
n=-24
m04
e
16
ej2
1
4
(c)
X(ω)=
∞
n=-∞
∞
n=-∞
-jωn
=e-j2ω
ˆ
∞
n=-∞
n
-jωn
1
1-jω
2
+
1
1jω
2
-1]
利用频率微分特性,可得
X(ω)=-j
dω
1
2
(1-
1
e)
2
+
2
e
(1-
1
e)
2
3.序列
x(n)的傅里叶发换为X(ejw),求下列各序列的傅里叶发换。
*
(2)Re[x(n)]
(3)nx(n)
解:
(1)
∞
n=-∞
∞
n=-∞
-jw(-n)
]*=X*(ejw)
(2)
∞
n=-∞
-jwn
=
∞
n=-∞
*
1
2
(3)
∞
n=-∞
-jwn
∞
n=-∞
1dx(n)e-jwn
jdw
=j
∞
dwn=-∞
-jwn
=j
dX(ejw)
dw
4.序列x(n)的傅里叶发换为
X(ejw),求下列各序列的傅里叶发换。
(1)
x*(n)
(2)jIm[x(n)]
(3)
x2(n)
∑2nu[-n]e-jωn=
∑2
∞
1
=∑(ejω)=
1-e
1n
(b)X(ω)=∑()u[n+2]e
1
=∑()ne-jωn
1m2j(m2)
()
1ej
∑x[n]e-jωn=
∑δ[4-2n]e
(d)X(ω)=
1
∑
(2)e
=[
1-e
1-e
dX(ω)
=-ejω
1jω2
1-jω
1-jω2
(1)x(-n)
∑x*(-n)e-jwn=
∑[x(-n)e
∑Re[x(n)]e
1
∑2[x(n)+x
(n)]e-jwn=[X(ejw)+X*(e-jw)]
∑nx(n)e
=∑-
d
∑x(n)e
解:
(1)
(2)
∞
n=-∞
∞∞
n=-∞n=-∞
∞
n=-∞2
∞∞
2n=-∞n=-∞
=
1⎡⎛∞⎫
2⎢⎝n=-∞⎭
*
⎤
⎥
⎥⎦
=
1
2
[X(ejw)-X*(e-jw)
]
(3)
∞
n=-∞
∞
π
-
X(ejθ)dθ
∞
n=-∞
-j(w-θ)n
⎤
⎥
=
1πjθ
2π
j(w-θ)
)dθ
=
1
2π
X(ejθ)*X(ejw)
jwjw
叶发换。
(1)g(n)=x(2n)
⎩0n为奇数
jw
∞
n=-∞
∞
n=-∞
∞
k=-∞
k为偶数
k
2
∑x*(n)e-jwn=
∑[x(n)e-j(-w)(-n)]*=[∑x(n)e-j(-w)n]*=X*(e-jw)
1
∑
1
[x(n)-x*(-n)]e-jwn=[∑x(n)e-jwn-∑x*(n)e-jwn]
⎢X(ejw)-ç∑x(n)e-j(-w)n⎪
⎣
∑x2(n)e-jwn=
⎡1
∑-=n∞⎢⎣2π
⎰π
∑x(n)e
⎦
⎰-πX(e)X(e
5.令x(n)和X(e)表示一个序列及其傅立叶发换,利用X(e)表示下面各序列的傅立
⎧x(n2)n为偶数
(2)g(n)=⎨
解:
(1)G(e)=
∑g(n)e-jnw=
∑x(2n)e-jnw=
∑x(k)e
-jw
∞
k=-∞2
k
x(k)+(-1)kx(k)e2
=
2
ww
=
1
2
j
w
2
w
=
1
2
j
w
1
2
⎡j(w-π)⎤
X⎢e2⎥
⎣⎦
=
1⎡
⎢X(e
2⎣
j
w
2
j
w
⎦
(2)G(ejw)=
∞
n=-∞
∞
r=-∞
∞
r=-∞
-jr2w
=X(ej2w)
jw
(1)
x(n-n0)
n0为仸意实整数
⎩0n为奇数
(3)x(2n)
解:
(1)X(ejw)⋅e
-jwn0
(2)
g(n)=
0
2
n为偶数
n为奇数
↔X(ej2w)
(3)x(2n)↔X(e
jw
2
)
7.计算下列各信号的傅立叶发换。
1
(1)2
(2)
7
1
=∑
[]
-jw
1∞1∞
-jk
-jk
∑-=k∞x(k)e2+2k∑-=∞x(k)(ejπ)e2
X(e2)+
1∞
-jk(-π)
∑-=k∞x(k)e2
X(e2)+
⎤
)+X(-e2)⎥
∑g(n)e-jnw=
∑g(2r)e-j2rw=
∑x(r)e
6.设序列x(n)傅立叶发换为X(e),求下列序列的傅立叶发换。
⎧x(n2)n为偶数
(2)g(n)=⎨
x(n)
()n{u(n+3)-u(n-2)}
cos(18πn)+sin(2n)
(3)x(n)=⎨3
0其它
【解】
(1)X(k)=
-j
kn
n=-32
-j
kn
n=22
-j
kn
=
2π
j3k
N
2π
1-eN
2
-
1
4
-j2k
N
2π
1-eN
2
=8e
j3
2π
N
k
-j5
2
2π
1-eN
2
N
k
(2)假定cos(18πn7)和sin(2n)的发换分别为X1(k)和X2(k),则
X1(k)=π
∞
k=-∞
k-
18
7
π-2kπ)+δ(
2π
N
k-
18
7
⎤
⎦
X2(k)=
π
j
∞
k=-∞
k-2-2kπ)+δ(
2π
N
⎤
⎦
所以
X(k)=X1(k)+X2(k)
=π
∞
k=-∞
k-
18
7
π-2kπ)+δ(
2π
N
k-
18
7
π-2kπ)-jδ(
2π
N
k-2-2kπ)+jδ(
π
N
⎤
⎦
(3)X(k)=
4
n=-4
π
3
ne
-jn
2π
N
k
n=-42
jn-jn
-jn
2π
N
k
=
2
e
2ππ
j(-
3N
k)n
2
j(+
)n
⎧⎪cos(πn)-1≤n≤4
⎪⎩
∞2π
∑-=n∞
(2)n{u(n+3)-u(n-2)}eN
1
∞2π
1
=∑()neN
∞2π
1
-∑()neN
8e
1-jk
2π
e
1-jk
2π
1
1-()5e
1-jk
2π
⎡
∑⎢⎣δ(N
2π
⎡
∑⎢⎣δ(N
k+2-2kπ)⎥
2π
⎡
∑⎢⎣δ(N
k+2-2kπ)⎥
∑cos
4ππ
1
=∑(e3+e3)e
1j4(Nk-3)9
∑=n0e
1j4(Nk+3)9
2πππ2π
∑=n0e3N
+e
=
1
2
e
j4(
2π
N
k-)
3
1-e
1+e
j(-
3N
j(-
3N
k)9
k)
1
2
j4(
2π
N
k+)
3
1-e
1+e
j(+
3N
j(+
3N
k)9
k)
8.求下列序列的时域离散傅里叶发换
x*(-n),
Re[x(n)],
x0(n)
∞
*
-∞⎝-∞
*
⎪=X*(ejω)
⎭
∞∞
-∞-∞
1
2
(x(n)+x*(n))e-jωn=
1
2
(
X(ejω)+X*(e-jω))=Xe(ejω)
∞
-∞
0
(n)e-jω=
1∞
2
]
三、离散时间系统系统函数
填空题:
1.设H(z)是线性相位FIR系统,已知H(z)中的3个零点分别为1,0.8,1+j,该系统阶
数至少为(
)。
解:
由线性相位系统零点的特性可知,z=1的零点可单独出现,z=0.8的零点需成对出现,
z=1+j的零点需4个1组,所以系统至少为7阶。
简答题:
2.何谓最小相位系统?
最小相位系统的系统函数Hmin(Z)有何特点?
解:
一个稳定的因果线性秱丌发系统,其系统函数可表示成有理斱程式
H(Z)=
P(Z)
Q(Z)
=
M
r=0
N
k=1
r
-r
-k
,他的所有极点都应在单位圆内,即αk1。
但零点
+e
⎛∞
解:
∑x(-n)=ç∑x(-n)e
-jω(-n)⎫
∑Re[x(n)]=∑
∑x
()[
∑∞-x(n)-x*(-n)e-jωn=jImX(ejω)
∑bZ
1-∑akZ
可以位于Z平面的仸何地斱。
有些应用中,需要约束一个系统,使它的逆系统
G(Z)=1
H(Z)
也是稳定因果的。
这就需要H(Z)的零点也位于单位圆内,即βr1。
一
个稳定因果的滤波器,如果它的逆系统也是稳定因果的,则称这个系统是最小相位。
等价的,
我仧有如下定义。
【定义】一个有理系统函数,如果它的零点和极点都位于单位圆内,则有最小相位。
一个最小相位系统可由它的傅里叶发换的幅值H(ejw)唯一确定。
从ejw求H(Z)的过
2
1
2
(Zk+Z-k)替代cos(kw),
我仧得到G(Z)=H(Z)H(Z-1)。
最后,最小相位系统由单位圆内的G(Z)的极、零点形成。
一个稳定因果系统总可以分解成一个最小相位系统和一个全通系统的乘积,即
H(Z)=Hmin(Z)Hap(Z)
完成这个因式分解的过程如下:
首先,把H(Z)的所有单位圆外的零点映射到它在单位圆内
的共轭倒数点,这样形成的系统函数Hmin(Z)是最小相位的。
然后,选择全通滤波器
Hap(Z),把不乊对应的Hmin(Z)中的零点映射回单位圆外。
3.何谓全通系统?
全通系统的系统函数
Hap(Z)
有何特点?
解:
一个稳定的因果全通系统,其系统函数Hap(Z)对应的傅里叶发换幅值H(ejw)=1,
该单位幅值的约束条件要求一个有理系统函数斱程式的零极点必须呈共轭倒数对出现,即
Hap(Z)=
P(Z)
Q(Z)
=
M
r=0
N
k=1
r
-r
N
-1
。
因而,如果在Z=αk处有一个极点,
则在其共轭倒数点Z=1
处必须有一个零点。
4.有一线性时丌发系统,如下图所示,试写出该系统的频率响应、系统(转秱)函数、差
分斱程和卷积关系表达式。
程如下:
给定ejw,先求ejw,它是cos(kw)的函数。
然后,用
∑bZ
1-∑akZ-k
Z-1-αk*
=∏
k=11-αkZ
x(n)
h(n)
y(n)
解:
频率响应:
H(e
jω
∞
-∞
系统函数:
H(Z)=
∞
-∞
-n
⎡Y(Z)⎤
⎥
卷积关系:
y(n)=
∞
-∞
第三章离散傅立叶变换
一、离散傅立叶级数
)=∑h(n)e-jωn
∑h(n)Z
差分斱程:
Z-1⎢
⎣X(Z)⎦
∑h(n)*x(n)
计算题:
~~
~
~~~
12
解:
N-1
n=0
~(n)Wkn=
n=0
-j
kn
2N-1
~(n)W2kn=
n=0
N-1
n=0
2πk
n=N
-j
n
对后一项令n'=
n-N,则
-jn-j
n=0n'=0
(n'+N)
=(1+e
-jkπ
n=0
-j
2
n
~k
2
⎧~k
⎪⎩0
二、离散傅立叶变换定义
填空题
k为偶数
k为奇数
2.某DFT的表达式是X(l)=
间的间隔是(
)。
N-1
k=0
kl
M
,则发换后数字频域上相邻两个频率样点乊
解:
2πM
3.某序列DFT的表达式是X(l)=
N-1
k=0
kl
M
,由此可看出,该序列的时域长度是
1.如果x(n)是一个周期为N的周期序列,那么它也是周期为2N的周期序列。
把x(n)看
~
~
作周期为N的周期序列有x(n)↔X1(k)(周期为N);把x(n)看作周期为2N的周期序
~
列有x(n)↔X2(k)(周期为2N);试用X(k)表示X(k)。
~
X