自动控制原理上课笔记(不定期更新)

前言

根据上课 ppt 做的笔记,不定期更新。主要涉及系统的时域和频域分析。

文章比较长可作平时查阅使用,如果需要笔记的 pdf 电子文档请在评论区留下电子邮箱。

自动控制的一般概念

自动控制指在没有人直接参与的情况下,利用外加的设备或装置(称控制装置或控制器),使机器、设备或生产过程(统称被控对象(plant))的某个工作状态或参数(即被控量)自动地按照预定的规律运行。


以反馈控制原理为基础,以传递函数为工具,主要研究单输入、单输出反馈控制系统的理论问题,又常称为反馈控制理论。

其控制又分为:

  • 开环控制
  • 反馈/闭环控制
  • 复合控制

经典控制理论和现代控制理论的区别

研究对象数学工具常用分析方法局限性
经典控制理论单输入-单输出、线性定常系统微分方程、传递函数时域分析法、频域分析法、根轨迹分析法对复杂多变量系统、时变和非线性系统无能为力
现代控制理论多输入-多输出变系数、非线性等系统线性代数、矩阵理论状态空间法比较繁琐(但由于计算机技术的迅速发展,这一局限性已克服)

自动控制的基本原理与方式

反馈:把取出输出量送回输入端,并与输入信号相比较偏差信号的过程(正反馈与负反馈)。

反馈控制(闭环控制):能够在存在扰动和不确定因素的情况下,力图减小系统的输出量与参考输入量(也称参据量)之间的偏差的控制

自动控制系统的工作原理

  1. 当被控量偏离给定值的要求时(一是给定值变化;二是给定值不变而由干扰作用引起被控量变化)。
  2. 测量元件测得被控量的大小,经变换后由比较元件将其与给定值比较得出偏差。
  3. 根据偏差的大小,经放大、调节、执行等元件后,产生控制作用;偏差大,控制作用大;反之亦然。
  4. 控制作用使被控量回复到或趋近于要求值,从而使偏差消除或减小。

典型反馈(闭环)控制系统的系统框图

在这里插入图片描述


例如:人取桌上书的过程。

控制目标(要求):人用手取到书

通过输入信号(理想书位置)和观测信号(看到的书位置)的偏差控制被控对象(手)使输出量(实际书位置)达到设定值(理想书位置)。

在这里插入图片描述


自动控制系统基本控制方式

  • 开环控制方式

    • 按给定量(输入)控制

      缺点:对可能出现的被控量偏离给定值的偏差没有任何修正能力,抗干扰能力差,控制精度不高。

      在这里插入图片描述

    • 按扰动(顺馈)控制

      定义:利用干扰信号产生控制作用,以及时补偿干扰对被控量的直接影响。

      缺点:只能对可测干扰进行补偿,对不可测干扰以及受控对象、各功能部件内部参数变化对被控量的影响,系统自身无法控制。

      在这里插入图片描述

  • 反馈/闭环控制方式

    • 基于偏差的控制

      定义:通过计算被控量和给定值的差值来控制被控对象

      特点:只要扰动和不确定因素的作用点被反馈所包围,无论扰动和不确定因素是什么形式、已知或未知的,均可有效地加以抑制。

      缺点:闭环控制系统的参数如果匹配得不好,会造成被控量的较大摆动,甚至系统无法正常工作。

      在这里插入图片描述

  • 复合控制方式

    在这里插入图片描述

对自动控制系统的基本要求

自动控制系统满足的标准或性能

  • 抑制扰动
  • 稳态误差
  • 瞬态(动态)特性
  • 对系统参数变化的灵敏度

线性系统满足叠加原理和齐次性(均匀性)

线性系统的数学描述:线性函数,如线性微分方程、线性差分方程、线性代数方程


在理想情况下,自控系统的被控量和给定值,在任何时候相等,没有误差,也不受干扰的影响,即 c ( t ) ≡ r ( t ) c(t){\equiv}r(t) c(t)r(t)

一般来说整个调节过程分为两个阶段:

  1. 动/暂态过程:反映系统的动态特性(输出量处于激烈变化之中)。
  2. 稳态过程:反映系统的稳态特性(输出量稳定在新的平衡状态并保持不变)。

系统的稳定性表现为时域响应的收敛性

稳定:设线性定常系统处于一平衡状态,若此系统在干扰作用下离开了原来的平衡状态,那么在干扰作用消失后,系统能否回到原来的平衡状态,这就是系统的稳定性问题。


工程上常用稳、快、准三个字来表达控制系统的要求

  • :指系统重新恢复平衡工作状态的能力(稳定性)和动态过程的振荡倾向(平稳性)。
  • :指动态过程进行的时间长短(快速性)。
  • :指系统过渡到新的平衡状态后或系统受扰重新恢复平衡之后(过渡过程结束),系统所保持的精度(准确性)。

总的来说,对反馈控制系统的要求:首先系统必须是稳定,其次希望实际的调节过程尽可能接近于理想的调节过程,考虑到动态过程在不同阶段中的特点,工程上常常从动态(稳、快)和静态(准)提出稳、快、准三个基本性能指标(要求)。

控制系统的数学模型

数学模型

  • 静态模型:代数方程在静态条件下(即变量各阶导数为零),描述变量之间关系的代数方程。
  • 动态模型:微分方程、差分方程、状态方程等。

控制系统属于动态系统

  • 描述连续系统运动输入输出关系的为微分方程
  • 描述离散系统运动输入输出关系的为差分方程
  • 描述输入、输出和内部各变量间关系的为状态方程
范围描述
时域微分方程(差分方程)
复数域传递函数(脉冲传递函数)
频域频率特性

微分方程

列写元件微分方程的步骤可归纳如下:

  1. 根据元件的工作原理及其在控制系统中的作用,确定其输入量和输出量。
  2. 分析元件工作中所遵循的物理规律或化学规律,列写相应的微分方程。
  3. 消去中间变量,得到输出量与输入量之间关系的微分方程,便是元件时域的数学模型。一般情况下,应将微分方程写为标准形式,即与输入量有关的项写在方程的右端,与输出量有关的项写在方程的左端方程两端变量的导数项均按降幂排列

  • 建立控制系统的微分方程时

    一般先由系统原理线路图画出系统方块图,并分别列写组成系统各元件的微分方程;然后,消去中间变量便得到描述系统输出量与输入量之间关系的微分方程。

  • 列写系统各元件的微分方程时

    • 应注意信号传送的单向性,即前一个元件的输出是后一个元件的输入,一级一级地单向传送。
    • 应注意前后连接的两个元件中,后级对前级的负载效应

微分方程的解 = 齐次微分方程的通解 + 非齐次微分方程的任一特解

用拉氏变换法求解线性定常微分方程的过程可归结如下:

  • 考虑初始条件,对微分方程中的每一项分别进行拉氏变换,将变量 t (时域)的微分方程转换为变量 s (复域)的代数方程。
  • 由代数方程求出输出量拉氏变换函数的表达式。
  • 对输出量拉氏变换函数求拉氏反变换,得到输出量时域表达式,即为所求微分方程的解。

在这里插入图片描述

小偏差法

将非线性微分方程线性化——小偏差法

  • 实际工作点在某一平衡点 (x0, y0) 附近

  • 非线性 y = f ( x ) y = f(x) y=f(x) 连续可导

  • x x x 在一个很小的范围内变化, Δ x \Delta{x} Δx 很小(小偏差)

  • 可列写线性化增量方程:

    Δ y = K Δ x \Delta{y}=K\Delta{x} Δy=KΔx K = d f ( x ) d x ∣ x = x 0 K=\frac{df(x)}{dx}|_{x=x_0} K=dxdf(x)x=x0

传递函数

线性定常系统传递函数的定义:零初始条件下,系统输出量的拉氏变换与输入量的拉氏变换之比。

  • 传递函数是复变量 s 的有理真分式函数,具有复变函数的所有性质。

  • 传递函数是一种用系统参数表示输出量与输入量之间关系的表达式它只取决于系统或元件的结构和参数,而与输入量的形式无关,也不反映系统内部的任何信息

  • 传递函数与微分方程具有相通性

    • 传递函数分子多项式系数及分母多项式系数,分别与相应微分方程的右端及左端微分算符多项式系数相对应。
    • 将微分方程的算符 d d t \frac{d}{dt} dtd 用复数 s 置换便得到传递函数;反之,将传递函数多项式中的变量s用算符 d d t \frac{d}{dt} dtd 置换便得到微分方程。

传递函数是在零初始条件下定义的。控制系统的零初始条件有两方面的含义:

  • 指输入量是在 t ≥ 0 时才作用于系统,因此,输入量及其各阶导数在 t = 0 时均为零。

  • 指输入量加于系统之前,系统处于稳定的工作状态,即输出量及其各阶导数在 t = 0 时的值也为零,现实的工程控制系统多属此类情况。因此,传递函数可表征控制系统的动态性能,并用以求出在给定输入量时系统的零初始条件响应,即由拉氏变换的卷积定理,有:
    c ( t ) = L − 1 [ C ( s ) ] = L − 1 [ G ( s ) R ( s ) ] = ∫ 0 t r ( τ ) g ( t − τ ) d τ = ∫ 0 t r ( t − τ ) g ( τ ) d τ c(t)=L^{-1}[C(s)]=L^{-1}[G(s)R(s)]=\int_{0}^{t}r(\tau)g(t-\tau)d\tau =\int_{0}^{t}r(t-\tau)g(\tau)d\tau c(t)=L1[C(s)]=L1[G(s)R(s)]=0tr(τ)g(tτ)dτ=0tr(tτ)g(τ)dτ

    式中, g ( t ) = L − 1 [ G ( s ) ] g(t)=L^{-1}[G(s)] g(t)=L1[G(s)] 是系统的脉冲响应。

  • 传递函数 G ( s ) G(s) G(s) 的拉氏反变换是脉冲响应 g ( t ) g(t) g(t)

  • 脉冲响应(也称脉冲过渡函数) g ( t ) g(t) g(t) 是系统在单位脉冲 δ ( t ) \delta(t) δ(t) 输入时的输出响应


传递函数的表达形式

  • 一般表达形式
    G ( s ) = C ( s ) R ( s ) = b 0 s m + b 1 s m − 1 + . . . + b m − 1 s + b m a 0 s n + a 1 s n − 1 + . . . + a n − 1 s + a n = M ( s ) N ( s ) G(s)=\frac{C(s)}{R(s)}=\frac{b_0s^m+b_1s^{m-1}+...+b_{m-1}s+b_m}{a_0s^n+a_1s^{n-1}+...+a_{n-1}s+a_n}=\frac{M(s)}{N(s)} G(s)=R(s)C(s)=a0sn+a1sn1+...+an1s+anb0sm+b1sm1+...+bm1s+bm=N(s)M(s)

  • 零-极点表达形式——最高阶次项系数为1
    G ( s ) = b 0 ( s − z 1 ) ( s − z 2 ) . . . ( s − z m ) a 0 ( s − p 1 ) ( s − p 2 ) . . . ( s − p n ) = K ∗ ∏ i = 1 m ( s − z i ) ∏ j = 1 n ( s − p j ) G(s)=\frac{b_0(s-z_1)(s-z_2)...(s-z_m)}{a_0(s-p_1)(s-p_2)...(s-p_n)}=K^*\frac{\prod_{i=1}^{m}(s-z_i)}{\prod_{j=1}^{n}(s-p_j)} G(s)=a0(sp1)(sp2)...(spn)b0(sz1)(sz2)...(szm)=Kj=1n(spj)i=1m(szi)
    系数 K ∗ = b 0 / a 0 K^*=b_0/a_0 K=b0/a0 称为传递系数或根轨迹增益。这种用零点和极点表示传递函数的方法在根轨迹法中使用较多。

    在复数平面上表示传递函数的零点和极点的图形,称为传递函数的零极点分布图。在图中一般用 “o” 表示零点,用 “×” 表示极点。传递函数零极点分布图可以更形象地反映系统的全面特性。

  • 典型环节形式——常数项系数为1
    G ( s ) = b m ( τ 1 s + 1 ) ( τ 2 s + 1 ) . . . ( τ m s + 1 ) a n ( T 1 s + 1 ) ( T 2 s + 1 ) . . . ( T n s + 1 ) = K ∏ i = 1 m ( τ i s + 1 ) ∏ j = 1 n ( T j s + 1 ) G(s)=\frac{b_m(\tau_1s+1)(\tau_2s+1)...(\tau_ms+1)}{a_n(T_1s+1)(T_2s+1)...(T_ns+1)}=K\frac{\prod_{i=1}^{m}(\tau_is+1)}{\prod_{j=1}^{n}(T_js+1)} G(s)=an(T1s+1)(T2s+1)...(Tns+1)bm(τ1s+1)(τ2s+1)...(τms+1)=Kj=1n(Tjs+1)i=1m(τis+1)
    K K K 称传递系数或增益。传递函数这种表示形式在频率法中使用较多。


传递函数的极点就是微分方程的特征根,极点/特征根决定了所描述系统自由运动(零输入响应)的模态

传递函数的零点并不形成自由运动的模态,但它们却影响各模态响应中所占的比重(系数),因而也影响响应曲线的形状。

