第五章 线性系统的频域分析法

第五章频率响应法5.1频率特性5.2典型环节与开环系统频率特性5.3频域稳定性判据5.4频域稳定裕度5.5闭环系统的频域性能指标5.6线性系统频域分析的MATLAB方法5.7频率特性的试验确定方法 小结 学校三防设施建设情况幼儿园教研工作小结高血压知识讲座小结防范电信网络诈骗宣传幼儿园师德小结 5.1频率特性5.1.1频率特性的定义给－个线性系统输入－个正弦信号，它的稳态输出是－个同频率的正弦信号。输入正弦信号的频率提高，则输出正弦的幅值大幅降低，相位滞后加大。系统r(t)c(t)w=0.1;t=0:pi/1000:9*pi/w;u=sin(w*t);lsim(tf(1,[111]),u,t)输入不相频率的正弦信号，测量稳态输出正弦的幅值和相位差。可绘出输入、输出幅值比和相位差与频率的关系曲线。输入、输出正弦幅值比与频率的关系定义为幅频特性A(ω)CM(ω)——输出正弦函数稳态时的幅值；RM(ω)——输入正弦函数的幅值。输入、输出正弦相位差与频率的关系，定义为系统的相频特性——输出正弦稳态时的相位角；——输入正弦的相位角。反之，已知系统的频率特性G(jω)，可求得系统的稳态输出正弦的复数形式：稳态输出正弦的复数形式与输入正弦函数的复数形式之比是－个复数，复数的幅值就是幅频特性，复数的幅角就是相频特性。因此，可定义系统的频率特性G(jω)：把输入正弦和稳态输出正弦写成复数形式则有：5.1.2频率特性的求取一、实验法实验测取频率特性的基本原理是输入正弦信号，测试输出正弦信号。只是，为保证准确性，采取了一些抗干扰、克服非线性等措施。例如，对于图示的电路,当ui(t)是正弦信号时,输出uo(t)的求取过程如下:设ui(t)=Usinωt,则其拉氏变换为二、解析法根据频率特性的定义，若已知系统的传递函数，输入正弦信号计算出稳态时的输出正弦信号，可求得系统的频率特性。电路的传递函数为对上式进行拉氏反变换,可得输出的象函数为第一项是输出的暂态分量,第二项是输出的稳态分量。当时间t→∞时,暂态分量趋于零,电路的稳态响应为输入、输出的复数形式为根据频率特性的定义，得电路的频率特性为另一方面，若令电路的传递函数中的s=jω,得到频率特性等于传递函数令s=jω。这一结论可推广到所有稳定的线性定常系统？设系统的传递函数为在零初始条件下，对应的微分方程为线性系统输入正弦函数时，稳态输出也是正弦函数，设输入输出的复数形式为，。代入上式由上式可得5.1.2频率特性的 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示方法1)直角坐标法乃奎斯特图nyquist(Gs)当频率ω从0→∞变化时，以jIm(ω)为纵坐标；以Re(ω)为横坐标。ω从-∞→0变化时的乃氏图与ω从0→∞变化时乃氏图关于实轴对称。式中Re(ω)称为系统的实频特性；Im(ω)称为系统的虚频特性。例如，－阶系统的传递函数G(jω)的轨迹上的任意一点到坐标原点的连线长度即为系统的幅频特性；连线与正实轴的夹角即为相频特性。当频率ω从0→∞变化时，可得到许多矢量，把矢量的端点连接起来，同样可得到G(jω)的轨迹，两种表示方法之间存在如下关系：2)极坐标法3)对数坐标法分别求出系统的幅频特性和相频特性，并对幅频特性取对数：对数幅频特性分贝（dB）波德图bode(Gs)以logω为横坐标，以L(ω)为纵坐标绘制对数幅频曲线；以logω为横坐标，以为纵坐标绘制对数相频曲线。请注意对数刻度和线性刻度的区别log2*L=0.3L10倍频程=一个单位长度L对数频率特性采用ω的对数分度实现了横坐标的非线性压缩,便于在较大频率范围反映频率特性的变化情况。对数幅频特性采用20lgA(ω),则将幅值的乘除运算化为加减运算。设两个环节的频率特性为：系统由这两个环节串联组成，系统频率特性为系统的对数频率特性为：即已知环节的对数幅频、相频特性曲线时，只需将其叠加，就可以得到系统的对数幅频、相频特性曲线。5.2典型环节与开环系统频率特性1.比例环节比例环节的频率特性为G(jω)=K相应的幅频特性和相频特性为开环频率特性往往是典型环节频率系统的乘积，先求出环节的频率特性，再开环频率特性就容易了。jIm(ω)0Re(ω)K乃氏图Bode图对数幅频特性和相频特性为2.积分环节积分环节的频率特性为(5.18)其幅频特性和相频特性为(5.19)由式(5.19)可见,它的幅频特性与角频率ω成反比,而相频特性恒为-90°。对数幅频特性和相频特性为(5.20)图5-8积分环节的奈氏图图5-9积分环节的伯德图3.微分环节微分环节的频率特性为(5.21)其幅频特性和相频特性为(5.22)由式(5.22)可见,微分环节的幅频特性等于角频率ω,而相频特性恒为90°。对数幅频特性和相频特性为(5.23)图5-10微分环节的奈氏图图5-11微分环节的伯德图4.惯性环节惯性环节的频率特性为(5.24)它的幅频特性和相频特性为(5.25)式(5.24)写成实部和虚部形式,即则有即所以,惯性环节的奈氏图是圆心在(0.5,0),半径为0.5的半圆(见图5-12)。对数幅频特性和相频特性为(5.26)图5-12惯性环节的奈氏图图5-13惯性环节的伯德图当ωT=1时,ω=1/T称为交接频率,或叫转折频率、转角频率。惯性环节对数幅频特性曲线的绘制方法如下:先找到ω=1/T,L(ω)=0dB的点,从该点向左作水平直线,向右作斜率为-20dB/dec的直线。在低频段和高频段,精确的对数幅频特性曲线与渐近线几乎重合。在ω=1/T附近,可以选几个点,把由式(5.26)算出的精确的L(ω)值标在图上,用曲线板光滑地连接起来,就得精确的对数幅频特性曲线。渐近线和精确曲线在交接频率附近的误差列于表5-2中。表5-2惯性环节对数幅频特性曲线渐近线和精确曲线的误差由表可知,在交接频率处误差达到最大值:一般来说,这些误差并不影响系统的分析与设计。在低频段,ω很小,Ωt<<1,φ(ω)=0°;在高频段,ω很大,ωT>>1,φ(ω)=-90°。所以,φ(ω)=0°和φ(ω)=-90°是曲线φ(ω)的两条渐近线,在交接频率处有表5-3惯性环节对数相频特性曲线角度值惯性环节对数相频特性曲线是一条中心点对称的曲线,这可以证明如下:取两个关于ω=1/T对称的频率ω1=α/T和ω2=1/(αT),则有因此有这表明φ(ω)是关于ω=1/T,φ(ω)=-45°这一点中心对称的。用MATLAB画出的惯性环节的伯德图如图5-14所示(T=1)。图5-14MATLAB绘制的惯性环节的伯德图5.一阶微分环节一阶微分环节的频率特性为(5.27)幅频特性和相频特性为(5.28)对数幅频特性和相频特性为(5.29)图5-15一阶微分环节的奈氏图图5-16一阶微分环节的伯德图6.二阶振荡环节二阶惯性环节的频率特性为(5.30)它的幅频特性和相频特性为(5.31)对数幅频特性和相频特性为(5.32)由式(5.31)得(ωT≤1)(ωT＞1)所以有(ω=0)(ω→+∞)图5-17二阶振荡环节的奈氏图画二阶振荡环节的伯德图时分析如下：在低频段,ω很小,ωT<<1,L(ω)=0dB；在高频段,ω很大,ωT>>1,L(ω)=-20lg(ωT)2=-40lg(ωT)dB。其对数幅频特性曲线可用上述低频段和高频段的两条直线组成的折线近似表示,如图5-18的渐近线所示。这两条线相交处的交接频率ω=1/T,称为振荡环节的无阻尼自然振荡频率。在交接频率附近,对数幅频特性与渐近线存在一定的误差,其值取决于阻尼比ζ的值,阻尼比越小,则误差越大,如表5-4所示。当ζ＜0.707时,在对数幅频特性上出现峰值。根据表5-5可绘制出不同阻尼比的相频特性曲线。二阶振荡环节的伯德图如图5-18所示。表5-4二阶振荡环节对数幅频特性曲线渐近线和精确曲线的误差(dB)表5-5二阶振荡环节对数相频特性曲线角度值图5-18二阶振荡环节的伯德图7.迟后环节迟后环节的频率特性为(5.33)幅频特性和相频特性为(5.34)可见,其奈氏图是一个以坐标原点为中心,半径为1的圆。对数幅频特性和相频特性为(5.35)图5-19迟后环节的奈氏图图5-20迟后环节的伯德图3、控制系统开环频率特性曲线的绘制1）开环频率特性奈氏图的绘制以后我们将会看到,在绘制奈氏图时有时并不需要绘制得十分准确,而只需要绘出奈氏图的大致形状和几个关键点的准确位置就可以了。因此,由以上典型环节奈氏图的绘制,大致可将奈氏图的一般作图方法归纳如下:(1)写出A(ω)和φ(ω)的表达式;(2)分别求出ω=0和ω=+∞时的G(jω);(3)求奈氏图与实轴的交点,交点可利用G(jω)的虚部Im［G(jω)］=0的关系式求出,也可利用∠G(jω)=n·180°(其中n为整数)求出;(4)如果有必要,可求奈氏图与虚轴的交点,交点可利用G(jω)的实部Re［G(jω)］=0的关系式求出,也可利用∠G(jω)=n·90°(其中n为正整数)求出;(5)必要时画出奈氏图中间几点;(6)勾画出大致曲线。例5-1试绘制下列开环传递函数的奈氏图:解该环节开环频率特性为ω=0,A(ω)=10,φ(ω)=0°,即奈氏图的起点为(10,j0)；ω=+∞,A(ω)=0,φ(ω)=-180°,即奈氏图的终点为(0,j0)。显然,ω从0变化到+∞,A(ω)单调递减,而φ(ω)则从0°到-180°但不超过-180°。奈氏图与实轴的交点可由φ(ω)=0°得到,即为(10,j0);奈氏图与虚轴的交点可由φ(ω)=270°(即-90°)得到,即得1-0.1ω2=0,ω2=10,则故奈氏图与虚轴的交点为(0,-j2.87)。其奈氏图如图5-21所示。用MATLAB绘制的奈氏图如图5-22所示。注意,一般手绘的奈氏图,其频率范围是0~+∞,而MATLAB绘制奈氏图时,则是从-∞~+∞。MATLAB绘制程序如下：nyquist(［10］,conv(［11］,［0.11］))图5-21例5-1的奈氏图图5-22MATLAB绘制例5-1的奈氏图例5-2已知系统的开环传递函数为试绘制其奈氏图。解该传递函数的幅频特性和相频特性分别为因此有ω=0,A(ω)=1,φ(ω)=0°和ω=+∞,A(ω)=0,φ(ω)=-∞。即奈氏图的起点为(1,j0),终点为(0,j0),随着ω的增大,曲线距离原点越来越近,相角越来越负,奈氏图与实轴和虚轴有无穷多个交点。系统的奈氏图如图5-23所示。图5-23例5-2的奈氏图例5-3设系统的开环传递函数为试绘制其奈氏图。解该传递函数的幅频特性和相频特性分别为所以有ω=0+,A(ω)=+∞,φ(ω)=-90°-Δ为正的很小量,故起点在第Ⅲ象限;ω=+∞,A(ω)=0,φ(ω)=-270°+Δ,故在第Ⅱ象限趋向终点(0,j0)。因为相角从-90°变化到-270°,所以必有与负实轴的交点。由φ(ω)=-180°得即上式两边取正切,得2ω=1/ω,即ω=0.707,此时A(ω)=0.67。因此,奈氏图与实轴的交点为(-0.67,j0)。系统的奈氏图如图5-24所示。用MATLAB绘制(-1,j0)点附近的奈氏图如图5-25所示，其程序如下：nyquist(［1］,conv(conv(［10］,［11］),［21］))例5-3中系统型次即开环传递函数中积分环节个数ν=1,若分别取ν=2,3和4,则根据积分环节的相角,将图5-24曲线分别绕原点旋转－90°,－180°和－270°即可得相应的奈氏图,如图5-26所示。图5-24例5-3的奈氏图图5-25MATLAB绘制例5-3的奈氏图图5-26ν=1,2,3,4时的奈氏图例5-4设系统的开环传递函数为其中K=0.1,T=1,T1=0.2,T2=0.5。试绘制系统的奈氏图。解该传递函数的幅频特性和相频特性分别为根据系统的幅频特性和相频特性有:ω=0+,A(ω)=+∞，φ(ω)=-90°+Δ,故奈氏图起点在第Ⅳ象限;ω=+∞,A(ω)=0,φ(ω)=-180°+Δ,故系统奈氏图在第Ⅲ象限趋向终点(0,j0)。因为相角范围为－90°~－180°,所以必有与负虚轴的交点。由φ(ω)=－90°得-90°+arctgω-arctg0.2ω-arctg0.5ω=-90°即arctgω=arctg0.2ω+arctg0.5ω上式两边取正切,得ω2=3,即ω=1.732,此时A(ω)=0.0825。所以,奈氏图与虚轴的交点为(0,-j0.0825)。系统奈氏图如图5-27所示。图5-27例5-4的奈氏图2）开环频率特性伯德图的绘制控制系统一般总是由若干环节组成的,设其开环传递函数为G(s)=G1(s)G2(s)…Gn(s)系统的开环频率特性为或(5.36)则系统的开环对数频率特性为(5.37)其中,Li(ω)=20lgAi(ω),(i=1,2,…,n)。可见,系统开环对数幅频特性和相频特性分别由各个环节的对数幅频特性和相频特性相加得到。例5-5绘制开环传递函数为的零型系统的伯德图。解系统开环对数幅频特性和相频特性分别为图5-28例5-5的伯德图实际上,在熟悉了对数幅频特性的性质后,不必先一一画出各环节的特性,然后相加,而可以采用更简便的方法。由上例可见,零型系统开环对数幅频特性的低频段为20lgK的水平线,随着ω的增加,每遇到一个交接频率,对数幅频特性就改变一次斜率。例5-6设Ⅰ型系统的开环传递函数为试绘制系统的伯德图。解系统开环对数幅频特性和相频特性分别为不难看出,此系统对数幅频特性的低频段斜率为－20dB/dec,它(或者其延长线)在ω=1处与L1(ω)=20lgK的水平线相交。在交接频率ω=1/T处,幅频特性的斜率由－20dB/dec变为－40dB/dec,系统的伯德图如图5-29所示。图5-29例5-6的伯德图通过以上分析,可以看出系统开环对数幅频特性有如下特点:低频段的斜率为－20νdB/dec,ν为开环系统中所包含的串联积分环节的数目。低频段(若存在小于1的交接频率时则为其延长线)在ω＝1处的对数幅值为20lgK。在典型环节的交接频率处,对数幅频特性渐近线的斜率要发生变化,变化的情况取决于典型环节的类型。如遇到G(s)＝(1+Ts)±1的环节,交接频率处斜率改变±20dB/dec;如遇二阶振荡环节,在交接频率处斜率就要改变－40dB/dec,等等。综上所述,可以将绘制对数幅频特性的步骤归纳如下:(1)将开环频率特性分解,写成典型环节相乘的形式;(2)求出各典型环节的交接频率,将其从小到大排列为ω1,ω2,ω3,…并标注在ω轴上;(3)绘制低频渐近线(ω1左边的部分),这是一条斜率为-20νdB/dec的直线,它或它的延长线应通过(1,20lgK)点;(4)随着ω的增加,每遇到一个典型环节的交接频率,就按上述方法改变一次斜率;(5)必要时可利用渐近线和精确曲线的误差表,对交接频率附近的曲线进行修正,以求得更精确的曲线。对数相频特性可以由各个典型环节的相频特性相加而得,也可以利用相频特性函数φ(ω)直接计算。例5-7已知系统的开环传递函数为试绘制系统的伯德图。解将开环传递函数写成如下典型环节乘积形式:可见,此系统由一个比例环节、一个积分环节、一个惯性环节、一个一阶微分环节和一个二阶振荡环节组成,且ω1=1.414,ω2=2,ω3=3。20lgK=20lg7.5=17.5。阻尼比ζ=0.354。在确定了各个环节的交接频率和20lgK的值以后,可按下列步骤绘制系统的伯德图:(1)通过点(1,17.5)画一条斜率为－20dB/dec的直线,它就是低频段的渐近线;(2)在ω1=1.414处,将渐近线的斜率从－20dB/dec改为－60dB/dec,这是考虑振荡环节的作用;(3)由于一阶惯性环节的影响,从ω2=2起,渐近线斜率应减少20dB/dec,即从原来的－60dB/dec变为－80dB/dec;(4)在ω3=3处,渐近线的斜率改变20dB/dec,形成斜率为－60dB/dec的线段,这是由于一阶微分环节的作用;(5)根据相频特性φ(ω),求出若干点的相频特性曲线角度值,如表5-6所示,将各点光滑连接,可以绘制系统的相频特性。开环系统的伯德图如图5-30所示(虚线为渐近线)。绘制程序如下：bode(［1030］,conv(conv(［10］,［12］),［112］))表5-6例5-7系统对数相频特性曲线角度值图5-30例5-7的伯德图3）最小相位系统在以上几个例子中,系统传递函数的极点和零点都位于s平面的左半部,这种传递函数称为最小相位传递函数;否则,称为非最小相位传递函数。具有最小相位传递函数的系统,称为最小相位系统;而具有非最小相位传递函数的系统,则称为非最小相位系统。对于幅频特性相同的系统,最小相位系统的相位迟后是最小的,而非最小相位系统的相位迟后则必定大于前者。当单回路系统中只包含比例、积分、微分、惯性和振荡环节时,系统一定是最小相位系统。如果在系统中存在迟后环节或者不稳定的环节(包括不稳定的内环回路)时,系统就成为非最小相位系统。对于最小相位系统,对数幅频特性与相频特性之间存在着唯一的对应关系。根据系统的对数幅频特性,可以唯一地确定相应的相频特性和传递函数,反之亦然。但是,对于非最小相位系统,就不存在上述的这种关系。实用的大多数系统为最小相位系统,为了简化工作量,对于最小相位系统的伯德图,可以只画幅频特性。例如有一最小相位系统,其频率特性为另有一非最小相位系统,其频率特性如下:(T2＞T1＞0)从图5-31不难看出,这两个系统的对数幅频特性是完全相同的,而相频特性却根本不同。前一系统的相角φ1(ω)变化范围很小,而后一系统的相角φ2(ω)随着角频率ω的增加却从0°变到趋于-180°。绘制程序如下：bode(［11］,［1001］)holdonbode(［-11］,［1001］)图5-31最小相位系统和非最小相位系统的伯德图例5-8绘制开环传递函数为的伯德图。解系统的幅频特性和相频特性分别为可见,此系统的幅频特性与惯性环节相同,而其相频特性却比惯性环节多了一项-τω。显然,它的迟后相角增加很快。开环系统的伯德图如图5-32所示。图5-32例5-8的伯德图5.3频域稳定性判据1、映射定理设有一复变函数为s为复变量,以s复平面上的s＝σ+jω表示。F(s)为复变函数,记F(s)＝U+jV。(5.38)设对于s平面上除了有限奇点之外的任一点s,复变函数F(s)为解析函数,那么,对于s平面上的每一解析点,在F(s)平面上必定有一个对应的映射点。因此,如果在s平面画一条封闭曲线,并使其不通过F(s)的任一奇点,则在F(s)平面上必有一条对应的映射曲线,如图5-33所示。若在s平面上的封闭曲线是沿着顺时针方向运动的,则在F(s)平面上的映射曲线的运动方向可能是顺时针的,也可能是逆时针的,这取决于F(s)函数的特性。我们感兴趣的不是映射曲线的形状,而是它包围坐标原点的次数和运动方向,因为这两者与系统的稳定性密切相关。图5-33s平面与F(s)平面的映射关系根据式(5.38),复变函数F(s)的相角可表示为(5.39)假定在s平面上的封闭曲线包围了F(s)的一个零点z1,而其他零极点都位于封闭曲线之外,则当s沿着s平面上的封闭曲线顺时针方向移动一周时,向量(s－z1)的相角变化-2π弧度,而其他各相量的相角变化为零。这意味着在F(s)平面上的映射曲线沿顺时针方向围绕着原点旋转一周,也就是向量F(s)的相角变化了-2π弧度,如图5-34所示。若s平面上的封闭曲线包围着F(s)的Z个零点,则在F(s)平面上的映射曲线将按顺时针方向围绕着坐标原点旋转Z周。图5-34封闭曲线包围z1时的映射情况用类似分析方法可以推论,若s平面上的封闭曲线包围了F(s)的P个极点,则当s沿着s平面上的封闭曲线顺时针移动一周时,在F(s)平面上的映射曲线将按逆时针方向围绕着原点旋转P周。综上所述,映射定理可以归纳如下:映射定理设s平面上的封闭曲线包围了复变函数F(s)的P个极点和Z个零点,并且此曲线不经过F(s)的任一零点和极点,则当复变量s沿封闭曲线顺时针方向移动时,在F(s)平面上的映射曲线按逆时针方向包围坐标原点P-Z周。2、奈奎斯特稳定判据设系统的开环传递函数为m≤n此系统的特征方程为图5-35奈氏回线由上式可见,复变函数F(s)的零点为系统特征方程的根(闭环极点)s1、s2、…、sn,而F(s)的极点则为系统的开环极点p1、p2、…、pn。闭环系统稳定的充分和必要条件是,特征方程的根,即F(s)的零点,都位于s平面的左半部。为了判断闭环系统的稳定性,需要检验F(s)是否有位于s平面右半部的零点。为此可以选择一条包围整个s平面右半部的按顺时针方向运动的封闭曲线,通常称为奈奎斯特回线,简称奈氏回线,如图5-35所示。奈氏回线由两部分组成,一部分是沿着虚轴由下向上移动的直线段C1,在此线段上s＝jω,ω由－∞变到＋∞；另一部分是半径为无穷大的半圆C2。如此定义的封闭曲线肯定包围了F(s)的位于s平面右半部的所有零点和极点。设复变函数F(s)在s平面的右半部有Z个零点和P个极点。根据映射定理,当s沿着s平面上的奈氏回线移动一周时,在F(s)平面上的映射曲线CF=1＋G(s)H(s)将按逆时针方向围绕坐标原点旋转N=P－Z周。由于闭环系统稳定的充要条件是,F(s)在s平面右半部无零点,即Z＝0。因此可得以下的稳定判据。奈奎斯特稳定判据如果在s平面上,s沿着奈氏回线顺时针方向移动一周时,在F(s)平面上的映射曲线CF围绕坐标原点按逆时针方向旋转N=P周,则系统是稳定的。根据系统闭环特征方程有G(s)H(s)=F(s)-1(5.41)这意味着F(s)的映射曲线CF围绕原点运动的情况,相当于G(s)H(s)的封闭曲线CGH围绕着(－1,j0)点的运动情况,如图5-36所示。图5-36奈氏曲线映射在F(s)平面和G(s)H(s)平面上绘制映射曲线CGH的方法是:令s＝jω代入G(s)H(s),得到开环频率特性G(jω)H(jω),按前面介绍的方法画出奈氏图,再画出其对称于实轴的、ω从0变到－∞的那部分曲线。至于映射曲线上对应于的部分,由于在实际物理系统中m≤n,当n＞m时G(s)H(s)趋近于零,n＝m时G(s)H(s)为实常数。因此,只要绘制出ω从－∞变化到＋∞的开环频率特性,就构成了完整的映射曲线CGH。综上所述,可将奈氏判据表述如下:闭环控制系统稳定的充分和必要条件是,当ω从－∞变化到＋∞时,系统的开环频率特性G(jω)H(jω)按逆时针方向包围(－1,j0)点P周,P为位于s平面右半部的开环极点数目。显然,若开环系统稳定,即位于s平面右半部的开环极点数P＝0,则闭环系统稳定的充分和必要条件是:系统的开环频率特性G(jω)H(jω)不包围(－1,j0)点。例5-9已知开环传递函数为试绘制(1)K=5,(2)K=10时的奈氏图,并判断系统的稳定性。解(1)当K=5时,开环幅频特性和相频特性分别为从而有ω=0+时,A(ω)=5,φ(ω)=0°;ω=+∞时,A(ω)=0,φ(ω)=-270°+Δ,故奈氏图在第Ⅱ象限趋向终点(0,j0)。因为相角范围为0°~－270°,所以必有与负实轴的交点。当ω=1.8时,φ(ω)=-177°,A(ω)=0.66;当ω=1.9时,φ(ω)=－181°,A(ω)=0.59。所以,当ω=ω1,1.8＜ω1＜1.9时,φ(ω)=－180°,A(ω)=A1,0.59＜A1＜0.66,因此与实轴的交点在(－1,j0)点的右侧。奈氏图如图5-37所示。因为s平面右半部的开环极点数P＝0,且奈氏曲线不包围(－1,j0)点,即N=0,则Z＝P－N=0,所以系统稳定。(2)当K=10时,奈氏图形状与(1)相同,只是以坐标原点为中心,向外“膨胀”而已。“膨胀”的倍数为10/5=2,故与实轴的交点的横坐标在(－0.59×2,－0.66×2)之间,即交点在(－1,j0)点的左侧。因为s平面右半部的开环极点数P＝0,且奈氏曲线顺时针包围(－1,j0)点2次,即N=－2,则Z＝P－N=2,所以系统不稳定,有两个闭环极点在s平面右半部。用MATLAB绘制的奈氏图如图5-38所示，其程序如下：nyquist(［5］,conv(conv(［10.5］,［11］),［12］))图5-37例5-9的奈氏图图5-38MATLAB绘制例5-9的奈氏图3、虚轴上有开环极点时的奈氏判据虚轴上有开环极点的情况通常出现在系统中有串联积分环节的时候,即在s平面的坐标原点有开环极点。这时不能直接应用图5-35所示的奈氏回线,因为映射定理要求此回线不经过F(s)的奇点。为了在这种情况下应用奈氏判据,可以选择图5-39所示的奈氏回线,它与图5-35中奈氏回线的区别仅在于,此回线经过以坐标原点为圆心,以无穷小量ε为半径的,在s平面右半部的小半圆,绕过了开环极点所在的原点。当ε→0时,此小半圆的面积也趋近于零。因此,F(s)的位于s平面右半部的零点和极点均被此奈氏回线包围在内，而将位于坐标原点处的开环极点划到了左半部。这样处理是为了适应奈氏判据的要求,因为应用奈氏判据时必须首先明确位于s平面右半部和左半部的开环极点的数目。图5-39虚轴上有极点的奈氏回线当s沿着上述小半圆移动时,有当ω从0-沿小半圆变到0+时,s按逆时针方向旋转了180°,G(s)H(s)在其平面上的映射为ν为系统中串联的积分环节数目。由以上分析可见,当s沿着小半圆从ω=0-变化到ω=0+时,θ角从－90°经0°变化到＋90°,这时在G(s)H(s)平面上的映射曲线将沿着半径为无穷大的圆弧按顺时针方向从90ν°经过0°转到－90ν°。例5-10绘制开环传递函数为的奈氏图,并判断系统的稳定性。解开环幅频特性和相频特性分别为从而有ω=0+时,A(ω)=∞,φ(ω)=-90°-Δ,Δ为正的很小量,故起点在第Ⅲ象限;ω=+∞时,A(ω)=0,φ(ω)=-270°+Δ,故在第Ⅱ象限趋向终点(0,j0)。因为相角范围从－90°到－270°,所以必有与负实轴的交点。由φ(ω)=－180°得即上式两边取正切,得0.5ω=1/ω,即ω=1.414,此时A(ω)=1.67。因此奈氏图与实轴的交点为(－1.67,j0)。系统开环传递函数有一极点在s平面的原点处,因此奈氏回线中半径为无穷小量ε的半圆弧对应的映射曲线是一个半径为无穷大的圆弧:ω:0-→0+;θ:－90°→0°→＋90°;φ(ω):＋90°→0°→－90°奈氏图如图5-40所示。因为s平面右半部的开环极点数P＝0,且奈氏曲线顺时针包围(－1,j0)点2次,即N=－2,则Z＝P－N=2,所以系统不稳定,有两个闭环极点在s平面右半部。用MATLAB绘制(－1,j0)点附近的奈氏图如图5-41所示，其程序如下：nyquist(［10］,conv(conv(［10］,［11］),［12］))图5-40例5-10的奈氏图图5-41MATLAB绘制例5-10的奈氏图例5-11绘制开环传递函数为的奈氏图,并判断系统的稳定性。解开环幅频特性和相频特性分别为从而有ω=0+时,A(ω)=∞,φ(ω)=-180°-Δ,Δ为正的很小量,故奈氏图起点在第Ⅱ象限;ω=+∞时,A(ω)=0,φ(ω)=-360°+Δ,故在第Ⅰ象限趋向终点(0,j0)。系统开环传递函数有2个极点在s平面的原点处,因此奈氏回线中半径为无穷小量ε的半圆弧对应的映射曲线是一个半径为无穷大的圆弧:ω:0-→0+;θ:－90°→0°→＋90°;φ(ω):＋180°→0°→－180°奈氏图如图5-42所示。因为s平面右半部的开环极点数P＝0,且奈氏曲线顺时针包围(－1,j0)点2次,即N=－2,则Z＝P－N=2,所以系统不稳定,有两个闭环极点在s平面右半部。用MATLAB绘制(－1,j0)点附近的奈氏图如图5-43所示，其程序如下：nyquist(［10］,conv(conv(［100］,［11］),［12］))图5-42例5-11的奈氏图图5-43MATLAB绘制例5-11的奈氏图4、对数频率稳定判据对数频率稳定判据实际上是奈氏判据的另一种形式,即利用开环系统的伯德图来判别系统的稳定性。系统开环频率特性的奈氏图(极坐标图)和伯德图之间有如下对应关系:奈氏图上以原点为圆心的单位圆对应于伯德图对数幅频特性的0分贝线;奈氏图上的负实轴对应于伯德图上相频特性的－180°线。伯德图上,φ(ω)从－180°线以下增加到－180°线以上,称为φ(ω)对－180°线的正穿越;反之,称为负穿越。图5-44例5-12的伯德图对数频率稳定判据可表述如下:闭环系统稳定的充分必要条件是,当ω由0变到∞时,在开环对数幅频特性L(ω)≥0的频段内,相频特性φ(ω)穿越－180°线的次数(正穿越与负穿越次数之差)为P/2。P为s平面右半部开环极点数目。注意,奈氏判据中,s沿着奈氏回线顺时针方向移动一周,故ω由－∞变到∞,所以伯德图中ω由0变到∞时,穿越次数为P/2,而不是P。对于开环稳定的系统,此时，P=0，若在L(ω)≥0的频段内，相频特性φ(ω)穿越-180°线的次数(正穿越与负穿越之差)为0则闭环系统稳定;否则闭环系统不稳定。例5-12系统开环传递函数为试用对数稳定判据判断其稳定性。解伯德图如图5-44所示。此系统的开环传递函数在s平面右半部没有极点,即P=0,而在L(ω)≥0的频段内,相频特性φ(ω)不穿越－180°线,故闭环系统必然稳定。5.4稳定裕度从奈氏判据可知,若系统的开环传递函数没有右半平面的极点,且闭环系统是稳定的,那么奈氏曲线G(jω)H(jω)离(－1,j0)点越远,则闭环系统的稳定程度越高;反之,G(jω)H(jω)离(－1,j0)点越近,则闭环系统的稳定程度越低;如果G(jω)H(jω)穿过(－1,j0)点,则意味着闭环系统处于临界稳定状态。这便是通常所说的相对稳定性,它通过G(jω)H(jω)对(－1,j0)点的靠近程度来度量,其定量表示为相角裕度γ和增益裕度Kg,如图5-45所示。图5-45相角裕度和增益裕度图5-45相角裕度和增益裕度1.相角裕度γ在频率特性上对应于幅值A(ω)＝1的角频率称为剪切频率,以ωc表示,在剪切频率处,相频特性距－180°线的相位差γ叫做相角裕度。图5-45(a)表示的具有正相角裕度的系统不仅稳定,而且还有相当的稳定储备,它可以在ωc的频率下,允许相角再增加(迟后)γ度才达到临界稳定状态。因此相角裕度也叫相位稳定性储备。对于稳定的系统,φ必在伯德图－180°线以上,这时称为正相角裕度,或者有正相角裕度,如图5-45(c)所示。对于不稳定系统,φ必在－180°线以下,这时称为负相角裕度,如图5-45(d)所示。故有(5.42)2.增益裕度Kg在相频特性等于－180°的频率ωg处,开环幅频特性A(ωg)的倒数称为增益裕度,记做Kg,即(5.43)在伯德图上,增益裕度改以分贝(dB)表示,Kg=-20lgA(ωg)。此时,对于稳定的系统,L(ωg)必在伯德图0dB线以下,这时称为正增益裕度,如图5-45(c)所示。对于不稳定系统,L(ωg)必在0dB线以上,这时称为负增益裕度,如图5-45(d)所示。以上表明,在图5-45(c)中,对数幅频特性还可上移Kg,即开环系统的增益增加Kg倍,则闭环系统达到稳定的临界状态。在奈氏图中,奈氏曲线与负实轴的交点到原点的距离即为1/Kg,它代表在频率ωg处开环频率特性的模。显然,对于稳定系统,1/Kg＜1,如图5-45(a)所示;对于不稳定系统有1/Kg＞1,如图5-45(b)所示。对于一个稳定的最小相位系统,其相角裕度应为正值,增益裕度应大于1。严格地讲,应当同时给出相角裕度和增益裕度,才能确定系统的相对稳定性。但在粗略估计系统的暂态响应指标时,有时主要对相角裕度提出要求。保持适当的稳定裕度,可以预防系统中元件性能变化可能带来的不利影响。为使系统有满意的稳定储备,以及得到较满意的暂态响应,在工程实践中,一般希望γ为45°~60°，Kg≥10dB,即Kg≥3。例5-13单位反馈系统开环传递函数为分别求取K1=10及K1=100时的相角裕度和增益裕度。解相角裕度可通过对数幅频特性用图解法求出。K1=10时,ω1=1,ω2=5。20lgK=20lg2=6dB。画出对数幅频特性曲线,如图5-46所示。图5-46例5-13的伯德图(幅频特性)由图可知:所以剪切频率。相角裕度为当K1从10变到100时,幅频特性上移20lg(100/10)＝20dB,如图5-46中虚线所示。所以K1=100时对应的剪切频率为。相角裕度为欲求增益裕度,则须先求出ωg,这里给出MATLAB计算的值,如图5-47所示，其程序如下：sys=tf(［100］,conv(conv(［10］,［11］),［15］));margin(sys);figuresys=tf(［10］,conv(conv(［10］,［11］),［15］));margin(sys)图5-47MATLAB绘制的例5-13的伯德图图5-47MATLAB绘制的例5-13的伯德图5.5闭环系统的频域性能指标1、由开环频率特性估计闭环频率特性对于图5-48所示的系统,其开环频率特性为G(jω)H(jω),而闭环频率特性则为因此,已知开环频率特性,就可以求出系统的闭环频率特性,也就可以绘出闭环频率特性曲线。这里介绍的是已知开环频率特性,定性地估计闭环频率特性。图5-48闭环系统设系统为单位反馈,即H(jω)=1,则一般实际系统的开环频率特性具有低通滤波的性质。所以低频时|G(jω)|>>1,则高频时|G(jω)|<<1,则图5-49闭环幅频特性例5-13中取K1=10,则单位反馈系统开环传递函数为而闭环传递函数为用MATLAB绘制其闭环频率特性的伯德图如图5-50所示，其程序如下：g1=tf(［10］,conv(conv(［10］,［11］),［15］));g2=tf(［1］,［1］)sys=feedback(g1,g2)margin(sys)图5-50MATLAB绘制的闭环频率特性伯德图2、频域性能指标图5-51闭环系统频域性能指标截止频率(带宽频率)ωb是指对数幅频特性的幅值下降到－3dB时对应的频率。带宽BW是指幅值不低于对应的频率范围,也即0~ωb的频率范围。带宽反映了系统对噪声的滤波特性,同时也反映了系统的响应速度。带宽愈大,暂态响应速度愈快。反之,带宽愈小,只有较低频率的信号才易通过,则时域响应往往比较缓慢。谐振频率ωr是指产生谐振峰值对应的频率,它在一定程度上反映了系统暂态响应的速度。ωr愈大,则暂态响应愈快。对于弱阻尼系统,ωr与ωb的值很接近。谐振峰值Mr是指闭环幅频特性的最大值。它反映了系统的相对稳定性。一般而言,Mr值愈大,则系统阶跃响应的超调量也愈大。通常希望系统的谐振峰值在1.1~1.4之间,相当于二阶系统的ζ为0.4＜ζ＜0.7。对于二阶系统,其幅频特性为由得谐振频率ωr为(0≤ζ≤0.707)(5.44)则谐振峰值Mr为(0≤ζ≤0.707)(5.45)由得截止频率(带宽频率)ωb为(0≤ζ≤0.707)(5.46)5.6频率特性的试验确定方法要想用频率特性分析或设计系统,首先要求出系统的频率特性。频率特性可用以下方法求取:(1)如果已知系统的微分方程,可将输入变量以正弦函数代入,求系统的输出变量的稳态解,输出变量的稳态解与输入正弦函数的复数比即为系统的频率特性。(2)如果已知系统的传递函数,可将传递函数中的s代之以输入变量jω,即得到系统的频率特性。(3)可通过实验的手段来求出。图5-52频率特性的实验求取用实验法测试系统的频率特性,需要解决影响测试精度的主要因素:(1)由于被测系统具有某些非线性因素或其他漂移等影响,尽管输入是正弦波信号,但输出可能含有直流分量及高次谐波;(2)由于随机干扰,主要是噪声,使输出畸变。相关分析法能从被测系统的输出信号中检出正弦波的一次谐波,同时抑制直流分量、高次谐波和噪声。线性系统频率特性G(jω)可表示为复数形式:(5.47)实频特性X(ω)、虚频特性Y(ω)与幅频特性A(ω)及相频特性φ(ω)之间有下列关系:(5.48)图5-53相关分析法测试频率特性原理示意图在输入信号ur(t)=Usinωt的作用下,被测系统的输出信号为(5.49)式中,A0——输出信号中的直流分量;u(t)——输出信号中的噪声分量;Asin(ωt+φ)——输出信号中的基波分量;Ansin(nωt+φ)——输出信号中的高次谐波分量。若以幅值为一个单位的基准信号sinωt和cosωt分别与输出信号uc(t)相乘,然后在基波的整倍数周期内积分并求平均值,则可得到基波分量的实部和虚部,而抑制掉其它分量,此即相关滤波原理。设输入信号正弦波的频率为f,周期为T,则ω＝2πf＝2π/T,取整数倍数周期NT求相关值,则有(5.50)考虑到并且从相关理论知,一个信号与另一随机信号之间的相关值,将随所取积分时间的增加而降低,即有故当N值取得较大时,式(5.50)可以写作或者写成同理可求得因此,计算相关值后,除了可将测试时输出信号中夹杂的直流分量、高次谐波分量都滤掉,噪声的影响也因NT取得足够大而忽略不计外,还能根据式(5.51)、(5.52)很方便地求得被测系统频率特性的实部和虚部。为得到一定的测试精度,通常取N＞5。由实验测出系统的频率特性后,进而可以求出系统的传递函数，下面举例说明。例5-14图5-54实线是某系统用实验测出的频率特性伯德图,试求系统的传递函数。解由幅频特性低频段可见,该系统为0型系统,且K=1。用折线(见图中虚线)作为渐近线逼近幅频特性曲线,其高频段为－40dB/dec,两个交接频率为ω1=1(rad/s),ω2=2.4(rad/s)。由此可知,该系统为二阶系统,且对于最小相位系统,二阶系统的相频特性不会小于－180°,但该系统在高频段已小于－180°,且呈现不断下降的趋势,故可断定该系统是非最小相位系统,存在迟后环节,系统的频率特性有如下形式：图5-54由实验测出的频率特性由图可见,φ(1)=－85°,故φ(1)=-τ1·1·180°/π-arctg1-arctg0.417=-85°解得τ1=0.303。另由图可见,φ(2.4)=－155°,故φ(2.4)=-τ2·2.4·180°/π-arctg2.4-arctg1=-155°解得τ2=0.310。取τ=(τ1+τ2)/2=0.307。所以即传递函数为小结频率响应法是控制理论的重要组成部分,又是研究控制系统的一种工程方法。它是一种常用的图解分析法,其特点是可以根据系统的开环频率特性去判断闭环系统的性能,并能较方便地分析系统参量对时域响应的影响,从而指出改善系统性能的途径。学习本章应掌握以下几个方面的基本 内容 财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容 :(1)频率特性的定义及其物理意义,典型环节的频率特性奈氏图和伯德图,进而绘制复杂系统的奈氏图和伯德图。虽然用MATLAB可以方便地绘制这两种图,但如果不甚明了其原理且不善于迅速地画出图像和进行实际分析,那么这种工程方法的优点也就失去了一大半。(2)若系统传递函数的极点和零点都位于s平面的左半部,则这种系统称为最小相位系统。反之,若系统的传递函数具有位于s平面右半部的极点或零点,则这种系统称为非最小相位系统。对于最小相位系统,幅频和相频特性之间存在着唯一的对应关系,即根据对数幅频特性,可以唯一地确定相应的相频特性和传递函数。而对非最小相位系统则不然。(3)奈氏稳定判据是频率响应法的核心,可以用系统的开环频率特性去判断闭环系统的稳定性。依据开环频率特性不仅能够定性地判断闭环系统的稳定性,而且可以定量地反映系统的相对稳定性,即稳定的程度。系统的相对稳定性通常用相角裕度和增益裕度来衡量。(4)时域分析中的性能指标直观反映系统动态响应的特征,属于直接性能指标。而系统频域性能指标可以作为间接性能指标。常用的闭环系统的频域性能指标有两个,一是谐振峰值Mr,反映系统的相对稳定性;另一个是频带宽度或截止频率ωb,反映系统的快速性。(5)频率响应法突出的优点是:物理意义明确并且可以用实验的方法测定出来。在难以用解析方法确定系统频率特性的情况下,这点具有特别重要的意义。一旦测出频率特性,系统的传递函数也就确定了。