自控理论与设计

第一章 绪论

自动控制理论通常可分为经典控制理论现代控制理论智能控制理论

控制理论发展和概括

一、经典控制理论

传递函数为基础,研究单输入-单输出控制系统的分析和设计。

1. 线性(连续)控制理论

方法核心
时域分析法微分方程求解,分析系统的时间响应
频率特性法基于频率响应(Bode 图、Nyquist 图)
根轨迹法闭环极点随增益变化的轨迹
为什么这三种方法要求「线性」?

线性系统满足叠加原理T[αx1+βx2]=αT[x1]+βT[x2]T[\alpha x_1+\beta x_2]=\alpha T[x_1]+\beta T[x_2]),这意味着响应可以拆解叠加。三种方法本质上都是对这个性质的利用:

时域分析法: 线性 ODE 的解 = 齐次解(特征根组合 eλte^{\lambda t})+ 特解。非线性方程没有这种分解——x1(t)x_1(t)x2(t)x_2(t) 分别是解,x1+x2x_1+x_2 不一定是解。

频率特性法: 只有线性系统在输入 sinωt\sin\omega t 时输出才保持同频(仅变幅值和相位)。非线性系统会产生 sin2ωt\sin 2\omega tsin3ωt\sin 3\omega t谐波——输出不是单一频率,Bode 图和 Nyquist 图就没法画了。

根轨迹法: 依赖传递函数 G(s)G(s),而传递函数定义为零初始条件下输出与输入的拉氏变换之比,只有对线性定常系统才成立。非线性系统根本不存在传递函数。

方法依赖的线性性质非线性系统的后果
时域分析法特征根叠加无法解耦,无通式
频率特性法正弦输入 → 同频正弦输出产生谐波,无法画 Bode
根轨迹法存在传递函数 G(s)G(s)传递函数无定义

这就是为什么非线性系统需要独立的分析工具(描述函数法、相平面法、李雅普诺夫直接法)。

2. 非线性控制理论

方法核心
相平面法二阶系统状态轨迹的几何分析
描述函数法谐波线性化,分析非线性系统的稳定性
李雅普诺夫直接法不求解方程,直接判稳
波波夫法频域非线性稳定性判据

3. 采样控制理论(要求线性定常)

  • Z 变换法:离散系统的分析与设计

辨析:线性控制理论 vs 采样控制理论

教材中”线性控制理论”准确的叫法是连续线性控制理论——它既要求线性,也假设时间连续。而”采样控制理论”研究的是离散时间系统,但 Z 变换同样以线性定常为前提。

