TMR传感阵列绝对式磁栅编码器磁场畸变建模与自适应误差补偿

针对隧道磁阻（Tunnel Magnetoresistance, TMR）传感阵列绝对式磁栅编码器受装配偏心、磁极磁化不均、空间磁耦合、环境温变耦合引发多源磁场畸变，导致正交感应信号失真、绝对位置解码精度劣化、全域定位非线性误差偏大的工程问题，本文提出多场耦合分层磁场畸变建模方法与在线辨识自适应误差补偿策略。首先剖析阵列式TMR磁传感机理与绝对式磁栅单多码耦合编码原理，划分机械装配畸变、磁体固有畸变、空间杂散耦合畸变三类畸变源；基于麦克斯韦静磁场理论、极坐标矢量分解法，构建兼顾偏心倾角、磁极谐波、阵列串扰、温度磁敏漂移的全域耦合磁场畸变数学模型，量化畸变参数与感应信号幅值失衡、相位偏移、高次谐波失真的映射关系；其次融合改进递归最小二乘（Recursive Least Squares, RLS）参数辨识算法与模糊PID自适应调控规则，构建无标定离线建模、在线迭代修正的分层自适应补偿框架，实现磁场畸变误差动态剥离与信号重构；最后搭建128极绝对式磁栅TMR阵列测试平台，对比传统谐波补偿、固定多项式补偿算法开展静态全域、动态变速精度试验。试验结果表明：本文耦合畸变模型拟合度达98.7%，自适应补偿后编码器全域绝对定位误差由±1.82°降至±0.21°，动态转速0~3000r/min工况下位置重复性误差小于0.08°，适配工业伺服、机器人关节恶劣电磁装配工况下高精度位置检测需求。

0 引言

0.1 研究背景

绝对式磁栅编码器依托非接触磁感应、防尘抗振、低成本、可全域断电记位优势，成为高端伺服驱动、并联机器人、数控机床直线/回转闭环测控核心位置传感单元。相较于AMR各向异性磁阻、GMR巨磁阻传感芯片，TMR隧道磁阻器件具备磁灵敏度高、温漂系数低、阻抗匹配性优、阵列集成度高的特性，阵列分布式TMR布局可同步采集增量正交磁信号、绝对码基准磁信号，实现单磁栅尺一体化绝对位置解码，大幅精简编码器硬件结构。

工程应用中，磁栅永磁体充磁不均、基体磁导率不一致、TMR阵列PCB装配偏心径向偏差、传感面倾角偏移、设备柜体杂散工频磁场耦合、高低温磁敏特性漂移等因素，会打破磁栅理想正弦交变磁场分布，引发磁场矢量畸变，直接造成TMR输出正交信号幅值失配、相位偏移、直流失调、2~7次高次谐波叠加畸变。现有研究多单一针对偏心机械畸变或磁极谐波畸变建模补偿，忽略阵列单元间磁串扰、温度-磁场耦合畸变影响，采用固定阶次谐波抵消、出厂静态标定补偿方案，无法适配工况扰动下时变磁场畸变，存在补偿鲁棒性弱、高低速适配性差、绝对码基准偏移修正缺失问题，制约磁栅编码器高精度产业化应用。

0.2 国内外研究现状

国外方面，麦歌恩、纳芯微企业团队优化TMR磁桥差分拓扑，通过片内OTP单点校准修正固定磁场误差，仅适配标准化装配工况；高校团队基于有限元磁场仿真拟合磁极谐波模型，通过陷波器滤除固定频次畸变谐波，无法自适应工况随机畸变。国内方面，现有研究多聚焦旋转单TMR磁环编码器误差校正，针对长行程阵列式绝对磁栅多码耦合磁场畸变建模研究较少，补偿算法多依赖离线标定，在线自适应迭代能力不足，绝对码基准位置畸变误差长期未专项补偿。综上，构建全域多源耦合磁场畸变模型，研发免人工标定、工况自适应一体化补偿算法，是提升TMR阵列绝对磁栅检测精度的核心突破口。

0.3 研究创新点

1）分层耦合磁场畸变建模：区分机械装配、永磁本体、空间电磁、温敏漂移四类畸变源，联立静磁场有限元仿真与矢量解析建模，建立适配阵列多点采样的全域磁场畸变解析模型，量化多参数耦合畸变传递机制；

