CN111884984A - 快变载波多普勒频移捕获系统 - Google Patents

快变载波多普勒频移捕获系统 Download PDF

Info

Publication number
CN111884984A
CN111884984A CN202010608950.8A CN202010608950A CN111884984A CN 111884984 A CN111884984 A CN 111884984A CN 202010608950 A CN202010608950 A CN 202010608950A CN 111884984 A CN111884984 A CN 111884984A
Authority
CN
China
Prior art keywords
frequency
carrier
carrier doppler
unit
doppler
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
Application number
CN202010608950.8A
Other languages
English (en)
Other versions
CN111884984B (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.)
Southwest Electronic Technology Institute No 10 Institute of Cetc
Original Assignee
Southwest Electronic Technology Institute No 10 Institute of Cetc
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 Southwest Electronic Technology Institute No 10 Institute of Cetc filed Critical Southwest Electronic Technology Institute No 10 Institute of Cetc
Priority to CN202010608950.8A priority Critical patent/CN111884984B/zh
Publication of CN111884984A publication Critical patent/CN111884984A/zh
Application granted granted Critical
Publication of CN111884984B publication Critical patent/CN111884984B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38—Demodulator circuits; Receiver circuits
    • H04L27/3845—Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier
    • H04L27/3854—Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier using a non - coherent carrier, including systems with baseband correction for phase or frequency offset
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2614—Peak power aspects
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • H04L27/2659—Coarse or integer frequency offset determination and synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • H04L27/266—Fine or fractional frequency offset determination and synchronisation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开的一种快变载波多普勒频移捕获系统，旨在提供能够在高动态场景下有效、精确、快速地测量快变载波多普勒频移的捕获系统。本发明通过下述技术方案实现：第一级滤波采样单元对接收信号进行第一级滤波采样，第一级乒乓缓存单元实现信号采样的准实时处理，第一级频率补偿单元采用多个粗略频率子槽对第一级乒乓缓存单元中的采样数据完成载波频率的粗略补偿，第二级滤波采样单元对粗略补偿后的采样数据进行第二级滤波采样，第二级频率补偿单元采用多个精细频率子槽对第二级乒乓缓存单元中的采样数据完成载波频率的精细补偿，经过非线性变换单元、快速傅里叶变换单元、峰值搜索单元和频率解算单元，完成快变载波多普勒频移的精确捕获。

Description

