李相平 王明泽* 但 波 李 蔚 马俊伟

①(海军航空大学岸防兵学院 烟台 264001)

②(海军航空大学航空作战勤务学院 烟台 264001)

③(海军航空大学航空保障专业兵训练基地 青岛 266109)

1 引言

区别于自由空间中的其他雷达，穿墙成像雷达(Through-the Wall Imaging Radar, TWIR)需要对墙后目标进行探测成像[1,2]。在这一过程中，由墙体反射造成的杂波信号会以较大的幅值来“遮蔽”目标信号，或者在时域上与目标信号交织重叠。因此，有效抑制墙体杂波[3]，是TWIR对墙后目标准确成像的重要前提。

诸如奇异值分解[4](Singular Value Decomposition, SVD)一类的经典杂波抑制算法，仅能实现对主墙体杂波的消除，而残余杂波会显著降低目标成像质量。随着对TWIR的广泛应用和深入研究，目标检测、识别等工作对前期成像的要求越发严格，现有的杂波抑制算法难以达到穿墙成像的实时性和准确性条件。

近年来，机器学习理论被逐渐引入TWIR领域，如压缩感知[5,6](Compressive Sensing, CS)、矩阵补全[7](Matrix Completion, MC)以及鲁棒主成分分析[8](Robust Principal Component Analysis, RPCA)。其中，作为高光谱图像去噪及视频监控中的前景提取等研究方向的常用工具，RPCA[9]可以将数据矩阵唯一分解为低秩杂波矩阵与稀疏目标矩阵，即同时实现对杂波与目标的准确分离。对应于数据矩阵为穿墙回波矩阵或穿墙图像矩阵，RPCA可以分离得到其中的目标回波信号或目标像素分量，进而实现杂波抑制。所以，在本文中，基于RPCA理论，提出一种新的杂波抑制算法。

为满足穿墙成像的实时性和准确性要求，分别就算法的速度和精度进行设计改进：一是提出光滑化快速交替线性化方法来缩短RPCA问题的求解时间，从而加快算法的速度；二是多域联合处理，同时在回波域和图像域应用RPCA理论，通过对其得到的目标图像进行指数加权联乘融合处理来提高算法精度。

综上所述，本文的章节安排如下：第2节介绍RPCA原理；对回波域和图像域建模的工作在第3节完成；第4节对所提杂波抑制算法进行理论分析；在第5节设计穿墙场景仿真，通过与多种算法的性能对比，验证所提算法的速度和精度；最后在第6节总结全文。

2 RPCA原理

其中，//A//∗表示矩阵A的核范数，//E//1表示矩阵E的l1范数。在这里，矩阵A和E要达到秩稀疏不相干的条件，以确保凸优化问题的正确求解。显然，在相当广泛的情况下，该条件都是成立的，即矩阵A和E可以被高概率地精确恢复[13]。

对于式(3)中RPCA问题的求解，有以下几种常用的1阶方法：一是加速近端梯度[14](Accelerated Proximal Gradient, APG)方法，以线性化函数来实现对目标函数的局部逼近，使得单次迭代的计算量较小，但收敛速度较慢；二是精确增广拉格朗日乘子[15](Exact Augmented Lagrange Multipliers,EALM)方法，通过交替迭代增广拉格朗日函数来得到矩阵A和E，可实现高阶线性收敛，不过每次迭代的SVD次数较多，大大影响了计算速度；三是非精确增广拉格朗日乘子[15](Inexact Augmented Lagrange Multipliers, IALM)方法，将EALM中的多轮交替最小化削减为1轮，可在保持原有收敛速度的基础上有效减少SVD次数。

3 联合低秩稀疏模型

3.1 回波域建模

首先，设定成像场景如下：TWIR以收发同置的均匀天线阵列对墙后目标进行探测成像，阵元数量为N，且天线阵列与匀质墙体表面平行，则第n个天线阵元所接收的回波信号为

其中，uw(n,t)为墙体杂波信号，utg(n,t)为目标回波信号，uno(n,t)为噪声信号，uant(n,t)为天线耦合信号。在本文中，认为天线耦合波经预处理后被消除；且为便于建模分析，假设噪声信号对目标成像的影响相对有限，则回波信号可进一步简化为

随后，将式(5)的结果扩展至整个天线阵列，有2维回波矩阵

其中，T为各天线通道信号的采样点数，Uw为杂波信号矩阵，Utg为目标信号矩阵。

在上述成像场景下，单就墙体杂波信号来讲，各天线阵元对应的墙体反射系数及信号传输路径是相同的，则其各自接收的墙体反射回波也是相同的[16–18]，即Uw是低秩矩阵。从空间上看，墙后目标相对于整个探测区域是稀疏的；类似地，在回波域中，目标采样点相对于所有采样点也是稀疏的，即Utg为稀疏矩阵。那么，根据式(3)，在回波域中建立联合低秩稀疏[19](Joint Low-Rank and Sparse,JLRS)模型，有

