发明名称 一种基于广域差分增强的卫星授时方法
摘要 本发明涉及一种基于广域差分增强的卫星授时方法,技术特征在于:各个差分站同时进行码相位和载波相位的伪距观测,且主差分站与卫星导航系统时间同步,均能测得伪距的码相位观测值和载波相位观测值位置。主要效果有两个:(1)提高卫星授时精度,利用本方法能使卫星授时达到1~3nsec。(2)简化地面监测系统的配置,降低差分网的建设成本。
申请公布号 CN101609143B 申请公布日期 2012.10.24
申请号 CN200910022301.3 申请日期 2009.04.30
申请人 中国科学院国家授时中心 发明人 蔡成林;李孝辉;吴海涛;王文利;常兴旺
分类号 G01S5/14(2006.01)I 主分类号 G01S5/14(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于广域差分增强的卫星授时方法,其特征在于:各个差分站同时进行码相位和载波相位的伪距观测,且主差分站与卫星导航系统时间同步,均能测得伪距的码相位观测值和载波相位观测值位置,具体步骤如下:步骤1:码相位的伪距观测量利用载波相位观测量平滑,得到平滑后的伪距观测量;步骤2:在主差分站进行星钟误差的第一次改正,得到星钟误差的第一次改正数<img file="FSB00000843222500011.GIF" wi="367" he="126" /><maths num="0001"><![CDATA[<math><mrow><msubsup><mi>&rho;</mi><mi>CLK</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&rho;</mi><mn>1</mn><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>R</mi><mn>1</mn><mi>k</mi></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中:<img file="FSB00000843222500013.GIF" wi="52" he="58" />为卫星k到主差分站的伪距,<img file="FSB00000843222500014.GIF" wi="138" he="66" />表示卫星k的解算位置到主差分站的几何距离,E表示利用卫星的广播星历解算获得的卫星位置,简称卫星的解算位置;步骤3:在副差分站进行星历误差的相对分离,得到该副差分站星历伪距的相对改正量<img file="FSB00000843222500015.GIF" wi="223" he="74" /><maths num="0002"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>&rho;</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><mi>&Delta;</mi><msubsup><mi>&rho;</mi><mi>i</mi><mn>1</mn></msubsup><mo>=</mo><msubsup><mi>&rho;</mi><mn>1</mn><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>R</mi><mn>1</mn><mi>k</mi></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mn>1</mn><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>R</mi><mn>1</mn><mn>1</mn></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>[</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>R</mi><mi>i</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>R</mi><mi>i</mi><mn>1</mn></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>其中:<img file="FSB00000843222500017.GIF" wi="91" he="61" />表示卫星k到副差分站i的站星星历伪距误差,<img file="FSB00000843222500018.GIF" wi="75" he="58" />表示基准星到副差分站i的站星星历伪距误差,<img file="FSB00000843222500019.GIF" wi="56" he="58" />表示卫星k到副差分站i的伪距,<img file="FSB000008432225000110.GIF" wi="55" he="63" />表示基准星到副差分站i的伪距,<img file="FSB000008432225000111.GIF" wi="47" he="69" />表示基准星到主差分站的伪距;<img file="FSB000008432225000112.GIF" wi="124" he="57" />表示基准星的解算位置到主差分站的几何距离,<img file="FSB000008432225000113.GIF" wi="130" he="58" />表示利用卫星k的解算位置到副差分站i的几何距离,<img file="FSB000008432225000114.GIF" wi="121" he="58" />表示基准星的解算位置到副差分站i的几何距离;步骤4:利用三个或三个以上的差分站得到的该副差分站星历伪距的相对改正量<img file="FSB000008432225000115.GIF" wi="233" he="74" />采用三维用户定位方程解算得到星历改正数;步骤5:在主差分站进行星钟误差的第二次改正,得到星钟第二次改正数<img file="FSB000008432225000116.GIF" wi="145" he="71" /><maths num="0003"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>&rho;</mi><mi>CLK</mi><mi>k</mi></msubsup><mo>=</mo><msubsup><mi>R</mi><mn>1</mn><mi>k</mi></msubsup><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>R</mi><mn>1</mn><mi>k</mi></msubsup><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><mi>&Delta;</mi><msub><mi>t</mi><mi>u</mi></msub><mo>-</mo><mi>&Delta;</mi><msub><mi>t</mi><msup><mi>u</mi><mo>&prime;</mo></msup></msub><mo>,</mo></mrow></math>]]></maths>其中<img file="FSB000008432225000118.GIF" wi="127" he="58" />表示卫星k的真实位置到主差分站的几何距离,Δt<sub>u</sub>表示用户定时结果的真值解,Δt<sub>u′</sub>表示用户定时结果的实际解;所述用户定时结果的真值解对于主差分站为Δt<sub>u</sub>=0;所述用户定时结果的实际解:将多颗卫星利用步骤2得到的星钟改正数和步骤4得到星历改正数,代入四维用户定位方程,得到主差分站的定时结果的实际解Δt<sub>u′</sub>;将主差分站的定时结果的实际解Δt<sub>u′</sub>作为星钟第二次改正数;步骤6:以第一次改正数<img file="FSB00000843222500021.GIF" wi="319" he="70" />和星钟第二次改正数<img file="FSB00000843222500022.GIF" wi="120" he="60" />修正卫星钟的星钟误差后进行授时。
地址 710072 陕西省西安市临潼区书院东路3号