1第三节控制系统的根轨迹分析法2利用根轨迹,可以对闭环系统的性能进行分析和校正由给定参数确定闭环系统的零极点的位置;分析参数变化对系统稳定性的影响;分析系统的瞬态和稳态性能;根据性能要求确定系统的参数;对系统进行校正。3一、条件稳定系统的分析[例4-11]:设开环系统传递函数为:)14.1)(6)(4()42()(22sssssssksGgk试绘制根轨迹并讨论使闭环系统稳定时的区值范围。gk开环极点:0,-4,-6,,零点:714.07.0j732.11j实轴上根轨迹区间:]0,4[),6,(渐进线:与实轴的交点:13.3324.164mnzpii,3)12(mnk倾角:[解]根据绘制根轨迹的步骤,可得:0464分离角(点):2d22)(,42)('2ssNsssNssssssD246.43394.11)(2345242.871176.455)(234'sssssD3.9497.4579.3758.805.97131.6280-4-3.5-3-2.5-2.0-1.5-1-0.50sgdk的最大值为9.375,这时s=-2.5,是近似分离点。gdk由:dsgdsNsDksDsNsDsN|)()(0)()()()(''''可以求得分离点。近似求法:分离点在[-4,0]之间。5入射角:1032与虚轴的交点(略)。这时的增益值:195,64,14gpk64gk64gk14gk14gk195gk195gk由图可知:当和时,系统是稳定的(为什么?);当时,系统是不稳定的。140gk19564gk6414195ggkk和左图是用Matlab工具绘制的。6条件稳定系统:参数在一定的范围内取值才能使系统稳定,这样的系统叫做条件稳定系统。具有正反馈的环节。下面的系统就是条件稳定系统的例子:开环非最小相位系统,其闭环系统的根轨迹必然有一部分在s的右半平面;7[例]非最小相位系统:,试确定使系统稳定时的增益值。)2)(1()(ssksGgk[解]:根轨迹如右:有闭环极点在右半平面,系统是不稳定的。显然稳定临界点在原点。该点的增益临界值为。gpk闭环特征方程为:,当s=0时,,所以,系统稳定的条件是:022gkss2gpk2gk8二、瞬态性能分析和开环系统参数的确定利用根轨迹可以清楚的看到开环根轨迹增益或其他开环系统参数变化时,闭环系统极点位置及其瞬态性能的改变情况。以二阶系统为例:开环传递函数为)2()(2sssGnk闭环传递函数为2222)(nnnsss共轭极点为:nnjs22,11在s平面上的分布如右图:nnj21闭环极点的张角为:1222cos,)()1(cosnnn所以称为阻尼角。斜线称为等阻尼线。9我们知道闭环二阶系统的主要的性能指标是超调量和调整时间。这些性能指标和闭环极点的关系如下:%100%100%21ctgee)(33为极点实部nst0306090%80604020的关系如下图和%若闭环极点落在下图中红线包围的区域中,有:3%sctgte和nnj2110上述结论也可应用于具有主导极点的高阶系统中。如下例:[例4-12]单位反馈系统的开环传递函数为:若要求闭环单位阶跃响应的最大超调量,试确定开环放大系数。)6)(4()(sssksGgk%18%[解]:首先画出根轨迹如右。由图可以看出:根轨迹与虚轴的交点为+j5,-j5,这时的临界增益当时,闭环系统不稳定。240gpk240gkAB11下面计算超调量和阻尼角的关系。由于:%,100%1ctge当时解得:%18%60这是一个三阶系统,从根轨迹上看出,随着的增加,主导极点越显著。所以可以用二阶系统的性能指标近似计算。gk在根轨迹图上画两条与实轴夹角为的直线,与根轨迹交与A、B两点。则A、B两点就是闭环共轭主导极点,这时系统的超调量为18%。通过求A、B两点的坐标,可以确定这时的根轨迹增益,进而求得开环放大系数k。60gk046A123j设A点坐标为:,则:j360tg(1)相角条件为:3216412011tgtg(2)12由(1),(2)式解得:1.2,2.1共轭主导极点为:。1.22.12,1js计算对应的根轨迹增益。由幅值条件:1)6)(4(1.22.1jsgsssk解得:44gknjjmiiksTsksG11)1()1()(开环传递函数以的形式表示时,k称为开环放大系数。显然的关系为:,式中不计0极点。gkk与jigpzkkjp所以,开环放大系数:83.16444k由于闭环极点之和等于开环极点之和,所以另一个闭环极点为:6.73p13[特别提示]:开环零、极点对根轨迹形状的影响是值得注意的。一般说,开环传递函数在s左半平面增加一个极点将使原根轨迹右移。从而降低系统的相对稳定性,增加系统的调整时间。)1()(ssksGgk)2)(1()(sssksGgk14)3)(2)(1()(ssssksGg15若在开环传递函数中增加一个零点,则原根轨迹向左移动。从而增加系统的稳定性,减小系统响应的调整时间。)1(2)(ssssGk)1()12)(12()(ssjsjssGk16Matlab参考书推荐:现代控制工程,[美]KatsuhikoOgats,卢伯英译,电子工业出版社MATLAB控制系统设计,欧阳黎明著,国防工业出版社三、用Matlab绘制根轨迹17num=[0001];%开环传递函数分子系数,降幂排列den=[1320];%开环传递函数分母系数,降幂排列r=rlocus(num,den);[例子]系统的开环传递函数为:,试利用Matlab画出系统的根轨迹。)2)(1(1)(ssssGk[解]打开Matlab,创建一个m文件,输入下列程序片段:执行之,可得到根轨迹。18例4-13.已知系统开环传递函数为(1)画出系统的根轨迹;(2)计算使系统稳定的k值范围;(3)计算系统对于斜坡输入的稳态误差。331sssksG1365.6732.1j732.1j[解]:(1)画根轨迹:19求出射角:,得。该系统有三条根轨迹,一条从原点起始,终止于开环零点-1处;另两条从原点以的出射角起始,分别终止于-3和无穷零点处。121801k180,60160会合分离点:由方程得解得在根轨迹上,因此是会合点。不在根轨迹上,舍去。0sDsNsDsN0982ss65.242,1s65.61s35.12s20求与虚轴交点系统特征方程为劳斯表为当时,由辅助方程,可求出根轨迹与虚轴的交点为。(2)由劳斯表可知当时,系统稳定。(3)系统含有三个积分环节,属Ⅲ型系统,Ⅲ型系统对于斜坡输入的稳态误差为零。03423kksksskskskksks3034341012343k032kks3j43k21画根轨迹分离(会合)点分别为-2.93和-17.07,分离(会合)角为90度。根轨迹为圆,如右图所示。45例4-14.已知单位反馈系统的开环传递函数为(1)画出系统的根轨迹;(2)计算当增益k为何值时,系统的阻尼比是,并求此时系统的闭环特征根;(3)分析k对系统性能的影响,并求系统最小阻尼比所对应的闭环极点。510sssksG2122当时,阻尼角,表示角的直线为OB,其方程为,代入特征方程整理后得:令实部和虚部分别为零,有解得由图可知当时直线OB与圆相切,系统的阻尼比,特征根为。2245450521052kjkk0520105kkk5,5k5k2155j23对于分离点,由幅值条件可知对于会合点,有由根轨迹图可知,当时,闭环系统有一对不等的负实数极点,其瞬态响应呈过阻尼状态。当时,闭环系统有一对共轭复数极点,其瞬态响应呈欠阻尼状态。当时,闭环系统又有一对不等的负实数极点,瞬态响应又呈过阻尼状态。93.2858.093.21093.2593.21k07.1714.2907.17100.17507.172k858.00k14.29858.0kk14.294524由坐标原点作根轨迹圆的切线,此切线就是直线OB,直线OB与负实轴夹角的余弦就是系统的最小阻尼比,由上可知,此时系统的闭环极点为。55j25[例4-15]:设系统A和B有相同的被控对象,且有相同的根轨迹,如下图所示。已知系统A有一个闭环零点,系统B没有闭环零点。试求系统A和B的开环传递函数和它们所对应的闭环方块图。26②系统A和B的闭环传递函数分别为:211211211222sssksssksksssksDsksA211212222sssksskskssksDksB212sssksG122sksssD[解]:①由于两系统的根轨迹完全相同,因而它们对应的开环传递函数和闭环特征方程式也完全相同。由上页图可知系统A和B的开环传递函数为:特征方程为:27由此可知,系统A是一单位反馈系统,前向通路的传递函数为:。系统B的前向通路传递函数为:,反馈通路传递函数为:。由于系统A和B有相同的被控对象,因此,系统的A的前向通路传递函数可写为:,闭环方块图如下图(a)所示,系统B的闭环方块图如下图(b)所示。212sssk22ssk1s212ssks1s)2(2ssk)(sC)(sR图(a)A系统1s)2(2ssk)(sC)(sR图(b)B系统根轨迹相同的系统,开环传递函数和闭环极点都相同,但闭环零点却不一定相同。28[例4-16]:已知单位反馈系统的根轨迹如下图所示。(1)写出该系统的闭环传递函数;(2)试用适当的方法使系统在任意K值时均处于稳定的状态。29[解]:①由根轨迹图知系统的开环传递函数为:62ssksGksskkssks6622单位反馈系统的闭环传递函数为:提示:①加入比例微分控制后,系统增加了开环零点。②在系统中加入零点后,将使根轨迹左移,有利于系统的稳定性。62ssasksG90026aa60a②当在系统中加入比例微分控制时,开环传递函数增加了一个零点,此时:这时渐近线与实轴的夹角为:,只要渐近线与负实轴相交,系统的根轨迹就在左半S平面。因此有:,所以。30从下图可以看出:a越小,根轨迹越左,稳定性越好。a6时,根轨迹全部在s左半平面。a=6时,根轨迹有一部分在虚轴上。a6时,根轨迹有一部分在s右半平面。clearall;num1=[0013];den1=[1600];num2=[0015];den2=[1600];num3=[0017];den3=[1600];h1=tf(num1,den1);h2=tf(num2,den2);h3=tf(num3,den3);rlocus(h1,h2,h3)作业:4-7,4-10,4-1131小结条件稳定系统的分析—临界稳定增益的确定;瞬态性能分析和开环系统参数的确定—阻尼角和等阻尼线;—超调量、调整时间与闭环极点的关系