发明名称 一种高精度确定载体动态加速度的方法
摘要 本发明公开了一种高精度确定载体动态加速度的方法,将卫星导航接收机固联安装在航空重力测量载体上,通过卫星导航接收机的观测信息直接解算出载体的动态加速度。本发明具有精度高、无需基站配合,作业简单、适用范围广等优点。
申请公布号 CN103529482B 申请公布日期 2016.05.11
申请号 CN201310511193.2 申请日期 2013.10.25
申请人 中国人民解放军国防科学技术大学 发明人 张开东;李显;吴美平
分类号 G01V7/00(2006.01)I;G01P15/00(2006.01)I 主分类号 G01V7/00(2006.01)I
代理机构 湖南兆弘专利事务所 43008 代理人 周长清
主权项 一种确定载体动态加速度的方法,其特征在于,将卫星导航接收机固联安装在航空重力测量载体上,通过卫星导航接收机的观测信息直接解算出载体的动态加速度;所述解算载体动态加速度的过程包括:(1.1)利用相位观测值的线性组合消除电离层散射的影响;(1.2)利用精密星历确定导航卫星的位置;(1.3)利用精密钟差产品改正卫星钟差的影响;(1.4)采用建模的方法消除对流层折射的影响;所述步骤(1.1)的具体流程为:(1.1.1)利用卫星导航接收机两个频点的伪距和载波相位观测值得到以下四个观测方程:<img file="FDA0000901533840000011.GIF" wi="1502" he="247" /><img file="FDA0000901533840000012.GIF" wi="1494" he="247" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mi>c</mi></mfrac><msubsup><mi>&rho;</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>dt</mi><mi>k</mi></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msup><mi>dt</mi><mi>p</mi></msup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mi>c</mi></mfrac><msubsup><mi>T</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msubsup><mi>I</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>c</mi><mo>&CenterDot;</mo><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup><mo>+</mo><msubsup><mi>e</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000901533840000013.GIF" wi="1510" he="151" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mi>c</mi></mfrac><msubsup><mi>&rho;</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>dt</mi><mi>k</mi></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msup><mi>dt</mi><mi>p</mi></msup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mi>c</mi></mfrac><msubsup><mi>T</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msubsup><mi>I</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>c</mi><mo>&CenterDot;</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mn>3</mn></mrow><mi>p</mi></msubsup><mo>+</mo><msubsup><mi>e</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000901533840000014.GIF" wi="1508" he="150" /></maths>其中,上标p表示导航卫星;下标k表示导航接收机;下标1,2分别表示导航信号频点,括号中的时间t指明观测时刻;<img file="FDA0000901533840000015.GIF" wi="45" he="60" />为载波相位观测值;R为伪距观测值;ρ为接收机和卫星之间的几何距离;f为导航信号载波频率;dt<sub>k</sub>为接收机r时钟误差;dt<sup>p</sup>为导航卫星k时钟误差;c为真空光速,且c=299792458m/s;T为接收机对卫星观测传播路径上的对流层延迟;I为接收机对卫星观测传播路径上的自由电子总含量;N为相位观测值的整周模糊度;ε为相位观测值的随机测量噪声;e为伪距观测值的随机测量噪声;b为不同通道之间的偏差;<img file="FDA0000901533840000016.GIF" wi="317" he="79" />为相位观测值初始的随机相位偏差,t<sub>0</sub>为信号开始锁定时刻;(1.1.2)通过对上述四个方程进行线性组合,消除一阶电离层传播误差的影响,得到的组合相位观测值称为“消电离层组合观测值”,组合方式为:<img file="FDA0000901533840000017.GIF" wi="1038" he="135" />和<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi><mi>o</mi><mi>n</mi><mi>o</mi><mo>-</mo><mi>f</mi><mi>r</mi><mi>e</mi><mi>e</mi></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000901533840000021.GIF" wi="1053" he="143" /></maths>其中,<img file="FDA0000901533840000022.GIF" wi="181" he="71" />为消电离层组合相位观测值;<img file="FDA0000901533840000023.GIF" wi="174" he="71" />为消电离层组合伪距观测值;将式(1)~(4)分别代入式(5)和(6),可得:<img file="FDA0000901533840000024.GIF" wi="1658" he="286" />和<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi><mi>o</mi><mi>n</mi><mi>o</mi><mo>-</mo><mi>f</mi><mi>r</mi><mi>e</mi><mi>e</mi></mrow><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>-</mo><msub><mi>f</mi><mn>2</mn></msub></mrow><mi>c</mi></mfrac><mo>&lsqb;</mo><msubsup><mi>&rho;</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>T</mi><mi>k</mi><mi>p</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mrow><mo>(</mo><msub><mi>f</mi><mn>1</mn></msub><mo>-</mo><msub><mi>f</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&lsqb;</mo><msub><mi>dt</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>dt</mi><mi>p</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mi>c</mi><mi>o</mi><mi>d</mi><mi>e</mi></mrow><mrow><mo>*</mo><mi>p</mi></mrow></msubsup><mo>+</mo><msubsup><mi>e</mi><mi>k</mi><mrow><mo>*</mo><mi>p</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000901533840000025.GIF" wi="1694" he="127" /></maths>其中:<img file="FDA0000901533840000026.GIF" wi="228" he="87" />分别表示吸收了初始相位偏差的模糊度;<img file="FDA0000901533840000027.GIF" wi="387" he="136" />为组合相位通道偏差;<img file="FDA0000901533840000028.GIF" wi="125" he="85" />为组合伪距通道偏差,且<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mi>c</mi><mi>o</mi><mi>d</mi><mi>e</mi></mrow><mrow><mo>*</mo><mi>p</mi></mrow></msubsup><mo>=</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup><mo>+</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mn>3</mn></mrow><mi>p</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000901533840000029.GIF" wi="647" he="140" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>&epsiv;</mi><mi>k</mi><mrow><mo>*</mo><mi>p</mi></mrow></msubsup><mo>=</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>&epsiv;</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>+</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>&epsiv;</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup></mrow>]]></math><img file="FDA00009015338400000210.GIF" wi="600" he="143" /></maths>为组合相位测量噪声;<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>k</mi><mrow><mo>*</mo><mi>p</mi></mrow></msubsup><mo>=</mo><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>e</mi><mrow><mi>k</mi><mo>,</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>+</mo><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>e</mi><mrow><mi>k</mi><mo>,</mo><mn>2</mn></mrow><mi>p</mi></msubsup></mrow>]]></math><img file="FDA00009015338400000211.GIF" wi="581" he="136" /></maths>为组合伪距测量噪声;通过方程(7)、(8)可知,观测值之间的线性组合消除了一阶电离层的影响。
地址 410073 湖南省长沙市砚瓦池正街47号中国人民解放军国防科学技术大学三院自动控制系