CN112710993B - 一种tdm-mimo雷达空间谱估计补偿方法 - Google Patents

一种tdm-mimo雷达空间谱估计补偿方法 Download PDF

Info

Publication number
CN112710993B
CN112710993B CN201911022681.0A CN201911022681A CN112710993B CN 112710993 B CN112710993 B CN 112710993B CN 201911022681 A CN201911022681 A CN 201911022681A CN 112710993 B CN112710993 B CN 112710993B
Authority
CN
China
Prior art keywords
radar
echo
signal
target
array element
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
CN201911022681.0A
Other languages
English (en)
Other versions
CN112710993A (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.)
Arkmicro Technologies Inc
Original Assignee
Arkmicro Technologies Inc
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 Arkmicro Technologies Inc filed Critical Arkmicro Technologies Inc
Priority to CN201911022681.0A priority Critical patent/CN112710993B/zh
Publication of CN112710993A publication Critical patent/CN112710993A/zh
Application granted granted Critical
Publication of CN112710993B publication Critical patent/CN112710993B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system

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

本发明提供一种TDM‑MIMO雷达空间谱估计补偿方法,用于提高雷达对目标的测角精度,该方法包括:确定目标与雷达之间的距离以及目标相对雷达的速度;将雷达锯齿波的回波信号与发射信号进行混频和滤波后得到差频信号;根据回波延迟时间差确定虚拟阵元需要补偿的相位,根据差频信号和所述相位确定补偿后的虚拟阵元回波差频信号。该方法可以有效提升TDM‑MIMO雷达空间谱估计的估计精度,对于小阵列毫米波雷达测角有重要的理论意义和工程价值。

Description

