基于MATLAB的线性系统的时域分析上机仿真题.pdfVIP

基于 MATLAB 的线性系统的时域分析上机仿真题 一、目的 1.观察学习控制系统的时域（阶跃、脉冲、斜坡）响应； 2.记录时域响应曲线；给出时域指标。 3.掌握时域响应分析的一般方法。 二、实验内容 1．二阶系统为 10/ （s2+2s+10 ） ? 计算系统的闭环根、阻尼比、无阻尼振荡频率，并作记录。 输入： p=[1,2,10]; v=roots(p) wn=10^(0.5) e=2/2/wn 输出： v = -1.0000 + 3.0000i -1.0000 - 3.0000i wn = 3.1623 e = 0.3162 ? 记算实际测取的峰值大小 Cmax(tp) 、峰值时间 tp 、过渡时间 ts，并与理论值相比较。 输入： p=[1,2,10]; v=roots(p) wn=10^(0.5) e=2/2/wn 输出： v = -1.0000 + 3.0000i -1.0000 - 3.0000i wn = 3.1623 e = 0.3162 2 ．试作出以下系统的阶跃响应，并比较与原系统响应曲线的差别与特点，作出相应的 实验分析结果。 (a)G1(s)= （2s+1 ）/ （s2+2s+10），有系统零点情况。 (b)G2(s)= （s2+0.5 ）/ （s2+2s+10 ），分子、分母多项式阶数相等。 (c)G3(s)=s/ （s2+2s+10 ），分子多项式零次项系数为零。 输入： G1=tf([2,1],[1,2,10]); G2=tf([1,0,0.5],[1,2,10]); G3=tf([1,0],[1,2,10]); G0=tf([10],[1,2,10]); Y1=step(G1); Y2=step(G2); Y3=step(G3); Y0=step(G0); subplot(4,1,1); plot(Y0) subplot(4,1,2); plot(Y1) subplot(4,1,3); plot(Y2) subplot(4,1,4); plot(Y3) 3. 3 2 3s 16s 41s 6 5 4 3 2 s 14s 110s 528s 1494s 2117 s 112 的稳定性并计算稳态误差。 输入： a=[1,14,110,528,1494,2117,112]; b=[3,16,41,0]; k=my_stablity_lianxu(a,b) G=tf(a,b); sys=feedback(G,1); step(sys); t=[0:0.001:120]; y=step(sys,t); subplot(211),plot(t,y);grid subplot(212),ess=1-y ; plot(t,ess),grid ESS=ess(length(ess)) 输出： 稳定 k = 1 ESS = 3.1185e-005 分析： e1 时，过阻尼系统的阶跃响应其调节时间最长，进入稳态很慢； e=1 时，临界阻尼系统的阶跃响应没有超调量，且响应速度比过阻尼时要快； e1 时，无阻尼系统