发明名称 一种导航卫星主备钟切换前后相位和频率变化测试方法
摘要 一种导航卫星主备钟切换前后相位和频率变化测试方法,用相位测量设备对卫星钟的输出信号进行相位测量,使用等间隔时间采集相位测量数据,长时间记录卫星钟输出切换开关动作前后的相位测量数据,并记录切换开关时刻,对切换前后的相位测试数据分别采用最小二乘法的曲线拟合,根据切换时刻对应的切换前后拟合值计算相位和频率的变化量。本方法可以降低卫星钟频率漂移和测量仪器对测试结果的影响,具有高测量精度。
申请公布号 CN103424755B 申请公布日期 2015.04.22
申请号 CN201310364546.0 申请日期 2013.08.20
申请人 北京空间飞行器总体设计部 发明人 崔小准;王璐;毕少筠;米红;李懿;刘崇华;刘天雄
分类号 G01S19/23(2010.01)I 主分类号 G01S19/23(2010.01)I
代理机构 中国航天科技专利中心 11009 代理人 陈鹏
主权项 一种导航卫星主备钟切换前后相位和频率变化测试方法,其特征在于步骤如下:(1)设定第一星载原子钟作为卫星钟主钟的频率参考,设定第二星载原子钟作为卫星钟备钟的频率参考,并设定初始状态卫星钟的输出信号由主钟提供,在进行时钟切换后卫星钟的输出信号改由卫星钟备钟提供;所述的第一星载原子钟和第二星载原子钟均置于真空环境中;(2)设定原子频标作为相位测试设备的频率参考,利用相位测试设备对卫星钟的输出信号进行相位测量,每秒一个测量值;首先,测量卫星钟主钟输出信号的相位,连续测量得到N个相位测量值<img file="FDA0000585807830000011.GIF" wi="431" he="77" />在测量得到<img file="FDA0000585807830000012.GIF" wi="76" he="77" />测量值时刻,进行卫星钟切换动作,继续测量卫星钟备钟输出信号的相位,连续测量得到N个相位测量值<img file="FDA0000585807830000013.GIF" wi="538" he="82" />(3)对切换前N个相位测量值进行二次曲线拟合,得到二次曲线前拟合系数,并利用得到的二次曲线前拟合系数计算得到切换时刻卫星钟相位的估计值<img file="FDA0000585807830000014.GIF" wi="95" he="74" />(4)对切换后N个相位测量值,预摘取卫星钟相位测量值<img file="FDA0000585807830000015.GIF" wi="422" he="80" /><img file="FDA0000585807830000016.GIF" wi="86" he="77" />进行二次曲线拟合,K=5~30,得到二次曲线后拟合系数,并利用得到的二次曲线后拟合系数计算得到切换时刻卫星钟相位的估计值<img file="FDA0000585807830000017.GIF" wi="96" he="75" />(5)将<img file="FDA0000585807830000018.GIF" wi="62" he="76" />与<img file="FDA0000585807830000019.GIF" wi="62" he="78" />求差后的绝对值作为卫星钟切换时的相位变化测量值Δφ;(6)对切换前的N个相位测量值<img file="FDA00005858078300000110.GIF" wi="404" he="78" />前后求差,得到N‑1个频差<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>y</mi><mn>1</mn><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>2</mn><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mn>1</mn><mo>*</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><msubsup><mi>y</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>n</mi><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>;</mo></mrow>]]></math><img file="FDA00005858078300000111.GIF" wi="760" he="80" /></maths>对切换后的N个相位测量值<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>N</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>N</mi><mo>+</mo><mn>2</mn></mrow><mo>*</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mrow>]]></math><img file="FDA00005858078300000112.GIF" wi="379" he="80" /></maths><img file="FDA00005858078300000113.GIF" wi="84" he="79" />前后求差,得到N‑1个频差<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>y</mi><mrow><mi>N</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mi>N</mi><mo>+</mo><mn>2</mn></mrow><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>N</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><msubsup><mi>y</mi><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mn>2</mn><mi>N</mi></mrow><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>;</mo></mrow>]]></math><img file="FDA00005858078300000114.GIF" wi="935" he="83" /></maths>(7)对卫星钟切换前的N-1个频差进行一次曲线拟合,得到一次曲线前拟合系数,并利用一次曲线前拟合系数计算得到切换时刻卫星钟频差的估计值<img file="FDA00005858078300000115.GIF" wi="134" he="71" />(8)对卫星钟切换后的N‑1个频差,预摘取卫星钟相位测量值<img file="FDA0000585807830000021.GIF" wi="147" he="81" /><img file="FDA0000585807830000022.GIF" wi="422" he="78" />进行一次曲线拟合,得到一次曲线后拟合系数,并利用一次曲线后拟合系数计算得到切换时刻卫星钟频差的估计值<img file="FDA0000585807830000023.GIF" wi="121" he="76" />(9)将<img file="FDA0000585807830000024.GIF" wi="100" he="76" />与<img file="FDA0000585807830000025.GIF" wi="94" he="76" />求差后的绝对值作为卫星钟切换时的频率变化测量值Δf。
地址 100094 北京市海淀区友谊路104号