数字信号处理第二章上机题作业
《数字信号处理》课后上机题#优选.
subplot(2,2,2);stem(n,sn1,'.')
title('(b)系统1的单位阶跃响应');
xlabel('n');ylabel('s(n)')
%系统2
xn=[1,zeros(1,30)];
%xn=单位脉冲序列,长度N=31
xi=filtic(B2,A2,ys);
实验报告
第一章:时域离散信号和时域离散系统
*16.已知两个系统的差分方程分别为
(1) y(n)=0.6y(n-1)-0.08y(n-2)+x(n)
(2) y(n)=0.7y(n-1)-0.1y(n-2)+2x(n)-x(n-2)
分别求出所描述的系统的单位脉冲响应和单位阶跃响应.
解:(可附程序)
(1)系统差分方程的系数向量为
subplot(2,2,1);stem(n,hn1,'.')
title('(a)系统1的系统单位脉冲响应');
xlabel('n');ylabel('h(n)')
xn=ones(1,30);
%xn=单位阶跃序列,长度N=31
sn1=filter(B1,A1,xn,xi);
%调用filter解差分方程,求系统输出信号sn1
%设差分方程(2)系数向量
%系统1
xn=[1,zeros(1,30)];
ys=0;
%xn=单位脉冲序列,长度N=31
xi=filtic(B1,A1,ys);
%由初始条件计算等效初始条件输入序列xi
hn1=filter(B1,A1,xn,xi);
《数字信号处理》(2-7章)习题解答
第二章习题解答1、求下列序列的z 变换()X z ，并标明收敛域，绘出()X z 的零极点图。
(1) 1()()2nu n (2) 1()()4nu n - (3) (0.5)(1)nu n --- (4) (1)n δ+(5) 1()[()(10)]2nu n u n -- (6) ,01na a <<解：(1) 00.5()0.50.5nn n n zZ u n z z ∞-=⎡⎤==⎣⎦-∑，收敛域为0.5z >，零极点图如题1解图(1)。
(2) ()()014()1414n nn n z Z u n z z ∞-=⎡⎤-=-=⎣⎦+∑，收敛域为14z >，零极点图如题1解图(2)。
(3) ()1(0.5)(1)0.50.5nnn n zZ u n z z --=-∞-⎡⎤---=-=⎣⎦+∑，收敛域为0.5z <，零极点图如题1解图(3)。
(4) [](1Z n z δ+=，收敛域为z <∞，零极点图如题1解图(4)。
(5) 由题可知，101010910109(0.5)[()(10)](0.5)()(0.5)(10)0.50.50.50.50.50.5(0.5)n n nZ u n u n Z u n Z u n z z z z z z z z z z z --⎡⎤⎡⎤⎡⎤--=--⎣⎦⎣⎦⎣⎦⋅=-----==--收敛域为0z >，零极点图如题1解图(5)。
(6) 由于()(1)nn n a a u n a u n -=+--那么，111()(1)()()()nn n Z a Z a u n Z a u n z z z a z a z a a z a z a ----⎡⎤⎡⎤⎡⎤=---⎣⎦⎣⎦⎣⎦=----=-- 收敛域为1a z a <<，零极点图如题1解图(6)。
(1) (2) (3)(4) (5) (6)题1解图2、求下列)(z X 的反变换。
数字信号处理 第二章习题
1 为因果序列，故收敛域为： z 2
8
（2） (n n0 ) n0 0
解：
X ( z)
n
x(n) z n
n
(n n0 ) z n
X ( z) z
n0
1 n n0 (n n0 ) 0 other
1 n0 z
z 0.5 左边序列 0.5 z 2 双边序列 右边序列 z 2
16
采用围线积分法求解：
3 2 X ( z) 1 1 0.5 z 1 2 z 1 3(1 2 z 1 ) 2(1 0.5 z 1 ) 5 7 z 1 1 1 (1 0.5 z )(1 2 z ) (1 0.5 z 1 )(1 2 z 1 )
z1 1, z2 2
X(z)的收敛域为
左边序列 z 1 1 z 2 双边序列 z 2 右边序列
24
F ( z) X ( z) z
n 1
z ( z 3) ( z 3) n 1 z zn ( z 1)( z 2) ( z 1)( z 2)
z 2
21
当收敛域为： z 2 0.5
1 n n 1 x(n) 3( ) u (n) 2 u (n 1) 2
22
收敛域为： z 2
右边序列
n 0 ，围线c内有2个1阶极点
x(n) Re s[( z 0.5) F ( z), 0.5] Re s[( z 2) F ( z), 2] ( z 0.5) 5z 7 zn ( z 0.5)( z 2) ( z 2)
双边序列
n 0 ，围线c内有1个1阶极点
数字信号处理课后习题答案(全)1-7章
第 1 章 时域离散信号和时域离散系统
（6） y(n)=x(n2)
令输入为
输出为
x(n－n0)
y′(n)=x((n－n0)2) y(n－n0)=x((n－n0)2)=y′(n) 故系统是非时变系统。 由于
T［ax1(n)+bx2(n)］=ax1(n2)+bx2(n2) =aT［x1(n)］+bT［x2(n)］
第 1 章 时域离散信号和时域离散系统
题4解图（一）
第 1 章 时域离散信号和时域离散系统
题4解图（二）
第 1 章 时域离散信号和时域离散系统
题4解图（三）
第 1 章 时域离散信号和时域离散系统
(4) 很容易证明： x(n)=x1(n)=xe(n)+xo(n)
上面等式说明实序列可以分解成偶对称序列和奇对称序列。 偶对称序列可 以用题中(2)的公式计算， 奇对称序列可以用题中(3)的公式计算。
（2） y(n)=x(n)+x(nN+1)k 0
（3） y(n)= x(k)
（4） y(n)=x(n－nn0)n0
（5） y(n)=ex(n)
k nn0
第 1 章 时域离散信号和时域离散系统
解：（1）只要N≥1， 该系统就是因果系统， 因为输出 只与n时刻的和n时刻以前的输入有关。
如果|x(n)|≤M， 则|y(n)|≤M， （2） 该系统是非因果系统， 因为n时间的输出还和n时间以 后（（n+1）时间）的输入有关。如果|x(n)|≤M, 则 |y(n)|≤|x(n)|+|x(n+1)|≤2M,
=2x(n)+x(n－1)+ x(n－2)
将x(n)的表示式代入上式， 得到 1 y(n)=－2δ(n+2)－δ(n+1)－0.5δ(2n)+2δ(n－1)+δ(n－2)
数字信号处理答案第二章
数字信号处理答案第⼆章第⼆章2.1 判断下列序列是否是周期序列。
若是，请确定它的最⼩周期。
（1）x(n)=Acos(685ππ+n ) （2）x(n)=)8(π-ne j（3）x(n)=Asin(343ππ+n )解 (1)对照正弦型序列的⼀般公式x(n)=Acos(?ω+n )，得出=ω85π。
因此5162=ωπ是有理数，所以是周期序列。
最⼩周期等于N=)5(16516取k k =。
（2）对照复指数序列的⼀般公式x(n)=exp[ωσj +]n,得出81=ω。
因此πωπ162=是⽆理数，所以不是周期序列。
（3）对照正弦型序列的⼀般公式x(n)=Acos(?ω+n )，⼜x(n)=Asin(343ππ+n )＝Acos(-2π343ππ-n )＝Acos(6143-n π)，得出=ω43π。
因此382=ωπ是有理数，所以是周期序列。
最⼩周期等于N=)3(838取k k =2.2在图2.2中，x(n)和h(n)分别是线性⾮移变系统的输⼊和单位取样响应。
计算并列的x(n)和h(n)的线性卷积以得到系统的输出y(n)，并画出y(n)的图形。
(a)1111(b)-1-1-1-1-1-1222222 3333 3444………nnn nnnx(n)x(n)x(n)h(n)h(n)h(n)21u(n)u(n)u(n)a n ===2 2knhkx)()(按照折叠、移位、相乘、相加、的作图⽅法，计算y(n)的每⼀个取样值。
(a) y(0)=x(O)h(0)=1y(l)=x(O)h(1)+x(1)h(O)=3y(n)=x(O)h(n)+x(1)h(n-1)+x(2)h(n-2)=4,n≥2(b) x(n)=2δ(n)-δ(n-1)h(n)=-δ(n)+2δ(n-1)+ δ(n-2)y(n)=-2δ(n)+5δ(n-1)= δ(n-3)(c) y(n)= ∑∞-∞=--kkn knuku a)()(=∑∞-∞=-aa n--+111u(n)2.3 计算线性线性卷积(1) y(n)=u(n)*u(n)(2) y(n)=λn u(n)*u(n)解：(1) y(n)= ∑∞-∞=-kknuku)(-)()(kknuku=(n+1),n≥0 即y(n)=(n+1)u(n) (2) y(n)=∑∞-∞=-kk knuku)()(λ=∑∞=-0)()(k kk n u k u λ=λy(n)=λλ--+111n u(n)2.4 图P2.4所⽰的是单位取样响应分别为h 1(n)和h 2(n)的两个线性⾮移变系统的级联，已知x(n)=u(n), h 1(n)=δ(n)-δ(n-4), h 2(n)=a n u(n),|a|<1,求系统的输出y(n).解ω(n)=x(n)*h 1(n) =∑∞-∞=k k u )([δ(n-k)-δ(n-k-4)]=u(n)-u(n-4)y(n)=ω(n)*h 2(n) =∑∞-∞=k kk u a )([u(n-k)-u(n-k-4)]=∑∞-=3n k ka,n ≥32.5 已知⼀个线性⾮移变系统的单位取样响应为h(n)=an-u(-n),0系统的单位阶跃响应。
(完整word版)数字信号处理答案第二章
第二章2.1 判断下列序列是否是周期序列。
若是,请确定它的最小周期.（1）x （n ）=Acos(685ππ+n ） （2)x （n)=)8(π-ne j（3）x （n)=Asin(343ππ+n ) 解 （1）对照正弦型序列的一般公式x （n ）=Acos （ϕω+n )，得出=ω85π。
因此5162=ωπ是有理数，所以是周期序列。
最小周期等于N=)5(16516取k k =。
（2）对照复指数序列的一般公式x （n ）=exp[ωσj +］n,得出81=ω。
因此πωπ162=是无理数,所以不是周期序列。
(3）对照正弦型序列的一般公式x （n)=Acos(ϕω+n ），又x （n)=Asin （343ππ+n ）＝Acos （-2π343ππ-n )＝Acos(6143-n π）,得出=ω43π.因此382=ωπ是有理数，所以是周期序列。
最小周期等于N=)3(838取k k =2.2在图2.2中，x （n ）和h(n)分别是线性非移变系统的输入和单位取样响应。
计算并列的x （n ）和h （n)的线性卷积以得到系统的输出y(n ）,并画出y(n)的图形。
(a)1111(b)(c)111110 0-1-1-1-1-1-1-1-1222222 33333444………nnn nnnx(n)x(n)x(n)h(n)h(n)h(n)21u(n)u(n)u(n)a n ===22解 利用线性卷积公式y(n ）=∑∞-∞=-k k n h k x )()(按照折叠、移位、相乘、相加、的作图方法,计算y(n)的每一个取样值。
（a ） y （0）=x （O)h （0）=1y （l ）=x （O ）h(1)+x （1)h （O)=3y （n)=x(O)h （n ）+x （1)h(n-1)+x(2)h （n —2)=4,n ≥2 （b) x(n ）=2δ(n ）-δ（n-1）h(n)=-δ（n)+2δ（n —1)+ δ(n —2）y(n ）=-2δ(n)+5δ（n —1）= δ(n-3） (c ） y （n ）=∑∞-∞=--k kn k n u k u a)()(=∑∞-∞=-k kn a=aa n --+111u （n ）2。
中国科学院刘艳老师现代数字信号处理第二章上机作业
一、上机作业要求
假设一个点目标在 x，y 平面上绕单位圆做圆周运动，由于外界干扰，其运 动轨迹发生了偏移。其中，x 方向的干扰为均值为 0，方差为 0.05 的高斯噪声； y 方向干扰为均值为 0，方差为 0.06 的高斯噪声。 1、产生满足要求的 x 方向和 y 方向随机噪声 500 个样本； 2、明确期望信号和观测信号； 3、试设计一 FIR 维纳滤波器，确定最佳传递函数：hopt Rxx 1Rxs ，并用该滤波器 处理观测信号，得到其最佳估计。 （注：自行设定误差判定阈值，根据阈值确 定滤波器的阶数或传递函数的长度） 。 4、要求 3 中，也可以选择 Kalman 滤波器进行滤波处理，采用哪种滤波器可以自 由选择。 5、分别绘制出 x 方向和 y 方向的期望信号、噪声信号、观测信号、滤波后信号、 误差信号的曲线图； 6、在同一幅图中绘制出期望信号、观测信号和滤波后点目标的运动轨迹。 7、实验报告要求： 给出求解思路和结果分析，给出 MATLAB 实现源程序和程序注 解。
四、实验结果及分析
图 1 x 方向上的期望信号、误差信号、观测信号、滤波后信号
图 2 y 方向上的期望信号、误差信号、观测信号、滤波后信号
ห้องสมุดไป่ตู้
图 3 滤波后的信号与原始信号的对比 图 3 中原始信号为红色， 滤波后信号为绿色。滤波后的结果与期望还是比较 接近的， 基本上实现了最优滤波功能。中心出现的点可能与滤波中的卷积运算有 关，卷积使信号长度增加，从而出现了中心点。
if
ey<1e-2 break;
数字信号处理 课后习题答案 第2章.docx
习题1.设X（e"。
）和r（e JC0）分别是印7）和）仞的傅里叶变换，试求下面序列的傅里叶变换:(1) x("-"o) (3) x(-n) (5) x(")y(")(7) x(2n)⑵ x*(〃)(4) x(") * v(«) (6) nx(n) (8) /(〃)解：⑴00 FT[X(/7-Z70)] = £x(〃一〃o)e—S令n r = n-n0,即〃=n' + n Q,贝!J00FT[x(n-n o y\=工》(〃')以"''*""="初。
乂(烈)00 00(2)FT[x («)] = £ x* (n)e*= [ £ 戏〃)攻以]* = X* (e「W=—00 w=—00(3)00FT[x(—")]= 〃)e*"令=一〃,则00FT[x(—”)]= Zx(〃')e" =X(e—〃")”'=—00(4)00 x(〃) *'(〃)= ^\x(jrT)y(n -m)W=-0000 00FT[x(n) * v(w)] = Z【Z x("y("-初)]e""' n=-<x> w=-oo k = n-m,贝U00 00FT[x(ri)*y(ri)]= £[ £x(初) k=—CD W=-0000 00k=-<x> m=—cc= X(e5(em)_00 00 1时[x(M)贝〃)]= Z》(〃)贝〃)e「9 = Zx(〃)[-Lf/(em'"'"d 渺]e-加""=—00 〃=—00 2l "1 00=—£ Y(e j0)')2l " n=—<x>1 伙=一L "口")*?®"、技或者FT[x{n)y{ny\ = —「171 »兀oo(6)因为X(e，")= »("初,对该式两边口求导，得到叫、)=-J £仗"如=-jFT[nx(n)]因此矶孙(〃)]=j至@3)dco00⑺ FT\x(2ri)\=加n=-(x)令n' = 2n ,则FT[X(2W)]= £x(z/)e 7 %W--00,且取偶数00 1 r r・l 八1°0 . 1 00 . 1£?kO + (T)“x(")厂=| 广伽+£ef ("广伽〃=—oo 匕匕〃=—oo 〃=—00=L「xa*+x(/*E)F7[x(2z?)] = | X(e‘2") + X(—e'尸)(8) F7[X2(»)]= J X2(77)6^»=-OO利用(5)题结果，令x{n) = y{n),则F巾2(”)] = _£x(em)*X(eS) = —「X®。
数字信号处理(第三版)第2章习题答案
第２章 时域离散信号和系统的频域分析
2.3
求信号与系统的频域特性要用傅里叶变换。 但分析频 率特性使用Z变换却更方便。 我们已经知道系统函数的极、 零点分布完全决定了系统的频率特性， 因此可以用分析极、 零点分布的方法分析系统的频率特性， 包括定性地画幅频 特性， 估计峰值频率或者谷值频率， 判定滤波器是高通、 低通等滤波特性， 以及设计简单的滤波器（内容在教材第5 章）等。
X e (e j ) FT[xr (n)]
Hale Waihona Puke 1 1 ej2 1 e j2 1 (1 cos 2)
24
4
2
因为 所以
Xe
(e j
)
1 2
[X
(e j
)
X
(e j
)]
X(ejω)=0π≤ω≤2π
X(e-jω)=X(ej(2π-ω))=0 0≤ω≤π
第２章 时域离散信号和系统的频域分析
当0≤ω≤π时，
用留数定理求其逆变换， 或者将z=ejω代入X(ejω)中， 得到X(z)函数， 再用求逆Z变换的方法求原序列。 注意收 敛域要取能包含单位圆的收敛域， 或者说封闭曲线c可取 单位圆。
第２章 时域离散信号和系统的频域分析
例如， 已知序列x(n)的傅里叶变换为
X
(e
j
)
1
1 ae
j
a 1
1 求其反变换x(n)。 将z=ejω代入X(ejω)中， 得到 X (z) 1 az 1
三种变换互有联系， 但又不同。 表征一个信号和系统 的频域特性是用傅里叶变换。 Z变换是傅里叶变换的一种推 广， 单位圆上的Z变换就是傅里叶变换。
第２章 时域离散信号和系统的频域分析
数字信号处理(方勇)第二章习题答案
2－1 试求如下序列的傅里叶变换： （1）)()(01n n n x -=δ （2）)1(21)()1(21)(2--++=n n n n x δδδ （3）),2()(3+=n u a n x n10<<a（4）)4()3()(4--+=n u n u n x（5）∑∞=-⎪⎭⎫⎝⎛=05)3(41)(k nk n n x δ（6）()6cos ,14()0,n n x n π⎧-≤≤=⎨⎩其他解： （1） 010()()j n j j nn X e n n ee ωωωδ∞--=-∞=-=∑（2） 2211()()122j j nj j n X e x n e e e ωωωω∞--=-∞==+-∑ωsin 1j +=（3） 2232()(2)1j j nj nn j nj n n a e X e a u n ea eaeωωωωω-∞∞---=-∞=-=+==-∑∑， 10<<a（4） []4()(3)(4)j j nn X e u n u n eωω∞-=-∞=+--∑∑-=-=33n nj e ω∑∑==-+=313n n j n nj e eωωωωωωωj j j j j e e e e e --+--=--111134=⎪⎭⎫ ⎝⎛⎪⎭⎫⎝⎛=----ωωωω21sin 27sin 1137j j j e ee（5） 3350011()(3)44n kj jn j k n k k X e n k e e ωωωδ∞∞+∞--=-∞==⎛⎫⎛⎫=-= ⎪ ⎪⎝⎭⎝⎭∑∑∑∑∞+=--⎪⎭⎫⎝⎛-=⎪⎭⎫⎝⎛=033411141k j kj e e ωω（6） 44336441()cos 32j j j jn jn n n X e nee e e ππωωωπ---=-=-⎛⎫==+ ⎪⎝⎭∑∑994()()4()()3333001122j j n j j n n n e e e e ππππωωωω--++===+∑∑ ()9()9334()4()33()()3311112211j j j j j j e e e e e e ππωωππωωππωω-+-+-+⎡⎤⎡⎤--⎢⎥⎢⎥=+⎢⎥⎢⎥++⎢⎥⎢⎥⎣⎦⎣⎦2－2 设信号}1,2,3,2,1{)(---=n x ，它的傅里叶变换为)(ωj e X ，试计算(1)0()j X e （２）()j X ed πωπω-⎰（３）2()j X e d πωπω-⎰。
《现代数字信号处理》第2章习题答案
∞
∞
1 1− z
1 2 −1
+
1 3 1 −1 = ⋅ 1 1 −1 1− 2 z 4 (1 − 2 z )(1 − 1 2 z)
−1 1 (1 − 1 3 1 3 1 2 z ) (1 − 2 z ) = ⋅ ⋅ ⋅ = ⋅ −1 1 −1 1 1 −1 1 1 4 (1 − 2 z )(1 − 2 z ) (1 − 3 z ) (1 − 3 z ) 4 (1 − 3 z )(1 − 1 3 z )
1 1− ∑ a (k ) z
k =1 2 v p
−k
2 2 , Px ( z ) =H ( z ) H * (1/ z * ) σ w =σw
1 1− ∑ a (k ) e
k =1 p
2
− jkω
(b) Pz ( z ) = Px ( z ) + σ
2.4 设给定一个线性移不变系统，其系统函数为 H ( z ) = (1 −
σ ∑⎢ ⎣
i =1
N
⎡
2 x
−
2 2 1 2⎤ σx + σx ⎥ N N ⎦
=
N −1 2 σx N
(b) E
{(σ
2
x
− E {σ x }
2
)}
2
⎧⎛ 2 N − 1 2 ⎞ 2 ⎪ ⎫ ⎧ N − 1 2 2 ( N − 1) 2 4 ⎫ ⎪ ˆx − = E ⎨⎜ σ σ x ⎟ ⎬ = E ⎨σ x4 − 2 σ xσ x + σx ⎬ 2 N N N ⎝ ⎠ ⎩ ⎭ ⎪ ⎪ ⎩ ⎭
{ }
N
( N − 1) 2 4 σx N2
− x)
（I）
数字信号处理（吴镇扬）课后习题答案（比较详细的解答过程）第二章测试训练题解
数字信号处理（吴镇扬）课后习题答案（比较详细的解答过程）第二章测试训练题解1.DFT和DTFT之间的关系是2.DFT和DFS之间的关系是3.对于一个128点的DFT，最先4个DFT相应于数字频率4.某滤波器的频响为H(ω) = 0.3cos2ω- 0.2cosω+ 0.05，相应于6点的DFT的H[k]为5.采样频率为22.05kHz的1024点DFT所对应的频率分辨率为6.采样率为8kHz的信号的256点DFT的第一个周期覆盖的频率范围是从0Hz至7.信号[ 1 0 2 ]的DFT每隔3个样点值重复，为8.以1600Hz对一220Hz的信号采样，进行64点DFT，最接近的DFT频率为9.以12kHz的信号对一4.25kHz的信号抽样，其256点DFT幅谱图的基带最大峰值点所对应的下标为10.采样频率为6kHz，1kHz信号的频率分辨率要达到50Hz，需11.采样频率为16kHz，1024点DFT的窗口长度为12.关于谱泄漏与窗口长度的关系是13.频谱图是展现信号的什么14.周期性方波的频谱图15.在FFT中的乘数因子是16.与512点的DFT相比，512点的FFT只需约几分之一的计算量17、一个长度为N的有限长序列可否用N个频域的采样值唯一地确定？18、计算两个Ｎ点序列的线性卷积，至少要做多少点的ＤＦＴ？19、x(2n)与x(n)的关系20、对于高斯序列x(n)=exp[-(n-p)2/q]，取16点作FFT，其幅度谱中低频分量最多的是21、一般地说按时间抽取基二FFT的_______序列是按位反转重新排列的。
22、信号x（n）＝sin(nπ/4) - cos(nπ/7)的数字周期为23、Ｎ＝２Ｌ点基二ＦＦＴ，共有______列蝶形，每列有____个蝶形。
24、信号s(t)=sin(4000πt)+sin(600πt)，则采样频率至少应为25、用按时间抽取法计算256点的FFT时，n=233的二进制位反转值是26、FFT之所以能减少DFT的运算量，是因为：，FFT减少DFT 运算量的基本处理思想是。
数字信号处理课后答案 第2章高西全
jω
n = −∞
∑
∞
x( n′)e − j2ωn′ = X (e j2ω )
| ω |< ω0
1, X (e ) = 0,
ω0 <| ω | ≤ π
求X(ejω)的傅里叶反变换x(n)。
解：
1 x ( n) = 2π
∫ωe
−
0
ω0
jωn
sin ω0 n dω = πn
3. 线性时不变系统的频率响应（频率响应函数） H(ejω)=|H(ejω)|ejθ(ω)， 如果单位脉冲响应h(n)为实序列， 试 证明输入x(n)=A cos(ω0n+ϕ)的稳态响应为
∑
1
∑
1
j ω j (ω − π ) 1 2 ) + X (e 2 = [ X (e )] 2
1
1
或者
FT[ x(2n)] =
（8）
1 j ω 1 [ X (e 2 )
2
+
1 j ω X (−e 2 )]
FT[ x 2 (n)] =
n = −∞
∑
∞
x 2 ( n ) e − jω n
利用（5）题结果， 令x(n)=y(n), 则
∑
∗
∞
x(n′)e − jω ( n + n0 ) = e − jωn0 X (e jω )
′
n = −∞
∑ x ( n )e
∞
− jωn
jωn = x ( n ) e = X ∗ ( e − jω ) n = −∞
∑
∞
∗
（3）
FT[ x( −n)] =
n = −∞
∑
《数字信号处理》课后上机题
title(’（a)系统1的系统单位脉冲响应’)；
xlabel('n');ylabel(’h(n)')
xn=ones（1，30)；
%xn=单位阶跃序列，长度N=31
sn1=filter（B1,A1，xn，xi)；
％调用filter解差分方程，求系统输出信号sn1
解：(可附程序)
Fs=80000；
T=1/Fs;
wp=2＊pi*4000/Fs;
ws=2＊pi＊20000/Fs；
rp=0。5；rs=45;
[N,wc］=buttord（wp/pi,ws/pi，rp,rs）
[B，A］=butter（N,wc)
clf；
mpplot(B,A,rs）;
function mpplot(B,A，rs）
if nargin<3 ymin=—80；else ymin=-rs—20；end；
[H,W]=freqz(B，A,1000）；
m=abs（H)；
p=angle(H)；
subplot(2，1,1)；
plot（W/pi，20*log10（m/max(m))）；grid on;
xlabel('\omega/\pi'）;ylabel（’幅度(dB)’）
（结果展示):
yn =
2 -1 —2 2 -2 —1 2
ycn =
2.0000 -1.0000 -2。0000 2。0000 -2。0000 —1.0000 2.0000
第四章:快速傅里叶变换(FFT)
＊6.按照下面的IDFT算法编写MATLAB语言IFFT程序，其中的FFT部分不用写出清单，可调用fft函数。并分别对单位脉冲序列，矩形序列，三角序列和正弦序列进行FFT和IFFT，验证所编程序.
数字信号处理答案第二章习题解答
————第二章————教材第二章习题解答1. 设()jw X e 和()jw Y e 分别是()x n 和()y n 的傅里叶变换，试求下面序列的傅里叶变换： （1）0()x n n -； （2）()x n -； （3）()()x n y n ； （4）(2)x n 。
解：（1）00[()]()jwnn FT x n n x n n e∞-=-∞-=-∑令''00,n n n n n n =-=+，则'00()'0[()]()()jw n n jwn jw n FT x n n x n e e X e ∞-+-=-∞-==∑（2）****[()]()[()]()jwnjwn jw n n FT x n x n ex n e X e -∞∞-=-∞=-∞===∑∑（3）[()]()jwnn FT x n x n e∞-=-∞-=-∑令'n n =-，则'''[()]()()jwn jw n FT x n x n eX e ∞-=-∞-==∑（4） [()*()]()()jwjwFT x n y n X e Y e = 证明： ()*()()()m x n y n x m y n m ∞=-∞=-∑[()*()][()()]jwnn m FT x n y n x m y n m e ∞∞-=-∞=-∞=-∑∑令k=n-m ，则[()*()][()()] ()() ()()jwk jwnk m jwkjwnk m jw jw FT x n y n x m y k eey k e x m eX e Y e ∞∞--=-∞=-∞∞∞--=-∞=-∞===∑∑∑∑2. 已知001,()0,jww w X e w w π⎧<⎪=⎨<≤⎪⎩求()jw X e 的傅里叶反变换()x n 。
解： 00sin 1()2w jwn w w nx n e dw nππ-==⎰3. 线性时不变系统的频率响应(传输函数)()()(),jw jw j w H e H e eθ=如果单位脉冲响应()h n 为实序列，试证明输入0()cos()x n A w n ϕ=+的稳态响应为00()()cos[()]jw y n A H e w n w ϕθ=++。
数字信号处理课后答案 第2章高西全
( −1) n x( n) = 2
n = −3
（4） 因为傅里叶变换的实部对应序列的共轭对称部分， 即
Re [ X (e jω )] =
n = −∞
∑
∞
x e ( n ) e − j ωn
1 xe (n) = ( x(n) + x(− n)) 2
按照上式画出xe(n)的波形如题5解图所示。
题5解图
2． 已知
jω
n = −∞
∑
∞
x( n′)e − j2ωn′ = X (e j2ω )
| ω |< ω0
1, X (e ) = 0,
ω0 <| ω | ≤ π
求X(ejω)的傅里叶反变换x(n)。
解：
1 x ( n) = 2π
∫ωe
−
0
ω0
jωn
sin ω0 n dω = πn
3. 线性时不变系统的频率响应（频率响应函数） H(ejω)=|H(ejω)|ejθ(ω)， 如果单位脉冲响应h(n)为实序列， 试 证明输入x(n)=A cos(ω0n+ϕ)的稳态响应为
∑
∗
∞
x(n′)e − jω ( n + n0 ) = e − jωn0 X (e jω )
′
n = −∞
∑ x ( n )e
数字信号处理第二章上机作业
第二章上机作业1、ljdt(A,B)函数定义function ljdt(A,B)p=roots(A);q=roots(B);p=p';q=q';x=max(abs([p q 1]));x=x+0.1;y=x;clfhold onaxis([-x x -y y])w=0:pi/300:2*pi;t=exp(i*w);plot(t)axis('square')plot([-x x],[0 0])plot([0 0],[-y y])text(0.1,x,'jIm[z]')text(y,1/10,'Re[z]')plot(real(p),imag(p),'x')plot(ral(q),imag(q),'o')title('pole-zero diagram for discrete system') hold off例2.26a=[3 -1 0 0 0 1];b=[1 1];ljdt(a,b)p=roots(a)q=roots(b)pa=abs(p)程序运行结果如下：P=0.7255+0.4633i0.7255+0.4633i-0.1861+0.7541i-0.1861-0.7541i-0.7455q=-1pa=0.86080.86080.77680.77680.7455例2.27b=[0 1 2 1];a=[1 -0.5 -0.005 0.3];subplot 311zplane(b,a);xlabel('实部');ylabel('虚部'); num=[0 1 2 1];den=[1 -0.5 -0.005 0.3];h=impz(num,den);subplot 312stem(h);xlabel('k');title('单位脉冲响应'); [H,w]=freqz(num,den);subplot 313plot(w/pi,abs(H));xlabel('频率\omega');title('频率响应')例2.28a=[1,-1];b=[1];subplot 321impz(b,a);a1=[1,-0.8];b1=[1];subplot 322impz(b1,a1,10);a2=[1,-2];b2=[1];subplot 323impz(b2,a2,10);a3=[1,-2*0.8*cos(pi/4),0.8^2];b3=[1];subplot 324impz(b3,a3,20);a4=[1,-2*0.8*cos(pi/8),1];b4=[1];subplot 325impz(b4,a4,20);a5=[1,-2*1.2*cos(pi/4),1.2^2];b5=[1];subplot 326impz(b5,a5,20);例2.29b=[1,0,-1];a=[1,0,-0.81];figure(1)subplot(2,1,1);dimpulse(b,a,50);ylabel('h(n)'); subplot(2,1,2);dstep(b,a,50);ylabel('g(n)'); figure(2)w=[0:1:500]*pi/500; freqz(b,a,w)例2.30b=[1,0,0,0,0,0,0,0,-1];a0=1;a1=[1,0,0,0,0,0,0,0,-(0.8)^8];a2=[1,0,0,0,0,0,0,0,-(0.9)^8];a3=[1,0,0,0,0,0,0,0,-(0.98)^8];[H,w]=freqz(b,a0);[H1,w1]=freqz(b,a1);[H2,w2]=freqz(b,a2);[H3,w3]=freqz(b,a3);subplot(4,2,1);zplane(b,a0);xlabel('实部');ylabel('虚部');title('FIR梳状滤波器零点图')subplot(4,2,2);zplane(b,a1);xlabel('实部');ylabel('虚部');title('IIR梳状滤波器零点图a=0.8')subplot(4,2,3);plot(w/pi,abs(H));title('FIR梳状滤波器幅频响应曲线') subplot(4,2,4);plot(w/pi,abs(H1));title('IIR梳状滤波器幅频响应曲线a=0.8')subplot(4,2,5);zplane(b,a2);xlabel('实部');ylabel('虚部');title('IIR梳状滤波器零极点图a=0.9')subplot(4,2,6);zplane(b,a3);xlabel('实部');yalbel('虚部')title('IIR梳状零极点图a=0.98')a4=[1,0,0,0,0,0,0,0,-(0.9)^8];a5=[1,0,0,0,0,0,0,0,-(0.98)^8];[H4,W4]=freqz(b,a4);[H5,W5]=freqz(b,a5);subplot(4,2,7);plot(w/pi,abs(H4));title('IIR梳状滤波器频幅响应曲线a=0.9')subplot(4,2,8);plot(w/pi,abs(H5));title('IIR梳状滤波器零极点图a=0.98')例2.31num=[0.45 0.4 -1];den=[1 -0.4 -0.45];x0=[1 2];y0=[0 1];N=50;n=[0:N-1]';x=0.8.^n;Zi=filtic(num,den,y0,x0);[y,Zf]=filter(num,den,x,Zi);plot(n,x,'r-',n,y,'b--');title('响应');xlabel('n');ylabel('x(n)-y(n)'); legend('输入x','输入y',1);grid;例2.32num=[18];den=[18 3 -4 -1];[r,p,k]=residuez(num,den) 程序运行结果如下：r =0.36000.24000.4000p =0.5000-0.3333-0.3333k =[]（1）f=sym('cos(a*k)');F=ztrans(f)程序运行结果如下：F =(z-cos(a))*z/(z^2-2*z*cos(a)+1) （2）F=sym('a*z/(z-a)^2');f=iztrans(F)程序运行结果如下：f =a^n*n例2.34(1)f=sym('a^n')F=ztrans(f)程序运行结果如下：f =a^nF =-z/(a - z)(2)f=sym('1');F=ztrans(f)程序运行结果如下F =z/(z - 1)b=1;a=poly([0.9 0.9 -0.9]);[r,p,k]=residuez(b,a)程序运行结果如下：r =0.25000.50000.2500p =0.90000.9000-0.9000k =[]例2.36x1=[1,2,3];n1=-1:1;x2=[2,4,3,5];n2=-2:1;[y,n]=conv_m(x1,n1,x2,n2);运行结果如下：y =2 8 17 23 19 15。
《数字信号处理》朱金秀第二章习题及参考答案
第二章 习题及参考答案 一、习题1、 序列x(n)的表达式如下：⎪⎩⎪⎨⎧≤≤-≤≤-+=其它,040,414,32)(n n n n x（1） 画出序列x(n)的波形，并标出各序列值。
（2） 请用延迟的单位抽样序列及其加权和表示序列x(n)。
（3） 令y(n)=2x(n-2),请画出y(n)的波形。
2、 判断下列序列是否是周期序列，并求出周期序列的周期。
（1）是常数αππα,)843sin()(-=n n x（2）是常数ββπ,)()16(-=n j en x3、设x(n)与y(n)分别表示系统的输入和输出，试判断下列差分方程所描述的系统是否是线性移不变的？ （1）y(n)=x(n)+5x(n-2) （2）y(n)=3x(n)+1（3）y(n)=x(n-n 0) , n 0为整数 （4）y(n)=x(-n) （5）y(n)=x 2(n) （6）y(n)=x(n 2) （7）∑==ni 0x(i)y(n)4、试判断下列差分方程所描述的系统是否具有因果性、稳定性，并说明理由。
（1）y(n)=x(n)+x(n+2) （2）∑+-==x(m)y(n)n n n n m（3）y(n)=x(n-m) （4）y(n)=e x(n) （5）∑-==20k)-x(n y(n)N k5、输入序列x(n)及线性时不变系统的单位脉冲响应h(n)如下所示：)3(2)1()2()(-+-++-=n n n n x δδδ )2(5.0)1()(2)(-+-+=n n n n h δδδ求该系统的输出序列y(n),并画出y(n)的波形。
6、设由下列差分方程描述的系统为因果系统， 3)1()(3)1()(-++-=n x n x n y n y 要求用递推法求系统的单位脉冲响应。
7、设)()(3n R n x =, 试求x(n)的共扼对称序列)(n x e 和共扼反对称序列)(n x o ，并分别用波形图表示。
8、根据系统的单位脉冲响应h(n),分析下列系统的因果稳定性： （1）)(0n n -δ （2）)(n u （3）)(3n u n （4）)(3n R N n（5）)(31n u n - （6）n n u )(（7）!)(n n u9、已知线性移不变系统的单位脉冲响应h(n)以及输入序列x(n)，求输出序列y(n),并画出y(n)的波形图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理作业实验题报告
第一章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)')
运行结果：
实验结论：
该系统稳定。