简介面向电力电子与电气工程领域的科研人员与控制算法开发者这份 PDF 文献探讨内置式永磁同步电机IPMSM基于模型的无传感器控制中的参数漂移问题提出在估计旋转 γ-δ 坐标系下使用仿射投影算法APA在线辨识定子电阻与 dq 轴电感并以双时间尺度策略分别完成快速与慢速参数更新使扩展反电动势观测器在温度变化和磁饱和工况下仍保持稳定。压缩包内仅含 1 个 PDF 文件约 1.51MB为 IEEE Transactions on Power Electronics 论文原文涵盖摘要、索引术语、算法推导、MATLAB/Simulink 仿真以及基于 TMS320F28335 DSP 实验平台的验证数据。目前已有 62 人学习。读者可据此复现 APA 参数辨识流程理解参数可辨识性分析与观测器自适应修正的衔接并延伸至状态监测、故障诊断等场景适合科研选题、论文精读与工程调试参考。1. 从一台低速重载 IPMSM 的观测器发散说起一台内置永磁同步电机IPMSM在低速大转矩工况下切到无传感器控制磁链观测器的位置估计开始周期性抖动转速环输出跟着晃电流波形却看不出明显异常。冷态标定的定子电阻 Rs、dq 轴电感 Ld/Lq、永磁磁链 ψf运行十几分钟后随铜温和磁钢温度偏移基于模型的无传感器控制完全依赖这些参数构造反电动势或磁链参数一偏位置估计就偏。在线参数辨识要做的是用控制器已经采到的电压、电流和转速信息在运行过程中把这些参数实时拉回真实值附近让观测器不靠重新停机标定也能撑住低速和温升场景。这套方法适合做伺服驱动器、车用电驱和工业变频器固件的人尤其是已经跑通无传感器控制、却在低速和温升下反复调观测器增益的团队。2. 内置永磁同步电机参数与基于模型无传感器控制的耦合关系在线辨识不是给电机控制器加一个独立模块就完事它和基于模型的无传感器观测器共享同一套 dq 轴方程。方程里哪些参数会漂、漂了以后观测器的哪一项先崩、哪些参数在同一次激励里根本不可辨这三件事决定了辨识器该在什么工况下启动、用什么激励、更新哪几个参数。下面把方程、敏感度和可辨识性拆开讲。2.1 IPMSM 的 dq 轴方程与四个会被在线辨识盯上的参数内置式永磁同步电机的 dq 轴电压方程写成ud Rs id Ld did/dt − ωe Lq iquq Rs iq Lq diq/dt ωe (Ld id ψf)方程里的 Rs、Ld、Lq、ψf 就是在线参数辨识的四个主要目标。它们不是常数Rs 随绕组温度近似线性上升Ld 和 Lq 随电流饱和下降ψf 随磁钢温度下降。下表给出典型量级和它们对无传感器观测器的影响路径方便判断当前工况下该优先辨识谁。参数典型量级漂移来源对基于模型观测器的影响Rs0.01–1 Ω铜绕组温度每升 1 ℃ 约 0.39%低速反电动势小Rs 误差直接进入磁链积分Ld0.1–10 mHd 轴磁路饱和id 增大时下降影响磁阻转矩和凸极位置估计Lq0.1–10 mHq 轴饱和与交叉饱和影响高频注入位置估计和 q 轴反电动势ψf0.05–0.3 Wb磁钢温度每升 1 ℃ 约 −0.1% 至 −0.15%决定反电动势幅值直接缩放位置误差从表里能看出一个工程顺序低速段 Rs 和 ψf 的辨识需求最急因为它们直接决定观测器能不能在反电动势很弱时稳住中高速段 Ld、Lq 的饱和误差开始影响转矩和效率但位置估计对 ψf 的绝对精度反而没那么敏感。2.2 参数失配在磁链观测器和反电动势观测器里怎样被放大用电压模型构造定子磁链时d 轴磁链写成 ψd ∫(ud − Rs id ωe Lq iq) dt。积分项里的 Rs 误差不会被衰减而是随时间累积成磁链偏差电机转速越低反电动势项越小这个积分误差在总磁链里的占比越大。反电动势观测器的 q 轴分量写成 Eq uq − Rs iq − Lq diq/dt − ωe Ld idψf 的误差直接进入 Eq再经过观测器增益放大成位置误差。滑模观测器的切换增益如果按额定 ψf 整定当 ψf 因温升下降后有效反电动势变小同等增益下抖振占比变大位置估计的波动随之增加。这就是为什么在线参数辨识不能只做一次观测器对参数的敏感度随转速和负载变化同一个 Rs 误差在额定转速下可能只造成零点几度偏差在 5% 额定转速下就能让磁链观测器发散。2.3 为什么四个参数不能在同一次激励里全部辨识出来把四个参数同时放进一个最小二乘问题会遇到秩亏。当控制器采用 id 0 控制时d 轴方程里 Ld 的系数是 did/dt稳态下该系数接近零Ld 不可辨低速时 ωe 很小ψf 在 q 轴方程里的系数 ωe 也很小ψf 和 Rs 难以分开。工程上常见的做法是分时、分频、分工况辨识而不是追求一次激励全辨。工况主导可辨识参数推荐激励方式零速或极低速Ld、Lq、Rs高频电压注入加直流衰减低速 5%–15% 额定转速Rs、ψfd 轴电流方波注入中高速Rs、ψf自然反电动势无需额外注入全速段ψf 温漂趋势观测器误差慢反馈把这张表当成辨识器的调度表控制器根据当前转速和 id 状态选择启动哪一路辨识避免在不可辨工况下强行更新参数。3. 在线参数辨识的三条工程主线去掉学术论文里的复杂激励工程上能落地的在线辨识基本围绕三条主线递推最小二乘RLS做慢变参数、扩展卡尔曼滤波EKF做状态与参数联合估计、高频注入做电感辨识。三条线不是互斥的而是按工况切换。3.1 用递推最小二乘在线辨识 Rs 与 ψf在 q 轴方程里把 diq/dt 当作可测量或可差分得到定义中间量 y uq − ωe Ld id − Lq diq/dt方程化为 y Rs iq ωe ψf。写成标准线性形式 y φᵀθ其中 φ [iq, ωe]ᵀθ [Rs, ψf]ᵀ。接下来用带遗忘因子的 RLS 递推import numpy as np class RLS: def __init__(self, theta0, P01e3, lam0.998): self.theta np.array(theta0, dtypefloat) # [Rs, psi_f] self.P np.eye(2) * P0 # 协方差矩阵 self.lam lam # 遗忘因子 def update(self, phi, y): phi np.array(phi, dtypefloat).reshape(2, 1) # [iq, omega_e] y float(y) # 增益向量 K self.P phi / (self.lam phi.T self.P phi) # 参数更新 err y - (phi.T self.theta.reshape(2, 1)).item() self.theta self.theta K.flatten() * err # 协方差更新 self.P (self.P - K phi.T self.P) / self.lam return self.theta.copy()逻辑说明每一拍用当前 iq 和 ωe 构造回归向量 φ用 q 轴电压方程算出的 y 作为观测值先算增益 K再修正 Rs 和 ψf最后按遗忘因子 λ 衰减协方差。参数说明λ 取 0.995–0.999越小追踪越快但参数噪声越大P0 取 1e2–1e4太大时初期参数跳变明显iq 不能长期为零否则回归向量第一列无效ωe 太小时第二列接近零ψf 不可辨此时应冻结 ψf 更新。3.2 扩展卡尔曼滤波把转速、位置和参数放进同一个状态向量EKF 的做法是把 id、iq、ωe、θe、Rs、ψf 全部扩进状态向量 x用 dq 轴电流方程做预测用电流采样做更新。离散预测骨架如下import numpy as np def ekf_predict(x, u, Ts, Ld, Lq): # x [id, iq, we, the, Rs, psi_f] id_, iq_, we, the, Rs, psi_f x ud, uq u did (ud - Rs * id_ we * Lq * iq_) / Ld diq (uq - Rs * iq_ - we * (Ld * id_ psi_f)) / Lq x_new x Ts * np.array([did, diq, 0.0, we, 0.0, 0.0]) # 工程实现里再按上述方程对 x 求偏导得到雅可比 F return x_new逻辑说明预测步只更新电流和电角度转速、Rs、ψf 的过程模型设为零漂移靠协方差传播允许缓慢变化更新步用实测 id、iq 纠正状态。参数说明参数项的 Q 取 1e-10–1e-8太大时参数会跟着电流噪声漂电流测量 R 取 1e-4–1e-2太小容易滤波发散。EKF 的短板是参数与状态耦合强当 ψf 和 Rs 同时可辨性差时需要人为限制参数协方差的下降速度。3.3 高频注入辨识 Ld 与 Lq 的参数窗口在 d 轴注入高频电压 udh Vh sin(ωh t)测量 d 轴高频电流幅值 Idh则 Ld ≈ Vh / (ωh Idh)。q 轴同理在 q 轴注入后测 Iqh 得 Lq。频率选择要避开电流环带宽和 PWM 开关频率附近的边带常见 500 Hz–2 kHz幅值取额定电压的 5%–10%太小信噪比不够太大带来转矩脉动和噪声。注入参数建议范围作用风险注入频率500 Hz–2 kHz区分基波与高频响应接近电流环带宽时振荡注入幅值5%–10% 额定电压保证电流信噪比过大引起转矩脉动注入轴d 轴测 Ldq 轴测 Lq分离交叉饱和同时注入会互相干扰持续时间50–200 ms完成一次辨识过长影响转矩平稳高频注入和 RLS 可以共用同一套电流采样但要注意注入期间基波控制器的带宽要适当降低避免电流环把高频响应当成扰动压掉。3.4 分时、分频、分工况的参数辨识时序把三条主线编成调度零速和极低速先做高频注入拿 Ld、Lq再用直流衰减拿 Rs低速段切到 RLS 辨识 Rs 和 ψf中高速段让 RLS 继续跑同时用观测器误差的慢反馈修正 ψf 温漂趋势。每个工况只更新当前可辨的参数其余参数冻结避免秩亏导致参数互相抵消。4. 把辨识器接进无传感器控制闭环辨识器本身能收敛不等于接进闭环就能用。它和电流环、无传感器观测器共享采样和电压重构时间尺度、死区补偿和抗饱和任何一处没对齐辨识结果都会偏。4.1 时间尺度分离观测器跑 kHz辨识器跑 Hz参数是慢变量位置和转速是快变量。常见做法是电流环 1–10 kHz无传感器观测器 1–5 kHz参数辨识 5–50 Hz。辨识更新太快会把电流噪声和 PWM 边带带进参数太慢又跟不上温漂。下面是一个按周期分频执行的骨架# 在电流环中断里调用N 为分频系数 N 100 # 10 kHz 电流环对应 100 Hz 辨识更新 counter 0 def current_loop_isr(i_dq, u_dq, we): global counter counter 1 if counter N: counter 0 y u_dq[1] - we * Ld * i_dq[0] - Lq * diq_dt rls.update([i_dq[1], we], y)逻辑说明只有计数器到 N 才执行一次 RLS其余中断只跑电流环和观测器。参数说明N 根据电流环频率和期望辨识带宽确定10 kHz 电流环取 100 对应 100 Hz但实际参数带宽还要受遗忘因子限制通常有效带宽在 5–50 Hz。4.2 电压重构、死区补偿与采样对齐在线辨识需要知道逆变器实际输出的电压。没有直接电压采样时用参考电压加重构补偿。死区误差写成 v_err sign(i) × (Td / Ts) × Vdc补偿代码如下def deadtime_compensation(i_abc, v_ref, Td, Ts, Vdc): # i_abc: 三相电流v_ref: 三相参考电压Td: 死区时间Ts: PWM 周期 sign np.sign(i_abc) v_comp sign * (Td / Ts) * Vdc return v_ref v_comp逻辑说明根据电流方向判断死区造成的电压损失方向再按占空比折算成电压补偿量。参数说明Td 常见 0.5–2 μsTs 对应 8–20 kHz 开关频率Vdc 用母线电压采样。不补死区低速时电压误差在参考电压里占比很大RLS 会把死区误差算进 Rs 和 ψf辨识结果整体偏移。4.3 一个可复现的仿真骨架RLS 辨识加磁链观测器把上面几块串起来下面是一个最小仿真循环用来在离线环境验证辨识器和观测器的配合import numpy as np Rs_true, psi_true 0.35, 0.12 # 真实参数 rls RLS(theta0[0.20, 0.08]) # 初始参数故意偏 for k in range(20000): we 30.0 # 低速工况rad/s iq 2.0 * np.sin(0.01 * k) # 模拟 q 轴电流 id_ 0.0 diq 2.0 * 0.01 * np.cos(0.01 * k) uq Rs_true * iq we * psi_true 0.5 * diq # 每 100 拍更新一次 if k % 100 0: y uq - we * 0.002 * id_ - 0.005 * diq theta rls.update([iq, we], y) # 用辨识参数构造磁链观测器观察位置误差趋势逻辑说明用真实参数生成 q 轴电压再故意给 RLS 一个偏低的初始 Rs 和 ψf观察它是否收敛回真实值。参数说明we 取 30 rad/s 对应低速iq 幅值 2 Adiq 用差分近似。实际调试时把这段换成实时采样数据并把辨识参数喂给磁链观测器看位置误差是否随参数收敛而减小。4.4 整定参数表与失败排查参数建议范围作用异常表现遗忘因子 λ0.995–0.999追踪速度与噪声权衡太小则参数抖动大P01e2–1e4初始收敛速度太大则初期参数跳变注入频率500 Hz–2 kHzLd/Lq 辨识接近电流环带宽会振荡注入幅值5%–10% 额定电压信噪比过大引起转矩脉动EKF 参数 Q1e-10–1e-8参数变化速率太大则参数漂EKF 电流 R1e-4–1e-2电流噪声抑制太小则滤波发散排查顺序通常是先看电压重构和死区补偿是否生效再看辨识触发工况是否满足可辨识条件最后才调遗忘因子和协方差。5. 温漂、饱和与在线辨识的验证手法5.1 温度与饱和下的参数漂移实测冷态到热态的 Rs 变化最容易测停机注入直流测电压电流比间隔记录。ψf 的温漂不好直接测常用反电动势法在中高速空载下反推。Ld、Lq 的饱和曲线用高频注入扫电流工作点得到注意交叉饱和——q 轴电流会改变 d 轴电感所以扫点时要固定另一个轴的电流。5.2 验证辨识结果的三个动作验证动作做法合格判据阶跃扰动人为改变母线电压或负载看参数是否朝物理方向收敛收敛方向与温升、饱和一致离线比对用 LCR 表或堵转测试得到 Rs、Ld、Lq在线值与离线值偏差小于 10%交叉验证用辨识参数重算观测器误差位置误差峰峰值下降第三个动作最直接把辨识前后的参数分别代入同一个磁链观测器比较低速段的位置误差峰峰值。如果辨识后误差反而变大说明辨识器被死区误差或采样偏置带偏了先退回去查电压重构。5.3 低速用高频注入、中高速用 RLS 的混合策略工程上很少只用一种辨识方法。低速段反电动势弱RLS 对 ψf 不可辨但高频注入对 Ld、Lq 有效中高速段反电动势足够RLS 对 Rs 和 ψf 的辨识条件好但高频注入会带来额外损耗和噪声。常见的混合做法是按转速加权def blend_weight(we, we_low20.0, we_high80.0): # 低速权重高中高速权重低过渡区线性切换 if we we_low: return 1.0 if we we_high: return 0.0 return (we_high - we) / (we_high - we_low)逻辑说明权重为 1 时完全用高频注入结果更新 Ld、Lq为 0 时切到 RLS 更新 Rs、ψf过渡区按权重融合两路结果并用一阶滞后平滑避免切换点参数跳变。参数说明we_low 和 we_high 按电机额定转速的 8%–12% 和 30%–40% 设置融合时间常数取 0.1–0.5 s。把切换点设在观测器换挡区之外能避开低速观测器和中高速观测器切换时的位置跳变。本文还有配套的精品资源点击获取
