CN110837081B - 基于p/d波段雷达信号融合处理的高速目标检测方法 - Google Patents

基于p/d波段雷达信号融合处理的高速目标检测方法 Download PDF

Info

Publication number
CN110837081B
CN110837081B CN201911089281.1A CN201911089281A CN110837081B CN 110837081 B CN110837081 B CN 110837081B CN 201911089281 A CN201911089281 A CN 201911089281A CN 110837081 B CN110837081 B CN 110837081B
Authority
CN
China
Prior art keywords
band
target
radar
waveform
processing
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
Application number
CN201911089281.1A
Other languages
English (en)
Other versions
CN110837081A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201911089281.1A priority Critical patent/CN110837081B/zh
Publication of CN110837081A publication Critical patent/CN110837081A/zh
Application granted granted Critical
Publication of CN110837081B publication Critical patent/CN110837081B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/06—Systems determining position data of a target
    • G01S13/08—Systems for measuring distance only
    • G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
    • G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
    • G01S13/341—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal wherein the rate of change of the transmitted frequency is adjusted to give a beat of predetermined constant frequency, e.g. by adjusting the amplitude or frequency of the frequency-modulating signal
    • 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
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S13/584—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
    • 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/292—Extracting wanted echo-signals
    • G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
    • 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
    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

本发明公开了一种基于P/D波段雷达信号融合处理的高速目标检测方法，包括以下步骤：选取P波段工作雷达的发射波形形式为伪随机二进制相移键控脉冲信号，D波段工作雷达的发射波形形式为线性调频连续波信号；确定P波段雷达的伪随机二进制相移键控脉冲发射波形参数和D波段雷达的线性调频连续波发射波形参数；对P波段和D波段对应的工作波形进行协同发射和接收；对P波段和D波段雷达的发射信号回波进行处理，将两种波形的信号处理结果进行融合处理，得到目标的距离和速度信息。本发明解决了高速运动目标在D波段雷达下易发生测速模糊的问题，提高了目标测距与测速的精度。

Description

基于P/D波段雷达信号融合处理的高速目标检测方法
技术领域
本发明属于雷达信息技术领域，特别是一种基于P/D波段雷达信号融合处理的高速目标检测方法。
背景技术
近年来隐身技术的迅速发展对战略和战术防御系统提出了严峻挑战,迫使人们考虑如何摧毁隐身兵器并研究反隐身技术。现役隐身飞机效果最佳的频段在1～40GHz之间，如果雷达工作频率超出此范围，就可以显著提高雷达探测能力，从而使隐身飞机的隐身能力大大降低。
P波段频率范围为0.23～1GHz，D波段频率范围为110～170GHz，两者都处于隐身频段外，因此P波段与D波段雷达都具有出色的反隐身能力。P波段选用二进制相移键控(Binary Phase Shift Keying,BPSK)工作波形，D波段选用线性调频连续波(LinearFrequency Modulated Continuous Wave,LFMCW)工作波形，两者都有较高的测距精度，但在测速性能方面，P波段雷达发射信号波长大，因此存在测速精度差的问题。D波段LFMCW雷达测速精度高，但其工作波形波长小，容易导致因多普勒模糊造成的测速模糊现象。
发明内容
本发明的目的在于提供一种基于P/D波段雷达信号融合处理的高速目标检测方法，以提高目标距离、速度测量精度，同时达到反隐身的效果。
实现本发明目的的技术解决方案为：一种基于P/D波段雷达信号融合处理的高速目标检测方法，包括以下步骤：
步骤1，选取P波段雷达的发射波形形式为伪随机二进制相移键控脉冲信号，D波段雷达的发射波形形式为线性调频连续波波形；
步骤2，根据最大测距范围、最大测速范围以及测距精度确定P波段雷达的发射波形参数，包括码元宽度、脉冲宽度、脉冲重复周期；确定D波段雷达的发射波形参数，包括扫频带宽、扫频周期；
步骤3，对P波段和D波段雷达对应的工作波形进行协同发射和接收；
步骤4，对P波段雷达的BPSK发射信号回波进行处理：发射信号回波经AD采样和数字下变频后得到数字信号，对采样得到的数字信号数据进行相关处理，将结果保存至存储器RAM，在连续发射了N个脉冲完成相关处理后，将N个脉冲相关处理结果中同一个距离单元结果进行动目标检测处理，获得距离-多普勒二维平面，进行恒虚警检测得到目标峰值点，得到目标在P波段BPSK工作波形下的距离和速度信息；
步骤5，对D波段雷达的LFMCW发射信号回波进行处理：发射信号回波与本振信号进行混频得到差拍信号，对采样得到的差拍信号数据进行加窗快速傅里叶变换处理得到第一维FFT结果，经过一个相参积累时间，即N个脉冲重复周期后，进行第二维的FFT，MTD完成后形成距离-多普勒二维平面，对MTD处理结果进行CFAR检测处理，获得目标点迹的距离和多普勒信息，得到目标在D波段LFMCW工作波形下的距离和速度信息；
步骤6，对P波段和D波段雷达回波信号处理结果进行融合处理：将P波段和D波段探测结果在距离维进行匹配处理，得到同一目标在两种工作模式下的距离和速度测量结果；对同一目标，使用P波段速度测量结果判断目标在D波段多普勒频率测量模糊倍数，进行解模糊处理，得到目标距离和速度信息。
进一步的，步骤2确定P波段雷达的BPSK发射波形参数，根据最大测距范围、最大测速范围、测距精度确定BPSK发射波形参数，所述BPSK发射波形参数包括码元宽度、脉冲宽度、脉冲重复周期；确定D波段雷达的LFMCW发射波形参数，根据最大测距范围、测距精度确定LFMCW发射波形参数，所述LFMCW发射波形参数包括扫频带宽、扫频周期，具体如下：
首先设置P波段雷达BPSK发射波形参数如下：波形的脉冲宽度th,P，脉冲重复周期Tr,P，码元宽度τP；
码元宽度设置与天线带宽相匹配，对应发射信号3dB带宽为1/τP；
脉冲重复周期Tr,P对应雷达的最大作用距离Rmax：
Figure BDA0002266374630000021
其中c为光速；
发射波形带宽为1/τP，对应雷达距离分辨力为：
Figure BDA0002266374630000031
目标最大相对运动速度为vmax，对应最长积累时间：
Figure BDA0002266374630000032
设置相干处理间隔TCPI,P不能超过最长积累时间，对应多普勒分辨力为：
Figure BDA0002266374630000033
对应速度分辨力：
Figure BDA0002266374630000034
其中f0为载波频率；
设置D波段雷达LFMCW发射波形参数如下：波形的扫频带宽BD，波形的扫频周期Tr,D；
在LFMCW发射波的扫频段中，根据系统要求的测距分辨力和测距精度ΔR确定扫频带宽，即
Figure BDA0002266374630000035
波形扫频周期Tr,D对应雷达的最大作用距离Rmax：
Figure BDA0002266374630000036
进一步的，步骤3对P波段和D波段雷达对应的工作波形进行协同发射和接收：两种工作波形在对应的相参积累时间内同时发射，P波段雷达在发射脉冲结束后接收目标回波，D波段雷达在LFMCW发射信号扫频周期的接收有效期间进行目标差拍回波的采样和处理。
进一步的，步骤4所述的对P波段雷达的BPSK发射信号回波进行相关处理，具体如下：
BPSK回波信号经AD采样和数字下变频后得到数字基带信号，数字基带信号是发射码字经过延迟调制到多普勒频率上的调幅信号，目标回波信号变成数字信号，然后将数字基带信号送入各个距离门，和参考码进行相乘累加，码字匹配的距离门上得到最大的输出；各个距离门上的参考码是由发射码字依次经过延迟一个码元时间得到的；每个距离门依次相差一个码元时间，故距离分辩力为cτP/2；这样经过一次相关积累处理，将长度为P的一组数据压缩成一个数据，完成一次积累，因此相关积累也叫脉冲压缩；
接着对下一组长度为P的数据进行压缩，依次反复，将相关处理结果送入存储器RAM，直到处理完所有数据；
最后，按距离门对脉压后的数据做FFT,获得距离-多普勒二维平面，然后进行恒虚警检测CFAR处理，获得在P波段BPSK波形下目标的距离和速度信息；
对相关器输出同一距离门下的数据Yr(1:N,m)作N点FFT处理，得到对应的多普勒维数据Fr(n,m),即：
Fr(1:N,m)＝fft[Yr(1:N,m)] (8)
其中Fr(n,m)表示回波通过MTD处理后第m个距离单元的第n个多普勒单元的结果。
进一步的，步骤5所述的对D波段雷达的LFMCW发射信号回波进行处理，具体如下：
LFMCW发射信号回波与本振信号进行混频得到差拍信号，对回波差拍信号进行AD采样；
目标最大距离为Rmax，对应回波时延最大值为：
Figure BDA0002266374630000041
对一个扫频周期后τmax～Tr期间的回波差拍信号进行采样，进行第一维FFT变换按距离门分别存入对应存储器RAM，经过一个相干处理间隔CPI后，对每个距离门存储数据作第二维FFT变换，即MTD处理后形成距离-多普勒二维平面，对MTD处理结果进行CFAR检测处理，获得目标在D波段LFMCW工作波形下的距离和多普勒信息；
通过第一维FFT处理测得目标差拍频率为：
Figure BDA0002266374630000042
其中R为目标距离，fd为目标多普勒频率。
通过MTD处理测得目标模糊多普勒频率记为fd,D，带入上式得到目标距离为：
Figure BDA0002266374630000051
其中μ＝BD/Tr,D，μ为D波段LFMCW波形的调频斜率。
进一步的，步骤6所述对P波段和D波段雷达回波信号处理结果融合处理，具体如下：
将P波段与D波段探测结果在距离维上进行匹配，获得属于同一目标在两种工作模式下的距离和速度测量结果；
对于同一个目标，使用P波段速度测量结果vP,判断目标在D波段雷达下多普勒频率测量模糊倍数ND，如下式所示：
Figure BDA0002266374630000052
Figure BDA0002266374630000053
其中vmax,D为D波段LFMCW波形的最大不模糊探测速度，λD为D波段LFMCW波形的波长；
根据已求的D波段多普勒频率模糊倍数，进行解模糊处理，得到D波段下目标的不模糊多普勒频率fd：
fd＝fd,D+ND/Tr,D (14)
其中Tr,D为D波段LFMCW波形的扫频周期，由此得到目标精确的速度为：
Figure BDA0002266374630000054
将D波段LFMCW目标差拍频率fB与不模糊多普勒频率fd相结合，获得目标精确的距离为：
Figure BDA0002266374630000055
其中μ＝BD/Tr,D，μ为D波段LFMCW波形的调频斜率。
本发明与现有技术相比，其显著优点为：(1)雷达工作波形处于隐身涂料工作频段外，可以达到频域反隐身的目的；(2)解决了高速目标在D波段LFMCW工作模式下易发生测速模糊的问题，以及高速目标回波的多普勒频移会降低P波段相位编码信号在雷达接收机的压缩效果的现象；(3)目标配对方式简单，运算量低且测距、测速精度高，原理简单，即使使用低成本的器件也能实现很好的实时性，在性能与成本方面具有竞争力。
附图说明
图1是本发明基于P/D波段雷达信号融合处理的高速目标检测方法流程图.
图2是本发明P波段和D波段工作波形协同发射示意图。
图3是本发明P波段回波信号相关处理的示意图。
图4是本发明P波段回波信号处理的流程图。
图5是本发明D波段LFMCW回波信号处理的流程图。
图6是本发明P波段和D波段雷达信号融合处理的流程图。
图7是本发明P波段信号处理、D波段信号处理与P波段和D波段雷达信号融合处理的检测精度对比图，其中图7(a)为测距误差对比示意图，图7(b)为测速误差对比示意图。
具体实施方式
本发明提供一种基于P/D波段雷达信号融合处理的高速目标检测方法，主要包括P/D波段波形结构设计、波形参数设计和信号处理设计三部分，如图1所示，该方法具体步骤为：
步骤一：设计P/D波段发射波形形式：选取P波段雷达的发射波形形式为伪随机二进制相移键控(Binary Phase Shift Keying,BPSK)脉冲信号，D波段雷达的发射波形形式为线性调频连续波(LinearFrequency ModulatedContinuous Wave,LFMCW)波形；
步骤二：设计P/D波段雷达发射波形参数：确定P波段雷达的发射波形参数，根据最大测距范围、最大测速范围以及测距精度确定波形的参数，所述波形参数包括码元宽度、脉冲宽度、脉冲重复周期；确定D波段雷达的发射波形参数，根据最大测距范围、测距精度确定波形参数，所述波形参数包括扫频带宽、扫频周期；P/D波段雷达发射波形参数设置应保证在最大测距范围内不发生测距模糊，以便于融合处理过程中的目标匹配。
首先设置P波段雷达BPSK发射波形参数如下：波形的脉冲宽度th,P，脉冲重复周期Tr,P，码元宽度τP；
码元宽度设置与天线带宽相匹配，对应发射信号3dB带宽为1/τP；
脉冲重复周期Tr,P对应雷达的最大作用距离Rmax：
Figure BDA0002266374630000071
其中c为光速；
发射波形带宽为1/τP，对应雷达距离分辨力为：
Figure BDA0002266374630000072
目标最大相对运动速度为vmax，对应最长积累时间：
Figure BDA0002266374630000073
设置相干处理间隔TCPI,P不能超过最长积累时间，对应多普勒分辨力为：
Figure BDA0002266374630000074
对应速度分辨力：
Figure BDA0002266374630000075
其中c为光速，f0为载波频率；
设置D波段雷达LFMCW发射波形参数如下：波形的扫频带宽BD，波形的扫频周期Tr,D；
在LFMCW发射波的扫频段中，根据系统要求的测距分辨力与测距精度ΔR确定扫频带宽，即
Figure BDA0002266374630000081
其中c为光速；
波形扫频周期Tr,D对应雷达的最大作用距离Rmax：
Figure BDA0002266374630000082
其中c为光速；
对P波段和D波段雷达对应的工作波形进行协同发射和接收，如图2所示，两种工作波形在对应的相参积累时间内同时发射，P波段雷达在发射脉冲结束后接收目标回波，D波段雷达在LFMCW发射信号扫频周期的接收有效期间进行目标差拍回波的采样和处理。
步骤三：设计P/D波段信号处理方式：如图3、图4所示，对P波段雷达的BPSK发射信号回波进行相关处理：BPSK回波信号经AD采样和数字下变频后得到数字基带信号，数字基带信号实际上就是发射码字经过延迟调制到多普勒频率上的调幅信号，目标回波信号变成数字信号，然后把数字基带信号送入各个距离门，和参考码进行相乘累加，码字匹配的距离门上就能得到最大的输出。各个距离门上的参考码是由发射码字依次经过延迟一个码元时间得到的。每个距离门依次相差一个码元时间，故距离分辩力为cτP/2。这样经过一次相关积累处理，就能把长度为P的一组数据压缩成一个数据,完成一次积累，因此相关积累也叫脉冲压缩。接着对下一组长度为P的数据进行压缩，依次反复，将相关处理结果送入存储器RAM，直到处理完所有数据。最后,按距离门对脉压后的数据做FFT,获得距离-多普勒二维平面，然后进行恒虚警检测CFAR处理，获得在P波段BPSK波形下目标的距离和速度信息。
对相关器输出同一距离门下的数据Yr(1:N,m)作N点FFT处理，得到对应的多普勒维数据Fr(n,m),即：
Fr(1:N,m)＝fft[Yr(1:N,m)] (8)
其中Fr(n,m)表示回波通过动目标检测(MovingTarget Detection,MTD)处理后第m个距离单元的第n个多普勒单元的结果。
对D波段雷达的LFMCW发射信号回波进行处理，如图5所示：
LFMCW发射信号回波与本振信号进行混频得到差拍信号，对回波差拍信号进行AD采样。
目标最大距离为Rmax，对应回波时延最大值为：
Figure BDA0002266374630000091
为了使回波差拍信号在采样期间连续，对一个扫频周期后τmax～Tr期间的回波差拍信号进行采样，进行第一维FFT变换按距离门分别存入对应存储器RAM,经过一个相干处理间隔CPI后，对每个距离门存储数据作第二维FFT变换，即MTD处理后形成距离-多普勒二维平面，对MTD处理结果进行CFAR检测处理，获得目标在D波段LFMCW工作波形下的距离和多普勒信息。
通过第一维FFT处理测得目标差拍频率为：
Figure BDA0002266374630000092
其中R为目标距离，fd为目标多普勒频率。
通过MTD处理测得目标模糊多普勒频率记为fd,D，带入上式得到目标距离为：
Figure BDA0002266374630000093
其中μ＝BD/Tr,D，μ为D波段LFMCW波形的调频斜率。
对P波段和D波段雷达回波信号处理结果融合处理，如图6所示，具体如下：
第1步：将P波段与D波段探测结果在距离维上进行匹配，获得属于同一目标在两种工作模式下的距离和速度测量结果；
第2步：对于同一个目标，使用P波段速度测量结果vP,判断在D波段多普勒频率测量模糊倍数ND，如下式所示：
Figure BDA0002266374630000094
Figure BDA0002266374630000095
其中vmax,D为D波段LFMCW波形的最大不模糊探测速度，λD为D波段LFMCW波形的波长；
第3步：根据已求的D波段多普勒频率模糊倍数，进行解模糊处理，得到D波段不模糊多普勒频率fd：
fd＝fd,D+ND/Tr,D (14)
其中Tr,D为D波段LFMCW波形的扫频周期，由此得到目标精确的速度为：
Figure BDA0002266374630000101
其中λD为D波段雷达LFMCW工作波形的波长。
第4步：将D波段LFMCW目标差拍频率与不模糊多普勒频率相结合，获得目标精确的距离为：
Figure BDA0002266374630000102
其中μ＝BD/Tr,D，μ为D波段LFMCW波形的调频斜率。
下面结合实施例对本发明进行详细说明。
实施例
本发明的效果可以通过以下仿真实验和外场测试进一步说明。
在仿真实验的指标设计中，测距范围为0到150m，测速范围为150m/s到3000m/s，雷达距离分辨力不大于1m，测距精度优于5m/s。根据具体实施的设计方法，P波段雷达BPSK发射波形参数设计如表1所示：
表1
Figure BDA0002266374630000103
D波段雷达LFMCW发射波形参数设计如表2所示：
表2
Figure BDA0002266374630000111
上述方案中，P波段雷达相干处理间隔CPI取512个脉冲，D波段雷达相干处理间隔CPI取256个调制周期。
为了评估本发明处理方法的性能，我们在相同指标的条件下比较P/D波段下信号处理和本发明融合处理的测距精度和测速精度性能。
图7(a)为Matlab仿真为P波段BPSK工作模式、D波段LFMCW工作模式与融合处理工作模式下的测距精度与目标距离关系曲线。在同一目标速度条件下，对目标距离从1m到150m，间隔1m的所有情况下进行仿真，将检测结果与实际情况相减计算误差，通过比较误差矩阵的标准差来确定两种体制的测量精度的好坏。
图7(b)为Matlab仿真为P波段BPSK工作模式、D波段LFMCW工作模式与融合处理工作模式下的测速精度与目标速度关系曲线。在同一目标距离条件下，对目标速度从0m/s到3000m/s，间隔10m/s的所有情况下进行仿真，将检测结果与实际情况相减计算误差，通过比较误差矩阵的标准差来确定两种体制的测量精度的好坏。
从图7(a)中可以看出，P波段BPSK工作模式与D波段LFMCW工作模式下都具有较好的测距精度，因此可以用于融合处理过程中的目标匹配。从图7(b)中可以看出，在速度小于D波段最大不模糊探测速度时，D波段LFMCW工作模式具有不错的测速精度，但目标速度超过最大不模糊探测速度时，D波段LFMCW工作模式的测速精度由于测速模糊而发生急剧恶化，P波段BPSK工作模式没有发生测速模糊现象，但具有较大测速误差。本发明通过融合处理，目标在整个动态范围内都具有很好的测速精度，测速精度优于5m/s。

Claims (5)

1.一种基于P/D波段雷达信号融合处理的高速目标检测方法，其特征在于，包括以下步骤：
步骤1，选取P波段雷达的发射波形形式为伪随机二进制相移键控脉冲信号，D波段雷达的发射波形形式为线性调频连续波波形；
步骤2，根据最大测距范围、最大测速范围以及测距精度确定P波段雷达的发射波形参数，包括码元宽度、脉冲宽度、脉冲重复周期；确定D波段雷达的发射波形参数，包括扫频带宽、扫频周期；
步骤3，对P波段和D波段雷达对应的工作波形进行协同发射和接收；
步骤4，对P波段雷达的BPSK发射信号回波进行处理：发射信号回波经AD采样和数字下变频后得到数字信号，对采样得到的数字信号数据进行相关处理，将结果保存至存储器RAM，在连续发射了N个脉冲完成相关处理后，将N个脉冲相关处理结果中同一个距离单元结果进行动目标检测处理，获得距离-多普勒二维平面，进行恒虚警检测得到目标峰值点，得到目标在P波段BPSK工作波形下的距离和速度信息；
步骤5，对D波段雷达的LFMCW发射信号回波进行处理：发射信号回波与本振信号进行混频得到差拍信号，对采样得到的差拍信号数据进行加窗快速傅里叶变换处理得到第一维FFT结果，经过一个相参积累时间，即N个脉冲重复周期后，进行第二维的FFT，MTD完成后形成距离-多普勒二维平面，对MTD处理结果进行CFAR检测处理，获得目标点迹的距离和多普勒信息，得到目标在D波段LFMCW工作波形下的距离和速度信息；
步骤6，对P波段和D波段雷达回波信号处理结果进行融合处理：将P波段和D波段探测结果在距离维进行匹配处理，得到同一目标在两种工作模式下的距离和速度测量结果；对同一目标，使用P波段速度测量结果判断目标在D波段多普勒频率测量模糊倍数，进行解模糊处理，得到目标距离和速度信息；具体如下：
将P波段与D波段探测结果在距离维上进行匹配，获得属于同一目标在两种工作模式下的距离和速度测量结果；
对于同一个目标，使用P波段速度测量结果vP,判断目标在D波段雷达下多普勒频率测量模糊倍数ND，如下式所示：
Figure FDA0003692629670000021
Figure FDA0003692629670000022
其中vmax,D为D波段LFMCW波形的最大不模糊探测速度，λD为D波段LFMCW波形的波长；
根据已求的D波段多普勒频率模糊倍数，进行解模糊处理，得到D波段下目标的不模糊多普勒频率fd：
fd＝fd,D+ND/Tr,D (14)
其中Tr,D为D波段LFMCW波形的扫频周期，fd,D为通过MTD处理测得目标的模糊多普勒频率，由此得到目标精确的速度为：
Figure FDA0003692629670000023
将D波段LFMCW目标差拍频率fB与不模糊多普勒频率fd相结合，获得目标精确的距离为：
Figure FDA0003692629670000024
其中μ＝BD/Tr,D，其中BD为波形的扫频带宽，μ为D波段LFMCW波形的调频斜率。
2.根据权利要求1所述的基于P/D波段雷达信号融合处理的高速目标检测方法，其特征在于，步骤2确定P波段雷达的BPSK发射波形参数，根据最大测距范围、最大测速范围、测距精度确定BPSK发射波形参数，所述BPSK发射波形参数包括码元宽度、脉冲宽度、脉冲重复周期；确定D波段雷达的LFMCW发射波形参数，根据最大测距范围、测距精度确定LFMCW发射波形参数，所述LFMCW发射波形参数包括扫频带宽、扫频周期，具体如下：
首先设置P波段雷达BPSK发射波形参数如下：波形的脉冲宽度th,P，脉冲重复周期Tr,P，码元宽度τP；
码元宽度设置与天线带宽相匹配，对应发射信号3dB带宽为1/τP；
脉冲重复周期Tr,P对应雷达的最大作用距离Rmax：
Figure FDA0003692629670000031
其中c为光速；
发射波形带宽为1/τP，对应雷达距离分辨力为：
Figure FDA0003692629670000032
目标最大相对运动速度为vmax，对应最长积累时间：
Figure FDA0003692629670000033
设置相干处理间隔TCPI,P不能超过最长积累时间，对应多普勒分辨力为：
Figure FDA0003692629670000034
对应速度分辨力：
Figure FDA0003692629670000035
其中f0为载波频率；
设置D波段雷达LFMCW发射波形参数如下：波形的扫频带宽BD，波形的扫频周期Tr,D；
在LFMCW发射波的扫频段中，根据系统要求的测距分辨力和测距精度ΔR确定扫频带宽，即
Figure FDA0003692629670000036
波形扫频周期Tr,D对应雷达的最大作用距离Rmax：
Figure FDA0003692629670000037
3.根据权利要求1所述的基于P/D波段雷达信号融合处理的高速目标检测方法，其特征在于，步骤3对P波段和D波段雷达对应的工作波形进行协同发射和接收：两种工作波形在对应的相参积累时间内同时发射，P波段雷达在发射脉冲结束后接收目标回波，D波段雷达在LFMCW发射信号扫频周期的接收有效期间进行目标差拍回波的采样和处理。
4.根据权利要求1所述的基于P/D波段雷达信号融合处理的高速目标检测方法，其特征在于，步骤4所述的对P波段雷达的BPSK发射信号回波进行相关处理，具体如下：
BPSK回波信号经AD采样和数字下变频后得到数字基带信号，数字基带信号是发射码字经过延迟调制到多普勒频率上的调幅信号，目标回波信号变成数字信号，然后将数字基带信号送入各个距离门，和参考码进行相乘累加，码字匹配的距离门上得到最大的输出；各个距离门上的参考码是由发射码字依次经过延迟一个码元时间得到的；每个距离门依次相差一个码元时间，故距离分辩力为cτP/2；这样经过一次相关积累处理，将长度为P的一组数据压缩成一个数据，完成一次积累，因此相关积累也叫脉冲压缩；
接着对下一组长度为P的数据进行压缩，依次反复，将相关处理结果送入存储器RAM，直到处理完所有数据；
最后，按距离门对脉压后的数据做FFT,获得距离-多普勒二维平面，然后进行恒虚警检测CFAR处理，获得在P波段BPSK波形下目标的距离和速度信息；
对相关器输出同一距离门下的数据Yr(1:N,m)作N点FFT处理，得到对应的多普勒维数据Fr(n,m),即：
Fr(1:N,m)＝fft[Yr(1:N,m)] (8)
其中Fr(n,m)表示回波通过MTD处理后第m个距离单元的第n个多普勒单元的结果。
5.根据权利要求1所述的基于P/D波段雷达信号融合处理的高速目标检测方法，其特征在于，步骤5所述的对D波段雷达的LFMCW发射信号回波进行处理，具体如下：
LFMCW发射信号回波与本振信号进行混频得到差拍信号，对回波差拍信号进行AD采样；
目标最大距离为Rmax，对应回波时延最大值为：
Figure FDA0003692629670000051
对一个扫频周期后τmax～Tr期间的回波差拍信号进行采样，进行第一维FFT变换按距离门分别存入对应存储器RAM，经过一个相干处理间隔CPI后，对每个距离门存储数据作第二维FFT变换，即MTD处理后形成距离-多普勒二维平面，对MTD处理结果进行CFAR检测处理，获得目标在D波段LFMCW工作波形下的距离和多普勒信息；
通过第一维FFT处理测得目标差拍频率为：
Figure FDA0003692629670000052
其中R为目标距离，fd为目标多普勒频率；
将目标模糊多普勒频率fd,D带入上式得到目标距离为：
Figure FDA0003692629670000053
其中μ＝BD/Tr,D，μ为D波段LFMCW波形的调频斜率。
CN201911089281.1A 2019-11-08 2019-11-08 基于p/d波段雷达信号融合处理的高速目标检测方法 Active CN110837081B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911089281.1A CN110837081B (zh) 2019-11-08 2019-11-08 基于p/d波段雷达信号融合处理的高速目标检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911089281.1A CN110837081B (zh) 2019-11-08 2019-11-08 基于p/d波段雷达信号融合处理的高速目标检测方法

Publications (2)

Publication Number Publication Date
CN110837081A CN110837081A (zh) 2020-02-25
CN110837081B true CN110837081B (zh) 2022-09-06

Family

ID=69574750

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911089281.1A Active CN110837081B (zh) 2019-11-08 2019-11-08 基于p/d波段雷达信号融合处理的高速目标检测方法

Country Status (1)

Country Link
CN (1) CN110837081B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308437B (zh) * 2020-02-27 2022-03-22 南京慧尔视智能科技有限公司 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN111722218B (zh) * 2020-06-03 2023-04-07 武汉大学 一种双频复合波形高频雷达系统
CN112630740B (zh) * 2020-12-08 2022-10-14 中国人民解放军陆军工程大学 密集多目标的分辨方法及装置
CN112285684B (zh) * 2020-12-29 2021-03-16 南昌大学 相对运动速度跟踪方法、装置、可读存储介质及电子设备
CN113030893A (zh) * 2021-03-01 2021-06-25 成都安智杰科技有限公司 一种随机采样方法及装置
CN113608205B (zh) * 2021-06-25 2024-04-09 西安空间无线电技术研究所 一种超高精度地外天体微波着陆雷达测距测速方法
CN113267751A (zh) * 2021-06-29 2021-08-17 珠海上富电技股份有限公司 一种车载毫米波雷达抗干扰的方法
CN114415167A (zh) * 2021-12-23 2022-04-29 深圳成谷科技有限公司 速度检测方法、装置、电子设备及存储介质
CN115963467B (zh) * 2023-03-16 2023-06-06 福思(杭州)智能科技有限公司 调频波形参数的处理方法、装置和计算机设备

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8159385B2 (en) * 2010-02-04 2012-04-17 Sensis Corporation Conductive line communication apparatus and conductive line radar system and method
CN103675783A (zh) * 2012-09-25 2014-03-26 中国航天科工集团第二研究院二〇七所 一种宽带多波段成像的相参处理方法
ITPI20120110A1 (it) * 2012-11-05 2014-05-06 I D S Ingegneria Dei Sistemi S P A Metodo e dispositivo per per determinare la traiettoria di un proiettile emesso da una piccola arma da fuoco, in particolare per localizzare una postazione di sparo
CN106990403B (zh) * 2017-04-28 2019-08-02 西安电子科技大学 基于多波段两级信息融合的低空目标跟踪方法
CN107886121A (zh) * 2017-11-03 2018-04-06 北京清瑞维航技术发展有限公司 基于多波段雷达的目标识别方法、装置及系统
CN110208804A (zh) * 2019-06-02 2019-09-06 南京理工大学 适用于汽车防撞雷达的发射波形设计方法及目标探测算法

Also Published As

Publication number Publication date
CN110837081A (zh) 2020-02-25

Similar Documents

Publication Publication Date Title
CN110837081B (zh) 基于p/d波段雷达信号融合处理的高速目标检测方法
US5808580A (en) Radar/sonar system concept for extended range-doppler coverage
CN106405541B (zh) 全相参连续波多普勒雷达及其测距测速方法
CN107688178A (zh) 一种基于77GHz毫米波雷达的锯齿波测距测速方法
CN102680975B (zh) 一种天基步进频分时测角雷达空间非合作目标成像方法
CN109613506B (zh) 一种随机跳频重频捷变雷达目标回波信号检测方法
Schroeder et al. X-band FMCW radar system with variable chirp duration
AU712338B2 (en) Radar/sonar system concept for extended range-doppler coverage
CN111580063B (zh) 基于广义解调频-楔形变换的雷达目标检测方法
JP5460290B2 (ja) レーダ装置
Zheng et al. Radar detection and motion parameters estimation of maneuvering target based on the extended keystone transform (July 2018)
CN116087942A (zh) 航空高度表调制信号的产生方法
CN102508212A (zh) 一种线性调频信号多通道扩展脉压技术
CN109581350B (zh) 基于时频积分插值的雷达测距测速方法及装置
CN111693981A (zh) 变重复周期的二次相关处理方法
CN116500620A (zh) 毫米波雷达的数据处理方法、装置、存储介质及无人车
Malanowski et al. Noise vs. deterministic waveform radar—Possibilities and limitations
CN115561728A (zh) 一种fmcw雷达的单帧速度解模糊方法和装置
Taylor Ultra wideband radar
CN112083405B (zh) 一种基于混合波形的目标检测方法及相关装置
CN110082748B (zh) 一种无源雷达目标检测方法及装置
CN113820678A (zh) 一种外辐射源雷达高重频跳频信号相参积累方法
JP2013113723A (ja) レーダ装置
Vignesh et al. Design of less-detectable RADAR waveforms using stepped frequency modulation and coding
CN113030886B (zh) 一种高速目标距离徙动校正方法

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