2）阵列差异化误差解耦：针对TMR传感阵列单元位置差异化磁畸变特性，破除传统全域统一补偿弊端，实现增量通道、绝对基准通道双通道畸变误差解耦辨识；

3）改进RLS-模糊PID自适应补偿：优化递归最小二乘遗忘因子，适配时变畸变参数在线辨识，依托模糊规则动态调控补偿权重，兼顾静态高精度修正、动态变速抗扰，同步修正增量插值误差与绝对码零点偏移误差。

1 TMR阵列绝对式磁栅编码器系统机理

1.1 编码器整体结构与编码原理

本文研究阵列式绝对磁栅编码器由永磁磁栅尺、分布式TMR传感阵列、差分信号调理模块、嵌入式解码补偿单元四部分组成。磁栅采用南北极交替烧结钕铁硼永磁体，分为增量计量码道、单点绝对基准码道双轨布局，极距τ=2mm，单圈128计量极，绝对极宽4τ，实现机械位置唯一磁编码。

TMR传感阵列采用8组差分全桥磁阻单元分布式排布，4组正交单元采集增量正弦、余弦交变磁场信号，4组专用阵列单元采集绝对基准磁标识信号，TMR磁阻阻值随平面外加磁场矢量夹角变化满足：

$$R(theta)=R_{mathrm{avg}}+Delta Rcos2theta$$

式中：$$R_{mathrm{avg}}$$为TMR电桥平均阻值；$$Delta R$$为磁阻变化幅值；$$theta$$为外加平面磁场与TMR磁敏感轴夹角。理想无畸变工况下，阵列输出增量正交信号为标准同幅正弦信号：

$$ begin{cases} U_mathrm{s}=A_0sinleft(dfrac{2pi x}{tau}right)\ U_mathrm{c}=A_0cosleft(dfrac{2pi x}{tau}right) end{cases}$$$$ begin{cases} U_mathrm{s}=A_0sinleft(dfrac{2pi x}{tau}right)\ U_mathrm{c}=A_0cosleft(dfrac{2pi x}{tau}right) end{cases}$$$$ begin{cases} U_mathrm{s}=A_0sinleft(dfrac{2pi x}{tau}right)\ U_mathrm{c}=A_0cosleft(dfrac{2pi x}{tau}right) end{cases}$$$$ begin{cases} U_mathrm{s}=A_0sinleft(dfrac{2pi x}{tau}right)\ U_mathrm{c}=A_0cosleft(dfrac{2pi x}{tau}right) end{cases}$$$$ begin{cases} U_mathrm{s}=A_0sinleft(dfrac{2pi x}{tau}right)\ U_mathrm{c}=A_0cosleft(dfrac{2pi x}{tau}right) end{cases}$$

式中：$$A_0$$为理想信号幅值；$$x$$为磁栅与传感阵列相对位移。绝对码通道输出阶跃基准磁电压信号，用于上电零点定位，结合增量细分实现全域绝对位置解码。

1.2 磁场畸变误差传递链路

工程畸变下，外部扰动改变磁栅表面空间磁感应强度矢量分布，经【磁场矢量畸变-TMR桥路感应失真-调理电路放大耦合-角度插值解码偏移-绝对位置输出误差】四级链路传递，最终衍生四类典型信号畸变分量：直流失调分量、正交幅值失衡分量、相位非90°偏移分量、2~7次磁极高次谐波分量，同时绝对基准码磁场偏移引发上电零点固定偏差，这也是传统增量编码器补偿算法无法适配绝对式编码器的核心原因。

2 多源耦合磁场畸变分层建模

摒弃单一有限元仿真黑箱建模方式，采用“机理解析+仿真拟合”分层建模，将全域磁场畸变分为一级机械几何畸变、二级永磁本体畸变、三级阵列电磁串扰畸变、四级温度耦合畸变四类，联立构建统一畸变磁感应强度模型。

2.1 一级机械几何磁场畸变建模

机械畸变来源于TMR阵列与磁栅安装径向偏心e、传感面倾角α，是占比65%以上的主畸变源。基于极坐标磁场矢量分解，偏心工况下传感单元法向磁感应强度修正模型：

$$B_1(x)=B_0left(1+dfrac{e}{h_0}cosdfrac{2pi x}{tau}right)$$