3.2 图像域建模

相较于回波域，图像域是对探测区域的直观表现，其应用RPCA理论的效果理应更好[20]。对于式(6)中的2维回波矩阵，以后向投影 (Back Projection, BP)算法进行成像处理，则图像中任一像素的幅度为

4 理论分析

4.1 光滑化快速交替线性化方法

根据上文对RPCA问题求解方法的分析，可以得到：线性化近似会减小单次迭代的计算量，而交替迭代会加速收敛，降低迭代复杂度。由此出发，为了加快RPCA问题求解速度，在这里提出一种光滑化快速交替线性化(Smoothing Fast Alternating Linearization, SFAL)方法。

4.2 指数加权联乘多域图像融合处理

根据之前所建立的JLRS模型，不难发现，杂波和目标在各域中并非总是保持严格低秩或严格稀疏的。在回波域中，杂波的低秩性较好，目标的稀疏性较弱，在图像域中则恰恰相反。二者在性质上呈现出一种“互补关系”，而这种关系在各域得到的目标图像中体现得尤为明显。为了提高算法精度、改善目标成像质量，对目标图像进行指数加权联乘多域图像融合处理。

首先，对回波域和图像域中的目标图像作如下分析：在回波域中，JLRS模型中的目标信号矩阵为实矩阵(对原始回波数据取模求得)，隐藏了目标的相位信息，由此得到的目标图像相当于非相干BP成像结果，再加上残余杂波的影响，其方位向分辨率和聚焦效果都有所降低，但墙体主杂波抑制相对彻底；在图像域中，对原始回波成像直接进行分解，此时的目标图像为相干BP成像结果，其目标聚焦效果相对较好，但杂波分离不够彻底，且对目标形成一定的“遮蔽”效应。

经上述分析可得，相较于多通道或多角度的子图像融合，回波域和图像域中的目标图像存在一定的差异性，基本的联乘融合[26]对目标成像质量的改善很是有限。由此提出如下的指数加权联乘融合思路，即

随着加权指数a的增大，回波域目标图像在融合图像中的贡献相应增加，杂波被逐渐消除；在杂波基本消失后，继续增大加权指数，便会造成目标像素的损失，因此找到杂波消失的节点显得尤为重要。然而，在这一过程中，像素均值持续下降，且下降速度并没有出现明显的变化点。

表1 SFAL方法

不过，受此启发，可考虑超均值像素数这一指标。相较于像素均值，超均值像素数能更加全面地描述图像中像素变化的整体情况和局部细节，尤其是杂波像素和目标像素的幅度变化过程。设图像中所有像素点的集合为S，则该指标的定义为

其中，幅度超过像素均值的像素点的集合B ⊆S,crad(B)表示集合B的元素个数。

同样是在加权指数增大的过程中，超均值像素数呈现出有规律的变化：从一开始，随着幅度较大的杂波像素被消除，超均值像素数以较快的速度持续减少；在杂波基本消失后，目标像素成为图像中待消除像素的“主体”，指标的减少速度会明显放缓。所以，以超均值像素数作为评价指标，并以其减速放缓的第一个点作为加权指数a。

4.3 算法思路整合

结合上文内容，对所提杂波抑制算法的思路作进一步整合，其整体流程如图1所示。

图1 杂波抑制算法流程图

同时，对算法的主要步骤可归纳如下：

第1阶段：利用SFAL方法得到各域中的目标图像，包括步骤1和步骤2。

步骤1 根据穿墙雷达回波信号，直接建立回波域JLRS模型；同时，以BP算法处理回波信号，得到原始图像，并以此建立图像域JLRS模型；

步骤2 利用SFAL方法依次求解各域中的JLRS模型，分别得到各域中的目标图像；

第2阶段：对各域中的目标图像进行指数加权联乘融合处理，包括步骤3和步骤4。

步骤3 给定b=1，并在a ∈[0,1]的范围内以0.05的步长绘制超均值像素数的变化曲线，以曲线减速明显放缓的第一个点作为加权指数a；

步骤4 根据选定的加权指数a和b进行指数加权联乘多域图像融合处理，得到最终的融合图像。

5 仿真验证

