CN109557537B - 雷达高度计频率基准漂移的在轨监视及测高数据修正方法 - Google Patents

雷达高度计频率基准漂移的在轨监视及测高数据修正方法 Download PDF

Info

Publication number
CN109557537B
CN109557537B CN201811469693.3A CN201811469693A CN109557537B CN 109557537 B CN109557537 B CN 109557537B CN 201811469693 A CN201811469693 A CN 201811469693A CN 109557537 B CN109557537 B CN 109557537B
Authority
CN
China
Prior art keywords
frequency reference
radar altimeter
time
orbit
counter
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
CN201811469693.3A
Other languages
English (en)
Other versions
CN109557537A (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 CN201811469693.3A priority Critical patent/CN109557537B/zh
Publication of CN109557537A publication Critical patent/CN109557537A/zh
Application granted granted Critical
Publication of CN109557537B publication Critical patent/CN109557537B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/882Radar or analogous systems specially adapted for specific applications for altimeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了雷达高度计频率基准漂移的在轨监视及测高数据修正方法,所述雷达高度计频率基准漂移的在轨监视方法用于利用守时计数器监视星载雷达高度计的频率基准漂移;所述方法包括:采集雷达高度计两个秒脉冲下降沿到来时刻的GNSS授时及守时计数器读数,由此计算守时计数器时钟的周期,再根据守时计数器时钟的标称周期计算雷达高度计频率基准的精度S,由此监视频率基准仪器的工作状态。通过本发明的方法,频率基准的性能可以在轨监测,并代入数据产品进行修正,不需要长期的定期在轨定标过程,节省大量的人力财力物力。

Description

雷达高度计频率基准漂移的在轨监视及测高数据修正方法
技术领域
本发明涉及微波遥感技术领域,特别涉及雷达高度计频率基准漂移的在轨监视及测高数据修正方法。
背景技术
雷达高度计是一种主要应用于海洋探测的微波遥感器,它可以测量仪器正下方海面的平均高度、海面的有效波高和后向散射系数,这些测量结果经过进一步处理可以获得海面高度、海面有效波高和海面风速等数据产品,这些产品在海洋学、地球物理学和军事上具有广泛而重要的应用。
雷达高度计使用的频率基准的准确度决定着雷达高度计的测高准确度。由于温度变化、器件老化等原因,频率基准的工作频率点会发生漂移,从而导致高度计的测高准确度发生变化。因此需要对频率基准的准确度进行长期在轨标定。
目前,雷达高度计频率基准的在轨标定包括利用海洋定标场或者有源定标器进行定期在轨标定等方法。
GNSS授时方法采用GNSS整秒时间和GNSS整秒脉冲对用户设备进行校时,该方法可使用户设备的时间与绝对时间的误差普遍小于100ns,该误差不随时间变化而改变。
发明内容
本发明的目的在于提出雷达高度计频率基准漂移的在轨监视及测高数据修正方法。
本发明的技术方案为:一种雷达高度计频率基准漂移的在轨监视方法,用于利用守时计数器监视星载雷达高度计的频率基准漂移;所述方法包括:采集雷达高度计两个秒脉冲下降沿到来时刻的GNSS授时及守时计数器读数,由此计算守时计数器时钟的周期,再根据守时计数器时钟的标称周期计算雷达高度计频率基准的精度S,由此监视频率基准仪器的工作状态。
作为上述方法的一种改进,所述方法具体包括:
步骤1)在雷达高度计某一秒脉冲下降沿到来时刻,记录GNSS授时整秒时间TN1,以及采用高度计守时计数器读数N1;在另一秒脉冲下降沿到来时刻,记录GNSS授时整秒时间TN2,以及对应的守时计数器读数N2;
步骤2)计算雷达高度计的守时计数器时钟的周期
Figure BDA0001890629690000021
Figure BDA0001890629690000022
步骤3)计算雷达高度计频率基准的精度S:
Figure BDA0001890629690000023
其中,Tp为雷达高度计守时计数器时钟标称周期;
步骤4)将步骤3)获得的精度S与雷达高度计频率基准的许可范围进行比较,若在许可范围内,则频率基准仪器在轨正常;否则,频率基准仪器在轨异常,将异常信息和精度S发送至卫星管理方。
本发明还提供了一种雷达高度计频率基准漂移的测高数据修正方法,当卫星管理方接收到异常信息和精度S时,所述方法包括:对雷达高度计测量的高度值进行修正:
H'=H(1+S) (17)
其中,H为雷达高度计地面重跟踪处理获得的高度值,H'为修正频率基准漂移后的高度值。
本发明的优点在于:
通过本发明的方法,频率基准的性能可以在轨监测,并代入数据产品进行修正,不需要长期的定期在轨定标过程,节省大量的人力财力物力。
附图说明
图1(a)为雷达高度计发射和接收线性调频脉冲的示意图;
图1(b)为雷达高度计单个脉冲测高示意图;
图2为本发明的频率性能监测原理的示意图。
具体实施方式
下面结合附图对本发明进行详细的说明。
首先给出频率基准精度与雷达高度计测量的高度的关系。
图1(a)和图1(b)给出了雷达高度计单个脉冲测高过程示意。图1(a)表示雷达高度计在E1时刻向海面发射了一个线性调频脉冲,然后在E2时刻产生一个接收参考线性调频脉冲作为回波本振参考;E1、E2时间差采用周期为T的参考计数时钟进行计数。图1(b)表示回波信号与图1(a)中的接收参考Chirp混频后,其差频Δf被ADC采样,然后经FFT运算获得具体数值ΔfFFT
由图1(b)可以看出回波信号与接收参考Chirp的相对时延满足时频关系对:
Δt=ΔfFFT/K (1)
式中,Δt表示相对时延,K表示线性调频斜率。
因此,雷达高度计的测量高度H可以表示为:
Figure BDA0001890629690000031
由公式(2)可知,参考源的准确度可能从以下三个方面决定着雷达高度计的测高精度:a)、参考计数时钟周期T的准确度;b)、ADC参考时钟的准确度,进而影响数值ΔfFFT的准确度;c)、线性调频斜率K的准确度。下面对其进行详细讨论。
以HY-2A雷达高度计为例,Chirp信号采用直接数字频率合成(Direct DigitalFrequency Synthesizer,DDFS)方法产生。DDFS产生Chirp信号的机理可以简述为:1、获得理想的以标称时宽TDDS为采样周期的数字Chirp信号;2、将数字Chirp信号以采样周期为TDDS进行数模DAC转换。设参与DDFS的数字Chirp信号样本值为S[1]、S[2]、……S[n],当以标称时宽TDDS合成Chirp信号时,其带宽为B,则当以实际时宽TD'DS合成时,其带宽为:
Figure BDA0001890629690000032
因此,实际调频斜率K’与标称调频斜率K的关系为:
Figure BDA0001890629690000033
设ADC工作标称频率为FADC,实际频率为FA'DC,则实测值ΔfFFT与标称值Δf的关系可以表示为:
Figure BDA0001890629690000034
又由于时频关系对Δf=K'.Δt,代入公式(5)可得:
Figure BDA0001890629690000041
因此,实测值ΔfFFT对应的标称相对时延Δt,实测时延Δt'分别为:
Δt=ΔfFFT/K (7)
Figure BDA0001890629690000042
设参考计数时钟周期实际值为T’,雷达高度计的参考计数时钟、DDS参考时钟、ADC时钟都通过频率基准派生获得,因此满足关系:
Figure BDA0001890629690000043
上式中,
Figure BDA0001890629690000044
为频率基准的相对精度。
联立公式(2)、(7)、(8)、(9)可得测量高度H与实际高度H’的关系:
Figure BDA0001890629690000045
公式(10)表明,频率基准的准确度直接决定着高度测量结果的准确度。若频率准确度为10-9,轨道高度为900公里左右,则由频率准确度引起的误差为1mm左右,合成孔径雷达高度计的测量精度约为7mm,两者处于同一数量级别,若对频率基准漂移不加以修正,则会影响合成孔径雷达高度计的测量精度。
雷达高度计通过GPS授时获得准确的协调世界时间。GPS授时的时间准确度是由授时方法决定的,不会产生时间累积误差;若在高度计内部设立一个守时计数器,该计数器采用频率基准进行计时且不参与校时,则计数器的误差不断累积。如图2所示,秒脉冲及广播来自于GPS授时系统,高度计守时计数时钟由频率基准派生。因此,守时计数时钟周期的准确度即可表示频率基准的准确度。
基于GPS授时的频率基准的频率准确度的统计监测方法具体做法如下:
在某一秒脉冲下降沿到来时刻,记录GPS授时整秒时间TN1,以及采用高度计守时计数器读数N1;在另一秒脉冲下降沿到来时刻,记录GPS授时整秒时间TN2,以及对应的守时计数器读数N2。设高度计守时计数器时钟标称周期为Tp,则由守时计数时钟测量获得两次纪录之间的时间差为ΔTpulse=(N2-N1)*Tp;对应由GPS授时测量获得的时间差为ΔTGPS=TN2-TN1。
GPS秒脉冲下降沿能够提供的精度的典型值为100ns。则两次纪录之间的GPS授时时间差的误差ξ(ΔTGPS)≤200ns。
记守时计数器时钟周期误差为ξTp,则两次纪录之间的守时计数累积误差为ξ(ΔTpulse)=(N2-N1)×ξ(Tp)。当ξ(ΔTpulse)远大于ξ(ΔTGPS)时,可以用下式估计守时计数器时钟的周期:
Figure BDA0001890629690000051
其中,
Figure BDA0001890629690000052
表示采用累积法对Tp的估计。
公式(11)的估计方法引入的绝对误差为:
Figure BDA0001890629690000053
公式(12)两端各除以TP,得到估计方法的相对误差为:
Figure BDA0001890629690000054
根据公式(13)可以获得两次纪录之间需要累积的时间:
Figure BDA0001890629690000055
令ηTP≤1×10-10,ξ(ΔTGPS)=200ns,可得需要的累积时间Tacc≥2000s,这表示当累积时间大于2000s时,采用以上累积对比监测估计守时计数时钟周期的准确度为1×10-10
采用公式(11)获得守时计数器时钟标称周期的估计
Figure BDA0001890629690000056
之后,就可以估计雷达高度计频率基准的精度S,监测S是否在许可范围内,实现对频率基准的在轨性能监测;
Figure BDA0001890629690000061
将获得的S与频率基准的许可范围进行对比,若在许可范围内,则表明频率基准仪器的在轨性能良好;若超出许可范围,则表明频率基准仪器的在轨性能恶化,需要报告卫星管理方,且对测高数据进行修正。具体修正方法为:
记H为地面重跟踪处理获得的高度,H'为修正频率基准漂移后的高度,则由公式(9)、(10)可得
Figure BDA0001890629690000062
公式(15)代入公式(10),获得最终的修正表达式:
Figure BDA0001890629690000063
修正高度H'还可以表示为:
H'=H(1+S) (17)
假定GNSS授时误差为100ns,两次整秒脉冲间隔时间为2天,则以上估计方法的精度可达:
Figure BDA0001890629690000064
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (2)