系统产生的运动(响应)分为:

  • 强迫运动(零初始条件响应)。

  • 自由运动(零输入响应),是系统“固有”的成分。

典型环节

  • 比例环节

  • 积分环节

  • 惯性环节(非周期环节)

  • 振荡环节
    G ( s ) = C ( s ) R ( s ) = w n 2 s 2 + 2 ζ w n s + w n 2 = 1 T 2 s 2 + 2 ζ T s + 1 G(s)=\frac{C(s)}{R(s)}=\frac{{w_n}^2}{s^2+2\zeta{w_n}s+{w_n}^2 }=\frac{1}{T^2s^2+2\zeta{T}s+1} G(s)=R(s)C(s)=s2+2ζwns+wn2wn2=T2s2+2ζTs+11

  • 微分环节

    • 理想微分环节(与积分环节对应)
    • 一阶微分环节(与惯性环节对应)
    • 二阶微分环节(与振荡环节对应)
  • 延迟环节

    • 时域表示 c ( t ) = r ( t − τ ) l ( t − τ ) c(t)=r(t-\tau)l(t-\tau) c(t)=r(tτ)l(tτ)
    • 复域表示 G ( s ) = C ( s ) R ( s ) = e − τ s G(s)=\frac{C(s)}{R(s)}=e^{-\tau{s}} G(s)=R(s)C(s)=eτs

控制系统的结构图

控制系统的结构图是由许多对信号进行单向运算的方框和一些信号流向线组成

它包含四种基本单元:

  • 信号线: 信号线是带有箭头的直线,箭头表示信号流向,在直线旁标记信号的时间函数或象函数。
  • 引出点(或测量点): 引出点表示信号引出或测量的位置,从同一位置引出的信号在数值和性质上完全相同。
  • 比较点(或综合点): 比较点表示对两个以上的信号进行加减运算,+号表示相加,-号表示相减,+号可省不写
  • 方框(或环节): 方框表示对信号进行的数学变换,方框中写入元部件/系统传递函数;方框的输出变量等于方框的输入变量与传递函数的乘积,即 C(s) = G(s)U(s)。因此,方框可视作单向运算的算子

系统结构图也是控制系统的一种数学模型。

建议:第1个方程最好能包含输入,相加点,输出


信号流图与结构图的关系

  • 节点表示系统的变量。
  • 支路相当于乘法器。

信号流图的相关概念

  • 前向通路:信号从输入节点到输出节点传递时,每个节点只通过一次的通路,叫前向通路。
  • 前向通路总增益:前向通路上各支路增益之乘积,称前向通路总增益,一般用 p k p_k pk 表示。
  • 回路:起点和终点在同一节点,而且信号通过每一节点不多于一次的闭合通路称为单独回路,简称回路。
  • 回路增益:回路中所有支路增益之乘积叫回路增益,用 L a L_a La 表示。
  • 不接触回路:回路与回路之间没有公共节点时,这种回路叫做不接触回路。可以有两个或两个以上不接触的回路。
  • 不接触回路增益 L b L c L_bL_c LbLc 为所有互不接触的单独回路中,每次取其中2个回路的回路增益的乘积; L d L e L f L_dL_eL_f LdLeLf为所有互不接触的单独回路中,每次取其中3个回路的回路增益的乘积。

梅森增益公式

求取从任意源节点到任意阱节点之间传递函数的梅森增益公式记为:
P = 1 Δ ∑ k = 1 n p k Δ k P=\frac{1}{\Delta}\sum_{k=1}^{n}p_k\Delta_k P=Δ1k=1npkΔk

  • P P P 为从源节点到阱节点的传递函数(或总增益)。
  • n n n 为从源节点到阱节点的前向通路总数。
  • p k p_k pk 为从源节点到阱节点的第k条前向通路总增益。
  • Δ \Delta Δ 为流图特征式, Δ = 1 − ∑ L a + ∑ L b L c − ∑ L d L e L f + … \Delta=1-\sum{L_a}+\sum{L_bL_c}-\sum{L_dL_eL_f} +… Δ=1La+LbLcLdLeLf+
    • ∑ L a \sum{L_a} La 为所有单独回路增益之和。
    • ∑ L b L c \sum{L_bL_c} LbLc 为所有互不接触的单独回路中,每次取其中两个回路的回路增益的乘积之和。
    • ∑ L d L e L f \sum{L_dL_eL_f} LdLeLf 为所有互不接触的单独回路中,每次取其中三个回路的回路增益的乘积之和。
  • Δ k \Delta_k Δk 为流图特征式的余因子式,它等于流图特征式 Δ \Delta Δ 中除去与第 k 条前向通路相接触的回路增益项(包括回路增益的乘积项)以后的余项式。

  • 为了研究有用输入作用对系统输出 C ( s ) C(s) C(s) 的影响,需要求有用输入作用下的闭环传递函数 C ( s ) / R ( s ) C(s)/R(s) C(s)/R(s)
  • 为了研究扰动作用 N ( s ) N(s) N(s) 对系统输出 C ( s ) C(s) C(s) 的影响,也需要求取扰动作用下的闭环传递函数 C ( s ) / N ( s ) C(s) / N(s) C(s)/N(s)
  • 此外,在控制系统的分析和设计中,还常用到在输入信号 R ( s ) R(s) R(s) 或扰动 N ( s ) N(s) N(s) 作用下,以误差信号 E ( s ) E(s) E(s) 作为输出量的闭环误差传递函数 E ( s ) / R ( s ) E(s)/R(s) E(s)/R(s) E ( s ) / N ( s ) E(s)/N(s) E(s)/N(s)

在一定条件下,系统的输出只取决于反馈通路传递函数 H ( s ) H(s) H(s) 及输入信号 R ( s ) R(s) R(s),既与前向通路传递函数无关,也不受扰动作用的影响。特别是当 H ( s ) = 1 H(s)=1 H(s)=1,即单位反馈时, C ( s ) ≈ R ( s ) C(s){\approx}R(s) C(s)R(s),从而近似实现了对输入信号的完全复现,且对扰动具有较强的抑制能力。

绝不允许将各种闭环传递函数进行叠加后求其输出响应

线性系统的时域分析法

线性系统的稳定性分析

线性系统的稳定性仅取决于系统自身的固有特性,而与外界条件无关。因此,设线性系统在初始条件为零时,作用一个理想单位脉冲,这时系统的输出增量为脉冲响应 c ( t ) c(t) c(t)。这相当于系统在扰动信号作用下,输出信号偏离原平衡工作点的问题。若 t → ∞ t\to\infty t时,脉冲响应:
lim ⁡ t → ∞ c ( t ) = 0 \lim_{t\to\infty}c(t)=0 tlimc(t)=0
即输出增量收敛于原平衡工作点,则线性系统是稳定的。


根据李雅普诺夫稳定性理论,线性控制系统的稳定性可叙述如下:

若线性控制系统在初始扰动的影响下,其动态过程随时间的推移逐渐衰减并趋于零(原平衡工作点),则称系统渐近稳定,简称稳定;反之,若在初始扰动影响下,系统的动态过程随时间的推移而发散,则称系统不稳定。


线性系统稳定的充分必要条件:闭环系统特征方程的所有根均具有负实部;或者说,闭环传递函数的极点均严格位于左半 s 平面 。

若特征根中具有一个或一个以上零实部根,而其余的特征根均具有负实部,则脉冲响应 c ( t ) c(t) c(t) 趋于常数,或趋于等幅正弦振荡,按照稳定性定义,此时系统不是渐近稳定的,处于稳定和不稳定的临界状态,常称为临界稳定情况


在经典控制理论中,只有渐近稳定的系统才称为稳定系统;否则,称为不稳定系统

线性系统稳定的充分且必要条件:劳思表中第一列各值为正。

如果劳思表第一列中出现小于零的数值,系统不稳定,且第一列各系数符号的改变次数,代表特征方程的正实部根的数目

线性系统的稳态误差计算

  • 控制系统的稳态误差,是系统控制准确度(控制精度)的一种度量,通常称为稳态性能指标。
  • 在控制系统设计中,稳态误差是一项重要的技术指标。
  • 对于一个实际的控制系统,由于系统结构、输入作用的类型(控制量或扰动量)、输入函数的形式(阶跃、斜坡或加速度)不同,控制系统的稳态输出不可能在任何情况下都与输入量一致,也不可能在任何形式的扰动作用下都能准确地恢复到原平衡位置。
  • 此外,控制系统中不可避免地存在摩擦、间隙、不灵敏区、零位输出等非线性因素,都会造成附加的稳态误差。
  • 可以说,控制系统的稳态误差是不可避免的,控制系统设计的任务之一,是尽量减小系统的稳态误差,或使稳态误差小于某一容许值。

只有当系统稳定时,研究稳态误差才有意义


把在阶跃函数作用下没有原理性稳态误差的系统,称为无差系统。

而把具有原理性稳态误差的系统,称为有差系统。

非线性因素所引起的系统稳态误差,则称为附加稳态误差,或结构性稳态误差。


相关概念

  • 稳态:时间 t 趋于无穷大时系统的输出状态。
  • 误差:系统实际输出与期望输出之差。
  • 稳态误差:对于稳定的系统,暂态响应随时间的增长而衰减,当衰减到可以忽略的程度时,则希望的稳态响应与实际的稳态响应之差称为稳态误差。稳态误差是系统达到稳态时,系统精度的度量。

稳态误差定义
e s s ( ∞ ) = lim ⁡ t → ∞ ε ( t ) = lim ⁡ t → ∞ ( c 0 ( t ) − c ( t ) ) e_{ss}(\infty)=\lim_{t\to\infty}\varepsilon(t)=\lim_{t\to\infty}(c_0(t)-c(t)) ess()=tlimε(t)=tlim(c0(t)c(t))
因系统的输入可分为给定和扰动输入,所以稳态误差也可分为给定稳态误差和扰动稳态误差

  • 对于随动系统,给定的参考输入是变化的,要求输出严格跟随输入(给定)的变化,则其响应的希望值就是给定的参考输入值。所以,一般以系统的给定稳态误差去衡量随动系统的稳态性能。
  • 对于恒值调节系统,一般给定不变,主要考虑扰动的作用。所以,一般以扰动稳态误差去衡量恒值调节系统的稳态性能。

  • 误差:是从系统输出端来定义,它定义为系统输出量的希望值与实际值之差:
    ε ( s ) = C 0 ( s ) − C ( s ) \varepsilon(s)=C_0(s)-C(s) ε(s)=C0(s)C(s)

  • 偏差:是在系统输入端来定义,定义为系统给定量与反馈量之差:
    E ( s ) = R ( s ) − B ( s ) = H ( s ) ε ( s ) E(s)=R(s)-B(s)=H(s)\varepsilon(s) E(s)=R(s)B(s)=H(s)ε(s)

偏差值与误差值的差别为H(s),对于单位反馈系统,偏差与误差的值是一样的。

一般地用偏差代替误差:
ε s s ( ∞ ) = e s s ( ∞ ) = lim ⁡ t → ∞ e ( t ) = lim ⁡ t → ∞ ( r ( t ) − b ( t ) ) \varepsilon_{ss}(\infty)=e_{ss}(\infty)=\lim_{t\to\infty}e(t)=\lim_{t\to\infty}(r(t)-b(t)) εss()=ess()=tlime(t)=tlim(r(t)b(t))

  • 在误差信号 e ( t ) e(t) e(t) 中,包含瞬态分量 e t s ( t ) e_{ts}(t) ets(t) 和稳态分量 e s s ( t ) e_{ss}(t) ess(t) 两部分。
  • 由于系统必须稳定,故当时间趋于无穷时,必有 e t s ( t ) e_{ts}(t) ets(t) 趋于零。
  • 因而,控制系统的稳态误差(终值)定义为:误差信号 e ( t ) e(t) e(t) 的稳态分量 e s s ( ∞ ) e_{ss}(\infty) ess(),常以 e s s e_{ss} ess 简单标志。

如果有理函数 s E ( s ) sE(s) sE(s) 除在原点处有惟一的极点外,在 s 右半平面及虚轴上处处解析(即不存在极点),即 s E ( s ) sE(s) sE(s) 的极点均位于 s 左半平面(包括坐标原点),则可根据拉氏变换的终值定理,方便地求出系统的稳态误差:
e s s ( ∞ ) = lim ⁡ t → ∞ e ( t ) = lim ⁡ s → 0 s E ( s ) e_{ss}(\infty)=\lim_{t\to\infty}e(t)=\lim_{s\to0}sE(s) ess()=tlime(t)=s0limsE(s)
上式算出的稳态误差是误差信号 e ( t ) e(t) e(t) 的稳态分量 e s s ( t ) e_{ss}(t) ess(t) 在 t 趋于无穷时的数值,故有时称为终值误差或稳态误差终值,它不能反映 e ( t ) e(t) e(t) 时间 t 的完整变化规律,具有一定的局限性

当输入信号形式一定时,系统是否存在稳态误差就取决于开环传递函数描述的系统结构。

