实验1 时域分析.docVIP

实验一 时域分析 一、 实验目的： 1、掌握一阶系统和二阶系统的非周期信号响应。 2、理解二阶系统的无阻尼、欠阻尼、临界阻尼和过阻尼响应。 3、掌握分析系统的稳定性、瞬态过程和稳态误差。 4、理解高阶系统的主导极点对系统特性的影响。 5、理解系统的零点对系统动态特性的影响。 二、 实验设备 PC机及MATLAB平台 三、实验原理及方法 1、系统的单位阶跃响应 用下列指令step(num,den)或step(num,den,t),就可求取系统的单位阶跃响应。前者指令中虽没有时间t的出现，但时间矢量会自动生成；后者指令中的t是由用户确定的时间。响应曲线图的x轴和y轴坐标也是自动标注的。 例如：已知一系统的闭环传递函数为，试求该系统的单位阶跃响应。 解： num=[16]; den=[1 4 16]; step(num,den); grid on; xlabel(’t’);ylabel(’c(t)’); title(’Unit-Step Response of G(s)=16/s^2+4s+16’); 若step指令等号左端含有变量时，即[y,x,t]=step(num,den,t)。使用该命令时，屏幕上不显示该系统的输出响应曲线。若要获得系统的响应曲线，则需加上plot绘图指令。 对上例： num=[16]; den=[1 4 16]; t=0:0.1:10; [y,x,t]=step(num,den,t); plot(t,y); grid on; xlabel(t);ylabel(c(t)); title(Unit-Step Response of G(s)=16/s^2+4s+16); 例：二阶系统闭环传递函数的标准形式为，令为一定值，则系统的瞬态响应只与参变量有关。试分析分别为0、0.3、0.5、0.7、1时系统的单位阶跃响应。 %matlab程序 t=0:0.1:12;num=[1]; Zeta1=0;den1=[1 2*Zeta1 1 ]; Zeta2=0.3;den2=[1 2*Zeta2 1 ]; Zeta3=0.5;den3=[1 2*Zeta3 1 ]; Zeta4=0.7;den4=[1 2*Zeta4 1 ]; Zeta5=2;den5=[1 2*Zeta5 1 ]; [y1,x,t]=step(num,den1,t); [y2,x,t]=step(num,den2,t); [y3,x,t]=step(num,den3,t); [y4,x,t]=step(num,den4,t); [y5,x,t]=step(num,den5,t); plot(t,y1,t,y2,t,y3,t,y4,t,y5); grid on; 2、判别系统的稳定性 利用Matlab求根指令，可方便地判定系统的稳定性。 roots(den) 例：已知某系统闭环特征方程式为，试判别系统稳定性。 %matlab程序 den=[1 3 2 24]; roots(den) 运行结果 ans = -4.0000 0.5000 + 2.3979i 0.5000 - 2.3979i 显然该系统不稳定。 3、系统的单位斜坡响应和单位脉冲响应 impulse(G) %绘制系统的脉冲响应 用下列指令step(num,den)或step(num,den,t),就可求取系统的单位阶跃响应。前者指令中虽没有时间t的出现，但时间矢量会自动生成；后者指令中的t是由用户确定的时间。响应曲线图的x轴和y轴坐标也是自动标注的。 例如：已知一系统的闭环传递函数为，试求该系统的单位脉冲阶跃响应。 解： num=[16]; den=[1 4 16]; impulse(num,den); grid on; xlabel(t);ylabel(c(t)); title(Unit-impulse Response of G(s)=16/s^2+4s+16); 对应更一般输入的响应也可以得到。命令为lsim ，其命令格式为： y = lsim ( num , den , u , t ) ; 输入信号为矢量u 。输入信号u 的行数决定了计算的输出点数。对于单输入系统，u 是一个列矢量。对于多输入系统，u 的列数等于输入变量数。 例如：计算斜坡响应，t 为输入矢量。 ramp = t ;y = lsim ( num , den , ramp , t ) ; 例如：已知一系统的闭环传递函数为，试求该系统的单位斜坡响应。 num=[16]; den=[1 4 16]; t = [0:0.1:10] ; ramp = t ; y = lsim ( num , den ,