发明名称 一种正反向组合导航结果融合的POS后处理方法
摘要 本发明属于一种POS数据的后处理方法,具体涉及一种正反向组合导航结果融合的POS后处理方法。目的是实现正、反向组合导航结果的最优融合,获得全局最优的高精度位置姿态后处理结果。正、反向组合导航结果融合方程的表达形式如下:N=(P1+P2)-1(P2N1+P1N2),P1表示正向组合导航的误差估计协方差阵,N1表示正向组合导航的结果;P2表示反向组合导航的误差估计协方差阵,N2表示反向组合导航的结果,N表示POS的后处理结果。本发明的优点是采用导航结果融合方程对正、反向的组合导航结果进行融合处理,可以获得全局优于两个单向过程的位置姿态后处理结果。
申请公布号 CN102997908A 申请公布日期 2013.03.27
申请号 CN201110273066.4 申请日期 2011.09.15
申请人 北京自动化控制设备研究所 发明人 扈光锋;周祖洋;周东灵;李文耀;尚克军;张勤拓;王子静;王黎斌;谢仕民
分类号 G01C21/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 核工业专利中心 11007 代理人 李烨
主权项 1.一种正反向组合导航结果融合的POS后处理方法,其特征在于:设P<sup>1</sup>表示正向组合导航的误差估计协方差阵,N<sup>1</sup>表示正向组合导航的结果;P<sup>2</sup>表示反向组合导航的误差估计协方差阵,N<sup>2</sup>表示反向组合导航的结果,N表示POS的后处理结果;正、反向组合导航结果融合方程的表达形式如下:N=F(P<sup>1</sup>,P<sup>2</sup>)N<sup>1</sup>+G(P<sup>1</sup>,P<sup>2</sup>)N<sup>2</sup>                                (1)对于互不相关的正、反向两个组合导航过程,采用正、反向组合导航的误差估计协方差阵构造融合方程的权值系数如下:F(P<sup>1</sup>,P<sup>2</sup>)=(P<sup>1</sup>+P<sup>2</sup>)<sup>-1</sup>P<sup>2</sup>                                    (2)G(P<sup>1</sup>,P<sup>2</sup>)=(P<sup>1</sup>+P<sup>2</sup>)<sup>-1</sup>P<sup>1</sup>                                    (3)组合导航结果融合方程的表达形式为:N=(P<sup>1</sup>+P<sup>2</sup>)<sup>-1</sup>P<sup>2</sup>N<sup>1</sup>+(P<sup>1</sup>+P<sup>2</sup>)<sup>-1</sup>P<sup>1</sup>N<sup>2</sup>                            (4)将融合方程进一步写作:N=(P<sup>1</sup>+P<sup>2</sup>)<sup>-1</sup>(P<sup>2</sup>N<sup>1</sup>+P<sup>1</sup>N<sup>2</sup>)                                   (5)从而,根据上述组合导航结果融合方程的表达形式,采用正、反向组合导航的误差估计协方差阵的相应对角线元素,给出POS的位置融合方程为:<img file="FDA0000091443720000011.GIF" wi="1789" he="81" /><maths num="0001"><![CDATA[<math><mrow><mi>&lambda;</mi><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&lambda;</mi><mn>1</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>&lambda;</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&lambda;</mi><mn>2</mn></msubsup><msup><mi>&lambda;</mi><mn>1</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>&lambda;</mi><mn>1</mn></msubsup><msup><mi>&lambda;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>h</mi><mn>1</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>h</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>h</mi><mn>2</mn></msubsup><msup><mi>h</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><msubsup><mi>P</mi><mi>h</mi><mn>1</mn></msubsup><msup><mi>h</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:<img file="FDA0000091443720000014.GIF" wi="267" he="53" />分别为融合、正向组合、反向组合的纬度,单位:度;λ、λ<sup>1</sup>、λ<sup>2</sup>分别为融合、正向组合、反向组合的经度,单位:度;h、h<sup>1</sup>、h<sup>2</sup>分别为融合、正向组合、反向组合的高度,单位:米;<img file="FDA0000091443720000021.GIF" wi="167" he="62" />分别为正向组合、反向组合的纬度误差估计协方差;<img file="FDA0000091443720000022.GIF" wi="167" he="56" />分别为正向组合、反向组合的经度误差估计协方差;<img file="FDA0000091443720000023.GIF" wi="167" he="56" />分别为正向组合、反向组合的高度误差估计协方差。POS的姿态角融合方程为:<maths num="0003"><![CDATA[<math><mrow><mi>&gamma;</mi><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&gamma;</mi><mn>1</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>&gamma;</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&gamma;</mi><mn>2</mn></msubsup><msup><mi>&gamma;</mi><mn>1</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>&gamma;</mi><mn>1</mn></msubsup><msup><mi>&gamma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>&psi;</mi><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&psi;</mi><mn>1</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>&psi;</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&psi;</mi><mn>2</mn></msubsup><msup><mi>&psi;</mi><mn>1</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>&psi;</mi><mn>1</mn></msubsup><msup><mi>&psi;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&theta;</mi><mn>1</mn></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>&theta;</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>&theta;</mi><mn>2</mn></msubsup><msup><mi>&theta;</mi><mn>1</mn></msup><mo>+</mo><msubsup><mi>P</mi><mi>&theta;</mi><mn>1</mn></msubsup><msup><mi>&theta;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:γ、γ<sup>1</sup>、γ<sup>2</sup>分别为融合、正向组合、反向组合的滚动角,单位:度;ψ、ψ<sup>1</sup>、ψ<sup>2</sup>分别为融合、正向组合、反向组合的航向角,单位:度;θ、θ<sup>1</sup>、θ<sup>2</sup>分别为融合、正向组合、反向组合的俯仰角,单位:度;<img file="FDA0000091443720000027.GIF" wi="167" he="62" />分别表示正向组合、反向组合的滚动角误差估计协方差;<img file="FDA0000091443720000028.GIF" wi="172" he="62" />分别表示正向组合、反向组合的航向角误差估计协方差;<img file="FDA0000091443720000029.GIF" wi="167" he="56" />分别表示正向组合、反向组合的俯仰角误差估计协方差。
地址 100074 北京市丰台区云岗北区西里1号院