发明名称 基于地形轮廓匹配的惯性导航系统速度累积误差修正方法
摘要 本发明涉及基于地形轮廓匹配的惯性导航系统速度累积误差修正方法,主要解决现有地形轮廓匹配方法不能修正惯性导航系统速度误差的问题。其实现步骤是:1)根据惯性导航系统指示在地形高程图中找到航行器的位置;2)根据航行器的位置,获取实测高程序列的起始坐标;3)根据实测高程序列的起始坐标,获取实测地形高程序列,并构造地形匹配起始点搜索区;4)根据地形匹配起始点搜索区的范围读取搜索区高程序列,计算搜索区各高程序列与实测高程序列的平方差,得到地形匹配位置坐标;6)根据实测高程序列起始坐标和地形匹配位置坐标对惯性导航系统进行速度修正。本发明减小了惯性导航系统的速度误差和地形匹配位置误差,可用于地形匹配辅助导航。
申请公布号 CN106052688A 申请公布日期 2016.10.26
申请号 CN201610644276.2 申请日期 2016.08.08
申请人 西安电子科技大学;陕西长岭电子科技有限责任公司 发明人 周宇;万俊;张林让;杨慧婷;蒙妍
分类号 G01C21/16(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华
主权项 基于地形轮廓匹配的惯性导航系统速度累积误差修正方法,包括:(1)读取当前机载惯性导航系统指示的航行器位置坐标,根据此位置坐标,在机载数字地形高程图中找到当前航行器的位置;(2)判断航行器是否处在地形匹配区域中:若不在地形匹配区域中,则返回(1);若处在地形匹配区域中,则调整航行器的航行姿态和速度v进行匀速直线航行,并把当前惯性导航系统指示位置的坐标作为实测高程序列的起始坐标;(3)根据实测高程序列的起始坐标,获取实测地形高程序列;(4)根据实测高程序列的起始坐标,构造地形匹配起始点搜索区;(5)判断地形起始点匹配搜索区是否完全在地形匹配区域中:若不完全在该地形匹配区域中,则不能进行本次地形匹配;若完全在地形匹配区域中,则执行(6);(6)分别计算搜索区各高程序列和实测高程序列的平方差值,并比较各平方差值,得到地形匹配位置坐标;(7)根据实测高程序列起始坐标和地形匹配位置坐标对惯性导航系统的速度进行修正:7a)计算k时刻实测高程序列起始点位置坐标和k时刻地形匹配位置坐标的差值矩阵Z<sub>k</sub>:<maths num="0001"><math><![CDATA[<mrow><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mrow><mi>I</mi><mi>N</mi><mi>S</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>T</mi><mi>E</mi><mi>R</mi></mrow><mi>k</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mrow><mi>I</mi><mi>N</mi><mi>S</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mrow><mi>T</mi><mi>E</mi><mi>R</mi></mrow><mi>k</mi></msubsup></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001072678570000011.GIF" wi="310" he="127" /></maths>其中,k为大于0的自然数,<img file="FDA0001072678570000012.GIF" wi="65" he="55" />为k时刻实测高程序列起始点位置的东向坐标值,<img file="FDA0001072678570000013.GIF" wi="70" he="55" />为k时刻实测高程序列起始点位置的北向坐标值,<img file="FDA0001072678570000014.GIF" wi="70" he="63" />为k时刻地形匹配位置的东向坐标值,<img file="FDA0001072678570000015.GIF" wi="70" he="61" />为k时刻地形匹配位置的北向坐标值;7b)将差值矩阵Z<sub>k</sub>作为卡尔曼滤波器的测量值输入给卡尔曼滤波器,并输入卡尔曼滤波器均方误差矩阵初始值p<sub>0</sub>和卡尔曼滤波器状态估计值初始值<img file="FDA0001072678570000016.GIF" wi="73" he="54" />进行卡尔曼滤波递推计算,输出卡尔曼滤波速度状态估计值<img file="FDA0001072678570000017.GIF" wi="67" he="69" />7c)根据卡尔曼滤波器输出滤波速度状态估计值<img file="FDA0001072678570000018.GIF" wi="270" he="87" />对惯性导航系统的速度进行修正,得到修正后的速度<img file="FDA0001072678570000019.GIF" wi="51" he="39" /><maths num="0002"><math><![CDATA[<mrow><mover><mi>v</mi><mo>^</mo></mover><mo>=</mo><mi>v</mi><mo>+</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0001072678570000021.GIF" wi="198" he="70" /></maths>其中,v为前时刻惯性导航系统的速度,(·)<sup>T</sup>表示矩阵的转置。
地址 710071 陕西省西安市太白南路2号