CN113556188A - 一种测控天线组阵精确频率偏差估计与补偿装置 - Google Patents
一种测控天线组阵精确频率偏差估计与补偿装置 Download PDFInfo
- Publication number
- CN113556188A CN113556188A CN202110839557.4A CN202110839557A CN113556188A CN 113556188 A CN113556188 A CN 113556188A CN 202110839557 A CN202110839557 A CN 202110839557A CN 113556188 A CN113556188 A CN 113556188A
- Authority
- CN
- China
- Prior art keywords
- unit
- residual
- doppler
- time delay
- cross
- 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.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种测控天线组阵精确频率偏差估计与补偿装置,适用于接收天线组阵中深空信号的高效合成、实时合成等相关应用领域。该装置对多天线接收信号进行下变频、滤波、积分清洗等预处理,对预处理后的信号进行FFT,利用频域互谱信号提取多天线的相位差,并根据相位差的变化情况拟合提取精确的残余多普勒值。本发明适用于低信噪比情况下多天线接收信号之间的时频相差对齐,具有频率偏差估计和补偿精度高等优点。
Description
技术领域
本发明涉及到深空测控通信中的接收天线组阵信号合成技术、分布式雷达收发相参合成技术、干涉测量数据处理技术等,特别适用于接收天线组阵中深空信号的高效合成、实时合成等相关应用领域。
背景技术
在深空测控通信系统中,由于深空航天器目标与地球测站相距遥远,从几千万到上亿公里,致使现有测控通信设备的接收信号强度非常微弱。接收天线组阵是实现提高接收信号SNR的有效手段,但实现天线组阵合成首先需要将各天线接收信号进行时延偏差、相位偏差、频率偏差修正对齐。传统的频率偏差采用时域载波跟踪或频域FFT技术,只适用于高信噪比情况,且估计精度差。
发明内容
本发明所要解决的技术问题在于本发明公开了测控天线组阵自动频率偏差估计与补偿装置,该装置通过多天线接收信号的频域FX相关计算,基于自适应FX相关器能够实现低信噪比情况下的天线阵间时延和相位偏差高精度的自动估计,并基于相位偏差的估计值进行多项式拟合得到频率偏差的估计值。
本发明采用的技术方案为:
一种测控天线组阵精确频率偏差估计与补偿装置,包括频偏模型建立单元、正交下变频和积分清洗单元、FFT单元、时延相位补偿单元、频域合成单元、互谱计算单元和精确残余多普勒估计单元;
频偏模型建立单元:用于根据外部输入目标轨位信息和各天线的位置坐标信息,建立天线1~天线N多普勒多项式预报值fdk(t),将多普勒多项式预报值送入正交下变频和积分清洗单元;其中,N为大于1的自然数,fdk(t)=a0+a1·t+a2·t2+…,t为时间,k=1,2,…,N为天线标识,a0、a1、a2…分别为多项式系数;
正交下变频和积分清洗单元:用于根据外部输入的中频载波标称值、频偏模型建立单元输入的多普勒多项式预报值和精确残余多普勒估计单元输入的精确残余多普勒估计值构造下变频本地信号slocal=exp(-j2π(f0+fdk(t)+Δfk0(t))t),利用构造的下变频本地信号对输入信号进行正交下变频,得到复数基带信号,对复数基带信号的I、Q支路分别进行积分清洗降低数据速率后,将积分清洗后的I、Q基带信号送入FFT单元;其中,f0为中频载波标称值,fdk(t)为多普勒多项式预报值,Δfk0(t)为精确残余多普勒估计值;
FFT单元:用于对积分清洗后的I、Q基带信号进行复数FFT运算,并将计算结果送入时延相位补偿单元;
时延相位补偿单元:用于利用互谱计算单元计算得到的残余时延差和残余相位差对FFT单元的计算结果进行时延相位精补偿,将精补偿后的结果送入频域合成单元和互谱计算单元;
频域合成单元:将精补偿后的信号进行幅度加权后频域合成,将频域合成信号进行IFFT后输出,并将频域合成信号送互谱计算单元;
互谱计算单元:根据精补偿后的信号和频域合成信号进行互谱计算,将互谱结果的互谱相位信息通过最小二乘线性拟合进行计算,得到残余时延差和残余相位差;将残余时延差和残余相位差送入时延相位补偿单元,同时将残余相位差送入精确残余多普勒估计单元,其中,第k路信号的残余相位差为:θk0为固定的相位差值,Δfk0(t)为精确残余多普勒频差;
精确残余多普勒估计单元:根据互谱计算单元得到的残余相位差,通过多项式拟合得到精确残余多普勒估计值,并将得到的精确残余多普勒估计值Δfk0(t)更新到正交下变频和积分清洗单元;其中,精确残余多普勒估计值Δfk0(t)=b0+b1·t+b2·t2+…,b0、b1、b2分别为多项式拟合系数。
本发明与背景技术相比具有如下优点:
本发明的测控天线组阵精确频率偏差估计与补偿的装置,基于自适应FX相关器能够实现天线阵间时延和相位偏差高精度的自动估计,并基于相位偏差的估计值进行多项式拟合得到频率偏差的估计值,具有估计和补偿精度高、适应信噪比低等优点,可以实现多测控天线低信噪比情况下的高精度对齐合成,提高合成信号信噪比。
附图说明
图1是本发明的原理框图。
具体实施方式
参照图1,一种测控天线组阵精确频率偏差估计与补偿装置,包括频偏模型建立单元、正交下变频和积分清洗单元、FFT单元、时延相位补偿单元、频域合成单元、互谱计算单元和精确残余多普勒估计单元;
频偏模型建立单元:用于根据外部输入目标轨位信息和各天线的位置坐标信息,建立天线1~天线N多普勒多项式预报值fdk(t),将多普勒多项式预报值送入正交下变频和积分清洗单元;其中,N为大于1的自然数,fdk(t)=a0+a1·t+a2·t2+…,t为时间,k=1,2,…,N为天线标识,a0、a1、a2…分别为多项式系数;
正交下变频和积分清洗单元:用于根据外部输入的中频载波标称值、频偏模型建立单元输入的多普勒多项式预报值和精确残余多普勒估计单元输入的精确残余多普勒估计值构造下变频本地信号slocal=exp(-j2π(f0+fdk(t)+Δfk0(t))t),利用构造的下变频本地信号对输入信号进行正交下变频,得到复数基带信号,对复数基带信号的I、Q支路分别进行积分清洗降低数据速率后,将积分清洗后的I、Q基带信号送入FFT单元;其中,f0为中频载波标称值,fdk(t)为多普勒多项式预报值,Δfk0(t)为精确残余多普勒估计值;
FFT单元:用于对积分清洗后的I、Q基带信号进行复数FFT运算,并将计算结果送入时延相位补偿单元;
时延相位补偿单元:用于利用互谱计算单元计算得到的残余时延差和残余相位差对FFT单元的计算结果进行时延相位精补偿,将精补偿后的结果送入频域合成单元和互谱计算单元;
频域合成单元:将精补偿后的信号进行幅度加权后频域合成,将频域合成信号进行IFFT后输出,并将频域合成信号送互谱计算单元;
互谱计算单元:根据精补偿后的信号和频域合成信号进行互谱计算,将互谱结果的互谱相位信息通过最小二乘线性拟合进行计算,得到残余时延差和残余相位差;将残余时延差和残余相位差送入时延相位补偿单元,同时将残余相位差送入精确残余多普勒估计单元,其中,第k路信号的残余相位差为:θk0为固定的相位差值,Δfk0(t)为精确残余多普勒频差;
精确残余多普勒估计单元:根据互谱计算单元得到的残余相位差,通过多项式拟合得到精确残余多普勒估计值,并将得到的精确残余多普勒估计值Δfk0(t)更新到正交下变频和积分清洗单元;其中,精确残余多普勒估计值Δfk0(t)=b0+b1·t+b2·t2+…,b0、b1、b2分别为多项式拟合系数。
Claims (1)
1.一种测控天线组阵精确频率偏差估计与补偿装置,其特征在于,包括频偏模型建立单元、正交下变频和积分清洗单元、FFT单元、时延相位补偿单元、频域合成单元、互谱计算单元和精确残余多普勒估计单元;
频偏模型建立单元:用于根据外部输入目标轨位信息和各天线的位置坐标信息,建立天线1~天线N多普勒多项式预报值fdk(t),将多普勒多项式预报值送入正交下变频和积分清洗单元;其中,N为大于1的自然数,fdk(t)=a0+a1·t+a2·t2+…,t为时间,k=1,2,…,N为天线标识,a0、a1、a2…分别为多项式系数;
正交下变频和积分清洗单元:用于根据外部输入的中频载波标称值、频偏模型建立单元输入的多普勒多项式预报值和精确残余多普勒估计单元输入的精确残余多普勒估计值构造下变频本地信号slocal=exp(-j2π(f0+fdk(t)+Δfk0(t))t),利用构造的下变频本地信号对输入信号进行正交下变频,得到复数基带信号,对复数基带信号的I、Q支路分别进行积分清洗降低数据速率后,将积分清洗后的I、Q基带信号送入FFT单元;其中,f0为中频载波标称值,fdk(t)为多普勒多项式预报值,Δfk0(t)为精确残余多普勒估计值;
FFT单元:用于对积分清洗后的I、Q基带信号进行复数FFT运算,并将计算结果送入时延相位补偿单元;
时延相位补偿单元:用于利用互谱计算单元计算得到的残余时延差和残余相位差对FFT单元的计算结果进行时延相位精补偿,将精补偿后的结果送入频域合成单元和互谱计算单元;
频域合成单元:将精补偿后的信号进行幅度加权后频域合成,将频域合成信号进行IFFT后输出,并将频域合成信号送互谱计算单元;
互谱计算单元:根据精补偿后的信号和频域合成信号进行互谱计算,将互谱结果的互谱相位信息通过最小二乘线性拟合进行计算,得到残余时延差和残余相位差;将残余时延差和残余相位差送入时延相位补偿单元,同时将残余相位差送入精确残余多普勒估计单元,其中,第k路信号的残余相位差为:θk0为固定的相位差值,Δfk0(t)为精确残余多普勒频差;
精确残余多普勒估计单元:根据互谱计算单元得到的残余相位差,通过多项式拟合得到精确残余多普勒估计值,并将得到的精确残余多普勒估计值Δfk0(t)更新到正交下变频和积分清洗单元;其中,精确残余多普勒估计值Δfk0(t)=b0+b1·t+b2·t2+…,b0、b1、b2分别为多项式拟合系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110839557.4A CN113556188B (zh) | 2021-07-23 | 2021-07-23 | 一种测控天线组阵精确频率偏差估计与补偿装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110839557.4A CN113556188B (zh) | 2021-07-23 | 2021-07-23 | 一种测控天线组阵精确频率偏差估计与补偿装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113556188A true CN113556188A (zh) | 2021-10-26 |
CN113556188B CN113556188B (zh) | 2022-07-19 |
Family
ID=78132724
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110839557.4A Active CN113556188B (zh) | 2021-07-23 | 2021-07-23 | 一种测控天线组阵精确频率偏差估计与补偿装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113556188B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115037404A (zh) * | 2022-08-11 | 2022-09-09 | 中国人民解放军战略支援部队航天工程大学 | 一种gpu平台下的阵列信号时延补偿方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104253774A (zh) * | 2014-09-16 | 2014-12-31 | 重庆邮电大学 | 一种高动态环境下多普勒频偏估计系统和方法 |
CN104407323A (zh) * | 2014-12-11 | 2015-03-11 | 中国工程物理研究院电子工程研究所 | 一种高动态低信噪比扩频信号伪码时延测量方法 |
CN108365872A (zh) * | 2017-12-29 | 2018-08-03 | 中国电子科技集团公司第五十四研究所 | 一种适用于异地多天线组阵的信号合成方法 |
CN109633574A (zh) * | 2018-10-25 | 2019-04-16 | 上海无线电设备研究所 | 一种用于深空探测的宽范围高精度多普勒测量方法 |
CN110545247A (zh) * | 2019-07-31 | 2019-12-06 | 成都天奥集团有限公司 | 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法 |
CN111510409A (zh) * | 2020-04-16 | 2020-08-07 | 清华大学 | 利用bpsk数据估计天基机会信号多普勒的方法和系统 |
-
2021
- 2021-07-23 CN CN202110839557.4A patent/CN113556188B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104253774A (zh) * | 2014-09-16 | 2014-12-31 | 重庆邮电大学 | 一种高动态环境下多普勒频偏估计系统和方法 |
CN104407323A (zh) * | 2014-12-11 | 2015-03-11 | 中国工程物理研究院电子工程研究所 | 一种高动态低信噪比扩频信号伪码时延测量方法 |
CN108365872A (zh) * | 2017-12-29 | 2018-08-03 | 中国电子科技集团公司第五十四研究所 | 一种适用于异地多天线组阵的信号合成方法 |
CN109633574A (zh) * | 2018-10-25 | 2019-04-16 | 上海无线电设备研究所 | 一种用于深空探测的宽范围高精度多普勒测量方法 |
CN110545247A (zh) * | 2019-07-31 | 2019-12-06 | 成都天奥集团有限公司 | 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法 |
CN111510409A (zh) * | 2020-04-16 | 2020-08-07 | 清华大学 | 利用bpsk数据估计天基机会信号多普勒的方法和系统 |
Non-Patent Citations (1)
Title |
---|
李乐天,王赛宇,王力男: "基于5G标准的低轨卫星通信初始同步技术", 《计算机测量与控制》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115037404A (zh) * | 2022-08-11 | 2022-09-09 | 中国人民解放军战略支援部队航天工程大学 | 一种gpu平台下的阵列信号时延补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113556188B (zh) | 2022-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107645465B (zh) | 对无线电信号进行定时和定位的方法和系统 | |
JP4462574B2 (ja) | パッシブコヒーレント探索アプリケーションのための、広帯域事前検出信号を処理するシステムおよび方法 | |
CN113472371B (zh) | 自适应阵列天线数字波束合成抗干扰处理方法 | |
CN101807977B (zh) | 基于波形特征的空时盲自适应抗干扰方法 | |
CN108769908B (zh) | 多径环境下基于doa/toa联合估计的车辆定位参数估计方法 | |
CN112505622B (zh) | 一种高精度单基站室内定位方法 | |
JP5089460B2 (ja) | 伝搬遅延時間測定装置及びレーダ装置 | |
CN113824517B (zh) | 一种基于数字波束合成的无线在轨自适应幅相校正系统 | |
CN109975754A (zh) | 联合角度、时延与多普勒信息的运动多站直接定位方法 | |
CN113556188B (zh) | 一种测控天线组阵精确频率偏差估计与补偿装置 | |
CN111190197A (zh) | 一种导航卫星信号质量在轨优化与保持方法 | |
CN107317781B (zh) | 一种水下通信中多普勒扩展抑制方法 | |
CN103869334B (zh) | Gnss空间信号畸变的自动识别与处理方法 | |
CN105954713A (zh) | 一种基于tdoa观测量定位算法的时延估计方法 | |
CN102508265A (zh) | 基于信号分离估计理论的卫星导航信号多径干扰抑制方法 | |
Wen et al. | Maximum likelihood and signal-selective TDOA estimation for noncircular signals | |
Zemmari | Reference signal extraction for GSM passive coherent location | |
CN103259638B (zh) | 一种异地本振一致性误差下的基带时差估计方法 | |
Falk et al. | Direction finding for electronic warfare systems using the phase of the cross spectral density | |
CN105388471B (zh) | 一种自适应电磁场时延估计方法及装置 | |
CN113589054A (zh) | 一种测控天线组阵自动时延偏差估计与补偿装置 | |
Xia et al. | Doppler rate-of-change estimation method for BPSK signal based on spectrum phase | |
Noschese et al. | A low-complexity approach for time of arrival estimation in OFDM systems | |
Neunteufel et al. | Bayesian CRLB for Blind Indoor Localization with Imperfect Receiver Synchronization | |
Önen et al. | A Time Difference of Arrival Estimation with Integration of Generalized Cross Correlation Samples for Continuous Wave Signals on Passive Positioning Systems |
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 |