CN114966564A - 一种宽带线性调频信号多通道并行去斜方法 - Google Patents

一种宽带线性调频信号多通道并行去斜方法 Download PDF

Info

Publication number
CN114966564A
CN114966564A CN202210558439.0A CN202210558439A CN114966564A CN 114966564 A CN114966564 A CN 114966564A CN 202210558439 A CN202210558439 A CN 202210558439A CN 114966564 A CN114966564 A CN 114966564A
Authority
CN
China
Prior art keywords
signals
paths
parallel
signal
frequency
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.)
Pending
Application number
CN202210558439.0A
Other languages
English (en)
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.)
Hunan Guoke Ruicheng Electronic Technology Co ltd
Original Assignee
Hunan Guoke Ruicheng Electronic Technology Co ltd
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 Hunan Guoke Ruicheng Electronic Technology Co ltd filed Critical Hunan Guoke Ruicheng Electronic Technology Co ltd
Priority to CN202210558439.0A priority Critical patent/CN114966564A/zh
Publication of CN114966564A publication Critical patent/CN114966564A/zh
Pending legal-status Critical Current

Links

Images

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
    • 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/28—Details of pulse systems
    • G01S7/285—Receivers
    • 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/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H2017/0072—Theoretical filter design
    • H03H2017/0081—Theoretical filter design of FIR filters
    • 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)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种宽带线性调频信号多通道并行去斜方法，方法包括：利用高采样率ADC器件对大带宽中频回波信号进行高质量带通采样，并根据FPGA时钟速率做并行分解成M路数字信号；对M路回波信号进行多通道并行数字下变频和低通滤波处理，得到M路零中频IQ信号；基于直接数字频率合成原理，分M路并行产生相同带宽、时宽且斜率相反的线性调频信号作为去斜参考信号；将M路零中频回波信号与M路去斜参考信号相乘，然后经过M倍多相抽取滤波及FFT处理，得到去斜后带有目标距离信息的单频点回波信号。本发明将高速率大带宽回波信号的数字去斜处理转换为多通道并行处理，实现简单、资源利用率高、计算延时小。

Description

一种宽带线性调频信号多通道并行去斜方法
技术领域
本发明属于雷达信号处理技术领域，尤其是涉及一种宽带线性调频信号多通道并行去斜方法。
背景技术
随着科学技术的不断发展，为了使雷达系统具有高的距离分辨力，往往需要雷达发射信号波形具有宽的信号带宽，其原因是雷达的距离分辨力与雷达发射信号带宽成反比关系。在此前提下，雷达信号设计人员一般采用对雷达信号进行调制来提高信号带宽，常用的调制方式有调频和调相两种方式。在调频方式中，主要有线性调频信号和非线性调频信号两种信号样式，相较于非线性调频信号，线性调频信号除了具有宽的工作带宽之外，也易于硬件实现，因此线性调频信号是现阶段高性能雷达系统常用的信号样式之一。线性调频信号的原理参见非专利文献1-3。
但随着雷达发射信号带宽的不断提升，对雷达系统的硬件性能也提出了更高的要求，例如需要更高的信号采样频率，更高地数据处理速率以及更高的数据传输速率等。为了应对高速信号采集及数据处理对系统带来的挑战，一般通过对宽带线性调频信号进行去斜处理来降低回波信号带宽，即降低了信号的采样率，同时又可以保留完整的目标信息。常见的去斜处理方式有两种，分别是模拟去斜和数字去斜，具体参见非专利文献4-10，下面简要介绍这两种去斜方式。
模拟去斜是将雷达接收到的模拟的目标回波信号与模拟参考信号相乘，然后经过滤波放大，ADC采样，快速傅里叶变换(FFT)之后得到单频正弦波形式的回波信号，如图1所示。此时，单频信号频率与目标的距离等信息相关联，表达式如下：
Figure RE-GDA0003745856920000021
式中K为调频信号斜率，R为目标距离，R0为目标参考中心距离，c为真空中光速。模拟去斜后得到了包含目标信息的单频点信号，降低了回波信号带宽，从而降低了雷达系统对采样率的需求，但是模拟去斜方法要求产生模拟的去斜参考信号，增加了系统的复杂性和实现难度。为了更方便的实现大带宽线性调频信号的去斜处理，可以在数字域对回波信号进行去斜处理，即数字去斜方法。处理流程如图2所示，即先对中频回波信号进行AD采样后变成数字信号，然后经过数字下变频和低通滤波得到基带回波信号，再将其与数字产生的去斜参考信号相乘，最后经快速傅里叶变换(FFT)后得到单频点信号。数字去斜同样可以将回波信号变换成带有目标信息的单频点信号，与此同时，数字去斜参考信号更容易生成，参数调整更加方便。
数字去斜方法要求信号处理系统以较高的采样速率对模拟中频信号进行采样，但随着AD芯片的采样速率越来越高，采样速率将不再限制数字去斜方法的应用和发展。与此同时，非专利文献11提出基于FPGA高速可编程逻辑器件，可以完成极其复杂的高速高密度数字电路设计，从而实现复杂的数字信号处理功能。虽然现阶段FPGA芯片运算速率已经高达几百兆赫兹，但是为适应AD芯片数吉赫兹的信号采样速率，信号处理算法的多通道并行实现方法被广泛应用于各种数字信号处理领域。
在先文献
非专利文献
1.丁鹭飞，耿富录，陈建春.雷达原理(第4版)[M].电子工业出版社,2009.
2.理查兹，邢孟道，王彤,et al.雷达信号处理基础:Fundamentals of radarsignal processing[M].电子工业出版社,2008.
3.万永伦.宽带雷达信号产生与处理技术研究[D],2006.电子科技大学.
4.李和平，王岩飞，张弛,等.雷达回波信号的去斜接收方法,2016.
5.骆飞，王岩飞，牛晓丽,等.一种SAR高效数字去斜方法的研究[J].电子测量技术,2016,39(9):165-171.
6.万永伦，姒强，吕幼新.宽带信号去斜脉冲压缩处理方法的研究[J]，航空学报，2006,27(5)：917-921.
7.Pollock B,Goodman N A.Structured de-chirp for compressive samplingof LFM waveforms[J].7th Sensor Array and Multichannel Signal ProcessingWorkshop(SAM)，2012:37-40.
8.骆飞,王岩飞,牛晓丽,等.一种SAR回波信号去斜方法[P]，中国科学院电子学研究所，2016-07-26.
9.汪欣，陈海红.宽带信号全数字去斜与脉压方法研究[J].现代雷达， 2011,33(9):34-37.
10.Qian-Qiang L,Peng-Fei T,Bin Y,et al.A new dechirp method forwideband radar direct IF sampling signal[C]，Signal Processing(ICSP),2012 IEEE11th International Conference on.2012:1920-1924.
11.高亚军.基于FPGA的数字信号处理.第2版[M].电子工业出版社, 2015.
发明内容
本发明公开了一种宽带线性调频信号多通道并行去斜方法，可以在FPGA 运算速率要求范围内，在数字域对高采样率大带宽的线性调频信号完成去斜处理，同时兼具实现简单、资源利用率高、计算延时小的特点。
本发明提出的宽带线性调频信号多通道并行去斜方法，包括以下步骤：
利用高采样率ADC器件对大带宽中频回波信号进行高质量带通采样，并根据FPGA时钟速率做并行分解成M路数字信号；
对M路回波信号进行多通道并行数字下变频和低通滤波处理，得到M路零中频IQ信号；
基于直接数字频率合成原理，分M路并行产生相同带宽、时宽且斜率相反的线性调频信号作为去斜参考信号；
将M路零中频回波信号与M路去斜参考信号相乘，然后经过M倍多相抽取滤波及FFT处理，得到去斜后采样率为fs/M且带有目标距离信息的单频点回波信号，完成对整个大带宽线性调频信号数字去斜处理的多通道并行实现。
进一步的，基于带通采样定理，ADC采样频率满足下式：
Figure RE-GDA0003745856920000041
式中f0为信号中心频率，B为信号带宽，n取fs≥2B的最大正整数；
利用满足上式的ADC采样频率fs对回波信号进行中频采样得到数字回波信号x[n]，且将x[n]依据FPGA时钟工作频率fclock并行分解为M路回波信号，此时满足fs＝M*fclock，然后对M路信号进行数字下变频以及多通道并行低通滤波，所述数字下变频过程将M路以f0为中心频率的中频信号变成零中频基带信号，计算过程如下：
Figure RE-GDA0003745856920000042
式中cos(2πf0nTs)与sin(2πf0nTs)均由M路DDS并行产生；
经过数字下变频之后，得到零中频M路I信号和M路Q信号，将此IQ信号通过多通道并行低通滤波器进行滤波，从而滤除高频分量，得到零中频回波信号。
进一步的，对于2路并行的4阶FIR滤波器，所述低通滤波器的滤波过程包括：
滤波器输出y[n]表示为：
Figure RE-GDA0003745856920000051
式中h(m)为滤波器系数，将上式分为偶数ye[n]与奇数yo[n]两部分之和：
Figure RE-GDA0003745856920000052
ye(n)＝h(0)x(n)+h(2)x(n-2)
yo(n)＝h(1)x(n-1)+h(3)x(n-3)
同理将y(n+1)分为：
Figure RE-GDA0003745856920000053
ye(n+1)＝h(0)x(n+1)+h(2)x(n-1)
yo(n+1)＝h(1)x(n)+h(3)x(n-2)
ye(n+1)与yo(n)输入数据相同，同时yo(n+1)与ye(n)输入数据相差两个时钟周期，因此通过FIR滤波器进行并行滤波。
进一步的，所述FIR滤波器进行并行滤波包括：在n时刻，第一路滤波输出结果为第一路n时刻数据滤波结果的偶数部分和第二路n-1时刻数据滤波结果的奇数部分延时一个时钟周期后的和；第二路滤波输出结果为第二路n时刻数据滤波结果的偶数部分与第一路n时刻数据滤波结果的奇数部分之和，依次类推，将此滤波器并行实现结构扩展到n阶M路。
进一步的，利用DDS并行产生M路大带宽数字去斜参考信号包括：
线性调频信号表达式如下：
Sr[n]＝A cos[2πf0nTs+πk(nTs)2]
式中n∈[0，N]，Ts＝1/fs，
n+1时刻的表达式如下：
Sr[n+1]＝A cos[2πf0(n+1)Ts+πk(n+1)2Ts 2]
当M路并行产生时，n＝Mi+j，其中
Figure RE-GDA0003745856920000061
j∈[0，M]，上述两式中相位部分变换为：
Figure RE-GDA0003745856920000062
Figure RE-GDA0003745856920000063
式中K＝kTs/2＝k/2fs，因此：
Figure RE-GDA0003745856920000064
其中，记
Figure RE-GDA0003745856920000065
基于直接数字频率合成原理，得到M路不同的相位步进值，从而将大带宽线性调频去斜参考信号分M路并行产生。
进一步的，将并行生成的M路线性调频去斜参考信号Irm-j*Qrm与经过处理后的M路零中频回波信号Ix+j*Qm相乘，得到去斜后的Idm与Qdm两路回波信号，如下式所示：
(Im+jQm)·(Irm-jQrm)＝-Im·Qrm+Qm·Irm+j(Im·Irm+Qm·Qrm)，m＝0,1,…M
其中，Irm为M路线性调频去斜参考信号中的I路信号，Qrm为M路线性调频去斜参考信号的Q路信号，Im为M路零中频回波信号中的I路信号，Qm为M 路零中频回波信号的Q路信号；
将上述去斜之后的两路回波信号经过M倍多相抽取滤波以及FFT变换之后，得到采样率为fs/M且带有目标距离等信息的单频点回波信号。
进一步的，对回波信号进行处理时，首先将其并行分解为16路，并对此 16路线性调频回波信号进行多通道并行混频滤波，过程如下：
产生M路混频参考信号：
Figure RE-GDA0003745856920000071
同理：
Figure RE-GDA0003745856920000072
因此将16回波信号与上述16路混频参考信号一一对应相乘，在此过程中只需做一次乘法运算；
混频之后的信号再通过多通道并行滤波处理，在经过混频滤波处理后，将16路零中频I信号与16路零中频Q信号与16路并行生成的去斜参考信号相乘，做复数乘法时，每路只需做两次乘法运算和一次加法运算即可完成去斜处理运算过程。
本发明的有益效果是适应FGPA等硬件的时钟计算速率，将高速率大带宽回波信号的数字去斜处理转换为多通道并行处理，实现简单、资源利用率高、计算延时小。
附图说明
图1为雷达回波信号模拟去斜处理流程；
图2为雷达回波信号数字去斜处理流程；
图3为本发明的雷达回波信号数字去斜的多通道并行实现流程；
图4为本发明使用2路4阶FIR低通滤波器多路并行分解结构。
图5为本发明多通道并行数字去斜处理过程中的信号示例。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的说明书和权利要求书的术语“第一”、“第二”、“第三”、“第四”等(如果存在)或“左”、“右”、“前”、“后”、“顶”、“底”是用于区别类似的对象或便于本发明的结构描述，而不必用于描述特定的顺序或先后次序以及限制本发明的结构技术特征。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例，例如能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例
为了使ADC采样之后的高速率回波信号能够在FPGA等硬件系统中进行快速低延时处理，本发明提出了一种基于FPGA的大带宽线性调频信号多通道并行的零中频数字去斜实现方法。处理流程如图3所示。
基于带通采样定理，ADC采样频率应满足下式：
Figure RE-GDA0003745856920000081
式中f0为信号中心频率，B为信号带宽，n取fs≥2B的最大正整数。
利用满足上述条件的ADC采样频率fs对回波信号进行中频采样得到数字回波信号x[n]，且将x[n]依据FPGA时钟工作频率fclock并行分解为M路回波信号，此时满足fs＝M*fclock，然后对M路信号进行数字下变频(DDC)以及多通道并行低通滤波，数字低通滤波器选用常用的FIR滤波器。数字下变频过程将M路以f0为中心频率的中频信号变成零中频基带信号。计算过程如下：
Figure RE-GDA0003745856920000091
式中cos(2πf0nTs)与sin(2πf0nTs)均由M路DDS并行产生。
经过DDC之后，我们可以得到零中频M路I信号和M路Q信号，将此IQ信号通过多通道并行低通滤波器进行滤波,从而滤除高频分量，并行滤波过程如下，以2路并行的4阶FIR滤波器为例，滤波器输出y[n]可表示为：
Figure RE-GDA0003745856920000092
式中h(m)为滤波器系数，将上式分为偶数ye[n]与奇数yo[n]两部分之和：
Figure RE-GDA0003745856920000093
ye(n)＝h(0)x(n)+h(2)x(n-2)
yo(n)＝h(1)x(n-1)+h(3)x(n-3)
同理将y(n+1)分为：
Figure RE-GDA0003745856920000094
ye(n+1)＝h(0)x(n+1)+h(2)x(n-1)
yo(n+1)＝h(1)x(n)+h(3)x(n-2)
由y(n)和y(n+1)分解之后的表达式可以看出，ye(n+1)与yo(n)输入数据相同，同时yo(n+1)与ye(n)输入数据相差两个时钟周期，因此可利用图4 形式的FIR滤波器结构进行并行滤波。在n时刻，第一路滤波输出结果为第一路n时刻数据滤波结果的偶数部分和第二路n-1时刻数据滤波结果的奇数部分延时一个时钟周期后的和；第二路滤波输出结果为第二路n时刻数据滤波结果的偶数部分与第一路n时刻数据滤波结果的奇数部分之和。依次类推，可以将此滤波器并行实现结构扩展到n阶M路。多通道并行低通滤波之后得到零中频回波信号。本实施例以2路4阶的FIR滤波器结构为例，在此并不限制使用其它形式的FIR滤波器结构。
下一步利用DDS并行产生M路大带宽数字去斜参考信号，线性调频信号表达式如下：
Sr[n]＝A cos[2πf0nTs+πk(nTs)2]
式中n∈[0，N]，Ts＝1/fs。n+1时刻的表达式如下：
Sr[n+1]＝A cos[2πf0(n+1)Ts+πk(n+1)2Ts 2]
当M路并行产生时，n＝Mi+j，其中
Figure RE-GDA0003745856920000101
j∈[0，M]，上述两式中相位部分可变换为：
Figure RE-GDA0003745856920000102
Figure RE-GDA0003745856920000103
式中K＝kTs/2＝k/2fs。因此：
Figure RE-GDA0003745856920000104
其中，记ΔXf0(n)＝Mf0Ts,ΔXk(n)＝K[M(Mi+j)+M(Mi+j+1)+M2]Ts。基于 DDS原理，可以根据设计累加值，得到M路不同的相位步进值，从而将大带宽线性调频去斜参考信号分M路并行产生。
将并行生成的M路线性调频去斜参考信号(Irm-j*Qrm)与经过处理后的M路零中频回波信号(Im+j*Qm)相乘，我们可以得到去斜后的Idm与Qdm两路回波信号，如下式所示：
(Im+jQm)·(Irm-jQrm)＝-Im·Qrm+Qm·Irm+j(Im·Irm+Qm·Qrm)，m＝0,1,…M 将上述去斜之后的回波信号经过M倍多相抽取滤波以及FFT变换之后，我们
可以得到采样率为fs/M且带有目标距离等信息的单频点回波信号。
至此完成对大带宽线性调频信号数字去斜处理的多通道并行实现，将该方法总结如下：
(1)利用高采样率ADC器件对大带宽中频回波信号进行高质量带通采样，并综合考虑FPGA时钟速率做并行分解成M路数字信号。
(2)对M路回波信号进行多通道并行数字下变频(DDC)和低通滤波处理，得到M路零中频IQ信号。
(3)基于DDS原理，分M路并行产生相同带宽、时宽且斜率相反的线性调频信号作为去斜参考信号。
(4)将M路零中频回波信号与M路去斜参考信号相乘，然后经过M倍多相抽取滤波及FFT处理，得到去斜后采样率为fs/M且带有目标距离等信息的单频点回波信号，从而完成对整个大带宽线性调频信号数字去斜处理的多通道并行实现。
下面结合实际工程数据，根据上述方法进行大带宽线性调频信号数字去斜处理。
设定系统工作中心频率为2.4GHz，信号带宽为1.3G，基于带通采样定理，确定ADC采样频率为3.2GHz(n＝1)，FPGA工作时钟为200MHz。由此可得：
M＝fs/fclock＝16
同时依据上述参数生成3目标线性调频回波信号，目标中心参考距离分别为2000m，2025m，2045m。
对回波信号进行处理时，首先将其并行分解为16路，并对此16路线性调频回波信号进行多通道并行混频滤波，过程如下：
产生M路混频参考信号：
Figure RE-GDA0003745856920000121
同理：
Figure RE-GDA0003745856920000122
因此将16回波信号与上述16路混频参考信号一一对应相乘，在此过程中只需做一次乘法运算。混频之后的信号再通过多通道并行滤波处理，可设滤波器通带截止频率为650MHz，阻带截止频率为800MHz，阶数为32阶，且并行分解到16路，每路2阶。滤波器阶数，通带截止频率和阻带截止频率可根据信号带宽与滤波后要求的信号质量灵活设置。在经过混频滤波处理后，将16路零中频I信号与16路零中频Q信号与16路并行生成的去斜参考信号相乘，做复数乘法时，每路只需做两次乘法运算和一次加法运算即可完成去斜处理运算过程。
在去斜处理之后，对16路I信号和16路Q信号分别做16倍多相抽取滤波，滤波器阶数为64阶，并行分解到16路，每路4阶，同时滤波器通带截止频率和阻带截止频率仍为650MHz与800MHz。由此得到200M采样频率的经过去斜处理后的IQ两路信号。对此信号再进行FFT运算可以得到带有目标距离信息的单频点信号。信号处理过程如图5所示，图5(a)为回波信号FFT 处理后的结构，图5(b)为回波信号混频滤波后FFT的结果，图5(c)为回波信号去斜处理后FFT的结果，图5(d)为回波信号去斜处理后抽取滤波FFT 的结果。
至此，本实施例完成了对带有3目标距离信息的大带宽回波信号的去斜处理，理论计算结果与FPGA硬件运行结果保持一致。本实施例也表明了本发明为适应FGPA等硬件的时钟计算速率，将高速率大带宽回波信号的数字去斜处理转换为多通道并行处理，具有实现简单、资源利用率高、计算延时小等特点，在200MHz时钟频率下，上述去斜过程总体延时不超过100个时钟周期，即500ns，并且通过多相抽取滤波器降低信号的采样率。
以上仅是本发明的优选实施例，并非对本发明作任何形式上的限制。虽然本发明已以优选实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员，在不脱离本发明技术方案范围的情况下，都可利用上述揭示的技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

Claims (7)

1.一种宽带线性调频信号多通道并行去斜方法，其特征在于，所述方法包括以下步骤：
利用高采样率ADC器件对大带宽中频回波信号进行高质量带通采样，并根据FPGA时钟速率做并行分解成M路数字信号；
对M路回波信号进行多通道并行数字下变频和低通滤波处理，得到M路零中频IQ信号；
基于直接数字频率合成原理，分M路并行产生相同带宽、时宽且斜率相反的线性调频信号作为去斜参考信号；
将M路零中频回波信号与M路去斜参考信号相乘，然后经过M倍多相抽取滤波及FFT处理，得到去斜后采样率为fs/M且带有目标距离信息的单频点回波信号，完成对整个大带宽线性调频信号数字去斜处理的多通道并行实现。
2.根据权利要求1所述的宽带线性调频信号多通道并行去斜方法，其特征在于，基于带通采样定理，ADC采样频率满足下式：
Figure RE-FDA0003745856910000011
式中f0为信号中心频率，B为信号带宽，n取fs≥2B的最大正整数；
利用满足上式的ADC采样频率fs对回波信号进行中频采样得到数字回波信号x[n]，且将x[n]依据FPGA时钟工作频率fclock并行分解为M路回波信号，此时满足fs＝M*fclock，然后对M路信号进行数字下变频以及多通道并行低通滤波，所述数字下变频过程将M路以f0为中心频率的中频信号变成零中频基带信号，计算过程如下：
Figure RE-FDA0003745856910000012
式中cos(2πf0nTs)与sin(2πf0nTs)均由M路DDS并行产生；
经过数字下变频之后，得到零中频M路I信号和M路Q信号，将此IQ信号通过多通道并行低通滤波器进行滤波，从而滤除高频分量，得到零中频回波信号。
3.根据权利要求1所述的宽带线性调频信号多通道并行去斜方法，其特征在于，对于2路并行的4阶FIR滤波器，所述低通滤波器的滤波过程包括：
滤波器输出y[n]表示为：
Figure RE-FDA0003745856910000021
式中h(m)为滤波器系数，将上式分为偶数ye[n]与奇数yo[n]两部分之和：
Figure RE-FDA0003745856910000022
ye(n)＝h(0)x(n)+h(2)x(n-2)
yo(n)＝h(1)x(n-1)+h(3)x(n-3)
同理将y(n+1)分为：
Figure RE-FDA0003745856910000023
ye(n+1)＝h(0)x(n+1)+h(2)x(n-1)
yo(n+1)＝h(1)x(n)+h(3)x(n-2)
ye(n+1)与yo(n)输入数据相同，同时yo(n+1)与ye(n)输入数据相差两个时钟周期，因此通过FIR滤波器进行并行滤波。
4.根据权利要求3所述的宽带线性调频信号多通道并行去斜方法，其特征在于，所述FIR滤波器进行并行滤波包括：在n时刻，第一路滤波输出结果为第一路n时刻数据滤波结果的偶数部分和第二路n-1时刻数据滤波结果的奇数部分延时一个时钟周期后的和；第二路滤波输出结果为第二路n时刻数据滤波结果的偶数部分与第一路n时刻数据滤波结果的奇数部分之和，依次类推，将此滤波器并行实现结构扩展到n阶M路。
5.根据权利要求4所述的宽带线性调频信号多通道并行去斜方法，其特征在于，利用DDS并行产生M路大带宽数字去斜参考信号包括：
线性调频信号表达式如下：
Sr[n]＝Acos[2πf0nTs+πk(nTs)2]
式中n∈[0，N]，Ts＝1/fs，
n+1时刻的表达式如下：
Sr[n+1]＝Acos[2πf0(n+1)Ts+πk(n+1)2Ts 2]
当M路并行产生时，n＝Mi+j，其中
Figure RE-FDA0003745856910000031
j∈[0，M]，上述两式中相位部分变换为：
Figure RE-FDA0003745856910000032
Figure RE-FDA0003745856910000033
式中K＝kTs/2＝k/2fs，因此：
Figure RE-FDA0003745856910000034
其中，记
Figure RE-FDA0003745856910000035
ΔXk(n)＝K[M(Mi+j)+M(Mi+j+1)+M2]Ts，基于直接数字频率合成原理，得到M路不同的相位步进值，从而将大带宽线性调频去斜参考信号分M路并行产生。
6.根据权利要求5所述的宽带线性调频信号多通道并行去斜方法，其特征在于，将并行生成的M路线性调频去斜参考信号Irm-j*Qrm与经过处理后的M路零中频回波信号Im+j*Qm相乘，得到去斜后的Idm与Qdm两路回波信号，如下式所示：
(Im+jQm)·(Irm-jQrm)＝-Im·Qrm+Qm·Irm+j(Im·Irm+Qm·Qrm)，m＝0，1，…M
其中，Irm为M路线性调频去斜参考信号中的I路信号，Qrm为M路线性调频去斜参考信号的Q路信号，Im为M路零中频回波信号中的I路信号，Qm为M路零中频回波信号的Q路信号；
将上述去斜之后的两路回波信号经过M倍多相抽取滤波以及FFT变换之后，得到采样率为fs/M且带有目标距离等信息的单频点回波信号。
7.根据权利要求1所述的宽带线性调频信号多通道并行去斜方法，其特征在于，对回波信号进行处理时，首先将其并行分解为16路，并对此16路线性调频回波信号进行多通道并行混频滤波，过程如下：
产生M路混频参考信号：
Figure RE-FDA0003745856910000041
同理：
Figure RE-FDA0003745856910000042
因此将16回波信号与上述16路混频参考信号一一对应相乘，在此过程中只需做一次乘法运算；
混频之后的信号再通过多通道并行滤波处理，在经过混频滤波处理后，将16路零中频I信号与16路零中频Q信号与16路并行生成的去斜参考信号相乘，做复数乘法时，每路只需做两次乘法运算和一次加法运算即可完成去斜处理运算过程。
CN202210558439.0A 2022-05-20 2022-05-20 一种宽带线性调频信号多通道并行去斜方法 Pending CN114966564A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210558439.0A CN114966564A (zh) 2022-05-20 2022-05-20 一种宽带线性调频信号多通道并行去斜方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210558439.0A CN114966564A (zh) 2022-05-20 2022-05-20 一种宽带线性调频信号多通道并行去斜方法

Publications (1)

Publication Number Publication Date
CN114966564A true CN114966564A (zh) 2022-08-30

Family

ID=82985110

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210558439.0A Pending CN114966564A (zh) 2022-05-20 2022-05-20 一种宽带线性调频信号多通道并行去斜方法

Country Status (1)

Country Link
CN (1) CN114966564A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116910456A (zh) * 2023-09-13 2023-10-20 北京坤驰科技有限公司 一种滤波方法、装置、电子设备及计算机可读存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116910456A (zh) * 2023-09-13 2023-10-20 北京坤驰科技有限公司 一种滤波方法、装置、电子设备及计算机可读存储介质
CN116910456B (zh) * 2023-09-13 2023-12-01 北京坤驰科技有限公司 一种滤波方法、装置、电子设备及计算机可读存储介质

Similar Documents

Publication Publication Date Title
WO2014130297A1 (en) System and method for reconstruction of sparse frequency spectrum from ambiguous under-sampled time domain data
CN110210101B (zh) 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法
CN110927680B (zh) 基于数字去斜和频域均衡的宽带接收数字波束形成方法
CN114966564A (zh) 一种宽带线性调频信号多通道并行去斜方法
CN107748354B (zh) 基于分析与综合的宽带数字波束形成装置
CN109639277A (zh) 一种基于iserdes和并行fir滤波的高速信号预处理方法
CN112526458A (zh) 一种基于提取参数分数时延的宽带nlfm发射波束形成方法
CN201555948U (zh) 子带合成实现超大带宽信号产生与接收系统
CN107144821A (zh) 宽带数字阵雷达中基于时延波束形成的高效接收通道
Jiang et al. Design and realization of FPGA-based DRFM with high instantaneous bandwidth
Tang et al. Design and simulation of DRFM system based on digital channelized receiver and transmitter
Wang et al. Design and FPGA implementation of digital pulse compression for HF chirp radar based on modified orthogonal transformation
Agarwal et al. A factorization method for FPGA implementation of sample rate converter for a multi-standard radio communications
CN108696464B (zh) 一种iq与4通道tiadc联合失真盲估计与修正方法
Datta et al. Implementation of fractional sample rate digital down converter for radio receiver applications
Pang et al. Wide-band digital down conversion technique based on poly-phase filter
Xu et al. Design and implementation of digital channelized receiver in multi-FPGA
Liu et al. Broadband DDC based on polyphase filter and its FPGA implementation
Li et al. Compensation method for the CIC filter in digital down converter
Kumar et al. Reconfigurable band pass filter using kaiser window for satellite communication
Jiang et al. An improved wideband digital baseband conversion method
CN112187215B (zh) 一种级联半带插值滤波器结构
CN107565967A (zh) 基于周期信号混频的信号处理与采样方法
Liu et al. Design of digital IF receiver based on ADCs and FPGAs
Zhang et al. Jamming Structure of LFM Pulse Compression Radar Based on Digital Channelization

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