CN103792519B - 一种基于有源定标器的星载雷达时钟漂移在轨校正方法 - Google Patents

一种基于有源定标器的星载雷达时钟漂移在轨校正方法 Download PDF

Info

Publication number
CN103792519B
CN103792519B CN201410013535.2A CN201410013535A CN103792519B CN 103792519 B CN103792519 B CN 103792519B CN 201410013535 A CN201410013535 A CN 201410013535A CN 103792519 B CN103792519 B CN 103792519B
Authority
CN
China
Prior art keywords
spaceborne radar
source scaler
distance
spaceborne
radar
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410013535.2A
Other languages
English (en)
Other versions
CN103792519A (zh
Inventor
王彩云
郭伟
万珺之
赵飞
何佳宁
谌华
蔡朋飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN201410013535.2A priority Critical patent/CN103792519B/zh
Publication of CN103792519A publication Critical patent/CN103792519A/zh
Application granted granted Critical
Publication of CN103792519B publication Critical patent/CN103792519B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G—PHYSICS
    • G04—HOROLOGY
    • G04R—RADIO-CONTROLLED TIME-PIECES
    • G04R20/00—Setting the time according to the time information carried or implied by the radio signal
    • G04R20/20—Setting the time according to the time information carried or implied by the radio signal the radio signal being an AM/FM standard signal, e.g. RDS
    • G—PHYSICS
    • G04—HOROLOGY
    • G04R—RADIO-CONTROLLED TIME-PIECES
    • G04R40/00—Correcting the clock frequency

Abstract

本发明提出一种基于有源定标器的星载雷达时钟漂移在轨校正方法，通过采用有源定标器对星载雷达脉冲信号进行跟踪、接收并转发，通过对比有源定标器距离抛物线和星载雷达距离抛物线，得出星载雷达的时钟漂移量，基于此来实现星载雷达时钟运行状态的实时定量监测。本发明首次采用一种独立的方法精确测量了星载雷达的时钟漂移量，所述方法在星载雷达波束覆盖区域内任意位置，星载雷达过境期间任意时刻均可进行测试，具有很强的时空灵活性，通过在地面使用有源定标器接收脉冲时间间隔作为分析数据，对有源定标器无绝对定时基准要求，同时消除了星载雷达信号在传输过程中由大气造成的路径误差，进而能够实时准确的测量出星载雷达的时钟漂移量。

Description

一种基于有源定标器的星载雷达时钟漂移在轨校正方法
技术领域
本发明涉及星载雷达遥感技术领域，具体涉及一种基于有源定标器的星载雷达时钟漂移在轨校正方法。
背景技术
星载雷达在轨运行中，向地面发射脉冲信号，并接收地物目标的反射回波，通过测量脉冲的往返时间来确定星载雷达到被测目标的距离。星载雷达高度计就是这样一种指向星下点的主动式星载雷达。在轨运行期间，星载雷达垂直向下发射脉冲信号，信号经地球表面（海、陆、冰）或其他地面目标反射后回到星载雷达接收天线。星载雷达通过测量脉冲往返时间可以确定其质心到星下点的距离，进而计算星下点的海平面高度。通过分析返回脉冲的波形和强度，可以获取海面有效波高和海面风速等信息。因此，海平面高度、海面有效波高和海面风速是星载雷达高度计测量的三个基本参数。
其中，海平面高度的测量需要精确获得星载雷达发射/接收脉冲的往返时间，由d=tc/2，c为真空中的光速，求得星载雷达至地面（或其他反射目标）的距离d，然后根据地球物理参数修正星载雷达信号传输过程中的各项误差，以获取星载雷达星下点的海平面高度。
根据上述，星载雷达发射/接收脉冲往返时间t的测量精度直接影响到星下点海平面高度的测量精度，因此，为获得准确的海平面高度值，要求星载雷达具有高精度高稳定性的时钟系统。任何一种时钟系统，恒温晶振或原子钟，由于器件自身老化等原因，其频率都会发生随时间缓慢漂移的现象，时钟频率的变化将直接影响到星载雷达对地测量结果的准确性和数据的稳定性。
星载雷达发射前，可通过采用频率计或频谱仪测量的方法直接测量其时钟频率；星载雷达在轨运行期间，无法对其进行直接测量。目前，在国内外尚未发现一种专门测量星载雷达在轨运行期间时钟漂移量的方法，国际上通常将时钟漂移造成的误差和星载雷达硬件系统固有偏差做综合处理，作为星载雷达在轨测量的总误差。能够独立测量出星载雷达时钟漂移量可以准确分离出总误差中的可变误差量和固有误差量，时钟漂移作为可变误差量，其准确测量对于有效分析星载雷达整机设备的在轨运行状态，实时监测时钟部件的定时准确度具有非常重要的意义，进而对于精确测量星下点的海平面高度也具有重要的意义。
发明内容
本发明的目的在于提供一种基于有源定标器的星载雷达时钟漂移校正方法。采用有源定标器对星载雷达脉冲信号进行跟踪，接收并转发，通过对比有源定标器距离抛物线和星载雷达距离抛物线，得出星载雷达的时钟漂移量，基于此来实现星载雷达时钟运行状态的实时定量监测。
本发明解决上述技术问题所采取的技术方案如下：
一种基于有源定标器的星载雷达时钟漂移在轨校正方法，包括以下步骤：
步骤（1）：通过有源定标器跟踪并接收星载雷达的脉冲信号，并记录每个接收脉冲的到达时刻，得到有源定标器接收脉冲时刻序列；
步骤（2）：通过有源定标器内部的数据处理模块对有源定标器接收脉冲时刻序列的相邻两点求差值，得到有源定标器接收脉冲时间间隔序列Rn+1-Rn；
步骤（3）：通过有源定标器从星载雷达存储记录的数据中读取星载雷达每一个发射脉冲的时刻，得到星载雷达发射脉冲时刻序列，然后将该发射脉冲时刻序列的相邻两点求差值，得到星载雷达发射脉冲时间间隔序列tn+1-tn；
步骤（4）：通过有源定标器内部的数据处理模块将步骤（2）和步骤（3）得到的时间间隔序列变换为星载雷达与有源定标器间的距离之差序列；
步骤（5）：对步骤（4）得到的距离之差序列进行数据拟合得到以时间t为变量的一次直线函数；
步骤（6）：对步骤（5）得到的一次直线函数进行积分得到有源定标器距离抛物线，星载雷达时钟漂移误差使所述有源定标器距离抛物线的对称轴位置较理想状态下距离抛物线的对称轴位置发生偏移；
步骤（7）：从星载雷达存储记录的数据中读取星载雷达每一个发射脉冲的时刻和对应的每一个接收脉冲的时刻，并用每一个接收脉冲时刻减去对应的发射脉冲时刻，得到星载雷达收发脉冲时间间隔序列rn-tn；
步骤（8）：通过有源定标器内部的数据处理单元将步骤（7）得到的星载雷达收发脉冲时间间隔序列变换为星载雷达与有源定标器的距离序列，并通过对该距离序列进行数据拟合得到星载雷达距离抛物线函数，星载雷达时钟漂移误差引起所述星载雷达距离抛物线较理想状态下的距离抛物线发生上下平移；
步骤（9）：在有源定标器内部的数据处理模块中，对比步骤（6）所得到的有源定标器距离抛物线和步骤（8）所得到的星载雷达距离抛物线的对称轴位置，根据其对称轴位置之差得出有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差；
步骤（10）：在地面使用频率计或频谱仪，直接测量所述有源定标器的时钟频率，并根据步骤（9）得到的有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差，求出星载雷达的时钟频率测量值；
步骤（11）：将步骤（10）得到的星载雷达时钟频率测量值与星载雷达时钟频率理论设计值相比较，得出星载雷达时钟漂移量，基于所述时钟漂移量校正星载雷达时钟系统。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中步骤（1）中，有源定标器通过其内部的高速数据采集模块和精密的GPS定时模块记录每个接收脉冲的到达时刻。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中步骤（4）具体包括：通过有源定标器内部的数据处理模块将步骤2得到的有源定标器接收脉冲时间间隔序列Rn+1-Rn减去步骤3得到的对应的星载雷达发射脉冲时间间隔序列tn+1-tn，得到一个差值序列（Rn+1-Rn）－（tn+1-tn），再将该差值序列乘以光速c，得到相邻两个发射脉冲所对应的星载雷达与有源定标器间的距离之差序列。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中步骤（5）中所述拟合得到的一次直线函数为：，其中，b表示由时钟漂移误差造成的直线截距，a中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中步骤（6）中所述的有源定标器距离抛物线函数为：，其中ρ为星载雷达至有源定标器的距离，，b为时钟漂移引入的误差，a中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中步骤（8）具体包括：通过有源定标器内部的数据处理单元将步骤（7）得到的星载雷达收发脉冲时间间隔序列乘以光速c再除以2，以得到星载雷达至有源定标器的距离序列，然后对该距离序列进行数据拟合，得到以时间t为变量的二次函数曲线，然后将该二次函数曲线与理想情况下星载雷达与有源定标器间的距离抛物线函数曲线进行拟合比对得到实测情况下的星载雷达距离抛物线函数。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中所述实测情况下的星载雷达距离抛物线函数为：，其中，为星载雷达至有源定标器的距离，，d表示由星载雷达时钟漂移误差造成的抛物线函数的上下平移量，a、C中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径，H为有源定标器相对星下点的高度。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中各步骤中的数据存储于有源定标器数据存储单元中。
进一步的根据本发明所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其中理想状态下星载雷达与有源定标器间的距离关系为：，式中ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；ρ为星载雷达至有源定标器的距离；H为有源定标器相对星下点的高度；Re为地球半径；t为星载雷达飞行时间。
通过本发明的技术方案至少具备以下技术优势和创新效果：
1）、本发明首次采用一种独立的方法精确测量了星载雷达的时钟漂移量，是一种全新的星载雷达时钟漂移测量技术；
2）、本发明所述方法在星载雷达波束覆盖区域内任意位置，星载雷达过境期间任意时刻均可进行测试，具有很强的时空灵活性；
3）、本发明通过在地面使用有源定标器接收脉冲时间间隔作为分析数据，对有源定标器无绝对定时基准要求；
4）、本发明通过对有源定标器相邻两个接收脉冲的达到时刻做差值求得有源定标器接收脉冲时间间隔，消除了星载雷达信号在传输过程中由大气造成的路径误差，有源定标器接收脉冲时间间隔仅取决于星载雷达发射脉冲时间间隔和星载雷达至有源定标器的距离变化。
附图说明
附图1为本发明所述星载雷达与有源定标器之间的空间几何关系图；
附图2为本发明所述星载雷达在轨运行时与有源定标器之间的空间几何和信号收发时序关系图。
具体实施方式
以下结合附图对本发明的技术方案进行详细的描述，以使本领域技术人员能够更加清楚的理解本发明的方案，但并不因此限制本发明的保护范围。
首先给出本发明所述基于有源定标器的星载雷达时钟漂移在轨校正方法的基本原理过程。
星载雷达时钟漂移量的在轨定量化监测需使用布设在地面上的有源定标器。有源定标器对星载雷达发射的脉冲信号进行跟踪，记录，并发射射频信号至星载雷达。通过分析有源定标器的接收数据和星载雷达的接收数据，能够准确计算出星载雷达时钟漂移量。具体的如附图1所示，有源定标器设置于地面，星载雷达处于卫星轨道，有源定标器接收星载雷达的对地发射信号，并向星载雷达发射射频信号，星载雷达在轨飞行中，与布设在地面上的有源定标器之间的几何关系满足附图1所示关系，根据附图1所示的几何关系，利用余弦定理，可得出星载雷达至有源定标器的距离ρ与星载雷达飞行时间t的关系表达式（1）：
图中ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；ρ为星载雷达至有源定标器的距离；H为有源定标器相对星下点的高度；d为有源定标器至星载雷达天底点的地表距离；Re为地球半径；θ为距离d对应的地心张角；t为星载雷达的飞行时间，在星载雷达至有源定标器的距离最近的时刻为时间t的零时刻。
。
根据表达式（1）可知，星载雷达至有源定标器的距离ρ(t)是一个以时间t为变量的抛物线函数。
星载雷达在轨运行期间一直与有源定标器收发脉冲信号，具体的如附图2所示，图中tn为星载雷达第n个发射脉冲的发射时刻；rn为星载雷达第n个接收脉冲的接收时刻；对应星载雷达第n个发射脉冲时星载雷达至有源定标器的距离；Rn为有源定标器接收星载雷达第n个发射脉冲的接收时刻；n为自然数。
根据实际检测过程，以下参数条件能够准确获知，从而作为本发明所述方法的初始基本条件：
（1）星载雷达收发脉冲时间间隔可准确获知，即若星载雷达的第n个脉冲的发射/接收时刻为(tn,rn)，则(rn-tn)能够准确获知；
（2）星载雷达发射脉冲时间间隔可准确获知，即若星载雷达的第n-1个脉冲发射时刻为tn-1，第n个脉冲发射时刻为tn,则（tn-tn-1）能够准确获知；
（3）有源定标器能够准确跟踪并记录下每一个接收脉冲的相对到来时刻，即若第n个雷达脉冲达到有源定标器的时刻Rn，第n+1个脉冲到达有源定标器的时刻Rn+1，则要求有源定标器能够准确测得(Rn+1-Rn)的时间间隔。
在上述基本条件下给出本发明所述校正方法的基本原理：
有源定标器跟踪星载雷达的发射脉冲，并记录接收到的每一个脉冲的到达时刻，将每相邻两个脉冲的到达时刻做差值，得到一个序列，称为有源定标器接收脉冲时间间隔序列，该序列的拟合函数为一条直线，将该序列值逐点累加并转换为距离，即可得到在有源定标器上所获知的有源定标器与星载雷达间的距离关系，记为有源定标器距离抛物线，有源定标器距离抛物线中包含有星载雷达时钟频率和有源定标器时钟频率之间的频率偏差。
同时，星载雷达跟踪有源定标器发射的射频脉冲，并记录接收到的每一个脉冲的到达时刻，根据星载雷达自身记录下的发射脉冲时刻，可得到每一对发射/接收脉冲的时间间隔序列，将该序列转换为距离并做拟合，即可得到在星载雷达上所获知的有源定标器与星载雷达间的距离关系，记为星载雷达距离抛物线，星载雷达距离抛物线中也包含有星载雷达时钟频率和有源定标器时钟频率之间的时钟频率偏差。
但是在有源定标器距离抛物线中所包含的星载雷达时钟频率和有源定标器时钟频率之间的频率偏差反应在其抛物线的对称轴上（以下详述），而在星载雷达距离抛物线中所包含的星载雷达时钟频率和有源定标器时钟频率之间的频率偏差则反应在其抛物线的上下平移量上（以下详述），这种频率偏差并不引起星载雷达距离抛物线形状和对称轴位置的变化，所以通过对比星载雷达距离抛物线和有源定标器距离抛物线的对称轴位置，可以精确测量星载雷达时钟频率和有源定标器时钟频率之间的频率偏差，在通过地面上使用频率计或频谱仪直接测量有源定标器的时钟频率，即可准确求得星载雷达的时钟频率，进一步与其理论标准值对比即可求出星载雷达的时钟漂移量。下面从各距离抛物线的确定入手给出具体过程。
获得有源定标器距离抛物线关系。
根据星载雷达与有源定标器之间的空间几何和时序关系，对于有源定标器，其所获知的有源定标器与星载雷达间的距离关系，有如下关系式：
由于实际中（Rn-tn）不可获知，因此需将上面方程组内相邻两等式两两相减，得到：
式中c为光速，并令：
，表示相邻两个发射脉冲所对应的星载雷达与有源定标器间的距离之差；
，表示有源定标器接收到的相邻两个星载雷达发射脉冲的时刻之差，记为：有源定标器接收脉冲时间间隔；
，表示星载雷达相邻两个发射脉冲的发射时刻之差，记为：星载雷达发射脉冲时间间隔；
则上述(4)方程式变换为：
为有源定标器接收脉冲时间间隔，通过有源定标器内部的数据处理单元测得。为星载雷达发射脉冲时间间隔，其真实值无法直接测量，在计算中使用星载雷达发射脉冲时间间隔的理论设计值。由于星载雷达时钟存在时钟频率漂移，所以星载雷达发射脉冲时间间隔真实值和星载雷达发射脉冲时间间隔理论设计值不同，二者之间存在误差量。
根据上述表达式(1)，在没有误差的情况下，星载雷达至有源定标器的距离满足：，是一个以时间t为变量的抛物线函数，而是相邻两个发射脉冲所对应的星载雷达至有源定标器的距离之差，故序列等效于对表达式(1)所示抛物线函数的导函数，理论上为一次函数，因此在没有误差的情况下，，当t=0时，。
实际中由星载雷达时钟漂移等因素造成上述误差量的存在，有源定标器计算得到的表达式（5）所示的序列存在恒定误差，根据基于测量可确知的有源定标器接收脉冲时间间隔和星载雷达发射脉冲时间间隔，对表达式（5）进行数据拟合，得到一条直线函数，该函数的直线斜率与没有误差情况下的直线斜率基本相同，因此通过数据拟合得到存在误差量的情况下，表达式（5）所示的有源定标器测量得到的序列满足：
其中，a表示序列的直线斜率亦即，b表示由误差量造成的直线截距。
对（6）式所示有源定标器测量拟合得到的序列进行逐点累加，等效于对求积分，得到基于实际测量拟合的有源定标器距离抛物线，表达式为：
其中，，误差量引起的偏差体现在公式（7）所示有源定标器距离抛物线的对称轴上。
获得星载雷达距离抛物线关系。
根据星载雷达与有源定标器之间的空间几何和时序关系，对于星载雷达，理想情况下其所获知的有源定标器与星载雷达间的距离关系，有如下关系式：
实际中由于星载雷达存在时钟漂移，造成星载雷达测量到的脉冲收发间隔（rn-tn）序列存在误差量，所以星载雷达实际得到的距离方程组为：
由误差量造成的星载雷达至有源定标器的距离误差用d表示。根据基于测量可确知的星载雷达收发脉冲时间间隔序列，对进行数据拟合，得到一抛物线函数，此为考虑误差的星载雷达距离抛物线，所拟合得到的抛物线与表达式（1）所示理想情况下的有源定标器和星载雷达间的距离抛物线的形状和对称轴位置基本相同，因此通过数据拟合得到的存在误差量的情况下，星载雷达距离抛物线的表达式可表示为：
其中，，d表示由误差量造成的星载雷达距离抛物线的上下平移量。
由此可见星载雷达时钟漂移量不改变星载雷达距离抛物线的形状和对称轴位置，仅影响其上下平移量。所以通过对比星载雷达距离抛物线和有源定标器距离抛物线的对称轴位置，可以获知的数值，而其中的b则蕴含了由星载雷达时钟频率和有源定标器时钟频率间频率偏差造成的误差量，因此通过这种对比即可精确测量星载雷达时钟频率和有源定标器时钟频率之间的频率偏差，然后在通过地面上使用频率计或频谱仪直接测量有源定标器的时钟频率，即可准确求得星载雷达的时钟频率，进一步与其理论标准值对比即可求出星载雷达的时钟漂移量，基于该漂移量即可对星载雷达时钟系统进行校正。
接着在上述基本原理的基础上，给出本发明所述基于有源定标器的星载雷达时钟漂移在轨校正方法的具体步骤，包括：
步骤1：通过有源定标器跟踪并接收星载雷达的脉冲信号，通过其内部的高速数据采集模块和精密的GPS定时模块记录每个接收脉冲的达到时刻，得到有源定标器接收脉冲时刻序列；
步骤2：通过有源定标器内部的数据处理模块对有源定标器接收脉冲时刻序列的相邻两点求差值，得到有源定标器接收脉冲时间间隔序列Rn+1-Rn，并存储于有源定标器数据存储单元中；
步骤3：从星载雷达存储记录的数据中，读取星载雷达每一个发射脉冲的时刻，得到星载雷达发射脉冲时刻序列，将该发射脉冲时刻序列的相邻两个时刻依次做差值，得到星载雷达发射脉冲时间间隔序列tn+1-tn；Rn+1-Rn
步骤4：通过有源定标器内部的数据处理模块将步骤2得到的有源定标器接收脉冲时间间隔序列减去步骤3得到的对应的星载雷达发射脉冲时间间隔序列，得到一个差值序列，再将该差值序列乘以光速c，得到相邻两个发射脉冲所对应的星载雷达与有源定标器间的距离之差序列，将该距离之差序列存储在数据存储单元中；
步骤5：然后对步骤4得到的距离之差序列进行数据拟合得到以时间t为变量的一次直线函数，该直线函数的斜率拟合为；
步骤6：对步骤5得到的一次直线函数进行积分得到二次函数曲线即为有源定标器距离抛物线（公式（7）），并存储在有源定标器数据存储单元中；
步骤7：从星载雷达存储记录的数据中，读取星载雷达每一个发射脉冲的时刻和对应的每一个接收脉冲的时刻，用每一个接收脉冲时刻减去对应的发射脉冲时刻，得到星载雷达收发脉冲时间间隔序列rn-tn，将该序列也存储在有源定标器数据存储单元中；
步骤8：通过有源定标器内部的数据处理单元将步骤7得到的星载雷达收发脉冲时间间隔序列乘以光速c，再除以2，得到星载雷达测得的星载雷达至有源定标器的距离序列，然后对该距离序列进行数据拟合，得到一个以时间t为变量的二次函数曲线，通过与没有误差情况下星载雷达与有源定标器间的距离抛物线（公式(1)）的形状位置进行对比拟合得到公式（9）所示存在误差情况下的星载雷达距离抛物线函数，并存储于有源定标器数据存储单元中；
步骤9：在有源定标器内部的数据处理模块中，对比步骤6所得到的有源定标器距离抛物线和步骤8所得到的星载雷达距离抛物线的对称轴位置，根据其对称轴位置的不同，得到两个对称轴的位置之差，得出有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差；
步骤10：在地面使用频率计或频谱仪，直接测量有源定标器的时钟频率，并根据步骤9得到的有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差，求出星载雷达的时钟频率测量值；
步骤11：将步骤10得到的星载雷达时钟频率测量值与星载雷达时钟频率理论设计值相比较，得出星载雷达时钟漂移量，基于所述时钟漂移量校正星载雷达时钟系统。
本发明提出一种基于有源定标器的星载雷达时钟漂移校正方法，在星载雷达在轨测试中得到多次成功应用。实际测试结果表明：该种校正方法能够实时准确的测量出星载雷达时钟漂移量，测量精度优于1Hz。通过对星载雷达时钟漂移量的长期监测，能够更深入的掌握星载雷达的在轨运行状态，并为后续星载雷达的研制提供数据参考与方案指导。
上面所述是用于实现本发明的优选方案，因此本发明的范围不应该由此来限定。本领域的技术人员应该理解，在不脱离本发明范围内的任何修改或局部替换，均属于本发明权利要求来限定的范围。

Claims (9)

1.一种基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，包括以下步骤：
步骤（1）：通过有源定标器跟踪并接收星载雷达的脉冲信号，并记录每个接收脉冲的到达时刻，得到有源定标器接收脉冲时刻序列；
步骤（2）：通过有源定标器内部的数据处理模块对有源定标器接收脉冲时刻序列中的相邻两点求差值，得到有源定标器接收脉冲时间间隔序列Rn+1-Rn；
步骤（3）：通过有源定标器从星载雷达存储记录的数据中读取星载雷达每一个发射脉冲的时刻，得到星载雷达发射脉冲时刻序列，然后将该发射脉冲时刻序列中的相邻两点求差值，得到星载雷达发射脉冲时间间隔序列tn+1-tn；
步骤（4）：通过有源定标器内部的数据处理模块将步骤（2）和步骤（3）得到的时间间隔序列变换为星载雷达与有源定标器间的距离之差序列；
步骤（5）：对步骤（4）得到的距离之差序列进行数据拟合得到以时间t为变量的一次直线函数；
步骤（6）：对步骤（5）得到的一次直线函数进行积分得到有源定标器距离抛物线，星载雷达时钟漂移误差使所述有源定标器距离抛物线的对称轴位置较理想状态下距离抛物线的对称轴位置发生偏移；
步骤（7）：从星载雷达存储记录的数据中读取星载雷达每一个发射脉冲的时刻和对应的每一个接收脉冲的时刻，并用每一个接收脉冲时刻减去对应的发射脉冲时刻，得到星载雷达收发脉冲时间间隔序列rn-tn；
步骤（8）：通过有源定标器内部的数据处理单元将步骤（7）得到的星载雷达收发脉冲时间间隔序列变换为星载雷达与有源定标器的距离序列，并通过对该距离序列进行数据拟合得到星载雷达距离抛物线函数，星载雷达时钟漂移误差引起所述星载雷达距离抛物线较理想状态下的距离抛物线发生上下平移；
步骤（9）：在有源定标器内部的数据处理模块中，对比步骤（6）所得到的有源定标器距离抛物线和步骤（8）所得到的星载雷达距离抛物线的对称轴位置，根据其对称轴位置之差得出有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差；
步骤（10）：在地面使用频率计或频谱仪，直接测量所述有源定标器的时钟频率，并根据步骤（9）得到的有源定标器时钟频率和星载雷达时钟频率之间的时钟频率偏差，求出星载雷达的时钟频率测量值；
步骤（11）：将步骤（10）得到的星载雷达时钟频率测量值与星载雷达时钟频率理论设计值相比较，得出星载雷达时钟漂移量，基于所述时钟漂移量校正星载雷达时钟系统。
2.根据权利要求1所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，其中步骤（1）中，有源定标器通过其内部的高速数据采集模块和精密的GPS定时模块记录每个接收脉冲的到达时刻。
3.根据权利要求1所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，其中步骤（4）具体包括：通过有源定标器内部的数据处理模块将步骤（2）得到的有源定标器接收脉冲时间间隔序列Rn+1-Rn减去步骤（3）得到的对应的星载雷达发射脉冲时间间隔序列tn+1-tn，得到一个差值序列（Rn+1-Rn）－（tn+1-tn），再将该差值序列乘以光速c，得到相邻两个发射脉冲所对应的星载雷达与有源定标器间的距离之差序列。
4.根据权利要求1所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，其中步骤（5）中所述拟合得到的一次直线函数为：，其中，b表示由时钟漂移引入的误差，t为时间变量，a中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径。
5.根据权利要求1所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，其中步骤（6）中所述的有源定标器距离抛物线函数为：，b为时钟漂移引入的误差，t为时间变量，a中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径。
6.根据权利要求1所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，其中步骤（8）具体包括：通过有源定标器内部的数据处理单元将步骤（7）得到的星载雷达收发脉冲时间间隔序列乘以光速c再除以2，以得到星载雷达至有源定标器的距离序列，然后对该距离序列进行数据拟合，得到以时间t为变量的二次函数曲线，然后将该二次函数曲线与理想情况下星载雷达与有源定标器间的距离抛物线函数曲线进行拟合比对得到实测情况下的星载雷达距离抛物线函数。
7.根据权利要求6所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，所述实测情况下的星载雷达距离抛物线函数为：，其中，
，d表示由星载雷达时钟漂移误差造成的抛物线函数的上下平移量，a、C中各参数的含义为：ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；Re为地球半径，H为有源定标器相对星下点的高度。
8.根据权利要求1-7任一项所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，各步骤中的数据存储于有源定标器的数据存储单元中。
9.根据权利要求1-7任一项所述的基于有源定标器的星载雷达时钟漂移在轨校正方法，其特征在于，理想状态下星载雷达与有源定标器间的距离关系为：，式中ν为星载雷达飞行速度；ρ0为星载雷达至星下点的距离；ρ为星载雷达至有源定标器的距离；H为有源定标器相对星下点的高度；Re为地球半径；t为星载雷达飞行时间。
CN201410013535.2A 2014-01-13 2014-01-13 一种基于有源定标器的星载雷达时钟漂移在轨校正方法 Active CN103792519B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410013535.2A CN103792519B (zh) 2014-01-13 2014-01-13 一种基于有源定标器的星载雷达时钟漂移在轨校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410013535.2A CN103792519B (zh) 2014-01-13 2014-01-13 一种基于有源定标器的星载雷达时钟漂移在轨校正方法

Publications (2)

Publication Number Publication Date
CN103792519A CN103792519A (zh) 2014-05-14
CN103792519B true CN103792519B (zh) 2016-03-09

Family

ID=50668379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410013535.2A Active CN103792519B (zh) 2014-01-13 2014-01-13 一种基于有源定标器的星载雷达时钟漂移在轨校正方法

Country Status (1)

Country Link
CN (1) CN103792519B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104020489B (zh) * 2014-06-18 2017-01-18 中国海洋石油总公司 一种实现海上地球物理勘探的综合控制方法和装置
CN107179534B (zh) * 2017-06-29 2020-05-01 北京北科天绘科技有限公司 一种激光雷达参数自动标定的方法、装置及激光雷达
CN109001971B (zh) * 2018-07-25 2019-09-06 中国科学院国家空间科学中心 一种基于fpga的星载干涉成像高度计的自守时系统和方法
CN113671483B (zh) * 2021-08-19 2024-02-20 上海无线电设备研究所 一种基于秒脉冲的星载复合数据融合方法
CN114430308A (zh) * 2021-12-09 2022-05-03 西安昆仑工业（集团）有限责任公司 一种软件定时发送时间漂移导致时序偶发错位纠错方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6731236B1 (en) * 2003-06-11 2004-05-04 Honeywell International Inc. Methods and apparatus for optimizing interferometric radar altimeter cross track accuracy
CN103454620A (zh) * 2013-09-23 2013-12-18 中国科学院空间科学与应用研究中心 一种基于转发信号重建体制的跟踪型有源定标器

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6731236B1 (en) * 2003-06-11 2004-05-04 Honeywell International Inc. Methods and apparatus for optimizing interferometric radar altimeter cross track accuracy
CN103454620A (zh) * 2013-09-23 2013-12-18 中国科学院空间科学与应用研究中心 一种基于转发信号重建体制的跟踪型有源定标器

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Design of the interface of a calibration transponder and an;Xu Xi-Yu.etc;《eoscience and Remote Sensing Symposium》;20110729;全文 *
ＨＹ－２卫星雷达高度计时标偏差估算;汪栋等;《海洋学报》;20130930;第35卷(第5期);全文 *
Ocillator Clock Drift Compensation in Bistatic Interferometric SAR;Michael Eineder;《Geoscience and Remote Sensing Symposium》;20030725;全文 *
星载雷达高度计有源定标器回波模型;杜岳恒等;《遥感技术与应用》;20120430;第27卷(第2期);全文 *

Also Published As

Publication number Publication date
CN103792519A (zh) 2014-05-14

Similar Documents

Publication Publication Date Title
CN103792519B (zh) 一种基于有源定标器的星载雷达时钟漂移在轨校正方法
US7623065B2 (en) Refractivity retrieval via direct measurement of GNSS bending angle
CN101017202B (zh) 一种雷达高度表及采用该表对飞行器位置的测量方法
CN102753991A (zh) 短基线和超短基线相位图
CN109782289A (zh) 一种基于基线几何结构约束的水下航行器定位方法
RU2411538C2 (ru) Способ определения ошибки измерения скорости ла инерциальной навигационной системой и бортовой навигационный комплекс для его реализации
CN105116372A (zh) 一种浅海海底水平直线水听器阵朝向校准的方法
Cario et al. Analysis of error sources in underwater localization systems
RU2617830C1 (ru) Способ пассивной однопозиционной угломерно-разностно-доплеровской локации перемещающегося в пространстве радиоизлучающего объекта и радиолокационная система для реализации этого способа
CN103823209B (zh) 用于轻小型合成孔径雷达系统中低成本运动误差测量装置
RU2388008C1 (ru) Способ определения углового положения летательного аппарата по данным приемников спутниковой навигационной системы
Bahcivan et al. Radio Aurora Explorer: Mission science and radar system
Łabowski et al. Inertial navigation system for radar terrain imaging
RU2660159C1 (ru) Способ определения угла сноса летательного аппарата бортовой радиолокационной станцией
Roussel et al. Multi-scale volumetric soil moisture detection from GNSS SNR data: Ground-based and airborne applications
RU2483324C1 (ru) Способ навигации летательного аппарата по радиолокационным изображениям земной поверхности
Moccia et al. Bistatic synthetic aperture radar
Yan et al. Doppler Parameter Estimation Model Using Onboard Orbit Determination and Inter-satellite Distance Measurement for Spaceborne Bistatic SAR Real-time Imaging
Liu et al. The preliminary results about positioning accuracy of GF-3 SAR satellite system
RU2525650C2 (ru) Способ наведения летательных аппаратов на наземные объекты
Dai et al. Uniform calibration of networked radar based on unmanned aerial vehicle and GPS
Liu et al. An high-efficiency positioning method for ground testing SAR system
Carreno-Luengo et al. Submeter ocean altimetry with GPS L1 C/A signal
Yoshihara et al. Airborne GPS down-looking occultation experiments
Ruols et al. Development of a drone-based ground-penetrating radar system for efficient and safe 3D and 4D surveying of alpine glaciers

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1

Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS

Address before: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1

Patentee before: Space Science & Applied Research Centre, Chinese Academy of Sciences

CP01 Change in the name or title of a patent holder