发明名称 一种基于双星时差频差无源定位的伪距差分方法
摘要 本发明提供了一种基于双星时差频差无源定位的伪距差分方法,依次测量两颗卫星到参考站的伪距差,计算卫星的伪距改正数,利用伪距改正数改正未知辐射源的伪距差,得到接近真值的到达时间差,测量参考站到两颗卫星的到达频率差为,计算参考站到两颗卫星的接近真值的到达频率差,利用卫星的频率改正数改正未知辐射源到两颗卫星的到达频率差,得到接近真值的到达频率差,利用地球椭球方程、接近真值的到达时间差和到达频率差来实现对未知辐射源的高精度定位。本发明可以对系统的定位误差进行改正,从而提高了双星无源定位系统的精度。
申请公布号 CN103645485B 申请公布日期 2016.01.20
申请号 CN201310519847.6 申请日期 2013.10.28
申请人 中国科学院国家授时中心 发明人 薛艳荣;李孝辉;许龙霞;任烨;李雨薇
分类号 G01S19/41(2010.01)I 主分类号 G01S19/41(2010.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 一种基于双星时差频差无源定位的伪距差分方法,其特征在于包括下述步骤:1)测量两颗卫星到参考站的伪距差<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>&rho;</mi><mn>0</mn><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mn>0</mn><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>-</mo><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>&Delta;&rho;</mi><mn>0</mn><mi>s</mi></msubsup><mo>+</mo><msub><mi>&Delta;&rho;</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000011.GIF" wi="1216" he="77" /></maths>式中<img file="FDA0000819814970000012.GIF" wi="178" he="72" />分别为两颗卫星到参考站的伪距,<img file="FDA0000819814970000013.GIF" wi="171" he="72" />分别为两颗卫星到参考站的几何距离,可以根据已知参考站的位置与卫星广播星历计算获得,即:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000014.GIF" wi="1437" he="91" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000015.GIF" wi="1454" he="98" /></maths>式中(x<sub>0</sub>,y<sub>0</sub>,z<sub>0</sub>)为参考站的坐标,(x<sub>si</sub>,y<sub>si</sub>,z<sub>si</sub>),i=1,2是两颗卫星的位置,<img file="FDA00008198149700000110.GIF" wi="81" he="60" />表示与地面参考站位置和卫星位置有关的测距误差,Δρ<sub>0</sub>表示与接收机有关的测距误差;2)计算卫星的伪距改正数:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&rho;</mi><mo>=</mo><mrow><mo>(</mo><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>r</mi><mn>0</mn><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mn>0</mn><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mn>0</mn><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000016.GIF" wi="1179" he="84" /></maths>3)利用伪距改正数Δρ改正未知辐射源的伪距差,得到接近真值的到达时间差Δt':<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>c&Delta;t</mi><mo>&prime;</mo></msup><mo>=</mo><msup><mi>&rho;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msup><mo>-</mo><msup><mi>&rho;</mi><mrow><mi>s</mi><mn>1</mn></mrow></msup><mo>+</mo><mi>&Delta;</mi><mi>&rho;</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>z</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>z</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000017.GIF" wi="1779" he="99" /></maths>式中(x,y,z)表示未知辐射源的位置,c为光速;4)测量参考站到两颗卫星的到达频率差为Δf<sub>r</sub>;5)计算参考站到两颗卫星的接近真值的到达频率差Δf<sub>0</sub>,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfrac><mrow><msub><mi>c&Delta;f</mi><mn>0</mn></msub></mrow><msub><mi>f</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>x</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>y</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>z</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>z</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000018.GIF" wi="1857" he="149" /></maths>式中(v<sub>xi</sub>,v<sub>yi</sub>,v<sub>zi</sub>),i=1,2是卫星的速度;其中f<sub>0</sub>为参考站的载波频率,r<sub>1</sub>,r<sub>2</sub>分别为卫星1和卫星2到未知辐射源的几何距离;6)计算卫星的频率改正数:ΔF=Δf<sub>0</sub>‑Δf<sub>Y</sub>           (7)7)未知辐射源到两颗卫星的到达频率差为Δf,利用频率改正数ΔF改正未知辐射源的到达频率差,得到接近真值的到达频率差Δf':<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfrac><mrow><msup><mi>c&Delta;f</mi><mo>&prime;</mo></msup></mrow><msub><mi>f</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mrow><mi>c</mi><mrow><mo>(</mo><mi>&Delta;</mi><mi>f</mi><mo>+</mo><mi>&Delta;</mi><mi>F</mi><mo>)</mo></mrow></mrow><msub><mi>f</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>x</mi></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>x</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>y</mi></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>y</mi><mn>2</mn></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>z</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>-</mo><mi>z</mi></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>x</mi></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>y</mi></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>z</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>z</mi></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac><msub><mi>v</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000019.GIF" wi="1759" he="138" /></maths>8)利用地球椭球方程:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfrac><msup><mi>x</mi><mn>2</mn></msup><msup><mi>a</mi><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><mi>y</mi><mn>2</mn></msup><msup><mi>a</mi><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><mi>z</mi><mn>2</mn></msup><msup><mi>b</mi><mn>2</mn></msup></mfrac><mo>=</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819814970000021.GIF" wi="1077" he="181" /></maths>式中a是地球半长轴,b是地球半短轴;9)由公式(5)、(8)和(9)解算未知辐射源的位置,进而实现高精度定位。
地址 710600 陕西省西安市临潼区书院东路3号