基于相关向量机的石英加速度计温度误差补偿方法.pdf

第45卷第2 1期2023年11月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.21Nov.,2023基于相关向量机的石英加速度计温度误差补偿方法宁云晖1，赵凤昊，王凯，李硕，尚岳(1.海军研究院，北京10 0 16 1；2.天津航海仪器研究所，天津30 0 131)摘要：加速度计测量误差决定了惯导系统中位置、姿态、速度等导航信息精度。针对环境温度变化对加速度计测量精度的影响展开讨论，分析了工作温度与测量误差间的关系；提出一种基于相关向量机(Relevancevectormachine，RVM)的加速度计温度误差建模及补偿方法，能够实现温补模型参数的准确与快速解算。结果表明，补偿后的加速度计在10 50 的温度范围内测量误差减小了50%以上，提高了其对环境温度变化的适应能力。关键词：加速度计；温度误差补偿；相关向量机中图分类号：U666.12文章编号：16 7 2-7 6 49(2 0 2 3)2 1-0 17 1-0 6RVM-based temperature error compensation method for quartz accelerometer(1.Naval Academy,Beijing 100161,China;2.Tianjin Institute of Nautical Instruments,Tianjin 300131,China)Abstract:The accuracy of the position,attitude,and speed in navigation system is largely determined by the accuracyof accelerometer.Discussion was made in this paper,which is about the relationship between the errors and the temperature.Furthermore,a relevance vector machine based modeling and compensation method was proposed,which could compensatethe temperature errors accuracy and rapidly.The result of verification experiment shows that,the measurement error of theaccelerometer reduced by 50%at a temperature period of 10 50,w h i c h c o u l d i m p r o v e t h e a d a p t a b i l i t y t o t h e t e m p e r a t-ure variations.Key words:quartz accelerometer components;temperature error compensation;relevance vector machine0引言加速度计作为惯性导航系统中的核心惯性器件之一，主要用于观测运载体在测量参考坐标系下的线加速度，进而推算运载体在导航解算坐标系下的位置、速度以及姿态信息，其测量精度直接决定了系统的导航精度水平1-3。目前，在中高精度导航系统中应用最为广泛的是石英挠性加速度计，其具有精度高、可靠性强、长期稳定性好等优点，相较于其他类型的加速计，在导航领域更受青睐。理想情况下，加速度计的输出应与其敏感的加速度呈线性关系。但受到制造工艺、内部结构、采集线路等因素影响，不同元件间的输出模型存在一定差异，一般需要在实验室条件下对包括标度因数、偏值等参数进行标校后方能使用。然而，在实际工作环境收稿日期：2 0 2 2 11-0 3作者简介：宁云晖（19 7 3-），女，硕士，副研究员，研究方向为潜艇作战系统及人机交互技术。文献标识码：ANING Yun-huil,ZHAO Feng-hao,WANG Kai?,LI Shuo,SHANG Yue?doi:10.3404/j.issn.1672-7649.2023.21.032下，温度的波动相较于实验室环境波动剧烈，其可能会引起加速度计内部结构发生改变，导致已标校的参数失效，最终导致测量精度下降，因此必须采用必要措施予以控制补偿4-6 。王淑娟等7 构建了石英挠性摆式加速度计的温度模型，并通过单轴十二位置翻滚实验及正交多项式拟合方法对模型参数进行了辨识与补偿后，在2 0 50 范围内加速度计的测量精度提高了一个数量级。翁海娜等8 利用最小二乘法建立了五阶加速度计的偏值和标度因数误差模型，并在30 7 0 范围内对模型参数进行了拟合补偿，补偿后水平状态下加速度计输出精度提高了一个数量级，垂直状态下也提升了6 0%以上。2 0 18 年，俞茂超等9 1分析了石英加速度计的基本原理与温度误差特性，建立了多参量联合多项式模型对加速度计进行补偿，补偿后的加速度计偏值误差降低了7 5%。张金云等10 提出172采用多项式样条函数回归方法对惯性平台系统内石英加速度计进行建模与补偿，补偿后加速计在2 0 50范围内的标度因数稳定性有2 10-4提升到8 10,偏值由1.5LSB/s提升至0.3LSB/s。总的来看，尽管基于多项式的温漂模型能够对加速度计参数与温度变化间的关系进行有效描述与补偿，但难以描述包括温漂滞环效应在内的多种非线性效应，导致补偿不完全，影响测量精度。近年来，随着机器学习研究的不断深人，支持向量机、神经网络、粒子群算法等先进算法被用于加速度计的温漂模型辨识与补偿，有效提升了现有加速度计在变温环境下的测量精度1-13。作为以一种衍生于支持向量机的稀疏概率模型，相关向量机(RelevanceVectorMachine)被Tipping提出，其基于贝叶斯统计学习理论，通过移除不相关的点将模型稀疏化，从而极大减少了核函数所需的计算量，同时也克服了核函数必须满足Mercer条件的限制，在非线性模型的拟合应用中呈现出极大的优越性14。本文在对加速度计温度模型进行分析的基础上，采用相关向量机对模型参数进行估计，并设计了实验对所述方法的有效性进行验证。实验结果表明，该方法能够有效提升变温条件下的加速度计输出精度，补偿后的温度误差相较于补偿前降低了50%以上。1石英加速度计温度误差特性分析典型的石英加速度计结构如图1所示，当加速度计的敏感方向存在加速度时，由石英摆片构成的检测质量将偏离平衡位置，进而引起平衡电容传感器输出发生改变，该信号经放大滤波后，一部分反馈至力矩器线圈，产生使检测质量回到平衡位置的电磁力，一部分经信号采集电路转换后作为加速度计的信号输出。因此，加速度计的输出精度受摆组件的摆性、力矩器电磁力矩系数以及结构材料的稳定性决定。根据IEEE1293-2018线性单轴非陀螺加速度计标准规范格式及推荐测试流程，加速度计的输出公式可表示为检测质量挠性支撑永磁体图1石英加速度计结构示意图Fig.1The structure of quartz accelerometer舰船科学技术KoE=Ki3Ko+2+8。n=2其中，E为加速度计的输出；a为加速度计敏感轴方向输入的加速度分量；Ko为偏值，Kn(n=1,2,3.)为n阶标度因数；Kn为对应的标度因数不对称性；signO为取符号运算函数；8 为包括轴间耦合效应、门摆态变化等因素引起的误差。在中低精度的惯性导航系统中，除一阶标度因数与偏值外的其他因素对加速度计输出的影响极小，因此式(1)一般简化为式(2)的形式使用：E=Ki(Ko+a)。(2)所以加速度值可根据式(3)计算：Ea=K1-Ko实际上，受加工装调精度以及线路特性影响，加速度计间的输出特性均会存在一定差异，为保证加速度计能够输出准确结果，需要预先对其标度因数与偏值进行精确标定，以构建式(2)的输出模型。然而在真实的应用条件下，加速度计的工作温度将随外界环境发生剧烈变化，特别是对于没有温控措施的捷联式系统，加速度计内部的元件特性，如永磁体的磁感应强度、挠性支撑的结构等均会发生改变，最终反映到加速度计的参数上。此时，在实验室条件下标定得到的参数将失效，在最终的测量结果中引人误差。图2 为某石英加速度计在变温条件下的加速度输出曲线，整个试验过程加速度计始终处于静止状态，因此加速度计敏感的加速度在测量过程中未发生变化。随着环境温度逐步升高，加速度计的输出随之改变，在整个温变过程中加速度计输出变化了1.2 mg，显然这将对导航精度造成严重影响。因此，对加速度计的一加速度计输出-温度0.91100.91080.9106磁路0.91040.9102力矩器线圈0.91000.90980.909605.000图2 变温条件下加速度计输出曲线Fig.2The output of the accelerometer at a varying temperature第45卷KTS19n(a)Ja+250454035302520151050005.0002025时间/s(1)(3)000(8)202第45卷温度误差特性进行分析研究，并进行合理补偿，对提升惯导系统的测量精度具有重要意义。考虑加速度计的输出仅与标度因数与偏值相关，在变温条件下，式(2)应改写为：E(T)=Ki(T)a+Ki(T)Ko(T)。将Ki(T)与Ko(T)关于To展开，表现为多项式形式：Ki(T)=Ki(To)+)(5)(To)(T)。0Ko(T)=Ko(To)+)之n！n=1其中，To为标定参数的参考温度；K(To)和K(To)为标度因数与偏值的n阶导数。加速度计温度误差a 与温度变化量T存在一定近似线性关系，如图3所示，一般仅考虑一、二阶导数项即可加大降低温度误差影响。然而在实际应用过程中，该补偿方法仍存在一定缺陷。首先，在工作状态下，惯性导航系统内部复杂的温度场和外界环境温度的剧烈变化加上自身的内部器件发热，导致直接获取加速度计核心位置温度存在一定困难；其次，多项式温补模型存在高阶舍入误差，非线性温度效应的影响无法完全消除。最后，实际应用环境下，加速度计温度是自然连续变化的，传统的温度梯度标定方法忽略了温度变化率对结果的影响。因此，为进一步提升加速度计测量精度，建立温度采样数据与温度误差数据之间的非线性模型关系，而模型及参数的获取则成为解决这一问题的关键。0.91100.91080.91060.91040.91020.91000.90980.90965101520253035404550温度/图3加速度计输出与温度变化间的关系Fig.3 The relationship between output of the accelerometerand the temperature variation2基于RVM的加速度计组件温度误差建模方法传统的非线性模型建模方法包括多元回归，支持向量机、神经网络等，各有利弊。RVM是一种衍生于支持向量机的回归和分类算法，其基于概率学习思想，通过在权值上定义超参数来控制Gaussian先宁云晖，等：基于相关向量机的石英加速度计温度误差补偿方法2.1相关向量机函数拟合方法(4)设训练样本输入目标集为(3xn,-1，即tn=y(xn;w)+8n式中，8 n为附加噪声。且满足高斯分布snN(0,2)，则K(n(To)(AT),n！n=1一加速度计输出:173验概率，在贝叶斯框架下进行机器学习，可以计算输出的概率分布、很自然的适用于复杂类温度误差模型，并利用判据来移除野值，从而获得稀疏化模型。(6)RVM的输出模型为：(7)i=1式中，K(x,xi)为核函数，;为权值，且d(xi)=K(x,xi)。假设数据集tn服从独立分布，则其条件概率为：NP(tl,o2)=L I Ntil(xi;w),2=0(2元0 2)-N/2 exp/式中，向量t=(t,.tn)T，w=(w 1,.n)T，为N(N+1)维矩阵,(xn)=1,K(xn,x1),K(xn,xN)T。因此，对于给定样本集x*，预测输出集t的条件概率为：P(t t)=f P(t/,2)P(w,o2t)dwdo2。在贝叶斯框架下，该似然概率可设计到无限大，但为防止过拟合，为每个权值定义高斯先验概率分布来约束参数：P(ula)=L I N(w;l0,ai)。=0式中，=(0,)T为由超参数组成的(N+1)维向量。因此，式(9)可转化为：P(tl)=Jf P(t lu,2)P(w,o-1)dwdado2。(11)依据贝叶斯准则，通过已定义的先验概率，得到后验概率：P(tl,a,o2)P(w,a,o2)P(u,0,02l)=PP(t)因无法直接计算后验概率P(tlw,2)，将其分解为：P(u,02lt)=P(lt,2)P(,2t)。(13)将式(13)代人式(11)有:(9)N(10)(12)174 P(tal)=JJf P(tla,a,02)P(ul,02)P(a,21)dwdado2。(14)则权重的后验概率分布为：-1/2P(tlo,o2)exp(-1(a-2(a-),后验协方差和均值分别为：Z=(-2 sT$+A)1,=9-2式中,A=diag(o,n)。由于后验概率P(,2It)不能直接分解计算得到，采用狄拉克(Dirac Delta)函数近似计算：P(a,0-1)=0(a MP,0Mp)其中，MP,0Mp为后验概率P(,o21)的最优解，即Mp=arg max(P(alt),0/iup=arg max(P(a-10),通过上述转化有：P(,2t)P(la,2)P(a)P(2)。因此有：P(tla,o2)=P(lw,o)P(ala)dw=(2m)-/-2+0-1l0T-1/2exp-(o21+0A-l0T)-1对式(2 1)右边进行极大似然估计，可得：NL=log(P(llog a,log。)+Zl o g P(a i)+l o g P(.2)。1=0将式(17)代人式(2 3)整理得：+1og/A-02/(t-dw)N(alogai)+Nlogo?2=0对式(2 4)求偏导得：舰船科学技术8Llogi式中：为式17)中矩阵所对应的第i行平均权重值；Z;为式(16)后验协方差矩阵Z的第i个对角元素。假设ri=1-iZi，令式(2 5)为0，则有：(15)同理得：(16)(17)(18)(19)(20)(21)(22)(23)(24)第45卷1(25)2Tit-dull02NN-2ri=0采用狄拉克(DiracDelta)函数可将式(14)近似转化为：P(tal)=J P(tal,a MP,0Mp)P(ul,a MP,0Mp)do。(2 8)又由于 P(al,a MP,0MpP)和P(al,MP,0Mp)均满足正态分布，因此：P(t l)=N(tl,02),Ju=T0(x),0,=0Mp+0(x.)TE0(x),0(x)=1,(,1),.,K(xx,xN)T。选择合适的核函数对样本数据进行训练，将特征向量映射到高维空间。给定，2 初值，通过式(26)和式(2 7)进行迭代更新，待所有参数均满足收敛条件，迭代结束。此时，RVM回归预测模型确认完成。2.2加速度计温度误差建模与补偿在RVM拟合方法中，核函数的选择对回归预测的计算复杂度和精确度产生直接影响。为兼顾精度与拟合速度，采用高斯核函数和多项式核函数构成的混合核函数，对加速度计的温度误差模型进行拟合：K(x,xi)=Aexp(-x-x:/203)+(1-M)(x xi+1)d。式中：2为高斯核函数的宽度；d为多项式核函数参数；入为权重参数；x为核函数中心。基于RVM函数拟合的加速度计温度误差补偿方法主要包括数据归一、参数寻优和模型训练3个步骤进行。1)数据归一为优化内存空间，加快参数寻优以及模型训练的(26)(27)(29)(30)(31)(32)(33)第45卷速度，首先需要对加速度计组件温度数据和组件输出数据进行归一化处理，并记录归一化模型，处理方法如下式：ymax-ymin(x-Xmin)+yminoy=Xmax-Xmin式中，y为归一化输出数据，ymax=1，y mi n=-1，x 为归一化输人数据，Xmax为输人数据的最大值，Xmin为最小值。2)参数寻优选择不同的核参数与权值参数将影响混合核函数的性能。采用网格寻优来优化混合核参数d和2以及权重参数几的方法，来保证混合核函数具有更好的学习能力。根据经验，d，人和12 的取值范围为2-10 2 0 ,取0.1为步长，将d，和2 的取值以横、纵排列方式构成类似网格阵列。针对每个网格单元参数可得到对应的拟合函数模型，选取拟合准确率最高的一组参数作为最优参数，对应最优参数为dbest，A b e s 和obest。3)模型训练将归一化数据及最优参数代人进行训练，得到最佳拟合函数。此时根据实测数据温度变化即可推算对应的标度因数与偏值改变量，从而消除温度误差对测量结果的影响。3实验为验证所述方法，设计了实验对其有效性进行评估。具体实验装置如图4所示，待测加速度计通过安装孔固连在二维转台中央温箱内的台面上。调整工装，使加速度计的敏感轴方向在二维转台的外环与中环框架角度输出为0 时，垂直指地。温箱的温度可通过外部控制电流在-30 +55的范围内随意调整，温度传感器测量温箱内的实时温度，并与加速度计输出的脉冲、二维转台的实时角度通过数据采集系统同步采集，发送到上位机。加速度计二维转台图4实验装置示意图Fig.4Schematic diagram of the experiment宁云晖，等：基于相关向量机的石英加速度计温度误差补偿方法温箱175根据式(4)，加速度计的偏值与标度因数均存在温度效应，二者的温度特性并不相同，因此在实验过程中，首先应考虑如何分离两项参数以避免耦合效应影响温度误差建模。当加速度计敏感轴方向位于水平面(34)内时，其敏感的加速度应该为0，此时加速度计的输出仅受到偏值的影响；而加速度计敏感轴垂直本地水平面时，其敏感的加速度达到最大，为本地重力加速度G，从输出结果中扣除偏值温度误差，即可对标度因数的误差模型进行建模，从而实现二者的分离。具体的实验过程如下：首先控制二维转台x轴正向旋转90，使加速度计敏感轴水平；控制温箱进行温度循环，温度循环范围设定为10 50，循环周期1h，以10 0 Hz的频率同步记录温度传感器及加速度计的输出数据集为(To,JO)-1，采用上位机软件通过相关向量机拟合加速度计偏值温度误差模型，并对加速度计的输出结果进行补偿。控制二维转台的x轴复位，此时加速度计敏感轴垂直指天，同样控制温箱进行相同的温度循环，并以10 0 Hz的频率同步记录温度及加速度计输出，并根据前面建立的加速度计偏值模型，扣除加速度计输出中偏值波动的影响，得到的数据集记作(T%,0)%1，采用相关向量机建立加速度计标度因数温度误差模型，并对输出结果进行补偿。补偿前后的加速度计输出如图5所示，补偿前水平状态下加速度计的输出随温度变化发生明显改变，最大波动幅度达到了0.49 mG，经补偿后该误差被有效抑制，最大波动幅度降低了约50%，为0.2 7 mG。垂直状态下的数据同样可得到相同结论，补偿前的输出随温度变化的最大幅值为1.4mG，补偿后该变化降低到0.2 8 mG，温度误差降低了一个数量级。4结语加速度计温度误差是惯性导航系统中不可忽视的误差之一，对导航精度有着重要影响，而自然工况下环境温度的变化可能引起加速度计标度因数以及零偏发生改变，造成测量精度严重下降，因此必须采取措施予以补偿。本文提出了一种基于相关向量机加速度计温度误差建模与补偿方法，并通过实验验证了该方法的有效性。经实验验证，补偿后的加速度计输出精度相较于补偿前提升了50%以上，温度误差得到了有效抑制。后续将进一步开展加速度计温度误差模型的分析研究，主要完成器件结构温度变形机理等工作，为加速度计的温度误差补偿提供理论支撑。176Fig.5 Result of the temperature errors compensation experiment参考文献：1 王巍.惯性技术研究现状及发展趋势.自动化学报，2 0 13,39(6):723729.WANG Wei.Status and development trend of inertial technol-ogyJ.Acta automatica sinica,2013,39(6):723-729.2严恭敏.捷联惯导算法及车载组合导航系统研究D.西安：西北工业大学，2 0 0 4.【3 崔立尉.惯性捷联导航系统的标定和数据对准技术1.舰船舰船科学技术0.2科学技术,2 0 2 2,44(12):16 5-16 8.0.1CUI Liwei.Research on calibration and data alignment tech-0nology of inertial strap-down navigation systemJJ.Ship sci-ence and technology,2022,44(12):165-168.0.1【4】郭祎，李威，田兴，等.温度测试中石英加速度计热传递方式0.2研究.新技术新工艺,2 0 2 2(7):：55-6 0.-0.3GUO Yi,LI Wei,TIAN Xing,et al.Research on heat transfer0.4mode of quartz accelerometer in temperature measurementJ.1020温度/(a)加速度为0 时未补偿的加速度计输出0.30.20.10-0.10.20.310(b)加速度为0 时补偿后的加速度计输出0.80.4Du/00.40.810(c)加速度为G时未补偿的加速度计输出0.40.200.2-0.410(d)加速度为G时补偿后的加速度计输出图5加速度计温度误差补偿实验结果foraccelerometer第45卷30402030温度/2030温度/2030温度/50405040504050New Technology&New Process,2022(7):55-60.【5】金鑫.集成式石英谐振加速度计及温度补偿研究D.北京：北京信息科技大学，2 0 2 2.【6】陈雪.一种基于惯性器件的温度补偿模型研究D.成都：电子科技大学，2 0 2 2.【7】王淑娟，黄显林，刘升才.加速度计温度模型的辨识.中国惯性技术学报,19 9 7(1):32-37.WANG Shujuan,HUANG Xianlin,LIU Shengcai.Identifica-tion for temperature model of accelerometerJ.Journal of chi-nese inertial technolgy,1997(1):32-37.【8 翁海娜，胡小毛，裴志，等.一种新的加速度计温度误差补偿方法.中国惯性技术学报，2 0 0 9,17(4):47 9-48 2.WENG Haina,HU Xiaomao,PEI Zhi,et al.Novel method oftemperature error compensation for accelerometerJJ.Journalof chinese inertial technolgy,2009,17(4):479-482.【9 俞茂超，王新宇，阳洪，等.石英挠性加速度计温度补偿方法研究.压电与声光,2 0 18,40(2):2 2 0-2 2 2.YU Maochao,WANG Xinyu,YANG Hong,et al.Study ontemperature compensation method of quartz flexible ac-celerometerJ.Piezoelectrics and acoustooptics,2018,40(2):220-222.10】张金云，姜欢，赵军虎，等.惯性平台系统石英加速度计温度建模补偿技术J.导航与控制,2 0 2 1,2 0(3):6 6-7 3.ZHANG Jinyun,JIANG Huan,ZHAO Junhu,et al.Researchon temperature modeling and compensation technique ofquartz accelerometer in inertial platform systemJ.Navigationand control,2021,20(3):66-73.11】于湘涛,张兰,郭琳瑞,等.基于小波最小二乘支持向量机的加速度计温度建模和补偿.中国惯性技术学报，2 0 11,19(1):95-98.12 惠进，郭栓运，尹剑.石英振梁加速度计的小波神经网络温度补偿研究J.兵器装备工程学报,2 0 16,37(7)：118-12 2.HUI Jin,GUO Shuanyun,YIN Jian.Temperature compensa-tion of quartz vibrating beam accelerometer based on waveletneural networkJ.Journal of Ordnance Equipment Engineer-ing,2016,37(7):118-122.13于东康，杨功流,吴宜荣,等.基于粒子群算法的石英挠性加速度计温度补偿方法研究.导航定位与授时,2 0 18,5(5):90-95.14】杨树仁,沈洪远.基于相关向量机的机器学习算法研究与应用.计算技术与自动化,2 0 10,2 9(1):43-47.YANG Shuren,SHEN Hongyuan.Research and application ofmachine learning algorithm based on relevance vector ma-:chine.Computing technology and automation,2010,29(1):43-47.