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

数字信号处理作业实验题报告第一章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=31xi=filtic(B1,A1,ys);hn1=filter(B1,A1,xn,xi);%求系统输出信号hn1n=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);%求系统输出信号sn1n=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);ifmax(pm)1disp('系统因果稳定'),else,disp('系统因果不稳定'),endun=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)')运行结果：实验结论：该系统稳定。