为了验证本文所提杂波抑制算法的速度与精度，利用MATLAB构建如下的穿墙场景：选择匀质混凝土墙体，其厚度为0.2m，相对介电常数为5.0，磁导率为4π×10−7H/m，电导率为5.0×10−2S/m。天线阵列为均匀线阵，阵列孔径为2 m，阵元间距为0.1 m，阵元数量N=21，其距离墙体的垂直距离为1.0 m。天线工作在收发同置模式，各阵元依次发射信号并由自身接收回波，且各天线通道信号的采样点数T=1024。发射信号为单位幅度的高斯脉冲2阶导信号，脉冲形成因子为0.5 ns，脉冲持续时间约为1.1 ns，对应的–3 dB频谱范围为0.8～2.6 GHz。探测区域为墙后2×2.8 m的范围，共放置4个半径为3.5 cm的金属小球，各金属球完全相同，在此场景的分辨率下可视作点目标处理。场景如图2所示。

图2 穿墙场景示意图

5.1 算法速度

在本文中，对算法速度的改进主要体现在对RPCA问题求解方法的优化。那么，不妨以图像域JLRS模型为例，验证SFAL方法的性能。

首先，在上述穿墙场景下进行目标探测，以BP算法对回波信号处理成像，成像网格设置为256×256，并建立JLRS模型。随后，依次以SFAL,APG, EALM, IALM方法求解该模型，得到目标图像。在参数设置方面，各方法采用默认参数，其迭代终止阈值统一设定为ζ=10−7。在运行环境方面，所用联想台式机配置Intel Core i7-9700 3.00 GHz的8核中央处理器，其安装内存为8 GB，仿真软件为Windows 7系统下的MATLAB R2018b软件。最后，以迭代时间和目标杂波比等指标对各方法的速度及准确度进行对比，具体情况如表2所示。为了不失一般性，表中各数据均为多次独立重复试验结果的平均值。

表2 各方法性能对比

在准确度方面，各方法所得目标图像的目标杂波比基本持平，说明其目标成像质量相近；在速度方面，SFAL方法的迭代次数最少，迭代时间最短，相较于迭代次数较少的EALM方法及迭代时间较短的IALM方法仍有较大优势。总的来看，在准确度相当的前提下，SFAL方法以其速度优势展现出了良好的性能，而这也进一步加快了杂波抑制算法的整体速度。

5.2 算法精度

为了验证算法的精度，同样以图2的穿墙场景进行探测成像，其结果如图3所示。可以看到，其中只有墙体位置附近成像清晰，隐约可见墙后目标阴影。以本文所提算法进行杂波抑制，建立求解回波域和图像域JLRS模型，得到各自的目标图像，其结果分别如图4、图5所示。前者的墙体主杂波抑制较为彻底，但目标散焦严重，难以准确定位，且部分目标丢失；后者的目标聚焦准确，但不够清晰，且墙体位置附近有残余杂波。显然，单一域中的目标图像均未达到对目标成像质量的要求。

图3 原始回波成像

图4 回波域目标图像

图5 图像域目标图像

在进行多域联合处理之前，对加权指数a绘制超均值像素数的变化曲线。由图6可知，整体上该曲线呈下降趋势，初期的上升峰值可解释为杂波像素的快速消除和像素均值的急速下降所导致的超均值像素数的短暂增加，另有回波域加权指数a=0.65。同时，给定图像域加权指数b=1。

图6 超均值像素数变化曲线

根据上面得到的加权指数，进行指数加权联乘多域图像融合处理，最终成像结果如图7所示。相较于图8的背景对消成像和图9的SVD算法成像，多域联合成像中对墙体杂波尤其是残余杂波的抑制更为彻底，目标成像较为清晰，聚焦效果更好。

图7 多域联合成像

图9 SVD算法成像

此外，以目标杂波比来定量分析算法的精度，相应数据如表3所示。可以看到，与SVD算法和背景对消相比，本文所提算法对原始成像中目标杂波比的提升异常明显。总的来看，本文所提算法精度较高，能够较为彻底地抑制杂波，有效提升目标杂波比。

表3 各情况下的目标杂波比(dB)

6 结束语

本文针对穿墙雷达杂波抑制问题提出一种基于RPCA理论的多域联合抑制算法，该算法通过多域图像融合使目标成像准确聚焦，并利用SFAL方法来增强这一过程的实时性。经仿真证明，该算法具有良好的速度和精度，可实现对杂波的充分抑制，能够快速地为目标检测、识别等后续处理提供准确的目标信息。

实际上，本文算法仍限于BP成像和收发同置的工作模式。为了更广泛地适用于穿墙雷达领域，可考虑将所提算法推广至更一般的穿墙场景，如工作在收发分置模式下的天线阵列。而这意味着要解决一个关键问题，即杂波低秩性和目标稀疏性的增强问题。JLRS模型越标准，杂波和目标的分离就越彻底，目标成像质量就越高。同时，从降低单次迭代复杂度的角度入手来优化RPCA求解方法，可使算法的运算速度迈向新的台阶。