资源描述：
《实验3 根轨迹分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验3根轨迹分析一、实验目的1、学习和掌握利用MATLAB绘制根轨迹图的方法。2、学习和掌握利用系统根轨迹图分析系统的性能。二、实验原理根轨迹分析的MATLAB实现根轨迹是指系统某一参数变化时，闭环特征根在s平面上运动的轨迹。在MATLAB中，提供了用于根轨迹分析的专门函数。1）rlocus函数该函数的使用方法如下：rlocus(sys)绘制单输入单输出LTI系统的根轨迹图。rlocus(sys,k)使用用户指定的根轨迹增益k来绘制系统的根轨迹图。[r,k]=rlocus(sys)返回根轨迹增益值和闭环极点值，不绘制根轨迹图。2）rlocfind函数该函数
2、的使用方法如下：[k,poles]=rlocfind(sys)计算鼠标选取点处的根轨迹增益值和闭环极点值，可在图形窗口根轨迹图中显示出十字光标，当用户选择其中一点时，相应的增益值和极点值记录在k和poles中。[k,poles]=rlocfind(sys,p)计算最靠近给定闭环极点p处的根轨迹增益。3）sgrid函数该函数的使用方法如下：sgrid可在连续系统根轨迹或零极点图上绘制出栅格线，栅格线由等阻尼系数和等自然频率线构成。sgrid(‘new’)先清除当前的图形，然后会指出栅格线，并将坐标轴属性设置成holdon。sgrid(z,Wn)指定阻尼系数z
3、和自然频率Wn。sgrid(z,Wn,’new’)指定阻尼系数z和自然频率Wn，在绘制栅格线之前清除当前的图形并将坐标轴属性设置成holdon。三、实验内容1、已知系统开环传递函数为GS=K(s+5)s+1s+3(s+12)（1）试用MATLAB绘制系统的根轨迹图。（2）求根轨迹的两条分支离开实轴时的K值，并确定该K值对应的所有闭环极点。（3）以区间[-40,-5]之间的值替代s=-12处的极点，重新绘制根轨迹图，观察其对根轨迹图的影响。实验代码及结果：(1)>>z=[-5];>>p=[-1-3-12];>>k=[1];>>sys1=zpk(z,p,k)s
4、ys1=(s+5)------------------(s+1)(s+3)(s+12)Continuous-timezero/pole/gainmodel.>>rlocus(sys1)(2)可知此时增益K为3.38。[a,pole]=rlocfind(sys)Selectapointinthegraphicswindowselected_point=-2.1213+0.0000ia=3.3812pole=-11.7575-2.1213+0.0000i-2.1213-0.0000i(3)>>z=[-5];>>p=[-1-3-40];>>k=[1];>>sys
5、1=zpk(z,p,k)sys1=(s+5)------------------(s+1)(s+3)(s+40)Continuous-timezero/pole/gainmodel.>>rlocus(sys1)>>p=[-1-3-30];>>p=[-1-3-20];>>p=[-1-3-10];>>p=[-1-3-5];2、已知系统开环传递函数为Gs=K(s+8)ss+2(s2+8s+32)（1）试用MATLAB绘制系统的根轨迹图。（2）计算两条分支进入右半平面和两条分支复数极点出发在实轴相交处的K值。（3）以区间[-20,-1]之间的值替代零点的位置，重新
6、绘制根轨迹图，观察其对根轨迹图的影响。实验代码及结果：(1)>>z1=[-8];>>p1=[0-2];>>k1=[1];>>sys1=zpk(z1,p1,k1)sys1=(s+8)-------s(s+2)Continuous-timezero/pole/gainmodel.>>num2=[1];>>den2=[1832];>>sys2=tf(num2,den2)sys2=1--------------s^2+8s+32Continuous-timetransferfunction.>>sys=sys1*sys2sys=(s+8)-------------
7、----------s(s+2)(s^2+8s+32)Continuous-timezero/pole/gainmodel.>>rlocus(sys)(2)进入右半平面K=53.7在实轴相交出的K=2070(3)>>[z,p,k]=zpkdata(sys,'v')z=-8p=0-2.0000-4.0000+4.0000i-4.0000-4.0000ik=1>>z=[-20];>>gzpk=zpk(z,p,k)gzpk=(s+20)-----------------------s(s+2)(s^2+8s+32)Continuous-timezero/pole
8、/gainmodel.>>rlocus(gzpk)>>z=[-10
