自动控制原理-控制系统灵敏度分析.ppt

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第五节 控制系统灵敏度分析,控制系统在参数变化时的灵敏度是一个非常重要的概念。在开环系统中，所有的变化都会导致系统的输出产生偏差，并且系统自身没有能力消除这一偏差，这是由于开环系统没有反馈的缘故。但是，闭环系统能够察觉到输出所产生的偏差，并试图修正输出，这正是闭环反馈控制系统的一个主要好处，就是具有减少系统灵敏度的能力。,对于闭环系统 的情况，如果在所关心的复数域内，都有：,(3.71),成立，则可得到：,(3.72),那么，输出仅受到,H(s),的影响，而且,H(s),有可能是一个常数。如果,H(s)=1,，得到的结果正是期望的输入值，那就是，输出等于输入。但是，在对闭环控制系统应用式,(3.72),这样一个近似之前，必须注意式,(3.71),这一前提条件，可能会导致系统的响应为剧烈振荡，甚至于不稳定。尽管如此，增加开环传递函数,G(s)H(s),的大小会导致,G(s),对输出影响减少的事实是一个极有用的概念。因此，反馈控制系统的最重要优势就是被控过程参数,G(s),变化的影响被减少了。,|,G,(,s,),H,(,s,)|1,为描述参数变化的影响，假设被控过程,G,(,s,),发生变化，新被控过程就是,G,(,s,)+,G,(,s,),。那么，在开环情况下，输出的变化为,(3.73),在闭环系统中，有,(3.74),考虑到,则输出的改变就是：,(3.75),通常情况下，有,G(s)H(s)G(s)H(s),，于是：,(3.76),C,(,s,)=,G,(,s,),R,(,s,),观察式,(3.76),可以看出，由于,1+,G,(,s,),H,(,s,),在所关心的复数域范围内常常远大于,1,，因此闭环系统输出的变化减少了。因子,1+,G,(,s,),H,(,s,),在反馈控制系统的特征中起到了非常重要的作用。,系统灵敏度定义为系统传递函数的变化率与被控过程传递函数变化率的比值。如果系统传递函数为,则，灵敏度定义为,(3.77),取微小增量的极限形式，则式,(3.77),成为,(3.78),G,B,(,s,)=C(,s,)/,R,(,s,),很明显，从式,(3.73),可以看出，开环系统的灵敏度等于,1,。闭环系统灵敏度可以从式,(3.78),容易得到。设闭环系统的系统传递函数为,因此反馈系统关于,G(s),的灵敏度为,即,(3.79),再次可以看到，在所关心的复数域范围内,GH,(,s,),增加时，闭环系统的灵敏度将会低于开环系统的灵敏度。,同样道理，可以考察闭环系统对反馈环节,H,(,s,),改变时的系统灵敏度，令,(3.80),即,(3.81),当,G(s)H(s),很大时，灵敏度约为,1,，也就是,H(s),的变化将直接影响到系统的输出。因此，使用不随环境改变或基本恒定的反馈器件是很重要的。,由此可见，控制系统引入反馈环节后能减少因参数变化而造成的影响，尤其是因被控过程参数变化所造成的影响，这是反馈控制系统的一个重要优点。,下面介绍一个利用反馈减少灵敏度的简单例子。运算放大器是一种被广泛使用在电子线路上的集成电路器件，它的基本应用电路是图,3-36,（,a,）所示的反相放大器电路。,通常，运算放大器的增益,A,远大于,10,4,。由于输入阻抗很高，所以运算放大器的输入电流可以忽略不计，因此在节点,n,，可写出电流关系式如下,(3.82),由于放大器的增益是,A,，并且是反相接法，所以,uc=-Aun,，因此,(3.83),将,(3.83),代入,(3.82),，得到,(3.84),解出输出电压,u,c,，有,(3.85),可重写式,(3.85),如下,(,a,),电路原理图,(b),结构图,图,3-36,反相放大器,当,A,1,时，可忽略,R,1/,Rf,项，则,(3.86),其中，,k=R1/Rf,。反相放大器电路结构图如图,3-36,（,b,），图中反馈环节是,H(s)=k,，前向通道的传递函数是,G(s)=-A,。进一步，当,A1,时，反相放大器电路的传递函数为,(3.87),当运算放大器处于开环状态（即无反馈电阻,Rf,）时，相对于增益,A,的开环灵敏度为,1,。在闭环时，相对于增益,A,的闭环灵敏度为,(3.88),如果,A,=,10,4,而且,k,=0.1,，有,(3.89),则灵敏度接近于,0.001,，是开环灵敏度的千分之一。,再来考虑闭环时相对于因子,k,（或者反馈电阻,Rf,）的灵敏度。处理方法同上，得,(3.90),相对于,k,的闭环灵敏度接近于,1,。,第六节应用,MATLAB,分析控制系统的性能,这一节将用两个例子描述反馈控制的优点，同时,说明如何利用,MATLAB,来分析控制系统。系统分析的,主要内容包括如何抑制干扰、如何减小稳态误差、如,何调节瞬态响应以及如何减少系统对参数变化的影响,等。,第一个例子是带有负载转矩干扰信号的电枢控制直流电动机。开环系统结构图如图,3-37(a),所示，为了改善系统性能，加入速度反馈如图,3-37(b),所示。系统的各元器件参数值在表,3.6,中给出。,从图中可以看出，系统有,Ua,(,s,),（或,Vr,(,s,),）和,ML,(,s,),两个输入。由于这是一个线性系统，按叠加定理可以分别考虑两个输入的独立作用结果。为了研究干扰对系统的作用，可令,Ua,(,s,)=0,（或,Vr,(,s,)=0,），此时只有干扰,ML,(,s,),起作用。相反地，为了研究参考输入对系统的响应，可令,ML,(,s,)=0,。如果系统具有很好的抗干扰能力，则干扰信号,ML,(,s,),对输出,w,(,s,),的影响就应该很小，下面就来验证此结论。,参数名,R,a,K,m,J,B,K,e,K,a,K,s,参数值,1,10,2,0.5,0.1,54,1,图,3-37,速度控制系统结构图,首先，考虑图,3-37(a),所示的开环系统，从,M,L,(,s,),到,w,o,(,s,),(,此处的下标“,o”,表示开环,),的传递函数为,假设干扰信号为单位阶跃信号，即,M,L,(,s,),=1/s,。利用,MATLAB,可以计算系统的单位阶跃响应如图,3-38(a),所示，而用于分析此开环控制系统的,MATLAB,程序文本,opentach,.m,示于图,3-38(b),。,在输入信号,Ua,(s)=0,的情况下，稳态误差就是干扰响应,w,o,(,t,),的终值。在图,3-38(a),的曲线中，干扰响应,w,o,(,t,),在,t=7,秒后已近似不变，所以近似稳态误差值为,w,o,(),-0.663(,弧度,/,秒,),同样，通过计算从,M,L,(,s,),到,w,c,(,s,),(,此处下标“,c”,表示闭环,),的闭环传递函数可分析图,3-37(b),所示闭环系统的抗干扰性能。对于干扰输入的闭环传递函数为,(a),开环速度系统对阶跃干扰的响应曲线,%,开环速度控制系统对干扰信号的单位阶跃应,:opentach,.m,Ra=1;Km=10;J=2;B=0.5;Ke=0.1;,num1=1;den1=J B;,num2=Km*Ke/Ra;den2=1;,num,den=feedback(num1,den1,num2,den2);,%,干扰信号为负,num=-num;,printsys,(num,den),%wo,为输出,“o”,表示开环,wo,x,t=step(num,den);plot(t,wo,),xlabel(Timesec),ylabel,(Speed),grid,%,显示稳态误差，即,wo,的最后一个值,wo,(length(t),闭环系统对单位阶跃干扰输入的响应曲线,w,(,t,),和,MATLAB,程序文本,closedtach,.m,分别示于图,3-39(a)(b),。,同前，稳态误差就是,w,(,t,),的终值，稳态误差的近似值为,在本例中，闭环系统与开环系统对单位阶跃干扰信号的输出响应的稳态值之比为,可见通过引入负反馈已明显减小了干扰对输出的影响，这说明闭环反馈系统具有抑制噪声特性。,(a),闭环系统对阶跃干扰的响应曲线,(b)MATLAB,程序文本,:closetach,.m,图,3-39,闭环速度控制系统分析,%,闭环速度控制系统对干扰信号的单位阶跃响应,:closetach,.m,Ra=1;Km=10;J=2;B=0.5;Ke=0.1;Ka=54;Ks=1,num1=1;den1=J B;num2=Ka*Ks;den2=1;,num3=Ke;den3=1;num4=Km/Ra;den4=1;numa,dena,=parallel(num2,den2,num3,den3);,numb,denb,=series(numa,dena,num4,den4);,num,den=feedback(num1,den1,numb,denb);,%,干扰信号为负,num=-num;,printsys,(num,den),%wc,为输出,“c”,表示闭环,wc,x,t=step(num,den);plot(t,wc,),xlabel(Timesec),ylabel,(Speed),grid,%,显示稳态误差，即,wc,的最后一个值,wc,(length(t),第二个例子是分析闭环控制系统的控制器增益,K,对瞬态响应的影响。图,3-40,是闭环控制系统的结构图。在参考输入,R,(,s,),和干扰输入,N,(,s,),同时作用下系统的输出为,图,3-40,反馈控制系统的结构图,如果单纯考虑增益,K,对参考输入产生的瞬态响应的影响，可以预计增加,K,将导致超调量增加、调整时间减少和响应速度提高。在增益,K=20,和,K=100,时，系统对参考输入的单位阶跃响应曲线以及相应的,MATLAB,程序文本,gain_kr,.m,示于图,3-41,。对比两条响应曲线，可以看出上述预计的正确性。,尽管在图中不能明显看出增大,K,能减少调整时间，但是这一点可以通过观察,MATLAB,程序的运行数据得以验证。这个例子说明了控制器增益,K,是如何改变系统瞬态响应的。根据以上分析，选择,K=20,可能是一个比较好的方案。尽管如此，在做出最后决定之前还应该考虑其他因素。,(a),阶跃响应曲线,(b)MATLAB,程文本,:gain_kr,.m,图,3-41,单位阶跃输入的响应分析,%K=20,和,K=100,时，参考输入的单位阶跃响应：,gain_kr,.m,numg,=1;deng=1 1 0;,K1=100;K2=20;,num1=11 K1;num2=11 K2;den=0 1;,%,简化结构图,na,da,=series(num1,den,numg,deng);,nb,db=series(num2,den,numg,deng);,numa,dena=cloop(na,da,);,numb,denb=cloop(nb,db);,%,选择时间间隔,t=0:0.01:2.0;,c1,x,t=step(numa,dena,t);,c2,x,t=step(numb,denb,t);,plot(t,c1,-,t,c2),xlabel(Timesec),ylabel,(Cr(t),grid,在对,K,做出最后选择之前，非常重要的是要研究系统对单位阶跃干扰的响应，有关结果和相应的,MATLAB,程序文本如图,3-42,所示。从中可以看到，增加,K,减少了单位干扰响应的幅值。对于,K,=20,和,K,=100,，响应的稳态值分别为,0.05,和,0.01,。对干扰输入的稳态值可按终值定理求得,如果仅从抗干扰的角度考虑，选择,K=100,更合适。,在本例中所求出的稳态误差、超调量和调整时间（,2%,误差）归纳于表,3.7,。,(a),阶跃响应曲线,(b)MATLAB,程序文本,:gain_kn.m,图,3-42,单位阶跃干扰的响应分析,%K=20,和,K=100,时，干扰输入的单位阶跃响应：,gain_kn.m,numg,=1;deng=1 1 0;,K1=100;K2=20;,num1=11 K1;num2=11 K2;den=0 1;,%,简化结构图,numa,dena,=feedback(numg,deng,num1,den);,numb,denb,=feedback(numg,deng,num2,den);,%,选择时间间隔,t=0:0.01:2.5;,c1,x,t=step(numa,dena,t);,c2,x,t=step(numb,denb,t);,plot(t,c1,-,t,c2),xlabel(Timesec),ylabel(Cn,(t),grid,表,3.7,K,=20,和,K,=100,时，控制系统的响应特性,K,值,K=20,K,=100,超调量,s,p,%,4%,22%,调整时间,t,s,1.0,秒,0.7,秒,稳态误差,e,ss,5%,1%,(a),被控对象变化时系统的灵敏度,%,系统灵敏度分析,K=20;num=1 1 0;den=1 12 K;,%,取,s,=,jw,w,的范围为,10-1,103,，共取,200,点,w=logspace(-1,3,200);s=w*j;,%S,为灵敏度，,S2,为灵敏度的近似值,n=s.2+s;d=s.2+12*s+K;S=n./d;,n2=s;d2=K;S2=n2./d2;,subplot(211),plot(real(S),imag(S),xlabel(Real(S),ylabel(imag,(S),grid,subplot(212),loglog(w,abs(S),w,abs(S2),xlabel(wrad/sec),ylabel,(abs(S),(b)MATLAB,程序文本,图,3-43,系统的灵敏度分析,最后来分析被控对象变化时系统的灵敏度。在本例中，被控对象的传递函数和闭环系统的传递函数分别为,系统的灵敏度可由式,(3.78),得出,利用上式可计算不同,s,值所对应的灵敏度,S,G,，并绘制出频率,灵敏度曲线。图,3-43(a),中给出的是,K=20,，,s,=,j,w,，,w,=10-1 10-4,时，系统的灵敏度相对于频率,w,的变化曲线，图,3-43(b),是相应的,MATLAB,程序文本。在低频段，系统的灵敏度可近似为,可见，增大,K,值，可以减少系统的灵敏度。,第七节 设计实例,(,一,),英吉利海峡海底隧道钻机,连接法国和英国的英吉利海峡海底隧道于,1987,年,12,月,开工建造，,1990,年,11,月，从两个国家分头开钻的隧道首次,对接成功。隧道长,23.5,英里，位于海平面以下,200,英尺。隧,道于,1992,年完工，共花费,14,亿美元，每天能通过,50,趟列车。,这个工程把英国同欧洲大陆连接起来，将伦敦到巴黎的火,车行车时间缩短为,3,个小时,钻机分别从海峡两端向中间推进，并在海峡的中间对接。为了使对接达到必要的精度，施工时使用了一个激光导引系统以保持钻机的直线方向。钻机的控制模型如图所示，其中,Y,（,s,）是钻机向前的实际角度，,R,（,S,）是预期的角度，负载对机器的影响用干扰,D,（,S,）表示。,3.44,钻机控制系统框图,设计的目标是选择增益,K,使得对输入角度的响应满足工程要求，并且使干扰引起的误差最小。应用,Mason,信号流图增益公式，对两个输入的输出为：,为减少干扰的影响，我们希望增益大于,10,。当设置增益,K=100,，并令,D,（,T,）,=0,时，可得到系统对单位阶跃输入,R,（,T,）的响应，如图所示。令,r(t)=0,则可以确定系统对单位阶跃干扰的响应,y(T),如图所示。由此可见，干扰的影响是很少的。如果设置增益,K=20,，可得到系统对单位阶跃,r(t),和,d(t),的响应曲线如图所示。由于此时系统响应的超调量较小（小于,4%,），且在,2s,之内即达到稳态，所以我们选择,K=20,。为便于比较，特将上述结果同时列于表中。,表,3.8,钻机系统在两种增益情况下的响应,系统对单位阶跃输入,R,（,s,）,=1/s,的稳态误差为：,当干扰为单位阶跃,D,（,s)=1/s,输入,r(t)=0,时，,y(t),的稳态值为：,于是当,K=100,和,200,时，干扰响应的稳态值分别为,0.01,和,0.05,。,最后考虑系统对,G,（,s,）的变化的灵敏度，得：,在低频频段（），灵敏度近似为：,其中,K20,。可见当增益,K,增加时，系统的灵敏度会降低，因此我们选择,K=20,是一种合理的折中。,(,二,),火星漫游车,图所示的是以太阳能作动力的“逗留者号”火星漫游车。由于地球上发出的路径控制信号,r(t),能对该装置实施遥控。本例中我们将斜坡信号,r(t)=t,t0,用作输入信号。该装置既能以开环方式工作（不带反馈），,图,3.45,火星漫游车,又能够以闭环方式工作（带反馈），分别如图所示。本例设计的目的是使漫游车受干扰（如岩石）的影响较小，而且对增益,K,变化的灵敏度也较小。开环的传递函数是,：,闭环的传递函数是,K=2,时，,这时我们可对开环系统和闭环系统的灵敏度进行比较。,开环系统的灵敏度为：,图,3.46,漫游车的控制系统,闭环系统的灵敏度为：,为了考察系统在低频时的灵敏度，令,s=j,，得：,图,3.47,火星漫游车闭环系统的灵敏度幅值,对,K=2,，在频率,0,时，由上式可得：,由单位阶跃干扰引起的稳态误差为,-1/K,，可见增大,K,可以减少由单位阶跃干扰引起的瞬态响应。因此，我们要寻找一个较大的,K,和较大的,K/K1,，以保证系统对斜坡输入信号具有较小的稳态误差。同时我们还要保持前面已经确定的,=0.6,，以减少超调量。,为了完成设计，我们需要选择,K,。系统的特征方程为（,=0.6,）：,于是，,由式的分母可知，应同时有,故,或,的比值为：,由此可见，当选取,K=100,时，得到的是一个很,好的系统。,如果选择,K=25,，将有,k1=6,k/k1=4.17.,如果选择,K=100,，则有,K1=12,k/k1=8.33.,在实际系统中，我们还必须限制,K,，以使系统工作在线性区。当取,K=100,时，我们得到的系统如图所示。该系统对单位阶跃输入和单位阶跃干扰的响应如图所示，可以看出干扰的影响并不明显。,最后，我们可以得到系统对斜坡输入的稳态误差为,：,(,四,),循序渐进设计示例：磁盘驱动读取系统,现代磁盘在每厘米宽度内有多达,5000,个磁道，每个磁道的典型宽度仅为,1um,，因此磁盘读取系统对磁头的定位精度和磁头在磁道间移动的精度有严格的要求。本章将考虑弹性支架的影响，并在此情况下，建立磁盘驱动系统的状态空间模型。,图,3.49,磁盘驱动系统模型,考虑图所示的磁头支架，为了满足快速运动的需要，磁头支撑臂和簧片都非常轻，但必须考虑簧片对系统的影响。该簧片是由弹簧钢制成的很薄支架。,读取系统的模型如图所示，本例的实际目的在于精确控制磁头的位移,y(t)(,这对应于图,1.12,给出的设计流程的第,2,步,本例将建立土,3.49(a),所示系统的模型，为此记电机质量为,M1,，磁头质量为,M2,，簧片弹性系数为,k,，直流电机对质量,M1,的驱动力为,u,（,t,）。如果簧片是绝对刚性的（毫无弹性），则可以得到图,3.49,（,b,）所示的简化模型。该,2,体系统的典型参数列于表,3.2,。,首先，我们推导图,3.36,（,b,）所示简化系统的传递函数模型（对应设计流程）第,5,步）。注意到,M=M1+M2=20.5g=0.0205kg,，可得：,于是传递函数模型为：,表,3.2,对于表,3.2,3.45,3.45,接下来，我们推导图,3.44,（,a,）所示的,2,体系统的状态空间模型。可列写微分方程为：,质量,M1:,质量,M2:,选定状态变量为 ，便可得系统的,状态空间模型，其矩阵形式为：,其中：,图,3.46,小结,本章对控制系统分析的基本内容进行了讨论，概括地讲，就是稳定性、瞬态性能和稳态性能。另外还对控制系统关于某个环节的灵敏度以及用,MATLAB,进行系统分析等内容作了初步探讨。本章的内容是以后各章的基础。,使用闭环反馈，必然带来设备及元器件的花销和随之而来的系统复杂问题，另外，对于原本稳定的开环系统，由于反馈的引入，完全可以造成闭环系统的不稳定。但是，尽管如此，反馈控制系统在各行各业得到了广泛的应用，这是由于反馈控制能够：,（,1,）减少被控过程,G(s),中参数变化时系统的灵敏度；,（,2,）有利于控制和调节系统的瞬态响应性能；,（,3,）提高系统对干扰的抑制力；,（,4,）减小或消除系统的静态误差。,