1、1第14讲程向红奈奎斯特稳定判据对数稳定判据和稳定裕度控制系统的校正25.3.5 极坐标图的一般形状ReIm0型系统0型系统1型系统2000)1()1)(1()()1()1)(1()(2121jTjTjTjjjjKjGnmmn 00型系统:极坐标图的起点0是一个位于正实轴的有限值 极坐标图曲线的终点位于坐标原点,并且这一点上的曲线与一个坐标轴相切。11型系统:90的相角是j极坐标是一条渐近于平行与虚轴的直线的线段幅值为零,且曲线收敛于原点,且曲线与一个坐标轴相切。在总的相角中项产生的032在总相角中180的相角是由2)(j项产生的2型系统:Re01mn2mn3mn图5-34b高频区域内的极坐标
2、图 如果)(jG的分母多项式阶次)(jG的轨迹将沿者顺时针方向收敛于原点时,)(jG轨迹将与实轴或虚轴相切高于分子多项式阶次,那么当45.5奈奎斯特稳定判据(Nyquist Stability Criterion)C(s)R(s)G(s)H(s)图3-35 闭环系统闭环传递函数为)()(1)()()(sGsHsGsRsC为了保证系统稳定,特征方程0)()(1sGsH的全部根,都必须位于左半s平面。)()(sGsH的极点和零点可能位于右半s平面,但如果闭环传递函数的所有极点均位于左半s平面,则系统是稳定的。虽然开环传递函数充要条件55.5.2影射定理设)(sF为两个s的多项式之比,并设P为)(s
3、F的极点数,Z为)(sF的零点数,它们位于s平面上的某一封闭曲线内,)(sF的任何极点和零点。于是,s平面上的这一封闭曲线影射到)(sF平面上,也是一条封闭曲线。当变量s顺时针通过封闭曲线时)(sF平面上,相应的轨迹顺时针包围)(sF原点的总次数R等于Z-P。且有多重极点和多重零点的情况。设上述封闭曲线不通过在6若R为正数,表示)(sF的零点数超过了极点数;)(sF的极点数超过了零点数。)()(sGsH很容易确定)()(1)(sGsHsF的P数。因此,如果,)(sF的轨迹图中确定了R,则s平面上封闭曲线内的零点数若R为负数,表示在控制系统应用中,由很容易确定。)()()()(sAsBsGsH)
4、()()()()(1)(sAsBsAsGsHsF两者的极点数相同75.5.3影射定理在闭环系统稳定性分析中的应用为了分析线性控制系统的稳定性,令s平面上的封闭曲线包围整个右半s平面。这时的封闭曲线由整个j轴(从到该封闭曲线为奈奎斯特轨迹(轨迹的方向为顺时针方向)。因为奈奎斯特轨迹包围了整个右半s平面,所以它包围了)和右半s平面上半径为无穷大的半圆轨迹构成)()(1sGsH的所有正实部的极点和零点。)()(1sGsH则不存在闭环极点,因而系统是稳定的。如果在右半s平面不存在零点,8平面sj0图5-37 s平面内的封闭曲线ReIm平面GH1)()(1jHjG10ReIm0)()(1jHjG)()(
5、jHjG1平面GH曲线对原点的包围,恰等于)()(jGjH)()(1jGjH轨迹对-1+j0点的包围9这一判据可表示为:PRZZ函数)()(1)(sGsHsF在右半s平面内的零点数R对-1+j0点顺时针包围的次数P函数)()(sGsH如果P不等于零,对于稳定的控制系统,必须0Z或PR,这意味着必须反时针方向包围-1+j0点P次。5.5.5关于奈奎斯特稳定判据的几点说明式中在右半s平面内的极点数如果函数)()(sGsH在右半s平面内无任何极点,则RZ 因此,为了保证系统稳定,)()(jHjG的轨迹必须不包围-1+j0点。105.5.6)()(sHsG含有位于j上极点和/或零点的特殊情况平面sj0
6、j0j j j1ABC平面GHReIm,FEDFEDABC00变量s沿着j轴从 j运动到0j,从0j到0j,变量s沿着半径为1)的半圆运动,再沿着正j轴从0j运动到j()1()()(TssKsHsG11对于包含因子,3,2,1s的开环传递函数)()(sGsH,当变量s沿半径为(1)的半圆运动时,)()(sGsH的图形中将有个半径为无穷大的顺时针方向的半圆环绕原点。例如,考虑开环传递函数:)1()()(2TssKsHsGjes jeseKsHsGj22)()(lim当s平面上的9090时,)()(sGsH的相角18018012平面sj0j0j j j1ABC平面GHReImFED001在右半s平
7、面内没有极点,并且对所有的正K值,轨迹包围01j点两次。所以函数)()(1sGsH在右半s平面内存在两个零点。因此,系统是不稳定的。13如果在s平面内,奈奎斯特轨迹包含)()(1sGsH和P个极点,并且当s变量顺时针沿奈奎斯特轨迹运动时,不)()(1sGsH通过的任何极点或零点,则在)()(sGsH平面上相对应的曲线将沿顺时针方向包围01j点PZR次(负R值表示反时针包围01j点)。5.6稳定性分析的Z个零点a)不包围-1+j0如果这时)()(sGsH在右半s平面内没有极点,说明系统是稳定的;否则,系统是不稳定的。点。如果反时针方向包围的次数,等于)()(sGsH在右半s平面内没有极点数,则系
8、统是稳定的;否则系统是不稳定的。点。系统是不稳定的。c)顺时针包围-1+j0b)反时针包围-1+j014例5-3 设闭环系统的开环传递函数为:)1)(1()()(21sTsTKsGsH)()(jGjH的轨迹如图5-41所示。)()(sGsH在右半s平面内没有任何极点,并且)()(jGjH的轨迹不包围01j,所以对于任何的值,该系统都是稳定的。15Nyquist DiagramReal AxisImaginary Axis-1-0.8-0.6-0.4-0.200.20.40.60.81-0.6-0.4-0.200.20.40.6图5-41 例5-3中的)()(jGjH极坐标图 16例5-4 设系
9、统具有下列开环传递函数:)1)(1()()(21sTsTsKsGsH试确定以下两种情况下,系统的稳定性:增益K较小增益K较大。平面GHReIm001000ZRP平面GHReIm001220ZRPjjjj00小K值时是稳定的 大K值时是不稳定的 17例5-5 设开环传递函数为:)1()1()()(122sTssTKsGsH该系统的闭环稳定性取决于1T和2T相对大小。试画出该系统的奈奎斯特图,并确定系统的稳定性。ReIm00121TT 平面GH平面GHReIm00121TT 点矢量穿过01)()(jjHjG 21TT)()(sGsH的轨迹不包围01j系统是稳定的21TT)()(sGsH的轨迹通过0
10、1j点,这表明闭环极点位于轴上j18平面GHReIm00121TT 21TT)()(sGsH的轨迹顺时针方向包围01j点两次,因此系统有两个闭环极点位于右半s平面,系统是不稳定的。始开20例5-6 设一个闭环系统具有下列)1()()(TssKsHsG试确定该闭环系统的稳定性。平面GHReIm100开环传递函数:)()(jGjH极坐标图 图5-44解)1(TjjK)(1(2TjK)(1()1(2TjTjK2)(1)1(TjTjKarctg90 00,00,2)(1)1(TjTjKarctg90)1(Tj)1(Tj)1(Tj21)()(sGsH在右半s平面内有一个极点Ts11P图5-44中的奈奎斯
11、特图表明,)()(sGsH轨迹顺时针方向包围-1+0点一次1R2PRZ这表明闭环系统有两个极点在右半s平面,因此系统是不稳定的。平面GHReIm100)()(jGjH极坐标图 图5-4422)1()()(TjjKjHjG)1(Tj)1(TjjK)1(Tj)(1(2TjK)(1()1(2TjTjK2)(1)1(TjTjKarctg90 00,00,2)(1)1(TjTjKarctg90)1(Tj23例5-7 设一个闭环系统具有下列开环传递函数试确定该闭环系统的稳定性。1,)1()3()()(KsssKsHsG图5-45)()(jGjH极坐标图平面GHReIm1K)1()3(4)1(332222j
12、KjjjKjKK)1()3()1(42223KKjGH)13(4)3(jKjGH4)0(jKjGH4)0(00)(jjGH00)(jjGH)1(3)()(jjjKjHjG)1)(1()1)(3(jjjjjK解0K40渐近线24-10-8-6-4-202-15-10-5051015Real AxisImag Axis125-10-8-6-4-202-80-60-40-20020406080Real AxisImag Axis渐近线26Nyquist DiagramReal AxisImaginary Axis-8-7-6-5-4-3-2-10-150-100-500501001500 dB渐近线
13、271,)1()3()()(KsssKsHsG图5-45)()(jGjH极坐标图)()(sGsH在右半s平面内有一个极点1s1P因此开环系统是不稳定的)()(sGsH轨迹逆时针方向包围-1+j0一次1R0PRZ说明)()(1sGsH没有零点位于右半s平面内,闭环系统是稳定的。这是一个开环系统不稳定,但是回路闭合后,变成稳定系统的例子。图5-45表明平面GHReIm001KK43继续例5-728例5-8 一单位反馈控制系统的开环传递函数为)1)(1()(32221sTsTsTTKsG式中312,TTTK和均为正值。为使系统稳定,开环增益K与时间常数312,TTT和之间满足什么关系?解:频率特性)
14、(23222221)(1)()1()(jeTTTTKjG322121)(arctgTTTTarctg0)0(jKjG00)(jjG)1(1)()(32221jTjTjTTKjG29-1-0.500.511.52-1.5-1-0.500.511.5Real AxisImag Axis2,3,2,1321KTTT30-1-0.500.511.52-1.5-1-0.500.511.5Real AxisImag Axis2,3,2,1321KTTT31)1(1)()(32221jTjTjTTKjG1)()()(32232213321jTTjTTTTjTTTKjTTTTTTTTK)()(12321322
15、312令虚部为零即可 0232132TTTTT32132TTTTTc与负实轴相交于312323122312)(1)(1)(TTTTTTTTKTTTKjGcc1)(1313231TTTTTTK 1)(313231KTTTTTT展开?与负实轴的交点32-2-1012345678-6-4-20246Real AxisImag Axis8,3,2,1321KTTT33-2-1012345678-6-4-20246Real AxisImag Axis8,3,2,1321KTTT345.7.1相位裕度和增益裕度ReIm01平面G大时K小时K图5-46)1()1)(1()()1()1)(1()(2121jT
16、jTjTjjjjKjGnmmn 的极坐标图对于大的K值,系统是不稳定的。当增益减小到一定值时,)(jG的轨迹通过-1+j0点。对于小的K值,系统是稳定的。)(jG的轨迹对-1+j0点点的靠近程度,可以用来度量稳定裕量(对条件稳定系统不适用)。在实际系统中常用相位裕量和增益裕量表示。)(jG5.7相对稳定性35相位裕度、相角裕度(Phase Margin)设系统的截止频率(Gain cross-over frequency)为c1)()()(cccjHjGjA定义相角裕度为)()(180ccjHjG相角裕度的含义是度,则系统将变为临界稳定。对于闭环稳定系统,如果开环相频特性再滞后当0 时,相位裕
17、量为正值;0为了使最小相位系统稳定,相位裕度必须为正。在极坐标图上的临界点为0分贝和180时,相位裕度为负值。当36LogLogLogLog90270180Positive Gain MarginPositive Phase MarginNegative Gain MarginNegative Phase MarginStable SystemUnstable System0dB902701800dBcxcx37增益裕度、幅值裕度(Gain Margin)h设系统的相位穿越频率(Phase cross-over frequency)12()()()(kjHjGxxx,1,0 k定义幅值裕度为)
18、()(1xxjHjGh幅值裕度h对于闭环稳定系统,如果系统开环幅频特性再增大h倍,则系统将变为临界稳定状态。)()(log20 xxjHjGh的含义是,若以分贝表示,则有当增益裕度以分贝表示时,如果1h0)(dBh增益裕度为正值;1h,则0)(dBh正增益裕度(以分贝表示)表示系统是稳定的;负增益裕度(以分贝表示)表示系统是不稳定的。如果增益裕度为负值。x38LogLogLogLog90270180Positive Gain MarginPositive Phase MarginNegative Gain MarginNegative Phase MarginStable SystemUnst
19、able System0dB902701800dBcxcx39ReImh1PlaneGPositive Gain MarginPositive Phase Margin-11ReImh1PlaneGNegative Gain MarginNegative Phase Margin-11Stable SystemUnstable System)(jG)(jG400)(dBh0判断系统稳定的又一方法)()(180ccjHjG)()(log20 xxjHjGh41一阶或二阶系统的增益裕度为无穷大,因为这类系统的极坐标图与负实轴不相交。因此,理论上一阶或二阶系统不可能是不稳定的。当然,一阶或二阶系统在
20、一定意义上说只能是近似的,因为在推导系统方程时,忽略了一些小的时间滞后,因此它们不是真正的一阶或二阶系统。如果计及这些小的滞后,则所谓的一阶或二阶系统可能是不稳定的。对于稳定的最小相位系统,增益裕度指出了系统在不稳定之前,增益能够增大多少。对于不稳定系统,增益裕度指出了为使系统稳定,增益应当减少多少。一阶或二阶系统的增益裕度为多少?42只用增益裕度和相位裕度,都不足以说明系统的相对稳定性。为了确定系统的相对稳定性,必须同时给出这两个量。6030 与增益裕度应当大于6分贝。5.7.2关于相位裕度和增益裕度的几点说明控制系统的相位裕度和增益裕度是系统的极坐标图对-1+j0点 靠近程度的度量。这两个
21、裕度可以作为设计准则。对于最小相位系统,只有当相位裕度和增益裕度都是正值时,系统才是稳定的。负的裕度表示系统不稳定。适当的相位裕度和增益裕度可以防止系统中元件变化造成的影响,并且指明了频率值。为了得到满意的性能,相位裕度应当在之间,43例5-9 一单位反馈系统的开环传递函数为)05.01)(2.01()(sssKsGK=1时系统的相位裕度和增益裕度。要求通过增益K的调整,使系统的增益裕度20logh=20dB,相位裕度 40解:180)()()(xxxjHjG18005.02.090)(xxxarctgarctg即9005.02.0 xxarctgarctg2121211)(tgtgtgtgt
22、gxxxx05.02.0105.02.0 005.02.01xx 10 x相位穿越频率x增益裕度44)()(log20)(xxjHjGdBh)05.01)(2.01(1log20 xxxjjj22)1005.0(1log20)102.0(1log2010log20dB281720 在x处的开环对数幅值为45根据K=1时的开环传递函数c1)()(ccjHjG)05.01)(2.01(1)(ccccjjjjG1)0025.01)(04.01(122ccc 1c10405.02.090)(cccarctgarctg 76104180)(180c相位裕度增益穿越频率截止频率 46Bode Diagra
23、mFrequency(rad/sec)Phase(deg)Magnitude(dB)-40-30-20-1001020100101-225-180-135-90)(dBhc1K5.2K2.5KcK)(dBh)(dBh47 由题意知10h1.0)(xjG1.0)0025.01)(04.01(22xxxK 5.225.0141101.0K验证是否满足相位裕度的要求。根据 40的要求,则得:1404018005.02.090)(cccarctgarctg5005.02.0ccarctgarctg2.105.02.0105.02.0cccc 4c1)0025.01)(04.01(22cccK 2.50
24、2.128.14K不难看出,5.2K就能同时满足相位裕度和增益裕度的要求。48例5-11 设一单位反馈系统对数幅频特性如图5-50所示(最小相位系统)。写出系统的开环传递函数判别系统的稳定性如果系统是稳定的,则求ttr)(时的稳态误差。解:由图得)51)(01.01()1.01()(jjjjKjG看对数幅频特性4910-310-210-1100101102-40-20020406080-20dB/dec-20dB/dec-40dB/dec-40dB/dec0.010.115rad/sdB)(L)51)(01.01()1.01()(jjjjKjG501lg20)51(1lg20)01.01(1l
25、g20)1.01(1lg20lg20222K1110010K10K)2.01)(1001()101(10)(sssssG由于是最小相位系统,因而可通过计算相位裕度是否大于零来判断系统的稳定性。由图可知1c在c处4.1065101.011.0190)(arctgarctgarctgc则得6.73)(180c单位斜坡输入时,系统的稳态误差为1.01011vssKe 0 系统稳定515.7.3 标准二阶系统中阶跃瞬态响应与频率响应之间的关系S(S+2n)n2R(s)C(s)图3-8 标准形式的二阶系统方块图_在图3-8所示的标准二阶系统中,单位阶跃响应中的最大超调量可以精确地与频率响应中的谐振峰值联
26、系在一起。因此,从本质上看,在频率响应中包含的系统动态特性信息与在瞬态响应中包含的系统的动态特性信息是相同的。)2()(2nnsssG)2()(2nnjjjG书上例5-13p203设截止频率14)(2222nccncjGc则有5222224nncc 422244ncnc24)4(44222222nnnc )214(2422nc24214(ncnarctg290)(根据相位裕度的定义 nccarctg290180)(18022149024arctg 242142 arctg上式说明相位裕度仅仅与阻尼比有关。14)(2222nccncjG 5300.20.40.60.811.21.41.61.82
27、0102030405060708090图5-51标准二阶系统的相位裕度与阻尼比之间的关系54相位裕度与阻尼比直接相关。图5-51表示了相位裕度与阻尼比的函数关系。对于标准二阶系统,当时,相位裕度与阻尼比之间的关系近似地用直线表示如下:100因此,相位裕度相当于阻尼比。对于具有一对主导极点的高阶系统,当根据频率响应估计瞬态响应中的相对稳定性(即阻尼比)时,根据经验,可以应用这个公式。550.10.20.30.40.50.60.70.800.511.522.533.5221nr21nd对于小的阻尼比,谐振频率与阻尼自然频率的值几乎是相同的。因此,对于小的阻尼比,谐振频率的值表征了系统瞬态响应的速度
28、。56的值越小rM和pM的值越大。rrM和pM与之间的函数关系如图5-52所示。可以看出,当4.0时,rM和pM之间存在相近的关系。对于很小的值rM将变得很大,而pM却不会超过1。0.10.20.30.40.50.60.70.800.511.522.533.5575.7.4截止频率与带宽(Cutoff frequency and bandwidth)dB)(L0带宽b33图5-53 截止频率与系统带宽参看图5-53,当闭环频率响应的幅值下降到零频率值以下3分贝时,对应的频率称为截止频率。dBjRjCjRjC3)0()0(lg20)()(lg20b对于的dBjRjC0)0()0(lg20系统dB
29、jRjC3)()(lg20b58闭环系统滤掉频率大于截止频率的信号分量,但是可以使频率低于截止频率的信号分量通过。闭环系统的幅值不低于-3分贝时,对应的频率范围称为系统的带宽。带宽表示了这样一个频率,从此频率开始,增益将从其低频时的幅值开始下降。因此,带宽表示了系统跟踪正弦输入信号的能力。对于给定的n,上升时间随着阻尼比的增加而增大。另一方面,带宽随着的增加而减小。因此,上升时间与带宽之间成反比关系。59带宽指标取决于下列因素:1、对输入信号的再现能力。大的带宽相应于小的上升时间,即相应于快速特性。粗略地说,带宽与响应速度成反比。2、对高频噪声必要的滤波特性。为了使系统能够精确地跟踪任意输入信
30、号,系统必须具有大的带宽。但是,从噪声的观点来看,带宽不应当太大。因此,对带宽的要求是矛盾的,好的设计通常需要折衷考虑。具有大带宽的系统需要高性能的元件,因此,元件的成本通常随着带宽的增加而增大。60一阶系统的带宽为其时间常数的倒数。二阶系统,闭环传递函数为2222)()()(nnnssssRsC2222)2()1(1)(nnj因为1)0(j,由带宽的定义得2)2()1(2222nbnb于是 1)21(21222nb61第六章 控制系统的校正在前面几章中。讨论了控制系统几种基本方法。掌握了这些基本方法,就可以对控制系统进行定性分析和定量计算。本章讨论另一命题,即如何根据系统预先给定的性能指标,去设计一个能满足性能要求的控制系统。基于一个控制系统可视为由控制器和被控对象两大部分组成,当被控对象确定后,对系统的设计实际上归结为对控制器的设计,这项工作称为对控制系统的校正。在实际过程中,既要理论指导,也要重视实践经验,往往还要配合许多局部和整体的试验。所谓校正,就是在系统中加入一些其参数可以根据需要而改变的机构或装置,使系统整个特性发生变化,从而满足给定的各项性能指标。工程实践中常用的校正方法,串联校正、反馈校正和复合校正。62谢谢!结束