发明名称 一种惯性信息辅助的卫星接收机基带信号处理方法
摘要 本发明公开了一种惯性信息辅助的卫星接收机基带信号处理方法,包括以下步骤:(1)建立适用于高动态4阶PLL载波跟踪环的状态空间方程;(2)使用4阶PLL载波跟踪环模型和预测估计方法估计载波相位及多普勒频移,用于计算载波NCO控制量<img file="DDA00002125508300011.GIF" wi="150" he="69" />(3)引入惯性辅助信息计算载体至卫星视线速度,得到当前多普勒频移<img file="DDA00002125508300012.GIF" wi="140" he="64" />用以辅助载波环路的跟踪,并计算载波NCO控制量<img file="DDA00002125508300013.GIF" wi="149" he="69" />(4)使用4阶PLL载波跟踪环估计的多普勒频移辅助码跟踪环,码跟踪环采用2阶DLL。本发明具有动态性高、跟踪误差小、鲁棒性好、适用范围广等优点。
申请公布号 CN102830408B 申请公布日期 2014.05.14
申请号 CN201210335344.9 申请日期 2012.09.11
申请人 中国人民解放军国防科学技术大学 发明人 唐康华;何晓峰;潘献飞;胡小平;郭瑶;罗兵;罗勇
分类号 H04L25/00(2006.01)I;G01S19/25(2010.01)I;G01S19/29(2010.01)I 主分类号 H04L25/00(2006.01)I
代理机构 湖南兆弘专利事务所 43008 代理人 周长清
主权项 1.一种惯性信息辅助的卫星接收机基带信号处理方法,其特征在于,包括以下步骤:(1)建立适用于高动态4阶PLL载波跟踪环的状态空间方程,其中以载波相位估计值<img file="FDA0000430740460000011.GIF" wi="98" he="82" />多普勒频移<img file="FDA0000430740460000012.GIF" wi="116" he="83" />多普勒频移率<img file="FDA0000430740460000013.GIF" wi="125" he="94" />多普勒频移加速度<img file="FDA0000430740460000014.GIF" wi="92" he="94" />作为状态变量,以载波相位鉴相器的输出<img file="FDA0000430740460000015.GIF" wi="68" he="78" />作为测量值;(2)使用4阶PLL载波跟踪环模型和预测估计方法估计载波相位及多普勒频移,用于计算载波NCO控制量<img file="FDA0000430740460000016.GIF" wi="169" he="90" />(3)引入惯性辅助信息计算载体至卫星视线速度,得到当前多普勒频移<img file="FDA0000430740460000017.GIF" wi="160" he="85" />用以辅助载波环路的跟踪,并计算载波NCO控制量<img file="FDA0000430740460000018.GIF" wi="167" he="90" />(4)使用4阶PLL载波跟踪环估计的多普勒频移辅助码跟踪环,码跟踪环采用2阶DLL;所述步骤(1)的4阶PLL状态空间方程的状态矢量为<img file="FDA0000430740460000019.GIF" wi="697" he="93" />其中,<img file="FDA00004307404600000110.GIF" wi="89" he="87" />为第k时刻的载波相位估计值,<img file="FDA00004307404600000111.GIF" wi="119" he="88" />为第k时刻的多普勒频移,<img file="FDA00004307404600000112.GIF" wi="116" he="99" />为第k时刻的多普勒频移率,<img file="FDA00004307404600000113.GIF" wi="118" he="100" />为第k时刻的多普勒频移加速度;所述步骤(3)中载波NCO控制量的计算公式为:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>f</mi><mo>^</mo></mover><mrow><mi>ca</mi><mo>,</mo><mi>NCO</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>f</mi><mi>IF</mi></msub><mo>+</mo><mi>&Delta;</mi><msubsup><mover><mi>f</mi><mo>~</mo></mover><mrow><mi>ca</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>+</mo><mfrac><mrow><msubsup><mover><mi>&phi;</mi><mo>~</mo></mover><mrow><mi>ca</mi><mo>,</mo><mi>k</mi></mrow><mo>+</mo></msubsup><mo>-</mo><msubsup><mover><mi>&phi;</mi><mo>~</mo></mover><mrow><mi>ca</mi><mo>,</mo><mi>k</mi></mrow><mo>-</mo></msubsup></mrow><mrow><mn>2</mn><mi>&pi;&tau;</mi></mrow></mfrac><mo>+</mo><mi>&Delta;</mi><msub><mover><mi>f</mi><mo>&OverBar;</mo></mover><mrow><mi>INS</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></math>]]></maths>其中,f<sub>IF</sub>为GNSS基带信号的中频频率,τ为环路更新时间,<img file="FDA00004307404600000115.GIF" wi="94" he="83" />为第k时刻的载波相位估计值,<img file="FDA00004307404600000116.GIF" wi="120" he="87" />为第k+1时刻的载波相位估计值,上标-和+分别表示更新前和更新后。
地址 410073 湖南省长沙市砚瓦池正街47号中国人民解放军国防科学技术大学三院