应用于太赫兹成像的弧形多收多发阵列的设计方法 下载: 899次
1 引言
太赫兹波成像以其非电离辐射、高分辨及可穿透衣物等特性,在人体安检成像、无损检测等领域受到了越来越多的关注[1-3]。为了实现近场高分辨成像,天线阵列要有足够大的孔径,且阵列中天线间隔距离需满足空间采样定理。随着成像频率的提高,所需的天线数目迅速增加,这在实际系统中难以实现,而多收多发(MIMO)成像的出现为太赫兹波近场成像提供了一种新的思路[4]。MIMO成像是一种利用收发天线单元的多重组合来实现以较少的天线数目生成图像的方法,能够在保证图像分辨率的同时显著降低系统的成本及复杂度。基于MIMO阵列的太赫兹波成像系统的关键是阵列设计,它对成像质量、参数估计和目标检测等都有直接的影响。
基于等效阵列概念进行MIMO阵列设计是远场应用中一种有效的方式。在远场条件下,MIMO阵列可以等效为单收单发(SISO)阵列,该SISO阵列被称为MIMO阵列的等效阵列,等效阵列中天线的空间位置矢量是MIMO阵列中发射天线的空间位置矢量和接收天线空间位置矢量的卷积[5]。国内外一些研究人员已对该方法在近场应用中的可行性进行了初步探究,典型的有Zhuge等[6-7]设计的4发8收稀疏阵列以及12发13收的正交阵列,Ge等[8]提出的应用于近场成像的MIMO稀疏阵列,Tan等[9-10]提出的二维MIMO稀疏阵列,以及Ahmed等[11-12]所采用的正方形稀疏阵列,但上述研究仅基于直线MIMO阵列或平面MIMO阵列。由于镜面反射的影响,人体等曲面成像目标侧面回波信号在平面阵列及直线阵列等平面扫描方式下难以被接收,得到的成像物体侧面图相对模糊。而对弧形阵列进行单次扫描即可获得成像目标各角度的信息,能在一定程度上克服平面扫描时侧面成像质量差的缺点[13]。
本文针对太赫兹安检成像的需求,基于等效阵列概念提出了一种弧形MIMO阵列的设计方法,并通过仿真对该方法在弧形MIMO阵列设计中的有效性进行验证。
2 成像模型
式中:
式中:
由任一发射天线
式中:
式中:
式中:
式中:
3 弧形MIMO阵列设计方法
基于等效阵列概念的弧形MIMO阵列设计方法的基本原理为:首先设计满足成像性能要求的等效弧形SISO阵列,然后将其分解为相应的弧形MIMO阵列。具体而言,根据发射信号的中心波长
式中:
式中:
综合(7)~(9)式可知,当
式中:成像距离
进一步整理(10)式和(11)式可得
依据成像的空间采样条件[15],等效阵列天线阵元之间的间隔Δ
式中:
利用已获得的等效阵列参数,分解得到相应的弧形MIMO阵列参数。其中,弧形MIMO阵列参数应满足:
式中:
式中:
由于分解得到的弧形MIMO阵列和等效阵列不是一一对应的关系,故最后需要根据系统要求从分解得到的多个弧形MIMO阵列中选择最为合适的弧形MIMO阵列。
4 设计实例及仿真验证
4.1 设计实例
考虑一种太赫兹成像的场景,成像目标对应的成像半径为0.1 m,其中心点与天线阵列的距离为0.35 m,要求系统的方位向角分辨率优于0.006 rad,发射信号的中心频率设定为140 GHz。则由(12)、(14)式可得
取
图 2. 6种50发60收弧形MIMO阵列及其等效阵列的结构示意图
Fig. 2. Topologies of arc MIMO arrays with 50 transmitters and 60 receivers as well as their equivalent arrays
4.2 仿真验证
弧形MIMO阵列2及其等效阵列点扩展函数(PSF)的仿真结果如
图 3. 弧形MIMO阵列2及其等效阵列PSF(d=0.35 m,f=140 GHz)
Fig. 3. PSF for MIMO array 2 and its equivalent array (d=0.35 m, f=140 GHz)
图 4. 弧形MIMO阵列多点目标方位向成像结果。(a)点目标位置(0.1,π/2),(0.1,π/2-0.012);(b)点目标位置(0.1,π/2),(0.1,π/2-0.009);(c)点目标位置(0.1,π/2),(0.1,π/2-0.006)
Fig. 4. Azimuthal imaging results of multiple points for arc MIMO array. (a) Points located at (0.1,π/2),(0.1,π/2-0.012); (b) points located at (0.1,π/2), (0.1,π/2-0.009); (c) points located at (0.1,π/2), (0.1,π/2-0.006)
接下来对比直线MIMO阵列和弧形MIMO阵列对曲面目标的仿真成像结果。
图 5. 弧形MIMO阵列及其等效阵列对位于不同成像半径的点目标仿真成像结果。(a)点目标位置(0.12,π/2);(b)点目标位置(0.16,π/2);(c)点目标位置(0.20,π/2)
Fig. 5. Imaging results of points located at different imaging radii for arc MIMO array and its equivalent array. (a) Point located at (0.12,π/2); (b) point located at (0.16,π/2); (c) point located at (0.20,π/2)
图 6. MIMO阵列成像系统布置图。(a)弧形MIMO阵列;(b)直线形MIMO阵列
Fig. 6. Layouts of MIMO imaging systems. (a) Arc MIMO array; (b) linear MIMO array
图 7. 两种MIMO阵列的FEKO仿真成像结果。(a)弧形MIMO阵列;(b)直线形MIMO阵列
Fig. 7. Simulation imaging results of two MIMO arrays based on FEKO. (a) Arc MIMO array; (b) linear MIMO array
5 结论
从弧形MIMO阵列的成像模型出发,提出了一种基于等效阵列概念的弧形MIMO阵列设计方法。在成像目标对应的成像半径为0.1 m、中心点与天线阵列的距离为0.35 m时,弧形MIMO阵列的方位向角分辨率
[1] 张镜水, 朱维文, 赵跃进, 等. 基于曲柄摇杆机构的被动太赫兹波成像系统[J]. 光学学报, 2013, 33(12): 1211006.
[2] 闫昱琪, 赵成强, 徐文东, 等. 太赫兹主动关联成像技术研究[J]. 中国激光, 2018, 45(8): 0814001.
[3] 李文军, 王天一, 周宇, 等. 多层胶接结构胶层空气缺陷的太赫兹无损检测[J]. 光学学报, 2017, 37(1): 0111002.
[5] 粟毅, 朱宇涛, 郁文贤, 等. 多通道雷达天线阵列的设计理论与算法[J]. 中国科学:信息科学, 2010, 40(10): 1372-1383.
[8] 葛桐羽, 经文, 赵磊, 等. 毫米波多输入多输出雷达稀疏阵列近场成像[J]. 强激光与粒子束, 2016, 28(11): 113101.
[13] Sheen D M. McMakin D L, Hall T E. Cylindrical millimeter-wave imaging technique and applications[J]. Proceedings of SPIE, 2006, 6211: 62110A.
[15] Sheen D M. McMakin D L, Hall T E. Three-dimensional millimeter-wave imaging for concealed weapon detection[J]. IEEE Transactions on Microwave Theory and Techniques, 2001, 49(9): 1581-1592.
Article Outline
于洋, 乔灵博, 游燕, 赵自然. 应用于太赫兹成像的弧形多收多发阵列的设计方法[J]. 中国激光, 2019, 46(6): 0614014. Yang Yu, Lingbo Qiao, Yan You, Ziran Zhao. Design Method of Arc Multiple-Input Multiple-Output Array for Terahertz Imaging[J]. Chinese Journal of Lasers, 2019, 46(6): 0614014.