CN100507599C - 一种海陆兼容的卫星雷达高度计跟踪器 - Google Patents
一种海陆兼容的卫星雷达高度计跟踪器 Download PDFInfo
- Publication number
- CN100507599C CN100507599C CNB2006101144960A CN200610114496A CN100507599C CN 100507599 C CN100507599 C CN 100507599C CN B2006101144960 A CNB2006101144960 A CN B2006101144960A CN 200610114496 A CN200610114496 A CN 200610114496A CN 100507599 C CN100507599 C CN 100507599C
- Authority
- CN
- China
- Prior art keywords
- βlv
- boqi
- processing unit
- discr
- error
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种海陆兼容的卫星雷达高度计跟踪器，包括第一FFT功率谱分析器、MLE误差鉴别器和第一α-β滤波器，还包括第二FFT功率谱分析器、OCOG误差鉴别器、第二α-β滤波器和数据合成处理单元；其中，第一FFT功率谱分析器电连接到MLE误差鉴别器、第一α-β滤波器和数据合成处理单元上，MLE误差鉴别器分别与第一α-β滤波器和数据合成处理单元电连接，第一α-β滤波器与数据合成处理单元电连接；第二FFT功率谱分析器电连接到OCOG误差鉴别器、第二α-β滤波器和数据合成处理单元上，OCOG误差鉴别器分别与第二α-β滤波器和数据合成处理单元电连接，第二α-β滤波器与数据合成处理单元电连接；第一α-β滤波器还与第二α-β滤波器电连接。
Description
技术领域
本发明涉及一种卫星雷达高度计跟踪器，特别是涉及一种最大似然跟踪算法(maximum likelihood estimation，MLE)和重心跟踪算法(off center of gravity,OCOG)两种算法并行工作的雷达高度计跟踪器。
背景技术
卫星雷达高度计是用于海洋动力环境测量的重要遥感器，其直接测量的产品包括卫星到海面的高度，海面有效波高和海面后向散射系数，雷达高度计测量数据的进一步反演结果，可以应用于海洋地球物理学、海洋动力学、海洋气候与环境、海冰监测等方面。除此以外，雷达高度计还可以测量海冰、海陆交界处地形和陆地地表高度。自1973年Skylab试验以来，以美国和欧洲为代表，雷达高度计卫星已经发展成为系列化的应用卫星(如GEOSAT、Topex/poseidon、Jason-1、ERS-1、ERS-2和Envisat等)，技术日趋成熟。
现有在轨运行的雷达高度计跟踪器在同一时刻只有一种跟踪算法在工作，如GEOSAT和Topex/poseidon高度计只采用了门跟踪算法。门跟踪算法测量精度较高但健壮性较差，比较适合海洋测量，而在海陆交界、陆地和海冰的测量方面能力较差，容易失锁。Jason-1高度计跟踪器是一个基于MLE算法的跟踪器，图2为基于MLE的跟踪器结构图，该跟踪器由FFT功率谱分析器、MLE误差鉴别器、α-β滤波器和时序生成器组成。其中，FFT功率谱分析器分别与MLE误差鉴别器和α-β滤波器电连接，MLE误差鉴别器与α-β滤波器之间电连接，α-β滤波器还与时序生成器电连接。该跟踪器从高度计的接收机接收I、Q信号，生成波形数据、MLE跟踪参数、信号带宽、高度值、浪高值、AGC和工作时序等多种数据，并产生距离精调和距离粗调值，用于对跟踪器本身的调节。
Envisat高度计只采用了OCOG算法。OCOG算法跟踪精度较差但健壮性较强，比较适合海陆交界、陆地和海冰的测量；用OCOG算法进行海洋测量，需要在地面进行二次跟踪，但是存在累积误差，降低了海洋测量精度。图3为基于OCOG的跟踪器结构图。该跟踪器由FFT功率谱分析器、OCOG误差鉴别器、α-β滤波器和时序生成器组成。
ERS-1和ERS-2高度计在海洋测量时采用了准最大似然算法，而在海陆交界、陆地和海冰测量时采用了OCOG算法。这是一种比较好的方式，问题是这两种跟踪算法不是同时工作的，两种算法相互切换，特别是由准最大似然算法向OCOG算法切换时，再次捕获回波会使跟踪器一段时间内失锁，从而造成测量的不连续性。
发明内容
本发明的目的是克服现有的雷达高度计跟踪器无法很好地兼容于海洋和陆地的缺陷，从而提供一种既可用于海洋测量，又可同时用于陆地测量，且测量精度高、健壮性强的雷达高度计跟踪器。
为了实现上述目的，本发明提供了一种海陆兼容的卫星雷达高度计跟踪器，包括第一FFT功率谱分析器1、MLE误差鉴别器3和第一α-β滤波器5，其特征在于，还包括第二FFT功率谱分析器2、OCOG误差鉴别器4、第二α-β滤波器6和数据合成处理单元7；其中，所述的第一FFT功率谱分析器1电连接到所述的MLE误差鉴别器3、所述的第一α-β滤波器5和所述的数据合成处理单元7上，所述的MLE误差鉴别器3分别与所述的第一α-β滤波器5和所述的数据合成处理单元7电连接，所述的第一α-β滤波器5与所述的数据合成处理单元7电连接；所述的第二FFT功率谱分析器2电连接到所述的OCOG误差鉴别器4、所述的第二α-β滤波器6和所述的数据合成处理单元7上，所述的OCOG误差鉴别器4分别与所述的第二α-β滤波器6和所述的数据合成处理单元7电连接，所述的第二α-β滤波器6与所述的数据合成处理单元7电连接；所述的第一α-β滤波器5还与所述的第二α-β滤波器6电连接；
所述OCOG误差鉴别器4生成第二鉴别器状态和第二信号带宽，当所述第二鉴别器状态表明OCOG处于锁定状态且所述第二信号带宽为320MHz时，所述卫星雷达高度计跟踪器位于第二工作模式，否则位于第一工作模式；在所述第一工作模式下，所述数据合成处理单元7生成的第一状态控制信号表示待机状态，生成的第二状态控制信号表示主动工作状态；在所述第二工作模式下，所述数据合成处理单元7生成的第一状态控制信号表示工作状态，生成的第二状态控制信号表示被动工作状态；
所述第一α-β滤波器5接收到所述数据合成处理单元7生成的第一状态控制信号后，如果所述第一状态控制信号是待机状态，此时所述第一α-β滤波器5接收来自所述第二α-β滤波器6的第二高度信号和第二AGC，并以此值作为工作状态的初始值，同时所述第一α-β滤波器5将输出的第一距离精调信号设定为0，并处于待机状态；当所述第一状态控制信号为工作状态时，所述第一α-β滤波器5接收所述MLE误差鉴别器3的误差信号，将滤波处理后输出的第一高度信号和第一AGC的值送到所述数据合成处理单元7和第二α-β滤波器6，将滤波处理后输出的浪高值送到所述数据合成处理单元7，将滤波处理后输出的第一高度信号分成第一距离粗调和第一距离精调两部分，其中所述第一距离粗调送到所述数据合成处理单元7，所述第一距离精调送到所述第一FFT功率谱分析器1；
所述第二α-β滤波器6接收所述数据合成处理单元7生成的第二状态控制信号，当所述第二状态控制信号为被动工作状态时，所述第二α-β滤波器6接收来自所述OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的第二高度信号和第二AGC送到所述数据合成处理单元7；滤波处理后输出的第二高度信号分成第二距离粗调和第二距离精调两部分，其中所述第二距离粗调送到所述数据合成处理单元7，所述第二距离精调送到所述第二FFT功率谱分析器2；完成上述操作后，所述第二α-β滤波器6接收来自所述第一α-β滤波器5的第一高度信号和第一AGC，并复位所述第二α-β滤波器6的当前输出，使得第二高度信号等于第一高度信号，且第二AGC等于第一AGC；
当所述第二状态控制信号为主动工作状态时，所述第二α-β滤波器6接收所述OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的第二高度信号和第二AGC送到所述数据合成处理单元7和所述第一α-β滤波器5；滤波处理后输出的第二高度信号分成第二距离粗调和第二距离精调两部分，其中所述的第二距离粗调送到所述数据合成处理单元7，所述第二距离精调送到所述第二FFT功率谱分析器2。
上述技术方案中，所述的第二FFT功率谱分析器2、所述的OCOG误差鉴别器4和所述的第二α-β滤波器6在DSP芯片上编程实现。
上述技术方案中，所述的第一FFT功率谱分析器1、所述的MLE误差鉴别器3和所述的第一α-β滤波器5在DSP芯片上编程实现。
上述技术方案中，所述的数据合成处理单元7在DSP芯片和FPGA电路上编程实现。
本发明的优点在于：
1、本发明的海陆兼容的卫星雷达高度计跟踪器既适合于海洋测量，又适合于陆地和海冰的测量。
2、本发明的海陆兼容的卫星雷达高度计跟踪器具备跟踪精度高、跟踪健壮性强的优点，对提高卫星雷达高度计跟踪器的跟踪性能具有重要意义。
3、本发明的海陆兼容的卫星雷达高度计跟踪器在测量海陆交界处时，由于MLE和OCOG并行工作，不存在由于两种算法切换造成的跟踪测量的不连续。
附图说明
图1本发明跟踪器的结构图；
图2为基于MLE算法的跟踪器的结构图；
图3为基于OCOG算法的跟踪器的结构图。
1 第一FFT功率谱分析器 2 第二FFT功率谱分析器
3 MLE误差鉴别器 4 OCOG误差鉴别器
5 第一α-β滤波器 6 第二α-β滤波器
7 数据合成处理单元 8 时序生成器
具体实施方式
下面结合附图和具体实施方式对本发明的海陆兼容的卫星雷达高度计跟踪器进行说明。
如图1所示，本发明的跟踪器由3部分组成。其中第一部分完成MLE跟踪功能，第二部分完成OCOG跟踪功能，第三部分完成两个跟踪器的数据合成处理。所述的第一部分包括第一FFT功率谱分析器1、MLE误差鉴别器3和第一α-β滤波器5，所述的第二部分包括第二FFT功率谱分析器2、OCOG误差鉴别器4和第二α-β滤波器6，所述的第三部分包括数据合成处理单元7。第一FFT功率谱分析器1与高度计的接收机相连，并电连接到MLE误差鉴别器3、第一α-β滤波器5和数据合成处理单元7上，MLE误差鉴别器3分别与第一α-β滤波器5和数据合成处理单元7电连接，第一α-β滤波器5与数据合成处理单元7电连接。第二FFT功率谱分析器2也与高度计的接收机相连，并电连接到OCOG误差鉴别器4、第二α-β滤波器6和数据合成处理单元7上，OCOG误差鉴别器4分别与第二α-β滤波器6和数据合成处理单元7电连接，第二α-β滤波器6与数据合成处理单元7电连接，所述的第一α-β滤波器5还与第二α-β滤波器6电连接。
其中第一FFT功率谱分析器1，接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第一α-β滤波器5传送来的距离精调值①对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据①。第一FFT功率谱分析器1将波形数据①分别传送给数据合成单元7和MLE误差鉴别器3。
第二FFT功率谱分析器2，接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第二α-β滤波器6传送来的距离精调值②对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据②。第二FFT功率谱分析器2将波形数据②分别传送给数据合成单元7和OCOG鉴别器4。
MLE误差鉴别器3接收数据合成处理单元7的状态控制①，当状态控制①为待机状态时，MLE误差鉴别器3处于待机状态。当状态控制①为工作状态时进行如下操作：MLE误差鉴别器3对来自第一FFT功率谱分析器1的波形数据①进行MLE估计处理，并将输出的误差信号发送到第一α-β滤波器5；MLE误差鉴别器3根据MLE估计处理的结果，判定MLE处于跟踪锁定状态还是失锁状态，并将这一鉴别器状态①传送给数据合成处理单元7；MLE误差鉴别器3根据MLE估计的结果选定高度计工作的信号带宽①(320MHz、80MHz和20MHz中的一种)，并将信号带宽①送给数据合成处理单元7。
OCOG误差鉴别器4接收来自第二FFT功率谱分析器2的波形数据②进行OCOG处理，并将输出的误差信号发送到第二α-β滤波器6；OCOG误差鉴别器4根据OCOG处理结果，判定OCOG处于跟踪锁定状态还是失锁状态，并将这一鉴别器状态②传送给数据合成处理单元7；OCOG误差鉴别器4将处理后选定的信号带宽②(320MHz、80MHz和20MHz中的一种)传送给数据合成处理单元7。
第一α-β滤波器5接收数据合成处理单元7的状态控制①，如果状态控制①是待机状态，此时第一α-β滤波器5接收来自第二α-β滤波器6的高度②和AGC②，并以此值作为工作状态的初始值，同时第一α-β滤波器5将输出的距离精调①设定为0，并处于待机状态。当状态控制①为工作状态时进行如下操作：第一α-β滤波器5接收MLE误差鉴别器3的误差信号，将滤波处理后输出的高度①和AGC①值分别送到数据合成处理单元7和第二α-β滤波器6；滤波处理后输出的浪高值送到数据合成处理单元7；滤波处理后输出的高度①分成距离粗调①和距离精调①两部分，其中距离粗调①送到数据合成处理单元7，距离精调①送到第一FFT功率谱分析器1。
第二α-β滤波器6接收数据合成处理单元7的状态控制②，当状态控制②为被动工作状态时进行入下操作：第二α-β滤波器6接收来自OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的高度②和AGC②送到数据合成处理单元7；滤波处理后输出的高度②分成距离粗调②和距离精调②两部分，其中距离粗调②送到数据合成处理单元7，距离精调②送到第二FFT功率谱分析器2；完成上述操作后，第二α-β滤波器6接收来自第一α-β滤波器5的高度①和AGC①，并复位第二α-β滤波器6的当前输出，使得高度②＝高度①，且AGC②＝AGC①。
第二α-β滤波器6接收数据合成处理单元7的状态控制②，当状态控制②为主动工作状态时进行入下操作：第二α-β滤波器6接收来OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的高度②和AGC②分别送到数据合成处理单元7和第一α-β滤波器5；滤波处理后输出的高度②分成距离粗调②和距离精调②两部分，其中距离粗调②送到数据合成处理单元7，距离精调②送到第二FFT功率谱分析器2。
上述的第一FFT功率谱分析器1、MLE误差鉴别器3、第一α-β滤波器5、第二FFT功率谱分析器2、OCOG误差鉴别器4和第二α-β滤波器6均可在一个DSP芯片上通过编程实现。
数据合成处理单元7接收来自OCOG误差鉴别器4的鉴别器状态②和信号带宽②，如果鉴别器状态②表明OCOG处于锁定状态且信号带宽②为320MHz，此时设定系统工作模式为模式2，否则系统工作模式为模式1。在工作模式1情况下，数据合成处理单元7送出的状态控制①为待机状态，而状态控制②为主动工作状态；在工作模式2情况下，数据合成处理单元7送出的状态控制①为工作状态，而状态控制②为被动工作状态。
数据合成处理单元7在DSP芯片和FPGA电路上通过编程实现。
在工作模式1情况下，数据合成处理单元7输出的波形数据为接收来自第二FFT功率谱分析器2的波形数据②；数据合成处理单元7输出的AGC值为来自第二α-β滤波器6的AGC②；数据合成处理单元7根据来自第二α-β滤波器6的距离粗调②生成对应的工作时序；数据合成处理单元7输出的跟踪高度为来自第二α-β滤波器6的高度②；数据合成处理单元7输出的OCOG跟踪参数为来自OCOG误差鉴别器4的鉴别器状态②；数据合成处理单元7输出的信号带宽为来自OCOG误差鉴别器4的信号带宽②。
在工作模式2情况下，数据合成处理单元7输出的波形数据为接收来自第一FFT功率谱分析器1的波形数据①；数据合成处理单元7输出的AGC值为来自第一α-β滤波器5的AGC①；数据合成处理单元7根据来自第一α-β滤波器5的距离粗调①生成对应的工作时序；数据合成处理单元7输出的信号带宽为来自MLE误差鉴别器3的信号带宽①；数据合成处理单元7输出的跟踪高度为来自第一α-β滤波器5的高度①；数据合成处理单元7输出的MLE跟踪参数为来自MLE误差鉴别器3的鉴别器状态①；数据合成处理单元7输出的跟踪浪高为来自第一α-β滤波器5的浪高；数据合成处理单元7输出的OCOG跟踪参数为来自OCOG误差鉴别器4的鉴别器状态②。
在系统处于工作模式2的时候，当数据合成处理单元7检测到鉴别器状态①为失锁状态时，系统将工作状态切换到工作模式1。
本发明涉及的海陆兼容跟踪器采用MLE和OCOG并行工作的方式，在海洋测量时采用MLE的跟踪结果，在海陆交界、陆地和海冰的测量方面采用OCOG的跟踪结果。由于OCOG跟踪器的健壮性强，所以该跟踪器参与全程跟踪。当OCOG进入海洋测量区域时，其工作带宽会调整到最大工作带宽(320MHz)并且跟踪器处于锁定状态，此时启动MLE跟踪器工作；当进入海陆交界、陆地和海冰的测量区域时，MLE跟踪失锁，此时以OCOG工作为主。具体的工作过程如下：
步骤1，工作模式1：
该模式下数据合成单元7输出的状态控制①为待机状态，状态控制②为主动工作状态。
该模式下第二FFT功率谱分析器2接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第二α-β滤波器6传送来的距离精调值②对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据②。第二FFT功率谱分析器2将波形数据②分别传送给数据合成单元7和OCOG鉴别器4。
第一FFT功率谱分析器1，接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第一α-β滤波器5传送来的距离精调①对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据①。第一FFT功率谱分析器1将波形数据①分别传送给数据合成单元7和MLE误差鉴别器3。
OCOG误差鉴别器4接收来自第二FFT功率谱分析器2的波形数据②进行OCOG处理，并将生成的输出误差信号送到第二α-β滤波器6；OCOG误差鉴别器4根据处理结果，判定OCOG处于跟踪锁定状态还是失锁状态，并将这一鉴别器状态②传送给数据合成处理单元7；OCOG误差鉴别器4将处理后选定的信号带宽②(320MHz、80MHz和20MHz中的一种)传送给数据合成处理单元7。
第二α-β滤波器6接收来自OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的高度②和AGC②分别送到数据合成处理单元7和第一α-β滤波器5；滤波处理后输出的高度②分成距离粗调②和距离精调②两部分，其中距离粗调②送到数据合成处理单元7，距离精调②送到第二FFT功率谱分析器2。
数据合成处理单元7输出的波形数据为来自第二FFT功率谱分析器2的波形数据②；数据合成处理单元7输出的AGC值为来自第二α-β滤波器6的AGC②；数据合成处理单元7输出的信号带宽为来自OCOG误差鉴别器4的信号带宽②；数据合成处理单元7根据来自第二α-β滤波器6的距离粗调②生成对应的工作时序；数据合成处理单元7输出的跟踪高度为来自第二α-β滤波器6的高度②；数据合成处理单元7输出的OCOG跟踪参数为来自OCOG误差鉴别器4的鉴别器状态②。
第一α-β滤波器5接收来自第二α-β滤波器6的高度②和AGC②，并以此作为工作状态的初始值，此时第一α-β滤波器5输出的距离精调①设定为0。
步骤2、工作模式2：
该模式下数据合成单元7输出的状态控制①为工作状态，状态控制②为被动工作状态。
该模式下第一FFT功率谱分析器1，接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第一α-β滤波器5传送来的距离精调①对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据①。第一FFT功率谱分析器1将波形数据①分别传送给数据合成单元7和MLE误差鉴别器3。
第二FFT功率谱分析器2，接收来自高度计接收机采集的正交信号I和Q，并构成复信号I-jQ。利用第二α-β滤波器6传送来的距离精调②对I-jQ信号进行相位旋转，对旋转后的I-jQ信号进行FFT变换求功率谱，这个功率谱即为波形数据②。第二FFT功率谱分析器2将波形数据②分别传送给数据合成单元7和OCOG鉴别器4。
MLE误差鉴别器3对来自第一FFT功率谱分析器1的波形数据①进行MLE估计处理，并将生成的输出误差信号送到第一α-β滤波器5；MLE误差鉴别器3根据处理结果，判定MLE处于跟踪锁定状态还是失锁状态，并将这一鉴别器状态①传送给数据合成处理单元7；MLE误差鉴别器3同时将当前信号带宽①(320MHz、80MHz和20MHz中的一种)送给数据合成处理单元7。
OCOG误差鉴别器4接收来自第二FFT功率谱分析器2的波形数据②进行OCOG处理，并将生成的输出误差信号送到第二α-β滤波器6；OCOG误差鉴别器4根据处理结果，判定OCOG跟踪器处于跟踪锁定状态还是失锁状态，并将这一鉴别器状态②传送给数据合成处理单元7；OCOG误差鉴别器4同时将当前信号带宽②(320MHz、80MHz和20MHz中的一种)送给数据合成处理单元7。
第一α-β滤波器5接收MLE误差鉴别器3的误差信号，进行滤波处理后输出的高度①和AGC①分别送到数据合成处理单元7和第二α-β滤波器6；第一α-β滤波器5输出的浪高送到数据合成处理单元7；第一α-β滤波器5输出的高度①分成距离粗调①和距离精调①两部分，其中距离粗调①送到数据合成处理单元7，距离精调①送到第一FFT功率谱分析器1。
第二α-β滤波器6接收来自OCOG误差鉴别器4输出的误差信号并进行滤波，滤波处理后输出的高度②和AGC②送到数据合成处理单元7；滤波处理后输出的高度②分成距离粗调②和距离精调②两部分，其中距离粗调②送到数据合成处理单元7，距离精调②送到第二FFT功率谱分析器2；完成上述操作后，第二α-β滤波器6接收来自第一α-β滤波器5的高度①和AGC①，并复位第二α-β滤波器6的当前输出，使得高度②＝高度①，且AGC②＝AGC①。
数据合成处理单元7输出的跟踪高度为来自第一α-β滤波器5的高度①；数据合成处理单元7输出的跟踪浪高为来自第一α-β滤波器5的浪高；数据合成处理单元7输出的MLE跟踪参数为来自MLE误差鉴别器3的鉴别器状态①；数据合成处理单元7输出的波形数据为来自第一FFT功率谱分析器1的波形数据①；数据合成处理单元7输出的AGC值为来自第一α-β滤波器5的AGC①；数据合成处理单元7根据来自第一α-β滤波器5的距离粗调①生成对应的工作时序；数据合成处理单元7输出的信号带宽为来自MLE误差鉴别器3的信号带宽①；数据合成处理单元7输出的OCOG跟踪参数为来自OCOG跟踪器的鉴别器状态②。
步骤3，工作模式的转换：
启始工作的时候，跟踪器工作于模式1。当鉴别器状态②表明OCOG跟踪处于锁定状态且信号带宽②为320MHz时，系统工作模式转换为模式2。转换过后系统的工作带宽仍为320MHz，此时第一α-β滤波器5接收来自第二α-β滤波器6的高度②和AGC②，并以此作为工作状态的初始值；在该模式下此后第一α-β滤波器5不再接收来自第二α-β滤波器6的高度②和AGC②，而是将生成的高度①和AGC①传送给第二α-β滤波器6。
当鉴别器状态①表明MLE跟踪处于失锁状态时，系统工作模式转换为模式1。转换过后OCOG的工作带宽为20MHz，此后第二α-β滤波器6不再接收来自第一α-β滤波器5的高度①和AGC①，相反第二α-β滤波器6将高度②和AGC②输送到第一α-β滤波器5，并以此作为第一α-β滤波器5工作状态的初始值。
Claims (4)
1、一种海陆兼容的卫星雷达高度计跟踪器，包括第一FFT功率谱分析器(1)、MLE误差鉴别器(3)和第一α-β滤波器(5)，其特征在于，还包括第二FFT功率谱分析器(2)、OCOG误差鉴别器(4)、第二α-β滤波器(6)和数据合成处理单元(7)；其中，所述的第一FFT功率谱分析器(1)电连接到所述的MLE误差鉴别器(3)、所述的第一α-β滤波器(5)和所述的数据合成处理单元(7)上，所述的MLE误差鉴别器(3)分别与所述的第一α-β滤波器(5)和所述的数据合成处理单元(7)电连接，所述的第一α-β滤波器(5)与所述的数据合成处理单元(7)电连接；所述的第二FFT功率谱分析器(2)电连接到所述的OCOG误差鉴别器(4)、所述的第二α-β滤波器(6)和所述的数据合成处理单元(7)上，所述的OCOG误差鉴别器(4)分别与所述的第二α-β滤波器(6)和所述的数据合成处理单元(7)电连接，所述的第二α-β滤波器(6)与所述的数据合成处理单元(7)电连接；所述的第一α-β滤波器(5)还与所述的第二α-β滤波器(6)电连接；
所述OCOG误差鉴别器(4)生成第二鉴别器状态和第二信号带宽，当所述第二鉴别器状态表明OCOG处于锁定状态且所述第二信号带宽为320MHz时，所述卫星雷达高度计跟踪器位于第二工作模式，否则位于第一工作模式；在所述第一工作模式下，所述数据合成处理单元(7)生成的第一状态控制信号表示待机状态，生成的第二状态控制信号表示主动工作状态；在所述第二工作模式下，所述数据合成处理单元(7)生成的第一状态控制信号表示工作状态，生成的第二状态控制信号表示被动工作状态；
所述第一α-β滤波器(5)接收到所述数据合成处理单元(7)生成的第一状态控制信号后，如果所述第一状态控制信号是待机状态，此时所述第一α-β滤波器(5)接收来自所述第二α-β滤波器(6)的第二高度信号和第二AGC，并以此值作为工作状态的初始值，同时所述第一α-β滤波器(5)将输出的第一距离精调信号设定为0，并处于待机状态；当所述第一状态控制信号为工作状态时，所述第一α-β滤波器(5)接收所述MLE误差鉴别器(3)的误差信号，将滤波处理后输出的第一高度信号和第一AGC的值送到所述数据合成处理单元(7)和第二α-β滤波器(6)，将滤波处理后输出的浪高值送到所述数据合成处理单元(7)，将滤波处理后输出的第一高度信号分成第一距离粗调和第一距离精调两部分，其中所述第一距离粗调送到所述数据合成处理单元(7)，所述第一距离精调送到所述第一FFT功率谱分析器(1)；
所述第二α-β滤波器(6)接收所述数据合成处理单元(7)生成的第二状态控制信号，当所述第二状态控制信号为被动工作状态时，所述第二α-β滤波器(6)接收来自所述OCOG误差鉴别器(4)输出的误差信号并进行滤波，滤波处理后输出的第二高度信号和第二AGC送到所述数据合成处理单元(7)；滤波处理后输出的第二高度信号分成第二距离粗调和第二距离精调两部分，其中所述第二距离粗调送到所述数据合成处理单元(7)，所述第二距离精调送到所述第二FFT功率谱分析器(2)；完成上述操作后，所述第二α-β滤波器(6)接收来自所述第一α-β滤波器(5)的第一高度信号和第一AGC，并复位所述第二α-β滤波器(6)的当前输出，使得第二高度信号等于第一高度信号，且第二AGC等于第一AGC；
当所述第二状态控制信号为主动工作状态时，所述第二α-β滤波器(6)接收所述OCOG误差鉴别器(4)输出的误差信号并进行滤波，滤波处理后输出的第二高度信号和第二AGC送到所述数据合成处理单元(7)和所述第一α-β滤波器(5)；滤波处理后输出的第二高度信号分成第二距离粗调和第二距离精调两部分，其中所述的第二距离粗调送到所述数据合成处理单元(7)，所述第二距离精调送到所述第二FFT功率谱分析器(2)。
2、根据权利要求1所述的海陆兼容的卫星雷达高度计跟踪器，其特征在于，所述的第二FFT功率谱分析器(2)、所述的OCOG误差鉴别器(4)和所述的第二α-β滤波器(6)在DSP芯片上编程实现。
3、根据权利要求1所述的海陆兼容的卫星雷达高度计跟踪器，其特征在于，所述的数据合成处理单元(7)在DSP芯片和FPGA电路上编程实现。
4、根据权利要求1所述的海陆兼容的卫星雷达高度计跟踪器，其特征在于，所述的第一FFT功率谱分析器(1)、所述的MLE误差鉴别器(3)和所述的第一α-β滤波器(5)在DSP芯片上编程实现。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006101144960A CN100507599C (zh) | 2006-11-10 | 2006-11-10 | 一种海陆兼容的卫星雷达高度计跟踪器 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006101144960A CN100507599C (zh) | 2006-11-10 | 2006-11-10 | 一种海陆兼容的卫星雷达高度计跟踪器 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101178436A CN101178436A (zh) | 2008-05-14 |
| CN100507599C true CN100507599C (zh) | 2009-07-01 |
Family
ID=39404763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2006101144960A Expired - Fee Related CN100507599C (zh) | 2006-11-10 | 2006-11-10 | 一种海陆兼容的卫星雷达高度计跟踪器 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100507599C (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105353369B (zh) * | 2015-10-20 | 2017-06-30 | 中国电子科技集团公司第三十八研究所 | 一种宽雷达波束下高山地测高方法 |
| CN107783131B (zh) * | 2016-08-25 | 2021-05-18 | 大连楼兰科技股份有限公司 | 基于组合波形的自动驾驶汽车防撞雷达系统信号处理方法 |
| CN109839623B (zh) * | 2019-02-14 | 2020-09-11 | 北京遥感设备研究所 | 一种地外天体着陆测量雷达面目标回波信号测距处理方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1989002605A1 (en) * | 1987-09-17 | 1989-03-23 | The Singer Company | Closed loop velocity/altitude sensor for fm-cw doppler radars |
| EP0473082A2 (en) * | 1990-08-28 | 1992-03-04 | Honeywell Inc. | Radar altimeter |
| GB2315628A (en) * | 1986-04-17 | 1998-02-04 | Plessey Co Plc | Radar altimeter |
| CN1301968A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 一种海陆兼容并三维成象的雷达高度计系统及设计方法 |
| WO2002079798A1 (en) * | 2001-03-28 | 2002-10-10 | The Johns Hopkins University | Bistatic delay doppler radar altimeter |
| EP1116046B1 (en) * | 1998-10-02 | 2004-12-01 | Honeywell Inc. | Interferometric synthetic aperture radar altimeter |
-
2006
- 2006-11-10 CN CNB2006101144960A patent/CN100507599C/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2315628A (en) * | 1986-04-17 | 1998-02-04 | Plessey Co Plc | Radar altimeter |
| WO1989002605A1 (en) * | 1987-09-17 | 1989-03-23 | The Singer Company | Closed loop velocity/altitude sensor for fm-cw doppler radars |
| EP0473082A2 (en) * | 1990-08-28 | 1992-03-04 | Honeywell Inc. | Radar altimeter |
| EP1116046B1 (en) * | 1998-10-02 | 2004-12-01 | Honeywell Inc. | Interferometric synthetic aperture radar altimeter |
| CN1301968A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 一种海陆兼容并三维成象的雷达高度计系统及设计方法 |
| WO2002079798A1 (en) * | 2001-03-28 | 2002-10-10 | The Johns Hopkins University | Bistatic delay doppler radar altimeter |
Non-Patent Citations (2)
| Title |
|---|
| 基于准最大似然估计的星载海洋雷达高度计跟踪器的实现及机载试验验证. 许可，刘和光.电子学报，第32卷第6期. 2004 * |
| 星载海洋雷达高度计中基于FPGA的FFT运算模块的实现及性能分析. 唐昌平，唐月英，许可，刘和光.遥感技术与应用，第19卷第3期. 2004 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101178436A (zh) | 2008-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TWI400469B (zh) | 校準複數雷達的方法及裝置 | |
| WON et al. | Performance comparison of different forms of Kalman filter approaches for a vector‐based GNSS signal tracking loop | |
| AU2012205602B2 (en) | Method and system for determining clock corrections | |
| US7936305B2 (en) | System and method for determining attitude using spatial shift key (SSK) modulation signatures | |
| US20090043511A1 (en) | Assisted Seismic Measuring System Including GPS Receivers | |
| CN102033236B (zh) | 一种卫星导航位置速度联合估计方法 | |
| JPH0688865A (ja) | 単一受信機使用ｇｐｓ指向又は姿勢システム | |
| EP1538455A1 (en) | A method of calibrating an adaptive antenna array of a satellite navigation system | |
| CN107656300B (zh) | 基于北斗/gps双模软件接收机的卫星/惯性超紧组合方法 | |
| CA2407943A1 (en) | Low cost system and method for making dual band gps measurements | |
| CA2746187A1 (en) | Synthetic aperture antenna device for receiving signals of a system comprising a carrier and means for determining its trajectory | |
| AU2005226022A1 (en) | Method for back-up dual-frequency navigation during brief periods when measurement data is unavailable on one of two frequencies | |
| CN112639530B (zh) | 包括定位信号的信号检测的方法和系统 | |
| US6128557A (en) | Method and apparatus using GPS to determine position and attitude of a rotating vehicle | |
| CN100507599C (zh) | 一种海陆兼容的卫星雷达高度计跟踪器 | |
| CN110018506A (zh) | 基于和差组合卡尔曼滤波器的gnss双频联合跟踪算法 | |
| CN110398764A (zh) | 一种北斗导航系统双频段融合定位方法及装置 | |
| CN112946711A (zh) | 一种gnss/ins组合导航系统的导航方法 | |
| US20130287067A1 (en) | Gps aided open loop coherent timing | |
| CN201021939Y (zh) | 一种卫星雷达高度计跟踪器 | |
| Sousa et al. | Performance analysis of a VDFLL GNSS receiver architecture under ionospheric scintillation and multipath conditions | |
| Hummon et al. | A direct comparison of two RDI shipboard ADCPs: A 75-kHz ocean surveyor and a 150-kHz narrow band | |
| JP5163474B2 (ja) | Ｇｎｓｓ受信装置及び測位方法 | |
| CN108107454B (zh) | 一种惯性信息辅助卫星深组合环路 | |
| US6297769B1 (en) | System and method to estimate carrier signal in global positioning systems (GPS) |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090701 Termination date: 20161110 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |