自动控制原理入门教程——从反馈到稳定性

概述

自动控制原理是研究自动控制系统一般规律的学科，是电气工程、机械工程、航空航天等领域的核心基础。本教程从控制系统的基本概念出发，依次介绍传递函数与方框图、时域分析、根轨迹法、频域分析以及系统校正方法，帮助读者建立完整的控制理论分析框架。

知识点一：控制系统的基本概念与传递函数

开环控制与闭环控制

• 开环控制：控制器输出不受系统实际输出影响。结构简单，但无法自动修正误差。
• 闭环控制（反馈控制）：将系统输出反馈到输入端，与期望值比较产生偏差信号，再驱动控制器。闭环控制能自动纠偏，抗干扰能力强。

传递函数

传递函数是线性定常系统在零初始条件下，输出的拉普拉斯变换与输入的拉普拉斯变换之比：

\(G(s) = \frac{Y(s)}{X(s)}\)

传递函数完全描述了系统的输入-输出特性，是经典控制理论的核心工具。

方框图与信号流图

复杂控制系统可通过方框图表示各环节的连接关系。基本连接方式：

• 串联：\(G(s) = G_1(s) \cdot G_2(s)\)
• 并联：\(G(s) = G_1(s) + G_2(s)\)
• 反馈连接：闭环传递函数 \(\frac{Y(s)}{R(s)} = \frac{G(s)}{1 + G(s)H(s)}\)

其中 \(G(s)\) 为前向通道传递函数，\(H(s)\) 为反馈通道传递函数。\(1 + G(s)H(s) = 0\) 称为系统的特征方程。

例子：系统前向通道 \(G(s) = \frac{K}{s(s+1)}\)，单位反馈 \(H(s)=1\)，求闭环传递函数。

\(\Phi(s) = \frac{G(s)}{1+G(s)} = \frac{K}{s^2 + s + K}\)

知识点二：时域分析与性能指标

典型输入信号

• 阶跃信号 \(r(t) = R \cdot u(t)\)：最常用，考察系统对突变输入的响应
• 斜坡信号 \(r(t) = Rt\)：考察跟踪能力
• 冲激信号 \(r(t) = R\delta(t)\)：考察系统动态特性

二阶系统的时域响应

典型二阶系统传递函数： \(\Phi(s) = \frac{\omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}\)

其中 \(\zeta\) 为阻尼比，\(\omega_n\) 为自然频率。

根据 \(\zeta\) 的不同，系统响应分为：

• \(\zeta = 0\)：无阻尼，等幅振荡
• \(0 < \zeta < 1\)：欠阻尼，衰减振荡（最常见）
• \(\zeta = 1\)：临界阻尼
• \(\zeta > 1\)：过阻尼，缓慢趋近

性能指标

欠阻尼二阶系统的阶跃响应性能指标：

• 上升时间 \(t_r = \frac{\pi - \beta}{\omega_d}\)，其中 \(\omega_d = \omega_n\sqrt{1-\zeta^2}\)
• 峰值时间 \(t_p = \frac{\pi}{\omega_d}\)
• 超调量 \(\sigma\% = e^{-\frac{\zeta\pi}{\sqrt{1-\zeta^2}}} \times 100\%\)
• 调节时间 \(t_s \approx \frac{3}{\zeta\omega_n}\)（取5%误差带）

例子：二阶系统 \(\Phi(s) = \frac{25}{s^2 + 4s + 25}\)，求 \(\zeta\), \(\omega_n\), 超调量。

对比标准形式：\(\omega_n^2 = 25 \Rightarrow \omega_n = 5\)，\(2\zeta\omega_n = 4 \Rightarrow \zeta = 0.4\)。

超调量 \(\sigma\% = e^{-0.4\pi/\sqrt{1-0.16}} \times 100\% = e^{-0.4\pi/0.917} \times 100\% \approx 25.4\%\)

知识点三：根轨迹法

根轨迹的概念

根轨迹是当系统开环增益 \(K\) 从 \(0\) 变化到 \(+\infty\) 时，闭环极点在 \(s\) 平面上移动的轨迹。通过根轨迹可以直观地分析增益变化对系统稳定性和动态性能的影响。

绘制规则

1. 起点与终点：根轨迹起于开环极点（\(K=0\)），终于开环零点（\(K\to\infty\)）
2. 分支数：等于开环极点数
3. 对称性：根轨迹关于实轴对称
4. 实轴上的根轨迹：实轴上某点右侧的开环零、极点个数之和为奇数，则该点在根轨迹上
5. 渐近线：与实轴交角 \(\theta_k = \frac{(2k+1)\pi}{n-m}\)，交点 \(\sigma_a = \frac{\sum p_i - \sum z_j}{n-m}\)
6. 分离点：满足 \(\sum \frac{1}{d-p_i} = \sum \frac{1}{d-z_j}\)

例子：开环传递函数 \(G(s) = \frac{K}{s(s+2)(s+4)}\)，绘制根轨迹。

• 开环极点：\(0, -2, -4\)（3个），无零点
• 渐近线：角度 \(60°, 180°, 300°\)；交点 \(\sigma_a = \frac{0-2-4}{3} = -2\)
• 分离点：\(\frac{1}{d} + \frac{1}{d+2} + \frac{1}{d+4} = 0\)，解得 \(d \approx -0.85\)
• 与虚轴交点：令 \(s = j\omega\) 代入特征方程，可求得临界增益 \(K = 48\)

知识点四：频域分析与稳定性判据

频率特性

将传递函数中的 \(s\) 替换为 \(j\omega\)，得到系统的频率特性： \(G(j\omega) = |G(j\omega)| e^{j\angle G(j\omega)}\)

