Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
2021-12-27 16:35阅读: 1764评论: 0推荐: 0

自控理论 第6章 I 奈奎斯特判据

6.1 引入

  • 频率响应定义:正弦输入下系统的稳态响应。

6.1.1 LTI系统的频率响应

假设系统稳定(若不稳定可以串入超前补偿器使系统稳定),根据定义,系统输入为

x(t)=Xsin(ωt)˜x(s)=Xωs2+ω2

输出为

˜y(s)=G(s)˜x(s)=X(as+b)s2+ω2+{只与G(s)极点有关的部分因式}

因为G(s)稳定,所以后边那些项会随时间衰减至0,稳态响应为

Yss(s)=X(as+b)s2+ω2=X(cs+jω+¯csjω)

其中c可以由留数定理确定

c=lim

带入稳态误差化简得

\begin{aligned} Y_{ss}(s)&=\frac{jX}{2}[\frac{G(-j\omega)}{s+j\omega}-\frac{G(j\omega)}{s-j\omega}]\\ &=\frac{jX}{2}\{\frac{s[G(-j\omega)-G(j\omega)]-j\omega[G(-j\omega)+G(j\omega)]}{s^2+\omega^2}\}\\ &=X\frac{\omega\mathrm{Re}[G(j\omega)]+s\mathrm{Im}[G(j\omega)]}{s^2+\omega^2}\\ &=X|G(j\omega)|\left[\frac{\omega}{s^2+\omega^2}\cos\varphi_G(\omega)+\frac{s}{s^2+\omega^2}\sin\varphi_G(\omega)\right]\\ \end{aligned}

拉普拉斯逆变换得时域稳态响应

\begin{aligned} y_{ss}(t)&=X|G(j\omega)|\left[\sin{(\omega t)}\cos{[\varphi_G(\omega)]}+\cos{(\omega t)}\sin{\varphi_G(\omega)}\right]\\ &=X|G(j\omega)|\sin{[\omega t+\varphi_G(\omega)]} \end{aligned}

该式说明线性时不变系统对正弦输入的响应是一个频率相同、但幅值和相位变化了的正弦量。

6.1.2 频率响应的图形表示方法一览

  • Nyquist图 / Polar Plot
  • Bode图 / Logarithmic Plot
    • 幅频特性|G(j\omega)|
    • 相频特性\angle G(j\omega)
  • Nichols图(略)

6.2 奈奎斯特稳定判据

6.2.1 辐角原理

s平面(一个点对应一个s的复平面)上不通过F(s)任何零、极点的封闭曲线C_S包围s平面上F(s)Z个零点和P个极点。当s顺时针方向沿封闭曲线C_S移动一周时,在F(s)平面(一个点对应一个F(s)的复平面)上映射的封闭曲线C_F将以顺时针方向绕原点旋转N圈。NZP的关系为

N=Z-P

6.2.2 奈奎斯特稳定判据

image-20211228160537293

表示闭环传函\hat G(s)、闭环传函的分母D(s)、和开环传函G_o(s)如下

