作者中心
专家审稿
责编办公
编辑办公
Space-time Two-dimensional Clutter Suppression Method Based on Subarray Beam Pattern Design
-
摘要: 机载雷达接收端采用子阵处理时面临栅瓣杂波导致的复杂空时耦合分布，使主瓣波束方向存在多个由栅瓣杂波导致的检测性能凹口，目标检测性能恶化严重。针对此问题，该文首先分析了子阵处理中栅瓣杂波空时谱分布特性，并在此基础上提出了基于接收子阵波束设计的空时二维杂波抑制方法。该方法采用重叠子阵构型方案，通过子阵方向图设计形成在子阵间栅瓣区域处的宽凹口，实现对子阵间栅瓣区域杂波的预滤波。进一步构建子阵级空时处理器，由于栅瓣杂波已经在子阵内完成预滤波，避免了栅瓣杂波在空时二维平面上的耦合扩散，从而提高了杂波抑制和动目标检测性能。仿真结果表明，所提方法显著改善了栅瓣杂波区的输出信杂噪比损失性能。Abstract: Airborne radar receivers that utilize subarray processing face challenges owing to the complex space-time coupling distribution caused by grating-lobe clutter. This results in multiple performance notches in the main beam, which severely affects target detection performance. To address this issue, we analyze the characteristics of grating-lobe clutter distribution in subarray processing and propose an approach for space-time clutter suppression based on the design of a receiving subarray beam pattern. Our approach leverages an overlapping subarray scheme to form wide nulls in the regions between subarrays where grating-lobe clutter is prevalent through beam pattern design. This design facilitates grating-lobe clutter pre-filtering between subarrays. Furthermore, we develop a subarray-level space-time processor that avoids the grating-lobe clutter coupling diffusion in the space-time two-dimensional plane by performing clutter pre-filtering within each subarray. This strategy enhances clutter suppression and moving-target-detection capabilities. Simulation results verify that the proposed method can remarkably improve the output signal to clutter plus noise ratio loss performance in grating-lobe clutter regions.
-
1. 引言
近年来，由于国际和国内公共安全形势的严峻性，火车站、机场、地铁等公共安全场所的安检显得尤为重要[1]。金属探测器、X光安检机等传统安检手段因其各自的局限性，不适合或者不能高效地用于人体安检[2,3]。毫米波具有穿透能力强、对人体辐射小且分辨率高等先天优势，尤其适用于隐藏物体检测[4,5]，因此毫米波成像为人体安检提供了新的选择，并且具有广阔的应用前景。主动式毫米波雷达3维成像基于合成孔径雷达(Synthetic Aperture Radar, SAR)成像技术[6,7]，实现对人体快速3维成像并进行异物检测，毫米波3维成像安检系统受到了国内外广泛重视。
毫米波3维成像系统有多种实现方式，国内外研究比较深入且具备实际工程应用能力的主要是平面扫描成像和圆柱扫描成像[7]，平面扫描成像虽然具备3维分辨能力，但受限于波束扫描范围，因此不能通过一次扫描实现全方位的观测角度成像[1]。圆周扫描成像系统采用阵列SAR与圆迹SAR相结合的方法，从理论上解决了毫米波成像系统的3维分辨问题，并弥补了平面扫描方式的不足，在较小的系统代价条件下实现高性能3维成像[8]。目前国内研究重点多集中于圆周扫描成像系统，文献[6]给出了一种主动式毫米波圆周扫描3维成像的系统模型以及重构算法，并且搭建出了基于该模型的实际成像系统。文献[9]基于圆柱面扫描给出了一种基于巴克码稀疏采样和频域压缩感知的毫米波人体安检3维成像方法，并设计了一种基于巴克码和收发分置模式的稀疏平面阵列，得到了较好的3维成像实时处理效果。本文针对圆周扫描3维成像系统的特点，结合上述研究背景，分析了3维ω-K算法和3维后向投影(Back Projection, BP)算法在工程应用中的限制条件，提出了一种将ω-K算法和BP算法相融合的3维图像重构算法。
2. 现有成像算法分析
如图1所示，在主动式圆柱扫描毫米波3维成像系统中，两个天线阵列以Z轴为轴心，以R为半径做圆周旋转，形成一个圆柱扫描面，其中一列线阵作为发射天线，另一列线阵作为接收天线，天线阵列位置与X轴的夹角定义为
θ ,θ∈[0,π] 。该系统在Y方向和X方向的分辨率依靠雷达的距离向分辨率和沿角度θ 方向的合成孔径获得。在天线阵列旋转的同时通过收发开关来控制天线的收发时序，从而在竖直方向形成合成孔径获得Z方向的高分辨，Z方向合成孔径的长度与波束的俯仰角宽正相关，假设波束俯仰角宽为S，则竖直方向的合成孔径长度为Lz=2Rtan(S/2) 。
图 1 圆柱扫描3维成像系统模型Figure 1. Three-dimensional imaging system model of cylindrical scanning对上述模式成像系统重构算法的选取，主要应该考虑算法对雷达运动模式的适用性和成像算法的聚焦性能。BP算法在成像过程中没有近似，成像精度较高，同时对雷达的运动方式没有特殊要求，并且不受转动角度的限制，因此非常适用于本系统的工作模式[10]。3维BP成像处理对于网格点数较多的3维场景成像运算量巨大，针对高实时要求的人体安检应用存在困难。文献[6,11]提供了一类基于3维频域(波数域)成像算法(3维ω-K算法)，可以快速实现3维成像，该类算法在计算3维逆傅里叶变换时要求在角度维均匀采样，非线性相位补偿比较复杂[12]。在柱面扫描条件下，要求天线阵列扫描轨迹为理想圆柱，一定程度限制了系统工作方式，对天线阵列运动轨迹要求严格。本文结合两种算法的优点，提出了ω-K算法和BP算法相结合的方式来实现快速柱面扫描方式的3维成像。在阵列(高度)方向与距离方向采用ω-K算法进行解耦合与聚焦，减少后向投影的运算量；对于每一个阵元，在距离方向和角度方向采用BP算法进行合成孔径处理，避免角度维插值处理，同时适用于非理想圆柱扫描方式。
3. 融合ω-K和BP 3维重构算法
设某一脉冲时刻，雷达天线位置为
(Rcosθ, Rsinθ,z′) ，成像区域第i 个目标点的位置为Pi(xi,yi,zi) ，散射系数为σi ，则成像区域第i 个目标点到雷达天线中心的斜距为：Ri=√(Rcosθ−xi)2+(Rsinθ−yi)2+(z′−zi)2 (1) 该时刻雷达天线发射的线性调频信号为：
st(t)=rect(ˆtTp)ej2π(fct+γˆt22) (2) 式中，
ˆt 为快时间，Tp 为脉冲宽度，fc 为信号载频，γ 为线性调频信号的调频率，则在(t,θ,z′) 域中表示的Pi 点的回波信号为：sr(t,θ,z′)=σirect(t−ΔtTp)ej2π(fc(t−Δt)+12γ(ˆt−Δt)2) (3) 式中，
Δt=2Ri/c 为雷达天线发射信号与成像区域目标点的双程时延。对式(3)进行解线频调(dechirp)处理，选取的参考信号为：sref(t,θ,z′)=rect(ˆt−2Rref/cTp)⋅ej2π(fc(t−2Rrefc)+12γ(ˆt−2Rrefc)2) (4) 式中，
Rref 为参考距离。参考信号与回波共轭相乘，即作差频处理，回波变成单频信号，且其频率与回波和参考信号的距离差成正比。若令RΔi=Ri− Rref ，则式(3)与式(4)共轭相乘后的差频信号为：srs(t,θ,z′)=sr(t,θ,z′)s∗ref(t,θ,z′)=σirect(ˆt−ΔtTp)e−j4πγc(ˆt−2Rrefc)RΔi⋅e−j4πcfcRΔiej4πγc2R2Δi (5) 实际上回波数据是由成像区域内的多个散射点的回波相叠加而成的，于是式(5)改写为：
s(t,θ,z′)=∑iσirect(ˆt−ΔtTp)e−j4πγc(ˆt−2Rrefc)RΔi⋅e−j4πcfcRΔiej4πγc2R2Δi (6) 式(6)即为实际采集的回波信号，式中相位项第1项为差频之后的单频项，对应频率为
fi= −γ(2RΔi/c) ,RΔi 的变化会使对应的距离项中的频率即fi 发生改变，对其沿着快时间做傅里叶变换，即可将不同距离处的目标分离出来，即脉冲压缩得到1维距离像，其数学表达为：S(fi,θ,z′)=∑iATpsinc[Tp(fi+2γcRΔi)]⋅e−j4πcfcRΔiej4πγc2R2Δi (7) 式(6)第2项的相位变化使回波产生多普勒，对应同一个目标不同脉冲时刻的相位变化。第3项为解线频调带来的剩余视频相位(Residual Video Phase, RVP)，可以通过相位补偿的方式将其补偿掉。而在本系统中由于在量级上
c≫RΔi ，所以此项也可以忽略不计。式(6)为回波数据的时域表达形式，如果用波数域的表达，同时补偿掉RVP项，其数学表达式为：S(Kr,θ,z′)=∑iσirect(ΔKr4πγTp)e−jΔKrRΔie−jKrcRΔi=∑iσirect(ΔKr4πγTp)e−jKrRΔi (8) 式中，
ΔKr=4πγˆt/c∈[−2πγTp/c,2πγTp/c] ,Krc= 4πfc/c ,Kr=Krc+ΔKr 为频率在空间的表达，即波数表达方式，为了得到回波信号关于距离-高度两个维度的2维频谱，对式(8)沿着Z轴做傅里叶变换到竖直方向的波数域可得：Sr(Kr,θ,Kz′)=∫S(Kr,θ,z′)e−jKz′z′dz′=∫∑iσi rect(ΔKr4πγTp)⋅e−jKrRΔie−jKz′z′dz′ (9) 为了下述表示方便，用
Rxyz=RΔi+Rref 表示在3维坐标系下的目标点到雷达天线的斜距，用Rxy=√(Rcosθ−x)2+(Rsinθ−y)2 表示竖直方向某一天线阵元到该阵元所在XOY平面上目标点的斜距，则，Rxyz=√R2xy+(z′−zi)2 对式(9)乘以e−jKrRref 来补偿由于去斜参考时间延迟带来的相位偏移。采用驻定相位原理导出式(9)的频谱解析式：Sr(Kr,θ,Kz′)=∑iσirect(ΔKr4πγTp)⋅e−j√K2r−K2z′Rxye−jKz′z (10) 式(10)的相位如下：
ϕsq(Kr,θ,Kz′)=−√K2r−K2z′Rxy−Kz′z (11) 式(11)包含两项，第2项线性相位项代表目标高度位置的响应；第1项是距离和Z方向耦合产生的RCM(Range Cell Migration)，代表目标的倾斜距离对RCM空间的变化。使用ω-K算法中的变量代换对式(11)中的第1项耦合相位解耦合：
√K2r−K2z′→Krc+Kxy (12) 式(12)叫做Stolt映射(Stolt Mapping, SM)，通常采用插值来实现。由于SM移除了距离方向与高度方向的耦合，同时使得在波数域
Kz′ 与Kxy 垂直[13]。经过SM，式(10)变为式(13)形式：Sr(Kr,θ,Kz′)=∑iσirect(ΔKr4πγTp)⋅e−jKxyRxye−jKz′ze−jKrcRxy (13) 式(13)通过距离IFT(Inverse Fourier Transform)处理可以将不同距离点目标在不同的距离单元区分开来，得到距离脉冲压缩的结果，再经过高度维IFT将信号变回慢时间域，实现高度维聚焦。对式(13)分别沿着距离和高度方向进行IFT，同时考虑到高度方向的窗函数，得到距离高度2维脉冲压缩结果具有式(14)形式：
\begin{align} {S_{\rm r}}\left( {\widehat t,\theta ,z'} \right) =& \sum\limits_i {{A_i}\operatorname{sin\,c} \left( {\widehat t - \tau } \right)\operatorname}\\ & \cdot {\sin {\rm c}} \left( {z' - {z_i}} \right){{\rm{e}}^{ - {\rm{j}}{K_{{\rm{rc}}}}{R_{{\rm{xy}}}}}} \end{align} (14) 式(14)中
Ai=Aσi , A表示信号处理过程中幅度系数。τ=2Rxy/c ，表征了目标散射点距离雷达相位中心的距离。通过上述ω-K算法，实现了距离与高度维的解耦合，此时可以将3维聚焦问题转化为在竖直方向的聚焦问题和在XOY平面内的投影问题。具体做法是：将3维成像区域沿竖直方向划分为N个XOY切片堆叠的形式，如图2所示。此时只需要对每个阵元对应的XOY切片都执行相同的聚焦操作，即可得到3维的聚焦图像，因此问题转化为对每个切片的聚焦问题。由于天线阵列的运动轨迹为圆周，因此选取对曲线运动适应性较好的BP算法来实现每个切片的2维SAR图像。BP算法的基本思想是先通过距离脉冲压缩，使不同距离的目标点分离开来，然后在方位向，通过对回波信号补偿相位进行相干积累，从而将处于同一距离的不同散射点区分开来，即实现方位向的聚焦。BP算法的具体实现方式如下，由于通过ω-K算法已经实现了竖直方向与距离方向的解耦合，因此对于每一个垂直于Z轴方向的2维XOY切片来说，式(14)可以简写为：
Sr(ˆt,θ)=∑iAisinc(ˆt−τ)e−jKrcRxy (15) 采用函数
f(x,y) 来表示2维XOY切片上像素点位置(x,y) 的散射函数，由时域相关原理可知，位于(xi,yj) 处点目标的散射函数可以由式(16)重建：f(xi,yj)=∫Sr(ˆt,θ)ejKrcRij(θ)dθ=∑θSr(ˆt,θ)ej4πfccRij(θ) (16) 式(16)说明了位于
(xi,yj) 处点目标的散射函数的重建过程，对于解线频调信号由式(7)和式(14)可知，后向投影之前要补充相位因子才能重建散射点的sinc响应并完成相干积累。Rij(θ) 即是成像网格点(xi,yj) 到雷达天线阵元的斜距，Rij(θ)= √(Rcosθ−xi)2+(Rsinθ−yj)2 ，该值是天线阵元位置即θ 的函数，对于每一个XOY切片，运用上述BP算法，即可重建该切片的2维SAR图像，然后将各切片沿Z轴方向堆叠，可以重建3维SAR图像。最后，为了便于对场景目标观察和检测，对重建的3维SAR图像按照不同的观测角度进行2维投影，即可得到全方位观测角度下的2维图像，完整的算法流程如图3所示。4. 算法仿真与验证
4.1 多角度成像仿真
按照上述算法流程，使用Matlab对3维人体模型进行回波仿真并重构3维图像，以此验证算法的可行性。仿真时，天线阵列以固定的角速度沿圆周运动，设圆周半径为1 m，扫描角度范围为
θ∈[π/6,5π/6] 。仿真所用3维人体模型及运动天线阵列3维轨迹如图4所示。
图 4 3维人体模型及阵列运动轨迹Figure 4. Three-dimensional human model and the trajectory of antenna array成像时，对所有回波数据进行ω-K聚焦之后，BP投影采用了数据复用的方式实现多个视角下的成像，图5为使用仿真数据，采用文中所述重构算法得到的3维人体模型仿真结果。
由图5可知，重构算法较好地完成了3维场景的重建，3维人体模型的轮廓清晰，直观上看，图像聚焦效果较好。并且不同角度下的3维人体模型重构图像质量没有损失，说明本文算法可以实现3维目标的多角度成像，完成了算法理论上可行性的验证。
4.2 非理想圆柱仿真
为了验证本文算法对天线阵列非理想圆柱扫描轨迹的适应性，仿真了一组阵列沿椭圆轨迹运动的3维人体回波数据，椭圆长轴设为1 m，短轴设为0.8 m。分别采用3维ω-K算法和本文算法进行对比，成像结果如图6所示。
由图6可知，在天线阵列圆周轨迹运动条件下，本文算法和3维ω-K算法均能较好地完成3维成像；但在天线阵列椭圆轨迹运动条件下，3维ω-K算法由于3维匹配滤波以及波数均匀化不准确等原因导致图像散焦和镜像，而本文算法依然能够完成精确聚焦，说明本文算法在非理想圆柱扫描模式下，依然能够较好地完成人体3维成像。
4.3 实测数据分析
为了验证本文算法在实际安检成像系统的工程可用性，对采集的人体回波数据进行了合成孔径成像处理，并分析了本文算法与3维频域算法成像结果的区别。所选实际安检成像系统工作在Ka波段，系统带宽为6 GHz, X, Y, Z 3个维度的理论分辨率分别为5 mm, 25 mm, 5 mm。圆周轨迹运动，扫描半径为0.8 m，脉冲重复频率为200 kHz，阵列圆迹运动采用先加速后减速方式，因此在等时采样条件下，角度维的采样间隔是非均匀的。成像结果对比如图7所示。
图 7 圆柱扫描实测数据成像结果对比Figure 7. Comparison of imaging results of measured data from cylindrical scanning从两种算法的对比结果可知，两种算法均能够实现图像聚焦，人体口袋里隐藏的刀具和枪支模型等都可以清晰呈现。相比于3维频域算法，由于BP算法不需要在角度维进行插值操作，也不存在斜距上的近似，因此保留了更多的细节，成像之后场景更为丰富。同一观测角度下，图像中人体腿部的散焦问题也得到明显改善。
从运算量上分析，不同于3维ω-K算法的运算量主要集中于FFT，本文算法的运算量主要集中于BP投影。若假设X, Y, Z 3维成像区域像素点数分别为Nx, Ny, Nz，距离、竖直、角度3个维度的采样点数分别为Nr, Ne, Na，则3维BP算法的运算复杂度为
Nx×Ny×Nz×Ne×Na 量级。本文算法BP投影部分运算复杂度为Nx×Ny×Nz×Na 量级。由于算法具备很好的并行结构设计，实测数据处理表明(3维成像网格大小为380×200×64 )，在NVIDIA GTX1070平台上采用CUDA(Compute Unified Device Architecture)编程成像计算时间约为0.18 s，满足实时成像的需求。5. 结语
主动式毫米波3维成像系统在安检中具有对人体安全，穿透性好，成像效果佳等优势，因此具备广泛的应用前景，也是目前3维成像系统的发展方向。本文结合主动式圆柱扫描毫米波3维成像系统的工作方式，提出了一种ω-K算法和BP算法相结合的3维场景重建算法，给出了其理论依据。对3维人体模型进行仿真验证，结果表明，本文算法能够完成3维场景的聚焦成像，具备理论上的可行性。同时对比了本文算法与3维ω-K算法在非理想圆柱条件下的成像效果，验证了本文算法对系统运动方式的普适性。然后进一步在实际的3维成像系统采集的回波数据进行了验证，确保了算法在工程应用中的可行性，为该体制下的3维成像系统提供了一种新的简单易行的场景重构算法。此外，将本文提到的算法进行了CUDA工程实现，实现了3维人体成像的实时处理。
-
图 4 不同阵元重叠数时杂波空时谱示意图
Figure 4. Space-time clutter spectrum diagram with different array overlapping numbers
图 5 不同阵元重叠数下的杂波谱分布
Figure 5. Space-time clutter spectrum with different array overlapping numbers
图 10 有无幅相误差情形下输出SCNR与输入CNR曲线
Figure 10. Output SCNR versus input CNR with or without amplitude and phase error
表 1 雷达系统仿真参数
Table 1. Simulation parameters of radar system
参数 数值 参数 数值 阵元总数 52 子阵数 10 子阵内阵元数 16 重叠阵元数 12 脉冲数 16 工作载频 1.2 GHz 波长 0.25 m 脉冲重频 2 kHz 信号带宽 5 MHz 平台高度 10 km 平台速度 125 m/s 目标方位角 82.79° 目标俯仰角 4.99° 目标斜距 115 km 空间频率展宽 0.1 多普勒滤波器带宽 0.06 幅度误差 0.03 相位误差 3° 信噪比 –5 dB 杂噪比 50 dB 
下载: 导出CSV
-
[1] BRENNAN L, MALLETT J, and REED I. Adaptive arrays in airborne MTI radar[J]. IEEE Transactions on Antennas and Propagation, 1976, 24(5): 607–615. doi: 10.1109/TAP.1976.1141412. [2] 谢文冲, 王永良, 熊元燚. 机载雷达空时自适应处理[M]. 北京: 清华大学出版社, 2024: 86–95.XIE Wenchong, WANG Yongliang, and XIONG Yuanyi. Space Time Adaptive Processing Technique for Airborne Radar[M]. Beijing: Tsinghua University Press, 2024: 86–95. [3] XU Jingwei, LIAO Guisheng, and SO H C. Space-time adaptive processing with vertical frequency diverse array for range-ambiguous clutter suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(9): 5352–5364. doi: 10.1109/TGRS.2016.2561308. [4] WEN Cai, HUANG Yan, PENG Jinye, et al. Slow-time FDA-MIMO technique with application to STAP radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(1): 74–95. doi: 10.1109/TAES.2021.3098100. [5] SONG Di, FENG Qi, CHEN Shengyao, et al. Space-time adaptive processing using deep neural network-based shrinkage algorithm under small training samples[J]. IEEE Transactions on Aerospace and Electronic Systems, 2023, 59(6): 9697–9703. doi: 10.1109/TAES.2023.3255840. [6] 阚庆云, 许京伟, 廖桂生. 和差天线空时自适应测角方法及性能分析[J]. 电子学报, 2023, 51(1): 42–49. doi: 10.12263/DZXB. 20211458.KAN Qingyun, XU Jingwei, and LIAO Guisheng. Angle estimation approach and performance analysis for STAP with sum-difference antenna[J]. Acta Electronica Sinica, 2023, 51(1): 42–49. doi: 10.12263/DZXB.20211458. [7] XIE Lei, HE Zishu, TONG Jun, et al. A recursive angle-Doppler channel selection method for reduced-dimension space-time adaptive processing[J]. IEEE Transactions on Aerospace and Electronic Systems, 2020, 56(5): 3985–4000. doi: 10.1109/TAES.2020.2983533. [8] HUANG Penghui, ZOU Zihao, XIA Xianggen, et al. A novel dimension-reduced space-time adaptive processing algorithm for spaceborne multichannel surveillance radar systems based on spatial-temporal 2-D sliding window[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 5109721. doi: 10.1109/TGRS.2022.3144668. [9] GOLDSTEIN J S and REED I S. Reduced-rank adaptive filtering[J]. IEEE Transactions on Signal Processing, 1997, 45(2): 492–496. doi: 10.1109/78.554317. [10] HERD J S and CONWAY M D. The evolution to modern phased array architectures[J]. Proceedings of the IEEE, 2016, 104(3): 519–529. doi: 10.1109/JPROC.2015.2494879. [11] HE Xiongpeng, LIAO Guisheng, ZHU Shengqi, et al. Range ambiguous clutter suppression approach with elevation time diverse array[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(1): 359–373. doi: 10.1109/TAES.2021.3101786. [12] ZHANG Tianfu, WANG Zhihao, XING Mengdao, et al. Research on multi-domain dimensionality reduction joint adaptive processing method for range ambiguous clutter of FDA-Phase-MIMO space-based early warning radar[J]. Remote Sensing, 2022, 14(21): 5536. doi: 10.3390/rs14215536. [13] LI Yuanshuai, CHANG Shaoqiang, LIU Zihao, et al. Range ambiguity suppression under high-resolution estimation using the MUSIC-AP algorithm for pulse-Doppler radar[J]. Signal Processing, 2024, 214: 109237. doi: 10.1016/j.sigpro.2023.109237. [14] WU Yong, TANG Jun, and PENG Yingning. On the essence of knowledge-aided clutter covariance estimate and its convergence[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(1): 569–585. doi: 10.1109/TAES.2011.5705692. [15] WANG Degen, WANG Tong, CUI Wenchen, et al. A clutter suppression algorithm via enhanced sparse Bayesian learning for airborne radar[J]. IEEE Sensors Journal, 2023, 23(10): 10900–10911. doi: 10.1109/JSEN.2023.3263919. [16] VARADARAJAN V and KROLIK J L. Joint space-time interpolation for distorted linear and bistatic array geometries[J]. IEEE Transactions on Signal Processing, 2006, 54(3): 848–860. doi: 10.1109/TSP.2005.862941. [17] WANG Yongliang, DUAN Keqing, and XIE Wenchong. Cross beam STAP for nonstationary clutter suppression in airborne radar[J]. International Journal of Antennas and Propagation, 2013, 2013: 276310. doi: 10.1155/2013/276310. [18] MAILLOUX R J, SANTARELLI S G, ROBERTS T M, et al. Irregular polyomino-shaped subarrays for space-based active arrays[J]. International Journal of Antennas and Propagation, 2009, 2009: 956524. doi: 10.1155/2009/956524. [19] 黄飞, 盛卫星, 马晓峰. 随机错位子阵阵列天线及其优化设计[J]. 电波科学学报, 2008, 23(5): 917–921. doi: 10.3969/j.issn.1005-0388.2008.05.023.HUANG Fei, SHENG Weixing, and MA Xiaofeng. Plane antenna arrays with randomly staggered subarrays and its optimal design[J]. Chinese Journal of Radio Science, 2008, 23(5): 917–921. doi: 10.3969/j.issn.1005-0388.2008.05.023. [20] XIONG Ziyuan, XU Zhenhai, CHEN Siwei, et al. Subarray partition in array antenna based on the algorithm X[J]. IEEE Antennas and Wireless Propagation Letters, 2013, 12: 906–909. doi: 10.1109/LAWP.2013.2272793. [21] KARADEMIR S, PROKOPYEV O A, and MAILLOUX R J. Irregular polyomino tiling via integer programming with application in phased array antenna design[J]. Journal of Global Optimization, 2016, 65(2): 137–173. doi: 10.1007/s10898-015-0354-8. [22] DONG Wei, XU Zhenhai, LIU Xinghua, et al. Irregular subarray tiling via heuristic iterative convex relaxation programming[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(4): 2842–2852. doi: 10.1109/TAP.2019.2955070. [23] ROCCA P, ANSELMI N, POLO A, et al. An irregular two-sizes square tiling method for the design of isophoric phased arrays[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(6): 4437–4449. doi: 10.1109/TAP.2020.2970088. [24] ZHOU Jian, WANG Yiqing, WANG Zhangling, et al. Irregular subarray tiling via rotational symmetry[J]. IEEE Antennas and Wireless Propagation Letters, 2023, 22(4): 903–907. doi: 10.1109/LAWP.2022.3228108. [25] MA Yankai, YANG Shiwen, CHEN Yikai, et al. High-directivity optimization technique for irregular arrays combined with maximum entropy model[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(7): 3913–3923. doi: 10.1109/TAP.2020.3044653. [26] YANG Feng, MA Yankai, LONG Weijun, et al. Synthesis of irregular phased arrays subject to constraint on directivity via convex optimization[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(7): 4235–4240. doi: 10.1109/TAP.2020.3044632. [27] JIANG Hailing, GONG Yu, ZHANG Junyi, et al. Irregular modular subarrayed phased array tiling by algorithm X and differential evolution algorithm[J]. IEEE Antennas and Wireless Propagation Letters, 2023, 22(7): 1532–1536. doi: 10.1109/LAWP.2023.3250260. [28] TOYAMA N. Optimization of aperiodic array patterns using GA and SDA[C]. 2005 IEEE Antennas and Propagation Society International Symposium, Washington, USA, 2005: 69–72. doi: 10.1109/APS.2005.1551737. [29] WANG Hao, FANG Dagang, and CHOW Y L. Grating lobe reduction in a phased array of limited scanning[J]. IEEE Transactions on Antennas and Propagation, 2008, 56(6): 1581–1586. doi: 10.1109/TAP.2008.923354. [30] 程乃平, 潘点飞. 大型阵列天线子阵划分及栅瓣抑制方法[J]. 信号处理, 2014, 30(5): 535–543. doi: 10.3969/j.issn.1003-0530.2014.05.007.CHENG Naiping and PAN Dianfei. Subarray partition method and grating lobe suppression for large array antenna[J]. Journal of Signal Processing, 2014, 30(5): 535–543. doi: 10.3969/j.issn.1003-0530.2014.05.007. [31] KRIVOSHEEV Y V, SHISHLOV A V, and DENISENKO V V. Grating lobe suppression in aperiodic phased array antennas composed of periodic subarrays with large element spacing[J]. IEEE Antennas and propagation Magazine, 2015, 57(1): 76–85. doi: 10.1109/MAP.2015.2397155. [32] XU Xiaomin, LIAO Cheng, ZHOU Liang, et al. Grating lobe suppression of non-uniform arrays based on position gradient and sigmoid function[J]. IEEE Access, 2019, 7: 106407–106416. doi: 10.1109/ACCESS.2019.2932123. [33] YANG Gongqing, ZENG Hui, and XU Zhenhai. Adaptive gradient search algorithm for displaced subarrays with large element spacing[J]. IEEE Antennas and Wireless Propagation Letters, 2021, 20(7): 1155–1159. doi: 10.1109/LAWP.2021.3074253. [34] ELAYAPERUMAL S and RAJAGOPAL S. Design of antenna array architecture with large inter element spacing and low grating lobes[C]. 2021 IEEE International Conference on Electronics, Computing and Communication Technologies, Bangalore, India, 2021: 1–6. doi: 10.1109/CONECCT52877.2021.9622535. [35] LIU Zhixin, ZHU Shengqi, XU Jingwei, et al. Range-ambiguous clutter suppression for STAP-based radar with vertical coherent frequency diverse array[J]. IEEE Transactions on Geoscience and Remote Sensing, 2023, 61: 5106517. doi: 10.1109/TGRS.2023.3291738. [36] GUERCI J R. Theory and application of covariance matrix tapers for robust adaptive beamforming[J]. IEEE Transactions on Signal Processing, 1999, 47(4): 977–985. doi: 10.1109/78.752596. 期刊类型引用(5)
1. 杨磊，霍鑫，申瑞阳，宋昊，胡仲伟. 可信推断近场稀疏综合阵列三维毫米波成像. 雷达学报. 2024(05): 1092-1108 . 
本站查看2. 杨磊，宋昊，申瑞阳，陈英杰，胡仲伟，霍鑫，邢孟道. 强稀疏低副瓣近场聚焦稀疏阵列三维成像. 电子与信息学报. 2024(12): 4471-4482 . 
百度学术3. 杨磊，王腾腾，陈英杰，盖明慧，许瀚文. 低秩矩阵补全高分辨SAR成像特征重建. 电子与信息学报. 2023(08): 2965-2974 . 百度学术4. 张建新，黄平平，孙莹. 透视万物的力量——解析毫米波技术在高安保级别领域的融合应用. 中国安防. 2022(04): 30-35 . 百度学术5. 陈曦阳，王鑫，魏超，邓仕发，强勇，骆睿，吴光胜. 主动式毫米波人体成像安检仪关键实用技术解析. 中国安全防范技术与应用. 2022(Z2): 23-30 . 百度学术其他类型引用(9)
-
下载:
百度学术
计量
- 文章访问数: 426
- HTML全文浏览量: 109
- PDF下载量: 233
- 被引次数: 14
