自动控制原理与系统-第5章 频域分析法

频率法分析和设计控制系统的优点（1）当控制系统的结构和参数变化时，很容易确定相应的频率特性的变化。再通过频率特性指标和时域性能指标之间的关系，就将系统的结构和参数与时域性能指标联系起来了。（2）不解系统的闭环特征方程，用系统开环频率特性曲线图就可以研究闭环系统的稳定性和相对稳定性。（3）频率特性有明确的物理意义，控制系统或元部件的频率特性都可用实验方法测定，这对于很难列写运动方程的元部件或系统很有实用意义。（4）频率法不仅适用于线性定常系统的分析研究，还可推广应用于某些非线性控制系统。（5）当系统在某些频率范围内存在严重的噪声时，应用频率法，可以设计出能够抑制这些噪声的系统。5.1频率特性的基本概念与表示方法频率特性的基本概念在一般情况下，系统的传递函数为

若输入信号为正弦函数

式中

——输入信号的幅值；——输入信号的角频率。输入信号的拉氏变换为输入信号的拉氏变换为

则输出信号的拉氏变换为（5-1-2）为求出系统对正弦输入信号的响应，将式（5-1-2）展开成部分分式的形式。若系统的全部极点都不相同，则（5-1-3）式中为待定常数，为待定共轭复数。

对式（5-1-3）取拉氏变换，则得系统的输出响应函数为

对于稳定的系统，必须具有负实部，因此，当t趋于无穷大时，各项都趋于零。所以在稳态时，系统输出响应仅由式（5-1-4）前两项决定，即

如果Y(s)中含有h重极点，则y(t)中包含有一些项。对于稳定系统，由于有负实部，所以各项随趋于无穷大时都趋于零仍得式（5-1-5）结果。式（5-1-5）中系数可按下式计算：同理是一个复数，则可以用幅值和相角表示，即

（5-1-6）因为（5-1-7）式中

是的幅值;是的相角（或幅角）。同理

将式（5-1-6）、（5-1-7）代入式（5-1-5）得（5-1-8）式中——输出信号稳态分量的振幅。由式（5-1-8）可以看出系统稳态输出与输出y(t)信号频率相同。其幅值和相角通常不等于输入信号的幅值和相角，是输入信号角频率的函数。通常，把系统对正弦输入信号的稳态响应称为频率响应，当输入信号的角频率在某一范围内改变时所得到的一系列频率的响应称为这个系统的频率特性。定义为系统的幅频特性，它描述系统对不同频率输入信号的稳态响应幅值衰减（或放大）的特性。定义为系统的相频特性，它描述系统对不同频率输入信号的稳态响应，相位迟后或超前的特性。幅频特性和相频特性可由一个表达式表示，即称为系统的频率特性从上述定义可以看出，系统的频率特性可由该系统的传递函数，以代替s求得。例5-1-1求图5-1-1，电路的频率响应。

解：电路的传递函数为

式中——时间常数。图5-1-1正弦输入信号为

电路频率特性以代替可得

幅频特性为

相频特性为

系统频率特性为

电路参数给定后，

随频率变化规律就完全确定。所以频率特性反映了电路本身性质，与外界因素无关。例5-1-2设系统的传递函数，求输入信号频率为振幅为时，系统的稳态输出。解：[解题步骤]（1）输出与输入频率相同，所以

（2）求输出与输入相位差惯性环节相位落后为（3）求输出幅值

（4）稳态输出

5.1.2频率特性的表示方法系统或环节的频率特性的表示方法很多，其本质都是一样的，只是表示的形式不同而已。最常用的有；频率特性幅相频率特性对数频率特性1.幅相频率特性幅相频率特性可以表示成代数或指数形式。（1）代数表示形式设系统或环节的传递函数为，以代替s可得系统或环节的频率特性为（5-1-9）式中是频率特性的实部，称为实频特性；是频率特新的虚部，称为虚频特性。这就是频率特性的代数表示形式。（2）指数表示形式（5-1-10）式中是复数频率特性的模，称幅频特性；是复数频率特性的相位移，称相频特性。两种表示方法的关系为在工程分析和设计中，通常把频率特性画成一些曲线。幅相频率特性曲线，是以为参变量把幅频特性和相频特性同时表示在复平面上。例如，指数形式表示的频率特性，可以在极坐标中以一个矢量来表示，如图5-1-2（a）。矢量的长度等于模，相位移等于相对于极轴坐标的转角。若将极坐标与直角坐标重合，如图5-1-2（b）所示，极坐标原点与直角坐标原点重合，取极坐标轴为直角坐标轴的实轴，当由0变到时，的矢量终端将在极坐标系，或直角坐标系中描绘出一条曲线如图5-1-2（c），这条曲线称为系统或环节的幅相频率特性曲线或奈氏图。(a)(b)(c)图5-1-2幅相频率特性表示法2．对数频率特性图5-1-3对数坐标对数频率特性是将频率特性表示在半对数坐标中。对式（5-1-10）两边取以10为底的对数，得习惯上，一般不考虑0.434这个系数，只用相位移本身。图5-1-3对数坐标logω增加一个十倍频程对数频率特性曲线用两条曲线表示，即对数幅频频率特性曲线和相频频率特性曲线。通常把幅频和相频特性组成的对数频率特性曲线称为波特图。对数幅频特性曲线绘制在半对数坐标中，半对数坐标如图5-1-3横轴为角频率，采用对数比例尺（或称对数标度），每变化10倍，横坐标就增加一个单位长度，这个长度代表10倍频的距离，称之为“十倍频”或“十倍频程”（dec）。横坐标一个单位长度的值见表5-1-1所示。纵坐标以分贝（dB）为单位，等分分度。分贝值与幅频值的关系为例如则相频特性曲线。横坐标与相频特性的横坐标相同，纵坐标表示相角位移，单位为“度”或“弧度”，采用等分分度。