在一般情况下,分子阶次为 m,分母阶次为 n 的开环传递函数 G k G_k Gk 可表示为:
G ( s ) = G ( s ) H ( s ) = k ∏ i = 1 m ( τ i s + 1 ) s v ∏ j = 1 n ( T j s + 1 ) ( n ≥ m ) G(s)=G(s)H(s)=\frac{k\prod_{i=1}^{m}(\tau_is+1)}{s^v\prod_{j=1}^{n}(T_js+1)}(n\ge{m}) G(s)=G(s)H(s)=svj=1n(Tjs+1)ki=1m(τis+1)(nm)
式中, K K K 为开环增益; τ i \tau_i τi T j T_j Tj 为时间常数; v v v 为开环系统在 s 平面坐标原点上的极点的重数。现在的分类方法是以 v v v的数值来划分的: v = 0 v= 0 v=0 ,称为0型系统; v = 1 v=1 v=1 ,称为 I 型系统; v = 2 v=2 v=2 ,称为 II 型系统。当 v > 2 v>2 v>2 时,除复合控制系统外,使系统稳定是相当困难的。因此除航天控制系统外,III 型及III型以上的系统几乎不采用。

这种以开环系统在 s 平面坐标原点上的极点数来分类的方法,其优点在于:可以根据已知的输入信号形式迅速判断系统是否存在原理性稳态误差及稳态误差的大小。它与按系统的阶次进行分类的方法不同,阶次 m 与 n 的大小与系统的型别无关,且不影响稳态误差的数值。

误差系数

系统的静态位置误差系数 K p K_p Kp k p = lim ⁡ s → 0 G k ( s ) k_p=\lim_{s\to0}G_k(s) kp=lims0Gk(s)

系统的静态速度误差系数 K v K_v Kv k v = lim ⁡ s → 0 s G k ( s ) k_v=\lim_{s\to0}sG_k(s) kv=lims0sGk(s)

系统的静态加速度误差系数 K a K_a Ka k a = lim ⁡ s → 0 s 2 G k ( s ) k_a=\lim_{s\to0}s^2G_k(s) ka=lims0s2Gk(s)

如果系统为非单位反馈系统,其 H ( s ) = K h H(s)=K_h H(s)Kh 为常数,那么系统输出量的希望值为 R ′ s ) = R ( s ) / K h R's)=R(s)/K_h Rs)=R(s)/Kh,系统输出端的稳态位置误差为: e s s ′ = e s s / K h e'_{ss}=e_{ss}/K_h ess=ess/Kh

在这里插入图片描述


给定稳态误差级数 (动态误差系数)

利用动态误差系数法,可以研究输入信号几乎为任意时间函数时的系统稳态误差变化,因此动态误差系数又称广义误差系数。为了求取动态误差系数,写出稳态误差信号的拉氏变换式: E ( s ) = Φ e ( s ) R ( s ) E(s)=\Phi_e(s)R(s) E(s)=Φe(s)R(s)

将误差传递函数在 Φ e ( s ) = 0 \Phi_e(s)=0 Φe(s)0 的邻域内展成泰勒级数,得:
Φ e ( s ) = 1 1 + G ( s ) H ( s ) = Φ e ( 0 ) + Φ ˙ e ( 0 ) s + 1 2 ! Φ ¨ e ( 0 ) s 2 + ⋯ \Phi_e(s)=\frac{1}{1+G(s)H(s)}=\Phi_e(0)+\dot\Phi_e(0)s+\frac{1}{2!}\ddot\Phi_e(0)s^2+\cdots Φe(s)=1+G(s)H(s)1=Φe(0)+Φ˙e(0)s+2!1Φ¨e(0)s2+
于是
E ( s ) = Φ e ( 0 ) R ( s ) + Φ ˙ e ( 0 ) s R ( s ) + 1 2 ! Φ ¨ e ( 0 ) s 2 R ( s ) + ⋯ + 1 l ! Φ e ( l ) ( 0 ) s l R ( s ) + ⋯ E(s)=\Phi_e(0)R(s)+\dot\Phi_e(0)sR(s)+\frac{1}{2!}\ddot\Phi_e(0)s^2R(s)+\cdots+\frac{1}{l!}\Phi_e^{(l)}(0)s^lR(s)+\cdots E(s)=Φe(0)R(s)+Φ˙e(0)sR(s)+2!1Φ¨e(0)s2R(s)++l!1Φe(l)(0)slR(s)+
上述无穷级数收敛于 s=0 的邻域,称为给定输入作用下的稳态误差级数,相当于在时间域内 t → ∞ t\to\infty t 时成立。因此,当所有初始条件均为零时,对上式进行拉氏反变换,就得到作为时间函数的稳态误差表达式:
e s s ( t ) = ∑ i = 0 ∞ C i r ( i ) ( t ) C i = 1 i ! Φ e ( i ) ( 0 ) i = 0 , 1 , 2 , ⋯ e_{ss}(t)=\sum_{i=0}^{\infty}C_ir^{(i)}(t)\\ C_i=\frac{1}{i!}\Phi_e^{(i)}(0)\\ i=0,1,2,\cdots ess(t)=i=0Cir(i)(t)Ci=i!1Φe(i)(0)i=0,1,2,

减小稳态误差

为了减小或消除系统在输入信号和扰动作用下的稳态误差,可以采取以下措施:

  • 增大系统开环增益或扰动作用点之前系统的前向通道增益

    增大系统开环增益 K K K 以后,对于 0 型系统,可以减小系统在阶跃输入时的位置误差;对于 I 型系统,可以减小系统在斜坡输入时的速度误差;对于 Ⅱ 型系统,可以减小系统在加速度输入时的加速度误差。

    增大扰动点之后系统的前向通道增益,不能改变系统对扰动的稳态误差数值。

  • 在系统的前向通道或主反馈通道设置串联积分环节

    扰动作用点之前的前向通道积分环节数与主反馈通道积分环节数之和决定系统响应扰动作用的型别,该型别与扰动作用点之后前向通道的积分环节数无关。

    反馈控制系统中,设置串联积分环节或增大开环增益以消除或减小稳态误差的措施,必然导致降低系统的稳定性,甚至造成系统不稳定,从而恶化系统的动态性能

  • 采用串级控制抑制内回路扰动

    当控制系统中存在多个扰动信号,且控制精度要求较高时,宜采用串级控制方式,可以显著抑制内回路的扰动影响。

  • 采用复合控制方法

    复合控制系统是在系统的反馈控制回路中加入前馈通路,组成一个前馈控制与反馈控制相结合的系统。

一阶系统的时域分析

线性定常系统在任意输入信号作用下的时间响应,均可分解为零输入响应与零状态响应之和
y ( t ) = y z i + y z s y(t)=y_{zi}+y_{zs} y(t)=yzi+yzs

  • y z i y_{zi} yzi零输入响应

    • 只与系统的极点类型和分布有关,与系统的零点无关,与系统的初始条件有关。

    • 反映系统的稳定性。

  • y z s y_{zs} yzs零状态响应

    • 在时域:表示为系统的单位脉冲响应与输入信号的卷积。
    • 在复数域:表示为传递函数与输入信号拉氏变换的乘积。
    • 与系统的零极点都有关,还与输入信号的性质有关。
    • 反映系统的稳定性、暂态(快速及平稳性)及稳态(准确)特性。

数学模型

  • 微分方程 T d c ( t ) d t + c ( t ) = r ( t ) T\frac{dc(t)}{dt}+c(t)=r(t) Tdtdc(t)+c(t)=r(t)
  • 开环传函 G k ( s ) = 1 T s G_k(s)=\frac{1}{Ts} Gk(s)=Ts1
  • 闭环传函 Φ ( s ) = C ( s ) R ( s ) = 1 T s + 1 \Phi(s)=\frac{C(s)}{R(s)}=\frac{1}{Ts+1} Φ(s)=R(s)C(s)=Ts+11
  • 特征方程 T s + 1 = 0 Ts+1=0 Ts+1=0

