干涉仪对频分正交MIMO雷达的测向性能分析

摘 要：      频分正交MIMO雷达的空间合成信号脉冲具有时变包络和空时耦合特性， 为主动对抗无源电子侦察系统的测向提供了新的可能。 本文针对相位干涉仪测向系统分析了电子侦察系统对频分正交MIMO雷达的测向误差。 首先建立远场侦察接收机MIMO雷达合成信号模型， 分析其时变包络特性和相位波前畸变特性及其对测向性能的影响， 并提出一种通过初相调制来增大干涉仪测向误差的可行路径。 仿真结果表明， 相位波前畸变区和瞬时信噪比影响的干涉仪DOA估计误差主导性区域不同。 利用MIMO雷达阵元相位调制可操控包络起伏特性， 改变合成信号峰均比（PAR）从而增大抗截获性能。 在干涉仪分析窗内， 与相控阵雷达合成信号相比， 当侦察接收机先验位置已知时， MIMO雷达合成信号功率可下降25 dB左右， 而距离未知时也可下降10 dB左右， 因此MIMO雷达能在一定程度上增大干涉仪的测向误差。

关键词：     MIMO雷达； 相位调制； 相位波前畸变； 相位干涉仪； 测向误差

中图分类号：       TJ760; TN958

文献标识码：    A

文章编号：     1673-5048（2024）01-0103-08

DOI： 10.12132/ISSN.1673-5048.2022.0207

0 引  言

现代战争中， 制电磁权的核心就是侦察与反侦察的对抗， 随着电子战技术的发展， 雷达的生存面临巨大的威胁。 雷达对抗侦察技术一般分为两种， 一种是防御性的， 即通过功率管理、 频率捷变、 复杂调制、 大带宽等低截获概率（LPI）技术使得非合作侦察接收机对信号的检测、 识别和参数估计变得困难； 另一种是主动性的， 即设置干扰源实施噪声干扰， 或者发射与雷达信号相似的假辐射源信号， 从而掩护己方雷达信号［1–3］。 近年来， 一种具备干扰、 侦察、 定位等多种功能的微型空射诱饵（MALD）快速发展， 战时可由防区外电子战飞机发射进入对方防空区域抵近侦察［4］， 雷达信号将不可避免被截获。 由于辐射源到达角是侦察系统信号分选和定位的一个重要参数， 而准确定位则是进一步完成电子软硬杀伤的前提， 如果能有效降低电子侦察系统的测向精度， 必然能提升我方雷达的战场生存能力。

MIMO雷达是一种新体制雷达， 各阵元可同时发射不同波形， 具有较高的设计自由度， 从主动对抗侦察机的角度看， 在MIMO雷达信号特性层面展开分析， 研究空间合成信号对电子侦察系统测向定位的误差影响， 进而通过波形设计和相位调制的方式增大电子侦察系统的测向定位误差， 同时不影响己方雷达探测性能， 具有一定的可行性。 现阶段， 主动对抗无源定位系统的研究主要是利用相干干扰进行的［5–8］， 通过雷达波形本身的辐射特性主动对抗电子侦察系统测向定位的研究较少。 文献［9–12］主要研究了侦察机对具有较小频偏的频控阵雷达的测向定位误差， 表明阵列雷达可以通过不同频控函数有效增大定位误差。 文献［13］研究了OFDM-MIMO雷达信号的反侦察性能， 分析了MIMO雷达的脉内信号时变特性对侦察机测向的影响。 文献［14-15］研究了MIMO雷达信号抗到达时差定位的效果， 表明MIMO雷达的反定位效果要优于相控阵雷达。 然而， 这些研究并没有深入研究MIMO雷达发射加权参数对侦察机测向欺骗的影响效果， 未考虑多阵元信号空间合成信号的相位波前畸变特性。

一般的测向方法主要针对单载频恒包络脉冲信号， 但频分正交MIMO雷达远场合成信号将具有时变包络特性和相位波前畸变特性， 这一特性将影响电子侦察系统对MIMO雷达信号DOA参数的测量结果。 因此， 本文以干涉仪对频分正交MIMO雷达DOA测量为例， 在发射信号能量约束下建立MIMO雷达远场合成信号及其幅度和相位的空时分布模型， 在此基础上进一步分析DOA测量误差规律以及通过MIMO雷达信号初相调制增大相位干涉仪测向误差的可能性。