ω12345678910logω00.3010.4770.6020.6990.7880.8450.9030.9541

表5-1-1横坐标单位长度的刻度值对数频率特性表示法的优点：(1)在研究频率范围很宽的频率特性时，缩小了比例尺，在一张图上能很清楚地画出低频段，中频段和高频段。(2)将乘除运算变成了加减运算，使绘制频率特性曲线大大简化、方便。因为控制系统一般都是由许多环节串联构成的。串联后开环频率特性为写成指数形式：所以在极坐标中绘制幅相频率特性时，要进行乘法运算。作乘法运算，不如作加法运算方便。在半对数坐标中绘制幅频特性曲线时，由于将乘法运算变成加法运算，给绘图带来很大方便。（3）可以采用简便方法绘制近似的对数幅频曲线。5.2典型环节的频率特性1.比例环节的频率特性比例环节的传递函数为G(s)=K(1)幅相频率特性G(jω)=K=P(ω)+jQ(ω)=A(ω)(5-2-1)式中，实频特性能P(ω)=K；虚频特性Q(ω)=0.

幅频特性A(ω)=K;相频特性φ(ω)=0。由式（5-2-1）可以看出，比例环节的幅频特性、相频特性均与频率无关，所以ω由0时，G(jω)在图中为实轴上一矢量。φ(ω)=0表示输出与输入同相位。如图5-2-1。(2）对数频率特性对数幅频特性为L(ω)=20logA(ω)=20logKL(ω)为常数是平行于横轴的一条直线。当K>1时，20logK>0;K=1时20logK=0;K<1时,20logK<0。相频特性为,与横轴重合。比例环节的对数频率特性如图5-2-2所示。图5-2-1比例环节的幅相频率特性图5-2-2比例环节的对数频率特性2.惯性环节的频率特性

惯性环节的传递函数为

式中T——环节的时间常数(1)幅相频率特性

(5-2-2)式中实频特性=

虚频特性

=

幅频特性=

相频特性

给出一个频率可以算出相应的P（）和Q（）或者A（）和，分别得出在直角坐标中的一个点或者是极坐标中的一个矢量。当由0变为时，可以算出一组P（)和Q（）或者A（）和值。例如几个特征点的值列于表5-2-1中。根据这些数值，可以绘出幅相频率特性曲线图如图5-2-3所示。很容易证明，惯性环节幅相曲线是个半圆，圆心为（0.5,j0）,半径为0.5。(2)对数频率特性惯性环节对数幅频特性为从上式可以看到：在低频段Tω<<1(或ω<<1/T)时，可以近似地认为Tω=0，则L(ω)≈-20log1=0dB故在低频时，对数幅频特性可以用零分贝线近似表示，称为低频渐近线，如图5-2-4①线段。在高频段，Tω>>1(或ω>>1/T)时，忽略1，可以近似的认为L(ω)≈-20logTω这是一条斜线，其斜率为-20dB/十倍频（dB/dec）故在高频时可用斜率为-20dB/十倍频的直线近似表示，称为高频渐近线如图5-2-4②线段。上述分析表明：惯性环节对数幅频特性曲线，可由两条渐近线来表示。在0＜ω＜<1/T的范围内，是L（ω）=0的直线；在1/T＜ω＜

的范围内，是L(ω)≈-20logTω的一条斜率为-20dB/十倍频的直线。这两条直线在ω=1/T处相交，相交点的频率ω=1/T，称为转折频率或交接频率。在绘制对数幅频特性时，转折频率是一个重要参数。转折频率将近似的对数幅频特性曲线分成两段，低频段和高频段。用对数渐近幅频特性代替实际的对数幅频特性，将存在误差。最大误差发生在转折频率ω=1/T处，其误差值为L（ω=1/T）=-20log=-3dB不同频率上的误差值列于表5-2-2中。由表中可以看出，在高于或低于转折频率十倍频程处的误差值近似等于-0.04dB。惯性环节相频特性为

3.积分环节的频率特性积分环节的传递函数为幅相频率特性则。积分环节的幅相频率特性如图5-2-5所示。在0ω范围内，频率特性为负虚轴。2）对数频率特性积分环节的对数幅频特性为

由上式可知，对数幅频特性是一条斜率为-20dB/十倍频的直线，如图5-2-6所示。它在=1这一点穿过零分贝线。如果在传递函数中有N个积分环节串联，则L()=20log=-N这是一条斜率为-20dB/十倍频的斜线，且在=1处穿过横轴。积分环节的相频特性为

相位迟后与频率无关，在0ω范围内，为平行于横轴的一条直线，如图5-2-6所示。当传递函数中有N个串联积分环节时，相频特性为

图5-2-5积分环节幅相频率特性性图5-2-6积分环节对数频率特4.微分环节的频率特性理想微分环节的传递函数这G(s)=s(1)幅相频率特性G()=

（5-2-7）幅频特性为A（ω）=ω;相频特性为所以在0ω与正虚轴重合，如图5-2-7所示。范围内,幅相特性图5-2-7理想微分环节幅相频率特性图5-2-8理想微分环节对数频率特性(2)对数频率特性微分环节的对数幅频特性为

L()=20logA（）=

（5-2-8）显然，理想微分环节对数幅频特性曲线是一条斜率为20dB/十倍频，通过=1点的直线。相频特性为；是一条平行横轴的直线如图5-2-8所示。5.一阶微分环节的频率特性传递函数为G(s)=1+(1)幅相频率特性（5-2-9）则P()=1

Q()=,

,A()=