快变载波多普勒频移捕获系统
技术领域
本发明属于无线通信领域，涉及一种快变载波多普勒频移捕获系统。
技术背景
在无线通信中，卫星、飞行器等目标常处于高速机动状态，运动载体存在速度、加速度,这将导致在通信过程中载波存在多普勒频移，由于接收载体通常也在快速运动，因此发送端与接收端处于复杂相对运动状态，使得通信信号出现“多普勒效应”的现象，则接收信号将存在一个快速变化的载波多普勒频移，载波多普勒一阶变化率和载波多普勒高阶变化率会使信号频谱出现“时变平台效应”。传统频域测量方法主要是采用基于FFT运算对载波频率进行估计检测，当不存在载波多普勒变化率或其较小时，基于FFT运算的算法虽然能较精确地得到载波多普勒频偏估计值，然而当多普勒变化率较大时，由于进行分段FFT运算产生的信号功率谱会发生“频域平移”现象，将大大减小非相干累加的有效性，降低接收信号的检测性能，若采用最大似然估计方法基于三维搜索的思想对载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率进行估计测量，由于计算量过于庞大，工程实现上不易实现。因此，传统捕获系统已无法精确测量快变载波多普勒频移。
传统的载波环在处理存在快变载波多普勒频移的接收信号时，如果不增加环路带宽，载波多普勒频移及其变化率会使载波超出载波环的捕获频带，而环路带宽的增加将会引入更多噪声导致准确度的降低，当引入噪声的电平接近或超过环路的门限电压时，将导致跟踪环路失锁，为了跟踪环路快速而稳定地入锁，且保证跟踪环路锁定过程中不易失锁，在开始信号跟踪之前进行信号捕获时，需要进一步提高载波频率的测量精度，使得信号捕获估算出的载波频率估计值与真实载波频率误差较小，如何有效、精确、快速地测量快变载波多普勒频移，已成为无线通信系统的核心技术。
传统捕获系统对接收信号进行滤波采样处理后得到采样数据，采样数据先直接进行非线性变换实现载波恢复，经过快速傅里叶变换(Fast Fourier Transform，FFT)完成时频变换，最后经过峰值搜索、频率解算得到载波频率。当采用傅立叶变换法估计载波的频偏时,其频谱分辨率与采样点数成反比,采样点数越多,频偏分辨率越小,跟踪精度越高。但是在采样间隔固定的条件下,采样点数的增加意味着数据的准备时间加大,远远大于采用DSP进行FFT运算所需的时间。因此采用傅立叶变换估计快变载波多普勒频移是不现实的。由于传统捕获系统只实现了载波多普勒频移的测量，没有测量载波多普勒一阶变化率和载波多普勒二阶变化率，无法校正由于目标高速机动、处理时间延迟造成的频率偏移，测频数据实时性低，频率测量精度差，无法适应实际接收系统的捕获性能需求。
发明内容
本发明针对传统捕获系统测频实时性低，频率测量精度差，无法精确测量快变载波多普勒频移的技术缺陷，提供一种处理实时性好、频率测量误差小、估测准确度高，能够在高动态场景下有效、精确、快速地测量各种调制信号载波频率的快变载波多普勒频移捕获系统。
本发明的上述目的可以通过以下介绍方案予以实现，一种快变载波多普勒频移捕获系统，包括：包括：顺次串联的第一级滤波采样单元、第一级乒乓缓存单元、第一级频率补偿单元、第二级滤波采样单元、第二级乒乓缓存单元、第二级频率补偿单元、非线性变换单元、快速傅里叶变换单元、峰值搜索单元和频率解算单元，其特征在于：第一级滤波采样单元根据载波多普勒频移范围对接收信号进行第一级滤波采样，通过第一级乒乓缓存单元实现信号采样的准实时处理；第一级频率补偿单元根据载波多普勒频移、载波多普勒一阶变化率范围划分多个粗略频率子槽，采用粗略频率子槽对第一级乒乓缓存单元中的采样数据完成载波频率的粗略补偿，粗略补偿后的采样数据经过第二级滤波采样单元进行第二级滤波采样，并进入第二级乒乓缓存单元；第二级频率补偿单元根据第一级频率补偿单元中载波多普勒一阶变化率搜索步进和载波多普勒二阶变化率范围划分为多个精细频率子槽，采用精细频率子槽对第二级乒乓缓存单元中的采样数据完成载波频率的精细补偿，精细补偿后的采样数据经过非线性变换单元、快速傅里叶变换单元和峰值搜索单元完成非线性变换、快速傅里叶变换和峰值搜索，得到积分峰值及其所在轮次信息；频率解算单元根据积分峰值及其所在轮次信息进行频率解算，结合信号采样时间和数据处理时间，对频率解算单元输出的捕获结果进行校正，得到当前时刻接收信号的载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率，完成快变载波多普勒频移的精确捕获。
本发明相对于传统捕获系统具有以下有益效果：
处理实时性好。本发明对接收信号的载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率进行高精度的测量，相比于传统捕获系统，解决了高动态场景的捕获技术问题，可适应高动态场景的捕获性能需求。
频率测量精度高。本发明采用频率解算单元，根据采样时间和数据处理时间，校正得到当前时刻接收信号的载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率，相比于传统捕获系统，进一步提高了载波频率的测量精度，有效地提升了多普勒频偏测量的可靠性。
适应各种调制信号。本发明采用非线性变换单元，根据接收信号的调制类型选择相对应的倍频模式，对采样数据进行非线性变换以实现载波恢复，相比于传统捕获系统，可适应各种不同调制类型的接收信号。
附图说明
图1是本发明快变载波多普勒频移捕获系统结构原理示意图。
图2是图1中第一级滤波采样单元、第二级滤波采样单元的结构原理示意图。
图3是图1中第一级频率补偿单元的结构原理示意图。
图4是图1中第二级频率补偿单元的结构原理示意图。
图5是图1中非线性变换单元的结构原理示意图。
下面结合附图和实施例对发明进一步说明。
具体实施方式
参阅图1。一种快变载波多普勒频移捕获系统，包括：包括：顺次串联的第一级滤波采样单元、第一级乒乓缓存单元、第一级频率补偿单元、第二级滤波采样单元、第二级乒乓缓存单元、第二级频率补偿单元、非线性变换单元、快速傅里叶变换单元、峰值搜索单元和频率解算单元，其中：第一级滤波采样单元根据载波多普勒频移范围对接收信号进行第一级滤波采样，通过第一级乒乓缓存单元实现信号采样的准实时处理；第一级频率补偿单元根据载波多普勒频移、载波多普勒一阶变化率范围划分多个粗略频率子槽，采用粗略频率子槽对第一级乒乓缓存单元中的采样数据完成载波频率的粗略补偿，粗略补偿后的采样数据经过第二级滤波采样单元进行第二级滤波采样，并进入第二级乒乓缓存单元；第二级频率补偿单元根据第一级频率补偿单元中载波多普勒一阶变化率搜索步进和载波多普勒二阶变化率范围划分为多个精细频率子槽，采用精细频率子槽对第二级乒乓缓存单元中的采样数据完成载波频率的精细补偿，精细补偿后的采样数据经过非线性变换单元、快速傅里叶变换单元和峰值搜索单元完成非线性变换、快速傅里叶变换和峰值搜索，得到积分峰值及其所在轮次信息；频率解算单元根据积分峰值及其所在轮次信息进行频率解算，结合信号采样时间和数据处理时间，对频率解算单元输出的捕获结果进行校正，得到当前时刻接收信号的载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率，完成快变载波多普勒频移的精确捕获。
参阅图2。第一级滤波采样单元根据接收信号的数据频率fad和采样频率fsamp采用公式fsamp×232/fad进行数据转换，将得到的采样频率控制字Ksamp1送入直接数字式频率合成器DDS，对输入数据进行清零脉冲和积分清零，结合第二级滤波采样单元的积分点数M和快速傅里叶变换单元的快速傅里叶变换点数N，得到信号采样时间tsamp为M×N/fsamp。
参阅图3。第一级频率补偿单元将载波多普勒频移范围划分为总数Ldopl+1个载波多普勒一阶变化率搜索子槽，载波多普勒频移搜索轮次ldopl∈(0,1，…，Ldopl)，载波多普勒频移搜索步进
Figure BDA0002560199660000041
粗略频率子槽内载波多普勒频移
Figure BDA0002560199660000042
为
Figure BDA0002560199660000043
并从最小载波多普勒频移
Figure BDA0002560199660000044
顺次搜索至最大载波多普勒频移
Figure BDA0002560199660000045
将载波多普勒一阶变化率范围划分成总数Lrate1+1个载波多普勒一阶变化率搜索子槽，载波多普勒一阶变化率搜索轮次lrate1∈(0,1，…，Lrate1)，载波多普勒一阶变化率搜索步进
Figure BDA0002560199660000046
粗略频率子槽内载波多普勒一阶变化率
Figure BDA0002560199660000047
为
Figure BDA0002560199660000048
并从最小载波多普勒一阶变化率
Figure BDA0002560199660000049
顺次搜索至最大载波多普勒一阶变化率
Figure BDA00025601996600000410
再采用公式
Figure BDA00025601996600000411
对粗略频率子槽内的载波多普勒频移
Figure BDA00025601996600000412
进行数据转换，得到载波多普勒频移控制字
Figure BDA00025601996600000413
同时采用公式
Figure BDA00025601996600000414
对粗略频率子槽内的载波多普勒一阶变化率
Figure BDA00025601996600000415
进行数据转换，得到载波多普勒一阶变化率控制字
Figure BDA00025601996600000416
最后对载波多普勒频移控制字
Figure BDA00025601996600000417
进行一次累加得到查询地址，通过地址映射、查表生成载波多普勒频移的本地载波，同时对载波多普勒一阶变化率控制字
Figure BDA00025601996600000418
进行两次累加得到查询地址，通过地址映射、查表生成载波多普勒一阶变化率的本地载波，利用两路生成的本地载波进行复乘运算得到本地复合载波，本地复合载波与第一级乒乓缓存单元输出的采样数据进行复乘运算完成载波频率的粗略补偿。
第二级滤波采样单元采用公式232/M得到第二级采样频率控制字为Ksamp2，对粗略补偿后的采样数据进行第二级滤波采样，将第二级滤波采样的采样数据存入第二级乒乓缓存单元，实现两级数据处理的准实时后，结合第一级频率补偿单元中载波多普勒频移搜索轮次总数Ldopl+1、载波多普勒一阶变化率搜索轮次总数Lrate1+1和两级数据处理时钟fsys，得到数据处理时间tdeal为(Ldopl+1)×(Lrate1+1)×M×N/fsys。
参阅图4。第二级频率补偿单元将第一级频率补偿单元中的载波多普勒一阶变化率搜索步进
Figure BDA00025601996600000419
划分成总数Lrate2+1个载波多普勒一阶变化率搜索子槽，载波多普勒一阶变化率搜索轮次lrate2∈(0,1，…，Lrate2)，载波多普勒一阶变化率搜索步进
Figure BDA00025601996600000420
精细频率子槽内载波多普勒一阶变化率
Figure BDA00025601996600000421
为
Figure BDA00025601996600000422
并从
Figure BDA00025601996600000423
顺次搜索至
Figure BDA00025601996600000424
将载波多普勒二阶变化率范围划分成总数Lacce+1个载波多普勒二阶变化率搜索子槽，载波多普勒二阶变化率搜索轮次lacce∈(0,1，…，Lacce)，载波多普勒二阶变化率搜索步进
Figure BDA0002560199660000051
精细频率子槽内载波多普勒二阶变化率
Figure BDA0002560199660000052
为
Figure BDA0002560199660000053
并从最小载波多普勒二阶变化率
Figure BDA0002560199660000054
顺次搜索至最大载波多普勒二阶变化率
Figure BDA0002560199660000055
再采用公式
Figure BDA0002560199660000056
对精细频率子槽内的载波多普勒一阶变化率
Figure BDA0002560199660000057
进行数据转换，得到载波多普勒一阶变化率控制字
Figure BDA0002560199660000058
同时采用公式
Figure BDA0002560199660000059
对精细频率子槽内的载波多普勒二阶变化率
Figure BDA00025601996600000510
进行数据转换，得到载波多普勒二阶变化率控制字
Figure BDA00025601996600000511
最后对载波多普勒一阶变化率控制字
Figure BDA00025601996600000512
进行两次累加得到查询地址，通过地址映射、查表生成载波多普勒一阶变化率的本地载波，同时对载波多普勒二阶变化率控制字
Figure BDA00025601996600000513
进行三次累加得到查询地址，通过地址映射、查表生成载波多普勒二阶变化率的本地载波，利用两路生成的本地载波进行复乘运算得到本地复合载波，本地复合载波与第二级乒乓缓存单元输出的采样数据进行复乘运算完成载波频率的精细补偿。
参阅图5。非线性变换单元将单载波信号、PM信号划为非线性变换单元的非线性变换参数ksq＝1的单倍频模式，双相相移键控(BPSK)信号划为ksq＝2的二倍频模式，正交相移键控(QPSK)信号划为ksq＝4的四倍频模式，八进制相移键控(8PSK)信号划分为ksq＝8的八倍频模式，根据接收信号的调制类型选择相对应的倍频模式，对精细补偿后的采样数据进行载波恢复，载波多普勒频移的测量精度
Figure BDA00025601996600000514
为fsamp/M/ksq/N。
快速傅里叶变换单元对载波恢复的数据进行N点快速傅里叶变换得到积分数据，峰值搜索单元对所有轮次的积分数据进行比较搜索，得到积分峰值及其所在轮次信息，轮次信息包括第一级数据处理中载波多普勒频移轮次
Figure BDA00025601996600000515
载波多普勒一阶变化率轮次
Figure BDA00025601996600000516
第二级数据处理中载波多普勒一阶变化率轮次
Figure BDA00025601996600000517
载波多普勒二阶变化率轮次
Figure BDA00025601996600000518
和快速傅里叶变换指数
Figure BDA00025601996600000519
频率解算单元根据峰值搜索单元输出的峰值信息分别得到载波多普勒频移测量值
Figure BDA00025601996600000520
为
Figure BDA00025601996600000521
载波多普勒一阶变化率测量值
Figure BDA00025601996600000522
为
Figure BDA00025601996600000523
载波多普勒二阶变化率测量值
Figure BDA00025601996600000524
为
Figure BDA00025601996600000525
结合信号采样时间tsamp和数据处理时间tdeal，校正得到当前时刻接收信号的载波多普勒频移
Figure BDA00025601996600000526
为
Figure BDA00025601996600000527
载波多普勒一阶变化率
Figure BDA00025601996600000528
为
Figure BDA00025601996600000529
载波多普勒二阶变化率
Figure BDA00025601996600000530
为
Figure BDA00025601996600000531
下面以范例具体分析：
设接收信号为BPSK信号，调制符号速率为10.0kbps，系统处理时钟为180MHz，载波多普勒频移范围为±120.0kHz，载波多普勒一阶变化率范围为±7.5kHz/s，载波多普勒一阶变化率范围为±0.5kHz/s2。
第一级滤波采样单元中采样频率fsamp为640.0kHz，第二级滤波采样单元的积分点数M为128，快速傅里叶变换单元的快速傅里叶变换点数N为1024，则信号采样时间tsamp为0.2048s；第一级频率补偿单元将载波多普勒频移以3.0kHz为搜索步进划分41个粗略多普勒子槽、载波多普勒一阶变化率以300.0Hz/s为搜索步进划分成67个粗略一阶变化率子槽，数据处理时间tdeal约为0.201s；第二级频率补偿单元将载波多普勒一阶变化率以5.0Hz/s为搜索步进划分成125个精细一阶变化率子槽，载波多普勒二阶变化率以75.0Hz/s2为搜索步进划分成8个精细二阶变化率子槽，后续数据处理采用8路并行；非线性变换单元选择二倍频模式；峰值搜索单元对所有轮次的积分数据进行比较搜索，得到积分峰值及其所在轮次信息；频率解算单元根据峰值搜索单元输出的峰值信息结合信号采样时间tsamp和数据处理时间tdeal，校正得到当前时刻接收信号的载波多普勒频移
Figure BDA0002560199660000061
载波多普勒一阶变化率
Figure BDA0002560199660000062
和载波多普勒二阶变化率
Figure BDA0002560199660000063
经过试验分析验证：在范例场景下，本发明的载波多普勒频移测量误差优于±5.0Hz、载波多普勒一阶变化率测量误差优于±10.0Hz/s、载波多普勒一阶变化率测量误差优于±50.0Hz/s2。
以上对本发明实施例进行了详细介绍，本文中应用了具体实施方式对本发明进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及设备；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种快变载波多普勒频移捕获系统，包括：顺次串联的第一级滤波采样单元、第一级乒乓缓存单元、第一级频率补偿单元、第二级滤波采样单元、第二级乒乓缓存单元、第二级频率补偿单元、非线性变换单元、快速傅里叶变换单元、峰值搜索单元和频率解算单元，其特征在于：第一级滤波采样单元根据载波多普勒频移范围对接收信号进行第一级滤波采样，通过第一级乒乓缓存单元实现信号采样的准实时处理；第一级频率补偿单元根据载波多普勒频移、载波多普勒一阶变化率范围划分多个粗略频率子槽，采用粗略频率子槽对第一级乒乓缓存单元中的采样数据完成载波频率的粗略补偿，粗略补偿后的采样数据经过第二级滤波采样单元进行第二级滤波采样，并进入第二级乒乓缓存单元；第二级频率补偿单元根据第一级频率补偿单元中载波多普勒一阶变化率搜索步进和载波多普勒二阶变化率范围划分为多个精细频率子槽，采用精细频率子槽对第二级乒乓缓存单元中的采样数据完成载波频率的精细补偿，精细补偿后的采样数据经过非线性变换单元、快速傅里叶变换单元和峰值搜索单元完成非线性变换、快速傅里叶变换和峰值搜索，得到积分峰值及其所在轮次信息；频率解算单元根据积分峰值及其所在轮次信息进行频率解算，结合信号采样时间和数据处理时间，对频率解算单元输出的捕获结果进行校正，得到当前时刻接收信号的载波多普勒频移、载波多普勒一阶变化率和载波多普勒二阶变化率，完成快变载波多普勒频移的精确捕获。
2.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：第一级滤波采样单元根据接收信号的数据频率fad和采样频率fsamp采用公式fsamp×232/fad进行数据转换，将得到的采样频率控制字Ksamp1送入直接数字式频率合成器DDS，对输入数据进行清零脉冲和积分清零，结合第二级滤波采样单元的积分点数M和快速傅里叶变换单元的快速傅里叶变换点数N，得到信号采样时间tsamp为M×N/fsamp。
3.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：第一级频率补偿单元将载波多普勒频移范围划分为总数Ldopl+1个载波多普勒一阶变化率搜索子槽，并且载波多普勒频移搜索轮次ldopl∈(0,1，…，Ldopl)，载波多普勒频移搜索步进
Figure FDA0002560199650000011
粗略频率子槽内载波多普勒频移
Figure FDA0002560199650000012
为
Figure FDA0002560199650000013
并从最小载波多普勒频移
Figure FDA0002560199650000014
顺次搜索至最大载波多普勒频移
Figure FDA0002560199650000015
同时将载波多普勒一阶变化率范围划分成总数Lrate1+1个载波多普勒一阶变化率搜索子槽，并且载波多普勒一阶变化率搜索轮次lrate1∈(0,1，…，Lrate1)，载波多普勒一阶变化率搜索步进
Figure FDA0002560199650000016
粗略频率子槽内载波多普勒一阶变化率
Figure FDA0002560199650000017
为
Figure FDA0002560199650000018
并从最小载波多普勒一阶变化率
Figure FDA0002560199650000019
顺次搜索至最大载波多普勒一阶变化率
Figure FDA0002560199650000021
再采用公式
Figure FDA0002560199650000022
对粗略频率子槽内的载波多普勒频移
Figure FDA0002560199650000023
进行数据转换，得到载波多普勒频移控制字
Figure FDA0002560199650000024
4.如权利要求3所述的快变载波多普勒频移捕获系统，其特征在于：第一级频率补偿单元采用公式
Figure FDA0002560199650000025
对粗略频率子槽内的载波多普勒一阶变化率
Figure FDA0002560199650000026
进行数据转换，得到载波多普勒一阶变化率控制字
Figure FDA0002560199650000027
最后对载波多普勒频移控制字
Figure FDA0002560199650000028
进行一次累加得到查询地址，通过地址映射、查表生成载波多普勒频移的本地载波，同时对载波多普勒一阶变化率控制字
Figure FDA0002560199650000029
进行两次累加得到查询地址，通过地址映射、查表生成载波多普勒一阶变化率的本地载波，利用两路生成的本地载波进行复乘运算得到本地复合载波，本地复合载波与第一级乒乓缓存单元输出的采样数据进行复乘运算完成载波频率的粗略补偿。
5.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：第二级滤波采样单元采用公式232/M得到第二级采样频率控制字为Ksamp2，对粗略补偿后的采样数据进行第二级滤波采样，将第二级滤波采样的采样数据存入第二级乒乓缓存单元，实现两级数据处理的准实时后，结合第一级频率补偿单元中载波多普勒频移搜索轮次总数Ldopl+1、载波多普勒一阶变化率搜索轮次总数Lrate1+1和两级数据处理时钟fsys，得到数据处理时间tdeal为(Ldopl+1)×(Lrate1+1)×M×N/fsys。
6.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：第二级频率补偿单元将第一级频率补偿单元中的载波多普勒一阶变化率搜索步进
Figure FDA00025601996500000210
划分成总数Lrate2+1个载波多普勒一阶变化率搜索子槽，并且载波多普勒一阶变化率搜索轮次lrate2∈(0,1，…，Lrate2)，载波多普勒一阶变化率搜索步进
Figure FDA00025601996500000211
精细频率子槽内载波多普勒一阶变化率
Figure FDA00025601996500000212
为
Figure FDA00025601996500000213
并从
Figure FDA00025601996500000214
顺次搜索至
Figure FDA00025601996500000215
同时将载波多普勒二阶变化率范围划分成总数Lacce+1个载波多普勒二阶变化率搜索子槽，并且载波多普勒二阶变化率搜索轮次lacce∈(0,1，…，Lacce)，载波多普勒二阶变化率搜索步进
Figure FDA00025601996500000216
精细频率子槽内载波多普勒二阶变化率
Figure FDA00025601996500000217
为
Figure FDA00025601996500000218
并从最小载波多普勒二阶变化率
Figure FDA00025601996500000219
顺次搜索至最大载波多普勒二阶变化率
Figure FDA00025601996500000220
7.如权利要求6所述的快变载波多普勒频移捕获系统，其特征在于：第二级频率补偿单元采用公式
Figure FDA00025601996500000221
对精细频率子槽内的载波多普勒一阶变化率
Figure FDA00025601996500000222
进行数据转换，得到载波多普勒一阶变化率控制字
Figure FDA0002560199650000031
同时采用公式
Figure FDA0002560199650000032
对精细频率子槽内的载波多普勒二阶变化率
Figure FDA0002560199650000033
进行数据转换，得到载波多普勒二阶变化率控制字
Figure FDA0002560199650000034
最后对载波多普勒一阶变化率控制字
Figure FDA0002560199650000035
进行两次累加得到查询地址，通过地址映射、查表生成载波多普勒一阶变化率的本地载波，同时对载波多普勒二阶变化率控制字
Figure FDA0002560199650000036
进行三次累加得到查询地址，通过地址映射、查表生成载波多普勒二阶变化率的本地载波，利用两路生成的本地载波进行复乘运算得到本地复合载波，本地复合载波与第二级乒乓缓存单元输出的采样数据进行复乘运算完成载波频率的精细补偿。
8.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：非线性变换单元将单载波信号、PM信号划为非线性变换单元的非线性变换参数ksq＝1的单倍频模式，双相相移键控(BPSK)信号划为ksq＝2的二倍频模式，正交相移键控(QPSK)信号划为ksq＝4的四倍频模式，八进制相移键控(8PSK)信号划分为ksq＝8的八倍频模式，根据接收信号的调制类型选择相对应的倍频模式，对精细补偿后的采样数据进行载波恢复，载波多普勒频移的测量精度
Figure FDA0002560199650000037
为fsamp/M/ksq/N。
9.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：快速傅里叶变换单元对载波恢复的数据进行N点快速傅里叶变换得到积分数据，峰值搜索单元对所有轮次的积分数据进行比较搜索，得到积分峰值及其所在轮次信息，其中，轮次信息包括第一级数据处理中载波多普勒频移轮次
Figure FDA0002560199650000038
载波多普勒一阶变化率轮次
Figure FDA0002560199650000039
第二级数据处理中载波多普勒一阶变化率轮次
Figure FDA00025601996500000310
载波多普勒二阶变化率轮次
Figure FDA00025601996500000311
和快速傅里叶变换指数
Figure FDA00025601996500000312
10.如权利要求1所述的快变载波多普勒频移捕获系统，其特征在于：频率解算单元根据峰值搜索单元输出的峰值信息分别得到载波多普勒频移测量值
Figure FDA00025601996500000313
为
Figure FDA00025601996500000314
载波多普勒一阶变化率测量值
Figure FDA00025601996500000315
为
Figure FDA00025601996500000316
载波多普勒二阶变化率测量值
Figure FDA00025601996500000317
为
Figure FDA00025601996500000318
结合信号采样时间tsamp和数据处理时间tdeal，校正得到当前时刻接收信号的载波多普勒频移
Figure FDA00025601996500000319
为
Figure FDA00025601996500000320
载波多普勒一阶变化率
Figure FDA00025601996500000321
为
Figure FDA00025601996500000322
载波多普勒二阶变化率
Figure FDA00025601996500000323
为
Figure FDA00025601996500000324
CN202010608950.8A 2020-06-29 2020-06-29 快变载波多普勒频移捕获系统 Active CN111884984B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010608950.8A CN111884984B (zh) 2020-06-29 2020-06-29 快变载波多普勒频移捕获系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010608950.8A CN111884984B (zh) 2020-06-29 2020-06-29 快变载波多普勒频移捕获系统

