从经典控制理论到现代控制理论，全方位梳理机械工程控制基础公式体系。本文档旨在为工程师与学生提供详实的数学推导、物理意义阐释及工程实践案例，解决复杂机械系统的建模、分析与设计难题。
在机械工程领域，控制理论是确保机械设备高精度、高稳定性运行的核心学科。无论是数控机床的伺服进给系统，还是航空发动机的燃油调节系统，其背后都依赖于严密的数学模型与算法支撑。机械工程控制基础公式并非孤立的数学符号，而是对物理世界能量转换、信号传递及反馈调节过程的抽象表达。
以传递函数为核心，主要研究单输入单输出（SISO）线性定常系统。重点掌握频率响应法、根轨迹法及劳斯判据，适用于传统机械传动系统的分析与校正。
基于状态空间法，能够处理多输入多输出（MIMO）、非线性及时变系统。通过矩阵运算描述系统内部状态，是现代机器人、智能车辆控制的基础。
将连续控制算法离散化，利用微处理器或 DSP 芯片实现实时控制。涉及 Z 变换、差分方程及采样定理，是嵌入式控制系统的核心。
传递函数（Transfer Function）是线性定常系统在零初始条件下，输出量的拉普拉斯变换与输入量的拉普拉斯变换之比。它是分析机械系统动态特性的最常用工具。
对于一般的线性微分方程：
aₙ(dⁿy/dtⁿ) + ... + a₁(dy/dt) + a₀y = bₘ(dᵐu/dtᵐ) + ... + b₀u
其传递函数 G(s) 定义为：
G(s) = Y(s)/U(s) = (bₘsᵐ + ... + b₀) / (aₙsⁿ + ... + a₀)
考虑一个典型的机械振动系统，质量为 m，阻尼系数为 c，弹簧刚度为 k。根据牛顿第二定律：
m(d²x/dt²) + c(dx/dt) + kx = F(t)
进行拉普拉斯变换（假设初始条件为零），得到传递函数：
G(s) = X(s)/F(s) = 1 / (ms² + cs + k)
标准化后，可识别出自然频率 ωₙ = √(k/m) 和阻尼比 ζ = c/(2√(mk))。
稳定性是控制系统正常工作的首要条件。对于闭环系统，其特征方程的根（即极点）必须全部位于 s 平面的左半平面（实部小于零）。
劳斯判据是一种代数判据，无需计算特征方程的根，仅通过特征方程系数即可判断稳定性。
设特征方程为：D(s) = a₀sⁿ + a₁sⁿ⁻¹ + ... + aₙ = 0 (其中 a₀ > 0)。
| sⁿ行 | a₀ | a₂ | a₄ | ... |
|---|---|---|---|---|
| sⁿ⁻¹行 | a₁ | a₃ | a₅ | ... |
| sⁿ⁻²行 | b₁ | b₂ | 0 | ... |
| sⁿ⁻³行 | c₁ | c₂ | 0 | ... |
其中，b₁ = (a₁a₂ - a₀a₃)/a₁，以此类推。如果劳斯表第一列所有元素均为正数，则系统稳定；若第一列元素符号发生改变，改变的次数等于特征方程在右半平面的根的个数。
基于复变函数的幅角原理，通过分析开环频率特性 G(jω)H(jω) 的极坐标图（奈氏图）来判定闭环稳定性。
判据公式：Z = P - 2N。其中 P 为开环传递函数在右半平面的极点数，N 为奈氏曲线逆时针包围 (-1, j0) 点的圈数，Z 为闭环系统在右半平面的极点数。若 Z=0，则闭环系统稳定。
当系统某个参数（通常是开环增益 K）从零变化到无穷大时，闭环特征根在 s 平面上移动的轨迹。通过分析根轨迹，可以直观地看出系统稳定性随参数变化的趋势。
随着机械系统日益复杂（如多关节机器人），经典控制理论的传递函数方法显得力不从心。状态空间法引入状态变量，将高阶微分方程转化为一阶微分方程组，非常适合计算机求解和多变量系统设计。
状态方程：ẋ = Ax + Bu
输出方程：y = Cx + Du
其中 x 是 n×1 状态向量，u 是 r×1 输入向量，y 是 m×1 输出向量。A, B, C, D 分别为系统矩阵、输入矩阵、输出矩阵和直接传输矩阵。
倒立摆是典型的非线性、强耦合、不稳定系统。将其线性化后，选取小车位置 x、摆杆角度 θ 及其角速度为状态变量，可建立如下状态方程：
[ ẋ ] [ 0 1 0 0 ] [ x ] [ 0 ]
[ ẍ ] = [ 0 0 -m g/M 0 ] [ v ] + [ 1/M ] u
[ θ̇ ] [ 0 0 0 1 ] [ θ ] [ 0 ]
[ θ̈ ] [ 0 0 (M+m)g/ml 0 ] [ ω ] [ -1/ml ]
此模型可直接用于 LQR（线性二次型调节器）等现代控制算法的设计。
尽管现代控制理论层出不穷，但PID 控制器（比例 - 积分 - 微分）因其结构简单、鲁棒性强，依然是机械工程中应用最广泛的控制策略。
输出与误差成正比。增大 Kp 可加快响应速度，减小稳态误差，但过大会导致系统振荡甚至失稳。
输出与误差的积分成正比。只要有误差存在，积分作用就会不断累积，直至误差为零。但会引入相位滞后，降低稳定性。
输出与误差的变化率成正比。具有“预见性”，能预测误差趋势并提前制动。对噪声敏感，需谨慎使用。
在实际工程中，常用的整定方法包括：临界比例度法（先使系统产生持续振荡，记录临界增益 Ku 和临界周期 Tu，再查表取值）、衰减曲线法以及Ziegler-Nichols 经验公式。对于复杂的机械系统，往往需要结合实验数据进行微调。
理论公式必须通过仿真验证才能指导工程实践。以下展示如何利用 MATLAB 快速验证上述控制策略。
% 定义传递函数 num = [1]; den = [1 2 1]; % (s+1)^2 sys = tf(num, den); % 定义 PID 控制器 Kp = 10; Ki = 5; Kd = 2; C_pid = pid(Kp, Ki, Kd); % 构成闭环系统 T_cl = feedback(C_pid sys, 1);
% 绘制阶跃响应
figure;
step(T_cl);
grid on;
title('闭环系统阶跃响应');
% 获取性能指标
S = stepinfo(T_cl);
disp(['上升时间：', num2str(S.RiseTime), ' s']);
disp(['超调量：', num2str(S.Overshoot), ' %']);
disp(['调节时间：', num2str(S.SettlingTime), ' s']);A: 这通常是因为系统中缺少积分环节（I 作用）。对于阶跃输入，只有 I 型或更高型的系统才能消除稳态误差。请检查你的控制器中是否包含积分项，或者被控对象本身是否具有积分特性。
A: 全零行表明特征方程存在关于原点对称的根（如纯虚根或对称实根）。此时应使用全零行的上一行构造辅助多项式，对其求导得到新的系数代替全零行，继续计算劳斯表。
A: 传递函数法只关注输入输出关系，无法反映系统内部状态，且难以处理多输入多输出系统。状态空间法能完整描述系统内部动态，适合计算机求解和现代控制设计，两者可以通过约当标准型相互转换。