,=当由0增至大时，A()由1增至，相位由增至,是一条由实轴1点为起点，平行正虚轴的直线，如图5-2-9所示。图5-2-9一阶微分环节幅相频率特性(2)对数频率特性对数频率特性为L()=20logA（）=20log相频特性为=显然，一阶微分环节的对数幅频特性和相频特性，分别与惯性环节的相应特性相差一个符号，即对称于横轴。按照与惯性环节相似的绘图方法，绘制对数频率特性如图5-2-10所示。图5-2-10一阶微分环节对数频率特性6.振荡环节的频率特性振荡环节的传递函数为式中——时间常数；——阻尼比0<<1。

幅相频率特性

（5-2-10）由上式可求得实频特性虚频特性幅频特性相频特性振荡环节幅相频率特性的形状与阻尼比有关，以为参变量，给出一个频率就可计算出相应的和，在极坐标图上依据计算出的和就可画出一个相对应的矢量，当由0变到时，可算出一组和的值。在极坐标图上就可画出一组对应矢量。例如几个特征点的值列于表5-2-3中。根据这些数值，可以绘出由0变到时矢量端点的轨迹，该轨迹即是振荡环节的幅相频率特性曲线。如图5-2-11所示。由图可以看出，当等于无阻尼自然角频率时，即(==)时，=，=即幅值与阻尼比成反比。相位迟后为。(2)对数频率特性对数频率特性为L()=20logA（）=-20log

（5-2-11）仿照惯性环节的做法，先求对数幅频特性的渐近线，然后再修正。在低频段。Tω<<1（或ω<<横轴重合的直线，如图5-2-12①段，是低频渐近线。）范围内时，A(ω)≈1，L(ω)≈20log1=0。这是一条与在高频段。Tω>>1(或ω>>)范围内时，A(ω)≈

L(ω)≈-20log

=-40logTω

这是一条在处过横轴斜率为-40dB/dec的直线。如图5-2-12②段，是高频渐近线。低、高频渐近线的交点（在横轴的交点）ω==称为振荡环节的转折频率。这个频率是一个重要参数。渐近对数幅频特性，没有考虑阻尼比的影响，实际上精确的对数幅频特性与阻尼比有关。渐近特性与精确特性之间的误差在ω=附近的值最大。在转折频率ω=处的误差计算方法如下：用渐近特性计算中心用精确特性计算中心故误差为只有在=0.5时，二者相等。不同时精确特性如图5-2-12所示。由图可以看出=0.4～0.7，误差不大。很小时，出现一个尖峰，故在应用时要注意这个问题。相频特性为当ω=0时，=0；当ω=时，当ω=时，=。相频特性线如图5-2-12所示。由相频特性曲线上可以看出，三个特殊点的相位与阻尼比无关。整个相频特性与阻尼比有关。相频特性变化的特点是以（ω=，）点斜对称，在转折频率附近变化速度快：阻尼比越小，变化速度越快。7.二阶微分环节的频率特性二阶微分环节的传递函数为(5-2-13)式中——不具有振荡环节阻尼系数那样的物理意义。幅相频率特性幅频特性；相频特性几个特征点A(ω),的值如表5-2-4所示。二阶微分环节的幅相频率性如图5-2-13所示。(2)对数频率特性对数幅频特性为L(ω)=20logA(ω)=20log二阶微分环节的幅相频率性如图5-2-13所示。相频特性为=显然：二阶微分环节和振荡环节对数频率特性对称于横轴如图5-2-14所示。8.延迟环节的频率特性延迟环节的传递函数为幅相频率特性(5-2-14)幅频特性A(ω)=1，相频特性，故幅相频率特性是一个以原点为圆心，半径为1的圆如图5-2-15所示。对数频率特性延迟环节的对数频率特性对数幅频特性为L(ω)=20logA(ω)=0dB相频特性为

对数频率特性曲线如图5-2-16所示。5.3系统开环频率特性的绘制用实频特性P(ω)和虚频特性Q(ω)绘制，即给出不同的ω值（ω由0）计算相应P(ω)和Q(ω)的值，

在直角坐标系中得出相应的点，把这些点连接起来，就得到ω由0时的系统开环幅相频特性曲线。

也可以用幅频特性A(ω)和相频特性，给出不同的ω值，计算相应的A(ω)，的值，在极坐标系中得出相应的点，把这些点连接起来就得到ω由0时的开环幅相频率特性曲线。在对数坐标系中绘制频率特性曲线时，先绘出各环节的频率特性，然后相加，就得到对数频率特性曲线。5.3.1系统的开环幅相频率特性1.0型系统的开环幅相频率特性0型系统的开环传递函数为(n＞m)频率特性为下面分析这一类型系统幅频特性的特点。(1)在=0时，A(0)=所以幅相频特性由实轴上（K，0）点开始。(2)在=时，由于n＞m，A()=0，为坐标原点。为了确定特性线以什么角度进入坐标原点，需要求出时的相角。因为每一个惯性环节在时相位移为，每一个一阶微分环节在时相位移为，所以总相位移为：

(3)在0<的区间，频率特性形状与环节及参数有关。例5-3-1

开环系统频率特性为相位滞后随增加以一个方向连续减小，由0减到270°幅相频率特性的形状如图5-3-1所示。若频率特性为因为在分子中存在，当ω由0时，使相位移由，这时不能按一个方向连续变化，其幅相频率特性如图5-3-2所示。如果，则幅相频率特性将为图5-3-1所示形状。

图5-3-10型系统幅相频率特性图5-3-20型系统幅相频率特性2.1型系统的开环幅相频率特性1型系统的开环传递函数为n＞m开环频率特性为当时有，即幅值趋于，而相位移为。当时，其幅相频率特性如图5-3-3（a）所示。例5-3-2概略绘制开环系统幅相频率特性，其开环传递函数为令，得开环系统频率特性为由传递函数可知，系统为1型系统时时当时,幅相频率特性曲线的起点渐近线不是虚轴，而是横坐标为时的平行虚轴的直线。的求法如下：求过横轴点的坐标值，令，可得代入中得，概略幅相频率特性如图5-3-4所示图5-3-4例5-3-2的幅相频率特性

