公富康 张顺生
(电子科技大学电子科学技术研究院, 四川成都 611731)
利用多个发射机和多个接收机同时进行目标观测的分布式无源雷达成像系统,相比于传统收发一体的单站雷达能够获取更多的目标信息,具有更好的成像性能,因而成为雷达成像的一个热门研究领域。在分布式无源雷达成像系统中,由于接收到的是空间频率非均匀分布的空间信号,致使基于传统的SAR/ISAR算法并不能取得良好的成像效果[1]。研究表明,可以利用目标回波空间频谱和目标散射点之间的傅里叶变换关系实现分布式无源雷达成像[2]。
在实际的成像系统中,成本的制约使得难以构造数量众多的发射机和接收机来满足目标网格空间划分,考虑到接收信号空间谱和目标散射点的稀疏特性,基于压缩感知的成像技术会具有更好的性能。文献[3]中指出接收信号的空间频率对成像结果有重要的影响,它们与信号频谱,发射机和接收机的位置直接相关,同时将互相关系数作为优化目标函数对固定发射机位置情况下的接收机布局进行优化,因此在本文第2部分讨论了采用成像分辨率作为优化目标函数,基于遗传算法配置接收机的位置以实现最佳成像问题。
通过对分布式无源雷达系统接收机的布局优化能够在一定程度上提高其稀疏成像性能,但是在信噪比较低的情况下,基于OMP的稀疏重构算法需要引入大量其他信号分量来降低对前面信号的错误估计造成的影响[4],而稀疏贝叶斯学习的方法假设信号的每一个分量独立和具有不同的先验概率分布,进而获得信号的稀疏表示,如果将其应用于信号重构,能有效改善低信噪比情况下的成像性能[5]。因此在本文第3部分讨论了一种基于 SBL和协方差稀疏表示的分布式无源雷达成像算法来提高成像精度。在接下来的第4部分通过仿真实验对前面讨论的算法进行验证,最后总结全文。
图1表示一个包含m个发射机和n个接收机的分布式无源雷达成像系统,假设发射机m发射的信号为:
图1 分布式雷达成像构型示意图Fig.1 Schematic diagram of distributed radar imaging configuration
(1)
式中um(t)和fm分别为信号的包络和载频,则在成像空间Ω中接收机n接收到经过目标反射的发射信号为:
τmn(r)]exp(j2πfm(t-τmn(r))]dr+nn(t)
(2)
式中A(r)为传播路径衰减,α(r)为目标散射系数,nn(t)为接收噪声,τmn(r)为传输时延。
对式(1)(2)做相关处理并在滤波后信号的峰值进行采样后得到:
(3)
其中kmnx、kmny和kmnz为空间谱域离散点所对应的点,将成像空间划分为G=KLQ离散空间网格,K、L和Q是X、Y和Z方向的网格点数,推导出矩阵表示的形式可得:
(4)
(5)
根据压缩感知理论[6],如果矢量σ是稀疏矢量,则可以使用OMP算法来求解。
分布式雷达的空间谱与其几何构型有关,当按照距离向和方位向对空间谱范围进行划分时,可推导出当前构型的理想成像分辨率为:
(6)
x和y分别是x和y方向的极限分辨率,kxmax和kymax分别表示x和y方向上空间采样点分布的最大范围。由于ρx和ρy为θ和φ的非线性函数,则可以将成像分辨率最高作为优化目标函数,通过最优化算法就可以求出使得距离向分辨率和方位向分辨率最高的接收机分布。
图2 基于遗传算法的布局优化流程图Fig.2 Flow chart of layout optimization based on genetic algorithm
在对空中静止或运动目标一次快拍下成像时,基于压缩感知的稀疏成像算法是以目标信号的线性组合与噪声之和的数据模型为基础,因此,不同的模型会由于参数估计的误差得到不同的结果[8],针对这个问题,提出一种以观测数据协方差稀疏表示为基础的分布式无源雷达成像算法。
当在稀疏度为K的空间成像时,由于发射信号受到高斯白噪声的影响,则可以表示为:
(7)
(8)
(9)
根据贝叶斯理论,观测信号的高斯似然函数为[10]:
(10)
其中σ2为等效噪声转换为白噪声后的协方差,我们通过加入目标成像空间稀疏化先验约束来避免数据的过学习,假定参数σ2和γ是已知的,则其稀疏系数ω的高斯后验概论分布为:
P(ω|y,γ,σ-2)=
(11)
图3 基于协方差稀疏表示的SBL算法Fig.3 SBL algorithm based on covariance sparse representation
图4 遗传算法优化前后接收机配置示意图Fig.4 Schematic diagram of receiver configuration before and after genetic algorithm optimization
为了验证基于协方差稀疏表示的SBL成像算法的有效性,首先对比该算法对于静止点目标二维场景的重构能力,仿真参数如表1所示,由图5所知,协方差域成像恢复结果与仿真场景一致,原始数据域已经无法成像。
表1 分布式无源雷达成像仿真参数
图5 基于协方差稀疏表示的SBL算法成像仿真结果对比Fig.5 Comparison of imaging simulation results of SBL algorithm based on covariance sparse representation
其次对比基于协方差稀疏表示的SBL算法和基于数据域表示的算法在恢复性能上的差异,设计不同信噪比下协方差域与数据域SBL算法归一化重构信号误差RMSE蒙特卡洛仿真如图6所示。
图6 数据域与协方差域重构误差对比Fig.6 Comparison of data domain and based on covariance domain reconstruction error
图6所示,当信噪比大于5dB时,二者相差不大,都能实现重构;当信噪比在-10dB~ 0dB的时候,数据域性能极具下降,但是协方差域重构误差变化不大,仍然能够得到较好的图像。仿真结果说明,在低信噪比的情况下,基于协方差域的成像算法重构性能要好于基于数据域的成像算法。
本文针对分布式无源雷达稀疏成像问题,提出以成像分辨率为优化目标函数,基于遗传算法的接收机配置优化方法,并通过仿真结果表明优化后对成像性能有明显的提升。同时针对OMP算法在低信噪比情况下成像精度较差的问题,提出基于协方差稀疏表示的分布式无源雷达成像算法,并将稀疏贝叶斯学习(SBL)方法应用于目标场景重构,并通过仿真验证在低分辨率下算法的良好恢复性能。