一种MIMO雷达更窄主瓣波束形成方法

李朝海，李 明，张 伟

(电子科技大学电子工程学院 成都 611731)

.通信与信息工程.

一种MIMO雷达更窄主瓣波束形成方法

李朝海，李 明，张 伟

(电子科技大学电子工程学院 成都 611731)

介绍了一种获得更窄主瓣虚拟收发波束的方法。该方法通过天线阵元发射正交编码的波形，然后在接收端对回波信号进行处理。MIMO雷达由于靠稀疏阵列避免旁瓣的不利影响，所以可提高角分辨率和多普勒分辨率。但受到如空间、成本或者其他实际因素的限制，须将阵列布置成紧凑阵列，而紧凑阵列对低速目标检测效果较差。通过提出的加窗方式，紧凑阵列利用虚拟收发波束形成可获得更窄的主波束，提高对低速目标的检测效果。通过仿真验证，对归一化频率为0.062 5的低速目标，信干噪比提升系数(SINR IF)约提升7.3 dB。

紧凑阵列; MIMO雷达; 稀疏阵列; 加窗方式

MIMO雷达相比于传统雷达有许多优势[1-3]，如具有更高的角分辨率和多普勒分辨率等。当MIMO雷达发射正交信号[4-7]，发射模式和接收模式均可进行阵列处理，在接收端通过匹配滤波进行信号提取，提取到的信号可用于提高杂波的空间分辨率。利用不同发射天线和不同接收天线的相位偏差可构成一个新的虚拟阵列导向矢量，合理设计天线位置，MIMO雷达利用少量的天线就可获得一个非常长的严格采样的导向矢量数组[8]。

MIMO雷达使用稀疏阵列可避免旁瓣带来的不利影响[9]。但设置成稀疏阵列需要更大空间，制约了在机载雷达上的应用，且发射和接收不能使用相同天线，增加了成本。若使用紧凑阵列则可以很好地解决这些问题。本文利用新的算法使用紧凑阵列获得比传统波束形成方案更窄的主波束。

1 MIMO雷达信号模型

假设有TN个发射单元，间隔均为Td，有RN个接收单元，间隔均为Rd。采用单站式MIMO雷达，在每个发射单元，一个相关处理间隔发射脉冲的个数为L个。假设发射信号相互正交，系统内部没有杂波干扰，雷达系统平台的运动速度为v[10]。发射天线和接收天线有相同的方位角θ和俯仰角ϕ。设为波长。在每个脉波重复间隔(PRI)rT，系统发射波束集合为：

式中，γ=dT/dR；fd=2vsinθcosϕ⋅Tr/λ是移动平台的归一化多普勒频率；fs=dRsinθcosϕ/λ是空间频率；ξ(θ)是在θ处的反射系数。假设dR设置合理，保证−0.5≤fs≤0.5，从而避免空间频率混叠。

本文假设发射导向矢量和接收导向矢量分别为：

定义多普勒矢量为：

对于第l个脉冲，接收单元的杂波回声为：

式中，Yl∈CNR×1。接收端信号处理是为了分离不同信号的相对时间延迟特性，从而为后续处理提供信号向量。每个接收到的信号包含所有发射波形，设sH为接收端的滤波器组，有：

式中，X∈CNR×NT。使用向量形式保持第l个脉冲的过滤杂波数据：

式中，x∈CNRNT×1；“⊗”表示Kronecker积。距离环l被分为NC(NC>>NTNRL)个杂波片，那么第l个脉冲的杂波数据分离形式为：

然后将第L个脉冲的整个杂波回波保存为向量形式：

2 更窄主瓣波束形成方法

接收波束通过相位变换叠加处理来利用天线阵元的输出实现延迟叠加波束形成，从而实现数字合成[12]。MIMO雷达的正交信号可在接收端通过匹配滤波进行提取，提高杂波的空间分辨率[13]。维数为NTNR的阵列导向矢量可以等效为发射导向矢量与接收导向矢量的卷积。假设NR=NT=N ，对于波束形成，只考虑等间距的线性阵列。已经证明在匹配滤波之后实际获得的天线波束相比于对应的传统相控雷达具有相同的天线增益和空间分辨率。同时数字滤波器的系数矢量按如下方式设置：