图5-3-5不同类型系统的幅相频率特性4.总结：为了概略绘制系统开环幅相频率特性，可用如下方法确定特性的几个关键部分。（1）幅相特性的低频段开环系统频率特性的一般形式为

当时，可以确定特性的低频部分。对于0型系统，特性曲线由(K,j0)点开始；对于1型系统，特性趋于一条与虚轴平行的渐近线，其横轴坐标由下式确定；（5-3-1）对于2型系统，由负实轴无穷远处开始。如图5-3-5(a)所示。（2）幅相特性的高频部分时（5-3-2）即特性总是顺时针方向终止于原点。如图5-3-5（b）所示。（3）幅相频率特性曲线与实轴和虚轴的交点与实轴交点按下述步骤求：令0,求出交点频率，代入式中，求出交点坐标。与虚轴交点，则令=0，求出交点频率，代入式中求出交点坐标。（4）如果在传递函数的分子中，无一阶微分环节，则当ω由0过程中，特性平滑的变化。如果有一阶微分环节，则视时间常数大小不同，特性可能出现凹部。n＞m5.3.2系统的开环对数频率特性

1.对数幅频特性设系统开环传递函数由n个典型环节串联组成，则令s=jω得开环频率特性为或写成：对数频率特性为

20logA(ω)=（5-3-3）（5-3-4）上式说明，如果由n个典型环节串联组成，则其对数幅频特性曲线和相频特性曲线可由典型环节对应曲线相加而得。下面举例说明。例5-3-3以开环传递函数为例介绍对数频率特性的画法。对数幅频率特性为先绘上式中各环节的对数幅频率特性，然后将各环节的对数幅频率特性的纵坐标相加，即可得到系统开环对数幅频特性曲线。如图5-3-6所示。图5-3-6开环对数频率特性第一个环节是比例环节，为平行横轴的一条直线。第二个环节是积分环节，斜率为-20dB/dec在=1处过横轴的一条直线。第三和第四个环节和均为惯性环节，转折频率分别为和。绘出各环节的、、和以后，将各环节的纵坐标相加，就得到开环对数幅频特性曲线，即图5-3-6中线。实际绘制开环对数幅频特性曲线时，不必将各环节的特性单独的绘出，再相加。而按如下步骤进行即可绘制出渐近特性曲线。(1)确定转折频率，，，，并标在横轴上。(2)在处，标出纵坐标等于20logK值的A点，其中K为开环放大系数。(3)通过A点作一条斜率为-20NdB/dec（N为积分环节数）的直线，直到第一个转折频率，若第一个转折频率的值小于1时，则该直线的延长线经过A点。(4)以后每遇到一个转折频率（含第一个转折频率），就改变一次渐近线斜率。每当遇到环节的转折频率时，渐近线的斜率增加-20dB/dec；每当遇到（1+）环节的转折频率时，斜率增加+20dB/dec；每当遇到环节的转折频率（）时，斜率增加-40dB/dec。(5)绘出的渐近线，如果需要修正，就进行修正。系统开环对数幅频特性通过横轴，即=0或A（）=1时，对应的频率称为幅值穿越频率，穿越频率是开环对数频率特性的一个重要参数。用这种方法可直接绘出如图5-3-6实线所示的幅频渐近线。2.对数相频特性已知最小相位系统开环对数频率特性曲线如图5-3-7所示。试写出开环传递函数。解：[解题分析]这是一道已知开环对数幅频特性曲线求传递函数的练习题。本题的目的是熟练掌握各典型环节的对数频率特性曲线。（1）

的低频段斜率为[-20]，故低频段为K/s。增至，斜率由[-20]转为[-40]，增加[-20]，

=所以应为惯性环节的转折频率，该环节为。

增至，斜率由[-40]转为[-20]，增加[+20]，

，所以

应为一阶微分环节的转折频率，该环节为。

增至，斜率由[-20]转为[-40]，该环节为

，斜率保持不变。故系统开环传递函数应由上述各典型环节串联组成，即（2）确定开环增益K当=时A（）=1。所以A（）=故K=3.最小相位系统在右半[s]平面上没有零点和极点的传递函数，称为最小相位传递函数，反之，则称为非最小相位传递函数。具有最小相位传递函数的系统称为最小相位系统。反之称为非最小相位系统。5.4用频率法分析控制系统的稳定性本节介绍两种频率稳定判据：用开环幅相特性曲线判断系统的稳定性（或称奈氏判据）和用开环对数频率特性曲线判断闭环系统的稳定性。下面讨论用频率法判断闭环系统的稳定性。1.一阶系统

为简单起见，我们先分析一阶系统，它的特征方程式为式中p——系统的特征根(极点)如果p为负数，则特征根位于复数

平面的左半面，系统是稳定的。如果p为正数，则特征根位于复数平面的右半平面，系统是不稳定的。假设特征根p为负实根，如图5-4-1（a）所示，令，那么可得可以把看作是一个由p点为起点到为终点的一个矢量。当沿着虚轴由0增大到时，这个矢量的矢端将沿虚轴滑动，矢量与横轴的夹角将逆时针（取正值）旋转角度。记为如果特征根是一个正实根，如图5-4-1（b）所示，则当由时，这个矢量将顺时针（取负值）旋转角度。记为图5-4-1特征向量由此得出结论：对于一阶系统，当由时矢量逆时针方向旋转角度，则系统是稳定的。否则系统不稳定。2.二阶系统它的特征方程为式中若为负值时，这一对共轭复数极点位于复数平面的左半面如图5-4-2（a）所示，若为正值时，则这一对共轭复数极点位于复数平面的右半平面。如图5-4-2（b）所示。令，则式中当由时，矢量将逆时针旋转角度。矢量将逆时针旋转角度。故矢量的相角变化量总共为图5-4-2共轭复数特征矢量记为此式就是二阶系统的稳定条件。如果共轭复数极点在平面右半平面，当由时，顺时针旋转角度记为，则系统不稳定。3.n阶系统

