Geomatics Science and Technology
Vol. 09  No. 04 ( 2021 ), Article ID: 45994 , 9 pages
10.12677/GST.2021.94018

分布式MIMO雷达参数估计与目标检测 性能研究

马鹏1,王石泉2,赵飞1

1陆军装甲兵学院作战保障系战斗支援教研室,安徽 蚌埠

2陆军装甲兵学院教学科研处,安徽 蚌埠

收稿日期:2021年3月8日;录用日期:2021年10月19日;发布日期:2021年10月26日

摘要

本文针对不同的目标模型对分布式MIMO雷达工作模式进行了相干处理和非相干处理的分类,推导出参数估计器与目标检测概率的最终表达式,论证了不同的布站方式以及收发站数对参数估计和目标检测性能的影响,支撑了分布式MIMO雷达相关理论的研究。仿真表明布站方式会较大影响参数估计与复杂目标检测的精度,并与传统相控阵雷达进行对比,证明了文章结论的有效性。

关键词

分布式,MIMO雷达,相干检测,非相干检测,参数估计

Parameter Estimation and Target Detection for Distributed MIMO Radar

Peng Ma1, Shiquan Wang2, Fei Zhao1

1Departmet of Combat Support, Research Section of Battle Assist, Army Armored Force Institute, Bengbu Anhui

2Section of Teaching and Research, Army Armored Force Institute, Bengbu Anhui

Received: Mar. 8th, 2021; accepted: Oct. 19th, 2021; published: Oct. 26th, 2021

ABSTRACT

Corresponding MIMO radar signal models are established according to the modality of signal processing in parameter estimation and detection, based on which the working patterns of MIMO radar can be classified as coherent processing and non-coherent processing. On this foundation, the performances of estimation and detection are analyzed in detail, primarily focusing on the difference of two models, which support the conclusion of the rule for configuring transmitters and receiver is provided based on the analysis of effect to performance. The simulation results show the effectiveness of the proposed algorithms by quantitatively analyzing the main factors affecting the capability of forensics.

Keywords:Distributed Radar, MIMO Radar, Coherent Detecting, Non-Coherent Detecting, Parameter Estimation

Copyright © 2021 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

作为新体制雷达的代表,MIMO (MIMO: Multiple Input Multiple Output)雷达具有诸多优点,一直是学术界研究的热点。分布式MIMO雷达采用大间距配置,使雷达对目标的观测经历角度扩展,可充分利用空间分集增益以提高雷达系统的检测和估计性能 [1],代表了未来分布式雷达系统的发展方向,具有更大的研究潜力。但是相对于集中式MIMO雷达而言,由于无法直接利用阵列信号处理相关理论,无疑增加了分布式MIMO雷达建模与算法分析的难度 [2] [3]。

对于分布式MIMO雷达参数估计问题,Lehigh大学的R.S. Blum、H. Godrich以及Q. He较早地研究了分布式MIMO雷达相干模式的目标定位问题 [4],对克拉美–罗界(CRB: Cramer Rao Bound)、最优线性无偏估计量和几何精度因子(GDOP: Geometric Dilution of Precision)进行了分析。J. Li在此基础上推导分析了非相干模式时速度估计的CRB,并给出了参数估计的似然函数 [5],结合最大似然(ML: Maximum Likelihood)估计算法可以实现有效估计,但是文中在参数估计时为简化模型,直接匹配时延因子,这并不够合理。此外,Youssef Fayad结合速度估计的CRB,分析了布站方式对于参数估计性能的影响 [6],张栋和李云成分别通过稀疏重构和非均匀采样实现联合参数估计,并通过仿真验证了均匀布站为相干模式下速度估计的最优布站形式,但是并没有从理论上论证结论的有效性和唯一性 [7] [8]。目标检测一直都是分布式MIMO雷达研究的重点,对于相关工作模式,其检测方法与传统雷达基本相同,难点在于非相干MIMO雷达的检测问题。非相干MIMO雷达的检测问题面临两大挑战,一个是目标回波信号的临界相关问题 [8],一个是雷达通道的信噪比失衡问题,即不同通道中目标回波信号的信噪比不同 [9]。黄勇等提出一种MIMO雷达融合式GLRT检测器 [10],旨在解决含未知参数的高斯杂波背景下的目标检测问题。张新勋等将单基地雷达的GLRT检测器应用于机载MIMO雷达,提出了MIMO-GLRT检测器 [11],但检测器是建立于未知参数的检验框架,解析式中需要已知假设参数,并不符合实际。

