Multi-rank Range-spread Target Detection Method for Space/Time Symmetric Array Radar under Non-Gaussian Clutter Background
-
摘要: 针对多通道阵列雷达从实际杂波中检测目标场景,该文提出了一种面向多通道阵列雷达非高斯杂波背景的多秩距离扩展目标检测方法。首先,利用秩大于1的子空间矩阵和相应距离单元的坐标向量,建立了多秩距离扩展目标模型;然后,利用雷达接收单元空间或时间中心对称探测场景下杂波协方差矩阵的反对称结构信息,通过酉变换,采取广义似然比、Rao、Wald检验准则,构建待解参数小样本估计策略,设计了面向非高斯杂波背景的多秩距离扩展目标检测方法。最后,通过理论推导证明了所提检测方法相对于杂波协方差矩阵具有恒虚警特性。基于仿真数据和实测数据的实验结果表明,所提检测方法能够保证对杂波协方差矩阵具有恒虚警特性,此外,相较于现有检测方法,所提检测方法提升了小训练支持的目标检测性能,并且在导向矢量失配条件下,有效地改善目标检测的稳健性。Abstract: This study proposes a multi-rank range-spread target detection method for multi-channel array radar under a non-Gaussian clutter background. The method aims to detect the target from real clutter using the multi-channel array radar. First, a multi-rank range-spread target model was formulated using a subspace matrix with a rank greater than one and the coordinate vectors of corresponding range bins. Then, by exploiting the persymmetric structure information of the clutter covariance matrix under the detection scenario, wherein the radar receiver units were central symmetric in space or time, a small sample estimation strategy for the parameters to be solved through the unitary transformation was constructed. Further, a non-Gaussian clutter background multi-rank range-spread target detection method was designed based on the generalized likelihood ratio, Rao, and Wald tests. Finally, a theoretical derivation proved that the proposed detection method has the constant false alarm rate property. The experimental results based on both the simulated and measured data showed that the proposed detection method can ensure the constant false alarm rate property of the clutter covariance matrix. Additionally, compared with the existing detection methods, the proposed detection method improves the target detection performance under small sample support. Besides, the proposed detection method effectively improves the robustness of target detection under the condition of steering vector mismatch.
-
Key words:
- Adaptive detection /
- Range-spread target /
- Multi-rank subspace /
- Persymmetry /
- Non-Gaussian
-
1. 引言
空时自适应处理(Space-Time Adaptive Processing, STAP)是机载雷达地/海杂波抑制和运动目标检测的关键技术[1,2]。为设计空时滤波器自适应抑制杂波,STAP方法一般需要利用一定数量的独立同分布(Independent Identically Distributed, IID)训练距离单元估计待测距离单元(Range cell Under Test, RUT)的杂波协方差矩阵(Clutter Covariance Matrix, CCM)。为保证输出信杂噪比相比理想条件的损失不超过3 dB,传统STAP方法所需IID训练单元的数量至少为系统自由度的2倍。然而,在实际非均匀杂波环境中,通常难以获得足够的IID训练单元。为解决这一问题,学者提出降维、降秩、直接数据域、知识辅助和稀疏恢复等STAP新方法[3-8]。其中,稀疏恢复空时自适应处理(SR-STAP)方法基于杂波在角度-多普勒域(即空时二维平面)的稀疏特性,利用少量训练距离单元即可获得杂波空时谱的准确估计,从而重构CCM或杂波子空间,构造空时滤波器对杂波进行抑制[9-16]。
在实际应用中,机载雷达不可避免地存在着阵列误差,包括幅度误差和相位误差。由于误差信息隐含于CCM的估计之中,传统STAP方法具有较强的误差自适应补偿能力。然而,由于SR-STAP方法通常利用理想空时导向矢量构建杂波空时谱估计模型,其性能受误差的影响较大。阵列误差会降低杂波空时谱及CCM的估计准确性,从而严重影响SR-STAP方法的杂波抑制和目标检测性能。针对这一问题,文献[17]提出了基于迭代交替下降(Iterative Alternating Descent, IAD)算法的SR-STAP方法,能够同时估计杂波空时谱和阵列误差参数,但该方法的运算复杂度较高;文献[18]提出了基于ADMM算法的杂波空时谱和阵列误差参数联合估计方法,相比IAD方法运算复杂度较低,但需要同时对多个迭代参数进行设置。
在构建杂波空时谱估计模型或杂波空时谱及阵列误差参数联合估计模型的前提下,现有SR-STAP方法的性能往往依赖于所采用的SR算法。目前,典型的SR算法均基于模型驱动实现,具有理论保证性高、可解释性强等优点。但是,模型驱动类SR算法通常需要设置一个或多个参数,例如正则化因子、迭代步长等。不恰当的参数设置会影响SR算法的收敛速度和精度,从而使得SR-STAP方法的运算复杂度升高、杂波抑制性能下降,限制了其在实际中的应用。针对模型驱动类SR算法存在的问题,受深度学习技术的启发,学者提出了DU方法[19-23]。DU方法将特定SR算法展开为深度神经网络,将算法的迭代次数作为网络的层数、算法的参数作为网络的学习参数,利用训练数据集对SR算法所涉及的迭代参数进行训练,获得最优参数,从而提高SR算法的收敛速度和精度。例如,Gregor等人[19]基于迭代软阈值算法(Iterative Soft Thresholding Algorithm, ISTA),提出了学习型ISTA(Learned ISTA, LISTA)算法;Borgerding等人[21]对近似消息传递(Approximate Message Passing, AMP)算法进行展开,提出了LAMP算法;Yang等人[22]基于近端算子方法(Proximal Operator Methods, POM),提出了LePOM算法。相比其对应的SR算法,DU方法将模型驱动和数据驱动相结合,能够有效降低算法复杂度、提高算法性能。
目前,尚未有研究将DU方法引入到机载雷达SR-STAP之中,且上述DU方法仅能用于杂波空时谱估计,无法同时估计阵列误差参数。因此,为解决现有SR-STAP方法存在的参数设置困难和运算复杂度高等问题,本文提出了机载雷达DU-STAP方法,以验证DU方法在机载雷达杂波抑制和目标检测中的适用性。首先,建立了阵列误差条件下的机载雷达回波信号模型,并利用ADMM算法对杂波空时谱和阵列误差参数进行联合估计;接着,根据对其迭代步骤和数据流图的分析,将ADMM算法展开为深度神经网络,构建具有正则化因子、迭代步长、二次惩罚因子和比例因子等可学习参数的AE-ADMM-Net;然后,定义网络损失函数,基于充足完备的数据集对AE-ADMM-Net进行训练,获得最优参数;最后,利用训练后的AE-ADMM-Net对训练距离单元数据进行处理,快速获得杂波空时谱和阵列误差参数的准确估计,从而设计空时滤波器进行杂波抑制和目标检测。仿真实验表明:相比基于稀疏贝叶斯学习(Sparse Bayesian Learning, SBL)算法[13]、欠定系统聚焦式求解算法(Focal Under-determined System Solver, FOCUSS)[10]和ADMM算法的典型SR-STAP方法,本文所提出的DU-STAP方法均能够在保持较低运算复杂度的同时提高杂波抑制性能。
2. 理论基础
2.1 信号模型
如图1所示,假设机载雷达以高度H、速度v沿y轴匀速飞行,正侧视均匀线阵的阵元个数为M,阵元间距为
d=λ/2 (λ 为波长),脉冲重复频率为fr ,在一个相干处理间隔内共有N个脉冲。不考虑距离模糊杂波的影响,假设每个距离单元所对应的距离环中共有
Nc 个杂波块均匀分布在方位角θ∈[0,π] 上,则包含运动目标的RUT空时回波信号可表示为y0=αTsdsT+Nc∑i=1αisdsi+ε0=yT0+yC0+ε0∈CNM×1 (1) 其中,
αT 和αi 分别表示目标和第i个杂波块的复幅度;sdsT=sdT⊗ssT∈CNM×1 和sdsi=sdi⊗ssi∈CNM×1 分别表示目标和第i个杂波块的空时导向矢量;sdi=[1,ej2πfdi,⋯,ej2π(N−1)fdi]T∈CN×1 和ssi=[1,ej2πfsi,⋯,ej2π(M−1)fsi]T∈CM×1 分别表示第i个杂波块的时域导向矢量和空域导向矢量,fdi=2vcosφicosθi/(λfr) 和fsi=dcosφicosθi/λ 分别表示第i个杂波块的归一化多普勒频率和空间频率,φi 和θi 表示第i个杂波块的俯仰角和方位角;⊗ 表示Kronecker积,[⋅]T 表示转置,j=√−1 ;yT0 ,yC0 和ε0 分别表示目标、杂波和噪声信号。假设各个杂波块之间相互独立,与噪声不相关,且噪声服从均值为0、协方差矩阵为
RN = σ2INM 的复高斯分布,则杂波加噪声协方差矩阵(Clutter plus Noise Covariance Matrix, CNCM)可表示为RC=E[(yC0+ε0)(yC0+ε0)H]=Nc∑i=1E[|αi|2](sdsi)(sdsi)H+RN∈CNM×NM (2) 其中,
E[⋅] 表示期望,[⋅]H 表示共轭转置,INM 表示NM×NM 的单位矩阵。STAP通过计算空时回波信号的加权组合实现对杂波和噪声的抑制以及对运动目标的检测。为使输出信杂噪比(Signal to Clutter plus Noise Ratio, SCNR)最大,空时滤波器的最优权值可通过式(3)计算得出:
wopt=R−1CsdsT/[(sdsT)HR−1CsdsT]∈CNM×1 (3) 其中,
(⋅)−1 表示对矩阵求逆。实际上,RUT的CNCM是未知的,一般需要一定数量的无目标训练距离单元对其进行估计。假设训练距离单元与RUT的杂波独立同分布,则RUT的CNCM可以通过采样协方差矩阵求逆(Sample Matrix Inversion, SMI)方法估计得到[1],表示为
ˆRC=(1/L)L∑l=1ylyHl (4) 其中,
l=1,2,⋯,L ,L表示IID训练距离单元个数,yl 表示第l个训练距离单元的空时回波信号。根据RMB准则[2],SMI方法确保输出SCNR损失小于 3 dB所需的 IID 训练距离单元数应至少为 2 倍的系统自由度。在实际非均匀环境中,该条件通常难以得到满足。此外,实际机载雷达不可避免地存在阵列幅相误差。此时,RUT空时回波信号、CNCM和最优空时权值可分别表示为
˜y0=αT˜sdsT+Nc∑i=1αisdi⊗(e⊙ssi)+ε0=αT˜sdsT+Nc∑i=1αi˜sdsi+ε0=˜yT0+˜yC0+ε0 (5) ˜RC=Nc∑i=1E[|αi|2]sdi⊗(e⊙ssi)[sdi⊗(e⊙ssi)]H+RN=Nc∑i=1E[|αi|2]sdi(sdi)H⊗(eeH⊙ssi(ssi)H)+RN (6) ˜wopt=˜R−1C˜sdsT/[(˜sdsT)H˜R−1C˜sdsT] (7) 其中,
⊙ 表示Hadamard积,˜sdsT=sdT⊗(e⊙ssT) 表示阵列误差下的目标空时导向矢量,e=[e1,e2,⋯,eM]T 表示阵列幅相误差矢量,em=(1+νm)ejϕm ,m=1,2,⋯,M ,νm∈R 和ϕm∈R 分别为第m个实际阵元与理想阵元之间的幅度误差和相位误差。2.2 SR-STAP
由式(1)可以看出,杂波信号可由不同空间频率和多普勒频率的空时信号叠加而成。如果分别将空间频率和多普勒频率离散化为
Ns=κsM 和Nd=κdN 个网格点(其中κs>1 和κd>1 表示尺度因子),则第l个无目标训练距离单元的空时回波信号可表示为yl=NsNd∑q=1αqlsdq⊗ssq+εl=NsNd∑q=1αqlsdsq+εl=Aαl+εl (8) 其中,
αql 为第q个网格点对应的复幅度,q=1,2,⋯,NsNd ,αl=[α1l,α2l,⋯,αNsNdl]T∈CNsNd×1 为所有网格点对应的复幅度矢量,即杂波空时谱;sdq=[1,ej2πfdq,⋯,ej2π(N−1)fdq]T 和ssq=[1,ej2πfsq,⋯,ej2π(M−1)fsq]T 为第q个网格点对应的时域导向矢量和空域导向矢量,fdq 和fsq 为第q个多普勒频率和空间频率,A=[sds1,sds2,⋯,sdsNsNd]∈CNM×NsNd 为空时导向矢量字典,εl 为噪声信号。根据杂波空时谱的稀疏性,可将欠定问题(8)转化为如下约束优化问题进行求解:
argminαl||αl||0s.t.||yl−Aαl||22≤ξ (9) 其中,
||⋅||0 和||⋅||2 分别表示向量的L0范数和L2范数,ξ 表示噪声电平。在存在L个训练距离单元的情况下,式(9)可扩展至多观测模型,表示为
argminΛ||Λ||2,0s.t.||Y−AΛ||2F≤Lξ (10) 其中,
Y=[y1,y2,⋯,yL]∈CNM×L ,Λ=[α1,α2,⋯,αL]∈CNM×L ,||⋅||2,0 表示先对矩阵各行取L2范数再对列取L0范数,||⋅||F 表示矩阵的Frobenius范数。利用L1凸优化算法、FOCUSS算法或SBL算法等稀疏恢复算法对式(9)或式(10)进行求解,可获得
αl 或Λ 的高分辨估计。然后,可通过式(11)计算CNCM,并根据式(3)设计空时滤波器:ˆRC=(1/L)L∑l=1NsNd∑q=1|αql|2sdsq(sdsq)H+RN (11) 同理,当存在阵列误差时,第l个训练距离单元的空时回波信号可表示为
˜yl = Ns×Nd∑q=1αqlsdq⊗(e⊙ssq)+εl=EAαl+εl (12) 其中,
E=IN⊗diag(e) ,IN 表示N×N 的单位矩阵,diag(⋅) 表示取对角矩阵。此时,需要同时估计杂波空时谱
αl 和阵列误差参数e,表示为argminαl,e||αl||0s.t.||˜yl−EAαl||22≤ξ (13) 在求解(13)的基础上,CNCM可通过式(14)进行计算,从而根据式(7)设计空时滤波器:
ˆ˜RC=NsNd∑q=1|αql|2sdq(sdq)H⊗(eeH⊙ssq(ssq)H)+RN (14) SR-STAP方法利用少量甚至单个训练距离单元即可获得CNCM的准确估计,从而实现对杂波的抑制,在实际非均匀环境中具有显著优势。为简便起见,本文仅考虑单个训练距离单元的情况,即L=1,多个训练距离单元的情况可对本文算法进行拓展处理。此外,需要说明的是:在存在距离模糊的情况下,仍然可以建立如式(9)或式(10)所示的优化模型,利用SR算法进行求解,获得距离模糊杂波空时谱的高分辨估计,具体可参考文献[24,25]。
3. DU-STAP方法
为降低运算复杂度、提高杂波抑制性能,本文拟利用DU方法对杂波空时谱和阵列误差参数联合估计模型(13)进行求解。由文献[21]可知,对于
y=Aα + ε 所示的稀疏恢复问题,大多迭代类SR算法的步骤可表示为αk+1=P(αk−γkAH(Aαk−y)) 。其中,αk 为第k次迭代估计结果,γk 为迭代步长,P(⋅) 为非线性算子。令Wk=INM−γkAHA ,Bk=γkAH ,则SR算法的第k次迭代等价于αk+1=P(Wk⋅αk+Bky) 。将Wk 和Bk 定义为深度神经网络的权重参数,P(⋅) 定义为深度神经网络的激活函数,αk 和αk+1 分别定义为深度神经网络第k层的输入和输出,则SR算法的第k次迭代等价于深度神经网络的第k层运算。因此,DU方法可以看作基于SR算法的迭代步骤对深度神经网络的结构和参数进行设计。理论上,LISTA, LAMP和LePOM等DU方法[19-21]均可以实现对杂波空时谱的估计,即对式(9)进行求解。然而,这些方法无法同时估计阵列误差参数,即无法对式(13)进行求解。针对这一问题,本文对ADMM算法[18]进行分析,将其展开为深度神经网络,构建AE-ADMM-Net,实现对杂波空时谱和阵列误差参数的快速准确估计。3.1 ADMM算法
定义
T=IN⊗diag(t) ,其中t=[t1,t2,⋯,tM]T ,tm=e−1m=(1+νm)−1e−jϕm ,则可将式(12)变换为Ty=TEAα+Tε=Aα+ε′ (15) 其中,
y=˜yl ,α=αl ,ε=εl 。由于本文仅考虑单个训练距离单元的情况,因此忽略下标l。式(15)将阵列幅相误差矢量e转化为参数矢量t,可通过式(16)进行求解:
argminα,t||α||1+1/(2ρ)||Ty−Aα||22 (16) 其中,
ρ>0 表示正则化因子。定义辅助变量
η=Ty−Aα ,则式(16)可等效为argminα,η||α||1+1/(2ρ)||η||22 s.t. Aα+η=Ty (17) 等式约束问题(17)的增广拉格朗日函数可表示为
argminα,η,λ,t||α||1+1/(2ρ)||η||22−R{ λH(Aα+η−Ty)} +γ/2||Aα+η−Ty||22 (18) 其中,
λ∈CNM×1 表示拉格朗日乘子,γ>0 表示二次惩罚因子,R{⋅} 表示取实部操作。为避免零解,引入凸约束
∑Mm=1tm=(δ+jw) ,则可将式(18)改写为argminα,η,λ,t||α||1+1/(2ρ)||η||22−R{β∗(M∑m=1tm−δ−jw)}−R{λH(Aα+η−Ty)}+γ/2||Aα+η−Ty||22 (19) 其中,
δ∈R ,w∈R 为比例因子,β 为辅助参数,(⋅)∗ 表示共轭。ADMM算法利用K次迭代交替求解以下4个子问题对式(19)进行求解[18]:
{η(k+1)=argminη 1/(2ρ)||η||22+γ/2||Aα(k)+η−T(k)y−λ(k)/γ||22α(k+1)=argminα ||α||1+γ/2||Aα+η(k+1)−T(k)y−λ(k)/γ||22t(k+1)=argmint γ/2||Aα(k+1)+η(k+1)−Ty−λ(k)/γ||22−R{β∗(M∑m=1tm−δ−jw)}λ(k+1)=argminλ R{λH(Aα(k+1)+η(k+1)−T(k+1)y)} (20) 其中,
T(k+1)=IN⊗diag(t(k+1)) ,α(k+1) ,λ(k+1) ,η(k+1) 和t(k+1) 分别为α ,λ ,η 和t在第k+1次迭代时的估计,k=0,1,⋯,K−1 。式(20)中4个子问题的解可表示为
{X(k+1):η(k+1)=ργ/(1+ργ)(λ(k)/γ−Aα(k)+T(k)y)O(k+1):α(k+1)=soft(α(k)+τAHη(k+1)/(ργ),τ/γ)Z(k+1):t(k+1)=[(b1+β)/a1,⋯,(bM+β)/aM]TM(k+1):λ(k+1)=λ(k)−γ(Aα(k+1)+η(k+1)−T(k+1)y) (21) 其中,
τ 为α 的迭代步长,soft(x,c)=max{|x|−c,0}⋅x/|x| 为软阈值算子[19],am=∑Nn=1|y(n−1)M+m|2,bm=∑Nn=1y∗(n−1)M+mz(k)(n−1)M+m,z(k)=Aα(k+1)+η(k+1)−λ(k)/γ, β=[δ+jw−∑Mm=1(bm/am)]/∑Mm=1(1/am)。 综上所述,利用ADMM算法对式(13)进行求解的步骤如表1所示。需要强调的是:当不存在阵列误差时,表1所示的ADMM算法同样可以对式(9)进行求解。此时,可跳过步骤4,并令
T(k+1)=T(0) 保持不变;也可令比例因子δ=M ,w=0 , ADMM算法将输出阵列误差的估计e≈1M ,即νm≈ϕm≈0 。表 1 ADMM算法Table 1. ADMM algorithm输入:A, y,迭代次数K,正则化因子ρ,二次惩罚因子γ,迭代步长τ,比例因子δ和w。 步骤1 初始化:α(0)=0NdNs(NdNs×1的全0列向量),λ(0)=0NM(NM×1的全0列向量),t(0)=1M(M×1的全1列向量),
T(0)=IN⊗diag(t(0)),k=0;步骤2 η(k+1)=ργ/(1+ργ)(λ(k)/γ−Aα(k)+T(k)y); 步骤3 α(k+1)=soft(α(k)+τAHη(k+1)/(ργ),τ/γ); 步骤4-1
z(k)=Aα(k+1)+η(k+1)−λ(k)/γ, bm=∑Nn=1y∗(n−1)M+mz(k)(n−1)M+m, am=∑Nn=1|y(n−1)M+m|2,
β=[δ+jw−∑Mm=1(bm/am)]/∑Mm=1(1/am)。步骤4-2 t(k+1)=[(b1+β)/a1,(b2+β)/a2,⋯,(bM+β)/aM]T; 步骤5 λ(k+1)=λ(k)−γ(Aα(k+1)+η(k+1)−T(k+1)y); 步骤6 令k←k+1,若k≤K−1,则返回步骤2,否则结束。 输出: α=αK, em=1/tKm, e=[e1,e2,⋯,eM]T。 ADMM属于模型驱动类算法,其正则化因子
ρ 、二次惩罚因子γ 、迭代步长τ 、比例因子δ 和w等参数均需提前给定。在实际应用中,参数的设置是比较困难的。不恰当的参数设置会影响ADMM算法的收敛速度和精度,从而使式(13)的求解复杂度升高、杂波空时谱和阵列误差参数的估计准确性下降。即使能够通过理论分析、交叉验证的方法选择合适的参数,固定的参数设置并不能保证ADMM算法获得最好的收敛效果。为解决上述问题,基于DU方法的思路,本文将ADMM算法展开为深度神经网络AE-ADMM-Net,利用学习的方法获得其最优迭代参数。为构建AE-ADMM-Net,将ADMM算法的迭代步骤映射为一个数据流图,如图2所示。图2所示数据流图主要由ADMM算法所对应的不同图节点和不同图节点之间表示数据流动的有向边组成。数据流图的第k+1层表示ADMM算法的第k+1次迭代,表1的迭代步骤2—步骤5对应4个图节点:辅助变量更新节点(
X(k+1) )、杂波空时谱更新节点(O(k+1) )、误差参数更新节点(Z(k+1) )和拉格朗日乘子更新节点(M(k+1) )。可以看出:ADMM算法的K次迭代可以映射为一个K层的数据流图,输入的空时回波信号沿此数据流图进行传递,将获得杂波空时谱和阵列幅相误差的估计结果。3.2 AE-ADMM-Net
对于式(13)所示的优化问题,当机载雷达参数给定且杂波复幅度、阵列误差和噪声均服从一定分布时,训练距离单元的空时回波信号
yl 也将具有一定分布。此外,给定空时导向矢量字典A,杂波空时谱αl 也将具有一定稀疏分布。此时,可假设存在一组最优的参数序列,使得对于所有服从一定分布的空时回波信号、杂波空时谱和阵列误差,ADMM算法均能够快速准确地求解式(13)。因此,为解决ADMM算法存在的问题,结合模型驱动算法的可解释性和数据驱动深度学习方法的非线性拟合能力,本节基于ADMM算法的迭代步骤和数据流图,构建AE-ADMM-Net,将其用于求解式(13)。基于充足完备的训练数据集对AE-ADMM-Net进行训练,能够获得最优的迭代参数,从而提高杂波空时谱和阵列误差参数的估计速度和性能。下面对AE-ADMM-Net的网络结构、数据集构建方法、网络初始化与训练进行具体描述。3.2.1 网络结构
根据表1所示的算法步骤和图2所示的数据流图,可将ADMM算法等效为一个如图3所示的K层网络AE-ADMM-Net,其输入为y, A,
α(0) ,λ(0) 和t(0) ,可学习参数为Θ={Θ(k+1)}K−1k=0={ρk+1,γk+1,τk+1,δk + 1,ωk+1 ,ςk+1}K−1k=0 ,输出为α(K) 和t(K) ,从而可得杂波空时谱α=αK 和阵列误差em=1/tKm ,e=[e1,e2,⋯,eM]T 。其中,AE-ADMM-Net的第k+1层运算可表示为{α(k+1),λ(k+1),t(k+1)}=Fk+1{y,A,α(k),λ(k),t(k),Θ(k+1)} (22) 其中,
Fk+1{⋅} 对应一个4层子网络,包括辅助变量更新层(X(k+1) )、杂波空时谱更新层(O(k+1) )、误差参数更新层(Z(k+1) )和拉格朗日乘子更新层(M(k+1) ),如图4所示,其中实箭头表示正向传播的方向,虚线箭头表示反向传播的方向,具体描述如下:(1) 辅助变量更新层(
X(k+1) ):将y, A以及AE-ADMM-Net第k层中O(k) ,Z(k) 和M(k) 的输出α(k) ,t(k) 和λ(k) 作为其输入,则X(k+1) 的输出为η(k+1)=ρk+1γk+1/(1+ρk+1γk+1)[λ(k)/γk+1−Aα(k)+(IN⊗diag(t(k)))y] (23) 其中,
ρk+1 和γk+1 为第k+1层可学习的正则化因子和二次惩罚因子。X(k+1) 的输出η(k+1) 将作为第k+1层中O(k+1) ,Z(k+1) ,M(k+1) 的输入。(2) 杂波空时谱更新层(
O(k+1) ):将A以及AE-ADMM-Net第k层中O(k) 的输出α(k) 和第k+1层中X(k+1) 的输出η(k+1) 作为其输入,则O(k+1) 的输出为α(k+1)=soft(α(k)+τk+1AHη(k+1)/(ρk+1γk+1),τk+1/γk+1) (24) 其中,
τk+1 为第k+1层可学习的迭代步长。O(k+1) 的输出α(k+1) 将作为第k+1层中M(k+1) 和Z(k+1) 以及第k+2层中X(k+2) 和O(k+2) 的输入。(3) 误差参数更新层(
Z(k+1) ):将y,A,AE-ADMM-Net第k层中M(k) 的输出λ(k) , AE-ADMM-Net第k+1层中X(k+1) 和O(k+1) 的输出η(k+1) 和α(k+1) 作为其输入,则Z(k+1) 的输出为t(k+1)=[(b1+β)/a1,(b2+β)/a2,⋯,(bM+β)/aM]T (25) 其中,
z(k)=Aα(k+1)+η(k+1)−λ(k)/γk+1 ,bm=∑Nn=1y∗(n−1)M+mz(k)(n−1)M+m ,am=∑Nn=1|y(n−1)M+m|2 ,β=[δk+1+ jwk+1−∑Mm=1(bm/am)]/∑Mm=1(1/am) ,δk+1 和wk+1 为第k+1层可学习的比例因子。Z(k+1) 的输出t(k+1) 将作为第k+1层中M(k+1) 以及第k+2层中X(k+2) 的输入。(4) 拉格朗日乘子更新层(
M(k+1) ):将y, A, AE-ADMM-Net第k层中M(k) 的输出λ(k) , AE-ADMM-Net第k+1层中X(k+1) ,O(k+1) 和Z(k+1) 的输出η(k+1) ,α(k+1) 和t(k+1) 作为其输入,则M(k+1) 的输出为λ(k+1)=λ(k)−ςk+1[Aα(k+1)+η(k+1)−(IN⊗diag(t(k+1)))y] (26) 其中,
ςk+1 为第k+1层可学习的乘子更新参数。M(k+1) 的输出λ(k+1) 将作为第k+2层中M(k+2) ,X(k+2) 和Z(k+2) 的输入。需要强调的是:相比利用γk+1 作为乘子更新参数(如式(21)所示),添加新参数ςk+1 是为了进一步增强网络的学习能力,提高AE-ADMM-Net的性能。3.2.2 数据集构建方法
与现有DU方法相同,本文AE-ADMM-Net是一种“模型+数据”联合驱动的SR方法,合理构建具有泛化能力的数据集是决定其有效性的关键。此外,DU方法大多采用监督学习的方式,按照提前给定的数据及其标签对网络进行训练。为了使空时回波信号、杂波空时谱和阵列幅相误差均具有一定的分布,本文构建数据集的方式可以概括为“设定雷达参数、设定杂波分布、设定阵列幅相误差分布、生成空时回波信号、划分训练与测试数据集、构造空时导向矢量字典、获得训练和测试标签集”,具体描述如下:
步骤1 对于机载雷达正侧视均匀线阵,设定载机高度H、载机速度v、阵元数M、脉冲数N、阵元间距d、波长
λ 、脉冲重复频率fr和距离范围[Rmin, Rmax]等参数;步骤2 根据雷达距离分辨率将距离范围划分为L个距离单元,将每个距离单元所对应的距离环在方位角
θ∈[0,π] 上划分为Nc个杂波块,杂波块之间相互独立且幅度服从复高斯分布;步骤3 令阵元幅度误差
νm 和相位误差ϕm 分别服从[−νmax,νmax] 和[−ϕmax,ϕmax] 上的均匀分布,随机产生P个阵元误差矢量{ep}Pp=1 ,其中ep=[ep1,ep2,⋯,epM]T ,epm=(1+νpm)ejϕpm ,νpm∈U(−νmax,νmax) ,ϕpm∈U(−ϕmax,ϕmax) ,νmax 和ϕmax 分别表示幅度误差和相位误差的最大值;步骤4 对于每个阵元误差矢量
ep ,根据yl,p=∑Nci=1αl,isdl,i⊗(ep⊙ssl,i)+εl 产生L个空时回波信号{yl,p}Ll=1 ,其中,sdl,i 和ssl,i 分别为第l个距离单元上第i个杂波块的时域导向矢量和空域导向矢量,对应的复幅度为αl,i ,εl 为复高斯白噪声,杂噪比为CNR;步骤5 将
PL 个空时回波信号{{yl,p}Ll=1} Pp=1 随机划分为包含O个空时回波信号的训练数据集{ytraino}Oo=1 和包含S=(PL−O) 个空时回波信号的测试数据集{ytests}Ss=1 ;步骤6 设定空间频率和多普勒频率范围
[fsmin,fsmax] 和[fdmin,fdmax] 、网格数Ns=κsM 和Nd=κdN ,构造空时导向矢量字典A=[sds1,sds2,⋯,sdsNsNd] ;步骤7 通过理论分析和交叉验证的方式设置ADMM算法参数
ρ=ρ0 ,γ=γ0 ,τ=τ0 ,δ=δ0 ,ω=ω0 和K=K0 ,对式(13)进行求解,获得训练标签集{αtraino,etraino}Oo=1 和测试标签集{αtests,etests}Ss=1 。具体步骤为:基于理论分析,得到算法收敛时ρ ,γ 和τ 需满足的条件[26,27];在满足收敛条件的数值范围内,设置多组不同的ρ ,γ 和τ ;由于假设实际阵列误差服从均匀分布,设置δ=M ,ω=0 不变;基于不同的参数组合对空时回波信号进行处理,终止迭代的条件设为第K次迭代结果相对第K–1次迭代结果的归一化误差小于10–6;得到对于所有空时回波信号杂波空时谱估计均较为准确、杂波抑制性能均较优的一组参数ρ0 ,γ0 ,τ0 ,δ0 ,ω0 和K0 ,作为ADMM算法的参数,并获得其对应的训练和测试标签集。3.2.3 初始化与训练
网络的初始化和训练方法对AE-ADMM-Net的性能具有一定的影响,较好的初始化和训练方法能够使网络更容易达到收敛,在一定程度上避免陷入局部最优。AE-ADMM-Net的参数可根据3.2.2节中的步骤(7)进行初始化,即令
ρ1:K=ρ0 ,γ1:K=γ0 ,τ1:K=τ0 ,δ1:K=δ0 ,ω1:K=ω0 和ς1:K=γ0 。与采用固定参数设置的ADMM算法相比,AE-ADMM-Net经过训练后将在保证收敛性能的基础上,大幅提高收敛速度(即减少迭代次数),缩短求解式(13)的时间。基于所构建的训练数据集
{αtraino,etraino,ytraino}Oo=1 ,给定网络层数K,定义归一化均方根误差(Normalized Mean Square Error, NMSE)作为网络损失函数,则AE-ADMM-Net的最优参数Θ∗={ρ∗k+1,γ∗k+1,τ∗k+1,δ∗k+1 ,ω∗k+1,ς∗k+1}K−1k=0 可通过后向传播(Back Propagation, BP)方法[28]求解下式得到:Θ∗=argminΘ1OO∑o=10.5Loα+0.5Loe (27) 其中
{Loα=||α(K)(Θ,A,α(0),t(0),λ(0),ytraino)−αtraino||22/||αtraino||22Loe=||e(K)(Θ,A,α(0),t(0),λ(0),ytraino)−etraino||22/||etraino||22 (28) α(K)(ytraino,A,α(0),λ(0),t(0),Θ) 表示以ytraino , A,α(0)=0NdNs ,λ(0)=0NM 和t(0)=1M 为输入、以Θ 为参数的AE-ADMM-Net第K层中杂波空时谱更新层的输出,e(K)(Θ,A,α(0),t(0),λ(0),ytraino) 对应第K层中误差参数更新层的输出t(K)(Θ,A,α(0),t(0),λ(0),ytraino) ,满足e(K)=[e(K)1,e(K)2,⋯,e(K)M] ,e_m^{(K)} = 1/t_m^{(K)} 。经过训练得到最优参数后,即可将AE-ADMM-Net应用于实际训练距离单元空时回波信号的处理。对于测试数据
\{ {\boldsymbol{y}}_s^{{\text{test}}}\} _{s = 1}^S ,其杂波空时谱和阵列误差参数的估计可以表示为\left\{ \begin{gathered} \hat {\boldsymbol{\alpha}} _s^{{{\rm{test}}}} = {{\boldsymbol{\alpha}} ^{(K)}}({\boldsymbol{y}}_s^{{{\rm{test}}}},{\boldsymbol{A}},{{\boldsymbol{\alpha}} ^{(0)}},{{\boldsymbol{t}}^{(0)}},{{\boldsymbol{\lambda}} ^{(0)}},{{\boldsymbol{\varTheta}} ^ * }) \\ \hat {\boldsymbol{e}}_s^{{{\rm{test}}}} = {{\boldsymbol{e}}^{(K)}}({\boldsymbol{y}}_s^{{{\rm{test}}}},{\boldsymbol{A}},{{\boldsymbol{\alpha}} ^{(0)}},{{\boldsymbol{t}}^{(0)}},{{\boldsymbol{\lambda}} ^{(0)}},{{\boldsymbol{\varTheta}} ^ * }) \\ \end{gathered} \right. (29) 其中,
{{\boldsymbol{\alpha}} ^{(K)}}({\boldsymbol{y}}_s^{{{\rm{test}}}},{\boldsymbol{A}},{{\boldsymbol{\alpha}} ^{(0)}},{{\boldsymbol{t}}^{(0)}},{{\boldsymbol{\lambda}} ^{(0)}},{{\boldsymbol{\varTheta}} ^ * }) 表示以{\boldsymbol{y}}_s^{{{\rm{test}}}} , A,{{\boldsymbol{\alpha}} ^{(0)}} = {{\bf{0}}_{{N_{\rm{d}}}{N_{\rm{s}}}}} ,{{\boldsymbol{\lambda}} ^{(0)}} = {{\bf{0}}_{NM}} 和{{\boldsymbol{t}}^{(0)}} = {{\bf{1}}_M} 为输入、以{{\boldsymbol{\varTheta}} ^ * } 为参数的AE-ADMM-Net第K层中杂波空时谱更新层,{{\boldsymbol{e}}^{(K)}}({\boldsymbol{y}}_s^{{{\rm{test}}}},{\boldsymbol{A}},{{\boldsymbol{\alpha}} ^{(0)}},{{\boldsymbol{t}}^{(0)}},{{\boldsymbol{\lambda}} ^{(0)}},{{\boldsymbol{\varTheta}} ^ * }) 对应第K层中误差参数更新层的输出{{\boldsymbol{t}}^{(K)}}({\boldsymbol{y}}_s^{{{\rm{test}}}},{\boldsymbol{A}},{{\boldsymbol{\alpha}} ^{(0)}},{{\boldsymbol{t}}^{(0)}},{{\boldsymbol{\lambda}} ^{(0)}}, {{\boldsymbol{\varTheta}} ^ * }) 。4. 仿真实验
本节通过仿真对基于AE-ADMM-Net的DU-STAP方法进行验证,并与基于SBL,FOCUSS和ADMM等算法的典型SR-STAP方法进行对比分析,仿真参数如表2所示。所有仿真均基于MATLAB R2020b实现,系统配置为Intel(R) Core(TM) i9-10900K CPU @ 3.70 GHz和NVIDIA GeForce RTX 2080 Ti GPU。
表 2 仿真参数Table 2. Simulation parameters参数 数值 参数 数值 载机高度H 3000 m 载机速度v 100 ms–1 阵元数M 10 个 脉冲数N 10 个 阵元间距d 0.1 m 工作波长λ 0.2 m 脉冲重复频率fr 2000 Hz 距离范围[Rmin, Rmax] [21,31] km 距离单元数L 100 个 杂波块数Nc 361 个 阵元误差数P 100 个 杂噪比CNR 60 dB 训练数据集大小O 7500 测试数据集大小S 2500 频率范围f s和f d [–0.5,0.5] 网格数Ns和Nd 50 个 为验证所提方法在不同阵列误差条件下的性能,令阵列幅相误差的最大值
({\nu_{\max }},{\phi _{\max }}) 分别等于(0,{0^\circ }) ,(0.1,{10^\circ }) ,(0.2,{20^\circ }) 和(0.3,{30^\circ }) ,按照3.2.2节步骤(1)—步骤(5)所述方法构建4组不同的数据集。然后,设置ADMM算法的迭代参数为{\rho _0} = 0.5 ,{\gamma _0} = 0.01 ,{\tau _0} = 0.04 ,{\delta _0} = M ,{\omega _0} = 0 和{K_0} = 3000 ,按照3.2.2节步骤(6)—步骤(7)所述方法构建标签集。图5给出了不同阵列误差条件下,利用ADMM算法对某一数据进行处理得到的杂波空时谱和阵列误差参数估计,其从左到右分别对应({\nu_{\max }},{\phi _{\max }}) 等于(0,{0^\circ }) ,(0.1,{10^\circ }) ,(0.2,{20^\circ }) 和(0.3,{30^\circ }) 的情况,从上到下分别为杂波空时谱、阵列幅度误差和阵列相位误差的估计。可以看出,基于上述固定参数,ADMM算法在不同条件下均能获得较为准确的估计结果,因此可利用所构建的数据集对AE-ADMM-Net进行训练。图 5 固定参数ADMM算法杂波空时谱和阵列误差参数估计结果(a1—a4:不同阵列误差参数下的空时谱估计结果,b1—b4:不同阵列误差参数下的幅度误差估计结果,c1—c4:不同阵列误差参数下的相位误差估计结果)Figure 5. Clutter space-time spectra and array error parameters estimated via ADMM algorithm with fixed parameters (a1—a4: Clutter space-time spectra estimation results in different array error parameters, b1—b4: Amplitude error estimation results in different array error parameters, c1—c4: Phase error estimation results in different array error parameters)4.1 网络收敛性
本节验证AE-ADMM-Net的收敛性,并与固定迭代参数的ADMM算法进行对比分析。设置不同的网络层数K,按照3.2.3节所述方法对AE-ADMM-Net进行初始化和训练(Adam算法,迭代次数为3000),所得结果如图6所示。其中,图6(a),图6(b)为AE-ADMM-Net在网络层数K=25时的训练NMSE和测试NMSE,图6(c)为AE-ADMM-Net和ADMM算法在网络层数(迭代次数)K=15~45时的NMSE,图6(d)为ADMM算法在迭代次数K=60~180时的NMSE。从图6(a),图6(b)可以看出,无论是否存在阵列误差,AE-ADMM-Net的训练和测试NMSE均随着训练次数的增加而逐渐下降,且在训练1500次后基本达到收敛。从图6(c)可以看出,随着网络层数(迭代次数)的增加,AE-ADMM-Net和ADMM算法的NMSE均逐渐下降,但前者的NMSM远小于后者。从图6(c),图6(d)可以看出,当ADMM算法的迭代次数为AE-ADMM-Net的4倍时,两者才具有相近的NMSE。因此,可以得出结论:无论是否存在阵列误差,AE-ADMM-Net均能够从所构建的数据集中学习得到最优迭代参数,获得更好的收敛性能。需要说明的是:当网络层数达到一定数值(35~40)时,AE-ADMM-Net就可以获得比较准确的杂波空时谱估计结果,进一步增加网络层数并不能显著提高杂波抑制性能,反而会增加运算复杂度。因此,在实际应用中,可基于不同的仿真条件对AE-ADMM-Net进行离线训练,确定可获得较好杂波抑制性能和较低运算复杂度的网络层数取值范围,再根据实际情况进行选择。
4.2 杂波空时谱
本节验证AE-ADMM-Net的杂波空时谱估计性能,并与FOCUSS算法、SBL算法和固定迭代参数的ADMM算法进行对比分析。图7给出了不同阵列误差条件下,利用不同算法对图5所对应的数据进行处理获得的杂波空时谱估计结果,其从左到右分别对应
({\nu_{\max }},{\phi _{\max }}) 等于(0,{0^\circ }) ,(0.1,{10^\circ }) ,(0.2,{20^\circ }) 和(0.3,{30^\circ }) 的情况,从上到下分别对应迭代25次的ADMM算法、迭代45次的ADMM算法、迭代200次的FOCUSS算法(正则化参数设为10–3)、迭代400次的SBL算法(噪声功率初始值设为10–6)、层数为25的AE-ADMM-Net和层数为45的AE-ADMM-Net。可以看出:(1)与图5相比,固定迭代参数的ADMM算法在迭代次数较少时难以获得准确的杂波空时谱估计;(2)在不存在阵列误差时,SBL算法和FOCUSS算法均能够获得杂波空时谱的准确估计,但存在阵列误差时估计性能急剧下降;(3)无论是否存在阵列误差,AE-ADMM-Net均能够基于少量网络层数(迭代次数),实现对杂波空时谱的准确估计。因此,可以得出结论:相比典型的SR算法,AE-ADMM-Net在不同条件下均能快速获得杂波空时谱的准确估计。图 7 不同条件下不同算法的杂波空时谱估计结果(a1—a4:ADMM算法在不同阵列误差参数下的迭代25次的估计结果,b1—b4:ADMM算法在不同阵列误差参数下的迭代45次的估计结果,c1—c4:FOCUSS算法在不同阵列误差参数下的迭代200次的估计结果,d1—d4:SBL算法在不同阵列误差参数下的迭代400次的估计结果,e1—e4:25层的AE-ADMM-Net 在不同阵列误差参数下的的估计结果,f1—f4:45层的AE-ADMM-Net 在不同阵列误差参数下的估计结果)Figure 7. Clutter space-time spectra estimated via different algorithms under different conditions (a1—a4: estimation results of ADMM algorithm with 25 iterations in different array error parameters, b1—b4: estimation results of ADMM algorithm with 45 iterations in different array error parameters, c1—c4: estimation results of FOCUSS algorithm with 200 iterations in different array error parameters, d1—d4: estimation results of SBL algorithm with 400 iterations in different array error parameters, e1—e4: estimation results of AE-ADMM-Net with 25 layers in different array error parameters, f1—f4: estimation results of AE-ADMM-Net with 45 layers in different array error parameters)4.3 阵列误差参数
本节验证AE-ADMM-Net的阵列误差参数估计性能,结果如图8所示。图8从左到右分别对应
({\nu_{\max }},{\phi _{\max }}) 等于(0,{0^\circ }) ,(0.1,{10^\circ }) ,(0.2,{20^\circ }) 和(0.3,{30^\circ }) 的情况,上图和下图分别为幅度误差和相位误差估计结果。可以看出:在不同条件下,AE-ADMM-Net均能获得阵列幅度误差和相位误差的准确估计。图 8 不同条件下AE-ADMM-Net的阵列误差参数估计结果(a1—a4:不同阵列误差参数下的幅度误差估计结果,b1—b4:不同阵列误差参数下的相位误差估计结果)Figure 8. Array error parameters estimated by AE-ADMM-Net under different conditions (a1—a4: Amplitude error estimation results in different array error parameters, b1—b4: Phase error estimation results in different array error parameters)4.4 SCNR损失
本节验证基于AE-ADMM-Net的DU-STAP方法的杂波抑制性能,并与基于FOCUSS算法、SBL算法和固定迭代参数ADMM算法的SR-STAP方法进行对比分析。需要说明的是:由于SBL算法和FOCUSS算法无法有效估计阵列误差参数,因此在进行性能对比分析时,不考虑阵列误差参数,仅对不同算法得到的杂波空时谱进行处理,即基于式(11)估计CNCM
{\hat {\boldsymbol{R}}_{{\rm{C}}}} 、基于式(3)计算空时滤波器最优权值{{\boldsymbol{w}}_{{{\rm{opt}}}}} 。然后,利用SCNR损失衡量不同方法的杂波抑制性能,表示为{{{\rm{SCNR}}}_{{{\rm{Loss}}}}} = \frac{{{\sigma ^2}|{\boldsymbol{w}}_{{{\rm{opt}}}}^{{\rm{H}}}{\boldsymbol{s}}_{{\rm{T}}}^{{{\rm{ds}}}}{|^2}}}{{NM{\text{ }}{\boldsymbol{w}}_{{{\rm{opt}}}}^{{\rm{H}}}{{\hat {\boldsymbol{R}}}_{{\rm{C}}}}{{\boldsymbol{w}}_{{{\rm{opt}}}}}}} (30) 假设目标的空间频率为0(即
{\boldsymbol{s}}_{\rm{T}}^{\rm{s}} = {{\bf{1}}_M} )、归一化多普勒频率在[–0.5,0.5]范围内变化,不同方法对应的SCNR损失曲线如图9所示,其从左到右分别对应({\nu_{\max }},{\phi _{\max }}) 等于(0,{0^\circ }) ,(0.1,{10^\circ }) ,(0.2,{20^\circ }) 和(0.3,{30^\circ }) 的情况,下图对应上图的局部放大结果。可以看出:(1)基于FOCUSS和SBL算法的SR-STAP方法仅在无阵列误差时有效,在存在阵列误差时杂波抑制性能急剧下降;(2)基于固定迭代参数ADMM算法的SR-STAP方法在迭代次数较多的条件下(ADMM-opt, K = 3000)能够有效抑制杂波,但在迭代次数较少的条件下(K = 25和45),由于杂波空时谱估计不准确,其杂波抑制性能较差;(3)基于AE-ADMM-Net的DU-STAP方法基于少量网络层数(迭代次数)即可获得杂波空时谱的准确估计,实现对杂波的有效抑制,网络层数为K=45时的性能与ADMM-opt相当。因此,可以得出结论:相比典型的SR-STAP方法,DU-STAP方法在不同条件下均能获得较好的杂波抑制性能。图 9 不同条件下不同方法对应的SCNR损失曲线(a1—a4:不同阵列误差参数下的SCNR曲线结果,b1—b4:不同阵列误差参数下的SCNR曲线局部放大结果)Figure 9. SCNR loss curves corresponding to different methods under different conditions (a1—a4: SCNR loss curves results in different array error parameters, b1—b4: SCNR loss curves results with enlarged scale in different array error parameters)4.5 运算复杂度
本节分析AE-ADMM-Net的运算复杂度,并与FOCUSS算法和SBL算法进行对比。需要强调的是:由于可以采用离线训练、在线应用的方法[25,29],本文对AE-ADMM-Net的运算复杂度分析不包括网络训练所需的运算量。此外,在进行训练获得最优参数后,AE-ADMM-Net与ADMM算法的运算完全相同,仅在迭代参数上具有差异。因此,在网络层数(迭代次数)相同的条件下,AE-ADMM-Net与ADMM算法在应用时将具有相同的运算复杂度。以乘法次数为指标,可得不同算法进行一次迭代所需的运算复杂度如表3所示。可以看出,AE-ADMM-Net的运算复杂度远小于FOCUSS算法和SBL算法。为了对此进行验证,基于MATLAB的TIC和TOC命令获得不同条件下AE-ADMM-Net, FOCUSS和SBL算法的运行时间如图10所示。其中,图10(a)对应M=N=10、Nd=Ns=50、迭代次数K=15~45;图10(b)对应M=N=4~16、Nd=Ns=50、迭代次数K=45;图10(c)对应M=N=10、Nd=Ns=20~80、迭代次数K=45;图10(d)对应M=N=Nd/5=Ns/5=4~16、迭代次数K=45。可以看出:在不同条件下,AE-ADMM-Net的运行时间均远小于FOCUSS算法和SBL算法。此外,需要指出的是:与ADMM算法相似,参数固定的FOCUSS算法和SBL算法通常也需要相比AE-ADMM-Net更多的迭代次数以达到收敛。因此,可以得出结论:相比基于FOCUSS和SBL算法的SR-STAP算法,基于AE-ADMM-Net的DU-STAP方法具有更低的运算复杂度。
表 3 不同算法的运算复杂度Table 3. Computational complexities of different algorithms算法 运算复杂度 FOCUSS O\left( {3NM{N_{{\rm{s}}} }{N_{{\rm{d}}} } + { {(NM)}^3} + 2{ {(NM)}^2}{N_{{\rm{s}}} }{N_{{\rm{d}}} } } \right) SBL O\left( {5NM{N_{{\rm{s}}} }{N_{{\rm{d}}} } + { {(NM)}^3} + 2{ {(NM)}^2}{N_{{\rm{s}}} }{N_{{\rm{d}}} } + NM + {N_{{\rm{s}}} }{N_{{\rm{d}}} } } \right) AE-ADMM-Net O\left( {2NM{N_{{\rm{s}}} }{N_{{\rm{d}}} } + { {(NM)}^2} + 2NM + {N_{{\rm{s}}} }{N_{{\rm{d}}} } } \right) 4.6 实测数据处理
本节基于Mountain Top实测数据[16]对所提DU-STAP方法的实际性能进行验证,并与基于固定迭代参数ADMM算法的SR-STAP方法进行对比分析,其中ADMM算法的参数设置与仿真实验一致,DU-STAP方法直接采用由仿真数据训练得到的AE-ADMM-Net。Mountain Top数据的阵元数为14、脉冲数为16,目标位于第147个距离单元,为与仿真相匹配,取10个阵元和10个脉冲所对应的数据进行处理。假设不存在阵元误差并设保护距离单元个数为4,基于ADMM和AE-ADMM-Net对第152个距离单元的空时回波信号进行处理,从而估计杂波空时谱、设计空时滤波器进行杂波抑制和目标检测,结果如图11所示。其中,前3个子图依次对应迭代3000次的ADMM算法、迭代45次的ADMM算法和网络层数为45的AE-ADMM-Net,第4个子图为目标检测结果。可以看出,本文所提DU-STAP方法对实测数据进行处理仍能获得较好的结果,在迭代次数相同的条件下,杂波空时谱估计和目标检测性能均优于基于固定迭代参数ADMM算法的SR-STAP方法。
5. 结语
本文提出了基于DU的机载雷达STAP方法。在存在阵列误差的条件下,对基于ADMM算法的杂波空时谱和阵列误差联合估计方法进行了分析,针对其存在的问题构建了深度神经网络AE-ADMM-Net,并对其网络结构、数据集构建方法、网络初始化与训练方法进行了介绍。通过仿真实验对基于AE-ADMM-Net的DU-STAP方法进行了验证,结果表明:相比典型的SR算法,AE-ADMM-Net能够从数据中学习得到最优迭代参数,在不同阵列误差条件下快速获得杂波空时谱和阵列误差参数的准确估计;相比典型的SR-STAP方法,DU-STAP方法能够获得较好的杂波抑制性能,且运算复杂度更低。下一步将对载机偏航、距离模糊、杂波内部运动和网格失配等非理想条件下的算法改进与分析进行深入研究。
-
-
[1] 韩金旺, 张子敬, 刘军, 等. 基于贝叶斯的高斯杂波背景下MIMO雷达自适应检测算法[J]. 雷达学报, 2019, 8(4): 501–509. doi: 10.12000/JR18090HAN Jinwang, ZHANG Zijing, LIU Jun, et al. Adaptive Bayesian detection for MIMO radar in Gaussian clutter[J]. Journal of Radars, 2019, 8(4): 501–509. doi: 10.12000/JR18090 [2] 王永良, 刘维建, 谢文冲, 等. 机载雷达空时自适应检测方法研究进展[J]. 雷达学报, 2014, 3(2): 201–207. doi: 10.3724/SP.J.1300.2014.13081WANG Yongliang, LIU Weijian, XIE Wenchong, et al. Research progress of space-time adaptive detection for airborne radar[J]. Journal of Radars, 2014, 3(2): 201–207. doi: 10.3724/SP.J.1300.2014.13081 [3] DE MAIO A. Rao test for adaptive detection in Gaussian interference with unknown covariance matrix[J]. IEEE Transactions on Signal Processing, 2007, 55(7): 3577–3584. doi: 10.1109/TSP.2007.894238 [4] PASCAL F, CHITOUR Y, OVARLEZ J P, et al. Covariance structure maximum-likelihood estimates in compound Gaussian noise: Existence and algorithm analysis[J]. IEEE Transactions on Signal Processing, 2008, 56(1): 34–48. doi: 10.1109/TSP.2007.901652 [5] 许述文, 石星宇, 水鹏朗. 复合高斯杂波下抑制失配信号的自适应检测器[J]. 雷达学报, 2019, 8(3): 326–334. doi: 10.12000/JR19030XU Shuwen, SHI Xingyu, and SHUI Penglang. An adaptive detector with mismatched signals rejection in compound Gaussian clutter[J]. Journal of Radars, 2019, 8(3): 326–334. doi: 10.12000/JR19030 [6] GRECO M, GINI F, and RANGASWAMY M. Statistical analysis of measured polarimetric clutter data at different range resolutions[J]. IEE Proceedings - Radar, Sonar and Navigation, 2006, 153(6): 473–481. doi: 10.1049/ip-rsn:20060045 [7] SANGSTON K J, GINI F, and GRECO M S. Coherent radar target detection in heavy-tailed compound-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2012, 48(1): 64–77. doi: 10.1109/TAES.2012.6129621 [8] HE You, JIAN Tao, SU Feng, et al. Adaptive detection application of covariance matrix estimator for correlated non-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(4): 2108–2117. doi: 10.1109/TAES.2010.5595620 [9] CONTE E, LOPS M, and RICCI G. Asymptotically optimum radar detection in compound-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 1995, 31(2): 617–625. doi: 10.1109/7.381910 [10] STINCO P, GRECO M, and GINI F. Adaptive detection in compound-Gaussian clutter with inverse-gamma texture[C]. 2011 IEEE CIE International Conference on Radar, Chengdu, China, 2011. [11] LI Na, CUI Guolong, KONG Lingjiang, et al. Rao and Wald tests design of multiple-input multiple-output radar in compound-Gaussian clutter[J]. IET Radar, Sonar & Navigation, 2012, 6(8): 729–738. doi: 10.1049/iet-rsn.2011.0376 [12] LIU Weijian, LIU Jun, HUANG Lei, et al. Performance analysis of reduced-dimension subspace signal filtering and detection in sample-starved environment[J]. Journal of the Franklin Institute, 2019, 356(1): 629–653. doi: 10.1016/j.jfranklin.2018.10.017 [13] YAN Linjie, HAO Chengpeng, ORLANDO D, et al. Parametric space-time detection and range estimation of point-like targets in partially homogeneous environment[J]. IEEE Transactions on Aerospace and Electronic Systems, 2020, 56(2): 1228–1242. doi: 10.1109/TAES.2019.2928672 [14] 顾新锋, 简涛, 何友, 等. 局部均匀背景中距离扩展目标的GLRT检测器及性能分析[J]. 电子学报, 2013, 41(12): 2367–2373. doi: 10.3969/j.issn.0372-2112.2013.12.008GU Xinfeng, JIAN Tao, HE You, et al. GLRT detector of range-spread target in local homogeneous background and its performance analysis[J]. Acta Electronica Sinica, 2013, 41(12): 2367–2373. doi: 10.3969/j.issn.0372-2112.2013.12.008 [15] 何友, 简涛, 苏峰, 等. 非高斯杂波协方差矩阵估计方法及CFAR特性分析[J]. 中国科学:信息科学, 2011, 41(1): 90–99.HE You, JIAN Tao, SU Feng, et al. CFAR assessment of covariance matrix estimators for non-Gaussian clutter[J]. Science in China:Information Sciences, 2011, 41(1): 90–99. [16] GAO Yongchan, LI Hongbin, and HIMED B. Adaptive subspace tests for multichannel signal detection in auto-regressive disturbance[J]. IEEE Transactions on Signal Processing, 2018, 66(21): 5577–5587. doi: 10.1109/TSP.2018.2869123 [17] XU D, ADDABBO P, HAO C, et al. Adaptive strategies for clutter edge detection in radar[J]. Signal Processing, 2021, 186: 108127. doi: 10.1016/j.sigpro.2021.108127 [18] CAROTENUTO V, DE MAIO A, ORLANDO D, et al. Adaptive radar detection using two sets of training data[J]. IEEE Transactions on Signal Processing, 2018, 66(7): 1791–1801. doi: 10.1109/TSP.2017.2778684 [19] LIU Weijian, ZHANG Zhaojian, LIU Jun, et al. Detection of a rank-one signal with limited training data[J]. Signal Processing, 2021, 186: 108120. doi: 10.1016/j.sigpro.2021.108120 [20] CONTE E, DE MAIO A, and RICCI G. CFAR detection of distributed targets in non-Gaussian disturbance[J]. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(2): 612–621. doi: 10.1109/TAES.2002.1008990 [21] CONTE E and DE MAIO A. Distributed target detection in compound-Gaussian noise with Rao and Wald tests[J]. IEEE Transactions on Aerospace and Electronic Systems, 2003, 39(2): 568–582. doi: 10.1109/TAES.2003.1207267 [22] 简涛, 何友, 苏峰, 等. 非高斯背景下基于动态阈值的距离扩展目标检测器[J]. 电子学报, 2011, 39(1): 59–63.JIAN Tao, HE You, SU Feng, et al. Range-Spread target detector with dynamic threshold for non-Gaussian clutter[J]. Acta Electronica Sinica, 2011, 39(1): 59–63. [23] HE You, JIAN Tao, SU Feng, et al. Novel range-spread target detectors in non-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(3): 1312–1328. doi: 10.1109/TAES.2010.5545191 [24] COLUCCIA A, FASCISTA A, and RICCI G. A novel approach to robust radar detection of range-spread targets[J]. Signal Processing, 2020, 166: 107223. doi: 10.1016/j.sigpro.2019.07.016 [25] XU Shuwen, XUE Jian, and SHUI Penglang. Adaptive detection of range-spread targets in compound Gaussian clutter with the square root of inverse Gaussian texture[J]. Digital Signal Processing, 2016, 56: 132–139. doi: 10.1016/j.dsp.2016.06.009 [26] XUE Jian, XU Shuwen, and SHUI Penglang. Near-optimum coherent CFAR detection of radar targets in compound-Gaussian clutter with inverse Gaussian texture[J]. Signal Processing, 2020, 166: 107236. doi: 10.1016/j.sigpro.2019.07.029 [27] LIU Weijian, LIU Jun, HAO Chengpeng, et al. Multichannel adaptive signal detection: Basic theory and literature review[J]. Science China Information Sciences, 2022, 65(2): 121301. doi: 10.1007/s11432-020-3211-8 [28] SHI Bo, HAO Chengpeng, HOU Chaohuan, et al. Parametric Rao test for multichannel adaptive detection of range-spread target in partially homogeneous environments[J]. Signal Processing, 2015, 108: 421–429. doi: 10.1016/j.sigpro.2014.10.007 [29] 刘维建, 王利才, 狄源水, 等. 自适应能量检测器及在失配信号检测中的应用[J]. 雷达学报, 2015, 4(2): 149–159. doi: 10.12000/JR14132LIU Weijian, WANG Licai, DI Yuanshui, et al. Adaptive energy detector and its application for mismatched signal detection[J]. Journal of Radars, 2015, 4(2): 149–159. doi: 10.12000/JR14132 [30] WANG Zhihang, HE Zishu, HE Qin, et al. Adaptive CFAR detectors for mismatched signal in compound Gaussian sea clutter with inverse Gaussian texture[J]. IEEE Geoscience and Remote Sensing Letters, 2021, 19: 3502705. doi: 10.1109/LGRS.2020.3047390 [31] LIU Weijian, LIU Jun, GAO Yongchan, et al. Multichannel signal detection in interference and noise when signal mismatch happens[J]. Signal Processing, 2020, 166: 107268. doi: 10.1016/j.sigpro.2019.107268 [32] LIU Weijian, LIU Jun, LI Hai, et al. Multichannel signal detection based on Wald test in subspace interference and Gaussian noise[J]. IEEE Transactions on Aerospace and Electronic Systems, 2019, 55(3): 1370–1381. doi: 10.1109/TAES.2018.2870445 [33] GAO Yongchan, JI Hongbing, and LIU Weijian. Persymmetric adaptive subspace detectors for range-spread targets[J]. Digital Signal Processing, 2019, 89: 116–123. doi: 10.1016/j.dsp.2019.03.007 [34] REED I S, MALLETT J D, and BRENNAN L E. Rapid convergence rate in adaptive arrays[J]. IEEE Transactions on Aerospace and Electronic Systems, 1974, AES-10(6): 853–863. doi: 10.1109/TAES.1974.307893 [35] NITZBERG R. Application of maximum likelihood estimation of persymmetric covariance matrices to adaptive processing[J]. IEEE Transactions on Aerospace and Electronic Systems, 1980, AES-16(1): 124–127. doi: 10.1109/TAES.1980.308887 [36] 丁昊, 薛永华, 黄勇, 等. 均匀和部分均匀杂波中子空间目标的斜对称自适应检测方法[J]. 雷达学报, 2015, 4(4): 418–430. doi: 10.12000/JR14133DING Hao, XUE Yonghua, HUANG Yong, et al. Persymmetric adaptive detectors of subspace signals in homogeneous and partially homogeneous clutter[J]. Journal of Radars, 2015, 4(4): 418–430. doi: 10.12000/JR14133 [37] LIU Jun, ORLANDO D, ADDABBO P, et al. SINR distribution for the persymmetric SMI beamformer with steering vector mismatches[J]. IEEE Transactions on Signal Processing, 2019, 67(5): 1382–1392. doi: 10.1109/TSP.2019.2892027 [38] 闫林杰, 郝程鹏, 殷超然, 等. 部分均匀环境下适用于空间对称线阵的修正广义似然比检测方法[J]. 雷达学报, 2021, 10(3): 443–452. doi: 10.12000/JR20140YAN Linjie, HAO Chengpeng, YIN Chaoran, et al. Modified generalized likelihood ratio test detection based on a symmetrically spaced linear array in partially homogeneous environments[J]. Journal of Radars, 2021, 10(3): 443–452. doi: 10.12000/JR20140 [39] LIU Jun, LIU Weijian, CHEN Bo, et al. Modified Rao test for multichannel adaptive signal detection[J]. IEEE Transactions on Signal Processing, 2016, 64(3): 714–725. doi: 10.1109/TSP.2015.2491892 [40] GERLACH K. Spatially distributed target detection in non-Gaussian clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 1999, 35(3): 926–934. doi: 10.1109/7.784062 期刊类型引用(6)
1. 徐朋豪,占群峰,熊抒豪,刘小梅. 利用高阶累积量相关性的地震波传感器阵列目标检测方法. 宜宾学院学报. 2024(06): 38-45 . 百度学术
2. 简涛,马颖亮,王海鹏,郭磊,魏广芬. 异常值个数未知下辅助数据自适应筛选方法. 雷达学报. 2024(05): 1049-1060 . 本站查看
3. 田晗,张宇,许姗姗,高永婵,许智文. 子空间不确定下多重假设AMF、Rao与Wald检测方法. 海军航空大学学报. 2024(05): 622-632 . 百度学术
4. 林志凯,刘维建,钱李昌,李槟槟,周必雷,张昭建,陈浩. 子空间信号失配下的加权自适应检测器. 系统工程与电子技术. 2023(07): 1974-1980 . 百度学术
5. 杨政,程永强,吴昊,黎湘,王宏强. 基于正交投影的子带信息几何雷达弱小目标检测方法. 雷达学报. 2023(04): 776-792 . 本站查看
6. 刘旭,魏广芬,简涛,罗沅,周战,朱智林. 基于平方根分解的干扰加噪声背景下距离扩展目标子空间检测器. 微电子学与计算机. 2023(07): 36-44 . 百度学术
其他类型引用(4)
-