一种带噪声抑制的三轴磁传感器定标方法技术

本发明专利技术公开了一种带噪声抑制的三轴磁传感器定标方法，具体步骤是：步骤1：确定噪声上限；步骤2：对噪声方差进行估计；步骤3：计算对应的最小特征向量，采用奇异值svd分解计算最小特征向量，并代入步骤2中估计噪声方差；步骤4：利用ALS方法求得椭球参数Ce和d，进而求得C；步骤5：利用公式hE＝C-1hM-d，修正三轴传感器输出数据，即成。本发明专利技术的方法，使系统既能工作在较强的噪声干扰条件下又能保证工程可实现性，有效抑制强噪声信号，明显改善三轴磁传感器抗噪性能，提高了磁测系统的测量精度。

【技术实现步骤摘要】
一种带噪声抑制的三轴磁传感器定标方法
本专利技术属于弱磁探测
，涉及一种带噪声抑制的三轴磁传感器定标方法。
技术介绍
高精度三轴传感器以三轴磁通门传感器为主，由于在加工过程中，受到加工工艺、安装工艺和材料特性等因素的限制，三轴磁传感器本身存在零偏误差、标度系数误差、非正交误差等误差；同时，在三轴磁传感器使用过程中存在剩磁、软铁误差、硬铁误差和噪声等干扰因素，极大限制了三轴磁传感器的测量精度，导致其测量误差达到几百甚至上千纳特，使得三轴磁传感器难以满足高精度磁测的需求。目前的国内外主流研究中，误差模型一般以输出磁场的椭球假设为基础，建立以校正参数为变量、传感器输出分量和姿态为系数的方程组；实测试验经常采用旋转法或翻转法；基于试验数据，选取优化算法进行参数估计，从而得到误差模型中的各类参数。而存在高斯噪声时，若单纯的通过测量数据平均的方法去除噪声，则根据椭球假设模型求得的误差参数为真实误差量的非一致性估计，在环境噪声较大且数据长度受限的高实时性约束下，上述误差估计方法就会失效。因此，现有的三轴磁传感器的定标方法通常无法处理大噪声情况下的参数估计问题，且极少数考虑噪声影响的定标方法中需要预先得知噪声的分布特性，这在实际系统中通常难以得到。
技术实现思路
本专利技术的目的是提供一种带噪声抑制的三轴磁传感器定标方法，解决了现有技术中在噪声环境下三轴磁传感器定标不准确的问题，尤其是噪声较大且噪声特性未知情况下的误差参数估计失效的问题。本专利技术所采用的技术方案是，一种带噪声抑制的三轴磁传感器定标方法，具体按照以下步骤实施：步骤1：确定噪声上限假定其中噪声上限un>0，引入目标函数式：当即使得获得最小特征值λmin所对应的，当趋于无穷时，趋于0，因此对un的约束是，当满足很小时，所对应的un即确定为噪声上限；步骤2：对噪声方差进行估计估计噪声方差如下式：其为对噪声的一致性估计，且其精确度由噪声上限的约束条件控制；步骤3：计算对应的最小特征向量采用奇异值svd分解计算最小特征向量，并代入步骤2中估计噪声方差；步骤4：利用ALS方法求得椭球参数Ce和d，进而求得C；步骤5：利用公式hE＝C-1hM-d，修正三轴传感器输出数据，即成。本专利技术的有益效果是：该方法在确保现场测量可实现的前提下，大大提高了系统的稳健性，扩大了三轴磁传感器基于椭球假设误差模型定标方法的适用范围。附图说明图1是本专利技术带噪声抑制的定标方法与现有方法的原理比较示意图；图2是本专利技术UALS定标方法与现有的OLS、ALS方法的性能仿真效果图；图3是本专利技术UALS定标方法实施例数据采集平台获取的第一组数据样本；图4是本专利技术UALS定标方法实施例数据采集平台获取的第二组数据样本；图5是本专利技术UALS定标方法与OLS基于实验数据的性能对比。具体实施方式下面结合附图和具体实施方式对本专利技术进行详细说明。如图1，现有的定标方法，如OLS(OrdinaryLS)一般不考虑噪声，直接对误差参数进行估计，在噪声较大时得到的误差参数偏差较大，参照图1中的上半部分；本专利技术方法中使用的ALS(AdjustedLeastSquare)可在已知噪声分布特性时得到准确的误差参数；在此基础上，本专利技术进一步改进了ALS算法，得到UALS(NoiseUnknownAdjustedLeastSquare)算法，能够在噪声特性未知的情况下进行三轴磁传感器误差参数的有效估计，参照图1中的下半部分。对于任意给定的测量数据R3表示三维实空间，忽略噪声因素的影响后，在传感器固有误差影响下，对于包含多种误差因素的多姿态观测数据服从椭球分布公式：F(Ce,d)＝{x(l)∈R3:(x(l)-d)TCe(x(l)-d)＝H}，(1)其中H为磁场幅度，Ce和d为椭球参数，其中Ce＝(C-1)TC-1>0，C为乘性误差矩阵，C-1是乘性误差矩阵的逆，d即为加性误差；通过hE＝C-1hM-d将实测磁场hM修正得到准确测量磁场hE，使得修正后的输出磁场幅度hE满足(hE)ThE＝H2。现有三轴磁传感器的定标过程多为求解该最优问题，且多忽略了噪声的影响。由于忽略了噪声因素的影响引入的误差dn表示如下：dn＝||C-1hM-d||2-||hE+C-1d||2＝nTCen+2hEC-1n，(2)n为测量过程中的环境噪声，能够证明dn是严格正定的，且随着n增加，dn增大，此时传统定标方法将失效。常规三轴磁传感器误差参数估计方法仅通过对测量数据进行采样点平均来抑制噪声，其得到的是噪声的非一致性估计。如果通过改进最小二乘法(ALS,AdjustedLeastSquare)得到噪声抑制下的误差参数估计方法，其间对噪声的处理为一致性估计，该改进方法简述如下：对于测量数据l＝1,…,m，其中上标“-”表示真实磁场值，“～”表示测量数据中包含的噪声，则椭球表面表示为：其中于是，最优估计描述为：其中，p＝[vec(Ce)uTv]T,算子vec(*)为将上三角矩阵表示为向量，定义参数：于是，误差参数的改进估计式表示如下：即有此时，假设噪声服从均值为0，方差为的独立同分布，即当l1,l2＝1,…,m且l1≠l2时，若已知，则得到如下的关系式：...

