发明名称 一种区域定位导航增强信息提取算法
摘要 本发明涉及导航信号增强信息提取算法,具体为一种区域定位导航增强信息提取算法,包括以下步骤:多系统历史中频数据经过卫星信号接收机处理,由卫星信号接收机输出各路卫星信号的伪距和导航电文;通过高度角先验定权确定初始权矩阵;初始权矩阵通过基于M残差的方差分量估计验后定权法来确定最佳权矩阵;各路卫星信号的伪距和导航电文与最佳权矩阵进行加权融合,得到定位导航增强信息;定位导航增强信息与三维坐标进行对应,完成区域定位导航增强信息的提取。本发明有效的提高卫星系统定位的精度和实时性,特别是在遮挡严重、多径效应严重、可见卫星数目不足的恶劣环境中,在本算法提取出来的区域增强信息的辅助下,可以实现伪距精密单点定位。
申请公布号 CN105182379A 申请公布日期 2015.12.23
申请号 CN201510651325.0 申请日期 2015.10.10
申请人 太原理工大学 发明人 李灯熬;吴刚;赵菊敏
分类号 G01S19/37(2010.01)I 主分类号 G01S19/37(2010.01)I
代理机构 太原科卫专利事务所(普通合伙) 14100 代理人 朱源
主权项 一种区域定位导航增强信息提取算法,其特征在于包括以下步骤:四个卫星系统历史中频数据经过卫星信号接收机捕获、跟踪、位同步和帧同步处理,由卫星信号接收机输出各路卫星信号的伪距测量值和导航电文;四个卫星系统各自的权矩阵分别为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mfrac><mrow><msub><mi>sin&theta;</mi><mn>1</mn></msub></mrow><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mi>i</mi></msub><mn>2</mn></msubsup></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><msub><mi>sin&theta;</mi><mn>2</mn></msub></mrow><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mi>i</mi></msub><mn>2</mn></msubsup></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>...</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><msub><mi>sin&theta;</mi><mi>j</mi></msub></mrow><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mi>i</mi></msub><mn>2</mn></msubsup></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000818341690000011.GIF" wi="739" he="558" /></maths>i=1,2,3,4,分别代表GPS、BDS、GLOANSS、Galileo卫星系统,θ<sub>j</sub>为第j颗卫星的高度角;<img file="FDA0000818341690000012.GIF" wi="76" he="84" />为单点定位时观测卫星在天顶时的伪距观测单位权方差,则初始权矩阵P=(P<sub>1</sub>,P<sub>2</sub>,P<sub>3</sub>,P<sub>4</sub>);不同系统n颗卫星的伪距测量值观测方程组化矩阵△ρ=A△x+ε;式中△x为位置偏差,A表示系统矩阵,ε表示系统噪声;对n颗卫星的伪距测量值观测方程组化矩阵利用初始权矩阵P进行加权,得到P△ρ=PA△x+Pε,利用序贯最小二乘原理求解得到位置偏差△x=(A<sup>T</sup>PA)<sup>‑1</sup>A<sup>T</sup>P△ρ,并结合伪距测量值观测方程组化矩阵得到误差V=A△x‑△ρ;选择基于M残差的方差分量估计验后定权法来对误差进行迭代,使得误差V最小并得到最佳权矩阵P<sub>0</sub>:基于M残差的方差分量估计验后定权法的公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>=</mo><mi>W</mi><mo>,</mo></mrow>]]></math><img file="FDA0000818341690000013.GIF" wi="210" he="86" /></maths>式中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>&theta;</mi><mo>^</mo></mover><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mn>1</mn></msub><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mn>2</mn></msub><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mn>3</mn></msub><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>&sigma;</mi><msub><mn>0</mn><mn>4</mn></msub><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000818341690000014.GIF" wi="585" he="119" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>W</mi><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>V</mi><mn>1</mn><mi>T</mi></msubsup><msub><mi>P</mi><mn>1</mn></msub><msub><mi>V</mi><mn>1</mn></msub></mrow></mtd><mtd><mrow><msubsup><mi>V</mi><mn>2</mn><mi>T</mi></msubsup><msub><mi>P</mi><mn>2</mn></msub><msub><mi>V</mi><mn>2</mn></msub></mrow></mtd><mtd><mrow><msubsup><mi>V</mi><mn>3</mn><mi>T</mi></msubsup><msub><mi>P</mi><mn>3</mn></msub><msub><mi>V</mi><mn>3</mn></msub></mrow></mtd><mtd><mrow><msubsup><mi>V</mi><mn>4</mn><mi>T</mi></msubsup><msub><mi>P</mi><mn>4</mn></msub><msub><mi>V</mi><mn>4</mn></msub></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000818341690000015.GIF" wi="940" he="85" /></maths>S<sub>ii</sub>=k<sub>1</sub>tr(N<sup>(‑1)</sup>N<sub>i</sub>N<sup>(‑1)</sup>N<sub>i</sub>)‑2k<sub>2</sub>tr(N<sup>(‑1)</sup>N<sub>i</sub>)+n<sub>i</sub>,i=1,2,3,4,S<sub>ij</sub>=k<sub>1</sub>tr(N<sup>(‑1)</sup>N<sub>i</sub>N<sup>(‑1)</sup>N<sub>i</sub>),i=1,2,3,4,j=1,2,3,4,<img file="FDA0000818341690000021.GIF" wi="305" he="83" />其中A<sub>i</sub>表示第i个卫星系统单独定位时的系统矩阵,n<sub>i</sub>表示第i个卫星系统的可见卫星数目,N=N<sub>1</sub>+N<sub>2</sub>+N<sub>3</sub>+N<sub>4</sub>,其中:<img file="FDA0000818341690000022.GIF" wi="264" he="160" /><img file="FDA0000818341690000023.GIF" wi="269" he="160" />第i个系统的第K次更新后的权矩阵<img file="FDA0000818341690000024.GIF" wi="276" he="173" />式中c为任一常数,迭代计算至<img file="FDA0000818341690000025.GIF" wi="419" he="91" />时,误差V最小,此时得到的权矩阵P<sub>1</sub><sup>(k)</sup>就是最佳权矩阵P<sub>0</sub>;卫星信号的伪距测量值和导航电文通过最佳权矩阵P<sub>0</sub>进行加权融合,得到定位导航增强信息;定位导航增强信息与相应的三维位置坐标进行对应,得到区域定位导航增强信息,完成区域定位导航增强信息的提取。
地址 030024 山西省太原市迎泽西大街79号