对于n阶系统，矢量可表示为式中为特征方程式的极点。如果n个极点都位于平面的左半平面，那么由时，将逆时针方向旋转。即如果n个极点中有一个不在左半平面。那么只有n-1个在左半平面，当由时，n-1个在左半平面的矢量所引起的相角变化为(n-1).。一个在右半平面的矢量所引起的相角变化为，所以总的相角变化量为因此，我们可以把系统的稳定条件“特征方程式的全部极点必须位于平面的左半平面”转化为下面的条件：当由时，如果矢量的相角变化为则系统稳定。否则系统不稳定。或者将稳定条件变为则系统稳定。否则系统不稳定。5.4.1用开环幅相频率特性判断闭环系统的稳定性闭环系统是否稳定，取决于闭环特征方程式的根是都在

平面的左半平面。所以用开环频率特性来研究系统的稳定性，首先就要找出开环频率特性和闭环频率特性之间的关系。闭环系统如图5-4-3所示，此系统的开环传递函数为

系统的闭环传递函数为故闭环系统的特征方程为(5-4-4)(5-4-3)由此可知闭环系统的特征方程就是开环传递函数的分母和分子之和。一个实际的系统，由于，因此与的阶次相同，均为n次。我们引入一个辅助函数为用代替s，可得辅助函数的频率特性为辅助函数也可写成如下形式：之中辅助函数的零点和极点（5-4-5）（5-4-6）由上面的分析可知，辅助函数有如下特征：1)零点为闭环传递函数的极点；

2)极点为开环传递函数的极点；

3)零点和极点的个数相等；

4)辅助函数和开换传递函数只差常数1.根据式（5-4-5）和式（5-4-6）可知，判断闭环控制系统的稳定性问题，变成辅助函数的频率特性判断闭环系统的稳定性。又因为与只差常数1，因此也就可用开环频率特性判断闭环系统的稳定性。判断闭环系统的稳定性。

下面分三种情况讨论：开环是稳定的系统开环系统是不稳定的开环有串联积分环节的系统开环是稳定的系统闭环稳定判断开环是稳定的系统如果开环系统是稳定的，那么它的n个极点应在平面左半平面，则由时，矢量的相角变化量为如果闭环系统也是稳定的，那么，它的n个极点也应在平面左半平面，则当由时，矢量的相角变化量为于是，根据式（5-4-5），当由时，矢量的相角变化量为上式表明，如果开环系统是稳定的，那么，当由时，矢量的相角变化为0，也就是说

矢端轨迹不包围原点，如图5-4-4（a）所示。

这说明闭环系统的全部极点都在平面左半平面，闭环系统是稳定的。否则，当时，矢量的相角不为0，也就是

矢端轨迹包围原点，如图5-4-4（b）所示。

这说明闭环系统有位于平面右半平面的极点，因此，闭环系统不稳定。

是开环系统的幅相频率特性。有（5-4-5）式可知，与矢量仅相差一个常数1。因此，把平面的虚轴向右移动1就变成平面的虚轴，如图5-4-5所示。这样平面中的原点，就是开环幅相频率特性的极坐标中的点，矢量矢端轨迹在平面上就是开环幅相频率特性曲线。因此，要求的轨迹不包围平面的坐标原点，就是要求开环幅频率特性曲线在平面中不包围点。由此，可以得出用开环幅相频率特性判断闭环系统稳定的判据：

图5-4-5平面与平面关系如果开环系统是稳定的，那么，当由变到时，开环频率特性曲线在平面上不包围点，则闭环系统是稳定的。否则是不稳定的。上述判据称为奈氏稳定判据。（2）开环系统是不稳定的如果开环系统是的n个极点中有p个极点在平面的右半平面，而个极点在平面的左半平面。则开环系统是不稳定的。这时，当由时，矢量的相角变化为若闭环系统的特征方程式的n个极点中，有Z个极点在平面的右半平面，个极点在平面的左半平面，则当由时，矢量的相角变化量为在这种情况下，矢量的相角变化量为式中N——表示矢量的相角变化圈数（逆时针取正）。上式表明：若开环系统有P个极点在右半平面，闭环系统有Z个极点在右半平面，则当由时，矢量的矢量轨迹在平面上逆时针围绕坐标原点转N圈。或者说，开环频率特性曲线在平面上逆时针围绕点转圈。若闭环系统是稳定的，则

由此，我们得出关于稳定判据的第二中叙述方法如下：如果开环系统不稳定，当由

时，开环频率特性曲线，逆时针方向围绕点转的圈数N等于开环系统在右半平面的极点数P，即，则闭环系统稳定，否则闭环系统不稳定。例5-4-1一个单位反馈系统，其开环传递函数为

，判断闭环系统的稳定性。

解：这是一个不稳定的惯性环节。，绘制

幅相频率特性曲线如图5-4-6所示。由图可知，由于，所以闭环系统稳定。

(3)开环有串联积分环节的系统以开环系统传递函数：为例讨论系统中有串联积分环节（即在坐标原点上有极点）时，闭环系统稳定判据的方法。频率特性为

