发明名称 一种双星座组合导航系统几何精度因子计算方法
摘要 本发明公开了一种双星座组合导航系统几何精度因子计算方法;双星座组合导航系统几何精度因子计算公式中当几何矩阵H为方阵时,将几何精度因子计算公式进行转化,通过计算det(HH<sup>T</sup>W)和tr[adj(HH<sup>T</sup>W)]得到双星座组合导航系统几何精度因子;本发明将传统的计算公式进行转换,通过分别计算det(HH<sup>T</sup>W)和tr[adj(HH<sup>T</sup>W)],克服了传统计算方法存在的计算量大及计算过程发散等问题,实现了双星座组合导航系统几何精度因子的快速、准确计算。
申请公布号 CN105589078A 申请公布日期 2016.05.18
申请号 CN201510920489.9 申请日期 2015.12.11
申请人 电子科技大学 发明人 滕云龙;黄琦;井实;李坚;卢有亮
分类号 G01S19/01(2010.01)I 主分类号 G01S19/01(2010.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 周永宏;王伟
主权项 一种双星座组合导航系统几何精度因子计算方法,几何精度因子计算公式为:<maths num="0001"><math><![CDATA[<mrow><mi>G</mi><mi>D</mi><mi>O</mi><mi>P</mi><mo>=</mo><msqrt><mrow><mi>t</mi><mi>r</mi><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msup><mi>H</mi><mi>T</mi></msup><mi>W</mi><mi>H</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&rsqb;</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000874785470000011.GIF" wi="549" he="149" /></maths>其中,GDOP为几何精度因子,几何矩阵<maths num="0002"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>h</mi><mn>1</mn></msub></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd></mtr><mtr><mtd><msub><mi>h</mi><mi>&alpha;</mi></msub></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>h</mi><mrow><mi>&alpha;</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd></mtr><mtr><mtd><msub><mi>h</mi><mrow><mi>&alpha;</mi><mo>+</mo><mi>&beta;</mi></mrow></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000874785470000012.GIF" wi="451" he="454" /></maths>h<sub>i</sub>表示与接收机初始位置及卫星坐标相关的方向余弦向量,i=1,…,(α+β),α和β表示双星座组合导航系统卫星个数,tr(·)表示矩阵求迹,W为加权矩阵,H<sup>T</sup>为转置矩阵;其特征在于,当几何矩阵H为方阵时,根据矩阵理论得到tr[(H<sup>T</sup>WH)<sup>‑1</sup>]=tr[(HH<sup>T</sup>W)<sup>‑1</sup>],将几何精度因子计算公式转化为:<maths num="0003"><math><![CDATA[<mrow><mi>G</mi><mi>D</mi><mi>O</mi><mi>P</mi><mo>=</mo><msqrt><mrow><mi>t</mi><mi>r</mi><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msup><mi>HH</mi><mi>T</mi></msup><mi>W</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&rsqb;</mo></mrow></msqrt><mo>=</mo><msqrt><mfrac><mrow><mi>t</mi><mi>r</mi><mo>&lsqb;</mo><mi>a</mi><mi>d</mi><mi>j</mi><mrow><mo>(</mo><msup><mi>HH</mi><mi>T</mi></msup><mi>W</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><mi>det</mi><mrow><mo>(</mo><msup><mi>HH</mi><mi>T</mi></msup><mi>W</mi><mo>)</mo></mrow></mrow></mfrac></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000874785470000013.GIF" wi="1021" he="199" /></maths>其中,det(HH<sup>T</sup>W)通过分别计算矩阵行列式det(W)和det(H)得到,tr[adj(HH<sup>T</sup>W)]通过计算矩阵乘积HH<sup>T</sup>得到。
地址 611731 四川省成都市高新区(西区)西源大道2006号