工程实验室：控制、滤波、ODE 与 LP 面板

Engineering Lab 是一个横跨 11 个面板的标签栏：Control、Control Design、Filter Design、Signal Processing、ODE / Phase、Optimization、Curve Fitting、PDE Solver、MATLAB Kernel、Python / ML 和 R Kernel。每个面板保留自己的路由和读数。Control 涵盖伯德图、奈奎斯特图、根轨迹以及带增益裕度和相位裕度的阶跃/冲激分析；Control Design 增加离散化、极点配置、LQR 和 PID 整定；Filter Design 输出双二阶系数和截止频率扫描；ODE 面板在相图上对平衡点分类；Optimization 求解线性规划和整数规划，并为二变量模型的可行域着色。面板之间交换带类型的结果：一个面板建立的传递函数会成为另一个面板的受控对象系数。

更新 2026-10-015 个来源约 5 分钟读完

本文是草稿,尚未经过复核。

如何运行面板

选择面板——从标签栏选择，或直接打开它的路由。
输入数据。Control 接收按降幂排列的分子和分母系数（1/(s(s+1)) 对应 1 和 1 1 0）；Filter 接收滤波器族、节数和 Fs/Fc；Optimization 接收每行一条约束的文本模型。
阅读卡片，而不只是曲线。稳定性、裕度、极点和零点位于绘图上方；阶跃指标出现在 Step 标签页；LP 状态、目标函数和变量出现在可行域上方。
切换视图标签——回路有 Bode、Nyquist、Root Locus、Step 和 Impulse；PDE 有 Heat 1-D、Wave 1-D 和 Laplace。
把结果传下去。从 Filter 或 Control Design 发送的传递函数会填入 Control 的系数字段。

计算示例

以下读数全部来自实时面板。

面板	输入	读数
Control	G(s) = `1` / `1 1 0`，K = 1（默认值）	Stable ✓；相位裕度 51.83° @ 0.786 rad/s；增益裕度 ∞（可证明不存在 −180° 穿越）；闭环极点 -0.5 + 0.866j 和 -0.5 − 0.866j；零点 —
Control · Step	同一闭环	OS 16.3%，上升时间 1.63s，整定时间 8.08s，终值 1
Filter Design · IIR	butterworth lowpass，2 节，Fs 48000 Hz，Fc 1000 Hz（默认值）	spec IIR butterworth lowpass, 2 biquad(s), Fs=48000, Fc=1000；第一节 `b0 0.00407407, b1 0.00814814, b2 0.00407407, a1 -1.88856, a2 0.904852`；扫描第 2 行读到 996.189 Hz
Optimization	`max: 3x + 2y` · `x + y <= 10` · `x <= 8` · `y <= 6`	Optimal ✓，目标函数 28，x 8，y 2；紧约束 `x + y <= 10 ; x <= 8`
ODE · Phase portrait	x' = `y`，y' = `-2*sin(x)`（单摆预设）	(-6.283, 0) 中心，(-3.142, 0) 鞍点，(4.424e-21, 0) 中心（`sin x` 的数值零点），(3.142, 0) 鞍点，(6.283, 0) 中心

默认开环在 0.786 rad/s 处穿越 0 dB，给出 51.83° 的相位裕度，闭环阶跃响应超调 16.3%，随后在 8.08s 内整定。

卡片描述的是哪个回路

Bode、Nyquist 和裕度卡片使用开环 K·G(s)；稳定性卡片和极点列表描述单位反馈闭环，Step/Impulse 默认使用该闭环，并提供 Open/Closed 切换。这就是为什么 1/(s(s+1)) 虽然受控对象在原点有极点，却显示 Stable ✓。K 缩放开环；可选的 PID 模块在同一分析之前放置 Kp、Ki 和 Kd。

需要记住的限制

数值读数带有方法限制。: 裕度来自扫描的伯德响应，当穿越可证明不存在时面板会说明——默认增益裕度显示 ∞。极点和零点是数值根；阶跃指标会在尾部尚未整定时给出警告。
优化器是线性求解器。: 目标函数和约束必须是线性的，整数用 int x y 声明。变量按非负处理：在 x <= 5 下 max: -x 返回 x = 0 处目标函数 0，因此带符号的量需要重新参数化。Optimal、infeasible 和 unbounded 描述的是模型，而不是现实过程。
它不做的事：电路、单位或符号运算。: 它只接收系数列表和普通数字；不模拟电路图、不换算单位、不推导公式。电路、RLC 阻抗、史密斯圆图和 FFT 在 Science Lab；符号微积分、矩阵和分步代数在 Math Tools；快速表达式请用 scientific calculator。

双二阶数量不是滤波器阶数

典型的陷阱是节数：Filter Design 的阶数字段统计的是二阶节数量，标签也写明了——Biquads (order = 2×)。默认值 2 得到的是 4 阶滤波器，而不是 2 阶。spec 行显示 IIR butterworth lowpass, 2 biquad(s), Fs=48000, Fc=1000，表格列出两节，第 2 行的扫描读到 996.189 Hz。类型 1 时 spec 缩小为 1 biquad(s)；类型 13 时面板设计 12 节，并警告 Order 13 exceeds the limit of 12; designed with 12。把设计对照教科书表格时，请比较 spec 行，而不是你输入的数字。

适用场景

核对控制教程

对于 G(s) = 1/(s(s+1))、K = 1，裕度卡片显示 51.83°，阶跃卡片显示 16.3% 超调——典型的弱阻尼二阶组合。根轨迹和奈奎斯特图显示 K 变化时极点的走向。

试制数字滤波器

设计 Fs = 48000 Hz、Fc = 1000 Hz 的低通滤波器，读取双二阶节和 1–12 的截止频率扫描，然后复制 C 或 Python 系数块，或把 H(z) 连同采样率一起发送出去。

模型、规划与动态

默认 LP 返回 x = 8、y = 2、目标函数 28；整数预设返回 x = 3、y = 2、目标函数 5。ODE 面板无需手动线性化即可对单摆的平衡点分类，拟合和 PDE 面板覆盖最小二乘拟合、系统辨识以及热、波动和拉普拉斯问题。

参考文献

Wikipedia, Bode plot, en.wikipedia.org （访问日期：2026-10-01）— 增益和相位裕度。
Wikipedia, Butterworth filter, en.wikipedia.org （访问日期：2026-10-01）— −3 dB 截止。
Wikipedia, Digital biquad filter, en.wikipedia.org （访问日期：2026-10-01）— 二阶节。
Wikipedia, Linear programming, en.wikipedia.org （访问日期：2026-10-01）— 非负变量。
Wikipedia, Phase plane, en.wikipedia.org （访问日期：2026-10-01）— 平衡点。