有限大底板微带天线分析软件包的研究与实现

刘汉 姜水桥 李有

摘   要：开发了一款用于有限大底板微带天线计算的软件包. 利用矩量法对微带天线进行三角形和四面体剖分，將阻抗矩阵分为导体对导体作用、导体对介质作用、介质对导体作用和介质对介质作用的4类阻抗矩阵;通过计算4类阻抗矩阵，求得电流系数，完成对天线性能的分析;最后基于矩阵实验室编程环境，实现软件包. 通过与FEKO软件的仿真结果对比表明，该软件包操作简单、可视性强、计算精度较高，可用于有限大底板微带天线的分析和计算.

关键词：有限大底板;微带天线;矩阵实验室;软件包

中图分类号：TN823                                     文献标志码：A

Research and Implementation of Software Package

for Analysis Microstrip Antenna with Finite Base Plate

LIU Han?，JIANG Shuiqiao，LI You

（College of Information and Communication，National University of Defense Technology，Wuhan 430010，China）

Abstract：A software package for the calculation of microstrip antenna with finite base plate is developed. The microtrip antenna is partitioned into triangle and tetrahedron using moment method. The impedance matrix is divided into four types：conductor-to-conductor，conductor-to-medium，medium-to-conductor and medium-to-medium. By calculating the four types of impedance matrix，the current coefficient is obtained and the performance of the antenna is analyzed. Finally，the software package is implemented based on matrix laboratory（MATLAB）  programming environment. Compared with the simulation results of FEKO software，the software package has the advantages of simple operation，high visibility and high calculation accuracy. It can be used for the analysis and calculation of microstrip antenna with finite base plate.

Key words：finite base plate;microstrip antenna;matrix laboratory（MATLAB）;software package

随着无线通信技术的发展，具有良好性能的天线吸引了科研人员的目光，微带天线[1-2]由于具有尺寸小、剖面低、易于集成等优点，逐渐成为了研究热点. 微带天线理论研究方法主要有传输线模型法、空腔模型法和全波分析法等. 对于传统的微带天线，在基板厚度h和波长λ满足h<0.01λ时，以上方法均适用，但是在h≥0.01λ情况下，必须考虑场在贴片垂直方向的变化，此时，传输线模型法和空腔模型法不适用，只能应用全波分析法. 全波分析法中的矩量法[3]是基于积分运算的算法，具有很高的精度，目前广泛地应用于天线的研究设计中. 本文将利用矩量法对微带天线进行研究，实现有限大底板微带天线分析的软件包，并将软件包仿真结果与FEKO[4-5]计算结果进行对比，验证软件包的正确性. 该软件包操作简单、可视性强、计算精度较高，并且可以与优化算法[6-7]结合，完成对微带天线的自动优化设计.

1   RWG-SWG 混合基函数微带天线的矩量

法分析

矩量法中所采用的格林函数[8]表征的常常是一些特殊的边界情形，如平面分层、柱面分层及球面分层天线等，且它们在推导时，理论上都认为这些分层结构是无限大的，故可以利用格林函数法进行研究. 对于有限大结构天线的分析，难以推导出其对应的并矢格林函数，故不能采用类似格林函数的分析方法，但能将微带天线的介质以四面体剖分，利用定义在四面体对上的基函数（SWG）[9]模拟介质中的电位移矢量和三角面对基函数（RWG）[10]模拟天线导体部分的面电流分布，那么就可直接利用自由空间中的格林函数实现有限大介质和金属混合体的仿真计算.

对于有限大底板的平面微带天线，仅仅依靠面剖分无法进行有效计算，需要进行体面混合剖分. 在计算中，可以选取RWG-SWG混合基函数. 其主要计算过程如下.

对于天线金属部分，利用导体表面切向电场边界条件，可以得到：

在介质体内，电位移矢量与入射电场和反射电场满足：

式（14）中：Zcc为导体对导体作用的阻抗矩阵;Zcd为介质对导体作用的阻抗矩阵;Zdc为导体对介质作用的阻抗矩陣;Zdd为介质对介质作用的阻抗矩阵，如图1所示.

通过上述阻抗矩阵的建立，通过体面积分方程[11]的MATLAB编程就可以分析天线的性能.

2   基于RWG-SWG的计算微带天线软件包

的实现

为了能够对微带天线进行直观的分析计算，根据上述矩量法理论及推导的计算公式，利用MATLAB软件，开发了一个可用于分析有限大底板微带天线的软件包，主界面如图2所示. 它包含的功能如下.

图2   开发的介质及底板有限大微带天线的软件包主界面

Fig.2   The main interface of the software package for

the finite microstrip antenna with dielectric and substrate

2.1   参数转换功能

它可以将Feko的.out文件直接读入，转换成element.txt，edge.txt，ted.txt文件，分别包含金属三角形面元信息，金属三角形边信息，介质四面体信息;也可以分3次分别加载转换好的txt文件.

2.2   参数的输入功能

参数输入功能包括天线计算参数设置部分，界面如图3所示，包括介电常数设置、起止频率设置、辐射计算频点设置.

此外还可以进行天线类型和馈电设置，如图4所示. 它包括数据转换部分，其功能与主界面的一致. 选择天线类型部分主要确定天线是平面、柱面还是球面[12]，给了形状后能够便于确定馈电点位置. 馈电边信号确定是通过选择馈电边所在的三角形实现，点击馈电面选择按钮，会出现图5界面，通过放大、平移方式设置好位置后，点击空格键鼠标会变为十字线，将十字线中心对准要选择的三角形，点击后变为白色即选中，一次可以选中多个，选完后按enter键，即退出回到图4界面，同时馈电编号和相应的正负三角形面号将显示在相应窗口内.