1.一种雷达高度计频率基准漂移的在轨监视方法,用于利用守时计数器监视星载雷达高度计的频率基准漂移;所述方法包括:采集雷达高度计两个秒脉冲下降沿到来时刻的GNSS授时及守时计数器读数,由此计算守时计数器时钟的周期,再根据守时计数器时钟的标称周期计算雷达高度计频率基准的精度S,由此监视频率基准仪器的工作状态;
所述方法具体包括:
步骤1)在雷达高度计某一秒脉冲下降沿到来时刻,记录GNSS授时整秒时间TN1,以及对应的守时计数器读数N1;在另一秒脉冲下降沿到来时刻,记录GNSS授时整秒时间TN2,以及对应的守时计数器读数N2;
步骤2)计算雷达高度计的守时计数器时钟的周期
Figure FDA0002436638090000011
Figure FDA0002436638090000012
步骤3)计算雷达高度计频率基准的精度S:
Figure FDA0002436638090000013
其中,Tp为雷达高度计守时计数器时钟标称周期;
步骤4)将步骤3)获得的精度S与雷达高度计频率基准的许可范围进行比较,若在许可范围内,则频率基准仪器在轨正常;否则,频率基准仪器在轨异常,将异常信息和精度S发送至卫星管理方。
2.根据权利要求1所述的雷达高度计频率基准漂移的在轨监视方法,其特征在于,当卫星管理方接收到异常信息和精度S时,所述方法还包括:对雷达高度计测量的高度值进行修正:
H'=H(1+S) (17)
其中,H为雷达高度计地面重跟踪处理获得的高度值,H'为修正频率基准漂移后的高度值。
CN201811469693.3A 2018-11-28 2018-11-28 雷达高度计频率基准漂移的在轨监视及测高数据修正方法 Active CN109557537B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811469693.3A CN109557537B (zh) 2018-11-28 2018-11-28 雷达高度计频率基准漂移的在轨监视及测高数据修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811469693.3A CN109557537B (zh) 2018-11-28 2018-11-28 雷达高度计频率基准漂移的在轨监视及测高数据修正方法

Publications (2)

Publication Number Publication Date
CN109557537A CN109557537A (zh) 2019-04-02
CN109557537B true CN109557537B (zh) 2020-08-04

Family

ID=65868767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811469693.3A Active CN109557537B (zh) 2018-11-28 2018-11-28 雷达高度计频率基准漂移的在轨监视及测高数据修正方法

Country Status (1)

Country Link
CN (1) CN109557537B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110471087A (zh) * 2019-08-13 2019-11-19 北京电子工程总体研究所 一种空间飞行器的时间漂移计算方法和系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101383653A (zh) * 2008-10-30 2009-03-11 航天东方红卫星有限公司 一种星上自主时统方法
CN101562451A (zh) * 2009-05-27 2009-10-21 西安华伟电力电子技术有限责任公司 二级频标的精密驯服保持方法
CN101799658A (zh) * 2010-02-24 2010-08-11 华中科技大学 一种gps校准的守时钟
CN102436174A (zh) * 2011-10-26 2012-05-02 东莞市泰斗微电子科技有限公司 一种守时设备晶振频率驯服方法及相应装置
CN103269262A (zh) * 2013-04-01 2013-08-28 柳州市达迪通信设备有限公司 一种时间同步装置的守时方法
CN107145061A (zh) * 2017-06-09 2017-09-08 广州北极瑞光电子科技有限公司 卫星授时信号的智能抗干扰技术在时间同步中的防御方法
JP6278643B2 (ja) * 2013-09-09 2018-02-14 三菱電機株式会社 子局装置、親局装置、制御装置、通信システムおよび時刻同期方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101383653A (zh) * 2008-10-30 2009-03-11 航天东方红卫星有限公司 一种星上自主时统方法
CN101562451A (zh) * 2009-05-27 2009-10-21 西安华伟电力电子技术有限责任公司 二级频标的精密驯服保持方法
CN101799658A (zh) * 2010-02-24 2010-08-11 华中科技大学 一种gps校准的守时钟
CN102436174A (zh) * 2011-10-26 2012-05-02 东莞市泰斗微电子科技有限公司 一种守时设备晶振频率驯服方法及相应装置
CN103269262A (zh) * 2013-04-01 2013-08-28 柳州市达迪通信设备有限公司 一种时间同步装置的守时方法
JP6278643B2 (ja) * 2013-09-09 2018-02-14 三菱電機株式会社 子局装置、親局装置、制御装置、通信システムおよび時刻同期方法
CN107145061A (zh) * 2017-06-09 2017-09-08 广州北极瑞光电子科技有限公司 卫星授时信号的智能抗干扰技术在时间同步中的防御方法

