CN113204017A - 一种基于lfmcw的mimo雷达信号处理方法 - Google Patents
一种基于lfmcw的mimo雷达信号处理方法 Download PDFInfo
- Publication number
- CN113204017A CN113204017A CN202110280931.1A CN202110280931A CN113204017A CN 113204017 A CN113204017 A CN 113204017A CN 202110280931 A CN202110280931 A CN 202110280931A CN 113204017 A CN113204017 A CN 113204017A
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 21
- 238000001914 filtration Methods 0.000 claims abstract description 16
- 238000005070 sampling Methods 0.000 claims description 43
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 2
- 238000010009 beating Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 16
- 238000000034 method Methods 0.000 abstract description 8
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 13
- 230000003111 delayed effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000017105 transposition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
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
本发明公开了一种基于LFMCW的MIMO雷达信号处理方法,包括:建立包括M个发射阵元和N个接收阵元的MIMO雷达;利用MIMO雷达发射一组正交频分LFMCW信号;根据正交频分LFMCW信号的回波信号,获得有效差拍信号;对有效差拍信号进行接收波束形成,获得接收波束形成输出;获得雷达距离扫描范围内的K个距离扫描点;根据K个距离扫描点获得一组扫描有效差拍信号;将接收波束形成输出与一组扫描有效差拍信号进行互相关,得到时空匹配滤波后的信号。本发明提出一种基于LFMCW的MIMO雷达,并给出其信号处理实现方法,利用扫描差拍信号对接收信号进行距离扫描,得到时空匹配滤波结果,实现雷达目标检测。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种基于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信号:
其中,si(t)为第i个阵元产生的信号波形,δi为第i个阵元发射信号的初相,为调频斜率,B为单个信号的带宽,T为信号时宽,t表示0到T内的时间,fi=f0+ciΔf,f0为中心载频,ci为频率编码,为步进频率,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,
S32:对所述差拍信号sc(t)中的每一个差拍信号scq(t)进行等间隔采样,获得采样后的差拍信号sc(n):
sc(n)=[sc1(n),sc2(n),...scq(n),...scN(n)]T,q=1,2,...N,
S33:将所述采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
在本发明的一个实施例中,所述S4包括:
根据有效差拍信号sce(n),得到接收波束形成输出:
u(n)=bT(θ0)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个距离扫描点的时延:
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离;
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT(θ0)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信号:
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到第k个扫描点的差拍信号hck(t):
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
其中,采样时间为信号时宽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
与现有技术相比,本发明的有益效果在于:
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)为:
其中,i为阵元序号,i=1,2,...,M,δi为第i个阵元发射信号的初相;为调频斜率,B为单个信号的带宽,T为时宽;fi=f0+ciΔf,其中f0为中心载频,ci为频率编码,为步进频率,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,
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,
S33:将采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
具体地,采样后的差拍信号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)=bT(θ0)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个距离扫描点的时延:
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离。
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT(θ0)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信号:
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t),进行混频、低通滤波后得到第k个扫描点的差拍信号hck(t):
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
其中,采样时间为信号时宽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
进一步地,本发明实施例的基于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表示矩阵的转置。
4.根据权利要求1所述的基于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,
S32:对所述差拍信号sc(t)中的每一个差拍信号scq(t)进行等间隔采样,获得采样后的差拍信号sc(n):
sc(n)=[sc1(n),sc2(n),...scq(n),...scN(n)]T,q=1,2,...N,
S33:将所述采样后的差拍信号sc(n)乘以矩形窗函数w(n),得到有效差拍信号sce(n):
sce(n)=[sce1(n),sce2(n),...sceq(n),...sceN(n)]T,q=1,2,...N,
5.根据权利要求4所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S4包括:
根据有效差拍信号sce(n),得到接收波束形成输出:
u(n)=bT(θ0)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。
7.根据权利要求6所述的基于LFMCW的MIMO雷达信号处理方法,其特征在于,所述S6包括:
S61:获得第k个距离扫描点的时延:
其中,c为光速,τk为第k个距离扫描点的时延,Rk为第k个距离扫描点与所述MIMO雷达之间的距离;
S62:根据发射导向矢量a(θ0)和第k个扫描点的时延τk,获得第k个扫描点的接收信号zk(t):
zk(t)=aT(θ0)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信号:
S63:根据第k个扫描点的接收信号zk(t)和第一个发射阵元发射的信号s1(t)进行混频和低通滤波处理后,得到第k个扫描点的差拍信号hck(t):
S64:对第k个扫描点的差拍信号hck(t)进行等间隔采样,获得采样后第k个扫描差拍信号hck(n):
其中,采样时间为信号时宽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。
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 true CN113204017A (zh) | 2021-08-03 |
CN113204017B 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) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798856A (zh) * | 2012-09-11 | 2012-11-28 | 重庆大学 | 一种基于小波包的24GHz LFMCW雷达测距方法 |
CN104931948A (zh) * | 2015-05-25 | 2015-09-23 | 西安电子科技大学 | 基于普通波束扫描的fda雷达第一种接收方案改进方法 |
JP2018072014A (ja) * | 2016-10-25 | 2018-05-10 | ソニーセミコンダクタソリューションズ株式会社 | レーダ装置、信号処理装置及び信号処理方法 |
CN108226871A (zh) * | 2018-01-18 | 2018-06-29 | 西安电子工程研究所 | 基于毫米波避撞雷达的分集相控阵设计方法 |
CN110456334A (zh) * | 2019-07-27 | 2019-11-15 | 南京理工大学 | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 |
WO2021022683A1 (zh) * | 2019-08-02 | 2021-02-11 | 南京慧尔视智能科技有限公司 | 一种mimo体制下的线性调频连续波波形优化方法 |
-
2021
- 2021-03-16 CN CN202110280931.1A patent/CN113204017B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798856A (zh) * | 2012-09-11 | 2012-11-28 | 重庆大学 | 一种基于小波包的24GHz LFMCW雷达测距方法 |
CN104931948A (zh) * | 2015-05-25 | 2015-09-23 | 西安电子科技大学 | 基于普通波束扫描的fda雷达第一种接收方案改进方法 |
JP2018072014A (ja) * | 2016-10-25 | 2018-05-10 | ソニーセミコンダクタソリューションズ株式会社 | レーダ装置、信号処理装置及び信号処理方法 |
CN108226871A (zh) * | 2018-01-18 | 2018-06-29 | 西安电子工程研究所 | 基于毫米波避撞雷达的分集相控阵设计方法 |
CN110456334A (zh) * | 2019-07-27 | 2019-11-15 | 南京理工大学 | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 |
WO2021022683A1 (zh) * | 2019-08-02 | 2021-02-11 | 南京慧尔视智能科技有限公司 | 一种mimo体制下的线性调频连续波波形优化方法 |
Non-Patent Citations (2)
Title |
---|
XIANG ZHAO等: "Performance analysis of airborne LFMCW-MIMO vitual array radar", 《2016 CIE INTERNATIONAL CONFERENCE ON RADAR (RADAR)》 * |
邹林等: "LFMCW车载雷达解速度模糊测角新方法", 《电子科技大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113204017B (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | MIMO-OFDM radar for direction estimation | |
CN109375213B (zh) | 基于子阵划分的频率分集阵列信号处理方法 | |
CN111090080B (zh) | 基于空时编码阵列的超宽带雷达单通道数字波束形成方法 | |
CN110456334A (zh) | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 | |
CN110109065B (zh) | 基于lfm信号的mimo雷达部分相关波形设计方法 | |
CN106569212A (zh) | 一种多通道sar‑gmti距离模糊杂波抑制方法 | |
CN104155647B (zh) | 基于机扫米波阵列雷达的估计目标方位角的方法 | |
CN110456314B (zh) | 基于主瓣展宽的集中式mimo雷达波形优化方法 | |
CN111257879B (zh) | 基于二范数解决毫米波mimo雷达目标分裂的方法 | |
CN112612013A (zh) | Fda-mimo雷达增量距离-角度二维波束形成方法 | |
Xu et al. | Range-angle matched receiver for coherent FDA radars | |
CN114895261A (zh) | 基于多频子脉冲编码阵列的杂波抑制方法 | |
Chen et al. | Transmit beamforming based on 4-D antenna arrays for low probability of intercept systems | |
Deng et al. | A virtual antenna beamforming (VAB) approach for radar systems by using orthogonal coding waveforms | |
CN108693526B (zh) | 用于mimo雷达的运动目标相位补偿方法 | |
CN113376590A (zh) | 一种空时编码阵列雷达中阵元间时延量的设计方法 | |
CN113238212B (zh) | 基于空时编码的频率分集阵列雷达距离分辨率增强方法 | |
CN106529440A (zh) | 相干频率分集阵雷达分段匹配滤波方法 | |
CN111812612A (zh) | 基于子阵正交lfm信号的mimo雷达部分相关波形设计方法 | |
CN109946655B (zh) | Mimo雷达lfm信号置零波形设计方法 | |
CN113204017B (zh) | 一种基于lfmcw的mimo雷达信号处理方法 | |
CN113204018B (zh) | 一种基于lfmcw的mimo雷达快速信号处理方法 | |
Gao et al. | Frequency diverse array antennas with random logarithmically increasing frequency offset | |
CN117148303A (zh) | 一种基于图信号处理的fda-mimo雷达二维参数估计方法 | |
CN117075030A (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 |