\left\{ \begin{aligned} G(s) &= \frac{N_1(s)}{D_1(s)}\\ H(s) &= \frac{N_2(s)}{D_2(s)} \end{aligned} \right.\\ \Rightarrow \left\{ \begin{aligned} \hat G(s) &=\frac{G(s)}{1+G(s)H(s)}=\frac{N_1(s)N_2(s)}{N_1(s)N_2(s)+D_1(s)D_2(s)}\\ G_o(s) &= G(s)H(s)=\frac{N_1(s)N_2(s)}{D_1(s)D_2(s)}\\ D(s)&=1+G_o(s)=\frac{N_1(s)N_2(s)+D_1(s)D_2(s)}{D_1(s)D_2(s)}\\ \end{aligned} \right.\\

  • D(s)的极点是G_o(s)的极点
  • D(s)的零点是\hat G(s)的极点

又根据辐角原理,若在s平面选取一个包围右半平面的闭合曲线,则对D(s)运用辐角原理可以得到

\begin{aligned} N&=D(s)右半零点数-D(s)右半极点数\\ &=\hat G(s)右半极点数-G_o(s)右半极点数\\ \Rightarrow \hat G(s)&右半极点数 = N+G_o(s)右半极点数 \end{aligned}

故当已知G_o(s)右半极点数和N时,便可得到\hat G(s)右半极点数,进而判断稳定性。

进一步,因为D(s)=1+G_o(s),可以很容易地使用G_o(s)代替以上表述中的D(s)来求取N,因为D(s)曲线绕原点转过的圈数等于G_o(s)曲线绕(-1,j0)点转过的圈数。最后得到:

奈奎斯特稳定判据:若系统的开环传递函数在右半平面上有P个极点,且开环频率特性曲线(s移动时对应G_o(s)在复平面上的轨迹)对(-1, j0)点包围的次数为NN>0顺时针,N<0逆时针),则闭环系统在右半平面的极点数为N + P。若N+P=0,则闭环系统稳定,否则不稳定。

6.2.3 0、1、2型系统中的奈奎斯特路径

要计算\hat G(s)右半极点数,需要首先确定s的闭合曲线,然后还要确定s的闭合曲线所映射出的G_o(s)的轨迹及其绕过(-1,j0)点的圈数N。本节要确定s的闭合曲线如何取才合适,这个闭合曲线又称作“奈奎斯特路径”。

0型系统

对应虚轴上没有零、极点,虚轴上的任意一点均可在奈奎斯特路径上。

_1_4vLeVjfTunrtXo5

大多数情况下

  • 上图中\mathrm{II}部分映射出的G_o(s)曲线会因为s\to\infty而缩成一个点,所以更多关注的是\mathrm{I}部分和第\mathrm{III}部分
  • \mathrm{I}部分和\mathrm{III}部分又共轭,所以只需求出0\to\infty\mathrm{I}部分所映射的G_o(s)曲线,就可以对称地得到\mathrm{III}部分所映射的。

因此在接下来得讨论中,更多地是在讨论G_o(j\omega)而不是G_o(s)

1型系统

此时原点处有G_o(s)的一重零点。幅角原理要求路径上不能有零极点,故选取下图所示的闭合曲线来满足原理的要求:

_2_NTGcaHrPlmVLIYC

现在讨论上图中新增的蓝色路径在G_o(s)平面上会映射出什么图形,此时有:

G_o(s)=\frac{K\prod\limits_{i=1}^m(\tau_is+1)}{s\prod\limits_{j=1}^n(T_js+1)}\\ \Rightarrow \left\{ \begin{aligned} \rho&=\lim\limits_{s\to0}|G_o(s)|=\infty\\ \theta_0&=\lim\limits_{\omega\to0^-}\angle G_o(j\omega)=-(-\frac{\pi}{2})=\frac{\pi}{2}\\ \theta_t&=\lim\limits_{\omega\to0^+}\angle G_o(j\omega)=-\frac{-\pi}{2} \end{aligned} \right.

又沿蓝色的\mathrm{IV}部分,\angle G_o(s)减小,故对应轨迹为顺时针方向旋转、半径为\infty的右半圆:

_3_aNCORjoY8vJPsiQ

右半圆的是基于开环传函写成\frac{K\prod\limits_{i=1}^m(\tau_is+1)}{s\prod\limits_{j=1}^n(T_js+1)}的形式而得到的,对于一般的系统要么化为标准形式,要么用同样的方法进行分析。下边的2型系统也是同样的道理。

2型系统

此时原点处有G_o(s)的二重零点。仍可选取1型系统的闭合曲线来满足原理的要求。此时有

G_o(s)=\frac{K\prod\limits_{i=1}^m(\tau_is+1)}{s^2\prod\limits_{j=1}^n(T_js+1)}\\ \Rightarrow \left\{ \begin{aligned} \rho&=\lim\limits_{s\to0}|G(s)|=\infty\\ \theta_0&=\lim\limits_{\omega\to0^-}\angle G_o(j\omega)=\pi\\ \theta_t&=\lim\limits_{\omega\to0^+}\angle G_o(j\omega)=-\pi \end{aligned} \right.

又沿蓝色的\mathrm{IV}部分,\angle G_o(s)减小,故对应轨迹为顺时针方向旋转、半径为\infty的整圆:

_4_ni2USRtwaf1XChd

6.2.4 绘制奈奎斯特曲线并判断稳定性

本节讨论s的闭合曲线所映射出的G_o(j\omega)的完整轨迹该如何绘制,进而求得N并判断系统的稳定性。G_o(s)的轨迹又称“开环频率特性曲线”。

  1. 确定G_o(s)右极点数,根据上一节的内容确定奈奎斯特路径

  2. G_o(j\omega)分解

    • 实部P(j\omega)
    • 虚部Q(j\omega)
    • 幅值A(\omega)
    • 相角\varphi(\omega)
  3. 求出几个特征点

    • G_o(0)或者\lim\limits_{\omega\to0}G_o(j\omega)\varphi(0)
    • \lim\limits_{\omega\to \infty}\varphi(j\omega)
    • \left.G_o(j\omega_1)\right|_{P(\omega_1)=0}
    • \left.G_o(j\omega_2)\right|_{Q(\omega_2)=0}
  4. 连接特征点,定性作出\omega0^+变化到\infty的曲线

    • 注意特征点附近的相角,虽然是定性作图但也别差得太离谱

    • 确定该段开环频率特性曲线的正穿越次数N'^+和负穿越次数N'^-

      image-20211122170330007

      \omega\to0^-时曲线刚好落在-1左侧的实轴上,这时要算0.5“次”穿越。

    • 计算N

      N=2(N'^--N'^+)

  5. 判断系统的稳定性

    • N=0,系统稳定
    • N>0,系统不稳定
    • N<0,算错了再检查检查。。。

6.3 几种特殊情况的奈奎斯特判据

不少其实是解题思路哈哈🤦‍♂️😂

6.3.1 条件稳定系统

对一些系统来说,只有一些系统参数处于一定范围内时它们才是稳定的,通过奈奎斯特判据可以求出能使系统稳定的参数区间。做法也很简单,就是先把参数作为代数画出奈奎斯特曲线,用这些代数表示曲线与实轴的交点,通过调整这些交点处于(-1,j0)点的左侧或右侧即可改变系统的稳定性。

6.3.2 多回路系统

首先应将内环系统视作一个独立的系统,用奈奎斯特判据求取其在右半平面的极点数,然后将内环闭环在右半平面的极点数与其它环节的右半极点数相加,得到原开环系统的右半极点数P,最后再对整体用奈奎斯特判据得到整体的闭环右半极点数。举个例子:

image-20211227154116930

首先把红框中的抽出来,用奈奎斯特判据得到其闭环右半极点数P_1,然后确定G_1G_3H_1、的右半极点数分别为P_1P_3P_4,则原开环系统的右半极点数为P=P_1+P_2+P_3+P_4,化简框图得到原系统的开环传递函数,使用奈奎斯特判据得到原系统的N,则原系统的右半极点数即为P+N

6.3.3 需要极点处于特定区域内

这个需要蛮有意义的,比如要求超调量小,则会对\zeta有要求,即主导极点的相角应小于某个值;再比如要求响应迅速,则会对\omega_n有要求,也即极点的实部应小于某个值。就拿这两点举例,则此时下图折线的右侧不应有极点:

image-20211227155506133

回忆奈奎斯特判据的推导过程,奈奎斯特曲线的左边界是虚轴求出的是在右半平面的极点个数,相应的如果把左边界换成上图的折线即可求出处于折线右端的极点的个数。不过要注意此时求特征点就复杂一些了。

6.3.4 最小相位系统

右半平面无零点和极点的系统称为最小相位系统,即有P=0,故要求N=0即可保证系统稳定。

  • 纯时延的最小相位系统

    纯时延的最小相位系统又多了一项e^{-T_ds},故其与负实轴的交点变为无穷多个,但因为前提是最小相位系统,保证所有交点都在(-1,j0)右侧即可保证系统稳定。

6.3.5 在虚轴除原点以外有开环零点、极点的系统

用无限小半圆路径扣掉极点,对应的开环频率特性曲线需要根据具体情况进行讨论。

6.3.6 根据开环右零点数判断开环右极点数

如果已知开环右零点数和奈奎斯特曲线,是可以直接得到开环右极点数的,举个例子来说明:

image-20211227162356622

这个图:小蓝色下半圆箭头标反了;虚线当成实线看就好。

假设A(-1,j0)点,现在已知G_o(s)的有零点数为Z_0=1。奈奎斯特曲线是G_o(s)的曲线,由辐角原理(注意不是奈奎斯特判据)知道该曲线顺时针(默认奈奎斯特轨迹是顺时针包围右半平面的)绕过原点的圈数N_0G_o(s)的右极点数P_0G_o(s)的右零点数Z满足N_0=Z_0-P_0,由上图知道N_0=-2,则P_0=Z_0-N_0=3

6.3.7 包围左半平面的奈奎斯特路径

相应的用奈奎斯特判据求出的是闭环左极点数,如果闭环左极点数等于系统阶次,那么仍然可以判断系统稳定;反之如果闭环左极点数小于系统阶次,那么说明系统在闭右半平面有极点,系统不稳定。

本文作者:Harold_Lu

本文链接:https://www.cnblogs.com/harold-lu/p/15736804.html

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   Harold_Lu  阅读(1764)  评论(0编辑  收藏  举报
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起