连续时间离散时间
线性线性控制理论(拉氏变换 G(s)G(s)采样控制理论(Z 变换 G(z)G(z)
非线性描述函数法 / 相平面法 / 李雅普诺夫法离散描述函数 / 离散李雅普诺夫
  • 采样控制理论也要求线性——Z 变换只对线性定常系统有定义
  • “线性控制理论”隐含连续——三种方法(时域、频域、根轨迹)都在连续时间域
  • 非线性离散系统是更高级的话题(研究生阶段),经典教材不展开

二、现代控制理论

状态空间法为基础,研究多输入-多输出时变非线性等控制系统的分析和设计。

分支核心
状态空间分析法用状态方程描述系统内部行为
最优控制理论在约束下求最优控制律
系统辨识从数据建立系统数学模型
最佳估计(Kalman 滤波等)从噪声中估计系统状态

三、智能控制理论

人工智能理论为基础,研究具有模糊性、不确定性、不完全性、偶然性的系统。

分支核心
专家系统基于规则的推理控制
模糊控制处理模糊集合和模糊逻辑
学习控制系统通过经验自主改进

三者关系

控制理论发展脉络:

经典控制理论 ──▶ 现代控制理论 ──▶ 智能控制理论
  (1940s-60s)     (1960s-80s)      (1980s-至今)
  传递函数         状态空间          人工智能
  SISO系统         MIMO系统          不确定性系统

自控系统分类

按信号传递路径{开环闭环按信号传递路径\begin{cases} 开环 闭环 \end{cases} 按控制作用{恒值控制系统(自动镇定)随动系统(自动跟踪)程序控制系统按控制作用\begin{cases} 恒值控制系统(自动镇定) 随动系统(自动跟踪) 程序控制系统 \end{cases}

对控制系统的要求

  • 稳定性
  • 动态特性
  • 稳态特性 三者不能兼容

典型测试信号

实际系统输入具有不确定性,无法用解析式表达。引入典型测试信号统一标准,便于对比分析不同系统的性能。

1. 阶跃输入函数(Step)

检测系统突变的跟踪能力。阶跃函数频带很宽,跃变等价于宽频域内无数正弦信号的合成。

r(t)={0,t<0R,t0(1-2)r(t) = \begin{cases} 0, & t < 0 R, & t \ge 0 \end{cases} \qquad (1\text{-}2)

单位阶跃:R=1R=1,记作 u(t)u(t)。 起始时刻 τ\tau 的阶跃记为 u(tτ)u(t-\tau)

2. 斜坡输入函数(Ramp / 速度函数)

检测系统匀速运动的跟踪能力。斜坡函数 = 阶跃函数对 tt 的积分。

r(t)={0,t<0Rt,t0(1-3)r(t) = \begin{cases} 0, & t < 0 Rt, & t \ge 0 \end{cases} \qquad (1\text{-}3)

单位速度函数:R=1R=1,记作 tu(t)t \cdot u(t)

3. 抛物线输入函数(Parabolic / 加速度函数)

表示匀加速信号,由速度函数对 tt 积分而得。

r(t)={0,t<0[4pt]12Rt2,t0(1-4)r(t) = \begin{cases} 0, & t < 0 [4pt] \dfrac{1}{2}Rt^{2}, & t \ge 0 \end{cases} \qquad (1\text{-}4)

单位加速度函数:R=1R=1,记作 12t2u(t)\frac{1}{2}t^{2} \cdot u(t)。实际很少用比这变化更快的信号。

4. 脉冲输入函数(Pulse)

模拟瞬时冲击。实际脉冲宽度 hh 应极小,一般要求 h<0.1Th < 0.1TTT 为系统时间常数)。

r(t)={0,t<0,th[4pt]Ah,0th(1-5)r(t) = \begin{cases} 0, & t<0, t\ge h [4pt] \dfrac{A}{h}, & 0 \le t \le h \end{cases} \qquad (1\text{-}5)

理想单位脉冲(A=1,h0A=1, h\to 0):

r(t)=δ(t),+δ(t)dt=1(1-6,1-7)r(t) = \delta(t), \quad \int_{-\infty}^{+\infty}\delta(t)dt = 1 \qquad (1\text{-}6, 1\text{-}7)

5. 正弦函数(Sinusoid)

检测系统不同频率下的响应特性。求得所有频率的正弦响应,可准确确定整个系统特性。

r(t)=Asin(ωt+ϕ)(1-8)r(t) = A\sin(\omega t + \phi) \qquad (1\text{-}8)
  • AA:振幅 ϕ\phi:相位移 ω\omega:振荡角频率

五种信号关系

信号数学关系检测性能
脉冲 δ(t)\delta(t)冲击(无穷短)系统全部动态特性
阶跃 u(t)u(t)脉冲的积分突变的跟踪能力
斜坡 tu(t)t \cdot u(t)阶跃的积分匀速运动的跟踪
抛物线 12t2u(t)\frac{1}{2}t^2 u(t)斜坡的积分匀加速跟踪
正弦 sinωt\sin\omega t不同频率稳态频率响应特性

阶跃→斜坡→抛物线逐级积分,变化逐渐加快。脉冲响应做拉氏变换即得传递函数 G(s)G(s)

第二章 控制系统的数学模型

  • 数学模型
{静态模型:在静态条件下(即变量各阶导数为零),描述变量间关系的代数方程。动态模型:用微分(差分)方程描述的变量在动态过程中的关系。\left\{ \begin{array}{l} \text{静态模型:在静态条件下(即变量各阶导数为零),描述变量间关系的代数方程。} \\ \text{动态模型:用微分(差分)方程描述的变量在动态过程中的关系。} \end{array} \right.
  • 数学模型表示形式
{图形表示:信号流图、方块图、频率特性图及 SIMULINK 框图描述等。数学表示{微分(差分)方程、代数方程、逻辑方程是基本的形式。传递函数或频率特性,适用于单输入单输出系统。状态空间描述,用于多变量系统及最优控制问题。数字计算机上的程序综合。\left\{ \begin{array}{l} \text{图形表示:信号流图、方块图、频率特性图及 SIMULINK 框图描述等。} \\[4pt] \text{数学表示}\left\{ \begin{array}{l} \text{微分(差分)方程、代数方程、逻辑方程是基本的形式。} \\[2pt] \text{传递函数或频率特性,适用于单输入单输出系统。} \\[2pt] \text{状态空间描述,用于多变量系统及最优控制问题。} \end{array} \right. \\[4pt] \text{数字计算机上的程序综合。} \end{array} \right.

LTS系统微分方程的一般形式(时域分析法)

andnc(t)dtn+an1dn1c(t)dtn1++a1dc(t)dt+a0c(t)=bmdmr(t)dtm+bm1dm1r(t)dtm1++b1dr(t)dt+b0r(t)(2-1)a_n \frac{\mathrm{d}^n c(t)}{\mathrm{d}t^n} + a_{n-1} \frac{\mathrm{d}^{n-1} c(t)}{\mathrm{d}t^{n-1}} + \cdots + a_1 \frac{\mathrm{d}c(t)}{\mathrm{d}t} + a_0 c(t) = b_m \frac{\mathrm{d}^m r(t)}{\mathrm{d}t^m} + b_{m-1} \frac{\mathrm{d}^{m-1} r(t)}{\mathrm{d}t^{m-1}} + \cdots + b_1 \frac{\mathrm{d}r(t)}{\mathrm{d}t} + b_0 r(t) \qquad (2\text{-}1)

式中:r(t)r(t)为输入量;c(t)c(t)为输出量。

在实际物理系统中,系数 aia_ibib_i 均为实数,并且一般有 nmn \ge m。这是因为一般物理系统均有质量、惯性或滞后的储能元件,所以 c(t)c(t) 的导数阶次一般比 r(t)r(t) 高(有时相等)。 对微分方程进行拉氏变化,可得到系统输出与输入之比,即系统的传递函数。

传递函数的性质

  • 系数均为实数
  • 只与系统或元件本身参数有关,与输入输出无关
  • 不能反映系统或元件的物理组成-不同元件可能有相同传递函数
  • 传递函数与系统单位脉冲响应为一组拉氏变换
  • 仅适用于线性定常系统

式 (2-7) 时间常数表达式中 ω\omega 的来源

式 (2-7) 是传递函数的时间常数表达式:

G(s)=Ki=1m1(τis+1)k=1m2[(sωnk)2+2ξnk(sωnk)+1]svj=1n1(Tjs+1)l=1n2[(sωpl)2+2ξpl(sωpl)+1](2-7)G(s) = \frac{K \prod\limits_{i=1}^{m_1}(\tau_i s + 1) \prod\limits_{k=1}^{m_2}\left[\left(\dfrac{s}{\omega_{nk}}\right)^2 + 2\xi_{nk}\left(\dfrac{s}{\omega_{nk}}\right) + 1\right]}{s^v \prod\limits_{j=1}^{n_1}(T_j s + 1) \prod\limits_{l=1}^{n_2}\left[\left(\dfrac{s}{\omega_{pl}}\right)^2 + 2\xi_{pl}\left(\dfrac{s}{\omega_{pl}}\right) + 1\right]} \qquad (2\text{-}7)

式中:KK 为稳态增益;m1+2m2=mm_1 + 2m_2 = mv+n1+2n2=nv + n_1 + 2n_2 = n

教材在此处跳跃了一步:没有解释 ωnk\omega_{nk}ωpl\omega_{pl} 从何而来。它们实际上来自零极点表达式 (2-8) 中共轭复根的标准化改写。


推导:共轭复根 \to ω\omegaξ\xi 的标准型

设系统有一对共轭复零点,其根为

zk=ξω±jω1ξ2-z_k = -\xi\omega \pm j\omega\sqrt{1-\xi^2}

ω\omega 为无阻尼固有频率,ξ\xi 为阻尼比)

还原为二阶因子:

(s+ξωjω1ξ2)(s+ξω+jω1ξ2)=s2+2ξωs+ω2(s + \xi\omega - j\omega\sqrt{1-\xi^2})(s + \xi\omega + j\omega\sqrt{1-\xi^2}) = s^2 + 2\xi\omega s + \omega^2

提取 ω2\omega^2

s2+2ξωs+ω2=ω2[(sω)2+2ξ(sω)+1]s^2 + 2\xi\omega s + \omega^2 = \omega^2\left[\left(\frac{s}{\omega}\right)^2 + 2\xi\left(\frac{s}{\omega}\right) + 1\right]

这个 ω2\omega^2 被吸进前面的增益 KK 中,因此式 (2-7) 中只保留中括号内的标准型。极点同理。


结构对照
零/极点类型零极点表达式 (2-8)时间常数表达式 (2-7)
实零点(s+zi)(s + z_i)(τis+1)(\tau_i s + 1)τi=1/zi\tau_i = 1/z_i
共轭复零点s2+2ξnkωnks+ωnk2s^2 + 2\xi_{nk}\omega_{nk}s + \omega_{nk}^2[(sωnk)2+2ξnk(sωnk)+1]\left[\left(\frac{s}{\omega_{nk}}\right)^2 + 2\xi_{nk}\left(\frac{s}{\omega_{nk}}\right) + 1\right]
实极点(s+pj)(s + p_j)(Tjs+1)(T_j s + 1)Tj=1/pjT_j = 1/p_j
共轭复极点s2+2ξplωpls+ωpl2s^2 + 2\xi_{pl}\omega_{pl}s + \omega_{pl}^2[(sωpl)2+2ξpl(sωpl)+1]\left[\left(\frac{s}{\omega_{pl}}\right)^2 + 2\xi_{pl}\left(\frac{s}{\omega_{pl}}\right) + 1\right]
原点处极点svs^vsvs^v

下标含义nn = zero(零点),pp = pole(极点);k,lk, l 为序号。


一句话总结

ωnk\omega_{nk}ωpl\omega_{pl} 分别是共轭复零点、共轭复极点对应的无阻尼固有频率。 式 (2-7) 本质上是把零极点按”实根 \to 一阶时间常数形式、复根 \to 二阶标准型”重新整理后的结果。


典型环节及其传递函数

出处:《自控理论与设计 徐薇莉》第二章 §2.2.4 · PDF第21–22页 · 经 OCR 识别

任何一个复杂系统都是由有限个典型环节组合而成的。典型环节有以下六种:

1. 比例环节

G(s)=K(2-15)G(s) = K \qquad (2\text{-}15)

式中:KK 为增益。 特点:输出量与输入量成比例,无失真和时间延迟。 实例:电子放大器、齿轮、电阻(或电位器)、感应式变送器等。

2. 惯性环节

G(s)=1Ts+1(2-16)G(s) = \frac{1}{Ts+1} \qquad (2\text{-}16)

式中:TT 为时间常数。 特点:含一个储能元件,对突变的输入其输出不能立即复现,输出无振荡。 实例:图 2-1 所示的 RCRC 网络,直流伺服电动机从输入电压到转速的传递函数也包含这种环节。

3. 微分环节

  • 理想微分:
G(s)=Ks(2-17)G(s) = Ks \qquad (2\text{-}17)
  • 一阶微分:
G(s)=τs+1(2-18)G(s) = \tau s + 1 \qquad (2\text{-}18)
  • 二阶微分:
G(s)=τ2s2+2ξτs+1(2-19)G(s) = \tau^2 s^2 + 2\xi\tau s + 1 \qquad (2\text{-}19)

特点:输出量正比于输入量变化的速度,能预示输入信号的变化趋势。 实例:测速发电机输出电压与输入角度间的传递函数即为微分环节。

4. 积分环节

G(s)=1s(2-20)G(s) = \frac{1}{s} \qquad (2\text{-}20)

特点:输出量与输入量的积分成比例,当输入消失后输出仍具记忆功能。 实例:电动机角速度与角度间的传递函数,模拟计算机中的积分器等。

5. 振荡环节

G(s)=ωn2s2+2ξωns+ωn2=1T2s2+2ξTs+1(2-21)G(s) = \frac{\omega_n^2}{s^2 + 2\xi\omega_n s + \omega_n^2} = \frac{1}{T^2 s^2 + 2\xi Ts + 1} \qquad (2\text{-}21)

式中:ξ\xi 称为阻尼比系数(0ξ<10 \le \xi \lt 1);ωn\omega_n 称为自然振荡角频率(无阻尼振荡角频率);T=1ωnT = \dfrac{1}{\omega_n}。 特点:环节中有两个独立的储能元件,并可进行能量交换,其输出出现振荡。 实例:RLCRLC 电路的输出与输入电压间的传递函数。

6. 纯时间延迟环节(时滞环节)

c(t)=r(tτ),G(s)=eτs(2-22)c(t) = r(t-\tau), \quad G(s) = e^{-\tau s} \qquad (2\text{-}22)

式中:τ\tau 为延迟时间。 特点:输出量能准确复现输入量,但须延迟一个固定的时间间隔。 实例:管道的温度、流量等物理量的控制问题的数学模型就包含有延迟环节。

拉普拉斯变换性质及使用

基本定义、常用变换对与核心性质速查;自控建模与微分方程求解常用。

一、基本定义

F(s)=L{f(t)}=0f(t)estdt,s=σ+jωF(s)=\mathcal{L}\{f(t)\}=\int_{0^-}^{\infty} f(t)\,e^{-st}\,\mathrm{d}t,\quad s=\sigma+j\omega

逆变换:

f(t)=L1{F(s)}=12πjσjσ+jF(s)estdsf(t)=\mathcal{L}^{-1}\{F(s)\}=\frac{1}{2\pi j}\int_{\sigma-j\infty}^{\sigma+j\infty} F(s)\,e^{st}\,\mathrm{d}s

二、常用变换对

序号原函数 f(t)f(t)t0t\geq 0象函数 F(s)F(s)收敛域
1单位脉冲 δ(t)\delta(t)11全平面
2单位阶跃 u(t)u(t)1(t)1(t)1s\dfrac{1}{s}Re(s)>0\mathrm{Re}(s)>0
3tt1s2\dfrac{1}{s^2}Re(s)>0\mathrm{Re}(s)>0
4tnt^nn=0,1,2,n=0,1,2,\ldotsn!sn+1\dfrac{n!}{s^{n+1}}Re(s)>0\mathrm{Re}(s)>0
5tnt^nn>1n>-1,一般实数)Γ(n+1)sn+1\dfrac{\Gamma(n+1)}{s^{n+1}}Re(s)>0\mathrm{Re}(s)>0
6eate^{-at}1s+a\dfrac{1}{s+a}Re(s)>a\mathrm{Re}(s)>-a
7teatt\,e^{-at}1(s+a)2\dfrac{1}{(s+a)^2}Re(s)>a\mathrm{Re}(s)>-a
8tneatt^n e^{-at}n!(s+a)n+1\dfrac{n!}{(s+a)^{n+1}}Re(s)>a\mathrm{Re}(s)>-a
9sin(ωt)\sin(\omega t)ωs2+ω2\dfrac{\omega}{s^2+\omega^2}Re(s)>0\mathrm{Re}(s)>0
10cos(ωt)\cos(\omega t)ss2+ω2\dfrac{s}{s^2+\omega^2}Re(s)>0\mathrm{Re}(s)>0
11eatsin(ωt)e^{-at}\sin(\omega t)ω(s+a)2+ω2\dfrac{\omega}{(s+a)^2+\omega^2}Re(s)>a\mathrm{Re}(s)>-a
12eatcos(ωt)e^{-at}\cos(\omega t)s+a(s+a)2+ω2\dfrac{s+a}{(s+a)^2+\omega^2}Re(s)>a\mathrm{Re}(s)>-a
13sinh(ωt)\sinh(\omega t)ωs2ω2\dfrac{\omega}{s^2-\omega^2}Re(s)>ω\mathrm{Re}(s)>\lvert\omega\rvert
14cosh(ωt)\cosh(\omega t)ss2ω2\dfrac{s}{s^2-\omega^2}Re(s)>ω\mathrm{Re}(s)>\lvert\omega\rvert
151eat1-e^{-at}as(s+a)\dfrac{a}{s(s+a)}Re(s)>0\mathrm{Re}(s)>0
161ba(eatebt)\dfrac{1}{b-a}\bigl(e^{-at}-e^{-bt}\bigr)1(s+a)(s+b)\dfrac{1}{(s+a)(s+b)}Re(s)>max(a,b)\mathrm{Re}(s)>\max(-a,-b)
171ba(bebtaeat)\dfrac{1}{b-a}\bigl(be^{-bt}-ae^{-at}\bigr)s(s+a)(s+b)\dfrac{s}{(s+a)(s+b)}Re(s)>max(a,b)\mathrm{Re}(s)>\max(-a,-b)
18ω1ζ2eζωtsin(ω1ζ2t)\dfrac{\omega}{\sqrt{1-\zeta^2}}e^{-\zeta\omega t}\sin\bigl(\omega\sqrt{1-\zeta^2}\,t\bigr)ω2s2+2ζωs+ω2\dfrac{\omega^2}{s^2+2\zeta\omega s+\omega^2}Re(s)>ζω\mathrm{Re}(s)>-\zeta\omega
19eat[cos(ωt)aωsin(ωt)]e^{-at}\bigl[\cos(\omega t)-\dfrac{a}{\omega}\sin(\omega t)\bigr]s(s+a)2+ω2\dfrac{s}{(s+a)^2+\omega^2}Re(s)>a\mathrm{Re}(s)>-a

三、核心性质

性质时域复频域
线性αf1(t)+βf2(t)\alpha f_1(t)+\beta f_2(t)αF1(s)+βF2(s)\alpha F_1(s)+\beta F_2(s)
时移f(ta)u(ta)f(t-a)u(t-a)a>0a>0easF(s)e^{-as}F(s)
频移eatf(t)e^{-at}f(t)F(s+a)F(s+a)
尺度变换f(at)f(at)a>0a>01aF(sa)\dfrac{1}{a}F\bigl(\dfrac{s}{a}\bigr)
时域微分f(t)f'(t)sF(s)f(0)sF(s)-f(0^-)
f(t)f''(t)s2F(s)sf(0)f(0)s^2 F(s)-s f(0^-)-f'(0^-)
f(n)(t)f^{(n)}(t)snF(s)k=0n1sn1kf(k)(0)s^n F(s)-\displaystyle\sum_{k=0}^{n-1}s^{n-1-k}f^{(k)}(0^-)
时域积分0tf(τ)dτ\displaystyle\int_0^t f(\tau)\,\mathrm{d}\tau1sF(s)\dfrac{1}{s}F(s)
tf(τ)dτ\displaystyle\int_{-\infty}^t f(\tau)\,\mathrm{d}\tau1sF(s)+1s0f(τ)dτ\dfrac{1}{s}F(s)+\dfrac{1}{s}\displaystyle\int_{-\infty}^{0^-}f(\tau)\,\mathrm{d}\tau
复频域微分tf(t)t\,f(t)dF(s)ds-\dfrac{\mathrm{d}F(s)}{\mathrm{d}s}
tnf(t)t^n f(t)(1)ndnF(s)dsn(-1)^n\dfrac{\mathrm{d}^n F(s)}{\mathrm{d}s^n}
复频域积分f(t)t\dfrac{f(t)}{t}sF(σ)dσ\displaystyle\int_s^{\infty} F(\sigma)\,\mathrm{d}\sigma
时域卷积f1f2=0tf1(τ)f2(tτ)dτf_1*f_2=\displaystyle\int_0^t f_1(\tau)f_2(t-\tau)\,\mathrm{d}\tauF1(s)F2(s)F_1(s)\cdot F_2(s)
复频域卷积f1(t)f2(t)f_1(t)\cdot f_2(t)12πjF1(s)F2(s)\dfrac{1}{2\pi j}F_1(s)*F_2(s)
初值定理f(0+)=limssF(s)f(0^+)=\displaystyle\lim_{s\to\infty}sF(s)
终值定理f()=lims0sF(s)f(\infty)=\displaystyle\lim_{s\to 0}sF(s)(需极点均在左半平面,虚轴上至多原点单极点)
周期函数f(t)=f(t+T)f(t)=f(t+T)11esT0Tf(t)estdt\dfrac{1}{1-e^{-sT}}\displaystyle\int_0^T f(t)e^{-st}\,\mathrm{d}t

四、记忆口诀

  • 阶跃倒 ss,斜坡倒 s2s^2nn 次幂配阶乘;
  • 指数平移分母加,正弦余弦分母二次;
  • 时移乘指数,微分乘 ss 减初值,积分除以 ss
  • 卷积变乘积,初值终值看 ss 极限。

常用概念

信号流图相关

概念含义
开通路与任一节点仅相遇一次的通路
前向通路从输入节点到输出节点的通路
通路增益开通路 / 前向通路上各支路增益之积
闭通路(回路)起始与终止为同一节点,且与其他节点仅相遇一次的通路
环增益回路上各支路增益之积

电路元件

电容

iC=CduCdt    {UC(s)=IC(s)C1s+uC(0)sIC(s)=sCUC(s)CuC(0)i_C = C\frac{\mathrm{d}u_C}{\mathrm{d}t} \;\Longleftrightarrow\; \begin{cases} U_C(s)=\dfrac{I_C(s)}{C}\cdot\dfrac{1}{s}+\dfrac{u_C(0^-)}{s} \\[6pt] I_C(s)=sC\cdot U_C(s)-C\cdot u_C(0^-) \end{cases}

电感

uL=LdiLdt    {UL(s)=sLIL(s)LiL(0)IL(s)=UL(s)L1s+iL(0)su_L = L\frac{\mathrm{d}i_L}{\mathrm{d}t} \;\Longleftrightarrow\; \begin{cases} U_L(s)=sL\cdot I_L(s)-L\cdot i_L(0^-) \\[6pt] I_L(s)=\dfrac{U_L(s)}{L}\cdot\dfrac{1}{s}+\dfrac{i_L(0^-)}{s} \end{cases}

梅森公式

亦称梅逊(Mason)增益公式。由信号流图(或方块图)直接求输入到输出的总增益(传递函数),不必逐步化简。

P=1Δk=1nPkΔkC(s)R(s)=kPkΔkΔP=\frac{1}{\Delta}\sum_{k=1}^{n} P_k\Delta_k \qquad\text{或}\qquad \frac{C(s)}{R(s)}=\frac{\sum_{k} P_k\Delta_k}{\Delta}
符号含义
PP从源节点(输入)到阱节点(输出)的总增益 / 传递函数
nn前向通路条数
PkP_kkk 条前向通路的通路增益(各支路增益之积)
Δ\Delta流图特征式(见下)
Δk\Delta_k余子式:把 Δ\Delta 中与第 kk 条前向通路相接触的回路项全部去掉后的余式;若该前向通路与所有回路都接触,则 Δk=1\Delta_k=1

特征式:

Δ=1La+LbLcLdLeLf+\Delta=1-\sum L_a+\sum L_b L_c-\sum L_d L_e L_f+\cdots
  • La\sum L_a:所有单独回路的环增益之和
  • LbLc\sum L_b L_c:所有两两互不接触回路,环增益乘积之和
  • LdLeLf\sum L_d L_e L_f:所有三三互不接触回路,环增益乘积之和
  • 符号规律:取 mm 个互不接触回路时,前面乘 (1)m(-1)^m

用法:数清前向通路与回路 → 算 Δ\Delta → 对每条前向通路求 PkP_kΔk\Delta_k → 代入求和。求扰动到输出时,把扰动当源节点即可;特征式 Δ\Delta 与选哪条前向通路无关(同一张图共用一个 Δ\Delta)。

第四章 控制系统的时域分析

劳斯判据

出处:徐薇莉《自控理论与设计》§4.5;胡寿松《自动控制原理》§3-5(劳斯–赫尔维茨)。判的是闭环特征根是否全在左半 ss 平面,不直接求根。

1. 用之前先看必要条件

闭环特征方程(系数均为实数)

D(s)=a0sn+a1sn1++an1s+an=0D(s)=a_0 s^n+a_1 s^{n-1}+\cdots+a_{n-1}s+a_n=0

无正实部根的必要条件(两本一致):各项系数同号,且无一为零(不缺项)。

  • 不满足 \Rightarrow 必不稳定(可直接判,不必列表)
  • 满足 \Rightarrow 还不够,必须再排劳斯表(或赫尔维茨行列式)

低阶时胡寿松给出的充要(系数已全正):

阶次再补一条
n=2n=2系数全正即可
n=3n=3a1a2a0a3>0a_1 a_2-a_0 a_3>0
n=4n=4a1a2a0a3>0a_1 a_2-a_0 a_3>0,且 Δ3>a1a4/a3\Delta_3>a_1 a_4/a_3

2. 劳斯表怎么排

前两行由系数直接填:奇数项一行、偶数项一行;空位补 0。以后各行由上面两行按「交叉相乘再除」得到。以第 3 行第 1 列为例:

b1=1a1a0a2a1a3=a1a2a0a3a1b_1=-\frac{1}{a_1}\begin{vmatrix} a_0 & a_2 \\ a_1 & a_3 \end{vmatrix}=\frac{a_1 a_2-a_0 a_3}{a_1}

其余元算法相同。一直排到 s0s^0 行。

  • 行数 =n+1=n+1;最后两行每行通常只有一个元素,s0s^0 行首元等于 ana_n
  • 一行可同乘 / 同除一个正数,不影响稳定性判断(上交 2021 填空)
  • 只用闭环特征方程;开环传函不能直接拿来排表

3. 判据(充分必要条件)

徐薇莉:首列非零且不改变符号。 胡寿松:在 a0>0a_0>0 约定下,第一列各元均为正

同时:

  • 第一列变号次数 = 右半平面根(正实部根)的个数
  • 第一列全正且无特殊情况 \Rightarrow 渐近稳定
第一列结论
全正、无 0渐近稳定
有负(变号)不稳定;变号几次就有几个右半平面根
出现 0表排不下去;系统不是渐近稳定(见下)

劳斯与赫尔维茨实质相同:劳斯表第一列各元与顺序赫尔维茨主子式 Δi\Delta_i 成比例。高阶算行列式麻烦时用劳斯表。

4. 两种特殊情况

出现 0 时先处理后才能数变号。徐薇莉:某行首元为 0 或全零行 \Rightarrow 系统必定不稳定(含临界稳定)

情况 1:某行只有第一列是 0,同行其余不全为 0

下一行会除以 0,出现 \infty

教材处理
徐薇莉用很小的正数 ε\varepsilon 代这个 0,继续排,再看首列变号
胡寿松用因子 (s+a)(s+a)a>0a>0)乘原特征方程,对新方程再列表(可避开该情况)

不表示一定有虚轴根,只表示原表算不下去。

情况 2:整行全为 0

存在对原点对称的根(可同时有几种):

  • 共轭纯虚根 ±jω\pm j\omega
  • 符号相反的实根 ±σ\pm\sigma
  • 对原点对称的共轭复根

处理(两本相同):

  1. 上一行系数写辅助方程 F(s)=0F(s)=0(次数一般为偶数)
  2. ss 求导,用 dF/ds\mathrm{d}F/\mathrm{d}s 的系数替换全零行
  3. 继续列表;对称根由 F(s)=0F(s)=0 解出

补充:

  • 全零行一般出现在 ss奇次行
  • 辅助方程解出的正根,已包含在首列变号次数里,不要重复加
  • 上交 2013 判断:「结束前有一行全 0 \Rightarrow 不稳定」——按「稳定 = 渐近稳定」算对(至少临界稳定)

5. 适用条件(徐薇莉注意点)

  • 特征方程必须是实系数代数方程
  • 系数中有复数,或含 eτse^{-\tau s} 等指数,不能直接用劳斯(延迟系统要先近似或改用奈氏)
  • 稳定性由结构和参数决定,与输入无关

6. 三个标准用法

(1) 判稳 / 数右半平面根 先看必要条件,再列表,处理特殊情况后看首列。

(2) 求使系统稳定的参数范围(上交常考) 特征式含 K,TK,T 等:令第一列各元 >0>0,并联立系数全正,画出 KKTT 稳定域。临界值常对应某行首元 =0=0(辅助方程 \Rightarrow 虚轴交点)。

(3) 相对稳定性 / 给定稳定度(两本都有) 劳斯本身看不出根离虚轴有多远。令

s=zσ(σ>0)s=z-\sigma\quad(\sigma>0)

代入原特征方程,对 zz 再排劳斯。若新表首列全正,则全部闭环极点在直线 s=σs=-\sigma(胡寿松称 σ\sigma 为给定稳定度)。

根轨迹与虚轴交点:令含 KK 的首元为 0,用辅助方程,与本条同一套路。

7. 做题顺序

  1. 写出闭环 D(s)=0D(s)=0,化成降幂、实系数
  2. 必要条件:同号、不缺项;否则已不稳定
  3. 列表;遇 0 先分清「仅首元 0」还是「全零行」
  4. 首列全正 \Rightarrow 渐近稳定;变号次数 = 右半平面根个数
  5. 含参:第一列 >0>0 解不等式;要根在 s=σs=-\sigma 左侧则先平移再列表

一句:劳斯不求根,只看闭环特征式系数——必要条件先筛,首列定稳与右根个数;第一列 0 先处理,全零行必有对称根。

稳态误差

出处:徐薇莉 §4.8;胡寿松 §3-6。只对已稳定的系统谈 esse_{ss};先劳斯/奈氏判稳,再算误差。

1. 定义与终值定理

开环(时间常数形式,vv 为型别 = 积分环节个数)

G(s)H(s)=Ki=1m(τis+1)svj(Tjs+1)G(s)H(s)=\frac{K\prod_{i=1}^{m}(\tau_i s+1)}{s^{v}\prod_{j}(T_j s+1)}

单位反馈(或一般 HH)时,误差传递函数与终值:

E(s)R(s)=11+G(s)H(s),ess=limte(t)=lims0sE(s)=lims0sR(s)1+G(s)H(s)\frac{E(s)}{R(s)}=\frac{1}{1+G(s)H(s)},\qquad e_{ss}=\lim_{t\to\infty}e(t)=\lim_{s\to 0}sE(s)=\lim_{s\to 0}\frac{sR(s)}{1+G(s)H(s)}

KK 为开环增益。型别 v=0,1,2,v=0,1,2,\ldots 分别称 0 型、I 型、II 型。胡寿松:阶跃下无原理性静差称无差系统;0 型有静差,I 型为一阶无差度,II 型为二阶无差度。

2. 三个静态误差系数

系数定义典型输入稳态误差
位置 KpK_plims0G(s)H(s)\displaystyle\lim_{s\to 0}G(s)H(s)r=R1(t), R(s)=R/sr=R\cdot 1(t),\ R(s)=R/sessp=R1+Kpe_{ssp}=\dfrac{R}{1+K_p}
速度 KvK_vlims0sG(s)H(s)\displaystyle\lim_{s\to 0}s\,G(s)H(s)r=Rt, R(s)=R/s2r=Rt,\ R(s)=R/s^2essv=RKve_{ssv}=\dfrac{R}{K_v}
加速度 KaK_alims0s2G(s)H(s)\displaystyle\lim_{s\to 0}s^2 G(s)H(s)r=12Rt2, R(s)=R/s3r=\dfrac{1}{2}Rt^2,\ R(s)=R/s^3essa=RKae_{ssa}=\dfrac{R}{K_a}

由指定允许误差反求系数(徐薇莉):

Kp=Resspessp,Kv=Ressv,Ka=RessaK_p=\frac{R-e_{ssp}}{e_{ssp}},\qquad K_v=\frac{R}{e_{ssv}},\qquad K_a=\frac{R}{e_{ssa}}

速度/加速度误差:不是速度、加速度本身有偏差,而是斜坡、抛物线输入下,稳态时位置上仍差一块。I 型跟斜坡:输出速度能跟上,但位置差 R/KvR/K_v

3. 型别速查

型别 vvKpK_p阶跃 esspe_{ssp}KvK_v斜坡 essve_{ssv}KaK_a抛物线 essae_{ssa}
0KKR/(1+K)R/(1+K)00\infty00\infty
I\infty00KKR/KR/K00\infty
II\infty00\infty00KKR/KR/K
\ge III\infty00\infty00\infty00

要阶跃无静差:至少 I 型。要斜坡无差:至少 II 型。型别不够 ess=\Rightarrow e_{ss}=\infty(与再调 KK 无关,上交 2025 选择题同理)。

4. 组合输入

线性系统可叠加。例如 r(t)=A+Bt+12Ct2r(t)=A+Bt+\dfrac{1}{2}Ct^2

ess=essp(A)+essv(B)+essa(C)e_{ss}=e_{ssp}(A)+e_{ssv}(B)+e_{ssa}(C)

某一项为 \infty,总误差就是 \infty

5. 输入误差与输出误差(非单位反馈)

比较点上的量才是 Kp,Kv,KaK_p,K_v,K_a 直接给出的(常见、默认):

E(s)=R(s)H(s)C(s)E(s)=R(s)-H(s)C(s)

H(s)=KhH(s)=K_h 为常数时,输出端希望值是 r/Khr/K_h(稳态比较点要平衡:Khc希望=rK_h c_{\text{希望}}=r),输出误差

E(s)=R(s)KhC(s)=E(s)KhE'(s)=\frac{R(s)}{K_h}-C(s)=\frac{E(s)}{K_h}
  • Kh=1K_h=1(单位反馈):E=EE'=E,两处是同一个量
  • Kh1K_h\neq 1:公式仍是 ess=ess/Khe'_{ss}=e_{ss}/K_h,不是换了一套;只是希望输出和 Kp=KKhK_p=K\cdot K_h 都随 KhK_h
  • 胡寿松例 3-13:G=10/(s+1)G=10/(s+1),阶跃下 ess=1/(1+10Kh)e_{ss}=1/(1+10K_h)Kh=1K_h=1ess=1/11e'_{ss}=1/11Kh=0.1K_h=0.1ess=5e'_{ss}=5

6. 注意点

  • 先稳定:不稳定谈 esse_{ss} 无意义(终值定理前提:极点在左半平面,虚轴至多原点单极点)。一阶、系数全正、参数给死了可一句话点出极点在左半平面;三阶及以上或 KK 待定要劳斯
  • 系数用 G(s)H(s)G(s)H(s)(开环);型别看开环积分个数;判稳看闭环 1+GH=01+GH=0
  • 扰动通道另写 En(s)/N(s)E_n(s)/N(s),不能直接套 Kp,Kv,KaK_p,K_v,K_a。令 R=0R=0 时希望输出增量为 0(扰动不该把 cc 拽走),故 En=CnE_n=-C_n
  • 测速、局部反馈会改等效开环,KvK_v 可能变小、斜坡误差变大(见下节)
  • 静态系数只管 tt\to\infty常值误差,不管超调、调节时间

7. 做题顺序

  1. 判稳
  2. 写成时间常数型,读出 vvKK
  3. 对输入类型查表,或直接
lims0sE(s)\lim_{s\to 0}sE(s)
  1. 组合输入则分行叠加
  2. 问输出端误差且 H=KhH=K_h 时,再除 KhK_h

一句:先稳再算;Kp,Kv,KaK_p,K_v,K_a 分别是 G(s)H(s)G(s)H(s)s0s\to 0 时乘 1,s,s21,s,s^2;型别够则对应误差为 0,不够则 \infty。非单位反馈先算比较点 E=RHCE=R-HC,输出端再除 KhK_h

动态误差

出处:徐薇莉 §4.10;胡寿松 §3-6 动态误差系数。静态系数 Kp,Kv,KaK_p,K_v,K_a 只能给出 tt\to\infty 时误差是有限、0 还是 \infty,看不出误差随时间怎么变,也不能方便地处理任意输入。动态误差用误差级数描述稳态阶段的 e(t)e(t)

1. 误差传递函数在 s=0s=0 邻域展开

Φe(s)=E(s)R(s)=11+G(s)H(s)\Phi_e(s)=\frac{E(s)}{R(s)}=\frac{1}{1+G(s)H(s)}

s=0s=0 处泰勒展开(该级数只在 s=0s=0 邻域收敛,对应时域 tt\to\infty):

Φe(s)=Φe(0)+Φe(0)s+12!Φe(0)s2+\Phi_e(s)=\Phi_e(0)+\Phi_e'(0)\,s+\frac{1}{2!}\Phi_e''(0)\,s^2+\cdots

反拉氏后得到误差级数(徐薇莉称它为动态误差的稳态项 ess(t)e_{ss}(t),一般不等于含暂态的完整 e(t)e(t)):

ess(t)=Φe(0)r(t)+Φe(0)r(t)+12!Φe(0)r(t)+e_{ss}(t)=\Phi_e(0)\,r(t)+\Phi_e'(0)\,r'(t)+\frac{1}{2!}\Phi_e''(0)\,r''(t)+\cdots

稳态误差是这条级数的极限:

ess=limtess(t)e_{ss}=\lim_{t\to\infty}e_{ss}(t)

完整误差可写成 e(t)=et(t)+ess(t)e(t)=e_t(t)+e_{ss}(t)(暂态 + 稳态项)。

2. 动态误差系数(徐薇莉记法)

名称定义与展开系数
动态位置 K0K_0K0=1/Φe(0)K_0=1/\Phi_e(0)Φe(0)=1/K0\Phi_e(0)=1/K_0
动态速度 K1K_1K1=1/Φe(0)K_1=1/\Phi_e'(0)Φe(0)=1/K1\Phi_e'(0)=1/K_1
动态加速度 K2K_2K2=1/(Φe(0)/2!)K_2=1\big/\bigl(\Phi_e''(0)/2!\bigr)Φe(0)/2!=1/K2\Phi_e''(0)/2!=1/K_2

于是

Φe(s)=1K0+1K1s+1K2s2+\Phi_e(s)=\frac{1}{K_0}+\frac{1}{K_1}s+\frac{1}{K_2}s^2+\cdots ess(t)=1K0r(t)+1K1r(t)+1K2r(t)+e_{ss}(t)=\frac{1}{K_0}r(t)+\frac{1}{K_1}r'(t)+\frac{1}{K_2}r''(t)+\cdots

胡寿松把展开系数直接叫 C0,C1,C2C_0,C_1,C_2(广义误差系数):C0=Φe(0)C_0=\Phi_e(0)C1=Φe(0)C_1=\Phi_e'(0)C2=Φe(0)/2!C_2=\Phi_e''(0)/2!。与徐薇莉的关系是 Ki=1/CiK_i=1/C_i。同一套级数,记号不同。

与静态系数的常见对应(单位反馈):0 型 K0=1+KpK_0=1+K_p;I 型 K0=K_0=\inftyK1=KvK_1=K_v。静态系数是动态系数在典型输入、tt\to\infty 时的特例。

3. 两个求法(徐薇莉)

方法一:定义直接算 Φe(0),Φe(0),Φe(0)\Phi_e(0),\Phi_e'(0),\Phi_e''(0)。只适合 Φe(s)\Phi_e(s) 很简单。

方法二(常用):把 Φe(s)\Phi_e(s) 的分子除以分母(在 s=0s=0 邻域作升幂除法),商式与

1K0+1K1s+1K2s2+\frac{1}{K_0}+\frac{1}{K_1}s+\frac{1}{K_2}s^2+\cdots

比较同次项,读出 K0,K1,K2K_0,K_1,K_2。再把 r,r,rr,r',r'' 代入误差级数。

4. 注意

  • 级数描述的是 tt\to\infty 附近的误差变化,不是过渡过程超调
  • 输入含随 tt\to\infty 增长的项(斜坡、抛物线等)时,用级数比只报一个 esse_{ss} 更清楚
  • 两个开环 Kp,Kv,KaK_p,K_v,K_a 相同的系统,静态误差可以一样,动态误差级数仍可不同(徐薇莉例 4-14 / 4-16)
  • 仍须系统稳定,泰勒展开才对应稳态项

5. 做题顺序

  1. 写出 Φe(s)=1/(1+GH)\Phi_e(s)=1/(1+GH)
  2. s=0s=0 展开或升幂相除,得 1/K01/K_01/K11/K_11/K21/K_2
  3. 代入 r(t)r(t) 及其导数,写出 ess(t)e_{ss}(t)
  4. 若只要终值,再 limtess(t)\lim_{t\to\infty}e_{ss}(t)

一句:把 Φe(s)\Phi_e(s)s=0s=0 展开,误差稳态项 = 各阶系数 ×\times 输入各阶导数;K0,K1,K2K_0,K_1,K_2 是这些系数的倒数。

动态参数

出处:徐薇莉 §4.2–4.4;胡寿松 §3-1–3-3。动态性能看单位阶跃 r(t)=1(t)r(t)=1(t) 下的 h(t)h(t);稳态性能看 tt\to\infty 的误差(见上节)。指标默认针对零初始、单位负反馈、无零点的标准型。

1. 阶跃响应上怎么读

单位阶跃响应动态指标
单位阶跃响应动态指标
指标定义反映什么
延迟时间 tdt_dh(t)h(t) 第一次到达 0.5h()0.5\,h(\infty) 的时刻反应快慢
上升时间 trt_r一阶 / 过阻尼:从 10%h()10\%h(\infty)90%h()90\%h(\infty)欠阻尼二阶:从 00 第一次到 h()h(\infty)(即 tr=(πβ)/ωdt_r=(\pi-\beta)/\omega_d响应速度
峰值时间 tpt_p超过 h()h(\infty) 后到达第一个峰值的时刻;无超调则没有 tpt_p
调节时间 tst_s进入并保持在 h()±Δh(\infty)\pm\Delta 误差带内的最短时间;Δ=5%\Delta=5\%2%2\%快速性(进入稳态要多久)
超调量 σ%\sigma\%σ%=h(tp)h()h()×100%\displaystyle\sigma\%=\frac{h(t_p)-h(\infty)}{h(\infty)}\times 100\%阻尼程度
  • 曲线从 00 上升,穿过 0.10.10.50.50.90.9 倍稳态值,冲过 h()h(\infty) 到峰,再衰减进误差带
  • 稳态误差是 tt\to\infty 时输出与希望值之差,不是动态指标

2. 一阶系统

微分方程与传函(时间常数 T>0T>0):

Tc˙+c=r,Φ(s)=C(s)R(s)=1Ts+1T\dot{c}+c=r,\qquad \Phi(s)=\frac{C(s)}{R(s)}=\frac{1}{Ts+1}

单位阶跃响应(单调上升,无超调、无 tpt_p):

h(t)=1et/T,t0h(t)=1-e^{-t/T},\quad t\ge 0
指标公式由来
tdt_d0.69T0.69Th=0.5t=Tln2h=0.5\Rightarrow t=T\ln 2
trt_r2.20T2.20T10%90%10\%\to 90\%
tst_s3T3TΔ=5%\Delta=5\%h=0.95h=0.95Δ=2%\Delta=2\% 时约 4T4T

TT 越大越慢。TT 就是闭环极点到原点距离的倒数。

3. 二阶无零点标准型

Φ(s)=C(s)R(s)=ωn2s2+2ζωns+ωn2\Phi(s)=\frac{C(s)}{R(s)}=\frac{\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}

单位负反馈时开环为

G(s)=ωn2s(s+2ζωn),H(s)=1G(s)=\frac{\omega_n^2}{s(s+2\zeta\omega_n)},\qquad H(s)=1

特征方程 s2+2ζωns+ωn2=0s^2+2\zeta\omega_n s+\omega_n^2=0,根

s1,2=ζωn±ωnζ21s_{1,2}=-\zeta\omega_n\pm\omega_n\sqrt{\zeta^2-1}
  • ωn\omega_n无阻尼自然振荡频率(极点到原点的距离)
  • ζ\zeta阻尼比。也可写成 ζ=f/fc\zeta=f/f_c(实际阻尼 / 临界阻尼)
  • 机电对象常见写法(不是定义本身):ωn=K/Tm\omega_n=\sqrt{K/T_m}ζ=1/(2TmK)\zeta=1\big/\bigl(2\sqrt{T_m K}\bigr)

分子必须等于分母常数项 ωn2\omega_n^2,且 Φ(0)=1\Phi(0)=1(无零点、单位阶跃稳态值为 1)。有零点或增益不是 ωn2\omega_n^2 时,下面 tr,tp,σ%t_r,t_p,\sigma\% 公式不能直接套

ζ\zeta名称阶跃表现
<0<0不稳定
=0=0无阻尼等幅振荡
0<ζ<10<\zeta<1欠阻尼超调、衰减振荡(最常考)
=1=1临界阻尼无超调,上升最快(在无超调一类里)
>1>1过阻尼无超调,上升更慢

定性:欠阻尼才超调振荡;临界无超调且最快;过阻尼无超调但慢。

4. 欠阻尼(0<ζ<10<\zeta<1,常考)

衰减系数、阻尼振荡频率、极点:

σ=ζωn,ωd=ωn1ζ2,s1,2=σ±jωd\sigma=\zeta\omega_n,\qquad \omega_d=\omega_n\sqrt{1-\zeta^2},\qquad s_{1,2}=-\sigma\pm j\omega_d

ss 平面:原点到极点为 ωn\omega_n;到虚轴水平距离为 σ\sigma;到实轴竖直距离为 ωd\omega_d。负实轴到极点的夹角 β\beta 满足

cosβ=ζ,β=arccosζ\cos\beta=\zeta,\qquad \beta=\arccos\zeta

单位阶跃:

h(t)=1eσt1ζ2sin(ωdt+β),t0h(t)=1-\frac{e^{-\sigma t}}{\sqrt{1-\zeta^2}}\sin(\omega_d t+\beta),\quad t\ge 0

单位脉冲(胡寿松记 k(t)k(t);零初始时 k(t)=dh/dtk(t)=\mathrm{d}h/\mathrm{d}t,即 Φ(s)\Phi(s) 的反拉氏):

k(t)=ωn1ζ2eσtsin(ωdt),t0k(t)=\frac{\omega_n}{\sqrt{1-\zeta^2}}e^{-\sigma t}\sin(\omega_d t),\quad t\ge 0

无阻尼 ζ=0\zeta=0 时退化为 k(t)=ωnsin(ωnt)k(t)=\omega_n\sin(\omega_n t)。临界:k(t)=ωn2teωntk(t)=\omega_n^2 t\,e^{-\omega_n t}

指标公式怎么变快 / 变小
tdt_d1+0.6ζ+0.2ζ2ωn1+0.7ζωn\dfrac{1+0.6\zeta+0.2\zeta^2}{\omega_n}\approx\dfrac{1+0.7\zeta}{\omega_n}ζ\zeta\downarrowωn\omega_n\uparrow td\Rightarrow t_d\downarrow(极点离原点或离虚轴更远)
trt_rπβωd\dfrac{\pi-\beta}{\omega_d}ζ\zeta\downarrowωd\omega_d\uparrow tr\Rightarrow t_r\downarrow
tpt_pπωd\dfrac{\pi}{\omega_d}阻尼振荡周期的一半
σ%\sigma\%eπζ/1ζ2×100%e^{-\pi\zeta/\sqrt{1-\zeta^2}}\times 100\%只与 ζ\zeta 有关ζσ%\zeta\uparrow\Rightarrow\sigma\%\downarrow,与 ωn\omega_n 无关
tst_sΔ=5%\Delta=5\%(默认):3.5ζωn\dfrac{3.5}{\zeta\omega_n}Δ=2%\Delta=2\%4.4ζωn\dfrac{4.4}{\zeta\omega_n}由衰减系数 σ=ζωn\sigma=\zeta\omega_n 决定

由超调量反求阻尼比(σ\sigma小数,如 16.3%σ=0.16316.3\%\Rightarrow\sigma=0.163,不要代入 16.3):

ζ=lnσπ2+(lnσ)2\zeta=\frac{\lvert\ln\sigma\rvert}{\sqrt{\pi^2+(\ln\sigma)^2}}

5. 注意

  • 先确认是不是无零点标准型;有零点会抬超调、改 tpt_p
  • trt_r 两种定义不要混:一阶用 10%90%10\%\sim 90\%;欠阻尼二阶用 0100%0\to 100\%
  • tst_s 必须写明误差带;胡寿松常用 3.5/σ3.5/\sigma(5%)和 4.4/σ4.4/\sigma(2%),有的书用 3344
  • 要减小超调:加大 ζ\zeta(或减小开环 KK,若 KKζ\zeta);要加快 tst_s:加大 ζωn\zeta\omega_n(极点左移)
  • ζ=0\zeta=0 等幅振荡谈不上 tst_sζ1\zeta\ge 1 无超调,不用 σ%\sigma\%tpt_p 那套公式

6. 做题顺序

  1. 化成 ωn2/(s2+2ζωns+ωn2)\omega_n^2\big/\bigl(s^2+2\zeta\omega_n s+\omega_n^2\bigr),读出 ζ,ωn\zeta,\omega_n
  2. ζ\zeta 判断欠 / 临界 / 过阻尼
  3. 欠阻尼:算 ωd\omega_dβ\beta,再套 tr,tp,σ%,tst_r,t_p,\sigma\%,t_s
  4. σ%\sigma\%ζ\zeta 时,超调先化成小数再代入反公式
  5. 一阶只记 td=0.69Tt_d=0.69Ttr=2.2Tt_r=2.2Tts=3Tt_s=3T(5%)

一句:动态指标全从单位阶跃曲线上读;一阶只由 TT 决定;二阶标准型由 (ζ,ωn)(\zeta,\omega_n) 决定——超调只看 ζ\zeta,调节时间看 ζωn\zeta\omega_n,振荡快慢看 ωd\omega_d

二阶系统性能改善

出处:胡寿松 §3-3(比例-微分与测速反馈);对象默认无零点标准型 G0(s)=ωn2/[s(s+2ζωn)]G_0(s)=\omega_n^2\big/\bigl[s(s+2\zeta\omega_n)\bigr],单位负反馈。816 不点名「测速」二字,但内环 KtsK_t s、先令 Kt=0K_t=0 再比 ζ\zeta 与斜坡 esse_{ss} 是固定套路(2009、2010、2013、2015、2020)。

原系统开环增益(速度误差系数)Kv=ωn/(2ζ)K_v=\omega_n/(2\zeta)。只调这个 KK 时,加大 KK \Rightarrow ζ\zeta 变小、超调变大,斜坡 esse_{ss} 变小,动态和稳态拧着,所以要另加阻尼通道。

1. 比例-微分(PD)

误差经 1+Tds1+T_d s 再进对象:不改 ωn\omega_n,加大阻尼

比例微分结构图
比例微分结构图
ζPD=ζ+12Tdωn\zeta_{\mathrm{PD}}=\zeta+\frac{1}{2}T_d\omega_n

闭环出现零点 s=1/Tds=-1/T_d,超调会比「同样 ζ\zeta、无零点」略大。开环 KvK_v 不变,斜坡 esse_{ss} 不因加微分而变差。

2. 测速反馈

输出是位置时,对其求导得速度,经 KtsK_t s内负反馈(外环仍为单位反馈)。不改 ωn\omega_n,加大阻尼,不引入闭环零点

测速反馈结构图
测速反馈结构图
Φ(s)=ωn2s2+2ζtωns+ωn2,ζt=ζ+12Ktωn,Kt=2(ζtζ)ωn\Phi(s)=\frac{\omega_n^2}{s^2+2\zeta_t\omega_n s+\omega_n^2},\qquad \zeta_t=\zeta+\frac{1}{2}K_t\omega_n,\qquad K_t=\frac{2(\zeta_t-\zeta)}{\omega_n}

化成时间常数型后,等效开环增益变小:

G(s)=Ks(Ts+1),K=ωn2ζ+Ktωn=ωn22ζωn+Ktωn2,T=12ζωn+Ktωn2G(s)=\frac{K}{s(Ts+1)},\qquad K=\frac{\omega_n}{2\zeta+K_t\omega_n}=\frac{\omega_n^2}{2\zeta\omega_n+K_t\omega_n^2},\qquad T=\frac{1}{2\zeta\omega_n+K_t\omega_n^2}

ζt,ωn\zeta_t,\omega_n 确定后,按欠阻尼公式算 σ%,tr,tp,ts\sigma\%,t_r,t_p,t_s

测速反馈压 σ%\sigma\%、缩短 tst_s,但 KvK_v 变小 \Rightarrow 斜坡 esse_{ss} 变大。同时加大前向增益,把 KvK_v 补回来,就可以既加阻尼又不恶化 esse_{ss}

PD测速反馈
ωn\omega_n不变不变
ζ\zeta增大增大
闭环零点
斜坡 KvK_v / esse_{ss}不变KvK_v 下降,esse_{ss} 变大(除非再加大前向 KK
输入端噪声微分会放大内反馈对输入噪声有滤波

3. 做题顺序

  1. 认出内环 KtsK_t s(或积分前引出经 k2k_2 反馈,等效测速)
  2. 化简内环,读出新的 ζt\zeta_tωn\omega_nωn\omega_n 一般不变)
  3. ζt\zeta_t 算动态指标;用化简后的开环算 KvK_vesse_{ss}
  4. 常考对比:先令 Kt=0K_t=0,再打开 KtK_t,看 ζ\zeta 变大、斜坡误差变大
  5. 给定 ζ\zetaesse_{ss} 联立:前向 KKKtK_t 一起解

一句:PD 和测速都能加大 ζ\zeta、不改 ωn\omega_n;测速不引进零点但会降低 KvK_v,要保 esse_{ss} 就再把前向增益抬上去。