【技术保护点】
一种带噪声抑制的三轴磁传感器定标方法，其特点在于，具体按照以下步骤实施：步骤1：确定噪声上限假定其中噪声上限un>0，引入目标函数式：QU(σn2)=|λmin(YALS(σn2))|,]]>当σn2∈(0,un2),---(7)]]>即使得获得最小特征值λmin所对应的，当趋于无穷时，趋于0，因此对un的约束是，当满足很小时，所对应的un即确定为噪声上限；步骤2：对噪声方差进行估计估计噪声方差如下式：σ^n2=argmin0≤σn2≤un2|λmin(YALS(σn2))|,---(8)]]>其为对噪声的一致性估计，且其精确度由噪声上限的约束条件控制；步骤3：计算对应的最小特征向量采用奇异值svd分解计算最小特征向量，并代入步骤2中估计噪声方差；步骤4：利用ALS方法求得椭球参数Ce和d，进而求得C；步骤5：利用公式hE＝C‑1hM‑d，修正三轴传感器输出数据，即成。

【技术特征摘要】
1.一种带噪声抑制的三轴磁传感器定标方法，其特征在于，具体按照以下步骤实施：步骤1：确定噪声上限假定其中噪声上限un＞0，引入目标函数式：当即使得获得最小特征值λmin所对应的当趋于无穷时，趋于0，因此对un的约束是，当满足时，所对应un上限值即确定为噪声上限；步骤2：对噪声方差进行估计估计噪声方差如下式：其中表示以为自变量代入YALS(x)，其为对噪声的一致性估计，且其精确度由噪声上限的约束条件控制；步骤3：计算对应的最小特征向量采用奇异值svd分解计算最小特征向量，并代入步骤2中估计噪声方差；步骤4：利用ALS方法求得椭球参数Ce和d，根据公式Ce＝(C-1)TC-1求得三轴磁传感器乘性误差矩阵C；所述的ALS方法具体过程是：改进最小二乘法(ALS,AdjustedLeastSquare)得到噪声抑制下的误差参数估计方法，其间对噪声的处理为一致性估计，该改进方法简述如下：对于测量数据其中上标“-”表示真实磁场值，“～”表示测量数据中包含的噪声，则椭球表面表示为：其中H为参考标量磁场，于是，最优估计描述为：

【专利技术属性】
技术研发人员：宋忠国，张金生，朱文奇，席晓莉，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：陕西;61