输入响应

  • 单位阶跃响应

    t = 3 T , c ( t ) = 0.95 t=3T,c(t)=0.95 t=3T,c(t)=0.95 —— 5% 误差带

    t = 4 T , c ( t ) = 0.98 t=4T,c(t)=0.98 t=4T,c(t)=0.98 —— 2% 误差带

    • 可以用时间常数 T 去度量系统的输出量的数值

    • 初始斜率为 1/T

    • 无超调

    • 稳态误差 e s s = 0 e_{ss}=0 ess=0

  • 单位脉冲响应

    • 可以用时间常数去度量系统的输出量的数值
    • 初始斜率为 − 1 / T 2 -1/T^2 1/T2
    • 无超调
    • 稳态误差 e s s = 0 e_{ss}=0 ess=0
  • 单位斜坡响应

    出现稳态误差( e s s = T e_{ss}=T ess=T

    一阶系统不能实现对加速度输入函数的跟踪

在这里插入图片描述

  • 系统对输入信号导数的响应,就等于系统对该输入信号响应的导数。
  • 系统对输入信号积分的响应,就等于系统对该输入信号响应的积分,积分常数由零输出初始条件确定。

二阶系统的时域分析

  • 微分方程

T 2 d 2 c ( t ) d t 2 + 2 ζ T d c ( t ) d t + c ( t ) = r ( t ) d 2 c ( t ) d t 2 + 2 ζ w n d c ( t ) d t + w n 2 c ( t ) = w n 2 r ( t ) T^2\frac{d^2c(t)}{dt^2}+2\zeta{T}\frac{dc(t)}{dt}+c(t)=r(t)\\ \frac{d^2c(t)}{dt^2}+2\zeta{w_n}\frac{dc(t)}{dt}+w^2_nc(t)=w^2_nr(t) T2dt2d2c(t)+2ζTdtdc(t)+c(t)=r(t)dt2d2c(t)+2ζwndtdc(t)+wn2c(t)=wn2r(t)

  • 开环传函

G k ( s ) = 1 T 2 s 2 + 2 ζ T s = w n 2 s 2 + 2 ζ w n s G_k(s)=\frac{1}{T^2s^2+2\zeta{T}s}=\frac{w^2_n}{s^2+2\zeta{w_n}s} Gk(s)=T2s2+2ζTs1=s2+2ζwnswn2

  • 闭环传函

Φ ( s ) = C ( s ) R ( s ) = w n 2 s 2 + 2 ζ w n s + w n 2 = 1 T 2 s 2 + 2 ζ T s + 1 \Phi(s)=\frac{C(s)}{R(s)}=\frac{{w_n}^2}{s^2+2\zeta{w_n}s+{w_n}^2 }=\frac{1}{T^2s^2+2\zeta{T}s+1} Φ(s)=R(s)C(s)=s2+2ζwns+wn2wn2=T2s2+2ζTs+11

  • 特征方程的根

s 1 , 2 = − ζ w n ∓ j w n 1 − ζ 2 s_{1,2}=-\zeta{w_n}\mp{jw_n\sqrt{1-\zeta^2}} s1,2=ζwnjwn1ζ2


阻尼比

阻尼比所属情况
ζ = 0 \zeta=0 ζ=0无(零)阻尼
0 < ζ < 1 0<\zeta<1 0<ζ<1欠阻尼
ζ = 1 \zeta=1 ζ1临界阻尼(重极点)
ζ > 1 \zeta>1 ζ>1过阻尼
  • 由于阻尼比 ζ \zeta ζ 为负,指数因子具有正幂指数,因此系统的动态过程为发散正弦振荡或单调发散的形式,从而表明 ζ < 0 \zeta<0 ζ<0 的二阶系统是不稳定的。

  • 如果 ζ = 0 \zeta=0 ζ0,则特征方程有一对纯虚根, s 1 , 2 = ± j w n s_{1,2}=\pm{jw_n} s1,2=±jwn,对应于 s 平面虚轴上一对共轭极点,可以算出系统的阶跃响应为等幅振荡,此时系统相当于无阻尼情况。

  • 如果 0 < ζ < 1 0<\zeta<1 0<ζ<1,则特征方程有一对具有负实部的共轭复根, s 1 , 2 = − ζ w n ∓ j w n 1 − ζ 2 s_{1,2}=-\zeta{w_n}\mp{jw_n\sqrt{1-\zeta^2}} s1,2=ζwnjwn1ζ2 ,对应于 s 平面左半部的共轭复数极点,相应的阶跃响应为衰减振荡过程,此时系统处于欠阻尼情况。

    • 系统特征方程既有实部也有虚部,其单位阶跃响应既有振荡成分也有衰减成分,是一个随时间 t 的增长,振幅按指数规律衰减的周期函数(衰减振荡曲线)。
  • 如果 ζ = 1 \zeta=1 ζ1,则特征方程具有两个相等的负实根, s 1 , 2 = − w n s_{1,2}=-w_n s1,2=wn,对应于s平面负实轴上的两个相等实极点,相应的阶跃响应非周期地趋于稳态输出,此时系统处于临界阻尼情况。

  • 如果 ζ > 1 \zeta>1 ζ>1,则特征方程有两个不相等的负实根 s 1 , 2 = − ζ w n ± w n ζ 2 − 1 s_{1,2}=-\zeta{w_n}\pm{w_n\sqrt{\zeta^2-1}} s1,2=ζwn±wnζ21 ,对应于 s 平面负实轴上的两个不等实极点,相应的单位阶跃响应也是非周期地趋于稳态输出,但响应速度比 1 临界阻尼情况缓慢,因此称为过阻尼情况。

二阶控制系统的设计,一般取 ζ = 0.4 ∼ 0.8 \zeta=0.4\sim0.8 ζ=0.40.8


二阶系统的单位阶跃响应

  • 根的实部反应在响应的指数部分(模态形状:是否衰减)

  • 根的虚部反应在响应的正弦部分(模态形状:是否振荡)

    • ζ = 0 \zeta=0 ζ=0 时为临界稳定状态

    • ζ ≈ 0.707 \zeta ≈ 0.707 ζ0.707 时调节时间最短(称为最佳阻尼比)

s 1 , 2 = − ζ w n ± j w n 1 − ζ 2 = − σ ± j w d s_{1,2}=-\zeta{w_n}\pm{jw_n\sqrt{1-\zeta^2}}=-\sigma\pm{jw_d} s1,2=ζwn±jwn1ζ2 =σ±jwd

  • 衰减系数 σ \sigma σ (实部)是闭环极点到虚轴之间的距离

  • 阻尼振荡频率 w d w_d wd (虚部)是闭环极点到实轴之间的距离: w d = w n ( 1 − ζ 2 ) w_d=w_n(\sqrt{1-\zeta^2}) wd=wn(1ζ2 )

  • 自然频率 w n w_n wn 是闭环极点到坐标原点之间的距离

    w n w_n wn 与负实轴夹角的余弦正好是阻尼比 ζ \zeta ζ,即 cos ⁡ φ = ζ \cos\varphi=\zeta cosφ=ζ

结论:调节时间与闭环极点的实部数值成反比。闭环极点距虚轴的距离越远,系统的调节时间越短。由于阻尼比值主要根据对系统超调量的要求来确定,所以调节时间主要由自然频率决定。若能保持阻尼比值不变而加大自然频率值,则可在不改变超调量的情况下缩短调节时间。

欠阻尼二阶系统阶跃响应

  • 峰值时间 t p = π w d = π w n 1 − ζ 2 t_p=\frac{\pi}{w_d}=\frac{\pi}{w_n\sqrt{1-\zeta^2}} tp=wdπ=wn1ζ2 π 与根(极点)的虚部成反比

  • 超调量 σ p = e − π ζ 1 − ζ 2 × 100 % \sigma_p=e^{-\frac{\pi\zeta}{\sqrt{1-\zeta^2}}}\times100\% σp=e1ζ2 πζ×100%

  • 上升时间 t r = π − φ w n 1 − ζ 2 t_r=\frac{\pi-\varphi}{w_n\sqrt{1-\zeta^2}} tr=wn1ζ2 πφ 为第一次达到稳态值时的时间

  • 调节时间 t s = 3.5 ζ w n ( Δ = 5 % ) t_s=\frac{3.5}{\zeta{w_n}}(\Delta=5\%) ts=ζwn3.5(Δ=5%) t s = 4.4 ζ w n ( Δ = 2 % ) t_s=\frac{4.4}{\zeta{w_n}}(\Delta=2\%) ts=ζwn4.4(Δ=2%) 与闭环极点的实部数值成反比

  • 延迟时间 t d ≈ 1 + 0.7 ζ w n t_d\approx\frac{1+0.7\zeta}{w_n} tdwn1+0.7ζ t d ≈ 1 + 0.6 ζ + 0.2 ζ 2 w n t_d\approx\frac{1+0.6\zeta+0.2\zeta^2}{w_n} tdwn1+0.6ζ+0.2ζ2 为第一次达到稳态值一半的时间,与震荡频率成反比,与 ξ 成正比

由于阻尼比值主要根据对系统超调量的要求来确定,所以调节时间主要由自然频率决定。

对于欠阻尼二阶系统,极点的阻尼角(阻尼比)决定响应的平稳性;阻尼比(阻尼角)一定时,极点与虚轴的距离( ζ w n \zeta{w_n} ζwn)决定响应的快速性。

欠阻尼二阶系统单位阶跃响应稳态误差为零

添加零点对典型二阶系统暂态响应特性影响

  • 添加闭环零点

    • 输出响应的变化率越大微分作用便越强,零点的影响就越大;闭环零点离虚轴越近,影响就越显著,若零点离虚轴越远则影响就越弱。
    • 一般而言,添加闭环零点,使响应加快,震荡加剧,超调增大。零点越靠近虚轴,作用越明显。
  • 添加开环零点

    增大阻尼比,不影响系统的稳定性,系统的动态性能可以得到改善,不改变系统的稳态精度。

    • 比例微分控制(在前向通道中添加零点)不影响系统的稳定性。
    • 比例微分控制(添加开环零点)可以增大系统的阻尼比,从而使阶跃响应的超调量下降、系统平稳性增强;可以使闭环增加同样的零点,从而使系统响应加快;比例微分控制(添加开环零点)可以全面改善系统的动态响应特性。
    • 比例( k p = 1 k_p=1 kp=1)微分控制(添加开环零点)不影响系统的稳态精度。
    • 比例微分控制(添加开环零点)主要用于改善系统的暂态响应(动态)性能

速度反馈(微分负反馈)

输出量的导数同样可以用来改善系统的性能。通过将输出的速度信号反馈到系统输入端,并与误差信号比较,其效果与比例一微分控制相似,可以增大系统阻尼,改善系统动态性能。

  • 微分负反馈控制器不影响系统的稳定性。
  • 微分负反馈控制器可以增大系统的阻尼比,从而使阶跃响应的超调量下降、系统平稳性增强、使系统调节时间加快;微分负反馈控制器可以改善系统的动态响应特性。
  • 微分负反馈控制器降低系统的稳态精度
  • 微分负反馈控制器主要用于改善系统的暂态响应(动态)性能,但会增大稳态误差。
  • 为了减小稳态误差,必须加大原系统的开环增益,而使 K t K_t Kt 单纯用来增大系统阻尼

比例-微分控制与测速反馈控制的比较

  • 附加阻尼来源:比例-微分控制的阻尼作用产生于系统的输入端误差信号的速度,而测速反馈控制的阻尼作用来源于系统输出端响应的速度,因此对于给定的开环增益和指令输入速度,后者对应较大的稳态误差值。
  • 使用环境比例-微分控制对噪声有明显放大作用,当系统输入端噪声严重时,不宜选用比例-微分控制。同时,微分器的输入信号为系统误差信号,其能量水平低,需要当大的放大作用,为了不明显恶化信噪比,要求选用高质量的放大器;而测速反馈控制对系统输入端噪声有滤波作用,同时测速发电机的输入信号能量水平较高,因此对系统组成元件没有过高的质量要求,使用场合比较广泛。
  • 对开环增益和自然频率的影响:比例-微分控制对系统的开环增益和自然频率均无影响;测速反馈控制虽不影响自然频率,但却会降低开环增益。因此,对于确定的常值稳态误差,测速反馈控制要求有较大的开环增益。开环增益的加大,必然导致系统自然频率增大,在系统存在高频噪声时,可能引起系统共振。
  • 对动态性能的影响:比例-微分控制相当于在系统中加入实零点,可以加快上升时间。在相同阻尼比的条件下,比例微分控制系统的超调量会大于测速反馈控制系统的超调量。

线性系统的根轨迹法

根轨迹简称根迹,它是开环系统某一参数从0变到 ∞ \infty 时,闭环系统特征方程式的根在 s 平面上变化的轨迹。

根轨迹的分类:180度根轨迹 、参数根轨迹、零度根轨迹、根轨迹簇。


当闭环系统没有零点与极点相消时,闭环特征方程式的根就是闭环传递函数的极点,我们常简称为闭环极点。

从已知的开环零、极点位置及某一变化的参数来求取闭环极点的分布,实际上就是解决闭环特征方程式的求根问题。

因为系统的稳定性由系统闭环极点唯一确定,而系统的稳态性能和动态性能又与闭环零、极点在 s 平面上的位置密切相关,所以根轨迹图不仅可以直接给出闭环系统时间响应的全部信息,而且可以指明开环零、极点应该怎样变化才能满足给定的闭环系统的性能指标要求。


根轨迹法应用前提条件是受控系统的传递函数没有零极点相消

  • 闭环系统根轨迹增益,等于开环系统前向通路根轨迹增益。对于单位反馈系统,闭环系统根轨迹增益就等于开环系统根轨迹增益。
  • 闭环零点由开环前向通路传递函数的零点和反馈通路传递函数的极点所组成。对于单位反馈系统,闭环零点就是开环零点。
  • 闭环极点与开环零点、开环极点以及开环根轨迹增益均有关。

根轨迹是系统所有闭环极点的集合
根轨迹方程  K ∗ ∏ j = 1 m ( s − z j ) ∏ i = 1 n ( s − p i ) = − 1 \text{根轨迹方程 }K^*\frac{\prod_{j=1}^{m}(s-z_j)}{\prod_{i=1}^{n}(s-p_i)}=-1 根轨迹方程 Ki=1n(spi)j=1m(szj)=1
模值条件
K ∗ = ∏ i = 1 n ∣ s − p i ∣ ∏ j = 1 m ∣ s − z j ∣ K^*=\frac{\prod_{i=1}^{n}|s-p_i|}{\prod_{j=1}^{m}|s-z_j|} K=j=1mszji=1nspi
相角条件:180度相角差


性质

  • 由于闭环特征方程中的某些系数是根轨迹增益 K ∗ K^* K 的函数,所以当 K ∗ K^* K 从零到无穷大连续变化时,特征方程的某些系数也随之而连续变化,因而特征方程式根的变化也必然是连续的,故根轨迹具有连续性。
  • 根轨迹必对称于实轴的原因是显然的,因为闭环特征方程式的根只有实根和复根两种,实根位于实轴上,复根必共轭,而根轨迹是根的集合,因此根轨迹对称于实轴
  • 根据对称性,只需做出上半 s 平面的根轨迹部分,然后利对称关系就可以画出下半 s 平面的根轨迹部分。

分离点

  • 一般情况下,常见的根轨迹分离点是位于实轴上的两条根轨迹分支的分离点。
  • 如果根轨迹位于实轴上两个相邻的开环极点之间,其中一个可以是无限极点,则在这两个极点之间至少存在一个分离点。
  • 同样,如果根轨迹位于实轴上两个相邻的开环零点之间,其中一个可以是无限零点,则在这两个零点之间也至少有一个分离点。

180度根轨迹绘制法则

在这里插入图片描述
在这里插入图片描述

对于特定 K ∗ K^* K 值下的闭环极点,可用模值条件确定。

在控制系统中,除根轨迹增益 K ∗ K^* K 以外,其他情形下的根轨迹统称为广义根轨迹。包括:

  • 参数根轨迹
  • 零度根轨迹
  • 非最小相位系统根轨迹

通常,将负反馈系统中 K ∗ K^* K 由零变到无穷大时的根轨迹叫做常规根轨迹(或180度根轨迹)。

以非开环增益/根轨迹增益为可变参数绘制的根轨迹称为参数根轨迹,以区别于以开环增益 K K K /根轨迹增益 K ∗ K^* K 为可变参数的常规根轨迹。

一般说来,零度根轨迹的来源有两个方面

  • 是非最小相位系统中包含 s 最高次幂的系数为负的因子。
  • 是控制系统中包含有正反馈内回路。

零度根轨迹绘制法则

在这里插入图片描述

如果系统的所有开环极点和零点都位于s平面的左半平面,则称为最小相位系统。

若系统的开环传递函数在右半 s 平面有零点或极点,则该系统称为非最小相位系统

  • 正反馈系统
  • 具有正反馈性质的系统
  • 具有右半s平面极零点

闭环零点的影响

闭环零点对系统动态性能的影响,相当于减小闭环系统的阻尼,从而使系统的过渡过程有出现超调的趋势,并且这种作用将随闭环零点接近坐标原点的程度而加强。

在开环系统中增加极点,会使根轨迹向右移动,从而降低系统的相对稳定性,增加系统响应的调节时间。

如果闭环零、极点相距很近,那么这样的闭环零、极点常称为偶极子。偶极子有实数偶极子和复数偶极子之分,而复数偶极子必共轭出现。不难看出,只要偶极子不十分接近坐标原点,它们对系统动态性能的影响就甚微,从而可以忽略它们的存在。

闭环零点的存在,将使系统的峰值时间提前,这相当于减小闭环系统的阻尼,从而使超调量加大,当闭环零点接近坐标原点时,这种作用尤甚。

一般说来,闭环零点对调节时间的影响是不定的。

闭环实数主导极点对系统性能的影响是:闭环实数主导极点的作用,相当于增大系统的阻尼,使峰值时间迟后,超调量下降

如果实数极点比共轭复数极点更接近坐标原点,甚至可以使振荡过程变为非振荡过程。

线性系统的频域分析法

频域特性

控制系统中的信号可以表示为不同频率正弦信号的合成。控制系统频率特性反映正弦信号作用下系统响应的性能。应用频率特性研究线性系统的经典方法称为频域分析法。频域分析法具有以下特点:

  • 控制系统及其元部件的频率特性可以运用分析法和实验方法获得,并可用多种形式的曲线表示,因而系统分析和控制器设计可以应用图解法进行
  • 频率特性物理意义明确。对于一阶系统和二阶系统,频域性能指标和时域性能指标有确定的对应关系;对于高阶系统,可建立近似的对应关系
  • 控制系统的频域设计可兼顾动态响应和噪声抑制两方面的要求
  • 频域分析法不仅适用于线性定常系统,还可以推广应用于某些非线性控制系统

在正弦信号作用下,输出稳态分量的幅值和相位的变化,称为幅值比 A ( w ) A(w) A(w) 和相位差 φ ( w ) \varphi(w) φ(w) , 且皆为输入正弦信号频率 w w w 的函数
{ A ( w ) = ∣ G ( j w ) ∣ φ ( w ) = ∠ G ( j w ) \left\{\begin{matrix} A(w)=\left | {G(jw)} \right | \\ \varphi(w)=\angle{G(jw)} \end{matrix}\right. {A(w)=G(jw)φ(w)=G(jw)
表明,对于稳定的线性定常系统,由谐波输入产生的输出稳态分量仍然是与输入同频率的谐波函数,而幅值和相位的变化是频率 w w w 的函数,且与系统数学模型相关。
为此,定义谐波输入下,输出响应中与输入同频率的谐波分量与谐波输入的幅值之比 A ( w ) A(w) A(w) 为幅频特性,相位之差 φ ( w ) \varphi(w) φ(w) 为相频特性,并称其指数表达形式:
G ( j w ) = A ( w ) e j φ ( w ) G(jw)=A(w)e^{j\varphi(w)} G(jw)=A(w)ejφ(w)
系统的频率特性

频率特性的定义既可以适用于稳定系统,也可适用于不稳定系统。稳定系统的频率特性可以用实验方法确定,即在系统的输入端施加不同频率的正弦信号,然后测量系统输出的稳态响应,再根据幅值比和相位差作出系统的频率特性曲线。频率特性也是系统数学模型的一种表达形式

对于不稳定系统,输出响应稳态分量中含有由系统传递函数的不稳定极点产生的呈发散或振荡的分量,所以不稳定系统的频率特性不能通过实验方法确定

频率特性定义(物理意义)

稳定系统的频率特性为零初始条件下,输出和输入的傅氏变换之比
G ( j w ) = C ( j w ) R ( j w ) = G ( s ) ∣ s = j w G(jw)=\frac{C(jw)}{R(jw)}=G(s)|_{s=jw} G(jw)=R(jw)C(jw)=G(s)s=jw
在工程分析和设计中,通常把线性系统的频率特性画成曲线,再运用图解法进行研究。常用的频率特性曲线有以下三种:

  • 幅相频率特性曲线

    • 又称为幅相曲线或极坐标图或奈奎斯特曲线(Nyquist Plot)

    • 以横轴为实轴、纵轴为虚轴,构成复数平面。对于任一给定的频率 w w w频率特性值为复数。若将频率特性表示为实数和虚数的形式,则实部为实轴坐标值,虚部为虚轴坐标值

    • 若将频率特性表示为复指数形式,则为复平面上的向量,而向量的长度为频率特性的幅值,向量与实轴正方向的夹角等于频率特性的相位

    • 由于幅频特性为 w w w 的偶数,相频特性为 w w w 的奇函数 w w w 从零变化至 + ∞ +\infty + w w w 从零变化至 − ∞ -\infty 的幅相曲线关于实轴对称,因此一般只绘制 w w w 从零变化至 + ∞ +\infty + 的幅相曲线

    • 在系统幅相曲线中,频率 w w w 为参变量一般用小箭头表示 w w w 增大时幅相曲线的变化方向

  • 对数频率特性曲线

    • 又称为伯德曲线或伯德图(Bode Diagram)

    • 对数频率特性曲线由对数幅频曲线和对数相频曲线组成

    • 对数频率特性曲线的横坐标按 lg ⁡ w \lg{w} lgw 分度,单位为弧度/秒(rad/s)

    • 对数幅频曲线的纵坐标按:
      对数幅频特性: L ( w ) = 20 lg ⁡ ∣ G ( j w ) ∣ = 20 lg ⁡ A ( w ) \text{对数幅频特性:}L(w)=20\lg{|G(jw)|}=20\lg{A(w)} 对数幅频特性:L(w)=20lgG(jw)=20lgA(w)
      线性分度,单位是分贝(dB)

    • 对数相频曲线的纵坐标按 φ ( w ) \varphi(w) φ(w) 线性分度,单位为度 ( ∘ ^\circ )。由此构成的坐标系称为半对数坐标系

  • 对数幅相曲线

    • 又称尼科尔斯曲线或尼科尔斯图(Nichols Plot)

    • 其特点是纵坐标为 L ( w ) L(w) L(w) ,单位为分贝 (dB);横坐标为 φ ( w ) \varphi(w) φ(w) ,单位为度 ( ∘ ^\circ ) ,均为线性分度,频率 w w w 为参变量

    • 在尼科尔斯曲线对应的坐标系中,可以根据系统开环和闭环的关系,绘制关于闭环幅频特性的等 M 簇线和闭环相频特性的等 α \alpha α 簇线,因而根据频域指标要求确定校正网络,简化系统的设计过程

典型环节与开环系统的频率特性

由于开环传递函数的分子和分母多项式的系数皆为实数,因此系统开环零极点或为实数或为共轭复数

根据开环零极点可将分子和分母多项式分解成因式,再将因式分类,即得典型环节

典型环节可分为两类:一类为最小相位环节;另一类为非最小相位环节

最小相位环节有下列七种:

  • 比例环节 K (K >0)

    • Nyquist 曲线: w 由 0 → ∞ ,始终为一点 k w\text{由}0\to\infty\text{,始终为一点 k} w0,始终为一点 k

    在这里插入图片描述

    • Bode 图

    在这里插入图片描述

    • 对数幅频特性曲线

      平行于频率轴的一条直线,值为 20 lg ⁡ k 20\lg{k} 20lgk

    • 对数相频特性曲线

      为零度直线

  • 惯性环节 1/(Ts+1) (T >0)

    • Nyquist 曲线

    在这里插入图片描述

    • Bode 图

      在这里插入图片描述
      在这里插入图片描述

  • 一阶微分环节 Ts + 1 (T >0)

  • 振荡环节 1 / ( s 2 / w n 2 + 2 ζ s / w n + 1 ) 1/(s^2/{w_n}^2+2\zeta{s}/w_n+1) 1/(s2/wn2+2ζs/wn+1) ( w n > 0 , 0 ≤ ζ < 1 ) (w_n>0,0\le\zeta<1) (wn>0,0ζ<1)

    • 谐振频率: w r = w n 1 − 2 ζ 2 , 0 < ζ ≤ 2 2 w_r=w_n\sqrt{1-2{\zeta}^2},0<\zeta\le\frac{\sqrt{2}}{2} wr=wn12ζ2 ,0<ζ22

    • 谐振峰值: M r = 1 2 ζ 1 − ζ 2 , 0 < ζ ≤ 2 2 M_r=\frac{1}{2\zeta\sqrt{1-\zeta^2}},0<\zeta\le\frac{\sqrt{2}}{2} Mr=2ζ1ζ2 1,0<ζ22

    • 由于 L a ( w ) L_a(w) La(w) 与阻尼比 ζ \zeta ζ 无关,用渐近线近似表示对数幅频曲线存在误差,误差的大小不仅和 w w w 有关,而且也和 ζ \zeta ζ 有关,误差曲线为一曲线簇,根据误差曲线可以修正渐近特性曲线而获得准确曲线。
      根据对数幅频特性定义还可知:
      (1)非最小相位振荡环节与振荡环节的对数幅频渐近特性曲线相同
      (2)二阶微分环节和非最小相位二阶微分环节与振荡环节的对数幅频渐近特性曲线关于0dB线对称

    • Nyquist 曲线

    在这里插入图片描述

    • Bode 图

    在这里插入图片描述

  • 二阶微分环节 s 2 / w n 2 + 2 ζ s / w n + 1 s^2/{w_n}^2+2\zeta{s}/w_n+1 s2/wn2+2ζs/wn+1 ( w n > 0 , 0 ≤ ζ < 1 ) (w_n>0,0\le\zeta<1) (wn>0,0ζ<1)

  • 积分环节 1/s

  • 微分环节 s

  • 延迟环节与延迟系统

    • 时域表达式 c ( t ) = 1 ( 1 − τ ) r ( t − τ ) c(t)=1(1-\tau)r(t-\tau) c(t)=1(1τ)r(tτ)
    • 传递函数 G ( s ) = e − τ s G(s)=e^{-\tau{s}} G(s)=eτs
    • 频率特性 G ( j w ) = e − j τ w G(jw)=e^{-j\tau{w}} G(jw)=ejτw

非最小相位环节共有五种:

  • 比例环节 K (K<0)
  • 惯性环节 1/(-Ts+1) (T >0)
  • 一阶微分环节 -Ts+1 (T >0)
  • 振荡环节 1 / ( s 2 / w n 2 − 2 ζ s / w n + 1 ) 1/(s^2/{w_n}^2-2\zeta{s}/w_n+1) 1/(s2/wn22ζs/wn+1) ( w n > 0 , 0 < ζ < 1 ) (w_n>0,0<\zeta<1) (wn>0,0<ζ<1)
  • 二阶微分环节 s 2 / w n 2 − 2 ζ s / w n + 1 s^2/{w_n}^2-2\zeta{s}/w_n+1 s2/wn22ζs/wn+1 ( w n > 0 , 0 < ζ < 1 ) (w_n>0,0<\zeta<1) (wn>0,0<ζ<1)

除了比例环节外,非最小相位环节和与之相对应的最小相位环节的区别在于开环零极点的位置

开环频率特性曲线的绘制方法

已知系统开环传递函数 G k ( s ) G_k(s) Gk(s) ,绘制开环辐相曲线

  • 由已知的 G k ( s ) G_k(s) Gk(s) G k ( j w ) = G k ( s ) ∣ s = j w , A ( s ) , ϕ ( w ) , P ( w ) , Q ( w ) G_k(jw)=G_k(s)|_{s=jw},A(s),\phi(w),P(w),Q(w) Gk(jw)=Gk(s)s=jw,A(s),ϕ(w),P(w),Q(w)
  • 分别对: w : 0 + → + ∞ , − ∞ → 0 − , 0 − → 0 + w:0^+\to+\infty,-\infty\to0^-,0^-\to0^+ w:0++,0,00+
    • 求N氏曲线的起点
    • 求N氏曲线的终点
    • 求一些特色点
  • w : 0 + → + ∞ , − ∞ → 0 − w:0^+\to+\infty,-\infty\to0^- w:0++,0曲线关于实轴镜像对称

频率域稳定判据

F ( s ) = ( s − z 1 ) ( s − z 2 ) … ( s − z m ) ( s − p 1 ) ( s − p 2 ) … ( s − p n ) = ∑ j = 1 m ( s − z j ) ∑ i = 1 n ( s − p i ) F(s)=\frac{(s-z_1)(s-z_2)\dots(s-z_m)}{(s-p_1)(s-p_2)\dots(s-p_n)}=\frac{\sum_{j=1}^{m}(s-z_j)}{\sum_{i=1}^{n}(s-p_i)} F(s)=(sp1)(sp2)(spn)(sz1)(sz2)(szm)=i=1n(spi)j=1m(szj)

如果 F(s) 是非奇异的(除有限个孤立奇点 p 1 p_1 p1 p 2 p_2 p2 p n p_n pn 外处处解析),则有:

  • 对于 s 平面上任意一点 s,通过复变函数 F(s) 的映射关系,在 F(s)平面上都可以确定一个关于 s 的象或映射 S F S_F SF
  • 对于 s 平面上任意一条不通过 F(s) 的任一零点和极点的封闭曲线 Γ \Gamma Γ ,都可以在F(s)平面上确定一条相应的封闭曲线 Γ F \Gamma_F ΓF Γ \Gamma Γ 的象或映射)

辐角原理

s 平面上不通过 F(s) 的任一零点和极点的封闭曲线 Γ \Gamma Γ 包围 F(s) 的P个极点和Z个零点,则 s 从封闭曲线 Γ \Gamma Γ上任一点 A 起,顺时针沿 Γ \Gamma Γ 移动一周,再回到 A 点时,

则相应地, 在F(s) 平面上,封闭曲线 Γ \Gamma Γ 的象或映射亦从点 F(A) 起,逆时针方向绕原点 R = P − Z R=P-Z R=PZ,回到 F(A) 点止,形成一条闭合曲线 Γ F \Gamma_F ΓF

  • R <0 和 R > 0 分别表示 Γ F \Gamma_F ΓF 顺时针包围和逆时针包围F(s) 平面的原点
  • R = 0表示不包围 F(s) 平面的原点

封闭曲线 Γ \Gamma Γ 及其象或映射 Γ F \Gamma_F ΓF 的形状对于R值的确定没有影响,也即,R值只与封闭曲线 Γ \Gamma Γ 包围 F(s) 零点、极点的个数有关,而与封闭曲线 Γ \Gamma Γ 的形状无关

F(s)具有以下特点:

  • F(s) 的零点为闭环特征方程的根,即闭环传递函数的极点
  • F(s) 的极点为开环传递函数的极点
  • F(s) 的零点数与极点数相同(均为n)(考虑开环传递函数分母多项式的阶次一般大于或等于分子多项式的阶次)
  • s 沿封闭曲线 Γ \Gamma Γ 在 s 平面运动一周所产生的象或映射:封闭曲线 Γ F \Gamma_F ΓF (闭环)和 Γ G H \Gamma_{GH} ΓGH (开环),在位置上只相差常数 1,即闭合曲线 Γ F \Gamma_F ΓF 可由 Γ G H \Gamma_{GH} ΓGH 沿实轴正方向平移一个单位长度获得

系统的闭环稳定性取决于系统闭环传递函数极点即F(s)的零点的位置

考虑选取 s 平面上可以围 住其整个右半平面的封闭曲线 Γ \Gamma Γ (称为N氏围线),若N氏围线 Γ \Gamma Γ 没有围住F(s)的零点,则闭环系统稳定;若围住了1个零点,系统就不稳定

设 Z 为 F(s) 零点在右半 s 平面中的个数,则 Z =0 表示F(s)在右半s平面无零点,闭环特征方程 F(s)=0 的根全部位于s平面的左半平面,系统稳定;Z≠0,系统不稳

设P为F(s)极点在右半s平面中的个数, 即为系统开环极点在右半s平面的个数;若已知 G k ( s ) G_k(s) Gk(s),则P已知

原始的Nyquist稳定判据步骤:

  1. 已知: G k ( s ) G_k(s) Gk(s)
    则P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)

  2. 在F(s) 平面上画出 F ( s ) = 1 + G k ( s ) F(s)=1+G_k(s) F(s)=1+Gk(s) 曲线
    则可知R(当s在s 平面N氏围线 Γ \Gamma Γ 上顺时针方向移动一周时,F(s)平面上 F ( s ) = 1 + G k ( s ) F(s)=1+G_k(s) F(s)=1+Gk(s) 曲线逆时针绕原点的圈数)

  3. 在自变量s平面(N氏围线 Γ \Gamma Γ) →因变量F(s) 平面(象或映射 Γ F \Gamma_F ΓF )上直接应用幅角原理,由R=P-Z,可得:Z=P-R

    若 Z =0 ,则系统稳定
    若 Z≠0,则系统不稳定

F(s)平面上的原点,就是G(s)H(s)平面上的(-1,j0)点

封闭曲线 Γ F \Gamma_F ΓF 包围 F(s) 平面原点的圈数,就等于封闭曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 包围 F(s) 平面点 (-1, j0) 的圈数

0型系统:当自变量 s 从 s 平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s) H(s)平面上曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 映射到 G(jω)H(jω) 平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) ,即为当
ω:-∞→+∞的Nyquist曲线(全闭合曲线)

