发明名称 一种提高高程数据稳定性的滤波方法
摘要 本发明属于定位、导航数据的处理方法,是一种提高高程数据稳定性的滤波方法。目的是设计出一种性能优良的高程数据处理方法,提高高程数据的精度和分辨率,对高程投影误差的进行抑制。本发明包括如下步骤:步骤一、通过车辆的俯仰角以及前进速度计算车辆在高度方向上的速度;步骤二、通过GPS高程数据对高程计数据进行校正;步骤三、构建对高度数据进行最优估计的卡尔曼滤波模型;步骤四、通过步骤三中构建的滤波模型,进行卡尔曼滤波;步骤五、最后将高程数据进行输出。本发明的优点是采用了卡尔曼滤波的数据处理方法,利用多个数据源对高程数据进行修正,较强地抑制了车载定位系统中的高程投影误差,使车载定位精度提高了18%。
申请公布号 CN101750050B 申请公布日期 2011.08.17
申请号 CN200810180520.X 申请日期 2008.11.28
申请人 北京航天发射技术研究所 发明人 辛廷慧;陈令刚;佟艳艳;白焕旭;何定洲
分类号 G01C5/00(2006.01)I 主分类号 G01C5/00(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 1.一种提高高程数据稳定性的滤波方法,包括如下步骤:步骤一、通过车辆的俯仰角以及行驶速度计算车辆在高度方向上的速度v<sub>z</sub>=vsinα≈vα,v车辆行驶速度;α为车辆俯仰角;步骤二、通过GPS高程数据对高程计数据进行校正,校正通过以下公式实现:h=h<sub>q</sub>+K<sub>h</sub>(hG-h),h为校正后的高程数据;h<sub>q</sub>为高程计输出;h<sub>G</sub>为GPS高程数据;K<sub>h</sub>为校正增益;步骤三、通过步骤一、步骤二中得到的h与v<sub>z</sub>建立状态变量,构建对高度数据进行最优估计的卡尔曼滤波模型:<maths num="0001"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>h</mi></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>&Phi;</mi><mo>=</mo><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></mrow></math>]]></maths>X<sub>k+1</sub>=ΦX<sub>k</sub>+W<sub>k+1</sub>,方差E[W<sub>k</sub>W<sub>j</sub>]=R·δ<sub>k,j</sub>;Z<sub>k</sub>=X<sub>k</sub>+V<sub>k</sub>,方差E[V<sub>k</sub>V<sub>j</sub>]=Q·δ<sub>k,j</sub>;式中:k表示当前时刻,X为状态变量,X<sub>k</sub>为第k时刻的状态;T为滤波周期;Φ为状态转移矩阵;Z<sub>k</sub>为第k时刻的量测量;δ为狄拉克函数;W,V为模型误差与量测误差;R,Q为模型误差与量测误差的方差;步骤四、通过步骤三中构建的滤波模型,进行卡尔曼滤波,滤波所依据的数学方程为:<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msub><mo>=</mo><mi>&Phi;</mi><msub><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi></msub></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msub><mo>+</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>Z</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>K<sub>k+1</sub>=P<sub>k+1</sub>/k(P<sub>k+1/k</sub>+R)<sup>-1</sup>P<sub>k+1/k</sub>=ΦP<sub>k</sub>Φ<sup>T</sup>+QP<sub>k+1</sub>=(I-K<sub>k+1</sub>)P<sub>k+1/k</sub>式中:P为方差矩阵;K为滤波增益;<img file="F200810180520XC00021.GIF" wi="50" he="61" />为X<sub>k</sub>的最优估计值;步骤五、最后将状态变量<img file="F200810180520XC00022.GIF" wi="47" he="60" />中的高程数据进行输出:<maths num="0005"><![CDATA[<math><mrow><msub><mi>h</mi><mi>c</mi></msub><mo>=</mo><msub><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi></msub><mrow><mo>(</mo><mn>1,1</mn><mo>)</mo></mrow></mrow></math>]]></maths>即为最后得到的k时刻的高程数据。
地址 100076 北京市丰台区南大红门路1号
您可能感兴趣的专利