数字信号处理上机实验实验一：系统响应及系统稳定性院（系、部）：信息工程学院姓名：学号：班级：专业：电气工程及其自动化、实验目的(1)掌握求系统响应的方法。(2)掌握时域离散系统的时域特性。(3)分析、观察及检验系统的稳定性。二、实验仪器电脑、MATLAB软件三、实验容(1)编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。(2)给定一个低通滤波器的差分方程为y(n)=0.05x(n)+0.05x(n一1)+0.9y(n一1)输入信号x(n)=R(n)18x(n)=u(n)2a)分别求出系统对x(n)=R(n)和x(n)=u(n)的响应序列，并画出其182波形。b)求出系统的单位冲响应，画出其波形。(3)给定系统的单位脉冲响应为h(n)=R(n)110求出系统的输出响应，并画出其波形。stem(y1n);title('(b)系统对R8(n)的响应yl(n)');boxony2n=filter(B,A,x2n);subplot(2,2,3);y=y2n;stem(y2n);title('(c)系统对u(n)的响应y2(n)');boxon容2：调用conv函数计算卷积程序：xln=[llllllll];hln=[ones(l,l0)zeros(l,l0)];h2n=[l2.52.5lzeros(l,l0)];y2ln=conv(hln,xln);y22n=conv(h2n,xln);figure(2)subplot(2,2,l);y=hln;stem(hln,y);title('(d)系统单位脉冲响应hl(n)');boxonsubplot(2,2,2);y=y21n;stem(y21n,y);title('(e)hl(n)与R8(n)的卷积y21(n)');boxonsubplot(2,2,3);y=h2n;stem(h2n,y);title('(f)系统单位脉冲响应h2(n)');boxonsubplot(2,2,4);y=y22n;stem(y22n,y);title('(g)h2h2(n)与R8(n)的卷积y22(n)');boxon容3：谐振器分析un=ones(1,256);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];y31n=filter(B,A,un);y32n=filter(B,A,xsin);figure(3)求出系统的输出响应，并画出其波形。subplot(2,1,1);y=y31n;stem(y31n,y);title('(h)谐振器对u(n)的响应y31(n)');boxonsubplot(2,1,2);y=y32n;stem(y32n,y);title('(i)谐振器对正弦信号的响应y32(n)');boxon五、实验程序运行结果及分析程序运行结果如图10.1.1所示。实验容(1)系统的单位冲响应、系统对x(n)二R(n)和x(n)=u(n)的响应序列182分别如图(a)、(b)和(c)所示；实验容(2)系统h](n)和h2(n)对x(n)二R(n)的输出响应分别如图(e)和(g)所示；实验容(3)系统对u(n)和x(n)=sin(0.014n)+sin(0.4n)的响应序列分别如图(h)和(i)所示。由图(h)可见，系统对u(n)的响应逐渐衰减到零，所以系统稳定。由图⑴可见，系统对x(n)=sin(0.014n)+sin(0.4n)的稳态响应近似为正弦序列sin(0.4n)，这一结论验证了该系统的谐振频率是0.4rad。伺系统单位脉冲响应h^)鮒系统对R88n的响酬M)⑹系统对对u(n的nn)图1容一输出波形0.50.5(dd橐统单位脉艸响应w)0o1(i)鰭谐振器对正弦信号的响应y32n)
