
Citation: | WANG Xu, HUANG Dongping, WEI Guohua, BAI Jiahao, ZHU Qinyuan. Sparse Uniform Array Grating Lobe Suppression Using Dual-carrier Frequency Pattern Multiplication[J]. Journal of Electronics & Information Technology, 2023, 45(2): 534-541. doi: 10.11999/JEIT211492 |
稀疏阵列具有布阵灵活经济和角分辨率高的优势,在雷达探测、成像等领域得到了广泛的应用。对稀疏阵列的研究一直是国内外学者的研究热点[1,2]。稀疏阵列天线稀疏分布时阵列方向图将出现大量栅瓣[3,4],这些栅瓣可能会导致在角度估计中出现模糊的情况,从而严重影响参数估计性能。因而对稀疏阵列栅瓣抑制方法的研究具有重要的意义。
针对天线间距达到数十倍波长稀疏阵列的栅瓣抑制问题,通过优化阵列中各天线的位置抑制方向图栅瓣的传统方法[5-10]会因为优化的自由度过低和搜索量过大而不再适用。文献[11]提出与自适应变化的填充因子机制相结合使用的迭代 Chirp-Z 变换技术,在线性阵列细化的快速低旁瓣方向图综合中可以获得十分优秀的旁瓣抑制性能。该算法避免了陷入局部最优问题的出现,加快了合成效率并促进了全局收敛,具有良好的鲁棒性。但针对天线数量以及位置确定的阵列该方法将不能起到很好的旁瓣抑制效果。为此,文献[12]在确定圆形阵列下通过相位模式分析计算天线的激励幅度和相位进行方向图综合可以有效降低旁瓣电平,但该方法仅适用稀疏度小阵列的旁瓣抑制。同时,有研究者提出了将多个不同稀疏阵列的方向图进行融合处理的栅瓣抑制方法。其中,文献[13]通过调整3个均匀稀疏线阵天线之间的间距,使得各阵列的栅瓣出现位置与零点出现位置接近,然后把3个不同稀疏阵列的方向图进行融合处理,从而达到抑制稀疏阵列栅瓣的目标。该方法对栅瓣虽有一定的抑制作用,但仍需要对子阵的间距进行优化,而且栅瓣位置随阵列位置变化过于随机,不利于选取合适的阵列组合。文献[14]在分布式星载雷达上利用频分多输入多输出(Multiple-Input Multiple-Output, MIMO)技术进行子阵方向图融合,达到了抑制稀疏阵列栅瓣的效果。文献[15]设计了由密集发射阵列与稀疏接收阵列构成的MIMO阵列,通过使用Kaiser窗对发射阵列进行加权,调整发射方向图的零点出现的位置,使其与接收方向图栅瓣重合,最后将发射方向图与接收方向图进行融合处理从而达到抑制阵列方向图栅瓣的目的。上述方法虽对稀疏阵列栅瓣能起到了抑制的效果,但对其研究均停留在定性分析与数值仿真层面,缺乏对影响算法性能关键参数的定量分析。
本文提出一种基于双载频方向图相乘的栅瓣抑制算法,可适用于天线间距达到数十倍波长的情况,在最优频率差下对稀疏阵列栅瓣能达到25 dB的抑制效果。多目标探测可以根据距离和速度等维度进行目标分辨并确定目标数量,因而可以将其分解为多个单目标探测任务。当稀疏阵列的栅瓣得到有效抑制后,单目标测角模糊的问题便得到了解决,随即多目标测角模糊的问题也得到了解决,因此本文重点在单目标情况下进行算法讨论。文中分析了影响该栅瓣抑制算法性能的因素,然后以峰值旁瓣比作为栅瓣抑制性能的评价指标,对影响该算法性能的频率差关键参数进行定量分析,推导出栅瓣抑制后峰值旁瓣比与频率差以及载频间的最大差异与阵列天线间距关系的表达式。
本文第2节介绍文中所使用的阵列和信号模型。第3节介绍本文所提栅瓣抑制算法的原理以及该算法对应信号处理的实现流程。第4节定量分析影响该算法栅瓣抑制性能的频率差因素,并给出峰值旁瓣比与频率差的关系表达式。第5节给出不同条件下的算法仿真。第6节陈述相应的结论。
本文考虑由M个天线组成的发射阵列和N个天线组成的接收阵列构成的均匀稀疏分布的1维分布式MIMO阵列。发射阵列天线之间的间距
假定发射天线将基带信号
z(t) = [z0(t),z1(t),⋯,zM−1(t)]T | (1) |
其中,
理想情况下接收天线的接收回波信号矢量
x(t)=αexp(2πf0t)ar0(θ)aTt0(θ)s(t)+αexp(2πf1t)ar1(θ)aTt1(θ)s(t) + n(t) | (2) |
其中,
ar0(θ) = [1,exp(−j2ρ0dRsinθ),⋯,exp(−j2(N−1)ρ0dRsinθ)]T | (3) |
at0(θ) = [1,exp(−j2ρ0dTsinθ),⋯,exp(−j2(M−1)ρ0dTsinθ)]T | (4) |
ar1(θ) = [1,exp(−j2ρ1dRsinθ),⋯,exp(−j2(N−1)ρ1dRsinθ)]T | (5) |
at1(θ) = [1,exp(−j2ρ1dTsinθ),⋯,exp(−j2(M−1)ρ1dTsinθ)]T | (6) |
其中,
将
x0(t)=[x(t)exp(−2πf0t)]∗h(t)=αar0(θ)aTt0(θ)s(t) + n0(t) | (7) |
x1(t)=[x(t)exp(−2πf1t)]∗h(t)=αar1(θ)aTt1(θ)s(t) + n1(t) | (8) |
其中,
载频1,2对应的接收向量经过MN个匹配滤波器滤波后的输出矢量
v0(t)=[x0(t)⊗IM]∗[IN⊗s∗(−t)]=αβ(t)a0(θ) + e0(t) | (9) |
v1(t)=[x1(t)⊗IM]∗[IN⊗s∗(−t)]=αβ(t)a1(θ) + e1(t) | (10) |
其中,
对载频1,2对应信号匹配滤波后结果
y0(t)=wH0v0(t) | (11) |
y1(t)=wH1v1(t) | (12) |
其中,
假设期望信号的来波方向为
F0(θ)=sin[MdTρ0(sinθ−sinθ0)]Msin[dTρ0(sinθ−sinθ0)]⋅sin[NdRρ0(sinθ−sinθ0)]Nsin[dRρ0(sinθ−sinθ0)] | (13) |
F1(θ)=sin[MdTρ1(sinθ−sinθ0)]Msin[dTρ1(sinθ−sinθ0)]⋅sin[NdRρ1(sinθ−sinθ0)]Nsin[dRρ1(sinθ−sinθ0)] | (14) |
分析式(13)和式(14)可知,载频1, 2对应的阵列方向图栅瓣出现的位置
θGL0=arcsin(iγ+sinθ0),i=±1,2,⋯ | (15) |
θGL1=arcsin(f0f1⋅jγ+sinθ0),j=±1,2,⋯ | (16) |
其中,
由式(15)和式(16)可知,阵列方向图栅瓣出现的位置与信号的载频有关。当阵列的主瓣位置指向一定,两组不同载频信号对应的阵列方向图具备主瓣位置相同但栅瓣位置存在差异的特性,利用该特性将两组不同载频对应的阵列方向图进行相乘处理。两组不同载频对应方向图相乘得到的归一化方向图
F(θ)=F0(θ)⋅F1(θ) | (17) |
为了便于后续的讨论,可将发射与接收天线间距系数表示为
由归一化阵列方向图特性可知,对
F0(θ)≤1F1(θ)≤1} | (18) |
然后结合式(17)和式(18)可知,对
F(θ)≤G(θ) | (19) |
其中,
基于上述原理提出了基于双载频方向图相乘的稀疏阵列栅瓣抑制算法,该算法在阵列信号处理过程中具体实现流程如图2所示,该算法的详细实现步骤如算法1所示。
算法1 栅瓣抑制算法 |
发射与接收端: |
(1) 发射双载频信号z(t); |
(2) 接收回波信号x(t)。 |
信号处理端: |
(1) 利用式(7)和式(8)对x(t)进行下变频和频域滤波分离出x0(t) 和x1(t)两路信号; |
(2) 利用式(9)和式(10)分别对x0(t)和x1(t)进行匹配滤波得到 v0(t)和v1(t); |
(3) 使用式(11)和式(12)对v0(t)和v1(t)信号分别进行波束形成 得到y0(t)和y1(t); |
(4) 对y0(t)和y1(t)进行归一化得到˜y0(t)和˜y1(t); |
(5) 将˜y0(t),˜y1(t)进行相乘处理到˜y(t); |
(6) 从˜y(t)选取检测单元,从˜y0(t)或˜y1(t)选取该检测单元对应 的参考单元进行恒虚警检测; |
(7) 调整波束指向转到步骤(3),完成空域的搜索,退出。 |
为评价栅瓣抑制算法对栅瓣的抑制性能,本文使用峰值旁瓣比(Peak SideLobe Ratio, PSLR)作为评价指标去衡量对栅瓣的抑制效果,其定义为
PSLR = 20lgUMSLUML | (20) |
其中,
为进一步讨论该算法对栅瓣抑制的性能,给出以下公式推导证明。均匀稀疏分布的1维分布式MIMO阵列两组不同载频方向图相乘后得到的方向图
F(θ)=|sin[MdTℜ0(θ)]Msin[dTℜ0(θ)]⋅sin[NdRℜ0(θ)]Nsin[dRℜ0(θ)]⋅sin[MdTℜ1(θ)]Msin[dTℜ1(θ)]⋅sin[NdRℜ1(θ)]Nsin[dRℜ1(θ)]| | (21) |
ℜ0(θ) = ρ0(sinθ−sinθ0) | (22) |
ℜ1(θ) = ρ1(sinθ−sinθ0) | (23) |
令
g(x)=20lg{|sin[MdTℜ0(x)]Msin[dTℜ0(x)]⋅sin[NdRℜ0(x)]Nsin[dRℜ0(x)]⋅sin[MdTℜ1(x)]Msin[dTℜ1(x)]⋅sin[NdRℜ1(x)]Nsin[dRℜ1(x)]|} | (24) |
当信号频率差在
N = {x|f0/f0qξf1qξf1+x0<x<1/1qξqξ+x0} | (25) |
其中,
那么求解
xmax=argmaxxg(x),s.t.x∈N | (26) |
令
g0(x)=g1(x),s.t.x∈N | (27) |
(1) 对称性
g0(2qξ+2x0−x)=20lg{|sin[Mπqa(2qξ+2x0−x−x0)]Msin[πqa(2qξ+2x0−x−x0)]⋅sin[Nπqb(2qξ+2x0−x−x0)]Nsin[πqb(2qξ+2x0−x−x0)]|}20lg{|sin[2Mπaξ−Mπqa(x−x0)]Msin[2πaξ−πqa(x−x0)]⋅sin[2Nπbξ−Nπqb(x−x0)]Nsin[2πbξ−πqb(x−x0)]|} | (28) |
因为
g0(2qξ+2x0−x)=20lg{|sin[Mπqa(x−x0)]Msin[πqa(x−x0)]⋅sin[Nπqb(x−x0)]Nsin[πqb(x−x0)]|}=g0(x) | (29) |
由此可知,
(2) 平移不变性
g1(x−Δx)=20lg{|sin[MdTρ1(x−x0−Δx)]Msin[dTρ1(x−x0−Δx)]sin[NdRρ1(x−x0−Δx)]Nsin[dRρ1(x−x0−Δx)]|}20lg10{|sin[Mπλ0dT(x−x0)−MπΔfcdT(x−x0−1qξ)]Msin[πλ0dT(x−x0)−πΔfcdT(x−x0−1qξ)]⋅sin[Nπλ0dR(x−x0)−NπΔfcdR(x−x0−1qξ)]Nsin[πλ0dR(x−x0)−πΔfcdR(x−x0−1qξ)]|} | (30) |
其中,
综上所述,当
PSLR=40lg{|sin[(2+β)Mπa2(1+β)ξ]sin[(2+β)Nπb2(1+β)ξ]MNsin[(2+β)πa2(1+β)ξ]sin[(2+β)πb2(1+β)ξ]|},0<β<1q(χ + 2) | (31) |
其中,频率差系数
为了验证本文所介绍的栅瓣抑制算法原理的正确性以及栅瓣抑制性能。本文设计了由10个发射天线与10个接收天线组成的均匀稀疏分布的1维分布式MIMO阵列的仿真实验。具体的仿真参数如表1所示。
发射天线 数量M | 接收天线 数量N | 发射天线 间距系数γT | 接收天线 间距系数γR | 天线间距 公约系数q | 基准信号 频率f0 |
10 | 10 | 20 | 25 | 5 | 10 GHz |
两组不同载频对应的阵列方向图相乘栅瓣抑制性能与两组信号的频率差有关。为了验证算法栅瓣抑制原理的正确性在表1所示的仿真条件基础上设计多组不同频率差的仿真实验。
在该仿真条件下最大频率差系数约为0.028,在该频率差系数范围内设计了频率差系数为0.013和0.026的两组仿真实验。频率差系数为0.013和天线波束最大值指向为0 °时两组不同载频对应的阵列方向图相乘得到的方向图及局部放大图如图3(a)所示;频率差系数为0.026和天线波束最大值指向为0 °时两组不同载频对应的阵列方向图相乘得到的方向图及局部放大图如图3(b)所示。频率差系数为0.013和天线波束最大值指向为20°时两组不同载频对应的阵列方向图相乘得到的方向图及局部放大图如图4(a)所示;频率差系数为0.026和天线波束最大值指向为20 °时两组不同载频对应的阵列方向图相乘得到的方向图及局部放大图如图4(b)所示。
由图3(a)和图4(a)可知,在频率差系数为0.013下使用该栅瓣抑制算法后,天线最大值指向为0°和20°方向图峰值旁瓣比分别–4.99 dB, –4.99 dB;由图3(b)和图4(b)可知,在频率差系数为0.026下使用该栅瓣抑制算法后,天线最大值指向为0°和20°方向图峰值旁瓣比分别为–22.08 dB, –22.08 dB;由此可知,天线波束进行扫描并不会影响算法的栅瓣抑制性能。
由图3(a)或图4(a)可知,当频率差系数为0.013时,使用该栅瓣抑制算法后方向图峰值旁瓣比为–4.99 dB;由图3(b)或图4(b)可知当频率差系数为0.026时,使用该栅瓣抑制算法后方向图峰值旁瓣比为–22.08 dB。由此可知,基于不同载频方向图相乘的栅瓣抑制算法可以有效抑制栅瓣,而且栅瓣的抑制性能与频率差有关。
由上节分析可知,该算法栅瓣抑制的性能与两组载频的频率差有关,在一定范围内频率差越大算法的栅瓣抑制性能越好。为验证本文所推导的该算法栅瓣抑制后峰值旁瓣比与频率差关系表达式的正确性,不失一般性地在表1所示仿真条件基础上设计了天线最大值指向为0°时多组不同频率差的仿真实验。
两组不同载频对应的方向图相乘后得到方向图峰值旁瓣对应角度与频率差的关系曲线(近似值与理论值),如图5所示。两组不同载频对应的阵列方向图相乘得到方向图的峰值旁瓣比与两组信号频率差的关系曲线(近似值与理论值),如图6所示。
由图5和图6可知,不同频率差下得到的近似峰值旁瓣比以及对应角度与理论峰值旁瓣比以及对应角度相吻合,由此验证该算法栅瓣抑制性能近似表达式的正确性。由图6可知,在一定频率差范围内峰值旁瓣比随频率差值增大而减小。因此,频率差越大算法栅瓣抑制性能越强。
为了抑制稀疏阵列因天线稀疏分布而产生的大量栅瓣,本文提出一种利用双载频方向图相乘的栅瓣抑制算法。该算法利用不同载频阵列方向图主瓣与栅瓣相对位置关系存在差异的特性实现对栅瓣的有效抑制。在最优频率差下,该算法对稀疏阵列栅瓣有25 dB的抑制效果,且存在无需进行大规模搜索的优势。文中推导得到的栅瓣抑制后峰值旁瓣比与频率差的关系表达式可以为雷达系统设计过程中快速选择最优频率差提供理论依据。
[1] |
DENG Weibo, WANG Hongyong, SUO Ying, et al. Grating lobe control of two-way pattern for linear distributed array[C]. 2020 IEEE International Symposium on Antennas and Propagation and North American Radio Science Meeting, Montreal, Canada, 2020: 621–622.
|
[2] |
陈志坤, 杜康, 彭冬亮, 等. 基于混合三角变异差分进化算法的平面稀疏阵列约束优化[J]. 电子与信息学报, 2020, 42(4): 895–901. doi: 10.11999/JEIT190705
CHEN Zhikun, DU Kang, PENG Dongliang, et al. Planar sparse array constraint optimization based on hybrid trigonometric mutation differential evolution algorithm[J]. Journal of Electronics &Information Technology, 2020, 42(4): 895–901. doi: 10.11999/JEIT190705
|
[3] |
WANG Zhe, WANG Wenqin, ZHENG Zhi, et al. Nested array sensor with grating lobe suppression and arbitrary transmit–receive beampattern synthesis[J]. IEEE Access, 2018, 6: 9227–9237. doi: 10.1109/ACCESS.2018.2804486
|
[4] |
ZHU Rongqiang, ZHOU Jianxiong, JIANG Ge, et al. Grating lobe suppression in near range MIMO array imaging using zero migration[J]. IEEE Transactions on Microwave Theory and Techniques, 2020, 68(1): 387–397. doi: 10.1109/TMTT.2019.2941188
|
[5] |
LU Baojun, SANG Lei, LI Xiangxiang, et al. Research on reducing grating lobe of adjoining regular subarray using genetic algorithms (GA)[C]. 2016 11th International Symposium on Antennas, Propagation and EM Theory (ISAPE), Guilin, China, 2016: 475–477.
|
[6] |
丛雯珊, 余岚, 沃江海. 基于粒子群算法的宽带真延时方向图栅瓣抑制方法[J]. 电子与信息学报, 2019, 41(7): 1698–1704. doi: 10.11999/JEIT180719
CONG Wenshan, YU Lan, and WO Jianghai. A grating lobe suppression method of wideband real time delay pattern based on particle swarm optimization algorithm[J]. Journal of Electronics &Information Technology, 2019, 41(7): 1698–1704. doi: 10.11999/JEIT180719
|
[7] |
RUDOLPH D L and BAROTT W C. Reduction of near-field grating lobes in sparse linear phased arrays[C]. 2014 IEEE Antennas and Propagation Society International Symposium (APSURSI), Memphis, USA, 2014: 1155–1156.
|
[8] |
LEMES D L, HECKLER M V T, GREDA L A, et al. Analysis of dual-band non-uniformly spaced arrays for mobile communications[C]. 2017 SBMO/IEEE MTT-S International Microwave and Optoelectronics Conference (IMOC), Aguas de Lindoia, Brazil, 2017: 1–5.
|
[9] |
CUI Can, LI Wentao, YE Xiutiao, et al. Hybrid genetic algorithm and modified iterative fourier transform algorithm for large thinned array synthesis[J]. IEEE Antennas and Wireless Propagation Letters, 2017, 16: 2150–2154. doi: 10.1109/LAWP.2017.2700865
|
[10] |
FU Yu, GUO Zhigui, WANG Haowen, et al. Optimization of planar thinned antenna array based on genetic and convex hybrid algorithm[C]. 2016 Progress in Electromagnetic Research Symposium (PIERS), Shanghai, China, 2016: 1299–1303.
|
[11] |
ZHOU Zhigang, ZENG Cao, and CHEN Baixiao. Fast low-sidelobe pattern synthesis for linear array thinning utilizing a modified iterative Chirp-Z transform technique[J]. IEEE Sensors Journal, 2021, 21(20): 23480–23491. doi: 10.1109/JSEN.2021.3108595
|
[12] |
SINGH M and GHOSH B. Antenna pattern synthesis using phase mode for circular array with reduced side-lobe level[C]. 2020 IEEE Asia-Pacific Microwave Conference (APMC), Hong Kong, China, 2020: 173–175.
|
[13] |
CUI Xiaowei, ZHANG Qingtai, and FENG Zhenming. Outage performance for maximal ratio combiner in the presence of unequal-power co-channel interferers[J]. IEEE Communications Letters, 2004, 8(5): 289–291. doi: 10.1109/LCOMM.2004.827401
|
[14] |
冯晶晶. 分布式星载雷达栅瓣抑制技术研究[D]. [硕士论文], 电子科技大学, 2010.
FENG Jingjing. Research on grating lobe suppression technology of distributed spaceborne radar[D]. [Master. dissertation], University of Electronic Science and Technology of China, 2010.
|
[15] |
ZHU Rongqiang, ZHOU Jianxiong, and FU Qiang. A grating and side lobe suppression method for near range MIMO array[C]. 2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Xiamen, China, 2019: 1–4.
|
[16] |
李道京, 侯颖妮, 滕秀敏, 等. 稀疏阵列天线雷达技术及其应用[M]. 北京: 科学出版社, 2014: 172–173.
|
1. | 刘海波,韩宇新,武兴恺,陈广茂. 圆柱共形阵低副瓣波束形成算法研究. 信号处理. 2023(12): 2190-2204 . ![]() |
算法1 栅瓣抑制算法 |
发射与接收端: |
(1) 发射双载频信号z(t); |
(2) 接收回波信号x(t)。 |
信号处理端: |
(1) 利用式(7)和式(8)对x(t)进行下变频和频域滤波分离出x0(t) 和x1(t)两路信号; |
(2) 利用式(9)和式(10)分别对x0(t)和x1(t)进行匹配滤波得到 v0(t)和v1(t); |
(3) 使用式(11)和式(12)对v0(t)和v1(t)信号分别进行波束形成 得到y0(t)和y1(t); |
(4) 对y0(t)和y1(t)进行归一化得到˜y0(t)和˜y1(t); |
(5) 将˜y0(t),˜y1(t)进行相乘处理到˜y(t); |
(6) 从˜y(t)选取检测单元,从˜y0(t)或˜y1(t)选取该检测单元对应 的参考单元进行恒虚警检测; |
(7) 调整波束指向转到步骤(3),完成空域的搜索,退出。 |
发射天线 数量M | 接收天线 数量N | 发射天线 间距系数γT | 接收天线 间距系数γR | 天线间距 公约系数q | 基准信号 频率f0 |
10 | 10 | 20 | 25 | 5 | 10 GHz |
算法1 栅瓣抑制算法 |
发射与接收端: |
(1) 发射双载频信号z(t); |
(2) 接收回波信号x(t)。 |
信号处理端: |
(1) 利用式(7)和式(8)对x(t)进行下变频和频域滤波分离出x0(t) 和x1(t)两路信号; |
(2) 利用式(9)和式(10)分别对x0(t)和x1(t)进行匹配滤波得到 v0(t)和v1(t); |
(3) 使用式(11)和式(12)对v0(t)和v1(t)信号分别进行波束形成 得到y0(t)和y1(t); |
(4) 对y0(t)和y1(t)进行归一化得到˜y0(t)和˜y1(t); |
(5) 将˜y0(t),˜y1(t)进行相乘处理到˜y(t); |
(6) 从˜y(t)选取检测单元,从˜y0(t)或˜y1(t)选取该检测单元对应 的参考单元进行恒虚警检测; |
(7) 调整波束指向转到步骤(3),完成空域的搜索,退出。 |
发射天线 数量M | 接收天线 数量N | 发射天线 间距系数γT | 接收天线 间距系数γR | 天线间距 公约系数q | 基准信号 频率f0 |
10 | 10 | 20 | 25 | 5 | 10 GHz |