信号系统实验报告

一,实验目的连续时间系统的单位冲激响应、单位阶跃响应在线性系统分析中的作用、地位及其MATLAB实现；握求解连续LTI系统响应的MATLAB实现方法；握连续时间信号的卷积在连续系统分析中的作用、卷积方法及其MATLAB的实现；二，实验原理设连续LTI(线性时不变)系统的激励为，响应为，则描述系统的微分方程可表示为(2-1)为了在MATLAB编程中调用有关函数，我们可以用向量和来表示该系统，即(2-2a)(2-2b)这里要注意，向量和的元素排列是按微分方程的微分阶次降幂排列，缺项要用0补齐。(2)系统的单位冲激响应单位冲激响应是指连续LTI系统在单位冲激信号激励下的零状态响应，因此满足线性常系数微分方程(2-1)及零起始状态，即，(2-3)按照定义，它也可表示为(2-4)对于连续LTI系统，若其输入信号为，冲激响应为，则其零状态响应为(2-5)可见，能够刻画和表征系统的固有特性，与何种激励无关。一旦知道了系统的冲激响应，就可求得系统对任何输入信号所产生输出响应。系统的单位阶跃响应单位阶跃响应是指连续LTI系统在单位阶跃信号激励下的零状态响应，它可以表示为(2-6)上式表明，连续LTI系统的单位阶跃响应是单位冲激响应的积分，系统的单位阶跃响应和系统的单位冲激响应之间有着确定的关系，因此，单位阶跃响应也能完全刻画和表征一个LTI系统。另外，对于二阶以上的高阶系统，系统的阶跃响应和冲激响应也能反映系统处于不同阻尼状态下的系统特性。连续信号的卷积积分两个连续时间信号和的卷积积分定义为(2-7)卷积是信号与系统分析的基本方法之一，是求系统零状态响应的重要工具。假设系统的输入信号为，单位冲激响应为，则系统的零状态响应为(2-8)此结果与式(2-5)完全相同，可见，信号卷积是分析系统特性的重要方法之一。关于卷积积分的物理含义、计算过程和计算结果的有关结论请参考[1]和[2]。三，实验内容下图所示为一RLC串联电路，已知R=5，L=1H，C=(1/6)F，请用MATLAB绘制出该系统的单位冲激响应和单位阶跃响应的波形，当输入信号时，请画出该系统的零状态响应波形图。vsvs(t)vc(t)RLCi(t)图2-1图2-1图2-1当电阻R分别为4、2、0.8、0.4时，观察它们的波形，并对波形的变化作出解释。已知某一连续LTI系统的单位冲激响应函数为，若系统的激励为，请用MATLAB的卷积方法，求该系统的零状态响应，并绘出的时域波形图，观察并说明不同的取样时间间隔p对波形的影响。课外练习设描述某连续时间系统的微分方程为试用MATLAB绘制出这个系统的冲激响应和阶跃响应波形。四，源代码及结果：figure(1)a=[1/6,5/6,1];b=[1];sys=tf(b,a);t=[0:0.1:10];subplot(3,1,1);impulse(sys,t);title('冲激响应')subplot(3,1,2);step(sys,t);title('阶越响应')u=0.5*exp(-2*t);subplot(3,1,3);lsim(sys,u,t);title('系统的零状态响应')figure(2)a1=[1/6,4/6,1];a2=[1/6,2/6,1];a3=[1/6,8/60,1];a4=[1/6,4/60,1];b=[1];sys1=tf(b,a1);sys2=tf(b,a2);sys3=tf(b,a3);sys4=tf(b,a4);t=[0:0.1:10];subplot(3,1,1);[y,t]=impulse(sys1,t);plot(t,y,'r')holdonsubplot(3,1,1);[y,t]=impulse(sys2,t);plot(t,y,'b')holdonsubplot(3,1,1);[y,t]=impulse(sys3,t);plot(t,y,'y')holdonsubplot(3,1,1);[y,t]=impulse(sys4,t);plot(t,y,'g')holdontitle('冲激响应图')subplot(3,1,2);[y,t]=step(sys1,t);plot(t,y,'r')holdonsubplot(3,1,2);[y,t]=step(sys2,t);plot(t,y,'b')holdonsubplot(3,1,2);[y,t]=step(sys3,t);plot(t,y,'y')holdonsubplot(3,1,2);[y,t]=step(sys4,t);plot(t,y,'g')holdontitle('阶跃响应图')u=0.5*exp(-2*t);subplot(3,1,3);[y,t]=lsim(sys1,u,t);;plot(t,y,'r')holdonsubplot(3,1,3);[y,t]=lsim(sys2,u,t);;plot(t,y,'b')holdonsubplot(3,1,3);[y,t]=lsim(sys3,u,t);;plot(t,y,'y')holdonsubplot(3,1,3);[y,t]=lsim(sys4,u,t);;plot(t,y,'g')holdontitle('系统的零状态响应图')figure(3)t=0:0.01:5;f=exp(-t);subplot(3,1,1);plot(f);h=t.*exp(-t)subplot(3,1,2);plot(h);y=conv(f,h);subplot(3,1,3);plot(y);结果分析：分别利用impulse(sys,t);step(sys,t);lsim(sys,u,t);得到了系统的冲激响应，阶跃响应和零状态响应。而阶跃响应是是在t=0开始作用，所以在t>0都有响应。而且可以看出冲激响应是阶跃响应的导数。五，实验体会由于长时间没