本文基于不同的目标模型,推导出MIMO雷达相干处理和非相干处理两种工作模式的数学模型,分析了非相干与相干目标检测与参数估计的性能,得出了非均匀布阵时不同的布站方式以及收发站数对性能的影响,论证了MIMO雷达相较于相控阵雷达对复杂目标检测具有一定优势;且在某一分集自由度时,非相干可逼近相干参数估计的性能,并降低了估计器的复杂度。

2. 分布式MIMO雷达的信号模型

分布式MIMO雷达–目标的空间配置,如图1所示:

Figure 1. Configuration of distributed MIMO radar and target

图1. 分布式MIMO雷达–目标的空间配置图

图中 T k = ( x k t , y k t ) , k = 1 , , M 表示第k部发射站的位置, R l = ( x l r , y l r ) , l = 1 , , N 表示第l部接收站的位置。假设目标由Q个散射体组成,其中第q个散射体的位置为 X q = ( x q , y q ) ,目标的中心记为 X 0 = ( x 0 , y 0 ) 。第k部发射站的发射信号 s k ( t ) 经第q个散射体散射后,被第l部接收站接收的回波信号

可表示为:

r l k q ( t ) = E / M s k ( t τ l k q ) ρ l k q + n l k q ( t ) (1)

式中 ρ l k q = ζ l k q exp ( j 2 π f c ( τ l k q ) ) ,且 τ l k q = ( ( x k t x q ) 2 + ( y k t y q ) 2 + ( x l r x q ) 2 + ( y l r y q ) 2 ) / c 表示的是路

径传输时延。

由于发射信号满足窄带假设,因此各收发路径中到Q个不同散射体的时延,均可近似等于到目标中

心的时延,即 s k ( t τ l k q ) s k ( t τ l k 0 ) = s k ( t τ l k ) ,则第k部发站射发射的信号 s k ( t ) 经目标散射后,被第l

部接收站接收的回波信号可记为:

r l k ( t ) = E / M σ l k s k ( t τ l k ) + n l k ( t ) (2)

式中 σ l k = q = 1 Q ρ l k q

因此,第l部接收站接收到的回波信号表示为所有发射站对应信号的混合,即:

r l ( t ) = E / M k = 1 M σ l k s k ( t τ l k ) + n l ( t ) (3)

假设各发射波形相互正交,对接收的回波信号进行匹配滤波、分离后,可记为:

y l = s l * ( t ) r l ( t ) d t = E / M σ l + w l (4)

式中 y l = [ y l 1 , , y l M ] T , s l ( t ) = [ s 1 ( t τ l 1 ) , , s M ( t τ l M ) ] T ,且 σ l , w l y l 同型。此处假设收发端保证了理

想同步,因而对不同波形进行匹配滤波时可采用不同的时延,以使得匹配滤波器的输出在时间上对齐。

将所有接收站的匹配输出表示为矢量形式,即:

y = H 1 / 0 E / M σ + w (5)

式中 H 1 / 0 = 1 或0,分别表示目标存在或不存在, y = [ y 1 T , y 2 T , , y N T ] T σ = [ σ 1 T , σ 2 T , , σ N T ] T w = [ w 1 T , w 2 T , , w N T ] T 。假设 σ , w 相互独立, w 服从零均值复高斯分布,即 w ~ C N ( 0 , σ w 2 I ) σ 由目标

RCS特性决定。

由于信号 r 可表示为 r = [ r 1 , r 2 , , r M N ] T 的列向量形式,为简化分析,这里省略时延 τ l k ,则相应的信

号模型可表示为:

r ( t ) = E / M Η s ( t ) + n (6)