• 幅频特性 \(|G(j\omega)|\)：系统对不同频率信号的放大/衰减
• 相频特性 \(\angle G(j\omega)\)：系统引入的相位滞后

Bode图

Bode图由幅频图（\(20\lg|G|\) vs \(\lg\omega\)）和相频图（\(\angle G\) vs \(\lg\omega\)）组成。

绘制要点：

• 每个零点在转折频率处使斜率增加 \(+20\) dB/dec
• 每个极点在转折频率处使斜率增加 \(-20\) dB/dec
• 一阶因子 \(\frac{1}{1+j\omega/\omega_c}\) 在 \(\omega_c\) 处相位变化 \(-90°\)

Nyquist稳定判据

Nyquist判据：闭环系统稳定的充要条件是，\(G(j\omega)H(j\omega)\) 的Nyquist曲线逆时针包围 \((-1, j0)\) 点的圈数等于开环传递函数在右半平面的极点数。

简化判据（开环稳定时）：Nyquist曲线不包围 \((-1, j0)\) 点，则闭环稳定。

稳定裕度

• 增益裕度 \(K_g\)：在相位穿越频率处，增益距0dB的距离
• 相位裕度 \(\gamma\)：在增益穿越频率处，相位距 \(-180°\) 的距离

一般要求 \(\gamma > 30°\sim60°\)，\(K_g > 6\) dB。

例子：\(G(s) = \frac{10}{s(s+1)(s+5)}\)，求相位裕度。

增益穿越频率：\(|G(j\omega_c)| = 1\)，近似求解得 \(\omega_c \approx 2.4\) rad/s。

相位：\(\angle G(j\omega_c) = -90° - \arctan(2.4) - \arctan(0.48) \approx -90° - 67.4° - 25.6° = -183°\)

相位裕度 \(\gamma = 180° - 183° = -3° < 0\)，系统不稳定。

知识点五：系统校正方法

PID控制器

PID控制器是最广泛使用的工业控制器：

\(u(t) = K_p e(t) + K_i \int_0^t e(\tau)d\tau + K_d \frac{de(t)}{dt}\)

• P（比例）：快速响应，减小偏差，但不能消除稳态误差
• I（积分）：消除稳态误差，但增大超调、减慢响应
• D（微分）：预测趋势，改善动态性能，减小超调

超前校正与滞后校正

• 超前校正：提供正相位，增大相位裕度，改善动态性能。适用于系统响应慢、需要增加稳定裕度的情况。
• 滞后校正：利用高频衰减提高低频增益，在不改变动态性能的前提下减小稳态误差。

校正设计的一般步骤：

1. 根据稳态误差要求确定开环增益
2. 绘制校正前的Bode图，确定相位裕度
3. 根据期望裕度确定校正网络参数
4. 验证校正后系统性能

例子：对 \(G(s) = \frac{4}{s(s+2)}\) 设计超前校正，使相位裕度 \(\geq 45°\)。

未校正系统在 \(\omega_c \approx 2\) rad/s 处的相位裕度约 \(18°\)。需要增加约 \(30°\) 相位。选择超前校正网络 \(G_c(s) = \frac{1+aTs}{1+Ts}\)，取 \(a = 3\)，在 \(\omega_m = 2\) rad/s 处提供最大超前角 \(\sin^{-1}\frac{a-1}{a+1} = 30°\)，可得 \(T = 0.29\)，校正网络为 \(G_c(s) = \frac{1+0.87s}{1+0.29s}\)。

练习题

题目1

单位反馈系统的开环传递函数为 \(G(s) = \frac{5}{s(s+1)}\)，求闭环传递函数。

答案： \(\Phi(s) = \frac{G(s)}{1+G(s)} = \frac{5}{s^2 + s + 5}\)

题目2

二阶系统 \(\Phi(s) = \frac{9}{s^2 + 2.4s + 9}\)，求阻尼比和自然频率。

答案：\(\omega_n = 3\) rad/s，\(2\zeta\omega_n = 2.4\)，\(\zeta = 0.4\)。系统为欠阻尼系统。

题目3

开环传递函数 \(G(s) = \frac{K}{s(s+3)}\)，求根轨迹的分离点。

答案：特征方程 \(1 + \frac{K}{s(s+3)} = 0 \Rightarrow K = -s(s+3) = -s^2 - 3s\)。令 \(\frac{dK}{ds} = -2s - 3 = 0\)，得 \(s = -1.5\)。分离点为 \(s = -1.5\)。

题目4

什么是增益裕度和相位裕度？为什么工程中通常要求相位裕度在 \(30°\sim60°\) 之间？

答案：增益裕度是Nyquist曲线穿越负实轴时，增益还能增加多少倍才到 \(-1\) 点。相位裕度是增益为1时相位距 \(-180°\) 的余量。\(30°\sim60°\) 的相位裕度折中了响应速度和超调量：太小则超调大、振荡剧烈；太大则响应迟缓。

题目5

简述PID控制器中P、I、D三个分量各自的作用。

答案：P（比例）快速减小偏差但有稳态误差；I（积分）消除稳态误差但增大超调；D（微分）预测误差变化趋势，改善动态性能、减小超调。

总结

自动控制原理的核心框架可以概括为：

1. 数学建模：用传递函数和方框图描述系统
2. 时域分析：通过阶跃响应的性能指标评价系统品质
3. 根轨迹法：图形化分析增益对闭环极点的影响
4. 频域分析：用Bode图和Nyquist判据分析稳定性，用稳定裕度量化系统鲁棒性
5. 系统校正：通过PID或超前/滞后校正改善系统性能

这些分析方法相辅相成，是控制系统设计与工程应用的理论基础。掌握这些工具，将为学习现代控制理论、数字控制等进阶课程奠定坚实基础。