1 MIMO雷达信号模型

MIMO雷达的每个阵元发射相互正交的信号， 不失一般性， 考虑由M个偶极子天线构成的MIMO雷达系统。 其一维均匀线阵（Uniform Linear Array， ULA）结构如图1所示， 图中， s1， s2， …， sM为各阵元的发射信号； d为阵元间距。

第m个阵元信号可表示为

发射信号之间满足

式中： TP为脉冲宽度； Es为信号能量。

以MIMO雷达为坐标原点， 考虑在雷达天线远区（即r≥2D2/λ， D为雷达孔径， λ为雷达波长）， 波阵面局部可用平面波表示。 在窄带条件下， MIMO雷达阵列孔径渡越时间远小于MIMO信号带宽的倒数， 因此可忽略不同阵元发射信号的包络延迟差异及其带来的附加相位差， 此时， 观测点（r， θ）处的MIMO雷达远场合成信号可表示为

式中： wm是第m个阵元的初相控制系数。

采用欧拉公式法［16］， 可得MIMO雷达远场合成信号的幅度分布函数和相位分布函数， 即

令每个阵元的信号功率为Pt=Es/TP， 归一化距离衰减因子后， 由式（1）～（3）可得， MIMO雷达信号的空间平均功率分布为

由式（7）可知， 正交MIMO雷达信号的空间平均功率不随角度变化， 表现为宽波束。

当各阵元的发射信号完全相同时， 即s1=s2=…=sM， 此时对应相控阵雷达的工作模式， 代入式（7）可得相控阵雷达平均辐射功率的空间分布：

由式（8）可知， 相控阵雷达发射的是高增益的窄波束， 其平均辐射功率是关于角度的函数。

2 频分正交MIMO雷达合成信号空时分布特性

MIMO雷达的正交波形可以从时域、 频域或码域上设计［17］。 本文考虑线性频分正交MIMO雷达信号， 即将同一信号包络调制到等间隔载频上的一组正交信号。 为简化分析， 考虑矩形包络脉冲， 记参考阵元的信号载频为fc， 相邻阵元的信号频差为fΔ（为了保证正交性， fΔ应为1/TP的整数倍）。 初相控制系数wm一般可取任意相位调制序列或如式（9）所示的相位系数。

式中： θa为参数， 类似相控阵雷达的波束扫描控制角。

2.1 频分正交MIMO雷达合成信号幅度特性

由式（10）可知， 频分正交MIMO雷达的幅度分布为辛克函数形状， 且具有时变特性。 进一步可得归一化的脉内峰值位置：

由式（11）可知， 幅度峰值的相对周期为1/fΔTP。 当fΔTP=1时， 幅度分布函数只有一个峰值。 通过参数θa可以灵活调控峰值位置， 例如， 需要将峰值置于信号中部， 令tpeak/TP=1/2便可解出所需的控制参数θa。

当MIMO雷达处于相控阵工作模式时， 即fΔ=0 Hz， 代入式（10）可得相控阵雷达远场合成信号的幅度分布：

由式（12）可知， 相控阵雷达的远场信号幅度只与观测角度有关， 且不同于正交MIMO雷达全向发射的工作方式， 相控阵雷达的发射波束需要通过相位控制实现角度扫描。 另外， 与式（10）对比可知， 由于式（12）中与时间参量t无关， 因此相控阵雷达远场合成信号脉冲包络表现为恒包络特性。

2.2 频分正交MIMO雷达合成信号相位特性

由式（5）和式（6）可得， 频分正交MIMO雷达合成信号的相位分布函数为

等相位面是某一时刻空间中相位相等点构成的平面。 类似于相干多点源干扰对抗单脉冲测角的原理。 多阵元矢量合成信号能够使等相位面在某点发生较大的畸变［18-19］， 称为相位波前畸变现象。 相位梯度垂直于等相位面， 其定义式为

