CN113204017B - 一种基于lfmcw的mimo雷达信号处理方法 - Google Patents

一种基于lfmcw的mimo雷达信号处理方法 Download PDF

Info

Publication number
CN113204017B
CN113204017B CN202110280931.1A CN202110280931A CN113204017B CN 113204017 B CN113204017 B CN 113204017B CN 202110280931 A CN202110280931 A CN 202110280931A CN 113204017 B CN113204017 B CN 113204017B
Authority
CN
China
Prior art keywords
signal
scanning
lfmcw
mimo radar
beat
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
CN202110280931.1A
Other languages
English (en)
Other versions
CN113204017A (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 CN202110280931.1A priority Critical patent/CN113204017B/zh
Publication of CN113204017A publication Critical patent/CN113204017A/zh
Application granted granted Critical
Publication of CN113204017B publication Critical patent/CN113204017B/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
    • G01S13/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/08Systems for measuring distance only
    • G01S13/32Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
    • G01S13/34Systems 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/341Systems 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

Abstract

本发明公开了一种基于LFMCW的MIMO雷达信号处理方法,包括:建立包括M个发射阵元和N个接收阵元的MIMO雷达;利用MIMO雷达发射一组正交频分LFMCW信号;根据正交频分LFMCW信号的回波信号,获得有效差拍信号;对有效差拍信号进行接收波束形成,获得接收波束形成输出;获得雷达距离扫描范围内的K个距离扫描点;根据K个距离扫描点获得一组扫描有效差拍信号;将接收波束形成输出与一组扫描有效差拍信号进行互相关,得到时空匹配滤波后的信号。本发明提出一种基于LFMCW的MIMO雷达,并给出其信号处理实现方法,利用扫描差拍信号对接收信号进行距离扫描,得到时空匹配滤波结果,实现雷达目标检测。

Description

一种基于LFMCW的MIMO雷达信号处理方法
技术领域
本发明属于雷达技术领域,具体涉及一种基于LFMCW的MIMO雷达信号处理方法,可用于基于LFMCW的MIMO雷达目标检测。
背景技术
多输入多输出(Multiple Input Multiple Output,MIMO)雷达系统采用多天线发射分集的信号波形,多天线接收目标回波信号的体制。MIMO雷达可以利用相位或频率编码等方法实现发射信号波形的分集,在波形设计上非常灵活。与传统雷达相比,MIMO雷达在抗杂波、抗干扰、低截获、角度分辨率等性能方面具有明显优势。
线性调频连续波(Line Frequency Modulated Continuous Wave,LFMCW)雷达的频率调制随时间线性变化,具有峰值功率低、测距分辨率高、无距离盲区、体积小、结构简单成本低等优点,近年来应用较为广泛。
现有技术已经报道了一种利用时间分集的所谓MIMO线性调频连续波雷达,即对发射信号在时间轴上进行区分,相同的发射信号在不同的时间周期内进行发射,这种方式容易造成资源浪费,且不是真正的波形分集MIMO雷达。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于LFMCW的MIMO雷达信号处理方法。本发明要解决的技术问题通过以下技术方案实现:
本发明提供了一种基于LFMCW的MIMO雷达信号处理方法,包括:
S1:建立包括M个发射阵元和N个接收阵元的MIMO雷达;
S2:利用所述MIMO雷达发射一组正交频分LFMCW信号;
S3:根据所述正交频分LFMCW信号的回波信号,获得有效差拍信号;
S4:对所述有效差拍信号进行接收波束形成,获得接收波束形成输出;
S5:获得雷达距离扫描范围内的K个距离扫描点;
S6:根据所述K个距离扫描点获得一组扫描有效差拍信号;
S7:将所述接收波束形成输出与所述一组扫描有效差拍信号进行互相关,得到时空匹配滤波后的信号。
在本发明的一个实施例中,所述发射天线与所述接收天线均为等距线阵,所述M个发射阵元发射的信号为:
s(t)=[s1(t),s2(t),...si(t),...sM(t)]T,i=1,2,...M,
其中,si(t)为第i个发射阵元发射的信号,(·)T表示矩阵的转置。
在本发明的一个实施例中,所述S2包括:
根据雷达发射阵列的M个发射阵元,得到每个发射阵元产生的正交频分LFMCW信号:
Figure BDA0002978734320000021
其中,si(t)为第i个阵元产生的信号波形,δi为第i个阵元发射信号的初相,
Figure BDA0002978734320000022
为调频斜率,B为单个信号的带宽,T为信号时宽,t表示0到T内的时间,fi=f0+ciΔf,f0为中心载频,ci为频率编码,
Figure BDA0002978734320000023
为步进频率,Te为有效时宽,j表示复数符号。
在本发明的一个实施例中,所述S3包括:
S31:将所述N个接收阵元接收的回波信号Y(t)与第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到N个差拍信号sc(t):
sc(t)=[sc1(t),sc2(t),...scq(t),...scN(t)]T,q=1,2,...N,
Y(t)=[y1(t),y2(t),...yq(t),...yN(t)]T,q=1,2,...N,
其中,
Figure BDA0002978734320000024
(·)*表示矩阵的共轭,δ1为第一个阵元发射信号的初相,f1=f0+c1Δf,c1为频率编码;
S32:对所述差拍信号sc(t)中的每一个差拍信号scq(t)进行等间隔采样,获得采样后的差拍信号sc(n):
sc(n)=[sc1(n),sc2(n),...scq(n),...scN(n)]T,q=1,2,...N,
其中,
Figure BDA0002978734320000025
L=fST为总采样点数,n为第n个采样点,采样时间为信号时宽T,fS为采样频率;
S33:将所述采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
其中,sceq(n)=w(n)scq(n),
Figure BDA0002978734320000031
L′=<0.1L>,<·>为取整函数。
在本发明的一个实施例中,所述S4包括:
根据有效差拍信号sce(n),得到接收波束形成输出:
u(n)=bT0)sce(n),
其中,θ0为形成的接收波束的指向,b(θ0)为接收导向矢量,b(θ0)=[1,exp(j2πdrsinθ0/λ),...exp(j2πdr(N-1)sinθ0/λ)]T,λ为信号波长,dr为接收阵元的间距。
在本发明的一个实施例中,所述S5包括:
根据所述MIMO雷达设定的参数和采样频率,确定雷达距离扫描范围和距离间隔,得到K个距离扫描点,其中,所述K个距离扫描点与所述MIMO雷达之间的距离分别为R1,R2,...Rk,...RK,K≥Le
在本发明的一个实施例中,所述S6包括:
S61:获得第k个距离扫描点的时延:
Figure BDA0002978734320000032
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离;
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT0)s(t-τk)
s(t-τk)=[s1(t-τk),s2(t-τk),...si(t-τk),...sM(t-τk)]T,i=1,2,...M
其中,a(θ0)=[1,exp(j2πdt sinθ0/λ),...exp(j2πdt(M-1)sinθ0/λ)]T,dt为发射阵元间距,s(t-τk)表示第k个扫描点带时延的LFMCW信号,si(t-τk)表示第i个发射阵元第k个扫描点带时延的LFMCW信号:
Figure BDA0002978734320000041
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到第k个扫描点的差拍信号hck(t):
Figure BDA0002978734320000042
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
Figure BDA0002978734320000043
其中,采样时间为信号时宽T,采样频率fS,L=fST为总采样点数,n为第n个采样点;
S65:对采样后第k个扫描差拍信号hck(n)乘以矩形窗函数w(n),得到第k个扫描有效差拍信号hcek(n):
hcek(n)=w(n)hck(n),
其中,n=L′,L′+1,...,L-1;
S66:重复步骤S61至S65,得到K个距离扫描点的一组扫描有效差拍信号:
hce1(n),hce2(n),...hcek(n),...hceK(n),k=1,2,...,K。
在本发明的一个实施例中,所述S7包括:
根据所述接收波束形成输出u(n)与所述一组扫描有效差拍信号hce1(n),hce2(n),...hcek(n),...hceK(n)进行互相关,得到时空匹配滤波后的信号p:
p=[p1,p2,...pk,...pK]T,k=1,2,...,K
其中,
Figure BDA0002978734320000051
与现有技术相比,本发明的有益效果在于:
1、本发明基于LFMCW的MIMO雷达信号处理方法,利用频率编码,提出一种发射阵元同时发射一组正交LFMCW信号的MIMO雷达,并给出其信号处理实现方法,利用扫描差拍信号对接收信号进行距离扫描,得到时空匹配滤波结果,实现雷达目标检测。
2、本发明基于LFMCW的MIMO雷达信号处理方法,在抗杂波、抗干扰、低截获、峰值功率低、无距离盲区、距离及角度分辨率等方面具有明显优势。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1为本发明实施例提供的一种基于LFMCW的MIMO雷达信号处理方法的流程示意图;
图2为发射信号频率顺序排列且初相为零时的信号处理全局距离扫描图;
图3为发射信号频率顺序排列且初相为零时的信号处理局部距离扫描图;
图4为发射信号频率顺序排列且初相为零时的信号处理波束扫描图;
图5为发射信号频率顺序随机排列,初相为零时的信号处理全局距离扫描图;
图6为发射信号频率顺序随机排列,初相为零时的信号处理局部距离扫描图;
图7为发射信号频率顺序排列,初相随机变化时的信号处理全局距离扫描图;
图8为发射信号频率顺序排列,初相随机变化时的信号处理局部距离扫描图。
具体实施方式
为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及具体实施方式,对依据本发明提出的一种基于LFMCW的MIMO雷达信号处理方法进行详细说明。
有关本发明的前述及其他技术内容、特点及功效,在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明,可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解,然而所附附图仅是提供参考与说明之用,并非用来对本发明的技术方案加以限制。
应当说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的物品或者设备中还存在另外的相同要素。
请参见图1,图1为本发明实施例提供的一种基于LFMCW的MIMO雷达信号处理方法的流程示意图。所述方法包括如下步骤:
S1:建立一个包括M个发射阵元和N个接收阵元的MIMO雷达。
在本实施例中,设定所述MIMO雷达的发射天线与接收天线均为等距线阵,且发射天线包括M个阵元,接收天线包括N个阵元。雷达的M个发射阵元发射的信号表达式为:
s(t)=[s1(t),s2(t),...si(t),...sM(t)]T,i=1,2,...M
其中,si(t)为第i个发射阵元发射的信号,(·)T表示矩阵的转置。
S2:利用所述MIMO雷达发射一组正交频分LFMCW信号。
由于雷达发射阵列有M个阵元,每个阵元的发射信号具有不同的中心频率,占据不同的频带,因此可以得到M个信号波形,其中,第i个阵元产生的信号波形si(t)为:
Figure BDA0002978734320000061
其中,i为阵元序号,i=1,2,...,M,δi为第i个阵元发射信号的初相;
Figure BDA0002978734320000062
为调频斜率,B为单个信号的带宽,T为时宽;fi=f0+ciΔf,其中f0为中心载频,ci为频率编码,
Figure BDA0002978734320000063
为步进频率,Te为有效时宽,有效时宽是指发射信号与接收信号在时间轴上的重叠部分,其值可通过雷达参数设置和实验确定,在本实施例中,扔掉十分之一无效段,取Te=0.9T。
S3:根据所述正交频分LFMCW信号的回波信号,获得有效差拍信号。
所述正交频分LFMCW信号通过发射天线向外辐射到空中,遇障碍物后返回,形成回波信号被接收天线接收,对所述回波信号进行混频、低通滤波处理,得到差拍信号后进行采样、截取,可以得到有效差拍信号。
具体地,所述S3包括:
S31:将所述N个接收阵元接收的回波信号Y(t)与第一个发射阵元发射的信号s1(t)进行混频、低通滤波处理后,得到N个差拍信号sc(t):
sc(t)=[sc1(t),sc2(t),...scq(t),...scN(t)]T,q=1,2,...N,
Y(t)=[y1(t),y2(t),...yq(t),...yN(t)]T,q=1,2,...N,
其中,
Figure BDA0002978734320000071
(·)*表示矩阵的共轭,δ1为第一个阵元发射信号的初相,f1=f0+c1Δf,c1为频率编码。
S32:对差拍信号sc(t)的每一行,即每一个差拍信号scq(t)进行等间隔采样,采样时间T为信号时宽,设定采样频率fS,总采样点数L=fST,采样后的差拍信号sc(n)表示为:
sc(n)=[sc1(n),sc2(n),...scq(n),...scN(n)]T,q=1,2,...N,
其中,
Figure BDA0002978734320000072
n为第n个采样点。
S33:将采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
其中,sceq(n)=w(n)scq(n),
Figure BDA0002978734320000073
L′=<0.1L>,<·>为取整函数,n=L',L'+1,...,L-1,Le=L-L'表示有效时宽内的采样点数。
具体地,采样后的差拍信号sc(n)长为L点,n=0,1,...L-1;但是此处只有后半部分n=L′,...L-1有效,因此需要对采样后的差拍信号进行截取,矩形窗函数w(n)在n=L′,...L-1是1,其余部分是0,因此采样后的差拍信号乘以矩形窗函数,就获得了截取后的信号,本实施例中有效时宽Te=0.9T,L′=<0.1L>,因0.1L未必是整数,所以对其取整。截取之后的信号是有效差拍信号,这个信号是从L′开始的,n=L′,...L-1。
S4:对所述有效差拍信号进行接收波束形成,获得接收波束形成输出。
设形成波束的指向为θ0,权系数为接收导向矢量b(θ0),根据有效差拍信号sce(n)得到接收波束形成输出,表示为:
u(n)=bT0)sce(n)
其中,b(θ0)=[1,exp(j2πdr sinθ0/λ),...exp(j2πdr(N-1)sinθ0/λ)]T,λ为信号波长,dr为接收阵元间距。
S5:获得雷达距离扫描范围内的K个距离扫描点。
根据所述MIMO雷达设定的参数和采样频率,确定雷达距离扫描范围和距离间隔,得到K个距离扫描点,其中,所述K个距离扫描点与所述MIMO雷达之间的距离分别为R1,R2,...Rk,...RK,K≥Le,Le表示有效时宽内的采样点数。
S6:根据所述K个距离扫描点,计算得到一组扫描差拍信号,再对所述一组扫描差拍信号进行采样和截取,得到一组扫描有效差拍信号。
具体地,S6包括:
S61:获得第k个距离扫描点的时延:
Figure BDA0002978734320000081
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离。
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT0)s(t-τk)
s(t-τk)=[s1(t-τk),s2(t-τk),...si(t-τk),...sM(t-τk)]T,i=1,2,...M
其中,a(θ0)=[1,exp(j2πdt sinθ0/λ),...exp(j2πdt(M-1)sinθ0/λ)]T,dt为发射阵元间距,s(t-τk)表示第k个扫描点带时延的LFMCW信号,si(t-τk)表示第i个发射阵元第k个扫描点带时延的LFMCW信号:
Figure BDA0002978734320000082
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t),进行混频、低通滤波后得到第k个扫描点的差拍信号hck(t):
Figure BDA0002978734320000091
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
Figure BDA0002978734320000092
其中,采样时间为信号时宽T,采样频率fS,L=fST为总采样点数,n为第n个采样点。
S65:对采样后的第k个扫描差拍信号hck(n)乘以矩形窗函数,得到第k个扫描有效差拍信号hcek(n):
hcek(n)=w(n)hck(n)
其中,n=L′,L′+1,...,L-1,hcek(n)的长度为Le=L-L′。
S66:重复步骤S61至S65,得到K个距离扫描点的一组扫描有效差拍信号:
hce1(n),hce2(n),...hcek(n),...hceK(n),k=1,2,...,K。
S7:接收波束形成的输出与一组扫描有效差拍信号进行互相关,得到时空匹配滤波结果。
根据步骤S4得到的接收波束形成输出u(n)与步骤S6得到的所述扫描有效差拍信号hce1(n),hce2(n),...hcek(n),...hceK(n)进行互相关,得到时空匹配滤波后的信号p:
p=[p1,p2,...pk,...pK]T,k=1,2,...,K
其中,
Figure BDA0002978734320000093
进一步地,本发明实施例的基于LFMCW的MIMO雷达信号处理方法的效果可通过以下仿真进一步说明:
(1)仿真条件
假设所述MIMO雷达的发射天线的阵元个数M=4,接收天线的阵元个数N=4,发射阵元间距dt=λ/2,接收阵元间距dr=λ/2;发射信号单个带宽B=25MHz,时宽T=300μs,有效时宽取0.9倍的时宽,即Te=270μs;差拍信号采样频率fS=2.5MHz,采样点数L=750点,有效时宽内的采样点数Le=675点;距离扫描点K=2048点,对应最大扫描距离RK=4500m。假设目标距离雷达2000m,目标与雷达的法线夹角为20°,信噪比SNR=10dB。
(2)仿真内容与结果
仿真1:发射LFMCW信号频率顺序排列,即频率编码[c1,c2,c3,c4]=[0,1,2,3],且初相为零,利用本发明实施例的方法对目标回波信号进行信号处理,即波束指向θ0对准仿真设置的目标角度时,距离扫描点改变,从0m到4500m进行扫描,全局距离扫描图如图2所示,局部距离扫描图如图3所示,图中横坐标是距离,纵坐标是归一化幅度。由图2可以看出:目标对应的距离是2000m,与仿真设置的目标距离相同。由图3可以看出:距离向主瓣较宽,而距离分辨率ΔR=c/2B=6m,仿真结果不符合距离分辨率的理论计算结果。
仿真2:发射LFMCW信号频率顺序排列,即频率编码[c1,c2,c3,c4]=[0,1,2,3],且初相为零,利用本发明实施例的方法对目标回波信号进行空域搜索,即距离扫描点Rk对准仿真设置的目标距离时,波束指向θ0改变,从-90°到90°进行扫描,接收导向矢量和发射导向矢量同时改变,得到双程扫描结果;常规相控阵雷达只能改变接收导向矢量,得到单程扫描结果,如图4所示,图中横坐标是角度,纵坐标是归一化幅度。由图4可以看出:目标对应的角度是20°,与仿真设置的目标角度相同;基于LFMCW的MIMO雷达的波束扫描图波束主瓣宽度比常规相控阵雷达的窄,说明MIMO雷达的角度分辨率更好,测角精度更高。
仿真3:发射LFMCW信号频率顺序随机排列,即频率编码设为[c1,c2,c3,c4]=[0,3,1,2],初相为零,利用本发明实施例的方法对目标回波信号进行信号处理,全局距离扫描图如图5所示,局部距离扫描图如图6所示,图中横坐标是距离,纵坐标是归一化幅度。由图5可以看出:目标对应的距离是2000m,与仿真设置的目标距离相同。由图6可以看出:距离向主瓣变窄,与距离分辨率的理论计算结果一致,距离分辨率提高。
仿真4:发射LFMCW信号频率顺序排列,即频率编码[c1,c2,c3,c4]=[0,1,2,3],初相随机变化,利用本发明实施例的方法对目标回波信号进行信号处理,全局距离扫描图如图7所示,局部距离扫描图如图8所示,图中横坐标是距离,纵坐标是归一化幅度。由图7可以看出:目标对应的距离是2000m,与仿真设置的目标距离相同。由图8可以看出:距离向主瓣变窄,与距离分辨率的理论计算结果一致,距离分辨率提高。
无论是仿真3中的改变发射信号频率顺序,还是仿真4中的改变发射信号初相,都不影响仿真2的波束扫描结果。本发明实施例提出的将LFMCW和MIMO结合的新体制雷达及其信号处理方法可用,雷达天线阵列为均匀线阵时,发射信号频率顺序随机排列,或者发射信号的初相随机变化,可有效提高距离分辨率。
综上,本发明实施例基于LFMCW的MIMO雷达信号处理方法,利用频率编码,提出一种发射阵元同时发射一组正交LFMCW信号的MIMO雷达,并给出其信号处理实现方法,利用扫描差拍信号对接收信号进行距离扫描,得到时空匹配滤波结果,实现雷达目标检测。在抗杂波、抗干扰、低截获、峰值功率低、无距离盲区、距离及角度分辨率等方面具有明显优势。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (8)

