
Citation: | ZHANG Shuangxi, QIAO Ning, XING Mengdao, et al. A novel clutter suppression approach for the space-borne multiple channel in the azimuth high-resolution and wide-swath SAR-GMTI system with an ambiguous Doppler spectrum[J]. Journal of Radars, 2020, 9(2): 295–303. doi: 10.12000/JR20005 |
作为信息的载体,电磁波可由多个维度调制信息,如幅度、相位、频率和极化等[1]。根据经典电动力学理论,Maxwell方程所描述的电磁波包含线动量和角动量[2],而角动量又包括自旋角动量(Spin Angular Momentum, SAM)和轨道角动量(Orbital Angular Momentum, OAM)[3]。理论上,携带轨道角动量的电磁波束具有无穷多种模态且各模态间相互正交,因此可实现模分复用从而提高频谱效率及信道容量以解决探测系统探测效率低的问题[4,5]。各模态的OAM波束的相位呈差异性分布,基于OAM波束的探测成像系统可从回波信号中获取更多的目标信息以实现高分辨率成像,因此,OAM波束在雷达探测领域也有着巨大的应用前景[6-10]。
针对OAM在雷达探测中的应用,相关学者进行了大量研究,主要包括基于OAM的目标重构及OAM波束的传播特性两个方面。针对基于OAM的目标重构,2013年,郭桂蓉团队[1]构建了基于OAM波束的成像模型,利用傅里叶变换的方法对回波信号进行处理,验证了OAM波束在雷达目标成像领域存在巨大的应用潜力。2017年,袁铁柱[6]分析了基于傅里叶变换电磁涡旋波成像方法,验证了目标处于小俯仰角条件下成像的可行性。2019年,Wang等人[7]提出一种涡流SAR(Synthetic Aperture Radar)成像方法,将OAM波束和条带SAR技术相结合,通过改进的距离多普勒(Range-Doppler, RD)算法来获取目标的二维图像。由仿真结果可知,该涡旋SAR成像方法的方位分辨率比常规SAR提高近8倍。2020年,Liu等人[8] 研究了双耦合OAM波束的电磁涡流成像,采用傅里叶变换的方法对目标进行重建,有效缩短了传统成像所需要的时间,大大提高了雷达成像效率。2020年,杜永兴等人[9]建立了多收多发电磁涡旋SAR成像模型,极大提高了方位向分辨率。2021年,Liu等人[10]提出一种基于分数阶OAM波束的电磁涡流增强成像方法,有效缓解了噪声对成像性能的影响。但上述研究均未考虑不同模态OAM波束的主瓣指向差异对成像系统性能的影响,本文对同主瓣指向的多模态OAM波束在目标特征识别的应用进行理论分析和仿真,进一步挖掘OAM波束在雷达探测中的潜力。针对OAM波束的传播特性,2017年,Yu等人[11]研究了远场目标对高阶贝塞尔涡旋波束的散射。2019年,相关学者将OAM波束进行分解,通过理论分析及数值仿真,归纳了OAM波束在两介质交界面及绝缘平板上的反射及折射规律[12]。2020年,Liu等人[13]研究了OAM波束与电大尺寸标准物体的相互作用,分析了金属球体和锥体的后向散射特性。为进一步分析多模OAM波束在目标探测及识别中的应用潜能,基于上述研究,本文设置多种观测场景及反射场景,通过理论分析及数值仿真对各场景下OAM波束呈现的传播特性及反射特性进行归纳总结并对同主瓣指向的多模态OAM波束在目标探测及识别方面的应用进行研究。
目标识别的主要流程包括:(1)从目标散射回波中提取目标特征;(2)由(1)获得的目标特征建立特征数据库;(3)获取待识别目标的特征数据;(4)将待识别目标的特征数据与数据库对比并输出判决。可知目标识别的关键是目标特征的提取,而卷积神经网络(Convolutional Neural Network, CNN)作为一种经典的深度学习算法是深度学习技术中特征提取的有效方法之一[14]。于1988年提出的LeNet-5被视为第1个真正意义上的卷积神经网络,因此也成为推动卷积神经网络发展的基石[15]。2012年,AlexNet[16]在ImageNet大规模视觉识别挑战赛上成功夺冠使得卷积神经网络再次引起相关学者的注意。随着对卷积神经网络研究的不断深入,VGGNet, GoogleNet, ResNet等网络被相继提出[17-19]。本文通过卷积神经网络建立目标反射信息与目标特征识别的关联,从而实现目标特征的智能识别。鉴于本实验的输入为数值矩阵且数据量较小,采用上述网络反而易引起过拟合、计算复杂度提高及梯度弥散等问题,因此本文基于LeNet-5网络搭建了一个3层的网络结构(其中卷积层和池化层作为一个layer)进行实验。
本文各章节内容介绍如下:第2节通过理论分析详细介绍了OAM波束特征的影响因素,并通过数值仿真获得在不同场景下OAM波束的传播及反射特性。第3节搭建卷积神经网络对第2节中各反射场景获得的回波信号进行训练,实现对未知场景的判断及场景内目标位置的识别。第4节对全文内容进行总结。
OAM波束的产生方式包括螺旋相位板、旋转抛物面、均匀圆形阵列及超表面等[20,21]。其他产生方式通过圆形阵列等效可以较容易地实现对OAM波束的幅相特性分析,本文以圆形阵列为基础,采用不同半径的多环结构产生发散角相同的多模态OAM,假设其阵元数目为
El(r,θ,φ)=−jωμN∑n=1∫Jnexp(−jk|r−rn|)4π|r−rn|dVn=−ωμWdNexp(jkr)j−l+14πr⋅exp(jlφ)Jl(kRsinθ) | (1) |
由式(1)不难发现,区别于传统电磁波束,OAM波束的幅度特征与贝塞尔函数保持一致,沿传播轴的电场幅度为0,且其相位与方位角呈线性关系,绕传播轴一周其相位变化
θ0=argmax{J2l(kRsinθ)} | (2) |
对于图1所示的实验场景,获得2, 3模态OAM波束的幅度和相位分布如图2所示,随模态的增大波束发散角增大。由式(2)可知通过调整圆环的半径可以改变波束的发散角,在发散角为
模态l | 阵列半径R (m) |
1 | 0.0450 |
2 | 0.0747 |
3 | 0.1026 |
4 | 0.1302 |
5 | 0.1569 |
6 | 0.1836 |
7 | 0.2097 |
8 | 0.2361 |
9 | 0.2619 |
本文以3种观测面为例,分析OAM波束在传播过程中的幅度及相位分布,如图3所示观测面1位于发射阵列的左上方且垂直于OAM波束的主瓣指向,观测面2位于发射阵列的正上方且垂直于
由图4可知,与观测面2的电场分布相比,观测面1与观测面3均无法展现OAM波束的幅相特性。同时对比图4(b)、图4(d)、图4(f)可以发现3种观测面的相位分布存在一个共性,如图5所示,沿红色箭头所指方向,相位变化速度越来越快。
为进一步验证图5所示规律,提取式(1)中的相位项
p(y)=−(π2+kr)=−(π2+k√y2+z2) | (3) |
其中,
|p′(y)|=ky√y2+z2 | (4) |
|p′(y)|′=kz2(y2+z2)3/2 | (5) |
以上分析了OAM波束传播到不同观测面的幅度和相位分布。下面将分析OAM波束经上述3种反射面获得的反射回波的幅度特征,为观察OAM波束与传统电磁波束反射特性的差异,本文以
−ωμWdNexp(jkrmn)j−l+14πrmnexp(jlφmn)Jl(kRsinθmn) |
假设入射波经反射面被全反射,结合惠更斯原理[22]则空间中某点
−ωμWdNexp(jkrmn)j−l+14πrmnexp(jlφmn)Jl(kRsinθmn) |
基于阵列合成的方法,可获得
Er,l(r,θ,φ)=−jωμM∑m=1N∑n=1∫−ωμWdNexp(jkrmn)j−l+14πrmn⋅exp(jlφmn)Jl(kRsinθmn)⋅exp(−jk|r−rn|)4π|r−rn|dVmn | (6) |
由于贝塞尔函数为特殊函数,难以将式(6)进一步简化为一个闭式表达式,因此只能通过仿真观察不同反射面对应的反射波特征,最终,获得OAM波束经3种反射面的反射回波如图7所示。由图7(a)和图7(b)不难看出,传统电磁波与2模态OAM波束经反射面2获得的反射幅度方向图有明显差异;而其余2种反射面对应的反射幅度方向图难以看出幅度差异。传统电磁波和OAM波束经反射面1和反射面3对应的反射波束差异甚微主要因为反射平面尺寸较小且距
l | θ (°) | φ (°) |
1 | 11.2523 | 89.9940 |
2 | 11.2523 | 89.9819 |
3 | 11.2523 | 89.9819 |
4 | 11.2523 | 89.9699 |
5 | 11.2523 | 89.9579 |
6 | 11.2523 | 89.9579 |
7 | 11.2523 | 89.9458 |
8 | 11.2523 | 89.9458 |
由表2不难看出,模态不同,主瓣指向相同的OAM波束经相同的反射平面获得的俯仰向反射角
本节设置3种反射场景,不同场景中的反射目标呈某一趋势动态变化,通过分析对应场景的反射回波讨论多模态OAM波束用于未知场景的判断及场景内目标位置识别的可行性,3种场景如图8所示。
场景1内的反射面初始位置与
为便于分析,本文将采集的回波幅度信息进行归一化。由图9可知,在反射场景1中,反射面处于不同位置时观测点处观测到不同模态OAM波束对应的反射波幅度变化规律基本保持一致;各模态OAM波束对应的反射波的相位变化趋势相同但奇数模态与偶数模态OAM波束对应的反射波束存在一定的相位差;对于场景2,随反射平面的位置变化,不同模态OAM波束对应的反射波束的幅度变化趋势存在明显差异且相位变化也存在差异,如图10所示;对于场景3,虽然不同模态OAM波束对应的回波信号的幅度保持一致,但相位具有明显差异。综合可知,3种场景间的回波信号特征各异,且各场景内不同位置反射目标的回波也存着幅度或相位的差异,因此对于多模OAM波束的探测系统可通过场景间的信息差异判断反射场景类别,通过提取场景内不同位置的差异信息识别场景内目标。
第2节介绍了基于圆形阵列的OAM波束产生方法,并通过仿真实验对OAM波束的传播特性及不同场景下的反射特性进行了深入研究。本节搭建卷积神经网络,将2.3节中的3种场景反射数据进行训练从而实现未知场景的判断,随后对3种反射场景的反射目标位置进行识别。
实验开始前先对原始数据进行预处理,将其转化为数值矩阵形式,而灰度图像也是一个数值矩阵,因此本网络的输入可视为单通道的灰度图像,且灰度图像的尺寸较小。因此,本文所研究的问题仅需一个结构简单的网络即可解决。本文所用网络由2个卷积层、2个池化层和1个全连接层构成,其中输入数据为28×28的实数矩阵,卷积核的尺寸为5×5,数目为6,池化窗口的尺寸为2×2,如图12所示。
(1) 参数初始化
训练开始前,首先采用随机数初始化各个参数权值。为避免网络因权值过大而进入饱和状态,一般选取较小的随机数作为初始值赋给各参数,且各权值的初始值应各不相同,否则网络丧失学习能力。本文采用的初始化方法为Xavier初始化,其思想是使前向传播和反向传播的输出与输入具有相似的分布,以便其方差保持一致,一定程度上可避免网络训练时的梯度发散问题[23]。Xavier初始化方法的表达式为
w=rand(−1,1)×√6m+n | (7) |
其中,
(2) 前向传播
将训练数据进入网络经卷积层conv1、池化层pooling1、卷积层conv2、池化层pooling2及全连接层获得输出向量。
(3) 反向传播
反向传播算法是基于梯度下降算法建立的,根据目标的负梯度方向调整参数,本文采用均方误差函数作为损失函数对结果进行优化,如式(8)所示。首先设定学习率,并计算出前向传播输出结果与目标的误差;然后,分别计算出输出层和前面隐层神经元的梯度项;最后,根据每一层的误差梯度更新各层权值及偏置。经多次迭代,损失函数逐渐减小从而提升网络训练效果。
Loss=N∑i=1(ei)22N | (8) |
其中,
2.3节分析了不同模态OAM波束在3种反射场景的反射特征,本节对3种反射场景中指定观测点处获得的反射数据进行训练从而实现未知场景的判断。已知2.3节采集的反射数据为复数形式,而所搭建网络的输入为
实验参数 | 参数值 |
频率 (GHz) | 10 |
阵元数目 (个) | 40 |
反射面尺寸 | 10λ |
迭代次数 | 100 |
反向传播学习率 | 1 |
各场景选取的训练数据数目和测试数据数目分别为300和40。对于场景1,反射平面的旋转范围
P=(1−mN)×100% | (9) |
其中,m表示错误识别数目,N表示测试数据集的总数目。
由实验结果可知,所搭建网络对于3种反射场景的场景判断准确率为
基于3.1节的实验场景,本节将对卷积神经网络用于场景判断的抗噪性能进行分析,以引入加性高斯白噪声为例,首先观察噪声引入后对OAM波束的影响,图14为信噪比为10 dB时2模态OAM波束的幅相特征,对比图4(a)和图4(b)可以看出,引入噪声后OAM波束的幅度及相位均发生一定畸变。
为分析网络抗噪性能,在上节仿真场景的基础上引入噪声。网络的训练数据为未被噪声污染的反射数据。测试数据采用120组引入噪声的反射数据,获得网络对场景的判断准确率与信噪比的关系如图15所示。其中,3种场景的准确判断个数随信噪比的关系如图16所示。
由图15可知,随着信噪比的增大所搭建的卷积神经网络的场景判断准确率增高。对120组测试数据的判断结果进行统计,如图15所示,随信噪比的改变,场景2和场景3的准确判断数目相对稳定,而场景1的判断准确率随信噪比的降低急速下降。造成这种现象的原因可由2.3节中3种反射场景回波的幅相特性进行解释,对于场景1,由图9可知随反射面的变化仅奇数模态OAM波束和偶数模态OAM波束对应反射波存在固定相位差,其余信息均相同,因此对于反射场景1的回波所携带的信息量较少使得其特征不具代表性,当噪声引入时很容易淹没其特征信息从而不易对该场景做出准确判断;反射场景2中,各模态的OAM波束对应的反射波束的幅度及相位均存在差异,因此其携带的特征信息更多,自身特征更加明显从而抗噪声干扰的能力更强,因此低信噪比条件下也能被准确判断;同理,对于场景3,虽然各模态OAM波束对应的反射波束的幅度变化规律相同,但其相位分布存在明显差异,因此,相比反射场景1,反射场景3获得的反射数据携带的信息更多,受噪声的影响更小。
本节将基于上节的卷积网络对某一反射场景中各位置的反射数据进行训练从而对场景间目标位置的识别进行讨论。以第2种应用场景为例对参数设置及数据选取进行说明,对于场景2,反射平面沿
最后,对所搭建网络用于场景内目标识别的抗噪性能进行分析,3种反射场景的可识别间隔分别设置为
为实现高效目标探测,本文将基于多模OAM波束的目标探测方法与卷积神经网络相结合实现对目标特征的智能识别。首先,本文对OAM的传输特性进行了理论分析,给出了满足多模OAM波束特征;在此基础上,对不同模态OAM波束在不同反射场景下的回波信号进行了研究,并采用卷积神经网络提取各应用场景的反射数据特征,从而实现未知场景的判断及场景内目标的识别,同时对整体系统进行了抗噪性能分析。实验结果表明,无噪条件下,目标场景判断的准确率可达
[1] |
GEBERT N. Multi-channel azimuth processing for high-resolution wide-swath SAR imaging[D]. [Ph.D. dissertation], Fakultät für Elektrotechnik und Informationstechnik, 2009.
|
[2] |
邢孟道, 林浩, 陈溅来, 等. 多平台合成孔径雷达成像算法综述[J]. 雷达学报, 2019, 8(6): 732–757. doi: 10.12000/JR19102
XING Mengdao, LIN Hao, CHEN Jianlai, et al. A review of imaging algorithms in multi-platform-borne synthetic aperture radar[J]. Journal of Radars, 2019, 8(6): 732–757. doi: 10.12000/JR19102
|
[3] |
KIM J H, YOUNIS M, PRATS-IRAOLA P, et al. First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(1): 579–590. doi: 10.1109/TGRS.2012.2201947
|
[4] |
LI Zhenfang, WANG Hongyang, SU Tao, et al. Generation of wide-swath and high-resolution SAR images from multichannel small spaceborne SAR systems[J]. IEEE Geoscience and Remote Sensing Letters, 2005, 2(1): 82–86. doi: 10.1109/LGRS.2004.840610
|
[5] |
SIKANETA I, GIERULL C H, and CERUTTI-MAORI D. Optimum signal processing for multichannel SAR: With application to high-resolution wide-swath imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10): 6095–6109. doi: 10.1109/TGRS.2013.2294940
|
[6] |
邓云凯, 赵凤军, 王宇. 星载SAR技术的发展趋势及应用浅析[J]. 雷达学报, 2012, 1(1): 1–10. doi: 10.3724/SP.J.1300.2012.20015
DENG Yunkai, ZHAO Fengjun, and WANG Yu. Brief analysis on the development and application of spaceborne SAR[J]. Journal of Radars, 2012, 1(1): 1–10. doi: 10.3724/SP.J.1300.2012.20015
|
[7] |
SJÖGREN T K, VU V T, PETTERSSON M I, et al. Suppression of clutter in multichannel SAR GMTI[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7): 4005–4013. doi: 10.1109/TGRS.2013.2278701
|
[8] |
ZHANG Shuangxi, XING Mengdao, XIA Xianggen, et al. Robust clutter suppression and moving target imaging approach for multichannel in azimuth high-resolution and wide-swath synthetic aperture radar[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(2): 687–709. doi: 10.1109/TGRS.2014.2327031
|
[9] |
WANG H S C. Mainlobe clutter cancellation by DPCA for space-based radars[C]. The 1991 IEEE Aerospace Applications Conference Digest, Crested Butte, USA, 1991: 1–128. doi: 10.1109/AERO.1991.154520.
|
[10] |
王肖洋, 高贵, 周石琳, 等. 一种基于双通道DPCA的SAR-GMTI杂波抑制方法[J]. 雷达学报, 2014, 3(2): 241–248. doi: 10.3724/SP.J.1300.2014.13121
WANG Xiaoyang, GAO Gui, ZHOU Shilin, et al. A clutter suppression approach for SAR-GMTI based on dual-channel DPCA[J]. Journal of Radars, 2014, 3(2): 241–248. doi: 10.3724/SP.J.1300.2014.13121
|
[11] |
KLEMM R. Introduction to space-time adaptive processing[J]. Electronics & Communication Engineering Journal, 1999, 11(1): 5–12. doi: 10.1049/ecej:19990102
|
[12] |
TENNETI S V and VAIDYANATHAN P P. iMUSIC: A family of MUSIC-like algorithms for integer period estimation[J]. IEEE Transactions on Signal Processing, 2019, 67(2): 367–382. doi: 10.1109/TSP.2018.2879039
|
[13] |
张兴良, 王可人, 樊甫华. 典型阵列快速MUSIC算法研究[J]. 雷达学报, 2012, 1(2): 149–156. doi: 10.3724/SP.J.1300.2012.20026
ZHANG Xingliang, WANG Keren, and FAN Fuhua. Study on fast MUSIC algorithm with typical array[J]. Journal of Radars, 2012, 1(2): 149–156. doi: 10.3724/SP.J.1300.2012.20026
|
[14] |
蒋柏峰, 吕晓德, 向茂生. 基于广义MUSIC 算法的低仰角估计新方法[J]. 雷达学报, 2103, 2(4): 422–429. doi: 10.3724/SP.J.1300.2013.13090
JIANG Baifeng, LÜ Xiaode, and XIANG Maosheng. A new low-elevation estimation method based on a general MUSIC algorithm[J]. Journal of Radars, 2103, 2(4): 422–429. doi: 10.3724/SP.J.1300.2013.13090
|
[15] |
ZHANG Shuangxi, XING Mengdao, XIA Xianggen, et al. A robust imaging algorithm for squint mode multi-channel high-resolution and wide-swath SAR with hybrid baseline and fluctuant terrain[J]. IEEE Journal of Selected Topics in Signal Processing, 2015, 9(8): 1583–1598. doi: 10.1109/jstsp.2015.2464182
|
1. | 潘时龙,章志健,王祥传,刘熙,杨思豪,曹欣荣,王立晗,程强,马丛,丁泽勇,郭德明,伍光新. 微波光子信号同步及其在分布式相参雷达中的应用. 现代雷达. 2024(02): 25-34 . ![]() | |
2. | 张广磊,顾新杨,张道明,莫修辞. 面向系统性能评估的机载微波光子雷达集成验证系统. 电光与控制. 2023(08): 89-95 . ![]() | |
3. | 刘兴华,王国玉,徐振海,汪连栋. 分布式孔径相参合成原理、发展与技术实现综述. 雷达学报. 2023(06): 1229-1248 . ![]() | |
4. | 刘泉华,张凯翔,梁振楠,曾涛,龙腾. 地基分布式相参雷达技术研究综述. 信号处理. 2022(12): 2443-2459 . ![]() | |
5. | 赵锋,徐志明,刘蕾,艾小锋. 弹道目标特征提取研究现状与展望. 信息对抗技术. 2022(03): 15-32 . ![]() | |
6. | 周宝亮. 分布式相参雷达LFM宽带去斜参数估计方法. 电子与信息学报. 2020(07): 1566-1572 . ![]() |
模态l | 阵列半径R (m) |
1 | 0.0450 |
2 | 0.0747 |
3 | 0.1026 |
4 | 0.1302 |
5 | 0.1569 |
6 | 0.1836 |
7 | 0.2097 |
8 | 0.2361 |
9 | 0.2619 |
l | θ (°) | φ (°) |
1 | 11.2523 | 89.9940 |
2 | 11.2523 | 89.9819 |
3 | 11.2523 | 89.9819 |
4 | 11.2523 | 89.9699 |
5 | 11.2523 | 89.9579 |
6 | 11.2523 | 89.9579 |
7 | 11.2523 | 89.9458 |
8 | 11.2523 | 89.9458 |
实验参数 | 参数值 |
频率 (GHz) | 10 |
阵元数目 (个) | 40 |
反射面尺寸 | 10λ |
迭代次数 | 100 |
反向传播学习率 | 1 |