-
摘要: 该文基于对称的均匀线阵提出了一种近场非圆信号参数的快速估计算法,算法基于信号的非圆特性以及阵列的对称性对近场导向矢量进行解耦,并利用多项式求根取代传统的谱峰搜索对近场源的角度及距离参数进行快速估计。基于给定的阵列结构,建立非圆信号参数估计的多项式数学模型,然后对其进行求根即可获得近场信号源位置参数。所提算法采用多项式求根的方法有效地降低了运算复杂度,同时利用信号的非圆特性提高了参数估计的自由度(DOF)。通过性能分析和计算机仿真实验可以看出该算法能够分辨更多的近场非圆信号,并且参数估计性能有所提升,更接近于近场源参数估计的克拉美罗界(CRB)。Abstract: In this paper, a fast parameter estimation algorithm for near-field non-circular signals is proposed, which is based on symmetric, uniform, and linear array, that decouple the near-field steering vector via non-circularity and symmetry. Then, the position parameters are estimated based upon polynomial rooting technique instead of traditional spectral search. After fixing the array structure, the parameters for estimation equation are ascertained. Then the values of position parameters are obtained by solving the polynomial equations; consequently, computational difficulty is alleviated effectively. In addition, the Degree-Of-Freedom (DOF) is increased by utilizing the non-circularity characteristics of the impinging sources. It can be seen from the performance analysis and computer simulation experiments that the proposed algorithm can resolve more near-field non-circular signals and can improve parameter estimation performance, which is closer to the Cramer–Rao Bound (CRB) estimation.
-
1. 引言
太赫兹波(THz waves)一般指电磁波频率从0.1 THz至10 THz之间的电磁波,具有光子能量低、穿透性好、高分辨等特性。太赫兹成像系统在无损探测和非接触式安检成像等领域有重大的应用价值[1]。现有的成像系统多采用机械扫描的方式实现2维波束扫描[2],这类天线体积大、加工难、成本高。此外,机械扫描体制使得成像时间长达数秒甚至数十秒,因此系统的实际应用价值大打折扣。为实现太赫兹高帧率成像,电控波束扫描体制是一种很有潜力的解决方案。通常采用基于移相器的相控阵天线实现电控波束扫描,但是在F波段以上频率,移相器很难实现180°以上的相移范围[3]。另一种实现电控波束扫描的方式是频率扫描天线[4]。本文提出将频率扫描天线应用于太赫兹成像系统,频率控制的波束扫描方式在太赫兹高帧率成像方面有巨大的应用潜力。
目前在高于0.1 THz的频段范围内,国内外关于频率扫描天线的研究逐渐增多[5–8],已有的研究中太赫兹波段频率扫描天线的实现方式主要有两种,一种是基于空间馈电的反射栅式频率扫描天线[5,6],这种天线通过将入射高斯波束转化成高阶衍射波模式,实现频率控制的波束扫描,该种天线具有高增益的优点,但是扫描角度较小。另一种实现方式是基于漏波体制的频率扫描天线,该种频率扫描天线可以克服扫描角度小的缺陷。文献[7]提出一种工作在212 GHz的介质波导加载周期金属栅条的漏波结构,介质导中传播的慢波主模,周期性金属栅条起到辐射电磁波的作用。文献[8]同样通过介质波导传播慢波,不同的是在介质波导侧壁周期性开槽实现漏波辐射,天线工作在97~103 GHz。基于复合左右手传输线理论的漏波天线可以实现法线两侧的波束扫描,目前该项技术可以实现X波段的大角度波束扫描[9];吴柯等人提出一种Ka波段的基于半模基片集成波导[10]的毫米波漏波天线,随后提出基于介质集成镜像波导的漏波天线[11],将工作频率提高到100 GHz[12]。西班牙学者提出通过控制电磁波在波导中传输的长度,实现了太赫兹频段的1维频控波束扫描[13],随后提出将波导结构在E面、H面同时控制传输路径,调节两个维度的相位差,实现了中心频率270 GHz的2维频率扫描天线[14]。由于电磁波在介质波导中传播不可避免地会引起损耗,基于介质波导的漏波天线效率比较低,采用金属波导传播太赫兹波可以达到提高效率的目的。
目前,波导缝隙天线已成功应用于从S波段到Ka波段的雷达系统中,波导缝隙天线通常采用矩形波导主模馈电,缝隙分布在波导宽边或者窄边[15],实现主模的漏波辐射。应用较多的波导缝隙天线多为谐振式固定波束的天线[16,17],具有波束扫描功能的波导缝隙天线实现难度较大,原因在于该种天线是一种行波结构,工作频带较宽,对于特定分布的缝隙结构,在整个宽频带内实现低副瓣性能是非常有挑战性的。2009年,美国加州理工学院喷气推进实验室的Rengarajan教授等人,提出一种Ka波段的波导缝隙天线用于“GLISTIN”空间干涉仪[18],该天线具有低副瓣、高辐射效率的优点。目前工作频段最高的波导缝隙阵列天线,为2012年Cullens等人提出的一种工作在130~180 GHz应用于行星探测器安全着陆雷达的波导缝隙阵列天线[19],该天线通过在金属波导宽边均匀开缝实现漏波辐射,天线的效率比较高;但是由于采用了均匀缝隙,天线的副瓣较高,不利于在成像和目标探测等领域的应用;另外该天线通过沉积工艺加工,工艺复杂且成本较高。该课题组2013年优化了波导阵列设计,提出一种2维波导缝隙阵列,天线增益大大提高[20]。
频率扫描天线可以应用于诸多领域,例如机场外来物探测(Foreign Object Debris, FOD)[21]、行星着陆器成像雷达[22]、机场跑道成像雷达用于雨雾天气飞机安全着陆[23]、人体目标跟踪雷达[24] 等领域。在毫米波太赫兹成像领域,西班牙学者研究了合成孔径技术在频率扫描体制雷达上的应用,并进行了距离向和方位向的2维成像[13]。德国夫琅禾费高频物理研究所的研究人员,首先提出一种基于频率扫描天线的3维成像雷达[25],该雷达工作频率为60 GHz。本文设计的频率扫描天线工作在200 GHz,扫描角度宽,为大场景扫描、高分辨率太赫兹成像打下良好基础。
本文提出一种基于波导缝隙阵列的太赫兹频率扫描天线,考虑到太赫兹波的大气衰减和天线小型化和轻量化的要求,选择工作频段为165~215 GHz,天线通过标准WR-5波导馈电,如图1所示,波导缝隙阵列由31个缝隙单元组成。为增大波束扫描角度,天线辐射段波导宽度比馈电段小。通过微加工车削技术实现波导渐变结构,缝隙阵列刻在0.1 mm厚的铜箔上,铜箔固定在波导槽上侧实现波导缝隙阵列。采用泰勒综合法设计缝隙阵列分布以降低副瓣电平,经仿真和实验测试,波导缝隙天线具有扫描角度大、副瓣电平低的优良特性。据我们所知,本文提出的太赫兹波导缝隙阵列天线为目前工作频率最高的波导缝隙天线。
2. 波导缝隙阵列设计
2.1 缝隙间距
波导缝隙天线的辐射方向图由缝隙分布来控制,波导内传输电磁波可以是行波也可以是驻波,传播驻波时为谐振式缝隙阵,谐振式缝隙阵在波导末端通过短路活塞将结构封闭化,通常谐振式缝隙阵的主瓣指向不变。采用非谐振式波导缝隙阵,在波导末端是开口结构,传播行波带宽较大,行波在传播的同时向外辐射,波束指向随频率变化。交错排列的相邻缝隙的相位差为
2πd/λg−π (其中d为相邻缝隙的间隔,a1为辐射段波导宽边长度,λg=λ/√1−(λ/2a1)2 为波导中电磁波的波长),则辐射波束指向为[20]:θm=arcsin(λλg−λ2d) (1) 其中,
λ 为自由空间中电磁波的波长,为使辐射波束只有一个主瓣,不出现栅瓣,则要满足[20]:d<λλgλ+λg (2) 根据矩形波导的色散特性,当工作频率靠近基模截止频率时,波导中电磁波波长变化范围大,由式(1)可知波束扫描角度大。因此为使工作频率靠近波导截止频率,辐射段波导宽边宽度比馈电段窄,选择a1=0.95 mm, d=0.9 mm,缝隙尺寸一般根据中心频率确定,缝隙宽度通常小于
λ/10 ,缝隙长度经验上为0.48λ 。2.2 缝隙电导分布
通常契比雪夫分布的阵列辐射特性比较好,在相同阵列长度下对给定的副瓣电平,其主瓣是最窄的,副瓣电平是最低的。但是它的等副瓣特性会造成能量的损耗较多,而且两端单元的激励电流的幅度往往比其相邻单元电流的幅度变化大得多,这对缝隙的馈电以及端口反射特性造成较大的损失,而且副瓣电平对激励电平的误差比较敏感。本文采用泰勒分布降低副瓣电平。根据太赫兹成像和目标探测要求,选定副瓣电平SLL=–30 dB,辐射缝隙个数为31时,根据功率传输法[26],对应于泰勒分布的缝隙电导为
gi=U2iq−i+11ηN∑j=1U2j−i∑j=1U2jq−j+1 (3) 其中,
Ui 为第i个缝隙的口径激励幅度,Pri∝U2i ,根据此式缝隙电导可以由给定的口径分布和天线效率计算出。为避免较大的反射,缝隙电导一般不宜太大,为此设定天线效率为80%,则对应于泰勒分布的缝隙电导分布如图2所示。2.3 缝隙偏移量分布
下面求缝隙电导与其尺寸的关系,即缝隙电导函数。Stevenson根据传输线理论和波导的格林函数及边界条件推导了波导缝隙的阻抗计算公式,但是实际中缝隙间互耦很严重,必须要考虑。传统的驻波阵可以通过Eilliot公式优化设计。对于行波阵,缝隙参数在较宽的频带内有所变化,缝隙阵列在波导内外均会产生互耦,互耦的精确计算非常困难,一般通过数值的方法确定缝隙电导函数。采用全波仿真软件Ansys HFSS建立均匀分布的缝隙阵,缝隙单元数为20,模型中缝隙间距,缝隙宽度都与实际天线相同,相邻缝隙交替排布在波导中心线的两侧,变化缝隙偏移量求解端口S参数。假设每个缝隙辐射能量相同,每个缝隙的电导近似相等,且该电导包含了缝隙间的互耦;根据行波近似,可求得不同偏移量对应的缝隙电导[18]。再通过此偏移量与缝隙电导的对应关系,以及上一节中所需要的缝隙电导分布,拟合出所需要的缝隙偏移量。
在仿真中,馈电段波导宽边长度1.3 mm,高度0.65 mm,辐射段漏波波导的宽边长度0.95 mm,高度0.65 mm,缝隙周期0.9 mm,偏移量为0~0.2 mm。求得不同频率的缝隙偏移量变化时对应的缝隙电导。图3给出了170 GHz, 190 GHz和210 GHz 3个频点的缝隙电导随偏移量变化的曲线,从图中可知,相同偏移量不同频率对应的电导不同,因此在某一频率最优的缝隙偏移分布,在其他频率可能不是最优的。由于不同频率的缝隙偏移量理论值不同,将得到的200 GHz的缝隙偏移量分布作为初始值,利用粒子群优化算法[27],优化偏移量分布使得副瓣抑制效果达到频带内最优。最终的优化得到的缝隙偏移量分布如图4所示。
由于此波导缝隙阵列只有一个维度的周期缝隙,因此只在扫描平面有内窄波束;为实现非扫描平面波束宽度的锐化,采用了E面喇叭的方法,使得辐射波束在垂直于扫描平面的维度有较窄的波束宽度;经过仿真对比,最后选择喇叭张角为15°。
3. 仿真与实验测试
3.1 天线加工及测试环境
如图5所示,实际加工的波导缝隙天线,缝隙通过微加工工艺在0.1 mm厚的铜箔上刻缝形成,E面喇叭与波导槽将铜箔夹起来形成波导缝隙天线结构,馈电段为标准波导端口,与太赫兹发射模块的波导端口直接相连。
如图6所示波导缝隙天线测试场景,太赫兹发射模块通过计算机控制输出不同频率的电磁波,波导缝隙天线的辐射波束在xy平面扫描;检波器固定在可以2维移动的扫描平台上,检波器2维运动测量天线的辐射方向图,测量信号经数据采集后传回计算机。
3.2 扫描角度对比
如图7所示,扫描角度仿真与测试比较吻合,实测的扫描角度为–3.8°~53.8°,仿真的扫描角度为–7.37°~–56.98°,在60 GHz带宽内扫描了50°。如图8所示,若以–20 dB副瓣带宽计,在165~215 GHz频带范围内,扫描角度为–4.8°~–45.2°,在50 GHz带宽内扫描了40°,说明该太赫兹波导缝隙天线具有宽角度扫描特性;实测扫描角度比仿真的扫描角度整体向边射方向偏移,误差来自天线加工误差以及测量中天线摆放角度的误差。
3.3 方向图对比
如图8所示实验测得xy平面内扫描波束的归一化方向图,主瓣方向图比较光滑,在扫描平面内波束比较窄,3 dB波束宽度为4°~7°,且波束宽度随着频率的升高而减小;除去两端的频率点外,在165~215 GHz频段内副瓣电平平均小于–20 dB。
如图9所示仿真与实验方向图对比,可以发现在175 GHz, 190 GHz, 210 GHz处,仿真与实测的扫描角度大致相同,说明加工尺寸比较精确,波导缝隙天线中的相位常数变化不大;实测的副瓣电平比扫描的稍低,同时实测的波束宽度比仿真宽,说明金属波导损耗比仿真稍大,缝隙激励电压分布比仿真更加尖锐,使得波束变宽副瓣降低。经过试验验证,波导缝隙天线的增益可达15 dB。
如图10所示仿真与实验副瓣电平对比,在165~215 GHz频段内,仿真副瓣电平在–20 dB附近变化,实验测试的副瓣电平优于–20 dB,在整个频段内在–24 dB附近变化。在频带的两端,实测的副瓣电平都比较高,原因是两端的频率对应的缝隙电导变化较大,副瓣电平抑制效果没有中间频段好。
如图6所示,检波器固定在可以2维移动的扫描平台,记录不同角度的辐射强度,检波器前端连接一个全向性的波导端口,保证对各个入射方向的电磁波增益相同。测量波导缝隙阵的增益时,首先将波导缝隙阵列位置替换为VDI公司的增益为20 dBi的标准圆锥喇叭天线,通过检波器记录标准喇叭天线做发射时的辐射强度;然后再通过检波器测试波导缝隙阵列的辐射强度。两次测试相比较即可得到波导缝隙阵列的增益。经测试在中心频率190 GHz处,波导缝隙阵列增益约为15 dB。
表1比较了太赫兹频段几种典型的频扫天线的性能。频率扫描反射栅天线的突出优点是增益高,但是扫描范围较小。基于漏波或相移原理的频率扫描天线,可以克服扫描角度小的缺陷。几种见报道的波导缝隙行波天线,扫描角度都可以达到40°,但是副瓣电平比较高。本文设计的波导缝隙天线,扫描角度大同时副瓣电平低,由于采用1维阵列,天线增益适中,可以通过扩展到2维阵列的方式提高天线增益。
表 1 太赫兹频率扫描天线性能对比Table 1. Performance comparison of the Terahertz frequency scanning antennas4. 结论
太赫兹成像系统在无损探测和非接触式安检成像等领域有重大的应用价值,但是现有的成像系统多为机械扫描体制,成像时间较长。本文提出将频率扫描天线应用于太赫兹成像系统中,有望实现高帧率成像,并且提出一种基于波导缝隙阵列的太赫兹频率扫描天线。本文重点描述了天线的设计方法,采用泰勒综合法降低副瓣电平,通过软件仿真结合功率传输法设计最优的缝隙分布。实际加工了太赫兹波导缝隙阵列天线,通过太赫兹准光测试系统对天线性能进行测试,实测天线扫描角度超过40°,增益约为15 dB,副瓣电平抑制优于20 dB。本文提出的太赫兹频率扫描天线实现波束在一个维度的快速扫描,可以替代目前太赫兹成像系统中一个维度的机械扫描,接下来研究实现2维快速波束扫描的方法。测试结果表明波导缝隙天线具有大扫描角度和低副瓣的优良特性,可以应用于太赫兹成像和目标探测系统中,在太赫兹实时成像方面有重要的应用价值。
-
表 1 矩阵M中的元素
Table 1. Elements of matrix M
元素位置 表达式 元素位置 表达式 m11(z) ∑i∑jL(ui,uj)zui−uj m12(z) ∑i∑jL(ui,vj)zui−vj m13(z) ∑iS(ui,ui)z2ui−4+∑i∑jS(ui,uj)⏟i≠j m14(z) ∑i∑jS(ui,vj)zui−vj m21(z) ∑i∑jL(vi,uj)zvi−uj m22(z) ∑i∑jL(vi,vj)zvi−vj m23(z) ∑i∑jS(vi,uj)z−vi+uj m24(z) ∑i∑jS(vi,vj)zvi−vj m31(z) ∑iY(ui,ui)z−2ui+4+∑i∑jY(ui,uj)⏟i≠j m32(z) ∑i∑jY(ui,vj)z−ui+vj m33(z) ∑i∑jW(ui,uj)z−ui+uj m34(z) ∑i∑jW(ui,vj)z−ui+vj m41(z) ∑i∑jY(vi,uj)zvi−uj m42(z) ∑i∑jY(vi,vj)zvi−vj m43(z) ∑i∑jW(vi,uj)z−vi+uj m44(z) ∑i∑jW(vi,vj)zvi−vj 表 2 算法运算复杂度比较
Table 2. Algorithm complexity comparison
算法 统计量矩阵 特征值分解 谱峰搜索 GESPRIT O(2N2/L) 8N3/3 KN2/Δr(Rmax NCGESPRIT O(4{N^2}/L) 32{N^3}/3 4K{N^2}/{\Delta _r}({R_{\max } } - {R_{\min } }) + 720{N^2}/{\Delta _\theta } 本文算法 O(4{N^2}/L) 32{N^3}/3 -
[1] JOHNSON R C. Antenna Engineering Handbook[M]. 3rd ed. New York, USA: McGraw-Hill, 1993: 36–39. [2] HUANG Y D and BARKAT M. Near-field multiple source localization by passive sensor array[J]. IEEE Transactions on Antennas and Propagation, 1991, 39(7): 968–975. doi: 10.1109/8.86917 [3] CHALLA R N and SHAMSUNDER S. High-order subspace-based algorithms for passive localization of near-field sources[C]. The 29th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, USA, 1995: 777–781. doi: 10.1109/ACSSC.1995.540806. [4] ZHI Wanjun and CHIA M Y W. Near-field source localization via symmetric subarrays[J]. IEEE Signal Processing Letters, 2007, 14(6): 409–412. doi: 10.1109/LSP.2006.888390 [5] XIE Jian, TAO Haihong, RAO Xuan, et al. Comments on “near-field source localization via symmetric subarrays”[J]. IEEE Signal Processing Letters, 2015, 22(5): 643–644. doi: 10.1109/LSP.2014.2364174 [6] CHARGÉ P, WANG Yide, and SAILLARD J. A non-circular sources direction finding method using polynomial rooting[J]. Signal Processing, 2001, 81(8): 1765–1770. doi: 10.1016/S0165-1684(01)00071-8 [7] SHEN Lei, LIU Zhiwen, and XU Yougen. Parameter estimation of rectilinear signals using multiple noncollocated and diversely polarized antennas[C]. 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 2016: 1–5. doi: 10.1109/RADAR.2016.8059546. [8] XIE Jian, TAO Haihong, RAO Xuan, et al. Real-valued localisation algorithm for near-field non-circular sources[J]. Electronics Letters, 2015, 51(17): 1330–1331. doi: 10.1049/el.2015.0454 [9] XIE Jian, TAO Haihong, RAO Xuan, et al. Passive localization of noncircular sources in the near-field[C]. 2015 16th International Radar Symposium (IRS), Dresden, Germany, 2015: 493–498. doi: 10.1109/IRS.2015.7226309. [10] XIE Jian, TAO Haihong, RAO Xuan, et al. Efficient method of passive localization for near-field noncircular sources[J]. IEEE Antennas and Wireless Propagation Letters, 2015, 14: 1223–1226. doi: 10.1109/LAWP.2015.2399534 [11] KUANG Meidong, WANG Ling, XIE Jian, et al. Real-valued near-field localization of partially polarized noncircular sources with a cross-dipole array[J]. IEEE Access, 2019, 7: 36623–36632. doi: 10.1109/ACCESS.2019.2904837 [12] CHEN Hua, LI Youming, ZHU Weiping, et al. Joint DOA and range estimation for mixed near-field and far-field strictly noncircular sources[C]. 2018 IEEE 23rd International Conference on Digital Signal Processing (DSP), Shanghai, China, 2018: 1–5. doi: 10.1109/ICDSP.2018.8631547. [13] CHEN Hua, LIU Wei, ZHU Weiping, et al. Noncircularity-based localization for mixed near-field and far-field sources with unknown mutual coupling[C]. 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, Canada, 2018: 3236–3240. doi: 10.1109/ICASSP.2018.8462266. [14] FERREOL A, BOYER E, and LARZABAL P. Low-cost algorithm for some bearing estimation methods in presence of separable nuisance parameters[J]. Electronics Letters, 2004, 40(15): 966–967. doi: 10.1049/el:20040537 [15] STOICA P, LARSSON E G, and GERSHMAN A B. The stochastic CRB for array processing: A textbook derivation[J]. IEEE Signal Processing Letters, 2001, 8(5): 148–150. doi: 10.1109/97.917699 期刊类型引用(2)
1. 张涛,袁涛,陈峰峰,张广磊,路成军. 近场条件下数字阵列雷达自适应波束形成技术研究. 电波科学学报. 2023(01): 152-158 . 百度学术
2. 刘振,苏晓龙,刘天鹏,彭勃,陈鑫,刘永祥. 基于矩阵差分的远场和近场混合源定位方法. 雷达学报. 2021(03): 432-442 . 本站查看
其他类型引用(1)
-