📐 控制数学工具箱:拉普拉斯、传递函数与 Bode 图

把"微分方程"翻译成"代数运算"的三件套——拖动 ζ 与 ωn,实时看阶跃响应与 Bode 图如何联动,读懂带宽/超调/谐振峰。

拉普拉斯变换传递函数二阶系统Bode图

🎯 本页学习目标

1. 说清拉普拉斯变换解决了什么问题(微分→代数、时域→s 域)
2. 掌握二阶系统标准式,理解 ζ(阻尼)与 ωn(自然频率)如何同时决定超调与速度
3. 会读 Bode 图:−3dB 带宽、谐振峰值 Mr、相位滞后,并能联系到三环整定的带宽法
建议用时:约 30 分钟

1 拉普拉斯变换:控制工程师的"翻译器"

线性微分方程直接求解要背齐次解+特解,叠加初始条件,繁琐易错。拉普拉斯变换

F(s) = \mathcal{L}[f(t)] = \int_0^\infty f(t)\,e^{-st}\,dt

把微分变成乘 s、积分变成除 s、卷积变成乘法,整个微分方程变成关于 s 的一元代数方程:解出 F(s),再查表反变换回时域。三个最常用的对:

时域 f(t)	s 域 F(s)	控制系统里是谁
单位阶跃 1(t)	1/s	阶跃指令/突加负载
e^−at	1/(s+a)	一阶衰减(极点在 −a)
sin ωt	ω/(s²+ω²)	振荡/正弦激励
f′(t)	sF(s) − f(0)	微分(初值为 0 时就是 ×s)
∫f dt	F(s)/s	积分(速度→位置)

💡 一句话本质:s 平面上的极点位置 = 系统的"性格"——负实轴上的极点=衰减,离虚轴越远衰减越快;带虚部分量=振荡;虚轴上=等幅振荡(临界不稳);右半平面=发散。整定 PID 本质上就是在搬极点的位置。

2 一阶环节:时间常数 τ 的直觉

惯性环节 G(s)=1/(τs+1):阶跃响应以指数爬升,每过一个 τ 走完剩余距离的 63.2%,3τ 到 95%,5τ 到 99%。RL 电路(τ=L/R)、电机机械加速(τ=J/b)、温度对象都是一阶。FOC 电流环带宽法整定(Kp=ωcL, Ki=ωcR)做完极点对消后,闭环就化成这个形状——"带宽 1/(2π) Hz"和"时间常数"说的同一件事。

3 二阶环节交互实验(拖滑杆看全部)

标准二阶系统是三环里"速度环+位置环"级联后的典型形状:

G(s) = \frac{\omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}, \qquad M_p = e^{-\pi\zeta/\sqrt{1-\zeta^2}} \times 100\% \;(\zeta<1)

阻尼比 ζ: 0.35 自然频率 ωn (rad/s): 6.0

💡 怎么读上图:左图阶跃响应(绿线)——ζ↓超调↑振荡↑,ωn↓整体变慢;虚线包络 e^{−ζωn t} 是"它最多能荡到哪"。右图 Bode 幅频(对数横轴)——低频 0dB,过了 ωn 以 −40dB/十倍频下坠;ζ 小时 ωn 附近鼓起的"谐振峰"就是机械系统共振的数学影子。

4 Bode 图:频域体检报告

给系统灌不同频率的正弦,稳态输出的幅值比与相位差随频率的曲线就是 Bode 图。三个必读指标:

−3dB 带宽 ωb:幅值跌到 0.707 的频率——"跟得动"的指令频率上限,高于它的指令被低通掉(这正是三环"内环带宽 ≥ 5~10 倍外环"的物理含义)。
谐振峰值 M_r = 1/(2ζ√(1−ζ²)):ζ<0.707 时出现;Mr>2(ζ<0.26)意味着阶跃超调>40%,工程上先把 ζ 提上去再谈带宽。
相位滞后:每个极点贡献最多 −90°;级联三个环节累计 −270°,若在增益>1 的频率上相位到 −180°,闭环自激振荡(奈奎斯特判据的直观版)。

⚠️ 与整定的连接:调参手册的带宽法就是在放置闭环极点:电流环极点对消后是一阶(τ=1/ωc);速度/位置环整定后的闭环近似二阶——你现在拖动的 ζ 和 ωn,就是整定时"超调 5% 对应 ζ≈0.7、响应快慢对应 ωn"的那两个旋钮。数学与工程在此合流。

5 从开环到闭环:闭环传函与带宽

前四节看的都是"单个环节",真实的关节控制要把控制器、被控对象、传感器接成负反馈回路才工作。设前向通道(控制器×对象)传函为 G(s),反馈通道(传感器)为 H(s):控制器拿到的是误差——指令减去量测,即 E = R − HC;输出由前向通道产生,即 C = GE。消去 E 就得到闭环传函:

T(s) = \frac{C(s)}{R(s)} = \frac{G(s)}{1 + G(s)\,H(s)}, \qquad \frac{E(s)}{R(s)} = \frac{1}{1 + G(s)\,H(s)}

为什么压误差:误差对指令的传函是 1/(1+GH)——回路增益 |GH| 越大,同样 1° 的偏差被压得越小;低频 |GH| 远大于 1 时 T≈1/H,闭环特性几乎只由传感器决定,控制器与对象自身的参数漂移也被同一因子压制。负反馈"压误差、抗扰动、降参数敏感"共用这一个分母:1+GH。
为什么抗扰动:负载扰动 D 从对象入口混入时,对输出的影响约为开环时的 1/(1+GH)——开环系统对扰动照单全收,闭环把它压掉 (1+GH) 倍。
比喻(空调恒温):设定温度是指令 R,室温是输出 C,空调按"差几度"(误差)持续修正:温差大猛吹、温差小细调,把室温按在设定值附近。GH 就是"修正力度"——力度越大残余温差越小,但修正过急会让室温来回过冲,这正是"带宽想高"与"裕度要留"的矛盾。

闭环带宽:闭环本身也是个低通滤波器——低频指令幅值≈1(跟得上),高频指令既被衰减又被滞后(跟不上)。工程定义:闭环幅值跌到低频值的 0.707(−3dB)处的频率为闭环带宽 ω_b,即"跟得上输入的最高频率"。它像耳朵能听清的音域范围(人耳约 20Hz~20kHz):音域之外的声音要么听不见、要么糊成一片;带宽 10Hz 的关节面对 15Hz 的手抖指令就是"听不清"——幅度跟不上,相位还迟到。呼应第 2 节:对一阶闭环,ω_b=1/τ——"时间常数"与"带宽"本是同一件事的时域、频域两种说法。

算例:一阶对象 + 比例控制的闭环带宽(给定条件→逐步代入→判定)。已知被控对象为惯性环节 G₀(s)=K/(Ts+1),K=0.8,T=0.02s;比例控制器 K_p=2.5,单位负反馈 H=1;要求闭环带宽 ≥10Hz。

第 1 步,写出开环传函,回路增益 K_pK=2.5×0.8=2:

G(s) = K_p\,G_0(s) = \frac{2.5\times 0.8}{0.02s+1} = \frac{2}{0.02s+1}

第 2 步,代入 T=G/(1+G),分子分母同除 (1+K_pK)=3 化成标准一阶:等效时间常数由 20ms 压到 6.67ms,负反馈把系统提速 3 倍:

T(s) = \frac{G(s)}{1+G(s)} = \frac{2}{0.02s+3} = \frac{0.667}{0.00667\,s+1}

第 3 步,求带宽:一阶闭环的 −3dB 点恰在 1/τ 处(此处幅值比 0.707、相位滞后 45°):

\omega_b = \frac{1}{\tau_{cl}} = \frac{1+K_pK}{T} = \frac{3}{0.02} = 150\ \mathrm{rad/s} \approx 23.9\ \mathrm{Hz}

第 4 步,判定:要求 ≥10Hz,实际 ≈23.9Hz——✓ 达标,余量 (23.9−10)/10 ≈ 139%;对应阶跃响应 3τ≈20ms 基本到位、5τ≈33ms 完全稳定。

💡 算例附注:①纯 P 控制的代价同框可见:闭环直流增益 K_pK/(1+K_pK)=0.667,阶跃稳态误差 1/(1+K_pK)≈33%,真实电流环要再加积分项消差;②若把 K_p 翻倍到 5:ω_b=(1+4)/0.02=250rad/s≈39.8Hz,余量更大,但稳态误差只降到 20%,且回路增益越高、穿越点越滑向相位滞后更深的频段,稳定裕度越薄——带宽不能无限抬(下一节);③带宽只是"跟得上"的必要条件,充分性还差相频那一半。数值为教学示例,实际以电机手册与实测整定为准。

6 幅频与相频各看什么

第 3 节实验台的右图只画了幅频,但 Bode 图是两张图:幅频回答"幅度还剩多少",相频回答"时间晚了多少",合起来才是完整体检报告。只看幅频会漏掉稳定性隐患:当某个频率上相位逼近 −180°,而回路增益在该频率仍大于 1(0dB 线以上),负反馈就在这一频段变成正反馈——误差不被修正反被放大,闭环自激振荡。两个系统的幅频曲线可以几乎一样"平缓健康",一个稳如泰山、一个嗡嗡自振,差别全在相频。带宽法整定的完整口诀:幅频定快慢,相频定稳不稳。