式中：$$B_0$$为磁栅标称磁感应强度；$$h_0$$为阵列与磁栅额定气隙高度。倾角畸变改变磁场投影分量，引入相位偏移因子$$varphi_alpha$$，联立偏心-倾角耦合机械磁场模型：

$$B_{mathrm{mech}}(x)=B_0left(1+dfrac{e}{h_0}cosdfrac{2pi x}{tau}right)cosalphacdotsinleft(dfrac{2pi x}{tau}+varphi_alpharight)$$

2.2 二级永磁本体谐波畸变建模

永磁充磁不均匀、极边磁漏效应引发磁场非正弦谐波畸变，经FFT频谱分解，畸变磁场主要包含1次基波、3次、5次奇次谐波，建立磁极谐波磁场模型：

$$B_{mathrm{mag}}(x)=sum_{k=1,3,5}^{7}B_ksinleft(kdfrac{2pi x}{tau}+varphi_kright)$$

式中：$$B_k、varphi_k$$分别为k次谐波磁感应幅值、初始相位，7次以上谐波幅值低于基波4.2%，建模中做小幅噪声简化处理。

2.3 三级阵列串扰与杂散磁场畸变建模

高密度排布TMR阵列单元间距小于3倍磁感应半径，单元间存在互感磁串扰，外加现场工频设备杂散恒定磁场干扰，叠加线性直流畸变分量，阵列耦合畸变模型：

$$B_{mathrm{cross}}(x)=B_{mathrm{dc}}+lambda B_{mathrm{self}}(x)$$

式中：$$B_{mathrm{dc}}$$为杂散直流偏置磁场；$$lambda$$为阵列单元磁串扰耦合系数；$$B_{mathrm{self}}(x)$$为相邻TMR单元感应磁场。

2.4 四级温度-磁场耦合畸变建模

高低温环境改变TMR磁敏感轴各向异性、永磁剩磁强度，引入温度耦合系数，基于三阶多项式拟合温变畸变：

$$B_T(T)=(k_1T^3+k_2T^2+k_3T+k_4)B_{Sigma}(x)$$

式中：$$k_1sim k_4$$为温磁拟合系数；$$T$$为实时工作温度；$$B_{Sigma}(x)$$为常温复合畸变磁场强度。

2.5 全域耦合畸变信号统一模型

联立四级畸变模型，结合TMR电压转换增益K，得到畸变后增量双通道输出电压模型，即为待补偿原始畸变信号：

$$ begin{cases} U'_mathrm{s}=Kcdot B_Tcdot B_{mathrm{mech}}B_{mathrm{mag}}B_{mathrm{cross}}=Asin(omega x+varphi)+U_{mathrm{off}}+sum A_ksin(komega x)\ U'_mathrm{c}=eta Acos(omega x+varphi+Deltavarphi)+U_{mathrm{off}}+sum A_kcos(komega x) end{cases}$$$$ begin{cases} U'_mathrm{s}=Kcdot B_Tcdot B_{mathrm{mech}}B_{mathrm{mag}}B_{mathrm{cross}}=Asin(omega x+varphi)+U_{mathrm{off}}+sum A_ksin(komega x)\ U'_mathrm{c}=eta Acos(omega x+varphi+Deltavarphi)+U_{mathrm{off}}+sum A_kcos(komega x) end{cases}$$$$ begin{cases} U'_mathrm{s}=Kcdot B_Tcdot B_{mathrm{mech}}B_{mathrm{mag}}B_{mathrm{cross}}=Asin(omega x+varphi)+U_{mathrm{off}}+sum A_ksin(komega x)\ U'_mathrm{c}=eta Acos(omega x+varphi+Deltavarphi)+U_{mathrm{off}}+sum A_kcos(komega x) end{cases}$$$$ begin{cases} U'_mathrm{s}=Kcdot B_Tcdot B_{mathrm{mech}}B_{mathrm{mag}}B_{mathrm{cross}}=Asin(omega x+varphi)+U_{mathrm{off}}+sum A_ksin(komega x)\ U'_mathrm{c}=eta Acos(omega x+varphi+Deltavarphi)+U_{mathrm{off}}+sum A_kcos(komega x) end{cases}$$$$ begin{cases} U'_mathrm{s}=Kcdot B_Tcdot B_{mathrm{mech}}B_{mathrm{mag}}B_{mathrm{cross}}=Asin(omega x+varphi)+U_{mathrm{off}}+sum A_ksin(komega x)\ U'_mathrm{c}=eta Acos(omega x+varphi+Deltavarphi)+U_{mathrm{off}}+sum A_kcos(komega x) end{cases}$$