开环频率特性在处，轨迹不连续，难以说明是否包围点。在这中情况下，可作如下处理，把频率沿虚轴变化的路线在原点处作一修改。以为圆心，为半径，在右半平面作很小的半圆。使频率的变化路线沿着这个小半圆的表达式为,将代入（5-5-11）式，则得当从变到时，角变化为。令则（5-4-8）根据（5-4-8）式就可以将在平面上的无限小半圆变换到平面上，幅相特性轨迹路线。在图5-4-7中的a点，a的幅值，相角=，变换到平面上，则，相角=。这表明，在无限小半圆上的a点，变换到平面上，在正虚轴上无穷远处的a点，如图5-4-8所示。同理，图5-4-7上的b点，c点，变换到如5-4-8上的b点，c点。这样变换后，从图5-4-8上可以看出1型系统是一个完整的开环频率特性曲线，根据这个曲线，就可以应用前面介绍的方法判断闭环系统的稳定性。如果在原点处有重根，则

式中N——重根数目。开环传递函数为2型和3型系统的频率曲线如图5-4-9所示。。

图5-4-7在原点有极点时的S轨迹图图5-4-81型系统环节幅相频率特性实际应用时，对开环有串联积分环节的系统，从

开始，以半径为无穷大，画一个逆时针转过

角的增补圆，得到一个完整的开环频率特性曲线。用奈氏判据判断闭环系统的稳定性。例5-4-2

系统开环传递函数为K>0试确定使系统稳定时，K的取值范围。

开环幅相频率特性绘于图5-4-10.当K值较大时，由图可以看出，当由时，顺时针包围点两圈，N=2，故，系统不稳定。减小K值，则时，达到临界稳定，这时

所以所以当时，不包围点，N=0，N=P系统稳定。解：由于没有极点位于右半平面，P=0。图5-4-10例5-3-2稳定判据由例5-3-1中求得5.4.2用开环对数频率特性判断闭环系统稳定性开环对数频率稳定判据和开环幅相频率特性稳定判据其本质是相同的，但是绘制开环幅相频率特性曲线比较麻烦，绘制开环对数频率特性比较简单，工作可大为简化。前面介绍的稳定判据是用由

时是否包围（-1，j0）的方法来判断，现在我们把这中判断方法转到由时

频率特性曲线在负实轴区间正，负穿越次数来判断。如图5-4-11所示。开环幅相频率特性，在负实轴区间内，由上部穿越负实轴到下部称为正穿越，正穿越次数用表示。正穿越时，的相角位移，将有正的增加。如果特性曲线逆时针包围点，则一定存在正穿越。而把在区间,由下部穿越负实轴到上部，称为负穿越，负穿越次数用表示，如图5-4-11（a）所示。负穿越时，的相角位移，将有负的增加。如果特性曲线顺时针包围点则一定存在负穿越。根据正、负穿越，可将奈氏稳定判据表述如下：当由0变到时，在复平面上，开环幅相频率特性正穿越和负穿越次数之差为，则闭环系统稳定，否则系统不稳定。即(5-4-9)

下面将这一结论用于对数频率特性：开环系统幅相频率特性与对数频率特性之间存在如下对应关系：在平面上，的单位圆，对应于对数幅频特性的0分贝线。在平面上，单位圆外的区域，对应于对数幅频特性0分贝线以上的区域。在平面上，单位圆内的区域，对应于对数幅频特性0分贝线以下的区域。在平面上的负实轴，对应于相频特性上的（）线。当幅相频率特性正穿越时产生的相位移，对应对数相频特性应由下部向上部穿越线，称为正穿越；幅相频率特性负穿越时，产生负的相位移，对应对数相频特性应由上部向下部穿越线，称为负穿越，如图5-4-11（b）所示。依据上述对应关系，对数频率稳定判据可以叙述如下：在开环对数频率特性曲线上，在大于0分贝的所有频段范围内，对数相频特性与线，正、负穿越次数之差等于，则闭环系统稳定，否则系统不稳定。当开环传递函数串联由积分环节，应从开始，补上相位移的虚线到处。应将补上的虚线看成是对数相频特性的一部分。5.4.3系统的稳定裕量在频域中，通常用相角裕量和幅值裕量这两个量来表示系统的相对稳定性。1．相角裕量对于开环传递函数P=0系统，若开环系统幅相频率特性曲线不包围点，则闭环系统稳定；若通过点，则处于临界稳定；若包围点，则系统不稳定。幅相频率特性曲线与平面上的单位圆相交处的频率称为穿越频率，如图5-4-12所示。当时，相频特性与相差的角度

称为相角裕量。相角裕量是从负实轴算起，逆时针方向规定为正，顺时针方向为负。由图5-4-12可看出，相角裕量和相角位移有如下关系：。对于最小相位系统，为了使系统稳定，相角裕量必须为正值。相角裕量的物理意义是表示使系统达到临界稳定时，还需要增加（或减少）的滞后相角量。相角裕量是设计控制系统时的一个重要指标，它和系统瞬态特性的有效阻尼比有密切的关系。通常选取2．幅值裕量相频特性时频率称为相位穿越频率；在

时，幅相频率特性的幅值的倒数称为系统的幅值裕量。即如果以分贝表示，则

（5-4-11)当时，幅值裕量的分贝数为正值，当

时，幅值裕量的分贝为负值。对于最小相位系统，幅值裕量的分贝数为正，系统稳定；幅值裕量分贝数为负，系统不稳定。

幅值裕量的物理意义表示幅值再增加（或减少）

倍，使系统达到临界稳定。通常要求 幅值裕量大6分贝。相角裕量和幅值裕量在对数标图上的表示如图5-4-13所示对最小相位系统要求，在穿越频率

处，对数幅频特性斜率通常为。例5-4-3

