天线仿真设计Matlab天线工具箱-电子发烧友网

学个Antenna是以天线仿真和调试为主，理论原理为辅的干货天线技术专栏，包括天线入门知识以及各类天线的原理简介、仿真软件建模、设计、调试过程及思路。如有想看到的内容或技术问题，可以在文尾写下留言。

‍‍‍‍‍‍‍‍ 摘要：

提到天线仿真设计，大部分人会想到HFSS、CST、FEKO这几种常用的电磁仿真软件。不过号称“除了生孩子之外就是万能”的MATLAB在2015版本就已推出了天线工具箱-Antenna Toolbox。经过几年的迭代优化，现在的工具箱已经可以实现设计、分析和可视化天线单元和天线阵列。

本文使用的软件为MATLAB 2019a版本，文中编写的程序在低版本环境下可能会与提供的仿真结果有一定偏差！

0 1 简单介绍

Antenna Toolbox 提供了设计、分析天线单元和阵列并使其可视化的功能和应用程序。用户可以使用具有参数化几何结构的预定义单元或任意平面单元设计独立的天线并建立天线阵列。

Antenna Toolbox 使用矩量法 (MoM) 计算端口属性（如阻抗）、表面属性（如电流和电荷分布）以及电场属性（如近场和远场辐射方向图）。您可通过 2D 和 3D 的形式展示天线的几何结构以及分析结果。

也可将天线和阵列集成到无线系统当中，并使用阻抗分析方法设计匹配网络。Antenna Toolbox 提供了用于仿真波束成形和波束控制算法的辐射方向图。可从您的设计生成 Geber 文件，用于生产印刷电路板 (PCB) 天线。您可以在汽车或飞机等大型平台上安装天线，并分析结构对天线性能的影响。借助站点查看器，您可以使用各种传播模型在 3D 地形图上可视化天线覆盖范围。

0 2 天线单元设计

MATLAB的天线工具箱有函数库支持脚本式开发，亦有基于APPdesigner开发的傻瓜式天线建模仿真软件。这里先简单科普下脚本开发，当“调包侠”有多爽！

%偶极子天线单元设计clear;clc;%设置工作频率（Hz）fre=2*1e9;%创建天线单元种类：偶极子subobject=design(dipole,fre);%设置激励端口阻抗subobject.Load.Impedance=100;str1=[];str2=[];for num=1:3%设置天线尺寸 subobject.Length=0.058+0.02*(num-1); %S参数计算 S_data=sparameters(subobject,0.5*fre1.5*fre,201); subplot(1,2,1);h=rfplot(S_data);hold on; str1=[str1,strcat("偶极子长度(mm)：",num2str(subobject.Length*1000))]; legend(str1); subplot(1,2,2);smith(S_data,1,1);hold on; legend(str1);endholdoff;

CST仿真结果：宽度1.5mm，长度50.8mm,70.8mm,90.8mm

对比MATLAB快速计算和CST建模仿真-时域求解器下的仿真结果，可以看出有一定频偏差异。精准化仿真还是得靠HFSS、CST这类商业电磁仿真专业软件，不过前期进行快速仿真、参量预估，这个工具箱还是能发挥一席之地的！

%在上面的代码后追加%设定偶极子长度为0.078msubobject.Length=0.078;%绘制3D方向图figure(2);pattern(subobject,fre);%绘制2D方向图figure(3);patternElevation(subobject,fre);%绘制电流分布图figure(4);current(subobject,fre);view(90,0);

单元快速仿真计算结果

下面介绍基于APPdesigner开发的傻瓜式天线建模仿真软件。在matlab命令窗口输入antennaDesigner即可。进入界面后点击New，常用的简单天线种类应有尽有，特别方便，程序比较简单、非常容易上手！

点击上图最右边的Export按钮，可以将设计的天线导出为matlab脚本文件，方便开发者二次修改、维护。

0 3 阵列设计

%设置频率fre=2*1e9;%创建天线单元种类subobject=design(dipole,fre);subobject.Load.Impedance=100;%设置线阵arrayobject=design(linearArray('Element',subobject),fre,subobject);%设置单元数量arrayobject.NumElements=20;%设置单元间隔0.6倍工作波长arrayobject.ElementSpacing=0.15*0.6;%设置单元激励相差arrayobject.PhaseShift=zeros(1,20);%设置单元激励幅度arrayobject.AmplitudeTaper=ones(1,20);figure(1);layout(arrayobject);figure(2);pattern(arrayobject,fre);figure(3);patternAzimuth(arrayobject,fre);

线阵快速仿真计算结果-等幅同相激励

点击任意小图，即可放大观看

从上面仿真结果的第三张图片可以看出，组阵方向的2D方向图主波束较窄，主副瓣比仅有13.26dB，与等幅同相激励的理论值比较符合，下面尝试模拟切比雪夫低副瓣阵列。

设定36dB的主副瓣比，得到的20单元线阵激励电流幅度分布为：

%替换上述程序的幅度分布%设置单元激励幅度arrayobject.AmplitudeTaper=... [0.1749 0.2065 0.3091 0.4268... 0.5529 0.6790 0.7955 0.8930... 0.9632 1.0000 1.0000 0.9632... 0.8930 0.7955 0.6790 0.5529... 0.4268 0.3091 0.2065 0.1749];