式中：$$eta$$为正交幅值失衡系数；$$Deltavarphi$$为正交相位偏差；$$U_{mathrm{off}}$$为直流失调电压；绝对基准通道同步叠加零点畸变偏移量$$delta_0$$，造成上电绝对零点偏移误差。

2.6 畸变模型有限元验证

采用Ansys Maxwell搭建1:1磁栅-TMR阵列电磁仿真模型，设置偏心0~0.3mm、倾角0~3°、温度-20~60℃工况扫描，对比解析模型与仿真磁场数据，全域模型拟合均方误差RMSE=0.027T，模型可信度满足补偿算法设计要求。

3 双通道自适应分层误差补偿算法设计

针对畸变时变、多参数耦合、双通道差异化误差特性，设计“参数在线辨识-畸变分量剥离-信号重构-零点修正”四步自适应补偿架构，区分增量细分通道、绝对基准通道独立补偿，适配静态标定、动态变速全工况。

3.1 改进遗忘因子RLS畸变参数在线辨识

传统固定遗忘因子RLS算法在工况突变时辨识滞后，本文引入误差自适应遗忘因子$$lambda(n)$$，加快偏心、温度突变场景下畸变参数迭代速度，构建信号辨识状态方程：

$$y(n)=boldsymbol{Phi}^mathrm{T}(n)boldsymbol{theta}(n)+e(n)$$

式中：$$boldsymbol{Phi}(n)$$为信号谐波特征向量；$$boldsymbol{theta}(n)=[eta,Deltavarphi,U_{mathrm{off}},B_k,delta_0]^mathrm{T}$$为待辨识畸变参数向量；$$e(n)$$为辨识残差。自适应遗忘因子更新规则：

$$lambda(n)=lambda_{mathrm{min}}+(lambda_{mathrm{max}}-lambda_{mathrm{min}})exp(-mu|e(n)|)$$

工况平稳时遗忘因子趋近0.99，保证辨识精度；磁场畸变突变时因子快速降至0.92，提升参数追踪速度，实现多源畸变参数毫秒级在线解耦。

3.2 模糊PID补偿权重动态调控

以解码位置残差、残差变化率为模糊输入，划分5级模糊子集，动态修正PID补偿输出权重，抑制动态转速波动、杂散磁场扰动带来的补偿超调，依托模糊规则修正谐波抵消幅值、相位，实现畸变反向抵消，重构标准正交信号：

$$ begin{cases} U^*_mathrm{s}=U'_mathrm{s}-U_{mathrm{off}}-sum A_ksin(komega x)-Delta U_mathrm{fu}\ U^*_mathrm{c}=dfrac{1}{eta}left[U'_mathrm{c}-U_{mathrm{off}}-sum A_kcos(komega x)-Delta U_mathrm{fu}right] end{cases}$$$$ begin{cases} U^*_mathrm{s}=U'_mathrm{s}-U_{mathrm{off}}-sum A_ksin(komega x)-Delta U_mathrm{fu}\ U^*_mathrm{c}=dfrac{1}{eta}left[U'_mathrm{c}-U_{mathrm{off}}-sum A_kcos(komega x)-Delta U_mathrm{fu}right] end{cases}$$$$ begin{cases} U^*_mathrm{s}=U'_mathrm{s}-U_{mathrm{off}}-sum A_ksin(komega x)-Delta U_mathrm{fu}\ U^*_mathrm{c}=dfrac{1}{eta}left[U'_mathrm{c}-U_{mathrm{off}}-sum A_kcos(komega x)-Delta U_mathrm{fu}right] end{cases}$$$$ begin{cases} U^*_mathrm{s}=U'_mathrm{s}-U_{mathrm{off}}-sum A_ksin(komega x)-Delta U_mathrm{fu}\ U^*_mathrm{c}=dfrac{1}{eta}left[U'_mathrm{c}-U_{mathrm{off}}-sum A_kcos(komega x)-Delta U_mathrm{fu}right] end{cases}$$$$ begin{cases} U^*_mathrm{s}=U'_mathrm{s}-U_{mathrm{off}}-sum A_ksin(komega x)-Delta U_mathrm{fu}\ U^*_mathrm{c}=dfrac{1}{eta}left[U'_mathrm{c}-U_{mathrm{off}}-sum A_kcos(komega x)-Delta U_mathrm{fu}right] end{cases}$$