I型或I型以上系统:当自变量s从s平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s)H(s) 平面上曲线 Γ G H \Gamma_{GH} ΓGH 映射到 G(jω)H(jω) 平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) ,即为当 w: 0+ →+∞、 -∞→ 0-、0- → 0+的Nyquist曲线(全闭合曲线)

具有等幅振荡环节的系统:当自变量 s 从 s 平面的原点出发,顺时针沿N氏围线 Γ \Gamma Γ 移动一周时,G(s)H(s) 平面上曲线 Γ G ( s ) H ( s ) \Gamma_{G(s)H(s)} ΓG(s)H(s) 映射到G(jω)H(jω)平面上的曲线 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw),即为当 w: 0 → w n − w_{n-} wn w n − w_{n-} wn w n + w_{n+} wn+ w n + w_{n+} wn+→+∞、 -∞→ w n + w_{n+} wn+
− w n + -w_{n+} wn+ − w n − -w_{n-} wn − w n − -w_{n-} wn → 0 的 Nyquist 曲线(全闭合曲线)

幅角原理(推广)

设 s 平面上不通过 F(s) 的任一零点和极点的N氏围线 Γ \Gamma Γ (围住整个右半s平面)包围 F(s) 平面的 P 个极点和 Z 个零点,则自变量 s 从原点顺时针沿 Γ \Gamma Γ 移动一周时, 在 G(jω)H(jω) 平面上, N 氏围线 Γ \Gamma Γ 的象或映射 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(全闭合 Nyquist 曲线w: 0+ →+∞、 -∞→ 0-、0- → 0+) ,将以逆时针方向绕 (-1,j0) 点 R=P-Z 圈

