一种基于双轴旋转的激光陀螺仪组合误差系数分离方法技术

一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，首先将捷联惯性组合安装在双轴旋转机构上；由捷联惯性组合及双轴旋转机构共同组成双轴旋转惯性系统；其次推导双轴旋转机构的运动学方程；建立双轴旋转惯性系统中陀螺仪组合的静基座误差模型，进而得出双轴旋转机构的角速度误差方程；调整双轴旋转机构的基座，使得基座坐标系与地理坐标系重合；按照十六位置转位方法依次旋转双轴旋转机构；采集转停全过程中陀螺仪组合的输出；最终对采集得到的数据使用卡尔曼滤波方法实现陀螺仪组合的误差系数的分离。本发明专利技术适用于双轴旋转惯性系统在静基座条件下自主标校捷联惯性组合中陀螺仪组合误差系数的场合，避免了系统定期拆卸带来的繁琐、复杂的操作。

【技术实现步骤摘要】
一种基于双轴旋转的激光陀螺仪组合误差系数分离方法
本专利技术涉及一种误差系数标定方法，尤其涉及一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，属于捷联惯性组合标定技术，可用于标定捷联惯性组合中陀螺仪组合的场合。
技术介绍
陀螺仪是惯性导航和惯性制导系统的基本测量元件之一，它安装在运载体内部，用于测量运载体的运动角速度，并通过对角速度的积分，求得载体运动的角度。陀螺仪组合的性能和精度直接影响导航和制导系统的精度，对惯性系统的性能起着关键的作用。捷联惯性组合是将陀螺仪组合和加速度计组合集成在一起并直接安装在运载体上的惯性测量装置。捷联惯性组合在结构安装，惯性仪表以及系统的工程实现中各个环节都不可避免的存在误差。通常由于存储、气候环境变化、运输及捷联惯性组合自身性能不稳定等因素的影响，组合误差模型的某些误差系数会发生变化。解决的方法有两种：一是改进设计，工艺及提高元器件的性能指标，对影响惯性系统精度的陀螺仪、加速度计、再平衡回路等采取相应的措施使其在寿命周期内的变化小到可以满足要求；二是对系统进行定期标定。实践证明，单纯通过第一种方法来提高惯性导航系统的精度变得越来越困难，成本也越来越高，而采用第二种来提高系统的精度要容易得多，这是一种成本低、见效快的方法。目前通常采取的是第二种方法，即对其进行重新标定。传统的重新标定方法是每隔一段时间后将惯性系统从运载体上拆卸下来，需要大型的精密测试设备和专业人员的参与，工作量大、成本高、整个过程繁琐、费时。在民用领域，其高成本影响了惯性系统的推广应用，在军事领域，笨重的地面设备和复杂的操作限制了武器的机动性，制约了武器性能的提高，不利于部队在技术阵地和战备值班时的日常保障和快速反应。在不拆卸惯导的前提下，为了确保系统的对准和导航精度，提高外场标定效率，需要研究一种新型的基于双轴旋转的陀螺仪组合误差系数的分离方法。
技术实现思路
本专利技术解决的技术问题是：克服现有技术的不足，提供一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，实现了惯性系统的定期自标定。本专利技术的技术解决方案是：一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，步骤如下：1)将捷联惯性组合安装在双轴旋转机构上；由捷联惯性组合及双轴旋转机构共同组成双轴旋转惯性系统；2)推导双轴旋转机构的运动学方程；建立双轴旋转惯性系统中陀螺仪组合的静基座误差模型，进而得出双轴旋转机构的角速度误差方程；3)调整双轴旋转机构的基座，使得基座坐标系与地理坐标系重合；按照十六位置转位方法依次旋转双轴旋转机构，旋转过程中，在十六个位置中的每一个位置处停留时间大于60s；采集转停全过程中陀螺仪组合的输出；4)对采集得到的数据使用卡尔曼滤波方法实现陀螺仪组合的误差系数的分离。步骤2)中双轴旋转机构的运动学方程为：式中ωxp、ωyp、ωzp分别为陀螺仪组合绕本体坐标系的三个轴X、Y、Z的绝对角速度分量；Ωn＝ωiecosL，Ωn为地球转速投影到地理坐标系中的北向分量；Ωu＝ωiesinL，Ωu为地球转速投影到地理坐标系中的天向分量；ωie为地球自转角速度；L为对捷联惯性组合进行测试的地点的地球纬度；α为双轴旋转机构内环轴的旋转角；β为双轴旋转机构外环轴的旋转角；为双轴旋转机构内环轴的角速度，为双轴旋转机构外环轴的角速度；双轴旋转惯性系统中陀螺仪组合的静基座误差模型为：式中Δωxp、Δωyp、Δωzp为陀螺仪组合的测量误差；ωx、ωy、ωz为陀螺仪组合经补偿得到的角速度输出值；δD0x、δD0y和δD0z为陀螺仪组合的零次项偏差；δkgx、δkgy和δkgz为陀螺仪组合的标度因数偏差；δEYX为X轴相对于Y轴的安装误差角偏差；δEZX为X轴相对于Z轴的安装误差角偏差；δEXY为Y轴相对于X轴的安装误差角偏差；δEZY为Y轴相对于Z轴的安装误差角偏差；δEXZ为Z轴相对于X轴的安装误差角偏差；δEYZ为Z轴相对于Y轴的安装误差角偏差；双轴旋转机构的角速度误差方程为：式中为双轴旋转机构内环轴的角速度误差，为双轴旋转机构外环轴的角速度误差；Δα为双轴旋转机构内环轴的旋转角误差；Δβ为双轴旋转机构外环轴的旋转角误差。其特征在于步骤3)中十六位置转位方法的位置分别为：位置1：基座坐标系与地理坐标系重合时，记录此时的内环轴旋转角α(1)和外环轴旋转角β(1)；位置2：调节旋转机构使内环轴旋转角度α(2)＝α(1)+90°，外环轴旋转角度β(2)＝β(1)；位置3：调节旋转机构使内环轴旋转角度α(3)＝α(1)+180°，外环轴旋转角度β(3)＝β(1)；位置4：调节旋转机构使内环轴旋转角度α(4)＝α(1)+270°，外环轴旋转角度β(4)＝β(1)；位置5：调节旋转机构使内环轴旋转角度α(5)＝α(1)，外环轴旋转角度β(5)＝β(1)+180°；位置6：调节旋转机构使内环轴旋转角度α(6)＝α(1)+90°，外环轴旋转角度β(6)＝β(1)+180°；位置7：调节旋转机构使内环轴旋转角度α(7)＝α(1)+180°，外环轴旋转角度β(7)＝β(1)+180°；位置8：调节旋转机构使内环轴旋转角度α(8)＝α(1)+270°，外环轴旋转角度β(8)＝β(1)+180°；位置9：调节旋转机构使内环轴旋转角度α(9)＝α(1)+90°，外环轴旋转角度β(9)＝β(1)+90°；位置10：调节旋转机构使内环轴旋转角度α(10)＝α(1)+270°，外环轴旋转角度β(10)＝β(1)+90°；位置11：调节旋转机构使内环轴旋转角度α(11)＝α(1)，外环轴旋转角度β(11)＝β(1)+90°；位置12：调节旋转机构使内环轴旋转角度α(12)＝α(1)+180°，外环轴旋转角度β(12)＝β(1)+90°；位置13：调节旋转机构使内环轴旋转角度α(13)＝α(1)+270°，外环轴旋转角度β(13)＝β(1)+270°；位置14：调节旋转机构使内环轴旋转角度α(14)＝α(1)+90°，外环轴旋转角度β(14)＝β(1)+270°；位置15：调节旋转机构使内环轴旋转角度α(15)＝α(1)+180°，外环轴旋转角度β(15)＝β(1)+270°；位置16：调节旋转机构使内环轴旋转角度α(16)＝α(1)，外环轴旋转角度β(16)＝β(1)+270°。所述步骤4)中卡尔曼滤波器状态方程如下其中：X＝[δD0x,δD0y,δD0z,δkgx,δkgy,δkgz,δEYX,δEZX,δEXY,δEZY,δEXZ,δEYZ,Δα,Δβ]；为tk-1时刻到tk时刻的一步转移矩阵，Wk为白噪声；k为正整数；F1＝[01×2,1,01×2,ωz,01×4,ωx,ωy,0,Ωncosβ+Ωusinβ]F2＝[cosα,-sinα,0,ωxcosα,-ωysinα,0,ωycosα,ωzcosα,-ωxsinα,-ωzsinα,01×2,-(ωxpsinα+ωypcosα),0]卡尔曼滤波器量测方程如下：Zk＝HkXk+Vk其中，Hk为量测矩阵，Vk为量测噪声序列；当双轴旋转机构处于静止位置时，量测向量为：Zk_gyro_1＝[Δα,Δβ]量测矩阵为：Hk_gyro_1＝[02×12,I2]I表示单位矩阵；双轴旋转机构处于转动过程时，量测向量为：Zk_gyro_2＝[δ本文档来自技高网...

【技术保护点】
一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，其特征在于步骤如下：1)将捷联惯性组合安装在双轴旋转机构上；由捷联惯性组合及双轴旋转机构共同组成双轴旋转惯性系统；2)推导双轴旋转机构的运动学方程；建立双轴旋转惯性系统中陀螺仪组合的静基座误差模型，进而得出双轴旋转机构的角速度误差方程；3)调整双轴旋转机构的基座，使得基座坐标系与地理坐标系重合；按照十六位置转位方法依次旋转双轴旋转机构，旋转过程中，在十六个位置处停留时间大于60s；采集转停全过程中陀螺仪组合的输出；4)对采集得到的数据使用卡尔曼滤波方法实现陀螺仪组合的误差系数的分离。

