4. 模型降阶核心概念

模型降阶不是简单地把方程删掉几项，而是要在保留系统关键动态行为的前提下，把一个高维模型压缩成低维模型。要做到这一点，必须先理解四个基础概念：状态空间描述、传递函数、模态分析，以及降阶误差度量与稳定性保持。它们分别回答“系统怎么描述”“输入输出关系如何”“哪些动态最重要”“降阶后到底差多少、还稳不稳”这四个问题。

状态空间：降阶的出发点

一个线性时不变系统通常写成状态空间形式：

ẋ(t) = A x(t) + B u(t)
y(t)  = C x(t) + D u(t)

其中 x 是 n 维状态向量，u 是输入，y 是输出。A 描述状态之间的耦合与演化，B 描述输入如何进入系统，C 描述状态如何映射到输出，D 是直通项。模型降阶的目标，就是构造一个 r 维（r ≪ n）的近似系统：

ẋ_r(t) = A_r x_r(t) + B_r u(t)
y_r(t)  = C_r x_r(t) + D_r u(t)

使 y_r 在关心的频段和时间内尽量逼近 y。这里要特别注意：降阶不是只把 A 矩阵截断，因为 B、C 也决定了输入输出通道。只保留 A 的一部分而忽略 B、C 的投影关系，往往会导致输出误差很大。

传递函数：从频域看降阶目标

对状态空间做拉普拉斯变换，在零初始条件下可得传递函数：

G(s) = C (sI − A)⁻¹ B + D

传递函数把系统看成“输入到输出”的频域映射。降阶在频域中的核心要求是：在关键频段内，G_r(jω) 与 G(jω) 的幅频和相频特性接近。例如一个柔性机械臂模型，前若干阶低频模态决定主要运动，高频模态影响较小；降阶时应优先保留低频段增益和相位，而不是平均地压缩所有频段。

实际判断时，常看三个量：幅值误差 |G_r(jω)| − |G(jω)|、相位误差 ∠G_r(jω) − ∠G(jω)，以及峰值频率是否偏移。若在 0.1～100 rad/s 内幅值误差小于 1 dB、相位误差小于 5°，通常可认为该频段降阶效果可接受——具体阈值仍要由应用决定。

模态分析：找出“谁在主导动态”

模态分析从 A 矩阵的特征值出发。特征值 λ_i = σ_i + jω_i 对应系统的第 i 个模态：实部 σ_i 决定衰减或发散速度，虚部 ω_i 决定振荡频率。特征向量则描述该模态在各状态上的分布。对结构或电路类模型，常见做法是：

• 计算 A 的特征值，按 |Re(λ_i)| 或参与因子排序；
• 保留主导低频、弱阻尼模态，截去高频、强阻尼、能量贡献小的模态；
• 用左右特征向量构造投影矩阵 V、W，得到 A_r = WᵀAV、B_r = WᵀB、C_r = CV。

这就是模态截断与投影降阶的基本思路。需要强调的是，保留哪些模态不能只看特征值大小，还要看该模态在输入 B 和输出 C 上的可观测、可控程度。一个模态即使衰减很慢，如果它既不被输入激励、也不出现在输出中，对输入输出误差影响就很小。

降阶误差度量：不能只看“像不像”

降阶后必须量化误差。常用的度量有三类。第一类是时域误差，给定典型输入 u(t)，比较原模型输出 y(t) 与降阶模型输出 y_r(t)，计算最大绝对误差或均方根误差：

e_max = max_t |y(t) − y_r(t)|
e_rms = sqrt( (1/T) ∫₀ᵀ (y(t) − y_r(t))² dt )

第二类是频域误差，在关心频段 [ω₁, ω₂] 内计算幅值误差和相位误差，或直接看 H∞ 范数 ‖G − G_r‖∞。第三类是相对误差，例如 e_rms / ‖y‖_rms，便于不同量纲之间比较。实际工程中，建议至少同时看时域和频域两类指标：时域看瞬态和稳态，频域看带宽和共振峰。

稳定性保持：降阶不能把系统“降坏”

原系统稳定，不代表降阶后一定稳定。模态截断如果截掉了关键阻尼模态，或投影基选取不当，可能使 A_r 出现右半平面特征值。因此降阶后必须检查 A_r 的特征值是否全部满足 Re(λ_i) < 0。若发现不稳定，通常有三种处理方式：增加保留模态数 r、改用保持稳定性的降阶方法（如平衡截断及其稳定变体）、或对 A_r 做稳定性修正后再重新校验误差。

把四个概念串起来看：状态空间给出完整模型，传递函数给出频域目标，模态分析找出主导动态，投影降阶构造低维模型，误差度量和稳定性检查则决定这个低维模型能不能真正投入使用。对实时仿真而言，降阶模型不仅要“小”，还要“准”和“稳”，三者缺一不可。