一种双轴旋转惯性导航系统的重要参数快速自标定方法.pdf

第1 0卷 第3期2 0 2 3年5月导航定位与授时N a v i g a t i o nP o s i t i o n i n ga n dT i m i n gV o l.1 0N o.3M a y2 0 2 3d o i:1 0.1 9 3 0 6/j.c n k i.2 0 9 5-8 1 1 0.2 0 2 3.0 3.0 0 9一种双轴旋转惯性导航系统的重要参数快速自标定方法尹俊杰,李志强,陈 雨,王 盛,宋天骁(北京航天发射技术研究所,北京1 0 0 0 7 6)摘 要:为了提高双轴旋转惯导重要参数标定的快速性和精度,提出一种快速自标定方法。通过设置不同的标定路径可以在1 0m i n内完成陀螺和加速度计的零偏以及标度因数误差的标定。该方法利用基于姿态误差观测的卡尔曼滤波完成陀螺零偏的估计。通过六位置翻滚并以速度误差作为观测量进行卡尔曼滤波,完成加速度计的零偏及标度因数误差的标定。使天向陀螺绕方位轴旋转4周,使水平陀螺绕水平轴转动4周,通过计算旋转前后的姿态误差完成陀螺标度因数误差的估计。仿真和试验结果表明,该方法可以实现双轴旋转惯导重要参数1 0m i n内完成自标定,且具有较高的精度。关键词:快速自标定;卡尔曼滤波;双轴旋转式惯导中图分类号:V 2 4 9.3 2+2 文献标志码:A 文章编号:2 0 9 5-8 1 1 0(2 0 2 3)0 3-0 0 8 0-0 9F a s t s e l f-c a l i b r a t i o no fk e yp a r a m e t e r sb a s e do nd u a l-a x i sr o t a t i o n i n e r t i a ln a v i g a t i o ns y s t e mY I NJ u n j i e,L IZ h i q i a n g,C h e nY u,WAN GS h e n g,S ON GT i a n x i a o(B e i j i n gI n s t i t u t eo fS p a c eL a u n c hT e c h n o l o g y,B e i j i n g1 0 0 0 7 6,C h i n a)A b s t r a c t:T or e d u c e t h ec a l i b r a t i o nt i m ea n d i m p r o v ec a l i b r a t i o na c c u r a c yo fd u a l-a x i sr o t a t i o ni n-e r t i a ln a v i g a t i o ns y s t e m(D R I N S),af a s ts e l f-c a l i b r a t i o nm e t h o do fk e yp a r a m e t e r s i sp r o p o s e d.T h eb i a s a n ds c a l e f a c t o r e r r o r o f g y r oa n da c c e l e r o m e t e r c a nb e c a l i b r a t e dw i t h i n1 0m i nb ys e t t i n gd i f f e r e n t c a l i b r a t i o np a t h s.T h ee s t i m a t i o no fg y r ob i a sc a nb eo b t a i n e db yu s i n gK a l m a nf i l t e rb a s e do no b s e r v a t i o no f a t t i t u d e e r r o r.B ya s i x-p o s i t i o nc a l i b r a t i o np a t h,t h e c a l i b r a t o no f b i a s a n ds c a l e f a c t o re r r o ro f a c c e l e r o m e t e rc a nb ea c h i e v e dw i t hK a l m a nf i l t e rb yo b s e r v i n gv e l o c i t ye r r o r.T om a k e t h eu p w a r dg y r or e v o l v e4c i r c l e sa r o u n da z i m u t ha x i sa n dt h eo t h e r t w og y r o s r e v o l v e4c i r c l e sa r o u n dp i t c ha x i s r e s p e c t i v e l yc a no b t a i nt h e s c a l e f a c t o r e r r o ro f g y r ob ya n a l y z i n g t h ea t t i-t u d ee r r o r sb e f o r ea n da f t e r r e v o l v i n g.T h e r e s u l t so f s i m u l a t i o na n de x p e r i m e n t s i n d i c a t e t h a t t h ef a s t s e l f-c a l i b r a t i o no fk e yp a r a m e t e r so f t h eD R I N Sc a nb ed o n ew i t h i n1 0m i nu s i n g t h ep r o p o s e dm e t h o da n d i sa l s oo fh i g ha c c u r a c y.K e yw o r d s:F a s ts e l f-c a l i b r a t i o n;K a l m a nf i l t e r;D u a l-a x i sr o t a t i o ni n e r t i a ln a v i g a t i o n s y s t e m(D R I N S)收稿日期:2 0 2 2-1 2-2 9;修订日期:2 0 2 3-0 3-0 7作者简介:尹俊杰(1 9 8 7-),男,硕士,工程师,主要从事惯性导航与组合导航技术方面的研究。第3期一种双轴旋转惯性导航系统的重要参数快速自标定方法0 引言惯性器件误差标定技术是惯导系统的关键技术之一,通过事先对系统误差参数进行标定,在实际导航过程中对误差进行补偿,从而消除系统误差参数对导航精度的影响1-3。双轴旋转惯导系统在出厂前会进行全参数(陀螺和加速度计的零偏、标度因数、安装偏差角、加速度计二次项系数、杆臂误差等)的标定4-5,再将其用于系统上完成导航。当双轴旋转惯导随系统交付后,系统惯性器件误差会随着服役时间的加长而增大,为了消除这种长期稳定性误差,用户一般会提出定期标定的保障需求6。文献6 提出了一种双轴旋转惯导系统误差自标定技术,结合设计的自标定路径和卡尔曼滤波算法,在1h内实现了绝大部分惯性器件误差的分离和标定。文献7 采用多位置标定方法可以在1.5h内完成激光陀螺捷联惯导系统的全参数自标定。目前国内研究机构的研究方向主要关注全参数系统级标定,标定时间一般在1h以上,而针对惯性器件重要参数快速自标定方法还鲜有研究。陀螺和加速度计的零偏、标度因数等重要参数随着时间推移变化较快,其他参数在惯导系统不拆卸的前提下变化比较缓慢甚至基本保持不变。为了进一步提高标定的快速性,本文提出了一种惯性器件的标度因数和零偏等重要参数的快速自标定方法。该方法通过设计不同的状态变量,通过观测不同的系统导航误差,利用卡尔曼滤波算法依次估计出相应的惯性器件重要参数。仿真分析和实物试验结果表明,本文提出的快速自标定方法在1 0m i n内可以实现重要参数的自标定,且具有易于实现、精度高的优点,是一种可以工程应用的重要参数快速自标定方法。1 基于双轴惯导的快速自标定原理坐标系选取如下:选取东-北-天坐标系作为导航坐标系(n系);选取惯组右-前-上坐标系作为载体坐标系(b系);记地心惯性坐标系为i系;记地球坐标系为e系。1.1 标定误差模型快速自标定在惯组已经完成出厂标定的基础上进行。通过将激光惯组转至特定位置后解析待标定参数误差方程实现。陀螺和加速度计的快速标定误差模型可写成8-1 1bi b=KGbi b+b(1)fb=KAfb+b(2)式中KG=KgxKgyKgzbi b=bi b xbi b ybi b zTb=bxbybzTKA=KaxKayKazfb=fbxfbyfbzT,b=bxbybzT其中,Kgi、Kai(i=x,y,z)分别表示陀螺和加速度计的标度因数误差;bi b j、fbj(j=x,y,z)分别表示陀螺和加速度计的输入;bi、bi(i=x,y,z)分别表示陀螺和加速度计的常值零偏。需要进行标定的重要参数为陀螺标度因数误差Kgi和加速度计的标度因数误差Kai,陀螺零偏bi以及加速度计零偏bi(i=x,y,z)。1.2 快速自标定原理分3步顺序实现快速自标定:第一步,对当前姿态与对准完成时刻姿态的姿态误差进行观测并设计卡尔曼滤波器实现陀螺零偏进行估计;第二步,按设定的标定路径通过对速度误差进行观测并设计卡尔曼滤波器实现加速度计的零偏及其标定因数误差的标定;最后采用数值计算的方式通过双轴惯导系统的两个旋转轴对陀螺的标度因数误差进行计算。1.2.1 陀螺零偏标定原理由于双轴惯导系统交付用户后已经完成出厂标定,因此静态条件下,由地球自转角速度和标度因数误差产生的角速度测量误差较小,估计陀螺零偏时,该误差可忽略。两位置自对准完成后,可以估计出b系与n系的姿态矩阵Cnb,然后以3个失准角信息和三轴陀螺零偏作为状态量,以滤波过程中的姿态信息与自对准完成时刻“凝固”的姿态信息的姿态角误差作为观测量,通过构造卡尔曼滤波器估计三轴陀螺的零偏。陀螺零偏估计原理如图1所示。状态方程和观测方程如式(3)所示1 2Xg=FgXg+GgWbgZg=HgXg+Vg(3)式中18 导航定位与授时2 0 2 3年5月Fg=Ma a-CnbO3 33 3,Gg=-CnbO3 3Hg=I3O3 3Ma a=-(ni n)ni n=ni e+ne n其中,Wbg表示陀螺仪角速度测量白噪声;Vg表示量测噪声。状态变量Xg取为Xg=ENUbxbybzT(4)观测量Zg取为Zg=ENUT(5)式中,姿态角误差的获取方法如下。将对准完成时刻的姿态矩阵Cnb“凝固”,记其转换后的四元数为qn b0,记滤波过程中实时更新的姿态四元数为qn b,则误差四元数qn be r r可表示为1 2qn be r r=qe 0qe 1qe 2qe 3T=qn b(qn b0)*(6)式中,“”表示四元数乘法;(qn b0)*为qn b0的共轭四元数。姿态角误差为=2ac o s(qe 0)/s i n(ac o s(qe 0)qe 1qe 2qe 3T(7)图1 陀螺零偏估计原理F i g.1 P r i n c i p l eo f g y r ob i a s e s t i m a t i o n1.2.2 加速度计标定原理本方案采用卡尔曼滤波器对加速度计的标度因数误差以及加速度计零偏进行估计,状态方程和观测方程如式(8)所示1 2-1 4Xa=FaXa+GaWbaZa=Vn=HaXa+Va(8)式中Fa=Ma aMavO3 3O3 3MvaMv vMv kaCnbO6 3O6 3O6 3O6 3Mav=0-1/(RM+h)01/(RN+h)00t a nL/(RN+h)00Ga=O3 3CnbO6 3,Mva=fnfn=CnbfbxfbyfbzTMv v=(Vn)Mav-(ni e)Mv ka=Cnb(1,1)fbxCnb(1,2)fbyCnb(1,3)fbzCnb(2,1)fbxCnb(2,2)fbyCnb(2,3)fbzCnb(3,1)fbxCnb(3,2)fbyCnb(3,3)fbz其中,L为当地纬度;h为当地高度;RM为子午圈主曲率半径;RN为卯酉圈主曲率半径;fbi(i=x,y,z)为载体系下的比力分量;状态量Xa取为Xa=E N U VE VN VU KaxKay Kaz x y zT(9)观测量Za取为Za=VE VN VUT(1 0)为了提高对加速度计标度因数误差和零偏的可观测性,实际进行滤波时,在完成陀螺零偏估计后,以对准结果作为初始状态,并将陀螺零偏补偿后,设定转位机构标定路径,使每个轴向的加速度计均能朝上朝下并静止一段时间,从而完成加速度计的参数标定。加速度计零偏和标度因数估计原理如图2所示。图2 加速度计零偏及标度因数误差估计原理F i g.2 P r i n c i p l eo fa c c e l e r o m e t e rb i a sa n ds c a l ef a c t o r e r r o r e s t i m a t i o n1.2.3 陀螺标定因数标定陀螺标度因数标定采用数值解析的方式进行解算。完成陀螺零偏和加速度计标度因数误差以及加速度计零偏标定后,进行陀螺标度因数误差的标定。将陀螺各轴恢复到旋转框架零位,然后绕方位轴电机旋转整数周后,计算当前方位角与旋转前方位角的差值,由于零偏已经补偿,因此该值即为方位轴标度因数带来的偏差,由此可算出方位轴陀螺的标度因数误差。同理,将另外两个陀螺均转到28 第3期一种双轴旋转惯性导航系统的重要参数快速自标定方法俯仰轴,分别绕俯仰轴旋转整周,根据旋转后的俯仰角与旋转前俯仰角的差值可依次标定出另外两个轴陀螺的标度因数误差。2 快速自标定流程根据前面的描述,标定顺序安排如下:系统上电完成自检后,收到快速自标定指令,首先进行两位置自对准,自对准完成后,可以得到姿态的四元数,通过四元数可以转换为姿态矩阵Cnb,然后以实时姿态角与自对准完成后“凝固”的姿态差值作为观测量进行滤波,从而完成陀螺零偏的标定。然后通过转位机构分别使三轴加速度计朝上朝下各采集一定时间后,利用卡尔曼滤波估计完成对加速度计的标度因数误差和零位的标定。最后先将天向陀螺旋转4周,利用导航得到的方位与旋转前的方位差可算得天向陀螺标度因数,再将水平面内的陀螺分别绕俯仰轴旋转4周,利用导航得到的俯仰角与旋转前的俯仰角差值可分别算得水平两个陀螺的标度因数误差。两位置自对准总时间为3 0 5s,对准完成后获得姿态阵Cnb,从第二个对准位置开始存储惯性器件数据,然后开始下述的标定流程。2.1 陀螺零偏标定为了缩短标定时间和充分利用历史数据,在实际操作过程中,可以将对准过程中第二个对准位置的惯性器件信息进行存储和复用。即,完成对准后,逆向导航至第二个对准位置开始时刻,并开始以纯惯导导航的姿态信息输出与“凝固”的姿态信息之差作为观测量进行卡尔曼滤波,同时继续采集惯性器件数据5 0s。此过程共耗时5 0s。2.2 加速度计零偏及标度因数标定1)完成陀螺零偏标定后,此时z轴加速度计指向天,静态采集z轴加速度计数据1 5s,记录数据。转动俯仰轴电机9 0,使y轴加速度计指天,转动时间为t1,采集y轴加速度计数据1 5s,记录数据。2)继续转动俯仰轴电机9 0,使z轴加速度计指向地,转动时间t1,采集z轴加速度计数据1 5s,记录数据。3)继续转动俯仰轴电机9 0,使y轴加速度计指向地,转动时间t1,采集y轴加速度计数据1 5s,记录数据。4)转动方位轴电机9 0,使x轴加速度计指向前,转动时间t1,再转动俯仰轴电机9 0 使x轴指向天,转动时间t1,静态采集z轴加速度计数据1 5s,记录数据。5)转动俯仰轴电机1 8 0,使x轴加速度计指向地,转动时间t2,静态采集x轴 加速度计数据1 5s,记录数据。表1 加速度计标定过程电机转动安排T a b l e1 A r r a n g e m e n t o fm o t o r r o t a t i o n i na c c e l e r o m e t e r c a l i b r a t i o n序号电机转动角度/()俯仰轴电机方位轴电机停留时间/s备注1001 52+9 001 53+9 001 54+9 001 550+9 006+9 001 57+1 8 001 5完成后转入下阶段将采集的数据,按照前述的方法,利用卡尔曼滤波估计可完成对加速度计零偏和标度因数误差的标定。在实际标定过程中,设计使t1=t2=5s。根据上述内容可知,加速度计标定总耗时为1 2 0s。2.3 陀螺标度因数误差标定将三轴陀螺恢复至初始位置并进行调平,1 0s内完成,记录此时的方位角,然后先将z轴陀螺旋转14 4 0,转动时间为3 6s,将此时的方位角信息与转动前的方位角输出作差,记转换后的方位角差值为1(单位为),则z轴标度因数为(1+Kgz)T0dt=14 4 0+1(1 1)式中,T0dt=14 4 0,代入式(1 1)可得z轴陀螺标度因数为1+Kgz=1+114 4 0(1 2)则Kgz=114 4 0(1 3)然后将水平面内的陀螺分别绕俯仰轴旋转一定圈数,记导航俯仰角信息与转动前的俯仰角输出之差分别为1、2(单位为),按照上述的计算方法可分别得Kgx=114 4 0(1 4)38 导航定位与授时2 0 2 3年5月Kgy=214 4 0(1 5)在将z轴转换到前向过程中,需要5s的调整时间,由上述可知,陀螺标度因数标定总耗时为(3 63+1 0+5)s=1 2 3s;综上,快速自标定共耗时(3 0 5+5 0+1 2 0+1 2 3)s=5 9 8s,标定流程如图3所示。图3 快速自标定流程F i g.3 F l o wo f f a s t s e l f-c a l i b r a t i o n3 仿真验证为了验证上述方案的正确性,使用MAT L A B进行了仿真验证,由于陀螺和加速度计是分段完成的估计,因此本章也分别对陀螺零偏和加速度计及其标度因数误差进行了仿真,仿真步长统一设置为0.0 0 5s。3.1 陀螺零偏根据快速自标定流程,估计陀螺零偏时,复用了1 5 0s的位置2的数据,再考虑到对准后采集的5 0s数据,因此等价于陀螺零偏估计共用2 0 0s数据,故仿真时长设置为2 0 0s。其他仿真条件设置如下。初始姿态真值:俯仰角,横滚角,航向角=0,0,0 ;初始对准误差:E,N,U=1 0,1 0,5 5 ;陀螺常值零偏:x,y,z=-0.0 1()/h,0.0 3()/h,0.0 2()/h;陀螺 标 度 因 数 误 差:Kgx,Kgy,Kgz=1 01 0-6,81 0-6,-1 11 0-6;陀螺角随机游走:0.0 0 04()/h12,0.0 0 04()/h12,0.0 0 04()/h12;加速度计零偏:x,y,z=3 01 0-6g,-5 0 1 0-6g,4 0 1 0-6g(其 中g=9.8 0 66 5m/s2);加速度计随机游走:11 0-6g/H z12,11 0-6g/H z12,11 0-6g/H z12;按照上述仿真条件进行MAT L A B仿真,仿真结果如图4所示。(a)陀螺零偏估计(b)陀螺零偏估计结果的均方误差图4 陀螺零偏估计及其均方误差仿真结果F i g.4 S i m u l a t i o nr e s u l t so f g y r ob i a sa n di t sm e a n-s q u a r e e r r o r48 第3期一种双轴旋转惯性导航系统的重要参数快速自标定方法从图4中可以看出,进入卡尔曼滤波估计后,陀螺零偏很快收敛并逐渐达到稳定状态,最终x y z三轴陀螺 零 偏 估 计 结 果 依 次 为-0.0 0 94()/h、0.0 2 91()/h、0.0 2 02()/h,与仿真设定初值之差优于0.0 0 1()/h。3.2 加速度计零偏及标度因数误差仿真条件设置如下。仿真时长:1 2 0s;初始姿态真值:俯仰角,横滚角,航向角=0,0,0 ;初始对准误差:E,N,U=1 0,1 0,5 5 ;陀 螺 标 度 因 数 误 差:Kgx,Kgy,Kgz=1 0 1 0-6,81 0-6,-1 11 0-6;陀螺角随机游走:0.0 0 04()/h12,0.0 0 04()/h12,0.0 0 04()/h12;加速度计零偏:x,y,z=3 01 0-6g,-5 0 1 0-6g,4 01 0-6g;加速度计标度因数误差:Kax,Kay,Kaz=1 51 0-6,-1 81 0-6,1 21 0-6;加 速 度 计 随 机 游 走:11 0-6g/H z12,11 0-6g/H z12,11 0-6g/H z12;按照 上 述 仿 真 初 始 条 件 和 标 定 流 程,采 用MAT L A B按表1标定路径生成了六位置翻滚的轨迹,得到的仿真结果如图5和图6所示。从图5和图6中可以看出,进入卡尔曼滤波估计后,加速度计零偏和标度因数误差的估计均方误差曲线均是分段达到收敛状态,这是因为,随着加速度计六位置翻滚,加速度计零偏和其标度因数误差的可观测性逐步达到全部可观测,并最终收敛到稳定状态。三轴加速度计零偏估计结果依次为3 2.2 31 0-6g、-4 8.3 2 1 0-6g、3 8.2 8 1 0-6g,与仿真设定(a)加速度计零偏估计结果(b)加速度计零偏估计结果的均方误差图5 加速度计零偏及其均方误差仿真结果F i g.5 S i m u l a t i o nr e s u l t so fa c c e l e r o m e t e rb i a sa n di t sm e a n-s q u a r e e r r o r(a)加速度计标度因数误差估计结果(b)加速度计标度因数误差估计结果的均方误差图6 加速度计标度因数误差及其均方误差仿真结果F i g.6 S i m u l a t i o nr e s u l t so fa c c e l e r o m e t e r s c a l e f a c t o re r r o ra n d i t sm e a n-s q u a r e e r r o r初值之差优于2 1 0-6g。三轴加速度计标定因数误差估计结果依次为1 5.1 11 0-6、-1 7.5 51 0-6、58 导航定位与授时2 0 2 3年5月1 0.3 9 1 0-6,与仿真设定初值之差优于21 0-6。3.3 陀螺标定因数误差陀螺标定因数误差估计采用数值解析的方式,通过MAT L A B设置转动路径轨迹数据,计算得到陀螺 标 度 因 数 误 差 分 别 为9.6 51 0-6、8.5 71 0-6、-1 0.2 61 0-6、与 仿 真 设 定 值 之 差 均 优 于11 0-6。按照标定流程中的标定路径安排,采用MAT-L A B分别生成相应的仿真轨迹及惯性仪表数据,通过卡尔曼滤波估计及数值解析得到上述结果,仿真结果证明了本文提出的快速自标定方法的正确性和可实现性,并且具有较高的精度。4 试验验证为了验证快速自标定方法在实际应用中的效果,采用某型号双轴激光惯组进行了实物试验验证,按照标定流程设定了电机转动路径和停留时间,同时使用上位机采集软件采集了陀螺和加速度计数据,最后通过离线仿真得到图7图9的仿真结果。(a)陀螺零偏估计(b)陀螺零偏估计结果的均方误差图7 陀螺零偏及其均方误差曲线F i g.7 C u r v eo fg y r ob i a sa n d i t sm e a n-s q u a r e e r r o r(a)加速度计零偏估计结果(b)加速度计零偏估计结果的均方误差图8 加速度计零偏及其均方误差曲线F i g.8 C u r v eo fa c c e l e r o m e t e rb i a sa n d i t sm e a n-s q u a r e e r r o r(a)加速度计标度因数误差估计结果(b)加速度计标度因数误差估计结果的均方误差图9 加速度计标度因数误差及其均方误差曲线F i g.9 C u r v eo fa c c e l e r o m e t e r s c a l e f a c t o re r r o ra n d i t sm e a n-s q u a r e e r r o r68 第3期一种双轴旋转惯性导航系统的重要参数快速自标定方法从图7中可以看出,陀螺零偏估计均方误差曲线在开始滤波后能够快速收敛并达到稳态,陀螺零偏估计曲线均在1 4 0s后进入收敛状态。从图8和图9可以看出,加速度计的零偏和标度因数误差估计均方误差收敛曲线相似,均是随着六位置翻滚,加速度计零偏和标度因数误差逐步达到全部可观测,且均收敛到稳定状态。为了说明标定结果的正确性,采用该惯组进行了6次重复性快速自标定测试,测试结果如表2和表3所示。从表中可知,陀螺零偏标定重复性优于0.0 0 1()/h,陀螺标度因数 误差标定重 复性优于11 0-6,加速度计零偏标定重复性优于21 0-6g,加速度 计标度因 数 误 差 标 定 重 复 性 优 于21 0-6,可见本文所述快速自标定方法正确有效,且有较高的标定精度。表2 陀螺零偏及标度因数误差快速自标定重复性T a b l e2 R e p e a t a b i l i t yo f g y r ob i a sa n ds c a l e e r r o ro f f a s t s e l f-c a l i b r a t i o n测试组别陀螺零偏/()/h陀螺标度因数误差xyzxyz10.0 1 50.0 2 10.0 2 12.91 0-65.81 0-6-3.51 0-620.0 1 70.0 2 20.0 2 13.11 0-65.71 0-6-3.41 0-630.0 1 60.0 2 10.0 2 23.01 0-65.91 0-6-3.31 0-640.0 1 50.0 2 00.0 2 03.11 0-65.61 0-6-3.41 0-650.0 1 80.0 2 30.0 2 23.21 0-65.51 0-6-3.61 0-660.0 1 40.0 2 20.0 2 33.31 0-65.41 0-6-3.21 0-6重复性0.0 0 10.0 0 10.0 0 10.1 41 0-60.1 91 0-60.1 41 0-6表3 加速度计零偏及标度因数误差快速自标定重复性T a b l e3 R e p e a t a b i l i t yo fa c c e l e r o m e t e rb i a sa n ds c a l ee r r o ro f f a s t s e l f-c a l i b r a t i o n测试组别加速度计零偏/(ms-2)加速度计标度因数误差xyzxyz1-4 1.91 0-6g-6 7.71 0-6g-2 2.71 0-6g3 3.21 0-6-1 1 5.11 0-6-1 2 11 0-62-4 1.61 0-6g-6 6.41 0-6g-2 31 0-6g3 5.11 0-6-1 1 8.21 0-6-1 2 2.21 0-63-4 2.81 0-6g-6 5.61 0-6g-2 3.61 0-6g3 2.51 0-6-1 1 7.51 0-6-1 2 0.81 0-64-4 2.51 0-6g-6 8.11 0-6g-2 1.21 0-6g3 2.31 0-6-1 1 3.61 0-6-1 2 1.91 0-65-4 3.11 0-6g-6 5.81 0-6g-2 2.21 0-6g3 1.61 0-6-1 1 4.31 0-6-1 2 31 0-66-4 1.31 0-6g-6 7.31 0-6g-2 3.11 0-6g3 0.21 0-6-1 1 6.51 0-6-1 1 7.71 0-6重复性0.7 11 0-6g1.0 31 0-6g0.8 41 0-6g1.6 41 0-61.8 31 0-61.8 51 0-6 通过仿真和试验验证,表明本文所述的基于双轴旋转惯导的重要参数快速自标定方法可以实现1 0m i n内完成陀螺和加速度计的零偏及标度因数误差的估计,且具有较高的精度。5 结论本文从实际使用的角度出发,提出了一种基于双轴旋转惯导系统的快速自标定方法。通过合理地安排标定路径,分3步依次标定出陀螺零偏、加速度计标度因数、加速度计零偏和陀螺标度因数共1 2个重要参数。数值仿真和实际系统试验结果表明,采用本文所述快速自标定方法,1 0m i n内可以完成陀螺和加速度计1 2个重要参数的标定,相比于传统的标定方法,大幅缩短了标定时间,且具有较高的标定精度。该方法在保证现场标定精度的前提下,实现了现场的快速自标定,有助于缩短武器系统准备时间和快速打击,具有较高的工程应用价值。参考文献1 董春梅.基于误差分离与抑制技术的惯性仪表标定方法D.哈尔滨:哈尔滨工业大学,2 0 1 9.D o n gC h u n m e i.I n e r t i a l i n s t r u m e n t c a l i b r a t i o nm e t h o db a s e do ne r r o rs e p a r a t i o na n ds u p p r e s s i o nt e c h n o l o g yD.H a r b i n:H a r b i nI n s t i t u t eo fT e c h n o l o g y,2 0 1 9(i nC h i n e s e).78 导航定位与授时2 0 2 3年5月2L uJ,L e iC,L i a n gS.e ta l.A na l l-p a r a m e t e rs y s t e m-l e v e l c a l i b r a t i o nf o rs t e l l a r-i n e r t i a ln a v i g a t i o ns y s t e mo ng r o u n dJ.I E E ET r a n s a c t i o n so nI n s t r u m e n t a t i o n&M e a s u r e m e n t,2 0 1 7,6 6(8).2 0 6 5-2 0 7 3.3W a n g B,R e n Q,D e n g Z.e t a l.A s e l f-c a l i b r a t i o nm e t h o df o rn o n o r t h o g o n a la n g l e sb e t w e e ng i m b a l so fr o t a t i o n a l i n e r t i a l n a v i g a t i o ns y s t e mJ.I E E ET r a n s-a c t i o n so nI n d u s t r i a lE l e c t r o n i c s,2 0 1 5,6 2(4).4W a n gL,W a n g W,Z h a n g Q.e ta l.S e l f-c a l i b r a t i o nm e t h o db a s e do nn a v i g a t i o ni nh i g h-p r e c i s i o ni n e r t i a ln a v i g a t i o ns y s t e m w i t hf i b e ro p t i cg y r oJ.O p t i c a lE n g i n e e r i n g,2 0 1 4,5 3(6).5G a oW,Z h a n gY,W a n gJ.R e s e a r c ho ni n i t i a la l i g n-m e n t a n ds e l f-c a l i b r a t i o no fr o t a r ys t r a p d o w ni n e r t i a ln a v i g a t i o ns y s t e m s.J.S e n s o r s(B a s e l,S w i t z e r-l a n d),2 0 1 5,1 5(2).3 1 5 4-7 1.6 张金刚,姜述明.双轴旋转惯导系统误差自标定技术研究J.导航定位与授时,2 0 1 9,6(6):1 2 4-1 3 0.Z h a n gJ i n g a n g,J i a n gS h u m i n g.R e s e a r c ho ns e l f-c a l i-b r a t i o nt e c h n o l o g yo ft w o-a x i sr o t a t i o ni n e r t i a ln a v i-g a t i o ns y s t e mJ.N a v i g a t i o nP o s i t i o n i n g&T i m i n g,2 0 1 9,6(6):1 2 4-1 3 0(i nC h i n e s e).7 谢波,秦永元,万彦辉.激光陀螺捷联惯导系统多位置标定方法J.中 国惯 性 技术 学报,2 0 1 1,1 9(2):1 5 7-1 6 2.X i e B o,Q i n Y o n g y u a n,W a n Y a n h u i.M u l t i p o s i t i o nc a l i b r a t i o nm e t h o do f l a s e rg y r oS I N SJ.J o u r n a lo fC h i n e s e I n e r t i a lT e c h n o l o g y,2 0 1 1,1 9(2):1 5 7-1 6 2(i nC h i n e s e).8 杨明.车载捷联惯导系统标定技术研究D.哈尔滨:哈尔滨工程大学,2 0 1 7.Y a n g M i n g.R e s e a r c h o n c a l i b r a t i o n t e c h n o l o g y o fs t r a p d o w ni n e r t i a ln a v i g a t i o ns y s t e m o nt h ev e h i c l eD.H a r b i n:H a r b i nE n g i n e e r i n g U n i v e r s i t y,2 0 1 7(i nC h i n e s e).9G a oP,L iK,W a n gL,e ta l.As e l f-c a l i b r a t i o nm e t h o df o ra c c e l e r o m e t e rn o n l i n e a r i t ye r r o r si nt r i a x i sr o t a-t i o n a l i n e r t i a ln a v i g a t i o ns y s t e mJ.I E E E T r a n s a c-t i o n so nI n s t r u m e n t a t i o n&M e a s u r e m e n t,2 0 1 7,6 6(2).2 4 3-2 5 3.1 0 袁鹏.激光陀螺捷联惯导系统动基座标定技术研究D.北京:中国航天科技集团公司第一研究院,2 0 1 8.Y u a nP e n