[发明专利]一种适合车载光纤航姿系统的姿态测量方法无效
| 申请号: | 201010574097.9 | 申请日: | 2010-12-06 |
| 公开(公告)号: | CN102095424A | 公开(公告)日: | 2011-06-15 |
| 发明(设计)人: | 陈贵金;王勇刚 | 申请(专利权)人: | 国营红峰机械厂 |
| 主分类号: | G01C21/28 | 分类号: | G01C21/28;G01S19/49 |
| 代理公司: | 北京市德权律师事务所 11302 | 代理人: | 周发军 |
| 地址: | 432000 *** | 国省代码: | 湖北;42 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | 本发明涉及属于惯性导航技术领域,具体涉及一种适合车载光纤航姿系统的姿态测量方法。本发明采用基于卡尔曼滤波器的低成本光纤航姿系统与GPS组合导航系统,在SINS初始对准前,对GPS故障信号进行检测与隔离,保证SINS/GPS组合系统的可靠性,当GPS输出有效后,选取SINS输出的航向角、北向、东向速度与GPS接收机输出相应信息的差值作为观测量,构成3维观测方程。本发明能够大大加快计算速度,实现实时滤波计算,而且与速度+位置作为观测量的组合系统相比,其导航精度尤其是姿态计算机的精度有很大的提高。 | ||
| 搜索关键词: | 一种 适合 车载 光纤 系统 姿态 测量方法 | ||
【主权项】:
1.一种适合车载光纤航姿系统的姿态测量方法,采用基于卡尔曼滤波器的低成本光纤航姿系统与GPS组合导航系统,其特征在于,按照系统解算流程,当GPS输出有效后,利用GPS输出的姿态、速度和SINS输出的姿态、速度之差作为观测量,具体包括如下步骤:首先建立以SINS误差方程为基础的组合导航系统状态方程,并在误差方程的基础上建立组合系统的量程方程,用一阶线性随机微分方程来描述系统的误差状态传播特性和系统的量测误差状态,方程如下:X · ( t ) = F ( t ) X ( t ) + G ( t ) W ( t ) - - - ( 1 ) ]]> Z(t)=H(t)X(t)+V(t)……………………………………(2)选取SINS导航参数计算误差、光纤陀螺仪以及加速度计误差作为系统的状态变量,共12维,即:X T = [ φ E , φ N , φ U , δV E , δV N , δL , δλ , ϵ bx , ϵ by , ϵ bz , ▿ bx , ▿ by ] - - - ( 3 ) ]]> 式中:φE,φN,φU为SINS姿态计算误差,δVE和δVN为SINS的速度误差;δL和δλ为S INS位置误差;εbx、εby和εbz为3个轴向陀螺漂移误差![]()
为加速度计测量误差;选取SINS输出的航向角、北向、东向速度与GPS接收机输出相应信息的差值作为观测量,构成3维观测方程,即:Z ( t ) = Z v Z ψ = H v ( t ) H ψ ( t ) X ( t ) + V v ( t ) V ψ ( t ) - - - ( 4 ) ]]> 式中,HV(t)是观测量为速度时的量测矩阵,Hψ(t)是观测量为姿态时的量测矩阵,其中速度量测矢量为:Z v = v SINS . N - v GPS . N v SINSE - v GPS . E = H v X v + V v - - - ( 5 ) ]]> 式中,vSINS.N、vSINS.E是SINS输出的北向和东向速度信息,vSINS.N=vN+δvN,vSINS.E=vE+δvE,其中vN、vE是载体沿地理坐标系统北向和东向的真实速度;vGPS.N、vGPS.E是GPS接收机输出的北向和东向的速度信息,vGPS.N=vN-MN、vGPS.E=vE-ME,其中MN、ME是GPS接收机测速误差,设为白噪声,Hv是观测量为速度时的量测矩阵;Hv=[diag[1 1]02×10],Vv是GPS接收机输出速度的量测噪声:Vv=[MN ME]T,取量测噪声为白噪声,Zψ=ψSINS.yaw-ψGPS.yaw=HψXψ+Vψ……………………………(6)式中,ψSINS.yaw为SINS输出的航向角,ψSINS.yaw=ψaw+δψSINS.yaw,其中ψaw为真实航向角,δψSINS.yaw为SINS航向角测量误差;ψGPS.yaw为双GPS输出的航向角,ψGPS.yaw=ψyaw+δψGPS.yaw,其中δψGPS.yaw是双GPS航向角测量误差。视双GPS接收机测量精度和基线长度而定;Hψ是观测量为姿态时的量测矩阵,Hψ=[01×2 1 01×9];Vψ为双GPS输出航向角的量测噪声,Vψ=δψGPS.yaw,取量测噪声为白噪声,设GPS航向测量白噪声方差为
则Rψ=(HDOP.δψ)2离散后的状态方程和观测方程为:XK=Φk,k-1Xk-1+Γk-1Wk-1……………………………….…(8)ZK=HkKk+Vk卡尔曼滤波算法如下,状态下一步预测:
状态估计:
滤波增益:
下一步预测均方差:
估计均方差:P k = ( I - K k H k ) P k / k - 1 ( I - K k H k ) T + K k R k K k T . ]]>
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于国营红峰机械厂,未经国营红峰机械厂许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201010574097.9/,转载请声明来源钻瓜专利网。