1.一种基于LFMCW的MIMO雷达信号处理方法,其特征在于,包括:
S1:建立包括M个发射阵元和N个接收阵元的MIMO雷达;
S2:利用所述MIMO雷达发射一组正交频分LFMCW信号;
S3:根据所述正交频分LFMCW信号的回波信号,获得有效差拍信号;
S4:对所述有效差拍信号进行接收波束形成,获得接收波束形成输出;
S5:获得雷达距离扫描范围内的K个距离扫描点;
S6:根据所述K个距离扫描点获得一组扫描有效差拍信号;
S7:将所述接收波束形成输出与所述一组扫描有效差拍信号进行互相关,得到时空匹配滤波后的信号。
2.根据权利要求1所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,发射天线与接收天线均为等距线阵,所述M个发射阵元发射的信号为:
s(t)=[s1(t),s2(t),...si(t),...sM(t)]T,i=1,2,...M,
其中,si(t)为第i个发射阵元发射的信号,(·)T表示矩阵的转置。
3.根据权利要求2所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S2包括:
根据雷达发射阵列的M个发射阵元,得到每个发射阵元产生的正交频分LFMCW信号:
Figure FDA0003834957680000011
其中,si(t)为第i个阵元产生的信号波形,δi为第i个阵元发射信号的初相,
Figure FDA0003834957680000012
为调频斜率,B为单个信号的带宽,T为信号时宽,t表示0到T内的时间,fi=f0+ciΔf,f0为中心载频,ci为频率编码,
Figure FDA0003834957680000013
为步进频率,Te为有效时宽,j表示复数符号。
4.根据权利要求3所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S3包括:
S31:将所述N个接收阵元接收的回波信号Y(t)与第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到N个差拍信号sc(t):
sc(t)=[sc1(t),sc2(t),...scq(t),...scN(t)]T,q=1,2,...N,
Y(t)=[y1(t),y2(t),...yq(t),...yN(t)]T,q=1,2,...N,
其中,
Figure FDA0003834957680000014
(·)*表示矩阵的共轭,δ1为第一个阵元发射信号的初相,f1=f0+c1Δf,c1为频率编码;
S32:对所述差拍信号sc(t)中的每一个差拍信号scq(t)进行等间隔采样,获得采样后的差拍信号sc(n):
sc(n)=[sc1(n),sc2(n),...scq(n),...scN(n)]T,q=1,2,...N,
其中,
Figure FDA0003834957680000021
L=fST为总采样点数,n为第n个采样点,采样时间为信号时宽T,fS为采样频率;
S33:将所述采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
其中,sceq(n)=w(n)scq(n),
Figure FDA0003834957680000022
L′=<0.1L>,<·>为取整函数。
5.根据权利要求4所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S4包括:
根据有效差拍信号sce(n),得到接收波束形成输出:
u(n)=bT0)sce(n),
其中,θ0为形成的接收波束的指向,b(θ0)为接收导向矢量,b(θ0)=[1,exp(j2πdrsinθ0/λ),...exp(j2πdr(N-1)sinθ0/λ)]T,λ为信号波长,dr为接收阵元的间距。
6.根据权利要求5所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S5包括:
根据所述MIMO雷达设定的参数和采样频率,确定雷达距离扫描范围和距离间隔,得到K个距离扫描点,其中,所述K个距离扫描点与所述MIMO雷达之间的距离分别为R1,R2,...Rk,...RK,K≥Le,Le表示有效时宽内的采样点数。
7.根据权利要求6所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S6包括:
S61:获得第k个距离扫描点的时延:
Figure FDA0003834957680000023
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离;
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT0)s(t-τk)
s(t-τk)=[s1(t-τk),s2(t-τk),...si(t-τk),...sM(t-τk)]T,i=1,2,...M
其中,a(θ0)=[1,exp(j2πdtsinθ0/λ),...exp(j2πdt(M-1)sinθ0/λ)]T,dt为发射阵元间距,s(t-τk)表示第k个扫描点带时延的LFMCW信号,si(t-τk)表示第i个发射阵元第k个扫描点带时延的LFMCW信号:
Figure FDA0003834957680000031
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到第k个扫描点的差拍信号hck(t):
Figure FDA0003834957680000032
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
Figure FDA0003834957680000033
其中,采样时间为信号时宽T,采样频率fS,L=fST为总采样点数,n为第n个采样点;
S65:对采样后第k个扫描差拍信号hck(n)乘以矩形窗函数w(n),得到第k个扫描有效差拍信号hcek(n):
hcek(n)=w(n)hck(n),
其中,n=L′,L′+1,...,L-1;
S66:重复步骤S61至S65,得到K个距离扫描点的一组扫描有效差拍信号:
hce1(n),hce2(n),...hcek(n),...hceK(n),k=1,2,...,K。
8.根据权利要求7所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S7包括:
根据所述接收波束形成输出u(n)与所述一组扫描有效差拍信号hce1(n),hce2(n),...hcek(n),...hceK(n)进行互相关,得到时空匹配滤波后的信号p:
p=[p1,p2,...pk,...pK]T,k=1,2,...,K
其中,
Figure FDA0003834957680000041
CN202110280931.1A 2021-03-16 2021-03-16 一种基于lfmcw的mimo雷达信号处理方法 Active CN113204017B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110280931.1A CN113204017B (zh) 2021-03-16 2021-03-16 一种基于lfmcw的mimo雷达信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110280931.1A CN113204017B (zh) 2021-03-16 2021-03-16 一种基于lfmcw的mimo雷达信号处理方法

Publications (2)

Publication Number Publication Date
CN113204017A CN113204017A (zh) 2021-08-03
CN113204017B true CN113204017B (zh) 2022-11-04

Family

ID=77025452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110280931.1A Active CN113204017B (zh) 2021-03-16 2021-03-16 一种基于lfmcw的mimo雷达信号处理方法

Country Status (1)

Country Link
CN (1) CN113204017B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102798856B (zh) * 2012-09-11 2014-08-20 重庆大学 一种基于小波包的24GHz LFMCW雷达测距方法
CN104931948B (zh) * 2015-05-25 2018-03-06 西安电子科技大学昆山创新研究院 基于普通波束扫描的fda雷达的一种接收方案的改进方法
JP2018072014A (ja) * 2016-10-25 2018-05-10 ソニーセミコンダクタソリューションズ株式会社 レーダ装置、信号処理装置及び信号処理方法
CN108226871B (zh) * 2018-01-18 2021-08-10 西安电子工程研究所 基于毫米波避撞雷达的分集相控阵设计方法
CN110456334B (zh) * 2019-07-27 2023-04-28 南京理工大学 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法
CN110412514B (zh) * 2019-08-02 2023-05-23 南京慧尔视智能科技有限公司 一种mimo体制下的线性调频连续波波形优化方法

