数字信号处理第三版课后答案讲义.ppt

　　hn1=filter(B1， A1， xn， xi)； 　　　%调用filter解差分方程， 求系统输出信号hn1 　　n=0： length(hn1)-1；  　　subplot(3， 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 　　n=0： length(sn1)－1；  　　subplot(3， 2， 2)； stem(n， sn1， ′.′) 　　itle(′(b) 系统1的单位阶跃响应′)； xlabel(′n′)； 　　 　　ylabel(′s(n)′) 　　%====================================== 　　%系统2 　　xn=［1， zeros(1， 30)］； 　　　　%xn=单位脉冲序列， 长度N=31 　　xi=filtic(B2， A2， ys)； 　　　　%由初始条件计算等效初始条件输入序列xi 　　hn2=filter(B2， A2， xn， xi)； 　　　　%调用filter解差分方程， 求系统输出信号hn2 　　n=0： length(hn2)－1；  　　subplot(3， 2， 5)； stem(n， hn2， ′.′) 　　title(′(a) 系统2的系统单位脉冲响应′)； 　 xlabel(′n′)； ylabel(′h(n)′) 　　xn=ones(1， 30)； 　　　　%xn=单位阶跃序列， 长度N=31 　　sn2=filter(B2， A2， xn， xi)； 　　　　%调用filter解差分方程， 求系统输出信号sn2 　　n=0： length(sn2)－1；  　　subplot(3， 2， 6)； stem(n， sn2， ′.′) 　　title(′(b) 系统2的单位阶跃响应′)； xlabel(′n′)； 　　ylabel(′s(n)′) 　　程序运行结果如题16*解图所示。 题16*解图 　　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）。 　　解： 调用MATLAB函数filter计算该系统的系统响应的程序ex117.m如下：  　　%程序ex117.m 　　%调用filter解差分方程， 求系统单位脉冲响应和单位阶跃响应 　　B=0.866； A=［1， －0.8， 0.64］；%差分方程系数向量 　　%====================================== 　　%（1）求解系统单位脉冲响应， 并画出h(n)  　　xn=［1， zeros(1， 48)］； 　　　　%xn=单位脉冲序列， 长度N=31 　　hn=filter(B1， A1， 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(′n′)；