发明名称 基于GEO星座粗定位的北斗接收机快速定位方法
摘要 针对传统GNSS信号发射时间的整毫秒恢复算法,在接收机概略位置未知时存在计算量剧增而无法应用的问题,本发明利用北斗系统混合星座中GEO卫星电文速率高的这一特点,提供了一种基于GEO星座粗定位的北斗接收机快速定位方法。在接收机概略位置完全未知的情况下,首先使用已获取完整信号发射时刻的GEO卫星进行粗定位,然后根据粗定位获得的概率位置对非GEO卫星发射时刻毫秒整数时间进行恢复,最后利用所有可见卫星进行精确定位。本发明能够在北斗系统仅完成GEO星座帧同步的情况下,对完成捕获的非GEO卫星进行毫秒恢复,减小卫星星历已知条件下的接收机首次定位时间。
申请公布号 CN104865587A 申请公布日期 2015.08.26
申请号 CN201510328521.4 申请日期 2015.06.15
申请人 中国人民解放军国防科学技术大学 发明人 吴鹏;靖守让;刘文祥;苏映雪;孙广富;黄龙;刘增军;龚航;黄仰博
分类号 G01S19/25(2010.01)I 主分类号 G01S19/25(2010.01)I
代理机构 北京中济纬天专利代理有限公司 11429 代理人 胡伟华
主权项 一种基于GEO星座粗定位的北斗接收机快速定位方法,其特征在于包括以下步骤:步骤S1,判断GEO卫星是否完成帧同步,获取卫星发射时刻;步骤S2,卫星根据GEO卫星信号发射时刻和各类卫星传输时延对本地时间进行粗略校正,这里粗略校正使用单颗GEO卫星进行粗校正或者使用多颗卫星进行粗校正;步骤S3,使用已经获取了完整伪距的GEO星座进行粗定位,粗定位方法直接使用最小二乘、岭估计或者多普勒辅助的岭估计;步骤S4,在概略用户位置获取后,结合非GEO卫星的毫秒以下伪距对其进行完整伪距恢复;在接收完成帧同步和仅完成码相位测量的伪距观测量可以分别写为如下形式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mrow><msub><mi>&rho;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>r</mi><mi>i</mi></msub><mo>+</mo><msub><mi>b</mi><mi>u</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>N</mi><mi>j</mi></msub><mo>&CenterDot;</mo><mi>c</mi><mo>/</mo><mi>le</mi><mn>3</mn><mo>+</mo><msubsup><mi>&rho;</mi><mi>j</mi><mrow><mi>f</mi><mi>r</mi><mi>c</mi><mi>a</mi></mrow></msubsup><mo>=</mo><msub><mi>r</mi><mi>j</mi></msub><mo>+</mo><msub><mi>b</mi><mi>u</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000737932850000011.GIF" wi="493" he="164" /></maths>上式为卫星i的完整伪距和卫星j仅完成码相位测量的伪距表示;其中,ρ<sub>i</sub>为经过伪距改正的卫星i到用户的伪距,b<sub>u</sub>为用户钟差产生的等效距离误差;<img file="FDA0000737932850000012.GIF" wi="114" he="86" />为毫米以下的小数伪距值,N<sub>j</sub>表示毫米伪距模糊数,c表示光速,r<sub>i</sub>和r<sub>j</sub>分别表示用户到卫星i和卫星j的空间距离;选择第i颗卫星作为参考卫星进行伪距星间单差并通过取整操作可以得到卫星j的伪距毫秒数的计算公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>N</mi><mi>j</mi></msub><mo>=</mo><mi>r</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><mo>(</mo><msub><mi>&rho;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>r</mi><mi>i</mi></msub><mo>+</mo><msub><mi>r</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>&rho;</mi><mi>j</mi><mrow><mi>f</mi><mi>r</mi><mi>c</mi><mi>a</mi></mrow></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>le</mi><mn>3</mn><mo>/</mo><mi>c</mi><mo>)</mo></mrow>]]></math><img file="FDA0000737932850000013.GIF" wi="721" he="94" /></maths>步骤S5,非GEO卫星完整伪距恢复后,即可利用混合星座所有捕获到的卫星直接进行最小二乘定位;步骤S6,根据定位后的伪距残差判断伪距恢复是否合理;由于使用GEO星座进行粗定位获取用户概略位置,如果概略位置偏差较大会导致非GEO卫星的伪距恢复错误,因此需要对快速定位结果进行检验;伪距恢复后利用最小二乘迭代可以估计出用户位置同时可以计算出各卫星的伪距残差值,残差值和残差矢量的RMS值计算如下;式中<img file="FDA0000737932850000016.GIF" wi="39" he="65" />表示估计获得的用户位置与卫星i的空间距离;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>v</mi><mi>i</mi></msub><mo>=</mo><msub><mi>&rho;</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>r</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>b</mi><mi>u</mi></msub></mrow>]]></math><img file="FDA0000737932850000014.GIF" wi="286" he="69" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>R</mi><mi>M</mi><mi>S</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>v</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mo>/</mo><mi>n</mi></mrow></msqrt></mrow>]]></math><img file="FDA0000737932850000015.GIF" wi="452" he="158" /></maths>由于伪距恢复针对整毫秒恢复,如果伪距恢复出现错误时定位后会出现较大的伪距残差,利用定位后的伪距残差值的均方根是否超限判断快速定位方法是否成功,残差门限值根据测量精度或/和误差修正精度设置。
地址 410073 湖南省长沙市开福区德雅路109号