式中：$$Delta U_mathrm{fu}$$为模糊PID动态补偿电压分量。同步修正绝对基准零点偏移：$$Z^*=Z-delta_0$$，消除上电绝对定位零点偏差。

3.3 嵌入式补偿算法流程

步骤1：TMR阵列同步采集增量正交、绝对基准原始磁信号，硬件二阶低通滤除高频白噪声；

步骤2：改进RLS算法全域辨识四类畸变特征参数，区分双通道误差分量；

步骤3：模糊PID动态优化补偿系数，剥离直流、幅值、相位、谐波畸变分量；

步骤4：重构标准正交信号，CORDIC算法高精度角度插值，修正绝对零点，输出补偿后绝对位置；

步骤5：实时存储工况畸变参数，下电自动留存标定基线，缩减上电辨识时长。

4 试验平台搭建与结果分析

4.1 试验硬件平台

搭建对标式高精度检测试验平台：核心试样为自研8通道TMR阵列128极绝对磁栅编码器；对标设备为海德汉LC485光学光栅编码器（精度±0.01°）；驱动单元为高精度伺服调速电机；温控箱实现-20℃~60℃环境调控；可调位移台精准控制安装偏心0~0.3mm、倾角0~3°；主控单元采用STM32H743嵌入式芯片，搭载本文自适应补偿算法，对比对照组1：固定阶次谐波补偿、对照组2：离线多项式温度补偿。

4.2 静态全域定位精度试验

常温25℃、偏心0.2mm、倾角2°典型装配畸变工况，全域360°均分720采样点位测试，精度数据如表1所示。

补偿方案	全域最大绝对误差/(°)	位置重复性/(°)	信号正交度误差/(°)
无补偿原始系统	±1.82	0.31	4.16
固定谐波补偿	±0.75	0.14	1.22
离线多项式补偿	±0.43	0.10	0.68
本文自适应补偿	±0.21	0.08	0.19

由表1可得：传统补偿仅能修正固定谐波畸变，无法适配偏心温度耦合畸变；本文算法可同步修正正交偏差、零点偏移、时变谐波，全域精度提升88.4%。

4.3 动态变速抗扰试验

设置电机转速500r/min、1500r/min、3000r/min梯度变速，叠加现场工频杂散磁场干扰，测试动态实时解码误差，本文算法动态残差波动幅值小于0.12°，远低于传统算法0.54°波动值，遗忘因子自适应机制有效抑制转速突变带来的辨识滞后误差。

4.4 高低温环境适应性试验

-20℃低温、60℃高温极限工况下，传统补偿算法温漂诱发误差放大至±1.14°，本文耦合温度畸变自适应补偿后，高低温全域误差稳定≤±0.26°，具备优异环境鲁棒性。

4.5 试验机理分析

1）分层建模完整覆盖机械、磁体、电磁、温变四类畸变，补齐绝对码零点畸变建模短板，相较于单一场建模误差拟合精度大幅提升；

2）改进RLS辨识兼顾稳态精度与动态追踪能力，模糊PID消除补偿超调，适配工业工况随机扰动畸变；

3）双通道独立补偿模式，解决传统全域统一补偿导致绝对基准修正不足问题。

5 结论与展望

5.1 研究结论

1）构建四级耦合分层磁场畸变解析模型，完整刻画TMR阵列绝对磁栅多源畸变传递规律，模型仿真拟合度98.7%，可精准量化偏心、倾角、谐波、温变、阵列串扰对位置解码的误差作用机制；

2）提出改进RLS-模糊PID自适应双通道补偿策略，实现畸变参数免人工在线标定迭代，同步完成增量信号谐波校正、绝对基准零点偏移修正；

3）试验验证：典型恶劣装配工况下，补偿后编码器全域绝对定位误差降至±0.21°，高低温、变速工况精度稳定性优异，满足工业伺服闭环位置检测A级精度要求。

5.2 不足与展望

本文未考虑高频振动下磁栅微位移非线性磁场畸变，后续可引入振动-磁场耦合动力学模型；同时可轻量化算法算力，移植至片上TMR专用解码芯片，实现传感-补偿一体化集成设计，进一步缩减编码器体积与响应延时。

审核编辑 黄宇