吉大(2021-2022)学期《数字信号处理》在线作业二答案3

根据零、 极点分布可定性画幅频特性。 当频率由0到2π 变化时， 观察零点矢量长度和极点矢量长度的变化， 在极点 附近会形成峰。 极点愈靠进单位圆， 峰值愈高； 零点附近形 成谷， 零点愈靠进单位圆， 谷值愈低， 零点在单位圆上则 形成幅频特性的零点。 当然， 峰值频率就在最靠近单位圆的 极点附近， 谷值频率就在最靠近单位圆的零点附近。
故
X (z)z 1zN z 1 N (z 1 1 )zN z 1 N (z 1 1 )z2 1 N 1 zz N 1 1 2
［例2.4.4］ 时域离散线性非移变系统的系统函数H(z)为
H(z) 1 , a和b为常数 (za)(zb)
（1） 要求系统稳定， 确定a和b的取值域。 （2） 要求系统因果稳定， 重复（1）。 解： （1） H(z)的极点为a、 b， 系统稳定的条件是收敛 域包含单位圆， 即单位圆上不能有极点。 因此， 只要满足 |a|≠1, |b|≠1即可使系统稳定， 或者说a和b的取值域为除单位圆 以的整个z平面。 （2） 系统因果稳定的条件是所有极点全在单位圆内， 所以a和b
采样间隔T=0.25 s， 得到 xˆ ( t ) ， 再让 xˆ ( t ) 通过理想低通
滤波器G(jΩ)， G(jΩ)用下式表示:
G(j)0.025
≤ 4π 4π
(1) 写出xˆ ( t )的表达式;
(2) 求出理想低通滤波器的输出信号y(t)。
解：(1)
x ˆ(t) [c2 o πn s)T (co 5πs n()T ](tn)T n
（3） 若y(n)=x(n)h(n), 则
Y(ej)1H(ej)X(ej) 2π
这是频域卷积定理或者称复卷积定理。
（4）
xe(n)12[x(n)x(n)]

数字信号处理（姚天任江太辉）第三版课后习题答案第二章2.1 判断下列序列是否是周期序列。

若是，请确定它的最小周期。

（1）x(n)=Acos(685ππ+n ) （2）x(n)=)8(π-nej（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)的图形。

解 利用线性卷积公式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.3 计算线性线性卷积 (1) y(n)=u(n)*u(n) (2) y(n)=λnu(n)*u(n)解：(1) y(n)=∑∞-∞=-k k n u k u )()(=∑∞=-0)()(k k n u k u =(n+1),n ≥0即y(n)=(n+1)u(n) (2) y(n)=∑∞-∞=-k k k n u k u )()(λ=∑∞=-0)()(k kk n u k u λ=λλ--+111n ,n ≥0即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 nu(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)=a n-u(-n),0<a<1 用直接计算线性卷积的方法，求系统的单位阶跃响应。

数字信号处理习题及答案数字信号处理作业(1)1、画出离散信号的波形（1）)2(3)3(2)(1++-=n n n x δδ（2）)2()(2+-=n u n x（3）)5()()(3--=n u n u n x（4）)()()(214n u n x n ?= （5）)()25.0sin(3)(5n u n n x ??=π2、设x (n )、y (n )分别为系统的输⼊、输出变量，根据定义确定系统是否为：（1）线性，（2）稳定，（2）因果① )()]([ )(2n ax n x T n y == ② b n x n x T n y +==)()]([ )( ③ )0( )()]([ )(00>-==n n n x n x T n y④ ∑+-=>=0)0( )( )(0n n n n m n m x n y3、已知:描述系统的差分⽅程为)()1(5- )(n x n y n y =-且初始条件为： 0)1(=-y 求：系统的单位冲激响应h (n )4、已知：线性时不变系统的单位脉冲响应为10 , )( )(<求：该系统的单位阶跃响应。

数字信号处理作业(1)解答1、画出离散信号的波形（1）)2(3)3(2)(1++-=n n n x δδ（2）)2()(2+-=n u n x（3）)5()()(3--=n u n u n x（4）)()()(214n u n x n ?= （5）)()25.0sin(3)(5n u n n x ??=π2、设x (n )、y (n )分别为系统的输⼊、输出变量，根据定义确定系统是否为：（1）线性，（2）稳定，（3）因果因果：输出只取决于当前和之前的输⼊。

线性移不变系统的因果的充要条件：h (n )＝0 ， n < 0稳定系统：有界输⼊产⽣有界输出。

线性移不变系统稳定的充要条件：∞<=∑∞-∞=n n n x n x T n y （线性，稳定，因果）④ )0( )( )(0>=∑+-=n m x n y n n n n m （线性，稳定，⾮因果）注意：⾮线性系统的稳定、因果只能按定义判断，不能按线性、移不变系统的h (n )特点判断。

数字信号处理实验课后答案数字信号处理实验课后答案【篇一：数字信号处理第三版课后实验程序(高西全)】txt>close all;clear all%======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== a=[1,-0.9];b=[0.05,0.05]; %系统差分方程系数向量b和ax1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=r8(n)x2n=ones(1,128); %产生信号x2(n)=u(n)hn=impz(b,a,58); %求系统单位脉冲响应h(n)subplot(2,2,1);y=h(n);tstem(hn,y); %调用函数tstem绘图title((a) 系统单位脉冲响应h(n));box ony1n=filter(b,a,x1n); %求系统对x1(n)的响应y1(n)subplot(2,2,2);y=y1(n);tstem(y1n,y);title((b) 系统对r8(n)的响应y1(n));box ony2n=filter(b,a,x2n); %求系统对x2(n)的响应y2(n)subplot(2,2,4);y=y2(n);tstem(y2n,y);title((c) 系统对u(n)的响应y2(n));box on%===内容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=r8(n)h1n=[ones(1,10) zeros(1,10)];h2n=[1 2.5 2.5 1 zeros(1,10)];y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)subplot(2,2,1);y=h1(n);tstem(h1n,y); %调用函数tstem绘图title((d) 系统单位脉冲响应h1(n));box onsubplot(2,2,2);y=y21(n);tstem(y21n,y);title((e) h1(n)与r8(n)的卷积y21(n));box onsubplot(2,2,3);y=h2(n);tstem(h2n,y); %调用函数tstem绘图title((f) 系统单位脉冲响应h2(n));box onsubplot(2,2,4);y=y22(n);tstem(y22n,y);title((g) h2(n)与r8(n)的卷积y22(n));box on%=========内容3：谐振器分析======================== un=ones(1,256); %产生信号u(n)n=0:255;xsin=sin(0.014*n)+sin(0.4*n); %产生正弦信号a=[1,-1.8237,0.9801];b=[1/100.49,0,-1/100.49]; %系统差分方程系数向量b和a y31n=filter(b,a,un); %谐振器对u(n)的响应y31(n) y32n=filter(b,a,xsin); %谐振器对u(n)的响应y31(n)figure(3)subplot(2,1,1);y=y31(n);tstem(y31n,y);title((h) 谐振器对u(n)的响应y31(n));box onsubplot(2,1,2);y=y32(n);tstem(y32n,y);title((i) 谐振器对正弦信号的响应y32(n));box on10.2.2 实验程序清单1 时域采样理论的验证程序清单% 时域采样理论验证程序exp2a.mtp=64/1000; %观察时间tp=64微秒%产生m长采样序列x(n)% fs=1000;t=1/fs;fs=1000;t=1/fs;m=tp*fs;n=0:m-1;a=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=a*exp(-alph*n*t).*sin(omega*n*t);xk=t*fft(xnt,m); %m点fft[xnt)]yn=xa(nt);subplot(3,2,1);tstem(xnt,yn); %调用自编绘图函数tstem绘制序列图box on;title((a) fs=1000hz);k=0:m-1;fk=k/tp;subplot(3,2,2);plot(fk,abs(xk));title((a) t*ft[xa(nt)],fs=1000hz);xlabel(f(hz));ylabel(幅度);axis([0,fs,0,1.2*max(abs(xk))])%===================================== ============ % fs=300hz和fs=200hz的程序与上面fs=1000hz完全相同。