式中，κ满足归一化的要求；s0是来自期望方向θ0处的目标信号，有：

通过式(12)和式(13)，可以获得具有最大信噪比(SNR)的滤波器输出。而收发波束的指向可通过数字滤波器获得。将式(12)改写为：

式中，κ=κ1κ2，均为常数；wT为发射波束的数字滤波器的系数向量；wR为接收波束的权向量。利用式(12)和式(14)中的最佳滤波器，通过虚拟收发波束形成，可获得虚拟天线增益：

上式为离散傅里叶变换(DFT)[14]的形式。本文设置dR=dT=d，d/λ=0.5，NR=NT=N ,可得：

从式(15)和式(16)知，G(θ,θ0)可以看做是两个长度均为N的方波离散傅里叶变换(DFT)的积。图1阐述了一维均匀线性阵列导向矢量的卷积性质。

MIMO雷达的虚拟收发波束形成可以看作是两个方波的离散傅里叶变换的积。两个方波的卷积是三角波，其长度为2N−1，如图1所示。图1中的振幅是具有相同的相位中心的阵元的数量。根据傅里叶理论，发射波束的主瓣宽度被扩宽了。为了更明确这一点，改写MIMO雷达接收端经过匹配滤波后的信号的矩阵如下：

式中，γ=dT/dR，vT为发射导向矢量；vR为接收导向矢量，它们均已在式(3)和式(4)中进行了定义。如果考虑到噪声，则匹配滤波后的输出信号为：

可以发现式(19)中有很多元素相等，如图2所示。

根据离散傅里叶变换的性质，相同长度的方波主瓣比三角波更窄。在直觉上可以通过权值反转变换，将三角波变换为方波，如图3所示。

与式(18)或式(19)对应的权值转换矩阵为：

d=stack(D)。对应于虚拟收发波束的数字滤波器系数向量为：

如果使用窗函数，由于权向量wnew与期望信号s不能完美匹配，所以SNR有一定的损失。根据式(18)中的一维波束形成加窗的信噪比损失因子，定义虚拟收发波束加窗的信噪比损失因子为ηeff，它是加窗后的信噪比(SNRtap)与最优输出信噪比(SNRopt)之比值。为便于计算，ηeff通常定义为：

式中，s0是期望信号，其定义如式(13)。考虑锥形的情况，wnew=s0*d，保证单位增益约束为wHw=1，vec(t)必须满足：

图4为信噪比损失因子ηeff与阵元数目变化关系。从图4可以发现，随着阵元数目的增加，ηeff的减小速率将会降低。对应于16阵元、32阵元、80阵元，其ηeff分别为−2.515、−3.193、−3.998 dB。更窄的主瓣会带来SNR的损失，可以根据需要合理选择。更窄的波束可以提高角分辨率和多普勒分辨率，抑制来自目标的主瓣杂波干扰。通常，接收端的剩余旁瓣会带来足够的干扰掩盖期望的返回信号，在波束形成中是不允许存在过高的旁瓣。本文也可以使用切比雪夫窗函数。直接将切比雪夫窗函数式(21)写为：

将式(25)对应的加窗方式命名为方式1。根据式(20)中提出的加权方式，该加权方式将会改变切比雪夫窗函数的形状，旁瓣可能得不到有效地抑制。可以通过改变切比雪夫窗函数来克服该问题。根据图4描述的相位中心，将wch,2N−1扩展到N2个点，wch,2N−1是2N−1个点的切比雪夫窗函数的权向量。即由相同相位的匹配滤波器有相同的切比雪夫窗权系数为：

将式(26)对应的加窗方式命名为方式2，将式(26)代入式(25)，替换

可以得到虚拟天线增益为：

为了评估减小的主瓣杂波，虚拟收发波束形成被配置成串联T-S模式。时间权向量首先在每个天线通道的每个匹配滤波器的输出上进行计算，有：

式中，wt,l是对应于离散傅里叶变换系数的权向量；下标l表示第l个多普勒通道；下标t表示时域。多普勒处理器的输出结果在空域中进行进一步处理以获得检验统计量：

式中，下标s表示时域，有：