Publications (2)

Publication Number Publication Date
CN111884984A true CN111884984A (zh) 2020-11-03
CN111884984B CN111884984B (zh) 2022-09-02

Family

ID=73158177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010608950.8A Active CN111884984B (zh) 2020-06-29 2020-06-29 快变载波多普勒频移捕获系统

Country Status (1)

Country Link
CN (1) CN111884984B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112910499A (zh) * 2021-02-10 2021-06-04 西南电子技术研究所（中国电子科技集团公司第十研究所） 扩频信号精确捕获系统
CN112953873A (zh) * 2021-02-10 2021-06-11 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态微弱8psk/16psk信号载波捕获方法
CN112968850A (zh) * 2021-02-10 2021-06-15 西南电子技术研究所（中国电子科技集团公司第十研究所） 微弱bpsk信号载波捕获方法
WO2022126844A1 (zh) * 2020-12-17 2022-06-23 重庆邮电大学 基于沃特森模型的短波航空移动信道建模方法
US20230228885A1 (en) * 2022-01-18 2023-07-20 General Dynamics Mission Systems, Inc. Doppler tracking for signals of opportunity using rate-line detection

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0429769A2 (en) * 1989-11-22 1991-06-05 Pioneer Electronic Corporation A GPS satellite signal tracking method for GPS receivers
US5432521A (en) * 1994-01-14 1995-07-11 Motorola, Inc. Satellite receiver system having doppler frequency shift tracking
US20080272960A1 (en) * 2007-05-01 2008-11-06 Nobuhiro Kishimoto Method and system for GPS position measuring and frequency error detecting method
CN104218972A (zh) * 2014-08-27 2014-12-17 西安空间无线电技术研究所 一种跳扩频码相位和载波多普勒三维快速捕获方法
CN104378316A (zh) * 2014-10-30 2015-02-25 深圳市国创新能源研究院 一种多普勒频偏估计方法和装置
CN107835035A (zh) * 2017-10-30 2018-03-23 北京理工大学 低信噪比短帧突发通信开环解调方法及装置
CN108055058A (zh) * 2017-11-22 2018-05-18 西南电子技术研究所（中国电子科技集团公司第十研究所） 载波多普勒及其变化率的高精度测量方法
CN108809354A (zh) * 2018-06-04 2018-11-13 南京邮电大学 一种高动态多普勒捕获的长时间累积方法
CN109495410A (zh) * 2018-09-28 2019-03-19 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态pcm/fm信号载波频率精确估计方法
CN110417697A (zh) * 2019-06-29 2019-11-05 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态微弱mpsk信号的精确测频方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0429769A2 (en) * 1989-11-22 1991-06-05 Pioneer Electronic Corporation A GPS satellite signal tracking method for GPS receivers
US5432521A (en) * 1994-01-14 1995-07-11 Motorola, Inc. Satellite receiver system having doppler frequency shift tracking
US20080272960A1 (en) * 2007-05-01 2008-11-06 Nobuhiro Kishimoto Method and system for GPS position measuring and frequency error detecting method
CN104218972A (zh) * 2014-08-27 2014-12-17 西安空间无线电技术研究所 一种跳扩频码相位和载波多普勒三维快速捕获方法
CN104378316A (zh) * 2014-10-30 2015-02-25 深圳市国创新能源研究院 一种多普勒频偏估计方法和装置
CN107835035A (zh) * 2017-10-30 2018-03-23 北京理工大学 低信噪比短帧突发通信开环解调方法及装置
CN108055058A (zh) * 2017-11-22 2018-05-18 西南电子技术研究所（中国电子科技集团公司第十研究所） 载波多普勒及其变化率的高精度测量方法
CN108809354A (zh) * 2018-06-04 2018-11-13 南京邮电大学 一种高动态多普勒捕获的长时间累积方法
CN109495410A (zh) * 2018-09-28 2019-03-19 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态pcm/fm信号载波频率精确估计方法
CN110417697A (zh) * 2019-06-29 2019-11-05 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态微弱mpsk信号的精确测频方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LEI JIN等: ""Research on inertial information aided GPS tracking loop in high dynamic environment"", 《IEEE》 *
LIMIN YU 等: ""Optimum Receiver Design for Broadband Doppler Compensation in Multipath/Doppler Channels With Rational Orthogonal Wavelet Signaling"", 《IEEE》 *
余怀瑾等: ""基于UKF算法的高机动机载平台多普勒频移估计"", 《北京大学学报（自然科学版）》 *
金磊: ""高动态环境下载波频率的精确估计算法"", 《航天控制》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022126844A1 (zh) * 2020-12-17 2022-06-23 重庆邮电大学 基于沃特森模型的短波航空移动信道建模方法
CN112910499A (zh) * 2021-02-10 2021-06-04 西南电子技术研究所（中国电子科技集团公司第十研究所） 扩频信号精确捕获系统
CN112953873A (zh) * 2021-02-10 2021-06-11 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态微弱8psk/16psk信号载波捕获方法
CN112968850A (zh) * 2021-02-10 2021-06-15 西南电子技术研究所（中国电子科技集团公司第十研究所） 微弱bpsk信号载波捕获方法
CN112953873B (zh) * 2021-02-10 2022-07-29 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态微弱8psk/16psk信号载波捕获方法
US20230228885A1 (en) * 2022-01-18 2023-07-20 General Dynamics Mission Systems, Inc. Doppler tracking for signals of opportunity using rate-line detection
US11841442B2 (en) * 2022-01-18 2023-12-12 General Dynamics Mission Systems, Inc. Doppler tracking for signals of opportunity using rate-line detection

