CN112114306B - 一种提高探测器测量精度的方法和装置 - Google Patents
一种提高探测器测量精度的方法和装置 Download PDFInfo
- Publication number
- CN112114306B CN112114306B CN201910534530.7A CN201910534530A CN112114306B CN 112114306 B CN112114306 B CN 112114306B CN 201910534530 A CN201910534530 A CN 201910534530A CN 112114306 B CN112114306 B CN 112114306B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- doppler
- module
- analysis
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims abstract description 72
- 238000001514 detection method Methods 0.000 claims abstract description 39
- 238000005259 measurement Methods 0.000 claims abstract description 38
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims description 43
- 230000003321 amplification Effects 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 3
- 239000003550 marker Substances 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000001131 transforming effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
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
- G01S13/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
-
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种提高探测器测量精度的方法和装置,方法包括:S1,将多普勒信号转换为解析信号,并将解析信号分成第一解析信号和第二解析信号;S2,将第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;S3,实现探测信号的复现,生成复现信号;S4,补偿第二解析信号,以使第二解析信号与复现信号保持时间同步;S5,第二解析信号进行混频,输出零中频信号;S6,将零中频信号和复现信号进行共轭相乘,得到相位差信号;S7,滤除相位差信号的高频噪声,输出低频相位差信号;S8,将多普勒频率初步估计值与低频相位差信号保持时间同步;S9,变换得到复信号,根据复信号获得多普勒信号频率的精确估计值,实现探测器的精准测量。
Description
技术领域
本发明涉及深空探测技术领域,尤其涉及一种提高探测器测量精度的方法和装置。
背景技术
探测器的无线电通信主要包括上行链路和下行链路,主要用于探测器的测控和载荷数据的下传,同时对下行链路信号的测量还可以应用到无线电科学的研究中,通过无线电测距和多普勒频移测速可以获得探测器的位置和速度信号,是实现航天器定轨的一种常用手段,利用下行信号测量到的无线电波的幅度、多普勒频率和极化特性等可以进行行星科学的研究。根据上行链路和下行链路的形式可以将台站对探测器的跟踪划分为单程、双程和三程模式,根据链路是否闭合分为闭环和开环模式。为满足高精度定轨和行星无线电科学研究的需要,对多普勒频率的测量精度提出了很高的要求,对于下行信号,一般情况下数传信号的发射功率要远大于测控信标的功率,所以理论上采用数传信号可以获得更高的多普勒测量精度,对于双程模式,由于下行信号是完全已知的,可以通过相干积分的方式实现载波频率的测量,而相对于单程和三程模式,由于数传下行信号是未知的,一般主要采用科斯塔斯环实现,而由于科斯塔斯环中噪声的相干特性,由于其平方运算的非线性、噪声自相乘以及噪声信号交叉项导致信噪比恶化,尤其在低信噪比条件下,将导致噪声方差急剧增大,从而造成多普勒估计精度的严重下降。
发明内容
(一)要解决的技术问题
基于上述问题,本发明提供了一种提高探测器测量精度的方法和装置,对于基于多普勒频移的探测器的多普勒信号进行解析、复现、滤波等操作实现多普勒信号的精准估计,进而提高探测器的测量精度。
(二)技术方案
第一方面,本发明提供了一种提高探测器测量精度的方法,探测器基于多普勒频移进行测量,探测器发送探测信号,探测信号经发射后由地面天线接收,由于探测器与所述地面天线存在相对速度,使得探测信号产生多普勒频移,得到多普勒信号,方法包括:S1,将多普勒信号转换为解析信号,并将解析信号分成两路,分别为第一解析信号和第二解析信号;S2,将第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;S3,根据解调数据和探测信号的模式实现探测信号的复现,生成复现信号;S4,根据解调和信号复现的信号延迟补偿第二解析信号,以使第二解析信号与复现信号保持时间同步;S5,对步骤S4得到的第二解析信号进行混频,输出零中频信号;S6,将零中频信号和复现信号进行共轭相乘,得到相位差信号;S7,滤除相位差信号的高频噪声,输出低频相位差信号;S8,将多普勒频率初步估计值与低频相位差信号保持时间同步;S9,利用低频相位差信号和所述步骤S2输出的多普勒频率初步估计值进行变换得到复信号,根据复信号获得探测信号的载波相位和多普勒信号频率的精确估计值,进而实现探测器的精准测量。
可选的,步骤S5具体为:对步骤S4得到的第二解析信号与步骤S9得到的复信号进行混频,输出零中频信号。
可选的,步骤S1包括:S11,将多普勒信号转换为数字信号;S12,将数字信号转换为解析信号。
可选的,步骤S11之前还包括:S0,将对多普勒信号进行放大,并将多普勒信号变频为中频多普勒信号。
可选的,还包括预设频标信号,所述步骤S11具体为:根据频标信号将多普勒信号转换为数字信号。
可选的,步骤S7具体为:对信号相位差信号进行积分,以提高信号相位差信号的的信噪比,滤除相位差信号的高频噪声,输出低频相位差信号。
第二方面,本发明提供了一种提高探测器测量精度的装置,探测器基于多普勒频移进行测量,探测器发射探测信号,探测信号经反射后得到多普勒信号,装置包括:转换模块,用于将多普勒信号转换为解析信号,并将解析信号分成两路,分别为第一解析信号和第二解析信号;解调模块,用于将所述第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;信号复现模块,用于根据解调数据和探测信号的模式实现探测信号的复现,生成复现信号;第一时间延迟模块,用于根据解调和信号复现的信号延迟补偿所述第二解析信号,以使所述第二解析信号与复现信号保持时间同步;混频模块,用于对所述第一时间延迟模块得到的第二解析信号进行混频,输出零中频信号;共轭相乘模块,用于将零中频信号和复现信号进行共轭相乘,得到相位差信号;环路滤波模块,用于滤除相位差信号的高频噪声,输出低频相位差信号;第二时间延迟模块,用于将所述多普勒频率初步估计值与低频相位差信号保持时间同步;数字控制振荡器模块,用于利用低频相位差信号和所述第二时间延迟模块输出的多普勒频率初步估计值进行变换得到复信号,根据所述复信号获得探测信号的载波相位和多普勒信号频率的精确估计值,进而实现探测器的精准测量。
可选的,转换模块包括模数转换子模块和解析信号转换子模块,其中:模数转换子模块,用于将多普勒信号转换为数字信号;解析信号转换子模块,用于将数字信号转换为解析信号。
可选的,还包括低噪声放大模块和下变频模块,多普勒信号在输入转换模块前输入所述低噪声放大模块进行放大,并输入下变频模块变频为中频多普勒信号。
可选的,还包括原子钟模块和频率综合器模块,其中:原子钟模块,用于提供10MHz频率标准和1pps信号;频率综合器模块,用于输出与原子钟模块提供10MHz频率标准相同频率的预设频标信号,以使下变频模块和模数转换子模块以预设频标信号作为频率标准。
(三)有益效果
本发明提供了一种提高探测器测量精度的方法和装置,通过对接收的多普勒频移信号进行解析、复现、滤波等操作,可以克服传统方法中由于科斯塔斯环中噪声的相干特性,平方效应会导致的多普勒估计精度的下降的问题。本方法将噪声相干转化为非相干相关,可以显著提高多普勒估计精度,从而显著提高探测器轨道的估计精度和行星无线电科学的研究应用效果,尤其是针对深空探测低信噪比条件下效果更加显著。
附图说明
图1示意性示出了本公开实施例的提高探测器测量精度的方法的步骤图;
图2示意性示出了本公开实施例的提高探测器测量精度的装置的工作流程图;
图3示意性示出了本公开实施例的图2所示提高探测器测量精度的装置中解调模块的工作流程图;
图4示意性示出了本公开实施例的提高探测器测量精度的方法与传统科斯塔斯环方法得到测量精度的对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
第一方面,本发明提供了一种提高探测器测量精度的方法,探测器基于多普勒频移进行测量,探测器发送探测信号,探测信号经发射传输后由地面天线接收,由于探测器与接收天线的存在相对速度,使信号产生多普勒频移,参见图1,方法包括:S1,将多普勒信号转换为解析信号,并将解析信号分成两路,分别为第一解析信号和第二解析信号;S2,将第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;S3,根据解调数据和探测信号的模式实现探测信号的复现,生成复现信号;S4,根据解调和信号复现的信号延迟补偿第二解析信号,以使第二解析信号与复现信号保持时间同步;S5,对步骤S4得到的第二解析信号进行混频,输出零中频信号;S6,将零中频信号和复现信号进行共轭相乘,得到相位差信号;S7,滤除相位差信号的高频噪声,输出低频相位差信号;S8,将多普勒频率初步估计值与低频相位差信号保持时间同步;S9,利用低频相位差信号和步骤S2输出的多普勒频率初步估计值进行变换得到复信号,根据所述复信号获得多普勒信号频率的精确估计值,进而实现探测器的精准测量。以下将以一具体实施例为例对该方法进行详细介绍。
S1,将多普勒信号转换为解析信号,并将解析信号分成两路,分别为第一解析信号和第二解析信号;
具体的,接收到多普勒信号后,首先将该多普勒信号输入低噪声放大器等中进行放大,并将放大后的信号输入下变频器中变频为中频多普勒信号,然后执行如下步骤S11和S22,具体的:
S11,将多普勒信号转换为数字信号;
具体的,将上述中频多普勒信号经模数转换模块(ADC)进行数模转换,将该中频多普勒信号转换为数字信号。
S12,将数字信号转换为解析信号。
将数字采样信号经解析信号转换模块等转换为解析复信号,具体实现时可以采用希尔伯变换或正交采样方式进行,本发明实施例中,若数字信号为x(t),则经转换后的解析信号z(t)为:
其中,为x(t)的希尔伯特变换,t为时间序列;
对于窄带信号,解析信号z(t)可表示为:
其中,fc为下变频后的载波频率;φ为初始相位;j为虚数单位;a(t)为调制信号,其中,a(nTc)为二进制序列,n为序号,Tc为码片宽度,q(t)为矩形窗函数,对于BPSK信号
a(nTc)=±1
并将生成的解析信号分成两路相同的信号,即第一解析信号和第二解析信号。
S2,将第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;
将上述步骤S1输出的第一解析信号进行解调得到该第一解析信号携带的信息,输出解调数据和多普勒频率初步估计值。
本发明实施例中,对上述z(t)信号进行解调,获得解调后的二进制序列a(nTc)。
S3,根据解调数据和探测信号的模式实现探测信号的复现,生成复现信号;
根据上述解调数据和探测器所发送的探测信号的模式对探测器所发送的探测信号进行复现,也即重新生成连续信号。
本发明实施例中,根据输入的解调后的序列a(nTc)实现探测发送的下行数传信号也即探测信号的复现,复现信号为:
其中,fp为探测器理论发射频率减去下变频参考频率后的中频频率,且不受多普勒频率的影响;φp为探测信号的初始相位;τp为由于解调和信号复现带来的时间延迟。
S4,根据解调和信号复现的信号延迟补偿第二解析信号,以使第二解析信号与复现信号保持时间同步;
在对第二解析信号进行处理时需要参考上述步骤S2~S3中的解调和信号复现的信号延迟,并根据解调和信号复现的信号延迟补偿第二解析信号以使第二解析信号与复现信号保持时间同步。
本发明实施例中,经补偿后的第二解析信号可表示为:
其中,φ为初始相位。
S5,对步骤S4得到的第二解析信号进行混频,输出零中频信号;
本方法为闭环控制,将步骤S4得到的第二解析信号以及以下步骤S9得到的复信号进行混频,输出零中频信号。
在本发明实施例中,若步骤S9得到的输出的复信号为:
其中,fd为测量的多普勒频率的精确估计值,φd为探测信号中载波相位测量值。则经混频后的零中频信号z2(t)为:
其中,为zd(t)的共轭。
S6,将零中频信号和复现信号进行共轭相乘,得到相位差信号;
将步骤S6得到的零中频信号和步骤S3得到的复现信号进行共轭相乘,得到相位差信号。
本发明实施例中,混频后的信号z2(t)经共轭相乘后输出的相位差信号为:
其中,Δf=fc-fd-fp,Δφ=φ-φd-φp-2πfcτp+2πfpτp。
由上可以看出,当步骤S9输出的信号频率fd和相位φd满足如下条件fd=fc-fp和φd=φ-φp-2πfcτp+2πfpτp时,可得到γ(t)=a2(t-τp),从而实现相位差信号γ(t)的相位为0,实现多普勒频率侧理想估计。但实际情况下,由于噪声的影响,γ(t)需通过积分和滤波等操作,去除相位差信号中的高频分量和噪声,提高信噪比。
S7,滤除相位差信号的高频噪声,输出低频相位差信号;
对上述相位差信号进行积分,提高信噪比,滤除相位差信号的高频噪声,输出低频相位差信号。此步骤可采用通用的环路滤波器实现。
S8,将多普勒频率初步估计值与低频相位差信号保持时间同步;
将上述步骤S2得到的多普勒频率初步估计值进行时间延迟以使的多普勒频率初步估计值与步骤S7得到的低频相位差信号保持时间同步。
S9,利用低频相位差信号和步骤S2模块输出的多普勒频率初步估计值进行变换得到复信号,根据所述复信号获得所述探测信号的载波相位和多普勒信号频率的精确估计值,进而实现探测器的精准测量。具体实现方法如下:
将低频相位差信号的相位累加至步骤S9的输出复信号,即:
其中,zd(tm)为tm时刻的步骤S9输出给混频器的复信号;zd(tm)的其初始值zd(t0)采用步骤S2模块输出的多普勒频率初步估计值;φγ(tm)为γ(t)经环路滤波后在tm时刻的相位;m为采样点序号,且m=0,1,2,…;tm为时间序列。
若信号zd(tm)的相位为则多普勒频率和载波的初始相位的最小二乘估计为:
式中,为带辨识的参数向量,A为系数矩阵,b为载波相位差观测向量,且
式中,fd为测量的多普勒频率的精确估计值,φd为探测信号中载波相位测量值。上式只给出了多普勒频率的一次线性估计模型,根据多普勒动态变化程度,实际应用时也可以采用二次或更高次的估计模型。
另一方面,本发明提供了一种提高探测器测量精度的装置,探测器基于多普勒频移进行测量,探测器发射探测信号,探测信号经反射后得到多普勒信号,如图2所示,装置包括:
转换模块,用于将多普勒信号转换为解析信号,并将解析信号分成两路,分别为第一解析信号和第二解析信号;
转换模块包括模数转换子模块(ADC)和解析信号转换子模块,其中:
模数转换子模块,用于将多普勒信号转换为数字信号;
解析信号转换子模块,用于将数字信号转换为解析信号。
本发明实施例中,解析信号转换子模块实现数字信号向解析信号的转换,具体可采用希尔伯变换或正交采样方式进行,本发明实施例中,若经模数转换子模块(ADC)后的数字信号为x(t),则经解析信号转换子模块转换后的解析信号为:
其中,为x(t)的希尔伯特变换,t为时间序列;
对于窄带信号,解析信号可表示为:
其中,fc为下变频后的载波频率;a(t)为调制信号,其中,a(nTc)为二进制序列,n为序号,Tc为码片宽度,q(t)为矩形窗函数,对于BPSK信号
a(nTc)=±1
并将生成的解析信号分成两路相同的信号,即第一解析信号和第二解析信号。
解调模块,用于将第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;
解调模块对上述z(t)信号进行解调,获得解调后的二进制序列a(nTc)。如图3所示,解调模块由匹配滤波、定时同步模块、频率捕获模块、载波同步模块和符号判别等组成。图2是本发明解调模块数据处理流程示意图,信号经匹配滤波提升信噪比后,分别将输出信号送入频率捕获和定时同步模块完成多普勒频偏估计和符号同步;频率捕获模块获得多普勒频率初始值,作为载波同步模块的搜索初始值;载波同步模块采用基于科斯塔斯环的载波锁定模式,实现载波相位的锁定;载波同步模块跟踪获得的多普勒频移值在经实践延迟后输出给数字控制振荡器模块(NCO)模块作为多普勒频率的初值,输出的同相信号经符号判别模块获得解调后的二进制序列a(nTc)。
信号复现模块,用于根据解调数据和探测信号的模式实现探测信号的复现,生成复现信号;
信号复现模块根据输入的解调后的序列a(nTc)实现探测发送的下行数传信号也即探测信号的复现,复现信号为:
其中,fp为探测器理论发射频率减去下变频参考频率后的中频频率,且不受多普勒频率的影响;φp为发射信号的初始相位;τp为由于解调和信号复现带来的时间延迟。
第一时间延迟模块,用于根据解调和信号复现的信号延迟补偿第二解析信号,以使第二解析信号与所述复现信号保持时间同步;
本发明实施例中,经第一时间延迟模块补偿后的第二解析信号可表示为:
其中,φ为初始相位。
混频模块,用于对第一时间延迟模块得到的第二解析信号进行混频,输出零中频信号;
若数字控制振荡器模块(NCO)输出给混频模块的复信号为:
其中,fd为测量的多普勒频率的精确估计值,φd为探测信号中载波相位测量值。则经混频模块混频后的零中频信号为:
共轭相乘模块,用于将所述零中频信号和复现信号进行共轭相乘,得到相位差信号;
本发明实施例中,混频模块混频后的信号z2(t)经共轭相乘后输出的相位差信号为:
其中,Δf=fc-fd-fp,Δφ=φ-φd-φp-2πfcτp+2πfpτp。
由上可以看出,当NCO输出的信号频率fd和相位φd满足如下条件fd=fc-fp和φd=φ-φp-2πfcτp+2πfpτp时,可得到γ(t)=a2(t-τp),从而实现相位差信号γ(t)的相位为0,实现多普勒频率的理想估计。但实际情况下,由于噪声的影响,γ(t)需通过积分和滤波等操作,去除相位差信号中的高频分量和噪声,提高信噪比。
环路滤波模块,用于滤除相位差信号的高频噪声,输出低频相位差信号;
环路滤波器可采用锁相环通用的线性环路滤波器,其输出信号用于实现对NCO输出信号的闭环控制。
第二时间延迟模块,用于将多普勒频率初步估计值与低频相位差信号保持时间同步;
数字控制振荡器模块,用于利用低频相位差信号和第二时间延迟模块输出的多普勒频率初步估计值变换得到复信号,根据所述复信号获得多普勒信号频率的精确估计值,进而实现探测器的精准测量。
将低频相位差信号的相位累加至NCO模块输出的复信号,即:
其中,zd(tm)为tm时刻的NCO输出给混频器的复信号;zd(tm)的初始值zd(t0)采用多普勒频率初步估计值;φγ(tm)为γ(t)经环路滤波后在tm时刻的相位;m为采样点序号,且m=0,1,2,…;tm为时间序列。
若信号zd(tm)的相位为则多普勒频率和载波的初始相位的最小二乘估计为:
式中,为带辨识的参数向量,A为系数矩阵,b为载波相位差观测向量,且
式中,fd为测量的多普勒频率的精确估计值,φd为探测信号中载波相位测量值。上式只给出了多普勒频率的一次线性估计模型,根据多普勒动态变化程度,实际应用时也可以采用二次或更高次的估计模型。
另该装置还包括噪声放大模块、下变频模块、原子钟模块以及频率综合器模块,其中,多普勒信号在输入转换模块前输入低噪声放大模块进行放大,并输入下变频模块变频为中频多普勒信号;原子钟模块,用于提供10MHz频率标准和1pps信号;频率综合器模块,用于输出与原子钟模块提供10MHz频率标准相同频率的预设频标信号,以使下变频模块和模数转换子模块以预设频标信号作为频率标准。
图4是传统方法与本申请中所公开方法多普勒估计信噪比对比示意图。
传统的基于科斯塔斯环的多普勒频率测量误差可表示为
其中,SNRi为输入信号信噪比,Bi为输入信号带宽,BL为锁相环噪声带宽,λ为信号波长,Tc为码位宽度。而本发明提出的多普勒测量方法,由于实现了信号的相干积分,其多普勒频率测量误差可表示为:
针对BPSK信号,在码速率为1Mbps时,基于科斯塔斯环和新方法的多普勒测量误差,如图4所示。从中可以看出,在同等信噪比条件下,新方法明显低于科斯塔斯环的测量误差。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种提高探测器测量精度的方法,所述探测器基于多普勒频移进行测量,所述探测器发送探测信号,所述探测信号经发射后由地面天线接收,由于所述探测器与所述地面天线存在相对速度,使得所述探测信号产生多普勒频移,得到多普勒信号,所述方法包括:
S1,将所述多普勒信号转换为解析信号,并将所述解析信号分成两路,分别为第一解析信号和第二解析信号;
S2,将所述第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;
S3,根据所述解调数据和探测信号的模式实现所述探测信号的复现,生成复现信号;
S4,根据解调和信号复现的信号延迟补偿所述第二解析信号,以使所述第二解析信号与所述复现信号保持时间同步;
S5,对所述步骤S4得到的第二解析信号进行混频,输出零中频信号;
S6,将所述零中频信号和复现信号进行共轭相乘,得到相位差信号;
S7,滤除所述相位差信号的高频噪声,输出低频相位差信号;
S8,将所述多普勒频率初步估计值与所述低频相位差信号保持时间同步;
S9,利用所述低频相位差信号和所述步骤S2输出的多普勒频率初步估计值进行变换得到复信号,根据所述复信号获得所述探测信号的载波相位和多普勒信号频率的精确估计值,进而实现所述探测器的精准测量。
2.根据权利要求1所述的方法,所述步骤S5具体为:
对所述步骤S4得到的第二解析信号与步骤S9得到的复信号进行混频,输出零中频信号。
3.根据权利要求1所述的方法,所述步骤S1包括:
S11,将所述多普勒信号转换为数字信号;
S12,将所述数字信号转换为解析信号。
4.根据权利要求3所述的方法,所述步骤S11之前还包括:
S0,将对所述多普勒信号进行放大,并将所述多普勒信号变频为中频多普勒信号。
5.根据权利要求3所述的方法,还包括预设频标信号,所述步骤S11具体为:
根据所述频标信号将所述多普勒信号转换为数字信号。
6.根据权利要求1所述的方法,所述步骤S7具体为:
对所述相位差信号进行积分,以提高所述相位差信号的信噪比,滤除所述相位差信号的高频噪声,输出低频相位差信号。
7.一种提高探测器测量精度的装置,所述探测器基于多普勒频移进行测量,所述探测器发射探测信号,所述探测信号经反射后得到多普勒信号,所述装置包括:
转换模块,用于将所述多普勒信号转换为解析信号,并将所述解析信号分成两路,分别为第一解析信号和第二解析信号;
解调模块,用于将所述第一解析信号进行解调,得到解调数据和多普勒频率初步估计值;
信号复现模块,用于根据所述解调数据和探测信号的模式实现所述探测信号的复现,生成复现信号;
第一时间延迟模块,用于根据解调和信号复现的信号延迟补偿所述第二解析信号,以使所述第二解析信号与所述复现信号保持时间同步;
混频模块,用于对所述第一时间延迟模块得到的第二解析信号进行混频,输出零中频信号;
共轭相乘模块,用于将所述零中频信号和复现信号进行共轭相乘,得到相位差信号;
环路滤波模块,用于滤除所述相位差信号的高频噪声,输出低频相位差信号;
第二时间延迟模块,用于将所述多普勒频率初步估计值与所述低频相位差信号保持时间同步;
数字控制振荡器模块,用于利用所述低频相位差信号和所述第二时间延迟模块输出的多普勒频率初步估计值进行变换得到复信号,根据所述复信号获得所述探测信号的载波相位和多普勒信号频率的精确估计值,进而实现所述探测器的精准测量。
8.根据权利要求7所述的装置,所述转换模块包括模数转换子模块和解析信号转换子模块,其中:
模数转换子模块,用于将所述多普勒信号转换为数字信号;
解析信号转换子模块,用于将所述数字信号转换为解析信号。
9.根据权利要求8所述的装置,还包括低噪声放大模块和下变频模块,所述多普勒信号在输入所述转换模块前输入所述低噪声放大模块进行放大,并输入下变频模块变频为中频多普勒信号。
10.根据权利要求9所述的装置,还包括原子钟模块和频率综合器模块,其中:
原子钟模块,用于提供10MHz频率标准和1pps信号;
频率综合器模块,用于输出与所述原子钟模块提供10MHz频率标准相同频率的预设频标信号,以使所述下变频模块和模数转换子模块以所述预设频标信号作为频率标准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910534530.7A CN112114306B (zh) | 2019-06-19 | 2019-06-19 | 一种提高探测器测量精度的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910534530.7A CN112114306B (zh) | 2019-06-19 | 2019-06-19 | 一种提高探测器测量精度的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112114306A CN112114306A (zh) | 2020-12-22 |
CN112114306B true CN112114306B (zh) | 2023-08-18 |
Family
ID=73795128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910534530.7A Active CN112114306B (zh) | 2019-06-19 | 2019-06-19 | 一种提高探测器测量精度的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112114306B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1986005279A1 (en) * | 1985-03-01 | 1986-09-12 | Gene Wylie Adams | Imaging doppler interferometer |
US4743110A (en) * | 1983-03-29 | 1988-05-10 | Thomson Csf | Laser telemetry and Doppler measurement system with pulse compression |
DE4406404A1 (de) * | 1994-02-26 | 1995-08-31 | Daimler Benz Ag | Verfahren und Anordnung zur Bestimmung der Lage und Bewegung eines Objekts |
CN104237877A (zh) * | 2014-09-19 | 2014-12-24 | 陕西长岭电子科技有限责任公司 | 机载型自主式测速测高雷达系统及测速测高方法 |
CN105548968A (zh) * | 2015-12-14 | 2016-05-04 | 桂林电子科技大学 | 一种多带Chirp在动目标检测时抑制多普勒扩展的方法 |
CN109633574A (zh) * | 2018-10-25 | 2019-04-16 | 上海无线电设备研究所 | 一种用于深空探测的宽范围高精度多普勒测量方法 |
-
2019
- 2019-06-19 CN CN201910534530.7A patent/CN112114306B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4743110A (en) * | 1983-03-29 | 1988-05-10 | Thomson Csf | Laser telemetry and Doppler measurement system with pulse compression |
WO1986005279A1 (en) * | 1985-03-01 | 1986-09-12 | Gene Wylie Adams | Imaging doppler interferometer |
DE4406404A1 (de) * | 1994-02-26 | 1995-08-31 | Daimler Benz Ag | Verfahren und Anordnung zur Bestimmung der Lage und Bewegung eines Objekts |
CN104237877A (zh) * | 2014-09-19 | 2014-12-24 | 陕西长岭电子科技有限责任公司 | 机载型自主式测速测高雷达系统及测速测高方法 |
CN105548968A (zh) * | 2015-12-14 | 2016-05-04 | 桂林电子科技大学 | 一种多带Chirp在动目标检测时抑制多普勒扩展的方法 |
CN109633574A (zh) * | 2018-10-25 | 2019-04-16 | 上海无线电设备研究所 | 一种用于深空探测的宽范围高精度多普勒测量方法 |
Non-Patent Citations (1)
Title |
---|
舰载二维超低空搜索雷达总体设计;张显满;舰船科学技术(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112114306A (zh) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113238261A (zh) | 低轨卫星扩频通信体制信号捕获跟踪系统 | |
US7146146B2 (en) | Systems and methods for coherent adaptive calibration in a receiver | |
US20140099907A1 (en) | Rugged iq receiver based rf gain measurements | |
US20070230612A1 (en) | Apparatus, system, and method for amplitude-phase synchronization in polar transmitter | |
CN102571137A (zh) | 一种全数字直接序列扩频通信系统及其伪码快速捕获方法 | |
JP7478674B2 (ja) | 無線受信器 | |
JPH0821961B2 (ja) | デイジタル復調装置及び差動位相偏移キーイング復調装置ならびに低信号対雑音比入力信号復調方法 | |
US12047105B2 (en) | Apparatus and methods for DC-offset estimation | |
JP2661534B2 (ja) | スペクトラム拡散受信方法及び受信機 | |
CN108345014B (zh) | 一种正交复用boc调制信号的接收方法 | |
CN112578422A (zh) | 复合导航信号接收方法及接收机 | |
CN102638319A (zh) | 一种导航卫星二进制偏移载波信号的调制性能测试方法 | |
CN104765052A (zh) | 一种geo导航卫星高灵敏度载波跟踪方法 | |
EP0757461B1 (en) | Method of estimating signal quality for a direct sequence spread spectrum receiver | |
US8144752B2 (en) | Method and device for receiving a BOC modulation radio-navigation signal | |
EP3182606B1 (en) | Blanking using signal-based thresholding schemes | |
US20250088404A1 (en) | Electronic device and operating method to compensate for in-phase/quadrature imbalance | |
CN100417148C (zh) | 高灵敏度接收机及改进接收机灵敏度的方法 | |
CN112114306B (zh) | 一种提高探测器测量精度的方法和装置 | |
US4513429A (en) | Sample data phase locked loop for adaptive doppler tracker | |
JP2001230702A (ja) | スペクトラム拡散通信同期捕捉回路 | |
US5343210A (en) | Satellite radionavigation receiver | |
CN101576612B (zh) | 估计gps信号载噪比的方法和gps接收机 | |
US20130257652A1 (en) | Apparatus and method | |
US9419834B2 (en) | MPSK demodulation apparatus and method |
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 |