Also Published As

Publication number Publication date
CN113204017A (zh) 2021-08-03

Similar Documents

Publication Publication Date Title
CN110456334B (zh) 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法
Wang et al. Nonuniform frequency diverse array for range-angle imaging of targets
CN111090080B (zh) 基于空时编码阵列的超宽带雷达单通道数字波束形成方法
US20070152871A1 (en) Time duplex apparatus and method for radar sensor front-ends
Wang et al. FDA radar ambiguity function characteristics analysis and optimization
CN109375213B (zh) 基于子阵划分的频率分集阵列信号处理方法
CN110412570B (zh) 基于空间脉冲相位编码的hrws-sar成像方法
CN110346764B (zh) 一种mimo雷达目标距离角度解耦方法
CN111257879B (zh) 基于二范数解决毫米波mimo雷达目标分裂的方法
CN114152916A (zh) 一种基于先导信息的mimo雷达波形正交方法
Wang et al. Range-azimuth decouple beamforming for frequency diverse array with Costas-sequence modulated frequency offsets
CN111352079B (zh) 一种基于频控阵mimo雷达的低截获系统的设计方法
Xu et al. Range-angle matched receiver for coherent FDA radars
CN113238212B (zh) 基于空时编码的频率分集阵列雷达距离分辨率增强方法
CN114895261A (zh) 基于多频子脉冲编码阵列的杂波抑制方法
Chen et al. The applications and future of synthetic impulse and aperture radar
CN113204018B (zh) 一种基于lfmcw的mimo雷达快速信号处理方法
JP2006145251A (ja) 電波到来方向推定装置
CN113204017B (zh) 一种基于lfmcw的mimo雷达信号处理方法
CN112612013A (zh) Fda-mimo雷达增量距离-角度二维波束形成方法
Calvary et al. Spatio-temporal coding for radar array processing
Takayama et al. Hybrid SIMO and MIMO sparse array radar
Landi et al. Time-orthogonal-waveform-space-time adaptive processing for distributed aperture radars
CN114764136A (zh) 基于多时间尺度耦合网络的雷达抗干扰波形产生方法
CN113376590B (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