幅频:幅度还剩多少(精度·带宽)	相频:时间晚了多少(滞后·稳定)
低频段 0dB 附近:直流增益≈1,输出忠实复现指令——决定稳态精度	低频段 ≈0°:输出与指令同步,几乎无延迟
−3dB 点 = 闭环带宽:跟得动的指令频率上限	带宽处相位通常已滞后 45°~90°:带宽越高,动态延迟越小
高频段衰减越快:对测噪、齿槽转矩、振动的抑制越强	相位逼近 −180° 的频段,增益必须已衰减到 1 以下,否则自激振荡
谐振峰 M_r:阻尼不足、机械共振的前兆	相位裕度 PM=180°+φ(ω_c):整定留下的"安全垫",工程常取 30°~60°

💡 对照实验台:回到第 3 节把 ζ 调到 0.35:右图在 ω_n 附近鼓起谐振峰,而此刻(实验台没画的)相频在该处恰为 −90°、高频段一路逼近 −180°——幅频上那只"鼓包"能不能留,要靠相频与增益的配合(裕度)裁决。读图请把"幅频看精度、相频看稳定"练成肌肉记忆。

7 本节自测

1. 拉普拉斯变换给控制系统带来的核心便利是?

💡 解析:拉氏变换是"线性"工具(叠加原理成立),它不改变系统性质、不除噪、更不保证稳定——它只是把时域微积分翻译成 s 域代数,让求解与频域分析(Bode/根轨迹)成为可能。

2. 二阶系统 ζ=0.2,保持 ωn 不变把 ζ 提到 0.7,阶跃响应会?

💡 解析:Mp=e^{−πζ/√(1−ζ²)}:ζ=0.2 超调约 52%,ζ=0.7 约 4.6%。阻尼增加以轻微牺牲上升速度为代价(阻尼比 0.707 是"最快无超调"的工程甜点);稳态误差由系统型别与积分器决定,与 ζ 无关。

3. "内环带宽必须 ≥ 外环 5~10 倍"的频域解释是?

💡 解析:外环发出的指令频率 ≤ 外环带宽;若内环带宽不足,在这个频段内内环幅值衰减+相位滞后,级联后回路额外损失相位裕度 → 振荡。5~10 倍保证内环在外环工作频段的相位滞后 <~10°,这是级联控制的频域设计准则。

4. 二阶系统 Bode 图上「−3dB 带宽 ωb」处,输出幅值比约为?

💡 解析:−3dB 由 20·lg(幅值比)=−3 解得幅值比≈0.707(1/√2),对应第 4 节「幅值跌到 0.707 的频率」。常见误区:把 −3dB 当成「衰减一半」(0.5)——那是 −6dB。

5. 关于一阶惯性环节的时间常数 τ,下列说法错误的是?

💡 解析:惯性环节是指数渐近,永远「接近但不恰好到达」稳态:1τ 只到 63.2%,5τ 才≈99%(见第 2 节)。「过 1 个 τ 就到稳态」是典型错误认知。

6. 二阶系统 Bode 幅频曲线出现「谐振峰 M_r」的条件是?

💡 解析:谐振峰 M_r=1/(2ζ√(1−ζ²)),只有当 ζ<0.707 时才有峰值(根号内 1−ζ²>0 且 M_r>1);ζ≥0.707 幅频单调下坠、无峰(见第 4 节)。

📌 本节要点速查:①拉氏变换把微积分翻译成 s 域代数(微分×s、积分÷s、卷积×),极点位置=系统性格;②一阶看时间常数 τ:1τ 走完剩余 63.2%、3τ≈95%、5τ≈99%;③二阶由 ζ(管超调)与 ωn(管快慢)二维定性,Mp=e^{−πζ/√(1−ζ²)}、ζ=0.707 为无超调甜点;④Bode 三指标:−3dB 带宽(幅值 0.707)、谐振峰 Mr=1/(2ζ√(1−ζ²))(ζ<0.707 出现)、相位滞后(每极点最多 −90°);⑤闭环传函 T=G/(1+GH):回路增益越大,误差与扰动被压得越狠;闭环带宽=−3dB 处"跟得上输入的最高频率";幅频看精度与带宽,相频看稳定与裕度。

📌 本节小结:①拉氏变换把微积分翻译成 s 域代数,极点位置决定系统性格;②一阶看时间常数 τ(63.2%/3τ/5τ),二阶由 ζ 与 ωn 二维定性格:ζ 管超调、ωn 管快慢;③Bode 三指标:−3dB 带宽(跟得动)、谐振峰(共振风险)、相位滞后(稳定性);④整定=搬极点,带宽法公式背后就是这里的数学。

🤔 思考题:1. 把 ζ 固定 0.707,ωn 从 6 拖到 20:超调变化吗?上升时间呢?这解释了"带宽与阻尼可以独立整定"吗?
2. 为什么机械谐振(柔性关节)会在 Bode 图上表现为高频谐振峰?此时提高速度环增益可能发生什么?