通过渐变模式得到的更窄的主瓣波束有很多优势，比如提升角分辨率和多普勒分辨率，具有更尖锐的机载雷达杂波凹口。这些优势都可以应用于机载雷达对地面低速运动目标的检测。通常用单阵元、单脉冲的SINR值来衡量空时处理算法的性能。SINRin作为单阵元、单脉冲的SINR。通常情况下杂波和干扰比较强，所以SINRin相对是一个非常小的量。SINR的提高因子IFSINR被定义为：

SINRout(ϖ)在串联T-S模式下的配置为：

式中，st,l(fd)是期望目标的多普勒导向矢量：

Rl是第l个多普勒通道的协方差矩阵，有：式中，I是单位矩阵；R是整个协方差矩阵，有：

3 仿真结果

考虑一八阵元发射接收阵列，NT=NR= N=8，dT=dR=d=λ/2，主瓣方向θ0=0°，目标位于主瓣所在方向。载波频率为f0=1.2 GHz ，脉冲重复频率PRF为2 000 Hz，一个相关处理间隔发射脉冲数为M=64。雷达运动速度为v=125 m/s ，对于每一个天线阵元，每一个脉冲内的SNR=−5 dB，CNR=45 dB。式(25)对应方式1的波束形成方式，式(27)中对应方式2的波束形成方式，以及只考虑在接收波束形成中使用30 dB切比雪夫加窗(方式3)的相控阵雷达波束。对比3种不同的波束形式如图5所示。

由图5可知，由于方式1的切比雪夫窗函数形式的改变并没能有效抑制旁瓣。而方式2的切比雪夫加窗过程由于采用虚拟的发射接收波束，旁瓣得到了有效的抑制。对于相控阵雷达系统，由于考虑到最大发射功率的要求，切比雪夫加窗只用在接收波束形成过程中，虽然旁瓣得到了抑制，但主瓣也相应的展宽了。相比于其他方法的信干噪比改善因子如图6所示。为方便比较，在归一化多普勒频率处的最小的改善因子也被归一化为0 dB。

图6中所加切比雪夫窗为30 dB，本文提出的基于渐变处理的波束形成方法主要在小多普勒频率范围优于没有进行渐变处理的波束形成方法，这一点对于检测低速运动目标非常有用。比如对于低速运动目标，在归一化多普勒频率为0.062 5时，新方法的信干噪比改善因子可以达到7.3 dB。没有进行切比雪夫加窗时，新方法的信干噪比改善因子比没有进行渐变处理时在旁瓣区域差，但可以通过简单地修正切比雪夫窗函数来克服不足。在小多普勒频率区域，用30 dB切比雪夫加窗形成的相控阵雷达接收波束效果比较差，特别是在归一化多普勒频率位于(−0.1,0.1)时。通过上面的分析可知，波束主瓣宽度越窄，对低速目标的检测效果越好。

4 总 结

本文提出的MIMO雷达匹配滤波器输出信号的幅度经过加窗处理之后，对低速运动物体的检测效果有一定程度的提高。通过本文所提出的加窗方式可以获得更窄的主瓣，显著提高杂波抑制效果，改善角度分辨率，提升探测低速运动目标性能，因此可以提高对低速运动物体的检测效果。对于其他阵元数目的MIMO雷达，采用该加窗方式之后检测效果类似。方便将该方式应用于级联时空处理。本文的方法也可以应用于其他波束域的STAP算法。这种新型加窗方式提高了对低速运动目标的检测效果，完善了MIMO雷达的性能，对于现代机载雷达具有重要意义。

[1] BROOKNER E, RAYTHEON C. MIMO radar demystified and where it makes sense to use[C]//Phased Array Systems&Technology. Waltham, USA: IEEE, 2013.

[2] MOO P W, ZHEN Ding. Tracking performance of MIMO radar for accelerating targets[J]. Signal Processing, 2013, 61(21): 5205-5216.

[3] RANKIN G A, TIRKEL A Z, BUI L Q, et al. Radar imaging: Conventional and MIMO[C]//Communications and Electronic(ICCE). Hue, USA: IEEE, 2012.