式中 Η = [ σ 1 , 1 σ 1 , M σ N , 1 σ N , M ] 表示信道矩阵。

因此,结合上式的矢量形式,在对回波信号进行匹配滤波后,且分别以 y w 表示滤波后的回波与噪声,则匹配后可表示为:

y = H 1 / 0 E / M v e c ( H ) + w (7)

式中 v e c ( · ) 表示矩阵按列拉直后的列向量。

值得注意的是,对于分布式MIMO雷达参数估计的匹配滤波,实际上是在估计过程中实现的,即:在估计中匹配以获得真实的目标参数。对于检测而言,为了简化问题,便于理论分析,假设在检测前已完成了匹配滤波,其信号模型如式(7)所示。

3. 分布式MIMO雷达的工作模式

MIMO雷达的回波在不同条件下可能具有相干或者非相干特性,其空间相干性与目标模型、雷达与目标间距所引起的时延、相移等因素有关。当目标为复杂扩展目标时,雷达收发天线位于不同的目标波束宽度内,对应不同通道的接收回波复振幅也相互独立,即:MIMO雷达工作在非相干模式;而当目标为简单目标模型或点目标时,所有天线位于同一目标波束宽度内,空间分集条件不再成立,各通道接收的回波复振幅相等,系统只需对各收发站与目标的间距所引起的时延、相移进行补偿,此时各通道的回波具有强相干性,即:MIMO雷达工作在相干模式。

3.1. MIMO雷达的参数估计性能

针对分布式MIMO雷达的参数估计问题,目前的研究分析了单目标的位置、速度估计性能,推导参数估计的CRB,并基于CRB分析如何优化雷达收发站的配置。

1) 非相干MIMO雷达

Q. He等 [4] 给出了非相干处理时目标速度估计的信号模型,即:

r l ( t ) = k = 1 M E / M ζ l k s k ( t τ l k ) exp ( j 2 π f l k t ) + n l ( t ) (8)

式中, n l ( t ) 表示接收端混入的零均值复高斯白噪声,噪声功率为 σ n 2 f l k 为多普勒频移,是目标位置 ( x , y ) 与速度 ( v x , v y ) 的函数。为便于分析,定义未知参量包括目标速度及散射系数,表示为 2 N M + 2 维的矢量 θ ,即 θ = { v x , v y , ζ 11 R , ζ 11 I , , ζ N M R , ζ N M I } T ,同时定义总回波信号矢量 r = [ r 11 ( t ) , , r N M ( t ) ] T

对于含未知参量 θ 的回波 r ,其概率密度函数记为 f ( r | θ ) 。由于参量满足高斯分布的假设,并取对数后可得到回波信号的似然函数

ln f ( r | θ ) = 1 σ n 2 l = 1 N T | r l ( t ) E / M k = 1 M ζ l k s k ( t τ l k ) exp ( j 2 π f l k t ) | 2 d t + C (9)

CRB是衡量估计均方误差(MSE: Mean Square Error)的标准之一。推导CRB前,需要先得到费歇尔

信息矩阵(FIM: Fisher Information Matrix), J ( θ ) = E { θ ln f ( r | θ ) [ θ ln f ( r | θ ) ] T } ,而CRB即为FIM的逆,即 C C R B = [ J ( θ ) ] 1

利用链式法则(Chain rule),通过矩阵和代数运算,取上式矩阵对角线的前 2 × 2 维分块矩阵,即未知参量 θ 中的速度分量,则 C C R B v e l 可表示为:

C C R B v e l = ( U D U H U G L 1 ( U G ) H ) 1 = [ C R B v x C R B v y ] (10)

式中矩阵的具体形式由多个参量确定,例如载波波长、发射站的能量、噪声的功率谱密度、目标反射系数、感兴趣的距离单元、发射信号以及隐含在多普勒函数中的收发站数目和位置等等。

2) 相干MIMO雷达

A. De Maio等 [3] 研究了相干MIMO雷达位置估计的信号模型,即:

r l ( t ) = k = 1 M E / M ζ s k ( t τ l k ) exp ( j 2 π f c τ l k ) + n l ( t ) (11)