单位反馈系统开环传递函数为试分析不同K值时系统的稳定性确定当时系统的幅值裕量。解：本题介绍粗略地绘制系统开环幅相频率特性曲线的方法，并用幅相特性曲线（奈氏）稳定判据，判断系统放大系数的变化对闭环系统稳定性的影响。图5-4-13对数坐标图上幅值裕量和相角裕量首先将幅相频率特性用实部和虚部表示。实部：虚部：2）求曲线与负实轴的交点。与负实轴交点处求出对应此点的即，将代入实部，若闭环稳定，则所以3）粗略作的幅相曲线当当与负实轴交点当即可粗略的绘出幅相特性曲线如图5-4-14所示4）讨论不同K值时系统对稳定性的影响。当，当当

5）求幅值裕量。依据幅值裕量定义：例5-4-4

系统开环传递函数为1)试绘制系统的开环对数频率特性并计算；2)判断闭环系统的稳定性。解：1）首先将分成几个典型环节。显见该系统由放大环节，积分环节，惯性环节，一阶微分环节组成。

2)分别做各环节的对数频率特性曲线。对数幅频特性：相频特性：其对数频率特性曲线如图5-4-15所示。其对数频率特性曲线如图5-4-15所示。计算4）判断闭环系统的稳定性。由图可知部分，线无穿越，故系统闭环稳定。5.5系统开环频率特性和时域特性的关系5.5.1开环对数频率特性低频段与稳态误差1.0型系统0型系统的开环频率特性为在低频段范围内的对数幅频特性为如图5-5-1所示。0型系统对数幅频特性的低频段有如下特点：在低频段（），斜率为0dB/dec低频段的幅值为20logdB由低频段的分贝值就可确定稳态位置误差系数图5-5-10型系统对数幅频特性的低频段2.1型系统1型系统的开环频率特性为在低频段范围内的对数幅频特性为

20log当=1时，L（1）=20logK；或L(K)=0时，故K=K。由于有1个积分环节，低频段斜率为-20dB/dec,转折频率＞K时，低频段渐近线将在K=K处穿过横轴，如图5-5-2（a）所示。若＜K时，这时低频段渐近线不穿过横轴，其延长线在K=K处与横轴相交如图5-5-2（b）所示。图5-5-21型系统对数幅频特性的低频段1型系统对数幅频特性的低频段有如下特点：在低频段的渐近线的斜率为-20dB/dec。低频渐近线（或延长线）与横轴相交，交点频率为。低频渐近线（或延长线）在=1时的幅值为20logKdB。故可由低频段确定系统的速度误差系数Kv==K。K=K。

K

3．2型系统2型系统的开环频率特性为在低频段范围内的对数幅频特性为

=20log低频段斜率为-40dB/dec。在=1时，

L（1）=20logK或者L（K）=0时，

=。其对数幅频特性曲线如图5-5-3所示。

图5-5-32型系统对数幅频特性的低频段2型系统的对数幅频特性低频段有如下特点：低频渐近线斜率为-40dB/十倍频低频渐近线（或延长线）与横轴交点频率 低频渐近线（或延长线）在=1时的幅值为20logKdB可由低频渐段确定系统的加速度误差系数=

5.5.2开环对数频率特性中频段与时域瞬态特性1．开环对数频率特性中频段的基本特性对数幅频特性与横轴交点左右频段范围称为中频段。系统的瞬态特性好坏取决于中频段的特性。下面介绍常用的两种中频段特性。1．-40/-20/-40特性的传递函数为频率特性为(5-5-1)对数频率特性曲线如图5-5-4所示。在时，，即由于。近似可得相频特性为令，则相角裕量为（5-5-2）以n为参量，绘出与之间的关系如图5-5-5所示图5-5-4-2/-1/-2幅频特性图5-5-5与的关系从图中可以看出①（或K）为某一值时，相角裕量有最大值。求导令等于0得可得即选择K使时相角裕量有最大值否则相角裕量将下降。对两边取对数得(5-5-3)即在对数幅频特性中频段的中点，相角裕量最大。②最大相角裕量与中频段长度有关，中频段越长，最大相角裕量越大。2．-40/-20/-60特性（简称-2/-1/-3）-2/-1/-3特性的传递函数为频率特性为(5-5-4)

对数幅频特性曲线如图5-5-6所示。图5-5-6-2/-1/-3幅频特性当时，则相频特性为令则相角裕量为求导，令等于0则解得如果,作近似计算时：将上式代入相角裕量，可得最大相角裕量为取不同的n值时，所得绘于5-5-7中。

（5-5-5）（5-5-6）图5-5-7与n的关系

根据以上讨论可知，设计一个合理的系统，在低频段，要满足稳态精度的要求。对于中频段，要满足瞬态特性的要求来确定中频段。选择中频段要考虑以下几个方面的问题：(1)对数幅频特性曲线过横轴的斜率以-20dB/cde为宜(2)低频段和高频段可以有较大的斜率，低频段有较大的斜率，可以减少稳态误差，高频段有较大的斜率，可以更好排除高频干扰。(3)中频段的穿越频率值选择，要满足过度过程的要求。(4)中频段的长度对相角裕量有很大的影响，中频段越长，相角裕量越大。2．瞬态特性和开环频率特性的关系1）相角裕量和超调量之间的关系相角裕量为

式中二阶系统的超调量和阻尼比的关系为

（5-5-10）

（5-6-12）

（5-5-14）

上式是二阶系统与

之间的关系。将此关系绘成曲线，如图5-5-10所示。2）相角裕量和调节时间之间的关系相角裕量相同的两个系统，它们的超调量大致相同，但是它们的过渡过程时间与穿越时间成反比，穿越频率大的系统，过渡时间