式中： r和θ分别表示距离维和角度维的方向单位矢量。

由于相位干涉仪实际测量的是相位梯度的方向。 按照相位梯度法的基本原理［20］， 考虑位于平面内距离MIMO雷达任意一点r， θ处， 可用该处的相位梯度矢量与距离矢量的夹角表征该点的相位波前畸变程度。 因此， 定义相位梯度偏差角eθ为

图2给出了相位波前畸变以及相位梯度偏差角的示意图。 当eθ不为0时， 表示相位梯度不平行于距离方向（入射方向）矢量， 即发生了相位波前畸变。 将式（13）代入式（15）， 可得

式（16）表明， 频分正交MIMO雷达存在相位波前畸变现象， 其相位梯度偏差角大小与阵元间距d、 阵元数M、 观测点角度θ呈正相关， 与fΔ/fc以及观测点距离r呈负相关。

需要注意的是， 当MIMO雷达处于相控阵工作模式且波束扫描角为θa时， 相控阵雷达的空间相位分布函数为

将式（17）代入式（15）， 得

对比式（16）和式（18）可知， 由于式（16）多了fΔ/fc项， 因此频分正交MIMO雷达的相位梯度偏差角小于相控阵雷达。

3 干涉仪对频分正交MIMO雷达的测向误差

3.1 干涉仪测向原理

机载无源侦察系统通常采用相位干涉仪法来估计辐射源DOA。 如图3所示， 相位干涉仪的基本结构具有2个天线， 其中l为干涉仪基线长度。

与干涉仪阵元输出相连的窄带接收信道， 其采样输出可表示为

式中： r=2πlsinα/λ是附加相移； x［n］是干涉仪处MIMO雷达合成信号的离散采样序列； v1［n］和v2［n］是均值为零、 方差为σ2n的加性高斯白噪声序列。

通过两接收信号内积可以得到归一化正交相位差信号， 即

则相位差序列为

3.2 干涉仪测向误差分析

对式（22）全变量微分， 可得

做变量代换λ=c/f。 假设干涉仪基线长度误差dl=0， 其他误差源相互独立， 由误差传递公式可得

当电子侦察系统对辐射源信号的采样点数为L时， 双阵元干涉仪的联合接收矢量可写为

对于高频窄带信号， 相对测频误差σf/f的值很小， 可忽略不计， 把式（26）代入式（24）， 得到DOA估计的CRLB为

式（27）表明， 在加性高斯白噪声背景下， 干涉仪DOA估计的CRLB与接收信噪比η、 信号采样点数L、 干涉仪基线长度l呈负相关关系， 与信号波长λ和真实到达角α呈正相关关系。

4 仿真实验分析

由理论分析可知， 干涉仪对频分正交MIMO雷达的测向误差源主要包括两部分： （1）相位波前畸变的影响。 由于干涉仪实际上测量的是相位梯度的方向， 频分正交MIMO雷达的相位波前畸变导致相位梯度出现偏差角， 从而增大干涉仪测向误差。 （2）瞬时信噪比的影响。 频分正交MIMO雷达合成信号的时变包络特性会造成干涉仪接收信号瞬时信噪比的起伏， 从而增大测向误差。

注意到控制参数θa可以调控包络峰值位置。 考虑到电子侦察系统在截获辐射源信号后， 为了满足实时测量要求， 通常会截取一段时长的信号进行分析， 如现代宽带侦察接收机一般截取脉冲前沿的一小段时间TΔ进行参数估计。 因此， 通过调整脉冲包络峰值的位置， 可调节并改变侦察接收机分析时段内合成信号的信噪比， 从而增大干涉仪的测向误差。

通过仿真实验分析频分正交MIMO雷达对相位干涉仪DOA测量误差的影响。 仿真实验以地面防空雷达对抗双阵元干涉仪测向的MALD为例。 考虑MALD头锥部尺寸在分米级， L波段MIMO雷达的阵列孔径通常在米级， 因此， 仿真实验MIMO雷达和MALD测向接收机的典型参数设置如表1和表2所示。