式中 f c 为发射信号的载频。定义未知参量 θ 包括位置、散射系数,即 θ = [ x , y , ζ r e , ζ i m ] T 。由于目标位置可由时延 τ l k 表示,利用链式法则,FIM可记为 J ( θ ) = ψ / θ J ( ψ ) ( ψ / θ ) T

对于含未知参量 ψ 的观测信号 r ,其概率密度函数服从高斯分布,即 f ( r | ψ ) 表示为:

f ( r | ψ ) exp { 1 σ n 2 l = 1 N T | r l ( t ) E / M k = 1 M ζ s k ( t τ l k ) exp ( j 2 π f c τ l k ) | 2 d t } (12)

结合FIM的定义,给出以 ψ 为参量的FIM解析形式为:

J ( ψ ) = [ Λ x y V V T Θ ζ ] (13)

式中 Λ x y = 8 π 2 S N R ( f c 2 + β 2 ) I Θ ζ = S N R N M | ζ | 2 I V = 4 π f c S N R | ζ | 2 [ ζ i m ζ r e ζ i m ζ r e ] 。信噪比定义为 S N R = | E / M ζ | 2 / σ n 2 ,有效带宽 β = w f 2 | S ( f ) | 2 d f / w | S ( f ) | 2 d f ,W为信号带宽。

由于参数估计中只对目标位置参量感兴趣,表示为 C R B x = C C R B ( 1 , 1 ) , C R B y = C C R B ( 2 , 2 ) ,即: C C R B 矩阵对角线的前 2 × 2 分块矩阵元素。

C C R B l o c = [ C R B x C R B y ] (14)

在加性噪声的问题中,当SNR足够大时,ML估计的均方误差可近似达到CRB。ML估计器可表示为:

θ ^ M L = arg { max θ { ln f ( r | θ ) } } (15)

3.2. MIMO雷达的目标检测性能

分布式MIMO雷达的检测与传统雷达检测的根本区别在于:传统的相控阵雷达和组网雷达(netted radar)是通过相干信号处理,最大化相干处理增益来获得较好的检测性能,而分布式MIMO雷达主要是通过非相干信号处理,获得空间分集增益来改善目标的检测性能。

1) 非相干MIMO雷达

E. Fishler等 [2] 给出了非相干检测时的充分统计量为:

T = y H y = i = 1 M N | y i | 2 (16)

式中 y i 是向量 y 的第i个分量,进一步得到充分统计量的分布特性,即:

{ H 0 : T ~ 1 2 σ w 2 χ 2 M N 2 H 1 : T ~ 1 2 ( σ w 2 + E M ) χ 2 M N 2 (17)

此时虚警概率和检测概率分别表示为:

{ P f a = P ( T γ / H 0 ) = P ( χ 2 M N 2 2 γ σ w 2 ) = Q χ 2 M N 2 ( 2 γ σ w 2 ) P d = P ( T γ / H 1 ) = Q χ 2 M N 2 ( 2 γ σ w 2 ( 1 + S N R / M ) ) (18)

式中 S N R = E / σ w 2 表示信噪比形式, γ 表示由虚警概率 P f a 确定的检测门限, χ 2 M N 2 表示自由度为 2 M N 的中心式 χ 2 分布, P ( ) 表示概率分布函数, Q χ v 2 ( x ) = x p ( t ) d t p ( t ) χ v 2 的概率密度。

2) 相干MIMO雷达

假设目标只有一个散射体(点目标),且各向同性,假设已知目标参数,即对接收的回波信号进行相位

补偿,即 r l k ( t ) = E / M ζ s k ( t τ l k ) + n l k ( t ) ,则第l部接收站接收到的回波可表示为:

r l ( t ) = E / M k = 1 M ζ s k ( t τ l k ) + n l ( t ) (20)

式中 ρ l k = ζ exp ( j 2 π f c ( τ l k ) ) ζ ~ C N ( 0 , 1 ) 。匹配滤波后,接收站输出可写成矢量形式,即:

y = E / M ζ 1 M N × 1 + w (21)

对上式求对数似然比,进一步简化后得对数似然比,即:

ln f ( y / H 1 ) f ( y / H 0 ) = η 0 + E / M | y H 1 M × 1 | 2 k (22)

式中 η 0 与k均为与观测数据无关的常数。设充分统计量为:

T = | y H 1 M N × 1 | 2 = | i = 1 M N y i | 2 (23)

此时分布式MIMO雷达的检验统计量是各通道输出相干求和后的能量,则分布特性可记为:

{ H 0 : T ~ 1 2 M N σ w 2 χ 2 2 H 1 : T ~ 1 2 M N ( σ w 2 + E M M N ) χ 2 2 (24)

虚警概率和检测概率可分别表示为:

{ P f a = P ( T γ / H 0 ) = P ( χ 2 2 2 γ M N σ w 2 ) = Q χ 2 2 ( 2 γ M N σ w 2 ) P d = P ( T γ / H 1 ) = Q χ 2 2 ( 2 γ σ w 2 M N ( 1 + N S N R ) ) (25)

4. 仿真实验与数据分析

仿真中,以发射和接收站数目 M = 2 , N = 4 M = 4 , N = 6 为例。

1) MIMO雷达非相干、相干处理的参数估计

为简化分析,这里假设收发站与目标的距离均近似相等,即:对任意 l , k ,都有 R l = R k = 7000 m 。在此基础上,系统随机布阵,即 φ t = [ 80 , 120 ] φ r = [ 199 , 220 , 245 , 300 ] φ t = [ 110 , 100 , 51 , 79 ] φ r = [ 325 , 60 , 300 , 265 , 115 , 85 ] 。假设空间有1个目标,坐标位置为 ( 50 , 30 ) m ,速度为 ( 15.5 , 97.1 ) m / s 。发射信号为正弦脉冲信号,其数学形式为 s k ( t ) = e j 2 π k Δ f t / T ,其中频率增量 Δ f = 200 Hz 以保证正交性,信号时间宽度 T = 10 ms ,载频 f c = 1 GHz ,采样频率 f s = 140.625 kHz

图2所示,随着天线数增加,参数估计的性能迅速提高,且估计中Y方向的效果要优于X方向,这是由于布站配置中天线相对于Y方向投影的延展更大,因此Y方向对估计效果更敏感,参数估计的性能更高。且由于检测器的性能发挥会直接作用于参数估计,因此相对于传统的参数估计值而言,本文算法的估计性能会略微下降,但是更符合实际情况。

Figure 2. CRB of parameters estimation of distributed MIMO radar

图2. 分布式MIMO雷达参数估计的CRB曲线

根据上面的分析,可得如下结论:非均匀布阵中,阵列的布阵方式会影响参数估计的精度,收发站沿某方向投影所得到的延展越大,则该方向对应的估计精度越高。

2) MIMO雷达非相干、相干处理与相控阵雷达的检测性能

这里设虚警概率 P f a = 10 6 ,如图3所示。由图3(a)可以看出,在高检测概率时(近似 P d 0 .9 )或高SNR时,MIMO雷达的非相干检测性能优于相控阵雷达,而在低检测概率时(或低SNR时)相控阵雷达的检测性能优于MIMO雷达的非相干检测,这说明在高SNR时,通过空间分集获得的增益优于相干处理增益;而在低SNR时,相干处理增益比空间分集增益对雷达检测性能的改善更多。

MIMO雷达的相干模式的检测性能比相控阵的检测性能差,这是由于在相同的发射功率和回波接收功率条件下,MIMO雷达采用了波形分集技术,使得接收端的通道数目增加,引入了更多的噪声,使得总SNR下降,因而MIMO雷达的相干检测性能变差。分别增加2个发射站和2个接收站,重新进行仿真,检测性能如图3(b)所示,上述结论仍然成立。

根据上面的分析,可得如下结论:相控阵的检测性能优于MIMO雷达的相干检测,但在高SNR时,MIMO雷达的非相干检测性能要优于相控阵雷达。

Figure 3. The Detection of distributed MIMO radar and phased array radar

图3. 分布式MIMO雷达和相控阵雷达的检测性能曲线

5. 结束语

本文首先介绍了分布式MIMO雷达的信号模型,包括参数估计与检测的具体信号形式,并根据不同的目标模型对MIMO雷达工作模式进行了分类,即:相干处理和非相干处理。构建了相应的数学模型,并推导出估计器与检测概率的最终表达式,论证了不同的布站方式以及收发站数对参数估计和目标检测性能的影响。通过仿真表明:非均匀布站的方式会影响参数估计精度,收发站沿某方向投影所获得的延展越大,该方向的估计精度越高;在高SNR时,MIMO雷达的非相干检测相对于相控阵具有更好的检测性能,即MIMO雷达在检测方面对复杂目标具有一定优势,这与参数估计的结论相反。

文章引用

马 鹏,王石泉,赵 飞. 分布式MIMO雷达参数估计与目标检测性能研究
Parameter Estimation and Target Detection for Distributed MIMO Radar[J]. 测绘科学技术, 2021, 09(04): 146-154. https://doi.org/10.12677/GST.2021.94018

参考文献

  1. 1. Bliss, D.W. and Forsythe, K.W. (2003) Multiple-Input Multiple-Output (MIMO) Radar and Imaging: Degrees of Free-dom and Resolution. Proceeding 37th Asilomar Conference on Signal, System and Computers, 1, 54-59.

  2. 2. Fishler, E., Haimovich, A.M., Blum, R.S., Cimini, L., Chizhik, D. and Valenzuela, R. (2004) Performance of MIMO Radar Systems: Advantages of Angular Diversity. Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems & Computers, 1-2, 305-309.

  3. 3. De Maio, A. and Lops, M. (2007) Design Principles of MIMO Radar Detectors. IEEE Transactions on Aerospace and Electronic Systems, 43, 886-898. https://doi.org/10.1109/TAES.2007.4383581

  4. 4. He, Q., Lehmann, N., Blum R.S. and Haimovich, A.M. (2010) MIMO Radar Moving Target Detection in Homogeneous Clutter. IEEE Transactions on Aerospace and Electronic Sys-tems, 46, 1290-1301. https://doi.org/10.1109/TAES.2010.5545189

  5. 5. Li, J. and Stoica, P. (2009) Concepts and Applications of a MIMO Radar System with Widely Separated Antennas. In: Li, J. and Stoica, P., Eds., MIMO Radar Signal Processing, John Wiley & Sons, Inc., Hoboken, 365-410. https://doi.org/10.1002/9780470391488

  6. 6. Fayad, Y., Wang, C.Y. and Cao, Q.S. (2017) Temporal-Spatial Sub-spaces Modern Combination Method for 2D-DOA Estimation in MIMO Radar. Journal of Systems Engineering and Electronics, 28, 34-39. https://doi.org/10.21629/JSEE.2017.04.09

  7. 7. 张栋, 白新有, 邱彦文, 郑桂妹. 基于联合稀疏重构模型的双基MIMO雷达off-grid目标参数快速估计[J]. 现代雷达, 2020, 41(11): 30-34.

  8. 8. 李云成, 崔琛, 龚阳. 基于平行因子分析的时空非均匀采样下双基地MIMO雷达目标参数联合估计算法[J]. 电子工程学院学报, 2019(3): 1-7.

  9. 9. Cui, G. Kong, L. and Yang. X. (2010) Multiple-Input Multiple-Output Radar Detectors Design in Non-Gaussian Clutter. IET Radar, Sonar & Navigation, 4, 724-732. https://doi.org/10.1049/iet-rsn.2009.0056

  10. 10. 黄勇, 黄涛, 刘宁波, 陈小龙. 基于均匀性判定规则的统计MIMO雷达多通道融合检测技术[J]. 海军航空工程学院学报, 2018, 33(1): 101-105+118.

  11. 11. 张新勋, 周生华, 刘宏伟. 频率分集极化MIMO雷达目标检测[J]. 华中科技大学学报(自然科学版), 2019, 47(6): 11-16+22.

期刊菜单