【技术特征摘要】
1.一种基于双轴旋转的激光陀螺仪组合误差系数分离方法，其特征在于步骤如下：1)将捷联惯性组合安装在双轴旋转机构上；由捷联惯性组合及双轴旋转机构共同组成双轴旋转惯性系统；2)推导双轴旋转机构的运动学方程；建立双轴旋转惯性系统中陀螺仪组合的静基座误差模型，进而得出双轴旋转机构的角速度误差方程；所述双轴旋转机构的运动学方程为：式中ωxp、ωyp、ωzp分别为陀螺仪组合绕本体坐标系的三个轴X、Y、Z的绝对角速度分量；Ωn＝ωiecosL，Ωn为地球转速投影到地理坐标系中的北向分量；Ωu＝ωiesinL，Ωu为地球转速投影到地理坐标系中的天向分量；ωie为地球自转角速度；L为对捷联惯性组合进行测试的地点的地球纬度；α为双轴旋转机构内环轴的旋转角；β为双轴旋转机构外环轴的旋转角；为双轴旋转机构内环轴的角速度，为双轴旋转机构外环轴的角速度；双轴旋转惯性系统中陀螺仪组合的静基座误差模型为：式中Δωxp、Δωyp、Δωzp为陀螺仪组合的测量误差；ωx、ωy、ωz为陀螺仪组合经补偿得到的角速度输出值；δD0x、δD0y和δD0z为陀螺仪组合的零次项偏差；δkgx、δkgy和δkgz为陀螺仪组合的标度因数偏差；δEYX为X轴相对于Y轴的安装误差角偏差；δEZX为X轴相对于Z轴的安装误差角偏差；δEXY为Y轴相对于X轴的安装误差角偏差；δEZY为Y轴相对于Z轴的安装误差角偏差；δEXZ为Z轴相对于X轴的安装误差角偏差；δEYZ为Z轴相对于Y轴的安装误差角偏差；双轴旋转机构的角速度误差方程为：式中为双轴旋转机构内环轴的角速度误差，为双轴旋转机构外环轴的角速度误差；Δα为双轴旋转机构内环轴的旋转角误差；Δβ为双轴旋转机构外环轴的旋转角误差；3)调整双轴旋转机构的基座，使得基座坐标系与地理坐标系重合；按照十六位置转位方法依次旋转双轴旋转机构，旋转过程中，在十六个位置中的每一个位置处停留时间大于60s；采集转停全过程中陀螺仪组合的输出；4)对采集得到的数据使用卡尔曼滤波方法实现陀螺仪组合的误差系数的分离；卡尔曼滤波器状态方程如下：其中：X＝[δD0x,δD0y,δD0z,δkgx,δkgy,δkgz,δEYX,δEZX,δEXY,δEZY,δEXZ,δEYZ,Δα,Δβ]；为tk-1时刻到tk时刻的一步转移矩阵，Wk为白噪声；k为正整数；F1＝[01×2,1,01×2,ωz,01×4,ωx,ωy,0,Ωncosβ+Ωusinβ]F2＝[cosα,-sinα,0,ωxcosα,-ωysinα,0,ωycosα,ωzcosα,-ωxsinα,-ωzsinα,01×2,-(ωxpsinα+ωypcosα),0]卡尔曼滤波器量...

【专利技术属性】
技术研发人员：魏宗康，黄超，
申请(专利权)人：北京航天控制仪器研究所，
类型：发明
国别省市：北京;11