删除三角形贴片功能与馈电面选择类似，它是为天线优化设置的，通过人工去除贴片，达到优化天线的目的.

2.3   参数显示功能

参数显示功能如图6所示. 显示内容包括输入阻抗[13]、S11、驻波比、三维立体方向图和天线几何结构. 在进行方向图计算时，显示方向图同时还显示天线的部分结构，便于方向判别.

2.4   加载计算结果功能

结果加载功能是将计算好的结果存储的数据调出来显示和分析，点击加载计算结果会弹出文件对话框，选择所要加载的文件，就可以显示计算结果.

该软件包的主要特点有：

●   提供了两种方法，可以将FEKOCAD中建立的几何结构剖分网格导入到计算程序中，并可以通过图形输出窗口三维显示.

●   理论上，可以实现介质及金属底板有限大平面及共形微带天线的仿真分析与计算.

●   可以通过界面设置馈电点和移除三角形贴片，操作方便.

●   提供了6种计算结果查看功能，可以直接加载已计算好的结果数据并进行结果查看;另外，实现三维场分布与天线模型相结合，直观性更强.

下面通过本文开发的软件包仿真一款天线，将其得到的仿真结果与FEKO计算结果进行对比，以证明该软件包计算的正确性.

图7为所研究天线的结构示意图. 该天线制作在相对介电常数εr = 2.6、厚度H = 1 mm的聚四氟乙烯基板上，基板尺寸为15 mm × 15 mm. 天线具体尺寸参数如表1所示. FEKO仿真与软件包计算所得到的S11对比如图8所示.

图8 为天线的FEKO仿真结果与软件包计算结果对比图. 从图8看出，FEKO仿真结果与软件包计算结果比较吻合，证明了该软件包的正确性和有效性. 两者之间的误差来源于算法之间的固有误差，此外在计算过程中采取了近似，所以也会导致误差.

3   结   论

本文首先利用RWG-SWG混合基函数的矩量法对微带天线的分析计算，将计算矩阵分为4类：导体对导体作用、导体对介质作用、介质对导体作用和介质对介质作用，得到阻抗矩阵和电流系数的表达式. 在上述分析的基础上，利用MATLAB软件开发了友好界面的软件包，该软件包操作简单、可视性强、计算精度较高，可用于分析介质及金属底板有限大微带天线，为以后结合优化算法设计微带天线提供了有效的方法.

参考文献

[1]    PRINCY M P，KRISHNAMOORTHY K，MOHAMMAD S S，et al. Dispersion-engineered transmission line loaded llot antenna for UWB applications [J]. IEEE Antennas and Wireless Propagation Letters，2019，18（2）：323—327.

[2]    YAO H，ZHAO Z，DING X，et al. Compact UWB slot antenna utilizing traveling-wave mode based on slotline transitions [J]. IEEE Transactions on Antennas and Propagation，2019，67（1）：140—150.

[3]    XU Y，YANG H，SHEN R，et al. Scattering analysis of multiobject electromagnetic systems using stepwise method of moment [J]. IEEE Transactions on Antennas and Propagation，2019，67（3）： 1740—1747.

[4]    CHAI S R，GUO L X，LI K，et al. Combining CS with FEKO for fast target characteristic acquisition [J]. IEEE Transactions on Antennas and Propagation，2018，66（5）：2494—2504.

[5]    STEVEN W. Radiated fields of patch antennas——A comparison of simulated and analytic results[C]//2019 International Applied Computational Electromagnetics Society Symposium（ACES）. Miami：IEEE，2019：1—2.

[6]    DAI C，SONG Q. Heuristic computing methods for contact plan design in the spatial-node-based Internet of Everything [J]. China Communications，2019，16（3）：53—68.

[7]    HICHAM C，MEHDY K，OKEZIE O，et al. Simplified speed control of permanent magnet synchronous motors using？genetic algorithms [J]. IEEE Transactions on Power Electronics，2019，34（4）：3563—3574.

[8]    EDA K K A，ALI Y. Numerical calculation of 2-D inhomogeneous media Green's function and some applications in electromagnetic scattering problems [J]. IEEE Transactions on Antennas and Propagation，2019，67（1）：369—377.

[9]    ZHANG L M，SHENG X Q. Discontinuous Galerkin volume integral equation solution of scattering from inhomogeneous dielectric objects by using the SWG basis function [J]. IEEE Transactions on Antennas and Propagation，2017，65（3）：1500—1504.

[10]  SEUNG M S. Implementation of loop co-tree basis functions for EFIE formulation at low frequency [J]. IET Microwaves，Antennas and Propagation，2019，13（1）：118—121.

[11]  LI X J，LIN L，CHEN Y P，et al. VSIE-based domain decomposition method with simplified prism vector basis functions for planar thin dielectric-conductor composite objects [J]. IEEE Antennas and Wireless Propagation Letters，2018，17（9）：1608—1612.

[12]  ZHU Y J，LI J，ZHANG X，et al. A 3-D printed spherical antenna with bandwidth enhancement under operation of dual resonance [J]. IEEE Access，2020（8）：19345—19352.

[13]  周逢道，連士博，徐飞，等. 近地表电磁探测发射系统控制技术 [J]. 湖南大学学报（自然科学版），2017，44（8）：100—105.

ZHOU F D，LIAN S B，XU F，et al. Control Technology of near-surface electromagnetic detection transmitting system [J]. Journal of Hunan University （Natural Sciences），2017，44（8）：100—105.（In Chinese）