Also Published As

Publication number Publication date
CN111884984B (zh) 2022-09-02

Similar Documents

Publication Publication Date Title
CN111884984B (zh) 快变载波多普勒频移捕获系统
CN109495410B (zh) 高动态pcm/fm信号载波频率精确估计方法
CN110417697B (zh) 高动态微弱mpsk信号的精确测频方法
CN103630888B (zh) 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN106546817B (zh) 一种具有反馈功能的频率估计和能量估计方法
CN107797099B (zh) 一种多通道数字接收机实时内定标处理方法及装置
CN104698459B (zh) 一种应用于缺失数据的条带sar压缩感知成像方法
CN112953873B (zh) 高动态微弱8psk/16psk信号载波捕获方法
CN103837740A (zh) 一种高精度数字瞬时测频方法及装置
CN112910499B (zh) 扩频信号精确捕获系统
CN109655847B (zh) 一种适于动态信号的快速捕获方法
CN111865865B (zh) 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法
CN111123312A (zh) 一种高动态低信噪比环境下的卫星弱信号捕获方法及系统
CN113253313A (zh) 一种改进pmf-fft系统多普勒频率估计精度的方法
CN110708267B (zh) 频偏信息估计值确定方法
CN111049772A (zh) 应用于矢量信号分析仪平台的实现5g信号同步处理的系统及其方法
CN112859123B (zh) 一种快速捕获gps信号的方法
CN113885003A (zh) 一种基于相参雷达的非相参积累方式的目标探测方法
CN112968850A (zh) 微弱bpsk信号载波捕获方法
CN105044459A (zh) 一种谐波分析方法
CN117214928A (zh) 一种伪码跟踪方法、装置、终端设备及存储介质
CN115469341A (zh) 一种gps l2c信号cl码的快速捕获方法
KR101629493B1 (ko) 다중 파일럿 신호간의 진폭 및 위상 차 추출 장치 및 방법
CN112883787B (zh) 一种基于频谱匹配的短样本低频正弦信号参数估计方法
JP2005506799A (ja) ディジタル信号受信器におけるタイミング・エラーを推定するよう、訓練列を用いる装置及び方法

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