数字信号处理第二章上机题作业

《数字信号处理第二章上机题作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数字信号处理作业实验题报告第一章16.(1)实验目的：求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。实验要求：运川matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。源程序：Bl=l;Al=[l,-0.6,0.08J;ys=2；%设差分方程xn=[l,zeros(l,20)];%xn二单位脉冲序列，长度N二31xi=filtic(Bl,Al,ys);hn1=filter(B1,Al,xn,xi);%求系统输出信号hnln=0:lcngth(hnl)-l;subplot(2,I,

2、1);stem(n,hnl,7)titl*单位脉冲响应J;xlabel('n');ylabel('h(n)^)xn=ones(1,20);sn1=filter(B1,Al,xn,xi);%求系统输出信号snln二0:length(snl)-l;Subplot(2,1,2);stemm,snl,titleC单位阶跃响应J;xlabel('n');ylabel('s(n)')2-ys=2;3%设差分方程4-xn=[1,zeros(lj20)];5二单位脉冲序列，长度6-xi=filtic(BljAl,ys);7-hnl=filter(Bl,Al,xn

3、3xi);8%求系绒输出信号hnl9-n=0:length(hnl)-1;10-subplot(2,lj1);11-stem(rijhnl,?)12-titleC单位脉冲响应’)；13-xlabel(n);ylabelCh(n)')14-xn=ones(1,20);15-snl=filter(Bl,Al,xn,xi);16%求系统输出信号srd17-n=0:length(snl)-l;18-subplot(2,1,2);19-stem(rijsnl,?.7)20-titleC单位阶跃响应’)；21-xlabel('n');22-ylabel('

4、s(n)*)单位阶跃响应运行结果:%Editor-C:MATLAB7wFileEditTextCellToolsDebugDesktopWindowHelpD盲V為(S厂c爭關FPl)l?l电悄Stack:1Base▼!QFigure1-n1-B1=1;A1=[1,-0.6,0.08];FileEditViewInsertToolsDesktopWindowHelp□HsffeH紗®庚[口匡1叵回单位脉冲响应3r<2工.I・01_1_1_1_*_A_*_*_A_A_*_A_A_A_A_«_*_A_*05101520n实验分析:单位脉冲响M逐渐

5、趋于0,阶跃响应保持不变，由此可见，是个稳定系统。第二章31题实验目的：用matlab判断系统是否稳定。实验要求：用matlabM出系统的极，零点分布图，输入单位阶跃序列u(n)检杳系统是否稳定。源程序：A=[2,298,0.17,2.341&-1.5147J;B=[0,0,1,5,-50];subplot(2,l,l);zplanc(B,A);%求H⑵的极点p=roots(A);%求H⑵的模pm=abs(p);ifmax(pm)

6、filter(B,A,un);n=0:length(sn)-l;subplot(2,1,2);plot(n,sn)xlabel('n');ylabel('s(n)')运行结果:FileEditTextCellToolsDebugDesktopW□3船豎®gc昌阿化Figure1FileEditViewInsertToolsDesktopWindowHelp□Qa題Q礬⑥晏□S□Q-2:,I:,-10-5051-A二[2,-2.98,0.17,2.3418,-1.5147];2-B=[0,0,1,5,-50];3-subplot(2,1,1);4-

7、zplane(B,A);5滋求H(z)的极点6-p=roots(A);1%求班2：)的模8-pjn=abs(p);9-ifmax(pm)<1disp「系统因果稳定')』10-else,dispC系统因果不稳定'),end11-un=ones(1,800);12-sn=filter(B,A,un);13-n=0:length(sn)-l;14-subplot(2,1,2);plot(rijsn)15-xlabel(5n);ylabelCs(n)?)2:IIy"r/•x09考V?3实验结论:该系统稳定。