CN110967676A - 基于epc-mimo编解码设计的目标距离解模糊方法 - Google Patents

基于epc-mimo编解码设计的目标距离解模糊方法 Download PDF

Info

Publication number
CN110967676A
CN110967676A CN201911169336.XA CN201911169336A CN110967676A CN 110967676 A CN110967676 A CN 110967676A CN 201911169336 A CN201911169336 A CN 201911169336A CN 110967676 A CN110967676 A CN 110967676A
Authority
CN
China
Prior art keywords
mimo radar
array element
target
pulse
transmitted
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
CN201911169336.XA
Other languages
English (en)
Other versions
CN110967676B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201911169336.XA priority Critical patent/CN110967676B/zh
Publication of CN110967676A publication Critical patent/CN110967676A/zh
Application granted granted Critical
Publication of CN110967676B publication Critical patent/CN110967676B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/418Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/282Transmitters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/411Identification of targets based on measurements of radar reflectivity
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

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

本发明公开了一种基于EPC‑MIMO编解码设计的目标距离解模糊方法,主要用于解决无杂波环境下的雷达发射‑接收联合处理中的目标距离解模糊问题。本发明方法的步骤包括:构造编码向量;对基带信号编码;对编码后的基带信号进行上变频处理;对回波信号进行下变频处理;对模拟中频信号进行匹配滤波处理;构造解码向量;对匹配滤波信号解码;生成无距离模糊的目标回波信号。本发明充分利用了基于Fourier基函数的EPC‑MIMO雷达的解距离模糊能力,进一步拓展了编码形式,采用更灵活的编码设计从而提高波形捷变能力。

Description

基于EPC-MIMO编解码设计的目标距离解模糊方法
技术领域
本发明属于信号处理技术领域,更进一步涉及雷达信号处理技术领域中的一种基于阵元脉冲编码EPC(Element-Pulse Coding)多输入多输出MIMO(Multiple-InputMultiple-Output)编解码设计的目标距离解模糊方法。本发明可用于解决无杂波环境下的雷达发射-接收联合处理中的距离解模糊问题。
背景技术
在脉冲雷达中,如果采用较高的发射脉冲重复频率,同样两个目标,就很难测量出两个目标的真实距离。因为很难判定这些回波是哪一个发射脉冲的接收回波,这就是测距模糊,此距离解模糊问题会影响雷达捕获的目标的检测性能。
洪兴勇在其发表的论文“基于快速余差查表法的脉冲多普勒雷达解距离模糊算法”(微型机与应用,10.19358/j.issn.1674-7720.2017.04.025)中提出了一种采用快速余差查表的目标参数解模糊编码方法。该方法的具体步骤是:首先在每个CPI中,雷达目标距离的最大模糊次数是有限的,以每个重频CPI的最大不模糊距离的整数步数为步进建立解模糊余差距离表,这是该方法的关键;其次根据各个CPI点迹原始信息,填写解模糊距离表;再次读取解模糊距离表,判断是否满足解模糊准则,满足则输出目标点迹的幅度和频道;最后判断解距离模糊是否结束,以距离单元和多普勒频道来表示目标点迹信息。该方法虽然能有效解决PD雷达距离模糊问题,但是,该方法仍然存在的不足之处是,该方法在雷达发射信号与最大距离探测目标之间的距离不大时能够快速匹配出目标的真实距离,由于目标回波对应的可视距离形成的同余方程组的过程是一维空间到多维有限空间变换,随着雷达目标距离增大其空间复杂度会迅速增加,建立以及读取解模糊距离表的运算量会过大,从而解决目标距离模糊问题时花费的时间过长。
西安电子科技大学在其申请的专利文献“雷达阵元-脉冲编码与处理方法”(专利申请号:201910093293.5,申请公布号:CN 109946654 A)中公开了一种雷达阵元-脉冲编码与处理方法。该方法的具体步骤是:首先通过在发射阵元和发射脉冲之间加载相位编码设计,对不同阵元、不同脉冲的基带信号的初始相位进行编码调制,得到多通道发射信号;其次在接收端采用相应的解码处理;最后实现不同发射脉冲对应的回波信号的区分,实现目标距离参数解模糊。但是,该方法仍然存在的不足之处是,该方法构造编码向量时并没有关注编码参数的设置方法,从而导致该方法在解决无杂波环境下的目标距离解模糊问题时编码形式单一,不具有灵活性。
发明内容
本发明的目的是针对上述现有技术的不足,提出了一种基于EPC-MIMO编解码设计的目标距离解模糊方法,用于解决无杂波环境下的雷达发射-接收联合处理中的目标距离解模糊问题。
实现本发明目的的思路是,在发射阵元和发射脉冲之间加载相位编码设计,以便区分不同脉冲照射到目标上收到的回波信号。利用阵元脉冲编码EPC加权相位系数,对MIMO雷达中每个发射阵元待发射的基带信号进行相位编码标记。虽然接收维的Fourier功率谱实现聚焦,但由于加入编码序列,发射维的Fourier功率谱散焦。根据设计的编码形式在接收端进行解码处理,最后通过波束形成使得频谱集中在联合-接收空间域,实现对目标的距离解模糊。从而使得雷达的量测值准确反映目标的实际距离,提高目标的检测跟踪性能。
本发明的具体步骤如下:
步骤1,按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的编码向量:
Figure BDA0002288283440000021
其中,ck表示MIMO雷达中发射阵元待发射的第k个脉冲对应的基带信号的编码向量,k表示发射阵元待发射的脉冲的序号,k=1,2,…K,K表示MIMO雷达中一次相干处理时间内的脉冲数,e(·)表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,γ1表示MIMO雷达中第1个发射阵元的基带信号的加权相位系数的编码参数,γm表示表示MIMO雷达中第m个发射阵元的基带信号的加权相位系数的编码参数,
Figure BDA0002288283440000022
M表示MIMO雷达的发射阵元总数,∈表示属于符号,z表示整数编码,1≤z≤M-1,m表示发射阵元的序号,T表示转置操作;
步骤2,对基带信号编码:
(2a)按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的加权相位系数:
Figure BDA0002288283440000031
其中,cm,k表示MIMO雷达中第m个发射阵元待发射的第k个脉冲对应的基带信号的加权相位系数;
(2b)利用构造的加权相位系数,对MIMO雷达中每个发射阵元待发射的基带信号进行相位编码得到MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号;
步骤3,对基带信号进行上变频处理:
利用混频器,对MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号进行上变频处理,得到MIMO雷达中每个发射阵元待发射的脉冲信号;
步骤4,对回波信号进行下变频处理:
利用混频器,对MIMO雷达中每个发射阵元发射的脉冲信号的回波信号进行下变频处理,得到MIMO雷达中每个接收阵元对应的模拟中频信号;
步骤5,对模拟中频信号进行匹配滤波处理:
利用匹配滤波器,对MIMO雷达中每个接收阵元对应的模拟中频信号进行匹配滤波处理,得到每个模拟中频信号对应的匹配滤波信号;
步骤6,按照下式,构造MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量:
Figure BDA0002288283440000032
其中,
Figure BDA0002288283440000033
表示MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量,ck-q+1表示MIMO雷达中发射阵元发射的第(k-q+1)个脉冲对应的基带信号的编码向量,q表示对应MIMO雷达捕获的目标所处的距离区的序号的解码参数,*表示共轭操作,⊙表示哈达马Hadamard积操作,
Figure BDA0002288283440000034
表示MIMO雷达中发射阵元发射的第(k-pl+1)个脉冲对应的基带信号的编码向量,pl表示MIMO雷达捕获的目标所处的距离区的序号,pl=1,2,…P,P表示MIMO雷达捕获的目标所处的最大距离区序号;
步骤7,对匹配滤波信号解码:
利用步骤6构造的解码向量,对每个模拟中频信号对应的匹配滤波信号进行解码,得到MIMO雷达捕获的目标所处的距离区的序号的解码参数q与MIMO雷达捕获的目标所处的距离区的序号pl相对应则表示解码成功,得到MIMO雷达中接收阵元对应的解码后的目标回波信号:
Figure BDA0002288283440000041
其中,
Figure BDA0002288283440000042
表示MIMO雷达中接收阵元对应的解码后的目标回波信号,
Figure BDA0002288283440000043
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的期望距离区的回波信号,所述的参数向量β0包含待检测目标所处的期望距离区、角度和归一化多普勒频率信息,
Figure BDA0002288283440000044
p0表示MIMO雷达捕获的待检测目标所处的距离区的序号,p0=1,2,…P,Ω0表示MIMO雷达捕获的目标相对于阵列法线方向的角度,
Figure BDA0002288283440000045
表示MIMO雷达捕获的目标的归一化多普勒频率,
Figure BDA0002288283440000046
Figure BDA0002288283440000047
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的第k个发射脉冲相应的期望距离区的回波信号,
Figure BDA0002288283440000048
ξ0表示MIMO雷达捕获的待检测目标对发射阵元发射的脉冲信号的反射系数,s(β0)表示MIMO雷达捕获的待检测目标的等效的联合发射-接收端的空时二维导向矢量,
Figure BDA0002288283440000049
st表示MIMO雷达捕获的待检测目标的等效的时间导向矢量,
Figure BDA00022882834400000410
表示克罗内克Kronecker积操作,b表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure BDA00022882834400000411
k表示MIMO雷达的波数矢量,pR,1表示MIMO雷达中第1个接收阵元位置,pR,n表示MIMO雷达中第n个接收阵元位置,n=1,2,…N,N表示MIMO雷达的接收阵元总数,
Figure BDA00022882834400000412
表示MIMO雷达捕获的待检测目标的等效的解码后的发射端导向矢量,
Figure BDA00022882834400000413
a表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure BDA00022882834400000414
pT,1表示MIMO雷达中第1个发射阵元位置,pT,m表示MIMO雷达中第m个发射阵元位置,m=1,2,…M,
Figure BDA0002288283440000051
表示求和操作,
Figure BDA0002288283440000052
表示与MIMO雷达捕获的干扰目标的参数向量βi有关的模糊距离区的回波信号,
Figure BDA0002288283440000053
表示高斯白噪声的矢量;
步骤8,生成无距离模糊的目标回波信号:
对解码后的目标回波信号进行波束形成处理,得到MIMO雷达中接收阵元对应的无距离模糊的目标回波信号。
本发明与现有技术相比具有以下优点:
第一,本发明构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的编码向量,克服了现有技术构造编码向量时并没有关注编码参数的设置方法,从而导致该方法在解决无杂波环境下的目标距离解模糊问题时编码形式单一,不具有灵活性问题。使得本发明在距离模糊信号理想分离的情况下,天线阵元的EPC序列中的元素可以任意排序,进一步拓展了二维相位编码形式,可以采用更灵活的编码设计从而提高波形捷变的能力,编码向量的形式灵活多样,增加了MIMO雷达发射信号的灵活性。
第二,本发明构造MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量,克服了现有技术由于目标回波对应的可视距离形成的同余方程组的过程是一维空间到多维有限空间变换,随着雷达捕获的目标距离增大其空间复杂度会迅速增加,建立以及读取解模糊距离表的运算量会过大,从而解决目标距离模糊问题时花费的时间过长。使得本发明基于编码向量的形式构造解码向量,实现对模糊信号的快速分离,同时增加了发射信号被截获的难度。
附图说明
图1是本发明的流程图;
图2是本发明仿真1中接收机对信号进行正确解码前的Fourier功率谱分布图;
图3是本发明仿真1中接收机对信号进行正确解码后的Fourier功率谱分布图;
图4是本发明仿真2中接收机对信号进行正确解码前的Fourier功率谱分布图;
图5是本发明仿真2中接收机对信号进行正确解码后的Fourier功率谱分布图。
具体实施方式
下面结合附图对本发明的效果作进一步详细描述。
参照附图1,对本发明的具体实现步骤作进一步详细描述。
步骤1,按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的编码向量:
Figure BDA0002288283440000061
其中,ck表示MIMO雷达中发射阵元待发射的第k个脉冲对应的基带信号的编码向量,k表示发射阵元待发射的脉冲的序号,k=1,2,…K,K表示MIMO雷达中一次相干处理时间内的脉冲数,e(·)表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,γ1表示MIMO雷达中第1个发射阵元的基带信号的加权相位系数的编码参数,γm表示表示MIMO雷达中第m个发射阵元的基带信号的加权相位系数的编码参数,
Figure BDA0002288283440000062
M表示MIMO雷达的发射阵元总数,∈表示属于符号,z表示整数编码,1≤z≤M-1,m表示发射阵元的序号,T表示转置操作;所述基带信号是指发射阵元的第m个满足正交波形条件的基带信号Φm(t),m=1,2,…M,M表示发射阵元数;所述正交波形条件是指每个发射信号的功率相同,互相之间没有影响,实现从不同的距离模糊区完全分离。
Figure BDA0002288283440000063
的生成过程如下:
根据对应两个不同距离模糊区的导向矢量的内积为零生成非线性方程组:
Figure BDA0002288283440000064
其中,
Figure BDA0002288283440000065
表示从m=1至m=M对(·)求和,P表示MIMO雷达捕获的目标所处的最大距离区序号。
按照z=1时的γm构建加权相位系数cm,k所得的形式是Fourier基,所述
Figure BDA0002288283440000066
是指上述非线性方程组基于Fourier基的解;所述
Figure BDA0002288283440000067
确定的编码序列中的编码参数可以任意排序,不一定按照Fourier基递增顺序;当给定z的取值,所述集合中包含的元素为M个,编码序列的解是M个元素的任意排列组合,因此,存在解的个数为M!。
步骤2,对基带信号编码。
按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的加权相位系数:
Figure BDA0002288283440000071
其中,cm,k表示MIMO雷达中第m个发射阵元待发射的第k个脉冲对应的基带信号的加权相位系数。
在编码参数确定的基础上,只要已知MIMO雷达中发射阵元的序号和该发射阵元待发射的脉冲的序号,就可以对基带信号进行相位编码。
利用构造的加权相位系数,对MIMO雷达中每个发射阵元待发射的基带信号进行相位编码得到MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号。
步骤3,对基带信号进行上变频处理。
利用混频器,对MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号进行上变频处理,得到MIMO雷达中发射阵元待发射的脉冲信号;所述的MIMO雷达中发射阵元待发射的脉冲信号如下:
Figure BDA0002288283440000072
其中,sm,k表示MIMO雷达中第m个发射阵元待发射的第k个脉冲信号,rect表示脉冲调制操作,t表示脉冲重复周期Tr的时间序号,t∈[0,Tr],Tp表示MIMO雷达中发射阵元待发射的脉冲信号的脉冲持续时间,Φm(t)表示MIMO雷达中发射阵元待发射的第m个脉冲对应的基带信号,f0表示第1个发射阵元待发射脉冲信号的载频,tk表示MIMO雷达中发射阵元待发射的脉冲信号的相干处理时间CPI内的第k个脉冲信号的时间序号,tk∈[0,W],W表示MIMO雷达中发射阵元待发射的脉冲信号的相干处理时间,W=K×Tr,在给定的相干处理时间内满足tk=(k-1)Tr+t。
步骤4,对回波信号进行下变频处理。
利用混频器,对MIMO雷达中每个发射阵元发射的脉冲信号的回波信号进行下变频处理,得到MIMO雷达中每个接收阵元对应的模拟中频信号;所述的MIMO雷达中每个接收阵元对应的模拟中频信号如下:
Figure BDA0002288283440000081
其中,Xn表示MIMO雷达中第n个接收阵元对应的模拟中频信号,∑表示求和操作,ξ0表示MIMO雷达捕获的待检测目标对发射阵元发射的脉冲信号的反射系数,R表示MIMO雷达捕获的目标相对于雷达的距离,c表示光速,pT,m表示MIMO雷达中第m个发射阵元位置,m=1,2,…M,pR,n表示MIMO雷达中第n个接收阵元位置,n=1,2,…N,N表示MIMO雷达的接收阵元总数。
步骤5,对模拟中频信号进行匹配滤波处理。
利用匹配滤波器,对MIMO雷达中每个接收阵元对应的模拟中频信号进行匹配滤波处理,得到每个模拟中频信号对应的匹配滤波信号;所述的每个模拟中频信号对应的匹配滤波信号如下:
Figure BDA0002288283440000082
其中,xkl)表示与MIMO雷达捕获的目标的参数向量βl有关的第k个发射脉冲相应的模拟中频信号对应的匹配滤波信号,所述的参数向量βl包含距离区、角度和归一化多普勒频率信息,
Figure BDA0002288283440000083
Ω0表示MIMO雷达捕获的目标相对于阵列法线方向的角度,
Figure BDA0002288283440000084
表示MIMO雷达捕获的目标的归一化多普勒频率,b表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure BDA0002288283440000085
k表示MIMO雷达的波数矢量,pR,1表示MIMO雷达中第1个接收阵元位置,
Figure BDA0002288283440000086
表示克罗内克Kronecker积操作,
Figure BDA0002288283440000087
表示MIMO雷达中发射阵元发射的第(k-pl+1)个脉冲对应的基带信号的编码向量,⊙表示哈达马Hadamard积操作,a表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure BDA0002288283440000088
pT,1表示MIMO雷达中第1个发射阵元位置。
步骤6,按照下式,构造MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量:
Figure BDA0002288283440000091
其中,
Figure BDA0002288283440000092
表示MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量,ck-q+1表示MIMO雷达中发射阵元发射的第(k-q+1)个脉冲对应的基带信号的编码向量,q表示对应MIMO雷达捕获的目标所处的距离区的序号的解码参数,*表示共轭操作,
Figure BDA0002288283440000093
表示MIMO雷达中发射阵元发射的第(k-pl+1)个脉冲对应的基带信号的编码向量,pl表示MIMO雷达捕获的目标所处的距离区的序号,pl=1,2,…P;所述的解码向量与MIMO雷达中发射阵元发射脉冲的序号无关。
步骤7,对匹配滤波信号解码。
利用步骤6构造的解码向量,对每个模拟中频信号对应的匹配滤波信号进行解码,得到MIMO雷达捕获的目标所处的距离区的序号的解码参数q与MIMO雷达捕获的目标所处的距离区的序号pl相对应则表示解码成功,得到MIMO雷达中接收阵元对应的解码后的目标回波信号:
Figure BDA0002288283440000094
其中,
Figure BDA0002288283440000095
表示MIMO雷达中接收阵元对应的解码后的目标回波信号,
Figure BDA0002288283440000096
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的期望距离区的回波信号,所述的参数向量β0包含待检测目标所处的期望距离区、角度和归一化多普勒频率信息,
Figure BDA0002288283440000097
p0表示MIMO雷达捕获的待检测目标所处的距离区的序号,p0=1,2,…P,
Figure BDA0002288283440000098
Figure BDA0002288283440000099
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的第k个发射脉冲相应的期望距离区的回波信号,
Figure BDA00022882834400000910
s(β0)表示MIMO雷达捕获的待检测目标的等效的联合发射-接收端的空时二维导向矢量,
Figure BDA00022882834400000911
st表示MIMO雷达捕获的待检测目标的等效的时间导向矢量,
Figure BDA00022882834400000912
表示MIMO雷达捕获的待检测目标的等效的解码后的发射端导向矢量,
Figure BDA0002288283440000101
Figure BDA0002288283440000102
表示求和操作,
Figure BDA0002288283440000103
表示与MIMO雷达捕获的干扰目标的参数向量βi有关的模糊距离区的回波信号,
Figure BDA0002288283440000104
表示高斯白噪声的矢量。
步骤8,生成无距离模糊的目标回波信号。
对解码后的目标回波信号进行波束形成处理,得到MIMO雷达中接收阵元对应的无距离模糊的目标回波信号;所述的波束形成处理是按照下式完成的:
Figure BDA0002288283440000105
其中,y表示MIMO雷达中接收阵元对应的无距离模糊的目标回波信号,w表示MIMO雷达中期望距离区的目标回波信号的权矢量,w=s(β0),H表示共轭转置操作。
综上,传统的脉冲相位调制方法中的基带信号仅与发射阵元有关,脉冲之间没有差异;而本发明的相位编码在不同的发射阵元,以及同一发射阵元的不同脉冲之间均有差异,实现对不同脉冲发射的信号的区分。
下面结合仿真对本发明做进一步的描述。
1.仿真实验条件:
本发明的仿真实验的硬件平台为:Intel(R)Core(TM)i5-8265U CPU@1.60GHz,频率为1.8GHz,Nvidia GeForce MX250。
本发明的仿真实验的软件使用matlab2016b。
本发明的仿真参数:
考虑一个基于地基L波段(f0=3GHz)EPC-MIMO雷达,其发射天线单元数为M=10,接收单元数N=10。为简单起见,发射和接收阵列都是均匀的线性阵列,且阵列间距取半波长。相干脉冲数为K=16,脉冲重复频率为20kHz,于是最大不模糊距离为7.5km。雷达的作用距离是30km,因此最大距离模糊是P=4。仿真时,假定解码参数为q=1。上述仿真参数如表1所示:
表1仿真参数一览表
Figure BDA0002288283440000106
2.仿真内容与结果分析:
本发明仿真实验是采用本发明的方法分别对阵元脉冲编码多输入多输出EPC-MIMO雷达在联合发射-接收空域中单个目标和多个目标在不同脉冲下的Fourier功率谱分布情况进行了仿真,获得接收机对信号进行正确解码前后的Fourier功率谱分布图。功率谱分布图中白色部分表示MIMO雷达捕获的目标对应的联合发射-接收二维Fourier功率谱,功率谱分布图中的横坐标轴表示归一化后的接收空间频率,纵坐标轴表示归一化后的发射空间频率。
仿真1,在上述表1的仿真参数下,采用本发明的方法对接收机对信号进行正确解码前后,对阵元脉冲编码多输入多输出EPC-MIMO雷达在联合发射-接收空域中单个目标在不同脉冲下的Fourier功率谱分布情况进行了仿真,结果分别如图2和图3所示。
图2是本发明仿真1中接收机对信号进行正确解码前的Fourier功率谱分布图,其中,图2(a)表示接收机对信号进行正确解码前,发射阵元发射第1个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图2(b)表示接收机对信号进行正确解码前,发射阵元发射第2个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图2(c)表示接收机对信号进行正确解码前,发射阵元发射第3个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图2(d)表示接收机对信号进行正确解码前,发射阵元发射第4个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。
图3是本发明仿真1中接收机对信号进行正确解码后的Fourier功率谱分布图,其中,图3(a)表示接收机对信号进行正确解码后,发射阵元发射第1个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图3(b)表示接收机对信号进行正确解码后,发射阵元发射第2个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图3(c)表示接收机对信号进行正确解码后,发射阵元发射第3个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。图3(d)表示接收机对信号进行正确解码后,发射阵元发射第4个脉冲时MIMO雷达捕获的单个目标对应的联合发射-接收二维Fourier功率谱分布图。
由图2(a)可以看出,发射阵元发射第1个脉冲时发射维Fourier功率谱聚焦在[-0.2,0.2]的范围内。但由于本发明加入编码序列,由图2(b)、图2(c)和图2(d)可以看出,发射阵元发射不同脉冲时发射维的Fourier功率谱散焦。由此可以看出,由于EPC-MIMO雷达发射机存在相位调制,联合发射-接收空域中目标的Fourier功率谱会随着编码过程的脉冲而变化。
由图3可以看出,接收机对信号进行正确解码后,发射阵元发射不同脉冲时,目标的Fourier功率谱集中在联合发射-接收空间域。由此可见,本发明可以解决在无杂波环境下的目标距离解模糊问题。
仿真2,在上述表1的仿真参数下,采用本发明技术,接收机对信号进行正确解码前后,对阵元脉冲编码多输入多输出EPC-MIMO雷达在联合发射-接收空域中多个目标在不同脉冲下的Fourier功率谱分布情况进行了仿真,结果分别如图4和图5所示。
图4是本发明仿真2中接收机对信号进行正确解码前的Fourier功率谱分布图,其中,图4(a)表示接收机对信号进行正确解码前,发射阵元发射第1个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图4(b)表示接收机对信号进行正确解码前,发射阵元发射第2个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图4(c)表示接收机对信号进行正确解码前,发射阵元发射第3个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图4(d)表示接收机对信号进行正确解码前,发射阵元发射第4个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。
图5是本发明仿真2中接收机对信号进行正确解码后的Fourier功率谱分布图,其中,图5(a)表示接收机对信号进行正确解码后,发射阵元发射第1个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图5(b)表示接收机对信号进行正确解码后,发射阵元发射第2个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图5(c)表示接收机对信号进行正确解码后,发射阵元发射第3个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。图5(d)表示接收机对信号进行正确解码后,发射阵元发射第4个脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱分布图。
由图4(a)可以看出,发射阵元发射第1个脉冲时发射维Fourier功率谱聚焦在[-0.2,0.2]和[0.4,0.6]的范围内,对应的接收维Fourier功率谱聚焦在[-0.2,0.2]和[0.4,0.6]的范围内。由于本发明加入编码序列,由图4(b)、图4(c)和图4(d)可以看出,发射阵元发射不同脉冲时MIMO雷达捕获的多个目标对应的联合发射-接收二维Fourier功率谱散焦。显然,联合发射-接收空域中多个目标的Fourier功率谱在解码之前会随着编码过程的脉冲而变化。
由图5(a)和图5(b)可以看出,发射阵元发射第1个脉冲和第2个脉冲时发射维Fourier功率谱聚焦在[-0.2,0.2]和[0.4,0.6]的范围内,对应的接收维Fourier功率谱聚焦在[-0.2,0.2]和[0.4,0.6]的范围内。由图5(c)和图5(d)可以看出,发射阵元发射第3个脉冲和第4个脉冲时发射维Fourier功率谱聚焦在[-0.2,0.2]、[0.2,0.4][0.6,0.8]和[0.4,0.6]的范围内,对应的接收维Fourier功率谱聚焦在[-0.2,0.2]、[-0.2,0.2]、[-0.2,0.2]和[0.4,0.6]的范围内。显然,接收机对信号进行正确解码后,期望距离区内的多个目标聚焦,而其他区域的目标的功率谱散焦且被削弱。当天线方向图主瓣对准期望信号的波达方向时,本发明可以自动削弱距离模糊的目标,有效解决了雷达捕获的目标距离解模糊问题。

Claims (5)

1.一种基于EPC-MIMO编解码设计的目标距离解模糊方法,其特征在于,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的编码向量,构造MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量;该方法的具体步骤包括如下:
步骤1,按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的编码向量:
Figure FDA0002288283430000011
其中,ck表示MIMO雷达中发射阵元待发射的第k个脉冲对应的基带信号的编码向量,k表示发射阵元待发射的脉冲的序号,k=1,2,…K,K表示MIMO雷达中一次相干处理时间内的脉冲数,e(·)表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,γ1表示MIMO雷达中第1个发射阵元的基带信号的加权相位系数的编码参数,γm表示表示MIMO雷达中第m个发射阵元的基带信号的加权相位系数的编码参数,
Figure FDA0002288283430000012
M表示MIMO雷达的发射阵元总数,∈表示属于符号,z表示整数编码,1≤z≤M-1,m表示发射阵元的序号,T表示转置操作;
步骤2,对基带信号编码:
(2a)按照下式,构造MIMO雷达中发射阵元待发射的每个脉冲对应的基带信号的加权相位系数:
Figure FDA0002288283430000013
其中,cm,k表示MIMO雷达中第m个发射阵元待发射的第k个脉冲对应的基带信号的加权相位系数;
(2b)利用构造的加权相位系数,对MIMO雷达中每个发射阵元待发射的基带信号进行相位编码,得到MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号;
步骤3,对基带信号进行上变频处理:
利用混频器,对MIMO雷达中每个发射阵元待发射的进行相位编码后的基带信号进行上变频处理,得到MIMO雷达中发射阵元待发射的脉冲信号;
步骤4,对回波信号进行下变频处理:
利用混频器,对MIMO雷达中每个发射阵元发射的脉冲信号的回波信号进行下变频处理,得到MIMO雷达中每个接收阵元对应的模拟中频信号;
步骤5,对模拟中频信号进行匹配滤波处理:
利用匹配滤波器,对MIMO雷达中每个接收阵元对应的模拟中频信号进行匹配滤波处理,得到每个模拟中频信号对应的匹配滤波信号;
步骤6,按照下式,构造MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量:
Figure FDA0002288283430000021
其中,
Figure FDA0002288283430000022
表示MIMO雷达中发射阵元发射的每个脉冲对应的基带信号的编码向量的解码向量,pl表示MIMO雷达捕获的目标所处的距离区的序号,pl=1,2,…P,P表示MIMO雷达捕获的目标所处的最大距离区序号,q表示对应MIMO雷达捕获的目标所处的距离区的序号的解码参数;
步骤7,对匹配滤波信号解码:
利用步骤6构造的解码向量,对每个模拟中频信号对应的匹配滤波信号进行解码,得到MIMO雷达捕获的目标所处的距离区的序号的解码参数q与MIMO雷达捕获的目标所处的距离区的序号pl相对应则表示解码成功,得到MIMO雷达中接收阵元对应的解码后的目标回波信号:
Figure FDA0002288283430000023
其中,
Figure FDA0002288283430000024
表示MIMO雷达中接收阵元对应的解码后的目标回波信号,
Figure FDA0002288283430000025
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的期望距离区的回波信号,所述的参数向量β0包含待检测目标所处的期望距离区、角度和归一化多普勒频率信息,
Figure FDA0002288283430000031
p0表示MIMO雷达捕获的待检测目标所处的距离区的序号,p0=1,2,…P,Ω0表示MIMO雷达捕获的目标相对于阵列法线方向的角度,
Figure FDA0002288283430000032
表示MIMO雷达捕获的目标的归一化多普勒频率,
Figure FDA0002288283430000033
Figure FDA0002288283430000034
表示与MIMO雷达捕获的待检测目标的参数向量β0有关的第k个发射脉冲相应的期望距离区的回波信号,
Figure FDA0002288283430000035
ξ0表示MIMO雷达捕获的待检测目标对发射阵元发射的脉冲信号的反射系数,s(β0)表示MIMO雷达捕获的待检测目标的等效的联合发射-接收端的空时二维导向矢量,
Figure FDA0002288283430000036
st表示MIMO雷达捕获的待检测目标的等效的时间导向矢量,
Figure FDA0002288283430000037
表示克罗内克Kronecker积操作,b表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure FDA0002288283430000038
k表示MIMO雷达的波数矢量,pR,1表示MIMO雷达中第1个接收阵元位置,pR,n表示MIMO雷达中第n个接收阵元位置,n=1,2,…N,N表示MIMO雷达的接收阵元总数,
Figure FDA0002288283430000039
表示MIMO雷达捕获的待检测目标的等效的解码后的发射端导向矢量,
Figure FDA00022882834300000310
a表示MIMO雷达捕获的目标的等效的接收端导向矢量,
Figure FDA00022882834300000311
pT,1表示MIMO雷达中第1个发射阵元位置,pT,m表示MIMO雷达中第m个发射阵元位置,m=1,2,…M,
Figure FDA00022882834300000312
表示求和操作,
Figure FDA00022882834300000313
表示与MIMO雷达捕获的干扰目标的参数向量βi有关的模糊距离区的回波信号,
Figure FDA00022882834300000314
表示高斯白噪声的矢量;
步骤8,生成无距离模糊的目标回波信号:
对解码后的目标回波信号进行波束形成处理,得到MIMO雷达中接收阵元对应的无距离模糊的目标回波信号。
2.根据权利要求1所述的基于EPC-MIMO编解码设计的目标距离解模糊方法,其特征在于,步骤3中所述的MIMO雷达中发射阵元待发射的脉冲信号如下:
Figure FDA0002288283430000041
其中,sm,k表示MIMO雷达中第m个发射阵元待发射的第k个脉冲信号,rect表示脉冲调制操作,t表示脉冲重复周期Tr的时间序号,t∈[0,Tr],Tp表示MIMO雷达中发射阵元待发射的脉冲信号的脉冲持续时间,Φm(t)表示MIMO雷达中发射阵元待发射的第m个脉冲对应的基带信号,f0表示第1个发射阵元待发射脉冲信号的载频,tk表示MIMO雷达中发射阵元待发射的脉冲信号的相干处理时间CPI内的第k个脉冲信号的时间序号,tk∈[0,W],W表示MIMO雷达中发射阵元待发射的脉冲信号的相干处理时间,W=K×Tr,在给定的相干处理时间内满足tk=(k-1)Tr+t。
3.根据权利要求2所述的基于EPC-MIMO编解码设计的目标距离解模糊方法,其特征在于,步骤4中所述的MIMO雷达中每个接收阵元对应的模拟中频信号如下:
Figure FDA0002288283430000042
其中,Xn表示MIMO雷达中第n个接收阵元对应的模拟中频信号,∑表示求和操作,R表示MIMO雷达捕获的目标相对于雷达的距离,c表示光速。
4.根据权利要求3所述的基于EPC-MIMO编解码设计的目标距离解模糊方法,其特征在于,步骤5中所述的每个模拟中频信号对应的匹配滤波信号如下:
Figure FDA0002288283430000043
其中,xkl)表示与MIMO雷达捕获的目标的参数向量βl有关的第k个发射脉冲相应的模拟中频信号对应的匹配滤波信号,所述的参数向量βl包含距离区、角度和归一化多普勒频率信息,
Figure FDA0002288283430000044
Figure FDA0002288283430000045
表示MIMO雷达中发射阵元发射的第(k-pl+1)个脉冲对应的基带信号的编码向量,⊙表示哈达马Hadamard积操作。
5.根据权利要求4所述的基于EPC-MIMO编解码设计的目标距离解模糊方法,其特征在于,步骤8中所述的波束形成处理是按照下式完成的:
Figure FDA0002288283430000051
其中,y表示MIMO雷达中接收阵元对应的无距离模糊的目标回波信号,w表示MIMO雷达中期望距离区的目标回波信号的权矢量,w=s(β0),H表示共轭转置操作。
CN201911169336.XA 2019-11-26 2019-11-26 基于epc-mimo编解码设计的目标距离解模糊方法 Active CN110967676B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911169336.XA CN110967676B (zh) 2019-11-26 2019-11-26 基于epc-mimo编解码设计的目标距离解模糊方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911169336.XA CN110967676B (zh) 2019-11-26 2019-11-26 基于epc-mimo编解码设计的目标距离解模糊方法

Publications (2)

Publication Number Publication Date
CN110967676A true CN110967676A (zh) 2020-04-07
CN110967676B CN110967676B (zh) 2023-03-10

Family

ID=70031590

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911169336.XA Active CN110967676B (zh) 2019-11-26 2019-11-26 基于epc-mimo编解码设计的目标距离解模糊方法

Country Status (1)

Country Link
CN (1) CN110967676B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111505600A (zh) * 2020-05-19 2020-08-07 西北大学 基于stpc的fda-mimo雷达信号处理方法、装置及介质
CN114442048A (zh) * 2021-12-23 2022-05-06 西安电子科技大学 Epc-mimo雷达接收信号处理的fpga实现方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8044846B1 (en) * 2007-11-29 2011-10-25 Lockheed Martin Corporation Method for deblurring radar range-doppler images
WO2018045566A1 (zh) * 2016-09-09 2018-03-15 深圳大学 基于压缩感知的随机脉冲多普勒雷达角度-多普勒成像方法
CN107966688A (zh) * 2017-11-09 2018-04-27 东南大学 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN109946654A (zh) * 2019-01-30 2019-06-28 西安电子科技大学 雷达阵元-脉冲编码与处理方法
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8044846B1 (en) * 2007-11-29 2011-10-25 Lockheed Martin Corporation Method for deblurring radar range-doppler images
WO2018045566A1 (zh) * 2016-09-09 2018-03-15 深圳大学 基于压缩感知的随机脉冲多普勒雷达角度-多普勒成像方法
CN107966688A (zh) * 2017-11-09 2018-04-27 东南大学 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN109946654A (zh) * 2019-01-30 2019-06-28 西安电子科技大学 雷达阵元-脉冲编码与处理方法
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王伟伟等: "基于频率分集阵列的机载雷达距离模糊杂波抑制方法", 《电子与信息学报》 *
王成浩等: "FDA-SAR高分辨宽测绘带成像距离解模糊方法", 《电子学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111505600A (zh) * 2020-05-19 2020-08-07 西北大学 基于stpc的fda-mimo雷达信号处理方法、装置及介质
CN114442048A (zh) * 2021-12-23 2022-05-06 西安电子科技大学 Epc-mimo雷达接收信号处理的fpga实现方法

Also Published As

Publication number Publication date
CN110967676B (zh) 2023-03-10

Similar Documents

Publication Publication Date Title
US11582305B2 (en) Vehicle radar system with a shared radar and communication system
JP6278961B2 (ja) レーダ装置
CN105874351B (zh) 编码孔径雷达(car)信号的处理方法及装置
US10024958B2 (en) Radar apparatus
CN109946654B (zh) 雷达阵元-脉冲编码与处理方法
JP2019527835A5 (zh)
CN110412568B (zh) 基于扩展方位角相位编码的距离模糊杂波抑制方法
JP2016035440A (ja) レーダ装置および物体検知方法
WO2013141924A9 (en) Coded aperture beam analysis method and apparatus
EP2798370A2 (en) Coded aperture beam analysis method and apparatus
CN113777577B (zh) 基于mimo雷达的目标检测方法、装置及电子设备
US20220308160A1 (en) Automotive radar device
NO155416B (no) Aktivt deteksjonssystem som anvender samtidige multitransmisjoner.
Chahrour et al. Hybrid beamforming for interference mitigation in MIMO radar
CN110471033A (zh) 基于相干fda的二维空时编码干扰抑制方法
JP2012181109A (ja) レーダ装置
CN110967676B (zh) 基于epc-mimo编解码设计的目标距离解模糊方法
US10200081B2 (en) Systems and methods for signal detection and digital bandwidth reduction in digital phased arrays
CN111856410B (zh) 恒模约束下机载mimo雷达的稳健发射接收联合优化方法
CN112612013B (zh) Fda-mimo雷达增量距离-角度二维波束形成方法
CN113093137A (zh) 一种基于fda-mimo雷达最优频偏估计的杂波抑制方法
CN114609595A (zh) 一种频分正交mimo雷达信号处理方法
Landi et al. Time-orthogonal-waveform-space-time adaptive processing for distributed aperture radars
Cattenoz et al. Adaptive processing for MIMO radar realistic non perfectly orthogonal waveforms
Weiner et al. High-throughput communications using constant-modulus waveforms with mitigation of range-sidelobe modulation

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