调试伺服速度环的时候，最磨人的不是电流环，而是那个写在驱动参数表里的“惯量比”。同一台永磁同步电机，空载跑得稳，装上负载之后超调、低频振荡全跑出来了。机械臂抓不同重量的工件，卷绕机构料卷直径不断变化，运行中的负载惯量根本不是一个固定值，离线测一次根本不够用。永磁同步电机（PMSM）在 Simulink 仿真里有一个很实用的在线辨识方案——模型参考自适应法（MRAS），它不干扰正常运行，不需要额外叠加激励信号，就能实时估计转动惯量。这篇文章会把从公式推导、模型搭建、参数整定到波形判读、坑点排查的完整过程摊开讲，照着搭一遍，你也能看到估计值在加速段迅速收敛到真实转动惯量的那一幕。适合正在做伺服自整定、或者刚接触惯量辨识算法验证的朋友。
1. 为什么伺服系统必须知道转动惯量：速度环整定背后的物理逻辑
1.1 速度环增益与惯量比：藏在参数表里的核心关系
速度环的本质，是通过调整电磁转矩让机械速度跟随指令。机械子系统的传递函数近似为 1/(Js+B)，这里 J 是转动惯量，B 是阻尼系数。J 越大，同样的转矩变化产生的加速度就越小，系统对速度指令的响应越迟钝。
速度环比例增益 Kp 的作用，简单说就是“把机械惯量补偿到目标带宽”。理想情况下，速度环开环增益与 J 成正比：Kp 的设计值大约等于 J × ωc / Kt，其中 ωc 是期望的闭环带宽，Kt 是转矩系数。这个关系意味着，如果负载惯量从 1 倍变成 4 倍，而 Kp 还保持原值，那实际闭环带宽直接缩水到原来的四分之一，动态响应慢半拍；反过来，如果惯量变小而 Kp 没变，相位裕度下降，系统就会在某个速度段开始振荡。这也是为什么各家伺服驱动器的自整定功能都把“惯量比识别”放在第一位。
以前我调一套设备，空载时速度环参数整得很好，阶跃响应干净利落。装上一个夹具后，同样参数直接低频啸叫，曲线来回甩。后来逐一排查，发现就是负载惯量翻了三倍多。那一刻我才真正理解，所谓“速度环整定”，本质上是在和转动惯量这个物理量打交道。
1.2 离线辨识的局限与在线辨识的不可替代
常规离线辨识方法并不少：加减速时间常数法、转矩阶跃响应法、扫频法，原理上都能得到 J。但这些方法有一个共同前提——系统要停下来，专门执行一段辨识工况。对产线上的设备、无人机云台、机器人关节来说，停机辨识意味着节拍损失，有时候根本不具备停机条件。
更关键的是，停机辨识只能得到“当前状态”的惯量。机械臂抓取工件前后惯性负载完全不同；卷绕设备的料卷半径随生产不断增大，惯量持续变化。离线一次根本覆盖不了这些场景。在线辨识的价值就在于：在设备正常动作过程中，利用现有的速度动态信息实时估计 J，不需要中断运行，甚至可以直接把辨识结果喂给速度环，实现参数自整定。这正是模型参考自适应法在工程中被高频使用的原因。
2. 环顾一圈选型：MRAS凭什么能在转动惯量辨识里站C位
2.1 常见在线辨识思路速览
在线辨识转动惯量，业界常用的思路大约有三种：递推最小二乘（RLS）、扩展卡尔曼滤波（EKF）、模型参考自适应（MRAS）。它们在原理复杂度、计算量、可靠性上的取舍差别很大，先列个对比表：
| 方法 | 原理复杂度 | 计算量 | 激励要求 | 调试难度 | 典型场景 |
|---|---|---|---|---|---|
| 递推最小二乘 | 中 | 中 | 高，需要充分激发 | 中，需调遗忘因子 | 参数在线辨识、离线数据拟合 |
| 扩展卡尔曼滤波 | 高 | 高 | 较低 | 高，协方差矩阵难调 | 高精度状态估计、无传感器控制 |
| 模型参考自适应 | 低 | 低 | 中 | 低，直观易调 | 电机参数辨识、速度环自整定 |
递推最小二乘的思路是把运动方程改写成线性回归模型，把 1/J 和负载转矩当作未知系数，利用历史数据递归更新。它的精度确实高，但非常依赖输入信号持续激励，遗忘因子设大了稳不住，设小了跟不上慢变参数，调试周期不短。
扩展卡尔曼滤波通过扩维状态向量来同时估计速度、位置、惯量，理论上对噪声有最优抑制，但协方差矩阵 Q、R 的整定完全靠经验，计算量在普通 MCU 上吃紧，一般用于对精度要求极高的场合。
MRAS 的核心优势在于结构简洁：一个可调模型加一个自适应律，主要的数学运算就是乘法和积分，单片机里写起来非常轻松。它的稳定性有 Popov 超稳定性理论兜底，工程上调试也直观——曲线不对，调两个增益就行。对 Simulink 仿真验证来说，MRAS 可以用最少的模块搭出完整闭环，特别适合理解算法机理。
2.2 为什么在 Simulink 里选 MRAS 做验证
Simulink 里搭 MRAS 辨识模块，本质上只需要积分器、乘法器、PI 控制器、限幅、倒数和滤波这几个基础模块，整个辨识回路不超过十个模块。这比搭 EKF 那一堆矩阵运算要轻松太多，而且每个中间信号都可以拉出来看波形，调试体验非常友好。
另一个现实原因是可迁移性。Simulink 里验证通过的算法，改成嵌入式 C 代码后结构几乎能一对一对应，控制器实现时不容易走样。很多伺服驱动器产品里自带的“惯量辨识”功能，底层就是 MRAS 或者它的变体。在仿真阶段把这套东西吃透，后面做实物调试会从容很多。
3. MRAS辨识运动方程与自适应律推导：从电机方程到能放进积分器的公式
3.1 参考模型与可调模型怎么构建
先建立 PMSM 的转矩和运动方程。表贴式永磁同步电机在 dq 坐标系下的电磁转矩为：
T_e = 1.5 × p × ψ_f × i_q
其中 p 是极对数，ψ_f 是转子磁链，i_q 是 q 轴电流。如果是内置式永磁电机，还要加上磁阻转矩项 1.5 × p × (L_d - L_q) × i_d × i_q，这个后面踩坑部分再细说。
机械运动方程：
J × dω_m/dt = T_e - T_L - B × ω_m
其中 ω_m 是机械角速度，T_L 是负载转矩，B 是粘滞阻尼系数。
MRAS 的思路是把真实电机当作“参考模型”——真实系统里 J 和负载决定了实际的 ω_m 响应，这个响应就是标准答案。然后我们搭一个“可调模型”：用估计值 Ĵ（实际估计的是 1/Ĵ）去重构一个 ω̂_m，输入端接入同样的 T_e。可调模型的方程为：
Ĵ × dω̂_m/dt = T_e - T̂_L - B̂ × ω̂_m
通过自适应律不断修正 1/Ĵ，让 ω̂_m 实时跟踪 ω_m。只要可调模型的输出与实际速度误差趋于零，估计值就收敛到了真实 J。
这里要专门提醒一个误区：很多人一搭 MRAS 就习惯性地“参考模型”也自己搭一套，其实完全没必要。在 Simulink 仿真里，PMSM 本体模块本身就是参考模型，它的速度输出已经是真实 J 作用下的结果，你只需要在旁边附加一个可调模型就够了。重复搭参考模型不仅费模块，还容易因为参数不一致引入额外误差。
3.2 自适应律推导：Popov超稳定性给出的“操作手册”
先做最简单的假设：忽略阻尼 B，且仿真工况负载转矩 T_L = 0。此时参考模型和可调模型分别简化为：
dω_m/dt = (1/J) × T_e dω̂_m/dt = (1/Ĵ) × T_e
定义速度误差 e = ω_m - ω̂_m，参数估计为 a = 1/J，â = 1/Ĵ，参数误差 φ = â - a。对误差求导：
de/dt = (1/J - 1/Ĵ) × T_e = -φ × T_e
要让误差收敛，同时参数误差也收敛，Popov 超稳定性理论给了一个工程上很好用的结论：只要把系统转化为一个前向通路严格正实、反馈通路满足积分不等式的结构，再选择比例积分形式的自适应律，就能保证误差渐进收敛。具体推导过程这里不展开了，直接给出结论形式的自适应律：
â(t) = â(0) + K_p × q(t) + K_i × ∫₀ᵗ q(τ) dτ
其中 q(t) = T_e(t) × e(t)，â(0) 是 1/Ĵ 的初值。最终估计的转动惯量为：
Ĵ = 1 / â
这个公式看着简单，信息量其实很大。q(t) 是误差与转矩的乘积，它决定了参数修正的方向和力度；K_p 负责对误差的即时响应，K_i 负责累积消除稳态偏差。比例积分结构保证了无论初值偏差多大，只要 T_e 提供了足够的激励，â 就会朝着真实值 1/J 移动。
3.3 符号检查：自适应律方向搞反会怎样
自适应律方向如果连反，辨识值不会缓慢漂移，而是直接爆掉。调试时有个非常实用的检查方法：先把自适应律断开，手动固定 â = 1/Ĵ0，给电机一个阶跃转矩，同时观察 ω_m 和 ω̂_m 的响应曲线。
如果可调模型的估计惯量设得比真实值大，ω̂_m 的加速度会明显慢于 ω_m，此时 q = T_e × e 为正；反之 q 为负。明确了 q 的符号含义之后，再去检查积分器和乘法的极性，确保自适应律让 q 朝“缩小误差”的方向走。我是建议第一次搭的时候故意把符号接反一次，观察 Ĵ 发散的过程，这样你以后看到曲线爆炸时能第一时间想到符号问题。
4. Simulink模型搭建逐块拆解：从PMSM本体到自适应律的最小闭环
4.1 仿真整体结构与电机参数示例
搭建之前先给一个推荐的仿真结构，从上到下分四层：
- 电机本体层：PMSM 模块、逆变器、SVPWM 或者理想占空比给定
- 电流环层：id = 0 控制，dq 轴电流 PI
- 速度环层：速度 PI 输出 q 轴电流给定
- MRAS 辨识层：输入 T_e 和 ω_m，输出 Ĵ 的估计值
为了让辨识算法能独立调试验证，我建了一个中等功率的 PMSM 仿真模型作为示例，参数如下：
| 参数 | 数值 |
|---|---|
| 额定功率 | 2 kW |
| 极对数 p | 4 |
| 定子电阻 R_s | 1 Ω |
| dq 轴电感 L_d = L_q | 5 mH |
| 转子磁链 ψ_f | 0.25 Wb |
| 真实转动惯量 J | 0.005 kg·m² |
| 阻尼系数 B | 0 |
| 母线电压 | 311 V |
| PWM 开关频率 | 10 kHz |
辨识工况设计为：速度给定从 0 斜坡上升到 1000 r/min，上升时间 3 秒，负载转矩保持 0。这种持续加速的工况能提供足够强的转矩激励，是 MRAS 收敛的理想环境。
4.2 可调模型模块的实现细节
可调模型在 Simulink 里的实现非常直接：转矩信号 T_e 除以当前估计的 Ĵ，得到一个角加速度，送入积分器积分得到 ω̂_m。
注意这里有两个细节容易忽略：
第一，积分器初值必须设成与实际电机初始速度一致。仿真从零速启动的话，ω̂_m(0) 直接设 0 即可。如果电机起始速度不为零，而积分器从 0 开始，初始那一瞬间误差 e 会很大，自适应律会被一个虚假的大误差冲击，导致估计曲线先乱跳一阵。
第二，可调模型里如果引入了 -Bω̂_m 项，B 的取值必须和参考模型里的真实阻尼一致，否则会引入系统性偏差。前面参数表里把 B 设成 0，就是为了在原理验证阶段排除这个干扰因素。工程上电机高速时粘滞摩擦不能忽略，那是后续进阶要处理的问题。
4.3 自适应律与信号流接线清单
自适应律部分建议拆成几个模块分别搭建，方便调试时看中间信号。整体连线关系如下表：
| 信号 | 来源 | 去向 | 注意事项 |
|---|---|---|---|
| ω_m | PMSM 模块输出 | 误差计算、示波器 | 真实速度，单位 rad/s |
| T_e | PMSM 模块输出或 1.5pψ_f i_q | 可调模型、q 计算 | 表贴式忽略磁阻项 |
| ω̂_m | 可调模型积分器输出 | 误差计算 | 初值必须与 ω_m(0) 一致 |
| e | 减法器 | q 计算 | e = ω_m - ω̂_m |
| q | 乘法器 | PI 控制器 | q = T_e × e |
| 1/Ĵ | PI 控制器输出 | 倒数模块、限幅 | 初值设 1/Ĵ0，必须限幅防除零 |
| Ĵ | 倒数模块输出 | 低通滤波、示波器 | 滤波后再做闭环自整定 |
PI 控制器输出的物理量是 1/Ĵ，初值应该设成估计惯量初值的倒数。比如我仿真时把 Ĵ(0) 设成 0.05，也就是真实值的 10 倍，那 1/Ĵ 的初值就是 20。这样设定之后，启动阶段积分器会从 20 出发，在加速过程中逐渐逼近真实值 1/0.005 = 200。
倒数模块前方必须加限幅。假设 PI 输出不小心越过零变成负值或接近零，倒数模块会产生巨大数值甚至报错。限幅范围我习惯设成 [1/J_max, 1/J_min]，对应估计惯量上下限，比如 [10, 1000]，这样既安全又不影响正常收敛。
5. 关键参数整定与辨识波形判读：别让坏参数毁掉好算法
5.1 从波形判断收敛质量：快、稳、准三个尺度
仿真跑起来之后，重点看四条曲线：ω_m 与 ω̂_m 的跟踪情况、误差 e、估计转动惯量 Ĵ、以及 1/Ĵ 的变化。理想的收敛过程有三个特征：
- 快：从启动加速到估计值进入真实值 ±5% 区间，时间一般在 0.5 秒以内，取决于 K_i 大小
- 稳：收敛过程中没有大幅振荡，超调不超过真实值的 20%，然后平滑逐步收敛
- 准：最终稳定值与真实 J = 0.005 一致，误差在 2% 以内
如果看到 Ĵ 曲线一直在降但降不动，多半是 K_i 太小，收敛时间太长；如果 Ĵ 曲线绕着真实值反复穿越、振荡收不了口，多半是 K_p 或 K_i 偏大，或者低通滤波器截止频率设太高。可以先用一个记录模块把 Ĵ 波形存下来，对比不同参数下的收敛曲线，比单纯看数字变化直观得多。
速度进入稳态之后，T_e 趋近于零（因为 T_L = 0），此时 q = T_e × e 也趋近于零，1/Ĵ 的积分器基本停止修正，Ĵ 会保持在一个平台值。这个平台值如果在合理范围，就说明加速段的辨识结果是可信的。如果平台段还在缓慢漂移，说明数值噪声正在通过积分器累积，这就是后面要讲的激励不足问题。
5.2 Kp与Ki的初选方法与调节顺序
MRAS 自适应律的增益整定，我总结了一套非常实用的顺序：
- 先令 K_p = 0，只保留纯积分环节。这样曲线最平滑、最不容易发散，先把 K_i 调出来。
- K_i 从比较小的值开始，比如 1。观察 Ĵ 收敛时间，如果 1 秒内没进入真实值附近，把 K_i 乘以 3~5 倍再试。反复几次找到“能收敛、速度可接受”的区间。
- 继续增大 K_i 直到 Ĵ 曲线出现等幅振荡或明显超调，然后退回振荡临界值的 1/3 左右，这就是一个比较稳妥的 K_i。
- 最后再引入 K_p。K_p 的作用是抑制误差突变，一般取 K_i 的 1/20 到 1/10。K_p 太小没效果，太大会放大噪声，让估计值带毛刺。
调完增益之后，再用不同初值 Ĵ(0) 做几组对比仿真。如果初值差 10 倍都能在 1 秒内收敛到同一真实值，说明自适应律增益设定合理。这个“初值鲁棒性”测试比单看一组波形更有说服力。
5.3 工况设计：为什么要给足够强的速度激励
MRAS 本质上靠输入信号激励辨识参数。T_e 为零的区域，q 始终趋近于零，误差信息无法传递给自适应律，辨识就会停止。所以仿真工况必须刻意设计成“有动态过程”的样子。
推荐的做法是让速度做持续 1~3 秒的斜坡或 S 型曲线上升，整个过程中 T_e 始终保持在正值附近，持续驱动可调模型修正参数。你也会发现，加速斜率越大，T_e 越大，收敛越快；加速时间太短，则辨识数据量不足，估计值可能来不及稳定。
如果你想在有负载的工况下辨识，建议让负载转矩在加速中段阶跃加入。这样速度动态中既包含 J 的信息，又包含 T_L 突变引起的加速度变化，辨识算法可以在一个动态过程中同时把两者分离出来。纯稳态工况下无法同时辨识 J 和 T_L，这是数学上不可观的问题，不是算法不好调。
6. 实测模拟中容易翻车的五个细节与排查策略
6.1 激励不足导致的漂移与冻结策略
匀速段出现的问题是典型的“激励不足”：加速度为零时，T_e 几乎等于负载转矩，速度接近恒定，此时速度误差 e 虽然存在，但 q 中包含的有效信息量极少。可调模型积分器在数值噪声下会缓慢累积，Ĵ 曲线呈现缓慢爬升或下降的漂移现象。
处理办法是给自适应律加一个使能逻辑。我习惯用速度微分的绝对值作为判据：当 |dω_m/dt| 大于某个阈值时，积分器正常工作；速度接近恒定时，把积分器冻结，保持上一次的辨识结果。Simulink 里实现起来很简单，算速度微分后过一个比较器和增益，再用 Enable 端口控制积分器即可。当然，速度微分会有噪声，阈值要留足余量，建议先滤波再取微分。
6.2 阻尼B与负载转矩的系统性偏差
如果真实电机存在粘滞阻尼 B，但可调模型里忽略了 -Bω_m 项，辨识结果会系统性偏小。原因很直观：实际运动方程里有额外的阻尼转矩消耗了一部分电磁转矩，可调模型把它全部算成惯性加速，为了匹配速度响应，模型就会高估转动惯量。
排查时有个特征可以参考：让电机正向加速和反向减速分别辨识一次，如果两次结果不对称，基本可以断定是阻尼或负载转矩项没有建模好。工程上我建议把 B 纳入可调模型，并配合负载转矩观测器先行估计 T_L，再交给 MRAS 修正 J。这样三个未知量都能在动态过程中逐步收敛，不过这就是一个完整的进阶方案了，仿真中先跑通无阻尼版本理解机理更重要。
6.3 转矩计算错误与相位滞后
仿真中最容易踩的坑不是算法本身，而是 T_e 算错。用库里的 PMSM 模块时，T_e 可以从模块输出端直接取，几乎不会错。但如果自己搭 dq 模型，很多人会忘记磁阻转矩项。表贴式电机 L_d = L_q 时那项确实是零；换成内置式电机，漏掉 1.5 × p × (L_d - L_q) × i_d × i_q 会让转矩估计产生明显偏差，辨识结果自然不可信。
相位滞后问题更隐蔽。电流环闭环、PWM 更新、传感器滤波都会让转矩与速度之间引入延迟。可调模型里的 T_e 是一拍无延迟的理想信号，而 ω_m 却带有各路延迟，误差 e 就包含了系统延迟带来的相位差。处理办法是在可调模型输入的 T_e 通道也加一个低通滤波器，让两条路径的延迟尽量匹配。仿真中先不加滤波器看差异，再加滤波器看变化，你能非常直观地感受到延迟对辨识精度的影响。
6.4 离散化与积分器数值问题
Simulink 里如果整体模型用离散求解器，MRAS 模块的采样时间建议与 PWM 周期同步，或者取一个固定的较低频率。积分器在高频离散下一般没有大问题，但倒数模块在输入接近零时会产生尖峰，这是我最想提醒的一点。
另外，速度信号如果通过数字微分来算加速度，微分会放大高频噪声，估计惯量的曲线会布满毛刺。MRAS 的好处恰恰在于它不需要显式微分，误差 e 由 ω_m 和 ω̂_m 直接相减得到，噪声耐受性比很多基于微分的方案好。因此搭模型时务必不要画蛇添足再去求 dω/dt，直接把速度信号接进去用就好。
6.5 让“辨识结果”滚起来：闭环自整定的做法
辨识出 J 之后，最自然的用法是把它喂回速度环，实现 PI 参数在线自整定。仿真里可以做一个简单的增益调度：把速度环 Kp 写成 Kp_set = J_hat × scale，其中 scale 是根据目标带宽算出来的常数；Ki 同理取 Kp 的几分之一。
这样 J_hat 在辨识初期从初值跳到真实值的过程中，速度环 PI 也会跟着动态调整。你会发现一个很有意思的现象：当辨识值偏大时，速度环 Kp 偏大，系统短暂变激进；辨识收敛后 Kp 回归合理区间，速度响应重新变得平滑。这个闭环动态本身就是一套完整的“自整定”功能雏形，也是把仿真成果接到真实伺服驱动器的必经练习。
最后聊一点仿真之外的感受。很多教程把 MRAS 辨识做得非常理想化，那是因为仿真里没有测速噪声、没有电流环延时、没有滤波器相位滞后。自己从 Simulink 移植到实际控制器时，我第一次上电就发现辨识值在低速段乱飘，后来把速度测量的低通和自适应律的 PI 带宽一起降下来才稳住。建议复现的时候也先冻结积分器、降低开环增益，确认低速小负载下系统不出问题再逐步放开。辨识算法本身不复杂，复杂的是把它驯化成一台能长期稳定运行的参数自整定功能，这才是真正值钱的地方。先跑通 Simulink，再一点点往模型里加“现实世界的杂质”，这条路走踏实了，后面做什么伺服项目都会受益。