[4] HE Qian, LEHMANN N H, BLUM R S, et al. MIMO radar moving target detection in homogeneous clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(3): 1290-1301.

[5] 牛朝阳, 张剑云. MIMO雷达正交波形集设计—线性调频-位编码混合波形[J]. 计算机工程与应用, 2012, 48(13): 133-137. NIU Zhao-yang, ZHANG Jian-yun. Orthogonal waveform set of hybird LFM and phase code for MIMO radar[J]. Computer Engineering and Application, 2012, 48(13): 133-137.

[6] RABIDEAU D J. Non-adaptive multiple-input, multipleoutput radar techniques for reducing clutter[J]. IET Radar Sonar Navig, 2009, 3(4): 304-313.

[7] 杜晓林, 王旭, 苏涛, 等. MIMO雷达正交相位编码信号的代数设计方法[J]. 西安电子科技大学学报(自然科学版), 2014, 41(2): 64-70. DU Xiao-lin, WANG Xu, SU Tao, et al. Orthogonal phase coding signal design for MIMO radar via an algebraic method[J]. Journal of Xidian University(Natural Science Edition), 2014, 41(2): 64-70.

[8] CHEN Chun-yang, VAIDYANATHAN P P. Beamforming issues in modern MIMO radars with doppler[C]//The 40th Asilomar Coference on Signals, Systems and Computers. Pacific Grove, CA, USA: IEEE, 2006.

[9] FORSYTHE K W, BLISS D W. MIMO radar signal processing[M]. New York, USA: Wiley, 2009.

[10] DENG Hai, BRAHAM H. A virtual antenna beamforming (VAB) approach for radar systems by using orthogonal coding waveforms[J]. IEEE Transactions on Antennas and Propagation, 2009, 57(2): 425-435.

[11] WANG Guo-hua, LU Yi-long. Clutter rank of mimo radar with a speical class of waveforms[C]//Waveform Diversity and Design Conference. Kissimmee, FL, USA: IEEE, 2009.

[12] MONDAL B, HEATH R W. Performance analysis of quantized beamforming MIMO systems[J]. IEEE Trans Signal Process, 2006, 54(12): 4753-4766.

[13] CHEN Chun-yang, VAIDYANATHAN P P. MIMO radar space-time adaptive processing using prolate spheroidal wave functions[J]. IEEE Transactions on Signal Processing, 2008, 56(2): 623-635.

[14] 周剑雄, 石志广, 吴京. 离散时间信号频域分析中的概念辨析[J]. 电气电子教学学报, 2014, 36(1): 49-51. ZHOU Jian-xiong, SHI Zhi-guang, WU Jing. Conceptual discrimination in frequency domain analysis of discrete time signals[J]. Journal of EEE, 2014, 36(1): 49-51.

编 辑 税 红

An Approach to Achieving Narrower Mainlobe for MIMO Radar

LI Chao-hai, LI Ming, and ZHANG Wei
(School of Electronic Engineering, University of Electronic Science and Technology of China Chengdu 611731)

An innovative approach is introduced to form narrower mainlobe of virtual transmitting-receiving beamforming by transmitting orthogonal coding waveforms from the antenna elements and digitally processing their echoes at the receiver. The benefits such as increasing angle resolution and increasing Doppler resolution of MIMO radar mainly depend on the sparse arrays without the adverse effects of sidelobes. However, due to the actual limitation, such as room limitation for airborne, cost or other factors, we have to set the array in dense arrays configuration which is not good at detecting slow moving target. With the proposed amplitude taper approach, a narrower main beam can be achieved with dense arrays configuration after virtual transmitting-receiving beamforming and therefore the performance of detecting slow moving target is improved. Simulation result shows that for the normalized Doppler frequency 0.062 5, i.e., slow moving target, the signal to interference and noise improvement factor (SINR IF) is increased about 7.3 dB.

dense arrays; MIMO radar; sparse arrays; tapering

TN911.3

A

10.3969/j.issn.1001-0548.2015.03.001

2014 − 04 − 25；

2015 − 01 − 06

国家自然科学基金(61371184, 61201280, 61401062, 61301262)

李朝海(1972− )，男，高级工程师，主要从事相控阵雷达与自适应信号处理方面的研究.