一种TDM-MIMO雷达空间谱估计补偿方法
技术领域
本发明属于雷达技术领域,尤其涉及一种TDM-MIMO雷达空间谱估计补偿方法。
背景技术
目前的时分复用多输入多输出(英文全称Time Division Multiplexing-Multiple Input Multiple Output,简称TDM-MIMO)雷达不但避开了发射正交波形的工程难题,而且扩展了孔径。时分复用多输入多输出雷达广泛应用于车载毫米波雷达。
把无线通信系统中的多个输入和多个输出技术引入到雷达领域,并和数字阵列技术相结合而产生的一种新体制雷达。MIMO雷达的概念在2003年由美国林肯实验室的Bliss和Forsythe首次提出,目前已成为国内外雷达界的一个研究热点。在MIMO雷达定义中,多输入是指同时发射多种雷达信号波形(一般是多个天线同时发射不同的波形),多输出是指多个天线同时接收并通过多路接收机输出以获得多通道空间采样信号。因此,MIMO雷达可以看成是相控阵雷达的扩展,具有可控自由度多、发射波束设置灵活、测角精度高以及抗干扰能力强等优点,在军事和民用领域都有着广阔的应用前景,是雷达下一步发展的重要方向。
MIMO雷达技术应用于车载毫米波雷达的好处主要是利用其虚拟孔径扩展能力提升测角精度,不同的天线收发位置得到的虚拟天线孔径不同,所带来的天线孔径扩展也不同。相对于传统的相控阵雷达体制能扩展孔径40%~100%,此时,其测角精度也将提升40%~100%。
但是对于毫米波雷达(比如车载毫米波雷达,频率在77~81GHz)的应用来说,由于频率高,同时发射正交波形在工程上面临困难,于是一些研究者提出了TDM-MIMO雷达结构,这样不但避开了发射正交波形的工程难题,而且扩展了孔径。由于引入了时分复用,导致了形成虚拟孔径的各阵元数据不是来自于同一时刻目标回波,不同脉冲周期的数据失去了相干性,此时,如果直接用不同时刻的回波数据进行角度估计将会带来很大的误差(误差取决于雷达速度和目标速度)。仿真发现,对于2发4收的TDM-MIMO雷达系统,当雷达速度20m/s,目标相对于雷达的速度10m/s时,空间谱估计测角误差较大。
发明内容
有鉴于此,本发明提供一种TDM-MIMO雷达空间谱估计补偿方法用于解决现有技术的不足。
为达到上述目的,本发明采用如下技术方案:
本发明实施例提供了一种TDM-MIMO雷达空间谱估计补偿方法,用于提高雷达对目标的测角精度,该方法包括:
步骤一、确定目标与雷达之间的距离以及目标相对雷达的速度;
步骤二、将雷达锯齿波的回波信号与发射信号进行混频和滤波后得到差频信号;
(1)雷达锯齿波发射信号频率f(t)随时间变化的数学表达式为:
f(t)=fc+αt,0≤t≤Tc,其中:fc表示发射信号的中心频率,α表示调制信号的斜率,Tc表示单个锯齿波周期;
因此,单个锯齿波的发射信号St(t)的数学表达式为:
其中:At表示发射信号的幅度,nt(t)表示发射信号中的噪声;
(2)单个锯齿波经目标反射的回波信号Sr(t)的数学表达式为:
其中:Ar表示回波信号的幅度,nr(t)表示回波信号中的噪声,τ表示回波延迟时间;回波延迟时间τ的数学表达式为:
其中:R表示目标与雷达之间的距离,v表示目标相对雷达的速度,c表示光速;
(3)将回波信号与发射信号进行混频和滤波得到差频信号S0(t):
其中:/>
步骤三、根据回波延迟时间差确定虚拟阵元需要补偿的相位,根据差频信号和所述相位确定补偿后的虚拟阵元回波差频信号:
(1)第n+1个脉冲形成的虚拟阵元相对第n个脉冲的接收阵元有一个回波延迟时间差,该回波延迟时间差Δτ的数学表达式为:
(2)由回波延迟时间差Δτ得到虚拟阵元需要补偿的相位其数学表达式为:
(3)将虚拟阵元的回波差频信号乘以复向量得到补偿后的虚拟阵元回波差频信号:
其中:Sc(t)表示补偿后的虚拟阵元回波差频信号,S0(t)表示补偿前的虚拟阵元回波差频信号。
进一步地,所述确定目标与雷达之间的距离以及目标相对雷达的速度具体包括:
由单个阵元的采样数据进行FFT处理得到目标的差频信号,由差频信号得到目标的距离信息R;
在目标距离单元对不同脉冲数据再进行FFT处理,得到目标相对于雷达的速度v。
进一步地,所述雷达阵元数目为2发4收,雷达频率为77~77.4GHz。
本申请提供的技术方案,可以有效提升TDM-MIMO雷达空间谱估计的估计精度,对于小阵列毫米波雷达测角有重要的理论意义和工程价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种TDM-MIMO雷达空间谱估计补偿方法的流程示意图;
图2为不补偿虚拟阵元的8阵元空间谱估计结果;
图3为补偿虚拟阵元的8阵元空间谱估计结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供了一种TDM-MIMO雷达空间谱估计补偿方法,用于提高雷达对目标的空间谱估计测角精度,该方法包括:
S101、确定目标与雷达之间的距离以及目标相对雷达的速度。
确定目标与雷达之间的距离以及目标相对雷达的速度具体如下:由单个阵元的采样数据进行FFT处理得到目标的差频信号,由差频信号得到目标的距离信息R;在目标距离单元对不同脉冲数据再进行FFT处理,得到目标相对于雷达的速度v。
S102、将雷达锯齿波的回波信号与发射信号进行混频和滤波后得到差频信号。
(1)雷达锯齿波发射信号频率f(t)随时间变化的数学表达式为:
f(t)=fc+αt,0≤t≤Tc (1)
其中:fc表示发射信号的中心频率,α表示调制信号的斜率,Tc表示单个锯齿波周期。
因此,单个锯齿波的发射信号St(t)的数学表达式为:
其中:At表示发射信号的幅度,nt(t)表示发射信号中的噪声。
(2)单个锯齿波经目标反射的回波信号Sr(t)的数学表达式为:
其中:Ar表示回波信号的幅度,nr(t)表示回波信号中的噪声,τ表示回波延迟时间;回波延迟时间τ的数学表达式为:
其中:R表示目标与雷达之间的距离,v表示目标相对雷达的速度,c表示光速。
(3)将回波信号与发射信号进行混频和滤波得到差频信号S0(t):
其中:
S103、根据回波延迟时间差确定虚拟阵元需要补偿的相位,根据差频信号和所述相位确定补偿后的虚拟阵元回波差频信号。
(1)第n+1个脉冲形成的虚拟阵元相对第n个脉冲的接收阵元有一个回波延迟时间差,该回波延迟时间差Δτ的数学表达式为:
对于2发4收的MIMO雷达系统,将第一个脉冲的回波和第二个脉冲的回波一起处理可以得到8个阵元的虚拟阵列,但是此时第二个脉冲形成的虚拟阵元相对第一个脉冲的接收阵元有一个回波延迟时间差,该回波延迟时间差如果发射单元数大于2时,只需将式(6)乘以相应的系数便可以得到相应的回波延迟时间差,比如,第三个脉冲形成的虚拟阵元相对第一个脉冲的接收阵元的回波延迟时间差是2Δτ,第四个脉冲形成的虚拟阵元相对第一个脉冲的接收阵元的回波延迟时间差是3Δτ,以此类推。
(2)由回波延迟时间差Δτ得到虚拟阵元需要补偿的相位其数学表达式为:
(3)将虚拟阵元的回波差频信号乘以复向量得到补偿后的虚拟阵元回波差频信号:
其中:Sc(t)表示补偿后的虚拟阵元回波差频信号,S0(t)表示补偿前的虚拟阵元回波差频信号。
不是所有数据都需要补偿。只需要补偿虚拟阵元数据,对于2发4收,就是补偿偶数脉冲的接收数据。
下面以一个具体仿真实例加以说明:
设定雷达基本参数。雷达体制选用TDM-MIMO雷达;阵元数目选用2发4收;频率选用77~77.4GHz;发射波形选用锯齿波FMCW(线性调频连续波);单个锯齿波周期Tc为20μs;速度计算周期为2.56ms;采样率为40MHz;调频斜率α为20THz/s;雷达对地速度30m/s。
设定目标基本参数。目标与雷达之间的距离R为50m;速度为20m/s(速度为正表示靠近雷达);角度为0°(雷达天线法向);信噪比为15dB。
图2和图3分别给出了2发4收TDM-MIMO模式下不补偿虚拟阵元和补偿虚拟阵元的8阵元空间谱估计仿真结果。其中左图为角度为1度的搜索结果,右图为角度为0.1度的搜索结果。从图2中可以看出,如果不补偿则不但角度估计偏差很大,而且会产生虚假的目标。而经本发明方法补偿后目标角度估计精度接近理想情况。
本发明针对TDM-MIMO雷达的这个问题,提出了补偿方法,该方法首先由单个阵元数据进行FFT处理得到目标距离和目标速度,然后由目标距离和目标速度计算得到补偿因子,对不同脉冲的数据进行补偿。计算公式仿真结果表明,对于2发4收的TDM-MIMO雷达系统,如果提供的速度值很精确,理想情况下可以接近完全补偿,当雷达速度和目标速度估计的速度误差不大于0.5m/s时,补偿不完全带来的测角误差也不会超过0.2°。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种TDM-MIMO雷达空间谱估计补偿方法,用于提高雷达对目标的测角精度,其特征在于,该方法包括:
步骤一、确定目标与雷达之间的距离以及目标相对雷达的速度;
步骤二、将雷达锯齿波的回波信号与发射信号进行混频和滤波后得到差频信号;
(1)雷达锯齿波发射信号频率f(t)随时间变化的数学表达式为:
f(t)=fc+αt,0≤t≤Tc,其中:fc表示发射信号的中心频率,α表示调制信号的斜率,Tc表示单个锯齿波周期;
因此,单个锯齿波的发射信号St(t)的数学表达式为:
其中:At表示发射信号的幅度,nt(t)表示发射信号中的噪声;
(2)单个锯齿波经目标反射的回波信号Sr(t)的数学表达式为:
其中:Ar表示回波信号的幅度,nr(t)表示回波信号中的噪声,τ表示回波延迟时间;回波延迟时间τ的数学表达式为:
其中:R表示目标与雷达之间的距离,v表示目标相对雷达的速度,c表示光速;
(3)将回波信号与发射信号进行混频和滤波得到差频信号S0(t):
其中:/>
步骤三、根据回波延迟时间差确定虚拟阵元需要补偿的相位,根据差频信号和所述相位确定补偿后的虚拟阵元回波差频信号:
(1)第n+1个脉冲形成的虚拟阵元相对第n个脉冲的接收阵元有一个回波延迟时间差,该回波延迟时间差Δτ的数学表达式为:
(2)由回波延迟时间差Δτ得到虚拟阵元需要补偿的相位其数学表达式为:
(3)将虚拟阵元的回波差频信号乘以复向量得到补偿后的虚拟阵元回波差频信号:
其中:Sc(t)表示补偿后的虚拟阵元回波差频信号,S0(t)表示补偿前的虚拟阵元回波差频信号。
2.根据权利要求1所述的方法,其特征在于,所述确定目标与雷达之间的距离以及目标相对雷达的速度具体包括:
由单个阵元的采样数据进行FFT处理得到目标的差频信号,由差频信号得到目标的距离信息R;
在目标距离单元对不同脉冲数据再进行FFT处理,得到目标相对于雷达的速度v。
3.根据权利要求1所述的方法,其特征在于,所述雷达阵元数目为2发4收,雷达频率为77~77.4GHz。
CN201911022681.0A 2019-10-25 2019-10-25 一种tdm-mimo雷达空间谱估计补偿方法 Active CN112710993B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911022681.0A CN112710993B (zh) 2019-10-25 2019-10-25 一种tdm-mimo雷达空间谱估计补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911022681.0A CN112710993B (zh) 2019-10-25 2019-10-25 一种tdm-mimo雷达空间谱估计补偿方法

