资源描述：
《实验四--基于MATLAB的控制系统时域分析以及根轨迹分析.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验四基于MATLAB的控制系统时域分析以及根轨迹分析一、实验目的掌握运用MATLAB进行控制系统的时域暂态及稳态分析方法以及根轨迹图的绘制方法。二、实验题目1.已知系统闭环传递函数为：(1)判定该系统的闭环稳定性；(2)绘出系统在阶跃信号为0.1´1(t)时的响应曲线，并求系统性能指标：稳态值、上升时间、调节时间、超调量；(3)绘制该系统的单位脉冲响应图；(4)绘制该系统的单位斜坡响应图，并分析该系统的稳态误差。设误差定义为输入量减去输出量；(5)绘制该系统在0.2t+sin(2t)作用下的响应曲线，设初始条件为零。sys=tf([20],[18
2、364020]);root=roots(sys.den{1});iffind(real(root)>0)disp('该系统不稳定')elsedisp('该系统稳定')endsubplot(2,2,1);%title('0.1*阶跃响应');t=linspace(0,10,200);u=0.1*ones(1,200);lsim(sys,u,t);%0.1倍阶跃响应[postrtstp]=stepchar(sys,0.02);%超调量，上升时间，调整时间，峰值时间subplot(2,2,2)title('单位脉冲响应');impulse(sys);%单
3、位脉冲响应subplot(2,2,3)title('斜坡响应');t=linspace(0,10,200);u=t;lsim(sys,u,t);%单位斜坡响应subplot(2,2,4)title('0.2*t+sin(2*t)');t=linspace(0,10,200);u=0.2*t+sin(2*t);lsim(sys,u,t);实验结果fori=1:dimtify(i)>=ysstr=t(i);%调节时间break;endendfori=1:length(y)ify(i)>=(1-delta)*yss
4、
5、y(i)<=(1+delta)*ys
6、sts=t(i);%调整时间endendendStepchar代码function[postrtstp]=stepchar(g0,delta)%超调量，上升时间，调整时间，峰值时间[y,t]=step(g0);[mp,ind]=max(y);%最大值dimt=length(t);yss=y(dimt);%稳态值pos=100*(mp-yss)/yss;%超调量tp=t(ind);%峰值时间1.一个系统的结构框图如下所示。（1）试绘出z=0.5，p=0.1时系统闭环的根轨迹图；并在根轨迹图上任选一点，试计算该点的增益K及其所有极点的位置；计算系统临界
7、稳定时的K值；说明该系统单位阶跃响应是否会出现振荡？与K的取值是否有关？clear;G=zpk([-0.5],[-0.10-2-10],[1]);%开环传递函数rlocus(G);%绘制根轨迹[x,y]=ginput(1);%任选一点pp=x+i*y;%p点的坐标k=rlocfind(G,p);%计算p点的增益值k[p1,k1]=rlocus(G,k);%计算增益值为k时的闭环极点值和增益值k2=rlocfind(G);%临界稳定的k值disp(k2);%临界稳定的k大概在200左右，即当增益值小于200时，系统稳定实验结果（2）绘制p=0.1，K
8、=20时该系统闭环极点随着z变化的根轨迹图；同时，绘制当z=0.5,1,5,15时系统的单位阶跃响应图，通过对比根轨迹图以及单位阶跃响应图的特性，进而说明该系统开环零点对系统暂态、稳态以及稳定性的影响。figure;subplot(3,2,1)G1=tf([20],[112.121.2220]);%随着零点Z变换的等效开环传递函数rlocus(G1);z=[000.51515];fori=3:6subplot(3,2,i);G2=zpk([-z(i)],[0-2-10-0.1],[20]);G22=feedback(G2,1);sys=1.5*G2
9、2;%闭环传递函数step(sys);end实验结果（3）绘制z=0.5，K=20时该系统闭环极点与随着p变化的根轨迹图；同时，绘制当p=0.3,1,5,15时系统的单位阶跃响应图，通过对比根轨迹图以及单位阶跃响应图的特性，进而说明该系统开环极点对系统暂态、稳态以及稳定性的影响。figure;subplot(3,2,1)G2=tf([112200],[112202010]);%随着极点p变换的等效开环传递函数rlocus(G2);p=[000.31515];fori=3:6subplot(3,2,i);G2=zpk([-0.5],[0-2-10-p
10、(i)],[20]);G22=feedback(G2,1);sys=1.5*G22;%闭环传递函数step(sys);en
