CN113624253A - 一种对三轴磁传感器的旋转体误差补偿及实验方法 - Google Patents
一种对三轴磁传感器的旋转体误差补偿及实验方法 Download PDFInfo
- Publication number
- CN113624253A CN113624253A CN202110840487.4A CN202110840487A CN113624253A CN 113624253 A CN113624253 A CN 113624253A CN 202110840487 A CN202110840487 A CN 202110840487A CN 113624253 A CN113624253 A CN 113624253A
- Authority
- CN
- China
- Prior art keywords
- error
- geomagnetic
- rotator
- rotating body
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Abstract
本发明公开了一种对三轴磁传感器的旋转体误差补偿及实验方法,根据旋转体实际所处环境分析误差来源以及建立包含温度与环境影响的误差模型;采用椭球拟合对旋转体环境误差进行补偿;采用无磁转台以及平面旋转采点的方法对实际地磁传感器进行补偿验证。本发明有效的解决了环境磁场以及安装位置对传感器输出的影响,可以实现在磁环境干扰下对三轴磁传感器的旋转体误差标定。
Description
技术领域
本发明涉及一种对含误差数据进行处理的方法,具体涉及一种对三轴磁传感器的旋转体误差补偿及实验方法。
背景技术
现在的智能导航、飞行体的制导控制等项目越来越多利用地磁来实现姿态测量、获取旋转体姿态信息。地磁场属于弱磁场,磁传感器极易受到外部磁场的干扰,且磁传感器本身也存在着工艺、安装等误差。为了得到可靠的,高精度的导航信息,必须对磁传感器进行有效的误差标定与补偿。椭球拟合标定方法和最小二乘法因其简单、易实现的特点可在实际工程中应用。
发明内容
本发明的目的在于提供一种对三轴磁传感器的旋转体误差补偿及实验方法,基于工程应用的角度,能够实现消除安装环境及制造工艺对磁传感器的干扰,降低测量误差。
实现本发明目的的技术方案为:一种对三轴磁传感器的旋转体误差补偿及实验方法,根据旋转体实际所处环境分析误差来源以及建立包含温度与环境影响的误差模型;采用椭球拟合对旋转体环境误差进行补偿;采用无磁转台以及平面旋转采点的方法对实际地磁传感器进行补偿验证。
进一步的,采用两个两轴的传感器正交放置来采集三轴地磁数据,磁传感器的误差主要来源自四个方面:旋转体环境误差、不正交误差、灵敏度误差、零漂误差。
进一步的,误差模型表现为:
[Hx(T) Hy(T) Hz(T)]T是三轴地磁理论输出值,是有环境干扰下并受传感器自身影响的地磁输出值;旋转体系数矩阵,考虑了传感器的灵敏度误差与安装误差,G=[gx0,gy0,gz0]T为旋转体偏置矩阵,体现的是零位误差。
进一步的,采用二次多项式拟合对温度影响进行补偿,采用基于最小二乘法的椭球拟合对旋转体环境和传感器自身误差进行补偿。
进一步的,将三轴磁传感器的理论输出值由于制造误差的存在,原本构成的球面畸变成原点平移的椭球面。而地球磁场在短时间内可视为一个恒定磁场,因此地磁场强度H可视为一常量,从而有HTH=(H′-G)TM(H′-G),其中M=ATA,HTH为常数,A为旋转体系数矩阵,G是旋转体偏置矩阵,
进一步的,采集数据时,将三脚架中心点对准初始设定的地磁测量原点,无磁转台水平放置在三脚架上,调整三脚架每个脚的伸缩程度保证旋转体支撑平面平行于水平面,三脚架和转动平台均由无磁材料制成。用四个滚珠轴承支撑旋转体,将装有三轴磁强计的组件安装在头部,外接数据线与采集设备连接。转台和旋转装置可以分别模拟组件的俯仰、偏航和滚转方向的角度变化。
进一步的,采集地磁数据时,无磁转台在水平面内从磁北方向开始顺时针均匀转动8个45°度,在每个角度上分别顺、逆时针转动旋转体一分钟。验证标定方法时,上位机将计算得出的旋转体补偿系数发送给旋转体,旋转体经过嵌入式计算机解算后输出实时采集的补偿后的三轴地磁值,上位机进行采集后验证补偿后地磁值的正确性。
与现有技术相比,本发明的有益效果为:本发明对三轴磁传感器在工程实际中产生的误差进行分析建模,采用椭球拟合方法对三轴磁传感器的输出误差进行补偿,并提出了一种标准的旋转体标定环境及方法,快速、有效的解决了在外界环境干扰与自身误差影响下地磁数据发散的问题,最大程度上保障了旋转体上地磁数据的准确性和可靠性。
附图说明
图1为本发明的示意流程图。
图2为本发明的旋转体标定环境示意图。
图3为本发明三轴磁传感器实际输出地磁值三维图。
图4为本发明预处理后三轴地磁值三维散点图。
图5为本发明椭球拟合补偿方法流程图。
图6为本发明误差补偿后旋转体指北三轴地磁值。
图7为本发明误差补偿后旋转体指东三轴地磁值。
图8为本发明误差补偿后三轴地磁值三维散点图。
图中:1、HMC1052磁强计,2、+5V稳压直流电源,3、STM32F7处理器,4、旋转体对外通信接口FGG.0b,5、MOXA卡串口线。
具体实施方式
本发明针对在工程应用中三轴磁传感器进行标定时操作不便、标定困难等问题,在考虑了温度对磁传感器影响的前提下,公开了一种简单有效的去除误差的标定方法。该方法分析了磁传感器的误差来源,建立了包含安装位置、载体环境等影响因素的误差模型,采用最小二乘法原理对磁传感器的输出进行了椭球拟合补偿以及补偿实验方式,有效的解决了环境磁场以及安装位置对传感器输出的影响,可以实现在磁环境干扰下对三轴磁传感器的旋转体误差标定。
如图1所示,一种三轴磁传感器的旋转体误差补偿及实验方法,实现步骤包括:采集数据、数据预处理、平滑滤波、椭球拟合补偿、数据仿真验证。
采集数据是将组件水平固定放置,保证在旋转体旋转和转动过程中磁传感器相对旋转体保持静止状态。旋转体平台在同一温度前提下,平稳水平转动一周360°,并对多个角度进行旋转体旋转采集地磁数据;
数据预处理是将采集到的实际地磁值按不同角度进行分类,并进行去跳点处理;
滤波是对各个角度的地磁数据进行平滑处理,根据椭球所需筛选出合理的数据;
椭球拟合补偿模是对筛选过后的地磁值根据椭球公式进行拟合近似,利用最小二乘法原理计算出椭球系数与偏置系数,最终将计算得出的补偿系数保存记录到本地;
数据仿真是根据椭球拟合补偿模块计算出的补偿系数与误差补偿模型,带入地磁传感器输出的原始值计算仿真,验证椭球补偿算法的正确性。
进一步的,组件旋转体标定环境包括稳压直流电源、MOXA串口线缆、某型号旋转体、无磁三脚架、无磁转台,如图2所示。旋转体标定方法是将包含三轴磁传感器的组件安装在旋转体头部,外接数据线与采集设备连接,转动无磁转台,在同一平面内采集不同角度下的椭球地磁值。
进一步的,所述三轴磁传感器由两个两轴磁传感器正交放置构成;磁传感器型号为HMC1052。
进一步的,滤波模块采用10阶滑动滤波。
进一步的,椭球拟合补偿模块采用的椭球公式为:
f(x,y,z)=x2+by2+cz2+2fxy+2gxz+2hyz+2px+2qy+2rz+d=0
令ξ=[b,c,f,g,h,p,q,r,d]T,其中b,c表示待求椭球的曲面的平方项系数,f,g,h表示二次项系数,p,q,r表示一次项系数,d为常数值;
Hk=[xk,yk,zk]T,k=1,2,...,n,表示k组三轴磁传感器实际输出值;
所述椭球拟合补偿模块改写椭球公式为:
HTEH+2FTH+G=0
可将椭球方程改成:
(H-H0)TK(H-H0)=1
将矩阵M进行奇异值分解,得到M=a2U∑UT
式中a∈R,∑=Λ(λ1,λ2,λ3)为特征值λ1,λ2,λ3构成的对角阵,U为特征向量矩阵。则标定系数矩阵A=U∑UT。
所述椭球拟合补偿模块最终标定方程为:
一种对三轴磁传感器的旋转体误差补偿与实验方法具体为:
旋转体标定环境配置好,将无磁转台水平指北放置;
无磁转台顺时针每转动45°固定一次,在水平面内共转动8个角度,固定每个角度下的转台并将旋转体顺时针转动旋转体至少三周后逆时针转动至少三周,采集地磁数据;
按转台角度分别读取地磁值并选择在旋转体旋转后半段时采集的数据,使采集的地磁值满足椭球标准;
初始化滑动滤波系数为10;
初始化当地地磁场强度,根据奇异值分解求得特征值λ和特征向量矩阵U,从而求得标定系数矩阵A与偏置矩阵H0;
根据最终标定方程用标定系数矩阵A与偏置矩阵H0计算标定后的三轴地磁值。
下面结合说明书附图对本发明一种三轴磁传感器的旋转体误差补偿及实验方法作进一步详细说明。
实施例
如图2所示,一种三轴磁传感器的旋转体误差补偿实验方法为:将组件旋转体水平指北放置于无磁转台上,转台顺时针每转动45°角为一个测量点,在每个测量点上沿旋转体x轴匀速转动旋转体,分别正反转动3圈以上,使地磁原始值的斜置轴y与纵轴z的地磁值均满足平滑的正弦波形,相位相差
嵌入式处理器采用STM32F7芯片,误差补偿实验之前上位机给处理器STM32F7初始化当地地磁值、温度系数矩阵C、温度零偏矩阵B;
实验时+5V直流稳压电源给组件供电,旋转体尾部通过9针的FGG.0b雷莫接口与MOXA UPORT1150串口线缆连接,实时将三轴磁传感器原始输出值通过串口发送给上位机;
上位机通过椭球拟合补偿将最终的标定系数通过串口与FGG.0b雷莫接口发送给STM32F7,STM32F7根据误差模型将地磁原始值与标定系数带入最终标定方程计算,得到误差补偿后的地磁数据,并将该值通过串口传送给上位机监控整个传输过程采用IEEE754协议;
上位机接收到数据后,首先对各个角度进行分类,形成8个通道数据,并根据传输协议取出三轴地磁原始值,如图3所示,三轴地磁原始值在空间上表现为一个完整椭球的各个椭圆横截面,最大的椭圆在旋转体指东和指西时获得;
基于上述装置及实验方法,一种对三轴磁传感器的旋转体误差补偿方法采用滤波与基于最小二乘法的椭球拟合的形式。
在获得原始地磁数据后,对各通道地磁值进行去跳点处理,排除人工转动旋转体时可能出现的旋转体滑动对量测地磁值的阶跃影响。采用10阶的滑动滤波进一步减小旋转体转动时人工与外界对旋转体产生的微小扰动和滑动。滤波后的数据选择旋转体旋转的后半段来确保拟合的地磁数据可以稳定的构成平滑椭球的一部分,可以进一步减小标定误差,如图4所示,地磁在经过滤波和筛选后,分散点呈现清晰的椭球分布;
图5为整个误差补偿算法的流程图。误差模型中的C矩阵,在保持标定过程中旋转体温度变化很小的情况下设置为单位阵,偏置矩阵B仅与旋转体环境及安装误差相关,与温度无关。椭球方程中旋转轴x的二次项系数设置为1;因为短时间内的地磁场可视为一常数,所以通过球谐计算设置地磁场强度为48375nT。将滤波与预处理后的各通道地磁值拼接成x、y、z坐标矩阵组成椭球截面,带入计算,由最小二乘原理解的椭球系数矩阵E、F;
将矩阵M进行奇异值分解,得到M的特征值λ与特征向量U,带入标定方程求解出误差补偿系数A矩阵与偏置矩阵H0。
旋转体在指北和指东两个角度下误差补偿后的地磁数据如图6、图7所示,可以看出在经过误差补偿后的地磁值峰峰值与零点偏置值十分接近,具体计算方法为:
将误差补偿系数矩阵A通过串口线发送至嵌入式计算器STM32F7,ARM芯片将椭球拟合得到的A、H0矩阵与初始化时的C、B矩阵带入标定方程计算得出当前补偿后的地磁值H。
补偿后地磁数据通过串口线发送到上位机进行监控,数据三维散点图如图8所示,每个位置上的三轴地磁数据呈现出较为标准的圆形,地磁补偿效果较好。
Claims (8)
1.一种对三轴磁传感器的旋转体误差补偿及实验方法,其特征在于,根据旋转体实际所处环境分析误差来源以及建立包含温度与环境影响的误差模型;采用椭球拟合对旋转体环境误差进行补偿;采用无磁转台以及平面旋转采点的方法对实际地磁传感器进行补偿验证。
2.根据权利要求1所述的对三轴磁传感器的旋转体误差补偿及实验方法,其特征在于,采用两个两轴的传感器正交放置来采集三轴地磁数据,磁传感器的误差主要来源自四个方面:旋转体环境误差、不正交误差、灵敏度误差、零漂误差。
5.根据权利要求3所述的对三轴磁传感器的旋转体误差补偿及实验方法,其特征在于,采用二次多项式拟合对温度影响进行补偿,采用基于最小二乘法的椭球拟合对旋转体环境和传感器自身误差进行补偿。
7.根据权利要求1所述的对三轴磁传感器的旋转体误差补偿及实验方法,其特征在于,采集数据时,将三脚架中心点对准初始设定的地磁测量原点,无磁转台水平放置在三脚架上,调整三脚架每个脚的伸缩程度保证旋转体支撑平面平行于水平面,三脚架和转动平台均由无磁材料制成;用四个滚珠轴承支撑旋转体,将装有三轴磁强计的组件安装在头部,外接数据线与采集设备连接;转台和旋转装置分别模拟组件的俯仰、偏航和滚转方向的角度变化。
8.根据权利要求7所述的对三轴磁传感器的旋转体误差补偿及实验方法,其特征在于,采集地磁数据时,无磁转台在水平面内从磁北方向开始顺时针均匀转动8个45°度,在每个角度上分别顺、逆时针转动旋转体一分钟;验证标定方法时,上位机将计算得出的旋转体补偿系数发送给旋转体,旋转体经过嵌入式计算机解算后输出实时采集的补偿后的三轴地磁值,上位机进行采集后验证补偿后地磁值的正确性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110840487.4A CN113624253A (zh) | 2021-07-24 | 2021-07-24 | 一种对三轴磁传感器的旋转体误差补偿及实验方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110840487.4A CN113624253A (zh) | 2021-07-24 | 2021-07-24 | 一种对三轴磁传感器的旋转体误差补偿及实验方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113624253A true CN113624253A (zh) | 2021-11-09 |
Family
ID=78380818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110840487.4A Pending CN113624253A (zh) | 2021-07-24 | 2021-07-24 | 一种对三轴磁传感器的旋转体误差补偿及实验方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113624253A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115468485A (zh) * | 2022-11-15 | 2022-12-13 | 国仪量子(合肥)技术有限公司 | 工具面角的计算方法及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6014610A (en) * | 1997-01-31 | 2000-01-11 | Greenfield Enterprises, Inc | Navigation system and method |
CN103630137A (zh) * | 2013-12-02 | 2014-03-12 | 东南大学 | 一种用于导航系统的姿态及航向角的校正方法 |
CN108761548A (zh) * | 2018-05-25 | 2018-11-06 | 哈尔滨工程大学 | 条件约束最小二乘椭球拟合的微型磁力计磁干扰补偿算法 |
CN111780786A (zh) * | 2020-08-08 | 2020-10-16 | 武汉利科夫科技有限公司 | 一种三轴tmr传感器在线校准方法 |
-
2021
- 2021-07-24 CN CN202110840487.4A patent/CN113624253A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6014610A (en) * | 1997-01-31 | 2000-01-11 | Greenfield Enterprises, Inc | Navigation system and method |
CN103630137A (zh) * | 2013-12-02 | 2014-03-12 | 东南大学 | 一种用于导航系统的姿态及航向角的校正方法 |
CN108761548A (zh) * | 2018-05-25 | 2018-11-06 | 哈尔滨工程大学 | 条件约束最小二乘椭球拟合的微型磁力计磁干扰补偿算法 |
CN111780786A (zh) * | 2020-08-08 | 2020-10-16 | 武汉利科夫科技有限公司 | 一种三轴tmr传感器在线校准方法 |
Non-Patent Citations (1)
Title |
---|
安亮亮等: "弹载三轴磁强计的椭球拟合补偿方法", 《传感器与微系统》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115468485A (zh) * | 2022-11-15 | 2022-12-13 | 国仪量子(合肥)技术有限公司 | 工具面角的计算方法及存储介质 |
CN115468485B (zh) * | 2022-11-15 | 2023-03-10 | 国仪量子(合肥)技术有限公司 | 工具面角的计算方法及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107024674B (zh) | 一种基于递推最小二乘法的磁强计现场快速标定方法 | |
CN110849403B (zh) | 一种定向传感器的标定方法 | |
CN103153790B (zh) | 使用运动传感器和附接至装置的磁力计的测量数据估计该装置在重力参照系中的偏航角的设备和方法 | |
CN106125026B (zh) | 一种不依赖于地磁场场量的三轴磁强计全误差参数辨识与校正方法 | |
CN107544042B (zh) | 一种磁力计阵列校正方法 | |
CN103591949B (zh) | 三轴姿态测量系统非正交性误差的正交补偿方法 | |
CN105806364B (zh) | 一种矿用回转钻机测斜仪探管的校准方法 | |
CN109974749B (zh) | 一种三轴转台综合指向误差的评定方法 | |
CN110146839A (zh) | 一种移动平台磁梯度张量系统校正方法 | |
CN104198973B (zh) | 一种矢量磁力仪校准装置 | |
CN106997035A (zh) | 一种基于磁梯度不变量的磁梯度计校正方法 | |
CN107121707A (zh) | 一种三轴磁传感器测量基准与结构基准的误差校正方法 | |
Gang et al. | Linear calibration method of magnetic gradient tensor system | |
Zongwei et al. | A low-cost calibration strategy for measurement-while-drilling system | |
CN111413750B (zh) | 一种用于电缆定位的磁通门传感器的误差校正方法 | |
Gheorghe et al. | Calibration optimization study for tilt-compensated compasses | |
CN113624253A (zh) | 一种对三轴磁传感器的旋转体误差补偿及实验方法 | |
Wu et al. | Simultaneous hand–eye/robot–world/camera–IMU calibration | |
CN113866688B (zh) | 一种小姿态角条件下的三轴磁传感器误差校准方法 | |
CN108088431B (zh) | 一种自校正电子罗盘及其校正方法 | |
Gao et al. | A calibration method for the misalignment error between inertial navigation system and tri-axial magnetometer in three-component magnetic measurement system | |
Huang et al. | Two-step complete calibration of magnetic vector gradiometer based on functional link artificial neural network and least squares | |
Pang et al. | A new misalignment calibration method of portable geomagnetic field vector measurement system | |
Liu et al. | A hybrid calibration method for a three-axis magnetometer in limited-range attitudes | |
CN110398702A (zh) | 一种基于多传感器融合的实时在线磁校准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211109 |
|
RJ01 | Rejection of invention patent application after publication |