Also Published As

Publication number Publication date
CN109557537A (zh) 2019-04-02

Similar Documents

Publication Publication Date Title
US7548600B2 (en) Apparatus and method for compensating the drift of a local clock used as sampling frequency
US8997550B2 (en) Method and system for correcting for temperature variations in ultrasonic testing systems
CA2425250A1 (en) Time-base generator with self-compensating control loop
CN109557537B (zh) 雷达高度计频率基准漂移的在轨监视及测高数据修正方法
JP2010019728A (ja) 測位方法および測位装置
US9958551B1 (en) Method and device for determining availability of ARAIM system
CN116125512A (zh) 顾及时钟频率时变特性的ppp自适应钟差模型估算方法
CN100516900C (zh) 一种长波传播时延修正量的测量方法
WO2022001368A1 (zh) 雷达及雷达系统
Shmaliy et al. A thinning algorithm for GPS-based unbiased FIR estimation of a clock TIE model
KR100498185B1 (ko) 아이지에스를 이용한 지피에스 위치오차 보정정보제공방법 및 그를 이용한 지피에스 위치오차 보정방법
US8842037B2 (en) High frequency surfacewave radar
CN110244332A (zh) 一种天线阵大气相位扰动修正方法及系统
Webb A comparison of SEASAT 1 altimeter measurements of wave height with measurements made by a pitch‐roll buoy
US20130013254A1 (en) Self temperature-compensated high precision event timer using standard time reference frequency and its method
JP2010249578A (ja) 変換係数導出装置および変換係数導出方法
Tran et al. Assessment of the cycle-to-cycle noise level of the Geosat Follow-On, TOPEX, and Poseidon altimeters
CN115856793A (zh) 一种雷达信号频率偏差估计补偿方法
CN211180593U (zh) 时频检测设备
CA1114480A (en) Method and apparatus for automatically calibrating a radio altimeter
CN113608221A (zh) 一种加速度辅助机车雷达测速滤波方法
US6831596B1 (en) Calibrating the sampling frequency of a GPS receiver
CN112485813A (zh) Glonass测站间非组合测距码频间偏差校正方法及系统
EP2084556B1 (en) Weather radar signal processing
Mahyun et al. The Radar-Rainfall Relationship for Northern Region of Peninsular Malaysia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant