数字信号系统课后上机题答案.doc

课后习题第一章 17*. 已知系统的差分方程为 　　　　y(n)=－a1y(n－1)－a2y(n－2)+bx(n) 其中, a1=－0.8, a2=0.64, b=0.866。 （1） 编写求解系统单位脉冲响应h(n)(0≤n≤49)的程序， 并画出h(n)(0≤n≤49)； （2） 编写求解系统零状态单位阶跃响应s(n)(0≤n≤100）的程序， 并画出s(n)(0≤n≤100）。 %调用filter解差分方程,求系统单位脉冲响应和单位阶跃响 B=0.866;A=[1,-0.8,0.64];%查分方程向量 %(1)求解系统单位脉冲响应,并画出h(n) xn=[1,zeros(1,48)]; %xn=单位脉冲序列,长度N=31 hn=filter(B,A,xn);%调用filter解差分方程,求系统输出信号hn n=0:length(hn)-1; subplot(3,2,1);stem(n,hn,.) title( (a)系统的单位脉冲响应);xlabel(n);ylabel(h(n)) %(2)求解系统单位阶跃响应,并画出h(n) xn=ones(1,100);%xn=单位阶跃序列,长度N=100 sn=filter(B,A,xn);%调用filter解差分方程,求系统单位阶跃响应sn n=0:length(sn)-1; subplot(3,2,2);stem(n,sn, .);axis([0,30,0,2]) title( (b)系统的单位阶跃响应);xlabel(b);ylabel(s(n)) 第二章 31*. 假设系统函数如下式： (1) 画出极、零点分布图， 并判断系统是否稳定； (2) 用输入单位阶跃序列u(n)检查系统是否稳定。 解：（1）%程序ex231.m %判断系统的稳定性 A=[2,-2.98,0.17,2.3418,-1.5147]; %H(z)的分母多项式系数 B=[0,0,1,5,-50]; %H(z)的分子多项式系数用极点分布判断系统是否稳定 subplot(2,1,1); zplane(B,A); %绘制H(z)的零极点图 p=roots(A); %求H(z)的极点 pm=abs(p); %求H(z)的极点的模 if max(pm)1 disp(系统因果稳定),else disp(系统不因果稳定); end %画出u(n)的系统输出波形进行判断 un=ones(1,700); sn=filter(B,A,un); n=0:length(sn)-1; subplot(2,1,2); plot(n,sn) xlabel(n);ylabel(s(n)) 系统因果稳定 第三章 25*.已知序列h(n)=R6(n),x(n)=nR8(n）. (1)计算⑧x(n); (2)计算 x(n)和y(n)=h(n)*x(n); (3)画出h(n)、x(n)、和y(n)的波形图，观察总结循环卷积与线性卷积的关系。 解：%程序ex325.m hn=[1 1 1 1];xn=[0 1 2 3]; %用DFT计算4点循环卷积yc4n: H4k=fft(hn,4); %计算h(n)的4点DFT X4k=fft(xn,4); %计算x(n)的4点DFT Yc4k=H4k.*X4k; yc4n=ifft(Yc4k,4); %用DFT计算8点循环卷积yc8n: H8k=fft(hn,8); %计算h(n)的8点DFT X8k=fft(xn,8); %计算x(n)的8点DFT Yc8k=H8k.*X8k; yc8n=ifft(Yc8k,8); yn=conv(hn,xn); %时域计算线性卷积yn 第四章 19*.假设滤波器的系统函数为 在单位圆上采样六点，选择r=0.95,试画出它的频率采样结构，并在计算机上用DFT求出频率采样结构中的有关系数。 解：= = 式中，分母分子多项式各有一个零点z=1，相互抵消，因此系统仍然稳定，属于FIR系统。由系统函数得到单位脉冲响应为 h(n)=5δ(n)+5δ(n-1)+5δ(n-2)+3δ(n-3)+3δ(n-4)+3δ（n-5） H(k)=DFT[h(n)] k=0,1,2,3, …,5 H(z)=(1- %程序ex519.m hn=[5,5,5,3,3,3];r=0.95; Hk=fft(hn,6); for k=1:3, hk(k)=Hk(k);Wk(k)=exp(-j*2*pi*(k-1)/6); end H0=Hk(1) H3=Hk(4) r0k=2*real(hk) rlk=-2*real(r*hk.*Wk) 运行结果： H0 =2