就短。所以对数频率特性中的穿越频率是一个重要参数。它不仅影响系统的相角裕量，也影响系统的瞬态响应时间。5.5.3闭环系统频率特性与时域性能的关系1.闭环系统频率特性与开环系统频率特性的关系在直角坐标系中就可绘制出闭环系统频率特性曲线。如图5-5-11所示。图5-5-11闭环系统频率特性2.性能指标（1）谐振峰值——谐振峰值是闭环系统幅频特性的最大值，即谐振频率处取得的最大值。通常，越大，系统超调量也越大。系统振荡剧烈，稳定性差。（2）谐振频率——谐振频率是闭环系统幅频特性出现谐振峰值时的频率。（3）频带宽——闭环系统频率特性幅值，由其初始值时的频率或由的增益降低3dB时的频率，称为频带宽。带宽越宽，上升时间越短，但对高频抗干扰的能力越差。（4）剪切速度——剪切速度是指在高频时频率特性衰减的快慢。在高频区衰减越快，对信号和干扰两者的分辨能力越强。但剪切速度越快，谐振峰值越大。5.6用MATLAB解决频率分析问题1.Nyquist()函数命令调用格式

nyquist(sys)nyquist(sys,w)[re,im,w]=nyquist(sys)LTI对象sys可以是由函数tf()、zpk()、ss()中任何一个函数建立的开环系统模型。当函数为无等式左边输出变量格式时,函数在当前窗口中直接绘制出Nyquist曲线。nyquist(sys,w)函数用于显示绘制的系统nyquist曲线，w用来定义绘制曲线的频率范围或者是频率点。若定义频率范围,则w必须是[wmin,wmax]格式；如果是定义频率点,则必须是频率点构成的向量。[re,im,w]=nyquist(sys)可计算系统在频率w处的频率响应输出数据,而正常工作绘制曲线,其中re为频率响应的实部，im是频率响应的虚部。Re与im都是三维向量。2.Bode()函数命令调用格式bode(sys)bode(sys,w)[mag,phase,w]=bode(sys)格式与前面一样。[mag,phase,w]=bode(sys)函数中，mag是Bode图的振幅值,phase是Bode图的相位值,mag与phase是三维向量。3.控制系统中延迟环节的处理延迟环节为纯延迟环节可以用控制工具箱中函数命令pade()来近似求取。则pade()函数调用格式如下：[np,dp]=pade(tau,n)式中输入参量tau与n分别是延迟常数与pade()函数近似的阶次。延迟环节近似的传递函数可用命令tf(np,dp)来实现。4Nyquist图与Bode图的绘制图5-6-1惯性环节Nyquist图例5-6-1

试绘制惯性环节，当T=1时的Nyquist曲线与Bode图。解：在命令行输入num=[1];den=[11];s=tf(num,den);nyquist(s)bode(s)输出结果如下：例5-6-2

某具有延迟环节的非最小相位系统的开环传递函数为试绘制系统的Bode图与Nyquist曲线。解：(1)绘制系统的Bode图。nl=[1];dl=[51];G1=tf(n1,d1);tau=20;[np,dp]=pade(tau,2);Gp=tf(np,dp);G=G1*Gp;bode(G)(2)绘制系统的Nyquist曲线。nl=[1];dl=[51];Gl=tf(nl,d1);tau=20;[np,dp]=pade(tau,2);Gp=tf(np,dp);G=G1*Gp;nyquist(G)结果如下图5-6-3、图5-6-4所示。5.6.2利用频率特性计算系统的参数例5-6-3某反馈系统中，。试确定闭环系统稳定时反馈参数Kn的临界值。解：(1)求闭环系统开环传递函数对应的实频与虚频特性。symssGHUV;symsKnomegaomegacreal;s=j*omega;G=10/(s*(s-10));H=1+Kn*s;GH=G*H;U=factor(real(GH))V=factor(imag(GH))结果如下：U=-10*(1+10*Kn)/(omega-10*i)/(omega+10*i)V=-10*(-10+omega^2*Kn)/(omega-10*i)/omega/(omega+10*i)(2)当闭环系统处于临界稳定时，开环系统的频率响应G(j即Nyquist曲线将通过[G(jω)H(jω)]平面上的点(-1，j0)，此时。那么有

symsKnomegacreal;[Kn,omegac]=solve('-10*(1+10*Kn)/(omegac-10*i)/(omegac+10*i)=-l','-10*(-10+omegac^2*Kn)/(omegac-10*i)/omegac/(omegac+10*i)=0',Kn,omegac)Kn=llomegac=10^(1/2)/l^(1/2)-10^(1/2)/l^(1/2)即闭环系统稳定时反馈参数Kn=l；w=例5-6-4某控制系统的开环传递函数为，试根据该系统的频率响应，确定剪切频率ωc=5rad／s时系统的开环增益K之值解：(1)延迟环节的模(2)当=5rad／s时，对数振幅频率特性L=20log，则计算symsomegacKGH;omegac=5;GH=1/(j*omegac*(j*omegac+1)*(0.1*j*omegac+1));GH=K*abs(GH)结果如下：GH=1/325*K*130^(1/2)求开环增益K。②symsK;[K]=solve('1/325*K*130^(1/2)=1',K);K=vpa(K,3)结果如下：K=28.5例5-6-5某单位反馈系统的开环传递函数为

。试确定使系统的相角稳定裕度γ=+45时的τ之值。解：(1)确定系统的幅频特性与相频特性。symssGphi;symstauomegaomegacreal;G=(tau*s+1)/(s^2);G=subs(G,s,j*omega);Gabs=abs(G)phi=-pi+atan(tau*omega)结果如下：Gabs=(1+tau^2*omega^2)^(1/2)/omega^2phi=-pi+atan(tau*omega)确定使系统的相角稳定裕度γ=+45时的ωc。

symsomegacgam