Publications (2)

Publication Number Publication Date
CN112710993A CN112710993A (zh) 2021-04-27
CN112710993B true CN112710993B (zh) 2024-05-14

Family

ID=75540656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911022681.0A Active CN112710993B (zh) 2019-10-25 2019-10-25 一种tdm-mimo雷达空间谱估计补偿方法

Country Status (1)

Country Link
CN (1) CN112710993B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015109870A1 (zh) * 2014-01-24 2015-07-30 深圳大学 Mimo雷达系统及其目标端相位同步方法
CN107688178A (zh) * 2017-08-25 2018-02-13 上海通趣科技有限公司 一种基于77GHz毫米波雷达的锯齿波测距测速方法
CN108693511A (zh) * 2018-05-25 2018-10-23 中国人民解放军国防科技大学 时分复用mimo雷达的运动目标角度计算方法
CN109975806A (zh) * 2019-04-15 2019-07-05 西安中电科西电科大雷达技术协同创新研究院有限公司 基于阵元相位差的时分mimo雷达运动补偿方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7158076B2 (en) * 2004-10-01 2007-01-02 Bae Systems Information And Electronic Systems Integration Inc. Method and apparatus for correcting velocity-induced range estimate phase errors in a two-tone monopulse CW radar

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015109870A1 (zh) * 2014-01-24 2015-07-30 深圳大学 Mimo雷达系统及其目标端相位同步方法
CN107688178A (zh) * 2017-08-25 2018-02-13 上海通趣科技有限公司 一种基于77GHz毫米波雷达的锯齿波测距测速方法
CN108693511A (zh) * 2018-05-25 2018-10-23 中国人民解放军国防科技大学 时分复用mimo雷达的运动目标角度计算方法
CN109975806A (zh) * 2019-04-15 2019-07-05 西安中电科西电科大雷达技术协同创新研究院有限公司 基于阵元相位差的时分mimo雷达运动补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
77GHz多载频MIMO汽车雷达信号处理方法的研究;季晓宇 等;微波学报;20171229(第S1期);全文 *
Motion compensation for TDM MIMO radar by sparse reconstruction;Xueyao Hu 等;ELECTRONICS LETTERS;20171123;第53卷(第24期);全文 *