实用的Nyquist稳定判据(全闭合Nyquist曲线)步骤:

  1. 已知: G k ( s ) G_k(s) Gk(s)
    则v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)

  2. 画出w: 0 →+∞、 -∞→ 0 开环(全闭合)Nyquist曲线

    • 0型系统:画w: 0 →+∞、 -∞→ 0 Nyquist曲线
    • I型及以上系统: 画w: 0 + 0^+ 0+→+∞的Nyquist曲线
      补画w: -∞→ 0 − 0^- 0、0- → 0+的Nyquist曲线
    • 具有 s = ± j w n s=±jw_n s=±jwn系统: 画0 → w n − w_{n-} wn w n − w_{n-} wn w n + w_{n+} wn+ w n + w_{n+} wn+ →+∞、 -∞→ - w n + w_{n+} wn+、 - w n + w_{n+} wn+ →- w n − w_{n-} wn − w n − -w_{n-} wn→ 0 的Nyquist曲线

    则可根据画好的全闭合Nyquist曲线,从图上获知Nyquist曲线逆时针绕(-1,j0)点的圈数,即可得R(当s在s 平面N氏围线 Γ \Gamma Γ 上顺时针方向移动一周时)

  3. 在s平面上N氏围线 Γ \Gamma Γ →G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(全闭合Nyquist曲线),应用幅角原理,由R=P-Z(+逆,-顺)可得:Z=P-R

若Z =0 ,则系统稳定
若Z≠0,则系统不稳定

半闭合Nyquist曲线

设s平面上不通过F(s) 的任一零点和极点的N氏围线 Γ \Gamma Γ (围住整个右半s平面)包围F(s)平面的P个极点和Z个零点,则自变量s 从原点顺时针沿 Γ \Gamma Γ 移动一周时, 在 G(jω)H(jω) 平面上,N氏围线 Γ \Gamma Γ 的象或映射 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合Nyquist曲线 w: 0→ + ∞ +\infty +),将以逆时针方向绕(-1,j0) 点 R=P-Z 圈( R = 2 ( N + − N − ) = 2 N R=2(N^+-N^-)=2N R=2(N+N)=2N

扩展的Nyquist稳定判据(半闭合Nyquist曲线)步骤:(常用)

  1. 已知: G k ( s ) G_k(s) Gk(s)
    则 v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)

  2. 画出w: 0 → + ∞ +\infty +开环(半闭合)Nyquist曲线;

    • 0型系统:画w: 0 → + ∞ +\infty + Nyquist曲线;

    • I型及以上系统

      • 画w: 0 + → + ∞ 0^+→+\infty 0++Nyquist曲线
      • 从G(j0+)H(j0+)点补画(逆时针方向补)半径为 ∞ \infty 角度变化 v π 2 \frac{v\pi}{2} 2vπ 的虚圆弧;(全闭合曲线为半径为∞,角度从0- vπ/2 → 0+ -vπ/2的圆弧)
    • 具有 s = ± j w n s=±jw_n s=±jwn系统

      • 0 → w n − 、 w n + → + ∞ 0 →w_{n-} 、w_{n+} →+\infty 0wnwn++ Nyquist曲线
      • G ( j w n − ) H ( j w n − ) G(jw_{n-})H(jw_{n-}) G(jwn)H(jwn) 点起补画(顺时针方向补)半径为 ∞ \infty 角度变化 v × π v\times\pi v×π 的虚圆弧至 G ( j w n + ) H ( j w n + ) G(jw_{n+})H(jw_{n+}) G(jwn+)H(jwn+)

    则可根据画好的半闭合 Nyquist 曲线获知,位于(-1,j0)点左侧,从上向下穿越(II → III)次数和 N+(逆),以及从下向上穿越(III → II)次数和 N-(顺),即可得 R = 2 ( N + − N − ) = 2 N R=2(N^+-N^-)=2N R=2(N+N)=2N(当s在s 平面N氏围线 Γ \Gamma Γ上顺时针方向移动一周时)

  3. 在 s 平面上N氏围线 Γ \Gamma Γ →G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合Nyquist曲线),应用幅角原理,由R=P-Z,可得: Z=P-R

    若Z =0 ,则系统稳定
    若Z≠0,则系统不稳定

对数频率稳定判据(Bode图上的Nyquist稳定判据)

奈氏判据基于复平面的半闭合曲线 Γ G H \Gamma_{GH} ΓGH 判定系统的闭环稳定性,由于半闭合曲线 Γ G H \Gamma_{GH} ΓGH 可以转换为半对数坐标下的曲线,因此可以推广运用奈氏判据,其关键问题是需要根据半对数坐标下的 Γ G H \Gamma_{GH} ΓGH 曲线确定穿越次数 N 或 N+和N-

复平面 Γ G H \Gamma_{GH} ΓGH 曲线一般由两部分组成:开环幅相曲线和开环系统存在积分环节和等幅振荡环节时所补作的半径为无穷大的虚圆弧。而N 的确定取决于 A(w)>1 时 Γ G H \Gamma_{GH} ΓGH 穿越负实轴的次数,因此应建立和明确以下对应关系:

穿越点确定
{ A ( w c ) = ∣ G ( j w c ) H ( j w c ) ∣ = 1 L ( w c ) = 20 lg ⁡ A ( w C ) = 0 \left\{\begin{matrix} A(w_c)=|G(jw_c)H(jw_c)|=1 \\ L(w_c)=20\lg A(w_C)=0 \end{matrix}\right. {A(wc)=G(jwc)H(jwc)=1L(wc)=20lgA(wC)=0

w c w_c wc 为截止频率

对于复平面的负实轴和开环对数相频特性,当取频率为穿越频率 w x w_x wx 时:
φ ( w x ) = ( 2 k + 1 ) π , k = 0 , ± 1 , … \varphi(w_x)=(2k+1)\pi,k=0,\pm1,\dots φ(wx)=(2k+1)π,k=0,±1,
设半对数坐标下 Γ G H \Gamma_{GH} ΓGH 的对数幅频曲线和对数相频曲线分别为 Γ L \Gamma_{L} ΓL Γ φ \Gamma_{\varphi} Γφ 。由于 Γ L \Gamma_{L} ΓL 等于L(w)曲线,则 Γ G H \Gamma_{GH} ΓGH 在A(w)>1时,穿越负实轴的点等于 Γ G H \Gamma_{GH} ΓGH 在半对数坐标下,对数幅频特性L(w)>0时对数相频特性曲线 Γ φ \Gamma_{\varphi} Γφ ( 2 k + 1 ) π ; k = 0 , ± 1 , … (2k+1)\pi; k=0, \pm1, \dots (2k+1)π;k=0,±1,, 平行线的交点

Γ φ \Gamma_{\varphi} Γφ 确定

  • 开环系统无虚轴上极点时, Γ φ \Gamma_{\varphi} Γφ 等于 φ ( w ) \varphi(w) φ(w) 曲线

  • 开环系统存在积分环节 1 s v , ( v > 0 ) \frac{1}{s^v},(v>0) sv1,(v>0),复数平面的 Γ G H \Gamma_{GH} ΓGH 曲线,需从w=0+的开环幅相曲线的对应点G( j0+) H( j0+) 起,逆时针补作 v × 9 0 ∘ v\times90^{\circ} v×90 半径为无穷大的虚圆弧。对应地,需从对数相频特性曲线 w 较小且 L(w)>0 的点处向上补作 v × 9 0 ∘ v\times90^{\circ} v×90 的虚直线, φ ( w ) \varphi(w) φ(w) 曲线和补作的虚直线构成 Γ φ \Gamma_{\varphi} Γφ

  • 开环系统存在等幅振荡环节: 1 ( s 2 + w n 2 ) v 1 , v 1 > 0 \frac{1}{{(s^2+w_n^2)}^{v_1}},v_1>0 (s2+wn2)v11,v1>0
    时,复数平面的 Γ G H \Gamma_{GH} ΓGH 曲线,需从 w = w n − w = w_{n-} w=wn的开环幅相曲线的对应点 G ( j w n − ) H ( j w n − ) G( jw_{n-})H(jw_{n-}) G(jwn)H(jwn) 起,逆时针补作 v 1 × 18 0 ∘ v_1\times180^{\circ} v1×180 半径为无穷大的虚圆弧至 w = w n + w = w_{n+} w=wn+的对应点 G ( j w n + ) H ( j w n + ) G( jw_{n+})H(jw_{n+}) G(jwn+)H(jwn+)。对应地,需从对数相频特性曲线 φ ( w n − ) \varphi(w_{n-}) φ(wn)点起向上补作 v 1 × 18 0 ∘ v_1\times180^{\circ} v1×180 的虚直线至 φ ( w n + ) \varphi(w_{n+}) φ(wn+) 处, φ ( w ) \varphi(w) φ(w)曲线和补作的虚直线构成 Γ φ \Gamma_{\varphi} Γφ

穿越次数计算

  • 正穿越一次: Γ G H \Gamma_{GH} ΓGH 由上向下穿越 (-1, j0) 点左侧的负实轴一次,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由下向上穿越 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次。

  • 负穿越一次: Γ G H \Gamma_{GH} ΓGH 由下向上穿越 (-1, j0) 点左侧的负实轴一次,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由上向下穿越 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次

  • 正穿越半次: Γ G H \Gamma_{GH} ΓGH 由上向下止于或由上向下起于 (-1, j0) 点左侧的负实轴,等价于:在L(w)>0时, Γ φ \Gamma_{\varphi} Γφ 由下向上止于或由下向上起于 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次

  • 负穿越半次: Γ G H \Gamma_{GH} ΓGH 由下向上止于或由下向上起于(-1, j0)点左侧的负实轴,等价于:在 L(w)>0 时, Γ φ \Gamma_{\varphi} Γφ 由上向下止于或由上向下起于 ( 2 k + 1 ) π (2k+1)\pi (2k+1)π 线一次

应该指出的是,补作的虚直线所产生的穿越皆为负穿越

Bode图(半对数坐标 Γ G H \Gamma_{GH} ΓGH曲线)上的Nyquist稳定判据步骤:

  1. 已知: G k ( s ) G_k(s) Gk(s)
    则v、P已知(P为 G k ( s ) G_k(s) Gk(s) 在右半s平面上开环极点的个数)

  2. 画出 w : 0 → + ∞ w: 0 →+\infty w:0+开环半对数坐标 Γ G H \Gamma_{GH} ΓGH 曲线(Bode图)

    • 0型系统:画 w : 0 + → + ∞ w: 0^+ →+\infty w:0++ 对数相频特性曲线 φ ( w ) \varphi(w) φ(w)
    • I型及以上系统
      • w : 0 + → + ∞ w: 0^+ →+\infty w:0++ 对数相频特性曲线 φ ( w ) \varphi(w) φ(w)
      • 从对数相频特性曲线 w 较小且 L(w)>0 的点处向上补作 v × 9 0 ∘ v\times90^{\circ} v×90的虚直线, φ ( w ) \varphi(w) φ(w)曲线和补作的虚直线构成 Γ φ \Gamma{\varphi} Γφ
    • 具有 s = ± j w n s=±jw_n s=±jwn系统
      • w : 0 + → + ∞ w: 0+ →+\infty w:0++对数相频特性曲线 φ ( w ) \varphi(w) φ(w)
      • 从对数相频特性曲线 φ ( w n − ) \varphi(w_{n-}) φ(wn) 点起向上补作 v × 18 0 ∘ v\times180^{\circ} v×180 的虚直线至 φ ( w n + ) \varphi(w_{n+}) φ(wn+)处, φ ( w ) \varphi(w) φ(w)曲线和补作的虚直线构成 Γ φ \Gamma{\varphi} Γφ

    从画好的对数幅频特性L(w)>0时,对数相频特性曲线 Γ φ \Gamma{\varphi} Γφ ( 2 k + 1 ) π ; k = 0 , ± 1 , … (2k+1)\pi; k=0, \pm1, \dots (2k+1)π;k=0,±1,, 平行线的交点,可得R=2(N±N-)=2N

  3. 在 s 平面上N氏围线 Γ \Gamma Γ → G(jω)H(jω)平面上的 Γ G ( j w ) H ( j w ) \Gamma_{G(jw)H(jw)} ΓG(jw)H(jw) 曲线(半闭合 Nyquist 曲线),应用幅角原理,由R=P-Z,可得: Z=P-R
    若 Z = 0 ,则系统稳定;若 Z ≠ 0,则系统不稳定

稳定裕度

当开环传递函数的某些系数发生变化时, Γ G H \Gamma_{GH} ΓGH 包围 (-1, j0) 点的情况亦随之改变

在稳定性研究中,称 (-1, j0)点为临界点,而闭合曲线 Γ G H \Gamma_{GH} ΓGH 相对于临界点的位置即偏离临界点的程度,反映系统的相对稳定性

频域的相对稳定性即稳定裕度常用相角裕度 γ \gamma γ 和幅值裕度 h h h 来度量

相角裕度

w c w_c wc 为系统的截止频率(固定模为1,比较相角离 (-1, j0)点的远近):
A ( w c ) = ∣ G ( j w c ) H ( j w c ) ∣ = 1 A(w_c)=|G(jw_c)H(jw_c)|=1 A(wc)=G(jwc)H(jwc)=1
定义相角裕度为:
γ = 18 0 ∘ + ∠ G ( j w c ) H ( j w c ) \gamma=180^{\circ}+\angle{G(jw_c)H(jw_c)} γ=180+G(jwc)H(jwc)
相角裕度 γ \gamma γ 的含义是,对于闭环稳定系统,如果系统开环相频特性再滞后 γ \gamma γ 度,则系统将处于临界稳定状态

幅值裕度

w x w_x wx 为系统的穿越频率(固定相角为-π,比较模离 (-1, j0)点的远近) :
φ ( w x ) = ∠ G ( j w x ) H ( j w x ) = ( 2 k + 1 ) π ; k = 0 , ± 1 , … \varphi(w_x)=\angle{G(jw_x)H(jw_x)}=(2k+1)\pi;k=0,\pm1,\dots φ(wx)=G(jwx)H(jwx)=(2k+1)π;k=0,±1,
定义幅值裕度为:
h = 1 ∣ G ( j w x ) H ( j w x ) ∣ h=\frac{1}{|G(jw_x)H(jw_x)|} h=G(jwx)H(jwx)1
幅值裕度 h h h 的含义是,对于闭环稳定系统,如果系统开环幅频特性再增大 h 倍,则系统将处于临界稳定状态

对数坐标下,幅值裕度按下式定义:
h = − 20 lg ⁡ ∣ G ( j w x ) H ( j w x ) ∣ ( d B ) h=-20\lg{|G(jw_x)H(jw_x)|}(dB) h=20lgG(jwx)H(jwx)(dB)
仅用相角裕度或幅值裕度,都不足以反映系统的稳定程度。 对非最小相位系统并不都适用

用频域法分析系统性能

开环频域性能指标

稳:相角裕度、幅值裕度
快:幅值穿越频率(截止频率、剪切频率)、相角穿越频率
准:稳态误差

在开环频率特性曲线上Nyquist曲线、Bode图

闭环频域性能指标

稳:谐振峰值
快:带宽频率、谐振频率
准:零频幅值

在闭环频率特性曲线上Nyquist曲线、Bode图、闭环幅频曲线、等M/N圆、尼科尔斯图

  • 零频幅值 M 0 M_0 M0
    w = 0 M ( 0 ) = M 0 = C ( 0 ) R ( 0 ) M ( 0 ) = { k 1 + k , v = 0 1 , v > 1 w=0\\ M(0)=M_0=\frac{C(0)}{R(0)}\\ M(0)=\begin{cases} \frac{k}{1+k},v=0 \\ 1,v>1 \end{cases} w=0M(0)=M0=R(0)C(0)M(0)={1+kk,v=01,v>1
    M ( 0 ) = 1 M(0)=1 M(0)=1 表示系统阶跃响应的终值等于输入,静差为零

    M ( 0 ) ≠ 1 M(0)\ne1 M(0)=1表示系统有静差

    M 0 M_0 M0 与1相差的大小,反映了系统的稳态精度, M 0 M_0 M0 越接近1,系统精度越高

  • 谐振频率 w r w_r wr

    闭环幅频特性出现峰值时的频率称为谐振频率 w r w_r wr,它在一定程度上反映了系统的快速性

    w r w_r wr 越大,系统瞬态响应越快

  • 谐振峰值 M r M_r Mr

    闭环幅频特性最大值 M r = M ( w r ) M_r=M(w_r) MrM(wr) 称为谐振峰值 M r M_r Mr

    峰值大,表明系统对某个频率的正弦输入信号反映强烈,有共振的倾向,即意味着系统的平稳性较差,阶跃响应将有过大的超调量

    一般要求: M r < 1.5 M 0 M_r<1.5M_0 Mr<1.5M0

  • 带宽频率 w b w_b wb 与系统带宽 0 ~ w b w_b wb

    w b w_b wb幅频特性 M ( w ) M(w) M(w) 的数值从 M ( 0 ) M(0) M(0) 衰减到0.707 M ( 0 ) M(0) M(0) 时所对应的频率

    w b w_b wb 高: M ( w ) M(w) M(w) M ( 0 ) M(0) M(0) 到 0.707 M ( 0 ) M(0) M(0) 所占据的频率区间0~ w b w_b wb 较宽,
    明系统复现快速变化的信号能力强,失真小,也即快速性好,阶跃响应的调节时间短

    w b w_b wb 低: M ( w ) M(w) M(w) M ( 0 ) M(0) M(0) 到 0.707 M ( 0 ) M(0) M(0) 所占据的频率区间0~ w b w_b wb 较窄,带 宽窄,系统反映迟钝,失真大,快速性差(复现快速信号能力差)

系统带宽与信号频谱的关系

( 0 , w b ) (0,w_b) (0,wb) 称为系统的带宽

带宽定义表明,对于高于带宽频率的正弦输入信号,系统输出将呈现较大的衰减

当控制输入为矩形波时,系统的跟踪能力取决于带宽覆盖矩形波频谱的范围,带宽大则处于较高频率范围的谱线衰减小,故失真小

当控制输入信号频谱特性具有收敛形式时,可按跟踪要求选定常数 ε \varepsilon ε ,若 w > w b w>w_b w>wb 时, ∣ F ( j w ) ∣ < ε |F(jw)|<\varepsilon F(jw)<ε,从而确定系统所需要的带宽 w b w_b wb

系统分析应区分输入信号的性质、位置,根据其频谱或谱密度、相应的传递函数及噪声频段选择带宽

将控制系统中的各个变量看成一些信号而这些信号又是由许多不同频率的正弦信号合成的,各个变量的运动就是系统对各个不同频率的信号的响应的综合

频域分析法进行系统分析的基本思路

  1. 建模

    建立系统(1阶、2阶、任意阶次)数学模型:微分方程/传递函数/频率特性

  2. 画系统频率特性曲线

    • 系统开环频率特性曲线(Nyquist曲线/波德图)

    • 系统闭环频率特性(Nyquist曲线/波德图/幅频曲线)

  3. 系统分析(在开/闭环频率特性曲线上)

    • 稳定性分析

    • 性能分析(求系统的性能指标)

    • 动态响应性能分析(稳、快)

    • 稳态响应性能分析

线性系统的校正方法

所谓校正,就是在系统中加入一些其参数可以根据需要而改变的机构或装置,使系统整个特性发生变化,从而满足给定的各项性能指标

在控制系统设计中,采用的设计方法一般依据性能指标的形式而定:

  • 如果性能指标以单位阶跃响应的峰值时间、调节时间、超调量、阻尼比、稳态误差等时域特征量给出,一般采用根轨迹法校正
  • 如果性能指标以系统的相角裕度、幅值裕度、谐振峰值、闭环带宽、静态误差系数等频域特征量给出,一般采用频率法校正

目前,工程技术界多习惯采用频率法,故通常通过近似公式进行两种指标的互换

时域性能指标频域开环性能指标频域闭环性能指标
σ % \sigma\% σ% γ \gamma γ M r M_r Mr
t s , t p , t r , t d t_s,t_p,t_r,t_d ts,tp,tr,td w c w_c wc w r , w b w_r,w_b wr,wb
K , v ( e s s ) K,v(e_{ss}) K,v(ess) K , v K,v K,v M ( 0 ) , v M(0),v M(0),v

无论采用哪种校正方式,都要求校正后的系统既能以所需精度跟踪输入信号,又能抑制噪声扰动信号

在控制系统实际运行中,输入信号一般是低频信号,而噪声信号则一般是高频信号

为了使系统具有较高的稳定裕度,希望系统开环对数幅频特性在截止频率 w c w_c wc 处的斜率为 -20dB/dec,但从要求系统具有较强的从噪声中辩识信号的能力来考虑, 却又希望 w c w_c wc 处的斜率小于 40dB/dec

通常, 一个设计良好的实际运行系统, 其相角裕度具有45º左右的数值:

  • 过低于此值,系统的动态性能较差,且对参数变化的适应能力较弱
  • 过高于此值,意味着对整个系统及其组成部件要求较高。因此造成实现上的困难,或因此不满足经济性要求,同时由于稳定程度过好,造成系统动态过程缓慢

要实现 45º 左右的相角裕度要求,开环对数幅频特性在中频区的斜率应为-20dB/dec,同时要求中频区占据一定的频率范围,以保证在系统参数变化时,相角裕度变化不大

过此中频区后,要求系统幅频特性迅速衰减,以削弱噪声对系统的影响。这是选择系统带宽应该考虑的一个方面。另一方面,进入系统输入端的信号,既有输入信号r(t),又有噪声信号 n(t),如果输入信号的带宽为 0 ∼ w M 0\sim{w_M} 0wM,噪声信号集中起作用的频带为 w 1 ∼ w n w_1\sim{w_n} w1wn则控制系统的带宽频率通常取为:
w b = ( 5 ∼ 10 ) w M w_b=(5\sim10)w_M wb=(510)wM
且使 w 1 ∼ w n w_1\sim{w_n} w1wn 处于 0 ∼ w b 0\sim{w_b} 0wb 范围之外

根据稳、快、准的要求,选择时注意:

  • 各元部件的非线性因素,如不是敏区、磨擦、间隙、滞环和漂移等应尽量小
  • 避免小阻尼的环节和非最小相位环节
  • 受控对象和各元部件之间的时间常数应适当拉开
  • 各级增益的配置要有利于提高系统的线性度和抗扰能力

按照校正装置在系统中的连接方式,控制系统校正方式可分为四种:

  • 串联校正装置一般接在系统误差测量点之后和放大器之前,串接于系统前向通道之中

  • 反馈校正装置一般接在系统局部反馈通路之中

  • 前馈校正又称顺馈校正,是在系统主反馈回路之外采用的校正方式

    前馈校正装置接在系统给定值(或指令、参考输入信号)之后及主反馈作用点之前的前向通道上。这种校正方式的作用相当于对给定值信号进行整形或滤波后,再送入反馈系统

    另一种前馈校正装置接在系统可测扰动作用点与误差测量点之间,对扰动信号进行直接或间接测量,并经变换后接入系统,形成一条附加的对扰动影响进行补偿的通道

    前馈校正可以单独作用于开环控制系统,也可以作为反馈控制系统的附加校正而组成复合控制系统

  • 复合校正方式是在反馈控制回路中,加入前馈校正通路,组成一个有机整体;
    复合校正装置与前馈校正一致,有 2 种接入方式:

    • 按扰动补偿的复合控制形式
    • 按输入补偿的复合控制形式

在控制系统设计中,常用的校正方式为串联校正和反馈校正两种

一般来说,串联校正设计比反馈校正设计简单,也比较容易对信号进行各种必要形式的变换

  • 在直流控制系统中,由于传递直流电压信号,适于采用串联校正
  • 在交流载波控制系统中,如果采用串联校正,一般应接在解调器和滤波器之后,否则由于参数变化和载频漂移,校正装置的工作稳定性很差

串联校正装置又分无源和有源两类

  • 无源串联校正装置通常由 Rc 无源网络构成,结构简单,成本低廉,但会使信号在变换过程中产生幅值衰减,且其输入阻抗较低,输出阻抗又较高,因此常需要附加放大器,以补偿其幅值衰减,并进行阻抗匹配。为了避免功率损耗,无源串联校正装置通常安置在前向通路中能量较低的部位上
  • 有源串联校正装置运算放大器和 Rc 网络组成,其参数可以根据需要调整,因此在工业自动化设备中,经常采用由电动(或气动)单元构成的 PID 控制器(或称 PID 调节器),它由比例单元、微分单元和积分单元组合而成,可以实现各种要求的控制规律

在实际控制系统中,还广泛采用反馈校正装置。一般来说,反馈校正所需元件数目比串联校正多。由于反馈信号通常由系统输出端或放大器输出级供给,信号是从高功率点传向低功率点,因此反馈校正一般无需附加放大器

此外,反馈校正尚可消除系统原有部分参数波动对系统性能的影响。在性能指标要求较高的控制系统设计中,常常兼用串联校正与反馈校正两种方式

PID 控制器(或称 PID 调节器)

确定校正装置的具体形式时,应先了解校正装置所需提供的控制规律,以便选择相应的元件

包含校正装置在内的控制器,常常采用比例P、微分D、积分I等基本控制规律,或者采用这些基本控制规律的某些组合,如比例微分PD、比例积分PI、比例-积分-微分PID等组合控制规律,以实现对被控对象的有效控制

比例 P 控制规律

  • 具有比例控制规律的控制器称为P 控制器, 其中 K p K_p Kp 称为 P 控制器增益
  • P 控制器实质上是一个具有可调增益的放大器
  • 在信号变换过程中,P 控制器只改变信号的增益而不影响其相位
  • 在串联校正中,加大控制器增益 K p K_p Kp ,可以提高系统的开环增益,减小系统稳态误差,从而提高系统的控制精度,但会降低系统的相对稳定性,甚至可能造成闭环系统不稳定
  • 因此,在系统校正设计中,很少单独使用比例控制规律

在这里插入图片描述

串联滞后-超前网络校正装置的特性(作用)

  • 使滞后部分发生在系统频率特性的低频段,以提高系统的稳态性能
  • 使超前部分发生在系统频率特性的中频段, 以改善系统的动态性能

在这里插入图片描述

线性离散系统的分析与校正

z 变换法是离散系统理论的数学基础。利用 z 变换法研究离散系统,可以把连续系统中的许多概念和方法,推广应用于线性离散系统

线性离散系统分析方法的基本思路

  • 建模

    建立线性离散系统数学模型:(前向或后向)差分方程、闭环脉冲传递函数

  • 求时域响应

    已知条件:

    • 系统数学模型:脉冲传递函数(或:差分方程→脉冲传递函数)

    • 输入(典型输入):单位阶跃信号(单位脉冲,单位斜坡,单位加速度)

      求: c ∗ ( t ) / c ( n T ) / c ( k ) c^*(t)/c(nT)/c(k) c(t)/c(nT)/c(k)

  • 系统分析

    • 稳定性分析

    • 性能分析(求系统的性能指标)

      • 动态响应性能分析 (稳、快) - 定性分析/近似计算
      • 稳态响应性能分析 (准)

离散系统的基本概念

  • 如果控制系统中的所有信号,都是时间变量的连续函数,则这样的系统称为连续时间系统,简称连续系统
  • 如果控制系统中有一处或几处信号,是一串脉冲或数码,换句话说,这些信号仅定义在离散时间上,则这样的系统称为离散时间系统,简称离散系统

采样控制系统

通常,把系统中的离散信号是脉冲序列形式的离散系统,称为采样控制系统或脉冲控制系统

一般说来,采样系统是对来自传感器的连续信息在某些规定的时间瞬时上取值。如果在有规律的间隔上,系统取到了离散信息,则这种采样称为周期采样

反之,如果信息之间的间隔是时变的,或随机的,则称为非周期采样,或随机采样

  • 通常,测量元件、执行元件和被控对象是模拟元件,其输入和输出是连续信号,即时间上和幅值上都连续的信号,称为模拟信号
  • 而控制器中的脉冲元件,其输入和输出为脉冲序列,即时间上离散而幅值上连续的信号,称为离散模拟信号

为了使两种信号在系统中能相互传递,在连续信号和脉冲序列之间要用采样器,而在脉冲序列和连续信号之间用保持器,以实现两种信号的转换

采样器和保持器,是采样控制系统中的两个特殊环节

信号采样和复现

在采样控制系统中,把连续信号转变为脉冲序列的过程称为采样过程,简称采样。实现采样的装置称为采样器或称采样开关。

用 T 表示采样周期,单位为s; f s = 1 / T f_s=1/T fs=1/T 表示采样频率,单位为1/s; w s = 2 π f s = 2 π / T w_s=2\pi{f_s}=2\pi/T ws=2πfs=2π/T 表示采样角频率,单位为rad/s

在实际应用中,采样开关多为电子开关,闭合时间极短,采样持续时间 τ \tau τ 远小于采样周期 T, 也远小于系统连续部分的最大时间常数。为了简化系统的分析,可认为 τ \tau τ 趋于零,即把采样器的输出近似看成一串强度等于矩形脉冲面积的理想脉冲 e ∗ ( t ) e^*(t) e(t)

在采样控制系统中,把脉冲序列转变为连续信号的过程称为信号复现过程。实现复现过程的装置称为保持器。

数字控制系统

数字序列形式的离散系统,称为数字控制系统或计算机控制系统

数字控制系统,是一种以数字计算机为控制器去控制具有连续工作状态的被控对象的闭环控制系统。因此,数字控制系统包括工作于离散状态下的数字计算机和工作于连续状态下的被控对象两大部分

计算机作为系统的控制器,其输入和输出只能是二进制编码的数字信号,即在时间上和幅值上都离散的信号,而系统中被控对象和测量元件的输入和输出是连续信号,所以在计算机控制系统中,需要应用A/D (模/数) 和 D/A (数/模) 转换器,以实现两种信号的转换

数字计算机在对系统进行实时控制时,每隔 T 秒进行一次控制修正,T 为采样周期

在每个采样周期中,控制器要完成对于连续信号的采样编码 (即A/D 过程) 和按控制律进行的数码运算,然后将计算结果由输出寄存器经解码网络将数码转换成连续信号(即D/A过程)

因此,A/D 转换器和 D/A 转换器是计算机控制系统中的两个特殊环节

Z 变换

在这里插入图片描述

z 变换性质

  • 线性定理
    Z [ e 1 ( t ) ± e 2 ( t ) ] = E 1 ( z ) ± E 2 ( z ) Z [ a e ( t ) ] = a E ( z ) \mathscr{Z}[e_1(t){\pm}e_2(t)]=E_1(z){\pm}E_2(z)\\ \mathscr{Z}[ae(t)]=aE(z) Z[e1(t)±e2(t)]=E1(z)±E2(z)Z[ae(t)]=aE(z)

  • 实数位移定理

    实数位移定理又称平移定理。实数位移的含意,是指整个采样序列在时间轴上左右平移若干采样周期,其中向左平移为超前,向右平移为滞后
    Z [ e ( t − k T ) ] = z − k E ( z ) Z [ e ( t + k T ) ] = z k [ E ( z ) − ∑ n = 0 k − 1 e ( n T ) z − n ] \mathscr{Z}[e(t-kT)]=z^{-k}E(z)\\ \mathscr{Z}[e(t+kT)]=z^k[E(z)-\sum^{k-1}_{n=0}e(nT)z^{-n}] Z[e(tkT)]=zkE(z)Z[e(t+kT)]=zk[E(z)n=0k1e(nT)zn]

  • 复数位移定理
    e ∓ a t e ( t ) = E ( z e ∓ a T ) \mathscr{e^{{\mp}at}e(t)}=E(ze^{{\mp}aT}) eate(t)=E(zeaT)

  • 终值定理

    如果函数e(t)的 z 变换为 E(z),函数序列 e(nT) 为有限值 (n=0,1,2,…),且极限存在,则函数序列的终值:
    lim ⁡ n → ∞ e ( n T ) = lim ⁡ z → 1 ( z − 1 ) E ( z ) 也可表示为  e ( ∞ ) = lim ⁡ n → ∞ e ( n T ) = lim ⁡ z → 1 ( 1 − z − 1 ) E ( z ) \lim_{n\to\infty}e(nT)=\lim_{z\to1}(z-1)E(z)\\ \text{也可表示为 }e(\infty)=\lim_{n\to\infty}e(nT)=\lim_{z\to1}(1-z^{-1})E(z) nlime(nT)=z1lim(z1)E(z)也可表示为 e()=nlime(nT)=z1lim(1z1)E(z)

  • 卷积定理
    g ( n T ) = x ( n T ) ∗ y ( n T ) G ( z ) = X ( z ) . Y ( z ) g(nT)=x(nT)*y(nT)\\ G(z)=X(z).Y(z) g(nT)=x(nT)y(nT)G(z)=X(z).Y(z)

z 反变换

离散系统的稳定性与稳态误差

K p = lim ⁡ z → 1 [ 1 + G ( z ) ] K_p=\lim_{z\to1}[1+G(z)] Kp=z1lim[1+G(z)]

称为静态位置误差系数

  • 若 G(z) 没有 z = 1 的极点,则 K p ≠ ∞ K_p\ne\infty Kp=,从而 e ( ∞ ) ≠ 0 e(\infty)\ne0 e()=0,这样的系统称为0型离散系统
  • 若 G(z) 有一个或一个以上 z = 1 的极点,则 K p = ∞ K_p=\infty Kp=,从而 e ( ∞ ) = 0 e(\infty)=0 e()=0,这样的系统相应称为 I 型或 I 到以上的离散系统

K v = lim ⁡ z → 1 ( z − 1 ) G ( z ) K_v=\lim_{z\to1}(z-1)G(z) Kv=z1lim(z1)G(z)
称为静态速度误差系数。0 型系统的 K v = 0 K_v=0 Kv=0 ,I 型系统的 K v K_v Kv 为有限值,Ⅱ型和 II 型以上系统的 K v = ∞ K_v=\infty Kv=

K a = lim ⁡ z → 1 ( z − 1 ) 2 G ( z ) K_a=\lim_{z\to1}(z-1)^2G(z) Ka=z1lim(z1)2G(z)
称为静态加速度误差系数。0 型及 I 型系统的 K a = 0 K_a=0 Ka=0 ,Ⅱ型系统的 K a K_a Ka 为常值,III 型及 III 型以上系统的 K a = ∞ K_a=\infty Ka=

在这里插入图片描述

  • 30
    点赞
  • 169
    收藏
    觉得还不错? 一键收藏
  • 22
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 22
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值