本发明专利技术公开了一种加速度计二次项误差的系统级拟合标定方法，它有八大步骤：将惯组安装到转台上，惯组初始朝向为地—东—南，采样周期dt=0.01s；第一位置旋转，采集惯组数据，旋转完毕静置一分钟，停止采数；第二位置旋转，采集惯组数据，旋转完毕静置一分钟，停止采数；第三位置旋转，将惯组安装到转台上，惯组初始朝向为西—北—地，采集惯组数据，旋转完毕静置一分钟，停止采数；第四位置旋转，将惯组安装到转台上，三个坐标轴分别朝向北—天—东，采集惯组数据，旋转完毕静置一分钟，停止采数，至此四组旋转完毕；对以上各组惯组采集的数据分别进行导航解算，计算各轴向速度误差；设计卡尔曼滤波器；解析计算，求出加速度计的二次项系数误差。
【技术实现步骤摘要】
本专利技术涉及，属于惯性导航
技术介绍
标定技术是惯性导航领域的核心技术之一，是一种对误差的辨识技术，即建立惯性元件和惯导系统的误差数学模型，通过一系列的试验求解出误差模型中的误差项，进而通过软件算法来对误差进行补偿。标定可以从软件方面来提高惯导系统的精度，常见的标定误差项有加速度计和陀螺的零偏误差、标度因数误差、安装误差等。按层次来分，标定可分为分立式标定和系统级标定。当前分立式标定方法的理论研究已较为成熟，而系统级标定方法是由20世纪80年代开始发展起来，目前正成为标定技术研究的热点。分立式标定是直接利用加速度计和陀螺仪的输出作为观测量来对其误差参数进行辨识。分立式标定通过设计合适的旋转编排路径以抵消地球自转和重力加速度对惯导系统标定的影响。该方法标定精度受随机噪声影响，通常需对同一系统进行多次测试取其平均值。系统级标定是利用惯性导航输出的系统误差作为观测量来辨识惯性器件的误差参数。与分立式标定相比，系统级标定的算法相对复杂，但系统级标定能抑制标定过程中的测量噪声，并缩短标定时间，提高标定精度。系统级标定又分为拟合标定和滤波标定。拟合标定是通过建立特定运动激励观测导航误差，拟合估计器件误差参数。拟合标定方法一般将导航速度变化率或比力作为观测量。滤波标定是指通过设计Kalman滤波器，将惯导系统各器件误差参数作为滤波器状态，观测导航误差滤波估计各器件误差参数。杨杰，吴文启所写的“激光陀螺捷联系统高精度加速度计非线性模型参数标定”(见中国惯性技术学报2010，18 (5))提出基于角位置精度为5"的三轴转台，应用最小二乘分步辨识的思想，通过正交六位置辨识出加速度计的二次项误差。本方法属于分立式标定方法。张红良写的“陆用高精度激光陀螺捷联惯导系统误差参数估计方法研究”.长沙:国防科技大学，2010，10.分别对拟合标定方法和滤波标定方法进行了说明，但该论文所述标定方法中没有考虑加速度计的二次项误差。
技术实现思路
1、目的: 本专利技术的目的是提供了，能够准确地估计出加速度计的二次项误差。2、技术方案:本专利技术，该方法具体步骤如下:步骤一:将惯组安装到转台上，惯组初始朝向为地一东一南。设置陀螺输出已预先标定并进行补偿。惯组通电预热三十分钟，采样周期dt = 0.0ls0步骤二:第一位置即地-东-南旋转。开始采集惯组数据。首先静置一分钟，然后使惯组以10° /s的角速度绕Y轴正向旋转180°，旋转完毕再静置一分钟，停止采数。步骤三:第二位置即天-东-北旋转。开始采集惯组数据，首先使惯组在步骤二旋转完毕位置静置一分钟，然后使转台以10° /s的角速度绕惯组Z轴正向旋转180°，旋转完毕再静置一分钟。停止采数。步骤四:第三位置即西-北-地旋转。将惯组安装到转台上，惯组初始朝向为西一北一地。开始采集惯组数据，首先静置一分钟，然后使转台以10° /s的角速度绕惯组Y轴正向旋转180°，旋转完毕再静置一分钟。停止采数。步骤五:第四位置即北-天-东旋转。将惯组安装到转台上，三个坐标轴分别朝向北一天一东。开始采集惯组数 据，首先静置一分钟，然后使转台以10° /s的角速度绕惯组Z轴正向旋转180°，旋转完毕再静置一分钟。停止采数，至此四组旋转完毕。步骤六:对以上各组惯组采集的数据分别进行导航解算，并计算各轴向速度误差。将各组惯组数据采集对应的起始时刻记为Ttl，转动前静置一分钟的末时刻记为T1，转动完静置一分钟的末时刻记为T2。由于惯组只有角运动没有线运动，系统速度误差可用下式表示:Sv1jQ) = v'j (t)其中i = 1，2，3，4,表示第i位置旋转；j = x, y, z,表示对应j轴向的速度误差。步骤七:设计卡尔曼滤波器。状态变量为权利要求1.，其特征在于:该方法具体步骤如下: 步骤一:将惯组安装到转台上，惯组初始朝向为地一东一南；设置陀螺输出已预先标定并进行补偿，惯组通电预热三十分钟，采样周期dt = 0.0ls ； 步骤二:第一位置即地-东-南旋转，开始采集惯组数据，首先静置一分钟，然后使惯组以10° /s的角速度绕Y轴正向旋转180°，旋转完毕再静置一分钟，停止采数； 步骤三:第二位置即天-东-北旋转，开始采集惯组数据，首先使惯组在步骤二旋转完毕位置静置一分钟，然后使转台以10° /s的角速度绕惯组Z轴正向旋转180°，旋转完毕再静置一分钟，停止采数； 步骤四:第三位置即西-北-地旋转，开始采集惯组数据，首先静置一分钟，然后使转台以10° /s的角速度绕惯组Y轴正向旋转180°，旋转完毕再静置一分钟；停止采数；步骤五:第四位置即北-天-东旋转，开始采集惯组数据，首先静置一分钟，然后使转台以10° /s的角速度绕惯组Z轴正向旋转180°，旋转完毕再静置一分钟，停止采数，至此四组旋转完毕； 步骤六:对以上各组惯组采集的数据分别进行导航解算，并计算各轴向速度误差，将各组惯组数据采集对应的起始时刻记为Ttl，转动前静置一分钟的末时刻记为T1，转动完静置一分钟的末时刻记为T2,由于惯组只有角运动没有线运动，系统速度误差用下式表示:Jv；(/) = v；(/) 其中i = 1，2，3，4,表示第i位置旋转；j = X，y, z,表示对应j轴向的速度误差； 步骤七:设计卡尔曼滤波器:状态变量为[而Uo Sv1^t)抑;(0了；使用标准卡尔曼方程进行迭代，滤波完毕得到力彳⑴； 系统状态方程:全文摘要本专利技术公开了，它有八大步骤将惯组安装到转台上，惯组初始朝向为地—东—南，采样周期dt=0.01s；第一位置旋转，采集惯组数据，旋转完毕静置一分钟，停止采数；第二位置旋转，采集惯组数据，旋转完毕静置一分钟，停止采数；第三位置旋转，将惯组安装到转台上，惯组初始朝向为西—北—地，采集惯组数据，旋转完毕静置一分钟，停止采数；第四位置旋转，将惯组安装到转台上，三个坐标轴分别朝向北—天—东，采集惯组数据，旋转完毕静置一分钟，停止采数，至此四组旋转完毕；对以上各组惯组采集的数据分别进行导航解算，计算各轴向速度误差；设计卡尔曼滤波器；解析计算，求出加速度计的二次项系数误差。文档编号G01C25/00GK103234554SQ201310117069公开日2013年8月7日 申请日期2013年4月7日 优先权日2012年5月7日专利技术者芦佳振, 张春熹, 李保国, 汤卓 申请人:北京航空航天大学本文档来自技高网...
【技术保护点】
一种加速度计二次项误差的系统级拟合标定方法，其特征在于：该方法具体步骤如下：步骤一：将惯组安装到转台上，惯组初始朝向为地—东—南；设置陀螺输出已预先标定并进行补偿，惯组通电预热三十分钟，采样周期dt＝0.01s；步骤二：第一位置即地?东?南旋转，开始采集惯组数据，首先静置一分钟，然后使惯组以10°/s的角速度绕Y轴正向旋转180°,旋转完毕再静置一分钟，停止采数；步骤三：第二位置即天?东?北旋转，开始采集惯组数据，首先使惯组在步骤二旋转完毕位置静置一分钟，然后使转台以10°/s的角速度绕惯组Z轴正向旋转180°,旋转完毕再静置一分钟，停止采数；步骤四：第三位置即西?北?地旋转，开始采集惯组数据，首先静置一分钟，然后使转台以10°/s的角速度绕惯组Y轴正向旋转180°,旋转完毕再静置一分钟；停止采数；步骤五：第四位置即北?天?东旋转，开始采集惯组数据，首先静置一分钟，然后使转台以10°/s的角速度绕惯组Z轴正向旋转180°,旋转完毕再静置一分钟，停止采数，至此四组旋转完毕；步骤六：对以上各组惯组采集的数据分别进行导航解算，并计算各轴向速度误差，将各组惯组数据采集对应的起始时刻记为T0，转动前静置一分钟的末时刻记为T1，转动完静置一分钟的末时刻记为T2，由于惯组只有角运动没有线运动，系统速度误差用下式表示：δvji(t)=vji(t)其中i＝1,2,3,4，表示第i位置旋转；j＝x,y,z，表示对应j轴向的速度误差；步骤七：设计卡尔曼滤波器：状态变量为δvji(t)δv·ji(t)δv··ji(t)T;使用标准卡尔曼方程进行迭代，滤波完毕得到系统状态方程：ddtδvji(t)δv·ji(t)δv··ji(t)=010001000δvji(t)δv·ji(t)δv··ji(t)+00wj观测方程：Zj=100δvji(t)δv·ji(t)δv··ji(t)+μj其中wj和μj分别为状态噪声和观测噪声；步骤八：解析计算，求出加速度计的二次项系数误差；首先令误差观测方程如下表1所示：表1?误差观测方程则二次项误差解算如下表2所示：表2?二次项误差解算表误差项X轴Y轴Z轴加表二次项误差K2ax＝(c1?a3)/2g2K2ay＝(a2?c4)/2g2K2az＝(c3?b1)/2g2表中的符号说明如下：K2ax表示x轴加表二次项误差，K2ay表示y轴加表二次项误差，K2az表示z轴加表二次项误差；K2ax、K2ay、K2az是目标误差项；g表示重力加速度，ai、bi、ci是求解目标误差项过程中的误差观测量，通过ai、bi、ci解析计算出K2ax、K2ay、K2az的详细求法见上表1，i=1,2,3,4。FDA00003014089800013.jpg,FDA00003014089800021.jpg,FDA00003014089800022.jpg,FDA00003014089800023.jpg...
【技术特征摘要】
...
【专利技术属性】
技术研发人员：芦佳振，张春熹，李保国，汤卓，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。