Also Published As

Publication number Publication date
CN112710993A (zh) 2021-04-27

Similar Documents

Publication Publication Date Title
CN111308437B (zh) 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN111537966B (zh) 一种适用于毫米波车载雷达领域的阵列天线误差校正方法
CN110412559B (zh) 分布式无人机mimo雷达的非相参融合目标检测方法
CN110412558A (zh) 基于tdm mimo的解车载fmcw雷达速度模糊方法
CN108693511B (zh) 时分复用mimo雷达的运动目标角度计算方法
CN109471064B (zh) 基于脉冲压缩技术的时间调制阵列测向系统
US20240219516A1 (en) Target detection method and apparatus, radar, and vehicle
CN102156279A (zh) 基于mimo的双基地雷达地面动目标检测方法
CN103454632B (zh) 一站固定式调频连续波双基地sar成像方法
CN103823217A (zh) 基于双频发射的双基地mimo雷达高速运动目标参数估计方法
CN103018727A (zh) 一种基于样本训练的机载雷达非平稳杂波抑制方法
CN111239721B (zh) 车载mimo雷达求熵解速度模糊的方法
CN110412570B (zh) 基于空间脉冲相位编码的hrws-sar成像方法
CN107356908B (zh) 一种捷变频信号相参积累方法
Zhu et al. A type of M 2-transmitter N 2-receiver MIMO radar array and 3D imaging theory
CN111257879A (zh) 基于二范数解决毫米波mimo雷达目标分裂的方法
CN113093141B (zh) 多载频lfmcw雷达信号合成处理方法
CN102540153A (zh) 一种基于站间直达波干扰的阵列幅相误差校准方法
CN114076945A (zh) 一种nlfm信号的电离层色散效应解析模型的建立方法
CN111736129B (zh) 空时调制双基地雷达系统及其实现方法和介质
CN112710993B (zh) 一种tdm-mimo雷达空间谱估计补偿方法
CN109557536B (zh) 一种测角方法、测角装置及测角系统
Anajemba et al. Efficient switched digital beamforming radar system based on SIMO/MIMO receiver
CN111934737B (zh) 一种宽带数字阵列的时延补偿方法
CN113126076B (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