辐射方向图解析

分类:电磁场解析 > 天线 | 综合版 2026-04-11
Antenna radiation pattern 3D polar plot showing main lobe, side lobes and null points with directivity contour
辐射方向图解析 — 远场电界分布用E面、H面极坐标图表示

辐射方向图的理论基础

概要

🧑‍🎓

我经常看到辐射方向图的3D绘图,那到底是什么意思?

🎓

那是天线电界强度角度分布的展示。通俗讲,就是"电磁波在哪个方向以多大强度辐射出去"的立体图。一眼就能看出主瓣的方向和宽度、旁瓣电平、零点位置。

🧑‍🎓

这在实际工作中有那么重要吗? 似乎只要知道"朝向哪里"就够了。

🎓

想得太简单了。比如卫星通信的抛物面天线,旁瓣电平必须控制在$-30\,\text{dB}$以下。旁瓣太高的话,会对相邻卫星造成干扰。军用雷达更严格,如果旁瓣泄漏信号,隐身性就完全没有了。FEM分析的精度直接影响设计品质。

🧑‍🎓

那用什么物理量来评估呢?

🎓

主要有4个指标。远场电界$\mathbf{E}(\theta,\phi)$指向性$D$增益$G$半功率波束宽度HPBW。我们逐个讲解。

远场电界表达式

🧑‍🎓

什么叫"远场"?有明确的距离定义吗?

🎓

有的。天线最大尺寸为$D_{\max}$,波长为$\lambda$,那么远场的边界是

$$ r > \frac{2D_{\max}^2}{\lambda} $$
🎓

在这个距离以外,电磁波接近球面波,$\mathbf{E}$和$\mathbf{H}$都与$r$成反比,相互垂直,波阻抗均匀约为$\eta_0 \approx 377\,\Omega$的横向均匀波。

🧑‍🎓

远场电界怎么数学表达?

🎓

从天线表面电流分布$\mathbf{J}(\mathbf{r}')$推导远场电界,得到下式:

$$ \mathbf{E}(\theta,\phi) = -j\omega\mu_0 \frac{e^{-jkr}}{4\pi r} \int_V \mathbf{J}(\mathbf{r}') \, e^{jk\hat{r}\cdot\mathbf{r}'} \, dV' \bigg|_{\text{far-field}} $$
🎓

用方向图函数$\mathbf{F}(\theta,\phi)$简化表示为

$$ \mathbf{E}(r,\theta,\phi) = \frac{e^{-jkr}}{r} \mathbf{F}(\theta,\phi) $$
🎓

关键是,$\mathbf{F}(\theta,\phi)$相当于电流分布的傅里叶变换。因此阵列天线的方向图合成可以看作"空间频率滤波"。

🧑‍🎓

傅里叶变换! 跟信号处理一样的思路。那均匀矩形电流分布的话,方向图是不是像sinc函数?

🎓

正是如此。等幅等相均匀分布主瓣最窄,但旁瓣只能到$-13.2\,\text{dB}$。用泰勒分布或切比雪夫分布给电流加渐进,主瓣稍微宽一点,旁瓣可以降到$-20$~$-40\,\text{dB}$。

辐射强度与指向性

🧑‍🎓

有了方向图函数,接下来要定量评估"有多强"吧?

🎓

对的。首先定义辐射强度$U(\theta,\phi)$,就是单位立体角的辐射功率:

$$ U(\theta,\phi) = \frac{r^2}{2\eta_0} \left| \mathbf{E}(r,\theta,\phi) \right|^2 = \frac{1}{2\eta_0} \left| \mathbf{F}(\theta,\phi) \right|^2 $$
🎓

第二个等号很重要,说明在远场$U$不依赖距离,只是纯粹的角度函数。指向性$D$就是与全向辐射体相比的"峰值聚集度":

$$ D = \frac{4\pi \, U_{\max}}{P_{\text{rad}}} = \frac{4\pi \, U_{\max}}{\displaystyle\int_0^{2\pi}\!\!\int_0^{\pi} U(\theta,\phi)\sin\theta\,d\theta\,d\phi} $$
🧑‍🎓

$P_{\text{rad}}$是向全方向辐射的总功率吧? 数值怎么算?

🎓

在环绕天线的球面$S$上对波印廷矢量进行数值积分。商用求解器比如HFSS,先把近场变换到远场(NFTF),再用辛普森法则或高斯求积在$(\theta,\phi)$网格上积分。通常$\theta$和$\phi$都取1°步长就足够精确了。

增益与HPBW

🧑‍🎓

指向性$D$明白了。但数据表上通常写的是"增益$G$"吧? 有什么不同?

🎓

增益$G$包含了天线的损耗(导体损、介质损、表面波损等),是实际性能指标。用辐射效率$\eta$($0 < \eta \leq 1$)表示为

$$ G = \eta \, D $$
🎓

比如指向性$D = 30\,\text{dBi}$的抛物面天线,如果开口效率(照明效率+溅射损+表面精度损等)$\eta = 0.55$($-2.6\,\text{dB}$),实际增益就约$27.4\,\text{dBi}$了。

🧑‍🎓

HPBW(半功率波束宽度)怎么求?

🎓

就是主瓣顶部功率下降$-3\,\text{dB}$处的两点的角度差:

$$ \text{HPBW} = \theta_{-3\text{dB},\,\text{upper}} - \theta_{-3\text{dB},\,\text{lower}} $$
🎓

比如半波长振子的E面HPBW约$78°$。而20波长口径的天线,$\text{HPBW} \approx 51\lambda/(D_{\max}) \approx 2.55°$,波束就很尖锐了。波束越尖锐,指向性越高,对远处目标的能量集中度就越好。

旁瓣与零点

🧑‍🎓

旁瓣电平(SLL)越低越好对吧? 为啥要这么严格?

🎓

有几个原因。首先是干扰回避。卫星通信的ITU规范(比如ITU-R S.580)规定了旁瓣包络,$G(\theta) \leq 32 - 25\log_{10}\theta$ [dBi]必须满足,否则认证通不过。其次是杂波抑制。气象雷达的旁瓣如果拾取地面杂波,降雨量推估精度就会大打折扣。

🧑‍🎓

零点(null)怎么用?

🎓

零点是辐射功率几乎为零的方向。用零点指向技术,可以把零点故意指向干扰源,从而消除干扰。5G基地站的Massive MIMO中,还有把零点对准人体来降低SAR(比吸收率)的技术已经商用了。

Coffee Break 小知识

为什么用"极坐标图"而不是直角坐标

初次接触辐射方向图的人,常问为什么不用直角坐标而要用极坐标。直角坐标确实更容易读取数值,但极坐标的妙处在于"以天线为中心,直观看出各方向有多强"。旁瓣、反向瓣、零点的位置一眼明白,设计评审时"这个方向图触犯了ITU规范"能瞬间看出来。实务做法是用直角坐标比较精确数值,用极坐标检查整体形状。

辐射方向图的数值计算方法

解析方法的分类

🧑‍🎓

方向图计算除了FEM,还有其他方法吧?

🎓

大概有3种主流方法,各有所长:

方法原理适用场景计算复杂度
FEM(有限元法)体积分割、边元素复杂形状、非均质介质$O(N^{1.5})$~$O(N^2)$
MoM(矩量法)表面积分方程金属构造体、导线$O(N^2)$~$O(N^3)$
FDTD(时域差分法)Yee网格时间演进宽带特性、过渡过程$O(N)$(单步)
🧑‍🎓

怎么选呢?

🎓

实务做法是这样:微带贴片天线、导波管槽阵等有密闭区域的结构用FEM。导线天线、金属板散射这类开放区域只需离散化金属表面的用MoM。超宽带或全频段特性一次计算的话用FDTD。现在混合方法(FEM-BIE、MLFMM加速的MoM)也很主流。

FEM定式化

🧑‍🎓

FEM求方向图时,要解什么方程? 跟结构分析的$[K]\{u\}=\{F\}$一样吗?

🎓

思路差不多,但要解的是矢量亥姆霍兹方程

$$ \nabla \times \left( \frac{1}{\mu_r} \nabla \times \mathbf{E} \right) - k_0^2 \varepsilon_r \mathbf{E} = -jk_0 Z_0 \mathbf{J}_{\text{imp}} $$
🎓

其中$k_0 = \omega/c$是自由空间波数,$Z_0 = \sqrt{\mu_0/\varepsilon_0}$是自由空间波阻抗。FEM用边元素(Nédélec元素)离散化弱形式。边元素自动保证电界切线分量的连续性,消除了节点元素会出现的虚假模式(非物理的寄生解)。

🧑‍🎓

边元素是指自由度在边上的吧。离散化后的矩阵方程是什么样?

🎓

用加权余量法离散后,最终得到:

$$ \left( [S] - k_0^2 [T] \right) \{E\} = \{b\} $$
🎓

$[S]$是旋度-旋度矩阵,$[T]$是质量矩阵,$\{b\}$是励振源项。形式跟结构分析$[K]\{u\}=\{F\}$一样,但矩阵是复数且非正定。因此不能用CG法,要用GMRES或MINRES这样的迭代求解器。

MoM定式化

🧑‍🎓

矩量法跟FEM有什么本质区别?

🎓

FEM是把整个体积都分割网格,MoM只把金属表面上的电流分布作为未知量。电场积分方程(EFIE)是

$$ \hat{n} \times \mathbf{E}^{\text{inc}} = \hat{n} \times \left( jk_0 Z_0 \int_S \left[ \mathbf{J}(\mathbf{r}') + \frac{1}{k_0^2}\nabla'\left(\nabla'\cdot\mathbf{J}(\mathbf{r}')\right) \right] G(\mathbf{r},\mathbf{r}') \, dS' \right) $$
🎓

$G(\mathbf{r},\mathbf{r}') = e^{-jk_0|\mathbf{r}-\mathbf{r}'|}/(4\pi|\mathbf{r}-\mathbf{r}'|)$是自由空间格林函数。用RWG(Rao-Wilton-Glisson)基函数离散化。MoM的强项是不用给空气区域网格。导线天线、偶极子这种开放结构效率高得多。缺点是得到密矩阵,大规模问题要用MLFMM(多层快速多极子法)加速到$O(N\log N)$。

近场到远场的变换

🧑‍🎓

FEM算出来的结果是有限计算域内的,怎么得到远场?

🎓

这很关键。FEM/FDTD的计算域是有限的(用ABC/PML边界封闭),直接得不到远场。因此要用近场-远场变换(NFTF)。在环绕天线的虚拟闭曲面$S'$上取出$\mathbf{E}$和$\mathbf{H}$,转换为等效电流、磁流,再做傅里叶积分。

$$ \mathbf{E}_{\text{far}}(\theta,\phi) = -\frac{jk_0 e^{-jk_0 r}}{4\pi r} \oint_{S'} \left[ Z_0 (\hat{r}\times\mathbf{J}_s)\times\hat{r} - \mathbf{M}_s\times\hat{r} \right] e^{jk_0 \hat{r}\cdot\mathbf{r}'} dS' $$
🎓

其中$\mathbf{J}_s = \hat{n}\times\mathbf{H}$、$\mathbf{M}_s = -\hat{n}\times\mathbf{E}$是等效电流、磁流。HFSS的这个变换是自动化的,用户只需设置"Radiation Boundary"就行。

网格与收敛

🧑‍🎓

方向图计算的网格跟结构分析感觉不一样?

🎓

完全不一样。结构分析在应力集中处细分就行,但电磁场单元大小由波长决定。经验法则:

  • 至少$\lambda/10$以下的单元(二阶元素)
  • $\lambda/20$以下最安心(特别是介质边界附近)
  • 曲面网格曲率半径的$1/5$以下
  • 馈电点附近做局部加密(约$\lambda/50$)

HFSS的自适应网格(AMR)功能会自动加密单元,直到$\Delta S$(S参数变化量)低于阈值(通常0.02)为止。

🧑‍🎓

收敛怎么判定?

🎓

方向图分析必须同时检查S参数和增益。S参数收敛了,增益(特别是旁瓣)还在动的情况常见。网格加倍后增益变化超过$0.3\,\text{dB}$就要警惕。最后要把变化控制在$0.1\,\text{dB}$以内才算安全。

辐射方向图的实务应用

解析流程

🧑‍🎓

实际做方向图分析,要怎么一步步来?

🎓

大概5个阶段:

  1. 几何建模:天线本体+馈电结构+周边构件(罩、外壳)用CAD建立。有对称面就充分利用降低计算量
  2. 材料定义:导体(PEC或有限电导率)、介质基板($\varepsilon_r$、$\tan\delta$)、铁氧体等。注意$\tan\delta$的频率依赖性
  3. 励振设置:波导口(导波管、微带)或集总口(同轴、分立)。指定口的模式次数
  4. 边界+网格:辐射边界(PML或吸收边界)、对称面、初始网格设定
  5. 计算+后处理:S参数、方向图(E/H面切线、3D)、增益、效率全部检查

边界条件设置

🧑‍🎓

辐射边界放在哪里? 紧挨天线可以吗?

🎓

千万别。近场的无源成分会反射回来搞坏结果。指导原则:

  • PML(完全匹配层):离天线表面至少$\lambda/4$。PML本身厚度也要$\lambda/4$~$\lambda/2$
  • 吸收边界条件(ABC):离天线$\lambda/4$~$\lambda/2$。但大入射角反射率会升高,需注意
  • 对称面:PEC对称($\hat{n}\times\mathbf{E}=0$)和PMC对称($\hat{n}\times\mathbf{H}=0$)要按偏波方向区分用
🧑‍🎓

PEC和PMC搞反了会怎样?

🎓

偏波分量会整个反转。比如振子的E面图如果PEC/PMC用反了,$E_\theta$和$E_\phi$互换,峰值方向就偏了90°。这是菜鸟最常犯的错误。正确做法是先拿简单的振子模型验证一遍,然后再上复杂模型。

网格设计要点

🧑‍🎓

网格品质用什么指标评估?

🎓
指标理想值允许范围电磁场特有注意
单元大小/λ< 1/10< 1/6介质中用λ_d = λ/√ε_r
长宽比1.0< 5.0PML层内可例外
雅可比比1.0> 0.3单元退化对边元素有害
曲率分割数≥ 8≥ 4影响圆柱贴片的曲面精度
🎓

关键一点:介质内波长变短。$\varepsilon_r = 10$的基板中波长只有自由空间的$1/\sqrt{10} \approx 0.316$倍,网格也要同比细分。忽略这点的话,微带贴片的共振频率会偏几个百分点。

验证与有效性确认

🧑‍🎓

分析结果怎么验证? 测方向图又很费时…

🎓

分阶段来:

  1. 理论值对比:半波振子($D = 2.15\,\text{dBi}$、$R_{\text{rad}} = 73.1\,\Omega$)或矩形开口(sinc方向图)验证求解器配置是否正确
  2. 交叉工具验证:同一模型用FEM和MoM都算一遍,增益和旁瓣在$0.5\,\text{dB}$以内就基本对
  3. 网格收敛:至少3个网格密度等级,增益变化$0.1\,\text{dB}$以内收敛
  4. 与测量对比:消声室或紧凑靶场测量。测量系统误差通常$\pm 0.5$~$\pm 1.0\,\text{dBi}$
Coffee Break 小知识

天线方向图测量室的秘密

仿真值和测量值差异大多数来自测量环境。理想消声室误差在±0.5dBi以内,但实际的EMC试验室会因吸音材料老化和地面反射导致±1~2dBi偏差。拥有大型紧凑靶场的测试设备国内不超50套,价格以亿元计。"仿真正确、测量有问题"的情况在天线业界是公开秘密。

辐射方向图的软件比较

主要工具比较

🧑‍🎓

方向图分析用什么商用工具比较好?

🎓

主流工具对比如下:

工具开发方主求解器方向图支持特点
Ansys HFSSAnsys Inc.FEM(频域)3D/2D切线,NFTF自动自适应网格强大。业界标准
CST Studio SuiteDassault (SIMULIA)FDTD、FEM、MoM远场监视器自动宽带分析强。EMC集成
FEKOAltairMoM、MLFMM、PO/UTD3D、大规模阵列混合方法强
COMSOL (RF模块)COMSOL ABFEM远场计算节点多物理场耦合容易
NEC (4nec2)开源MoM3D方向图免费。仅限导线天线
🧑‍🎓

HFSS和CST的结果能对得上吗?

🎓

主瓣增益通常在$0.2$~$0.5\,\text{dB}$以内。旁瓣可能差$1$~$2\,\text{dB}$。特别是$-30\,\text{dB}$以下的低旁瓣,网格质量和边界条件微妙差别就会显著影响。用多个工具交叉验证是最保险的。

选择指南

🧑‍🎓

按用途怎么选?

🎓
  • 微带贴片、集成天线:HFSS首选。自适应网格收敛稳定
  • 超大阵列(100素子以上):FEKO的MLFMM+区域分解。MoM开放域自然处理
  • UWB、超宽带脉冲:CST时域求解器(一次计算全频段)
  • 天线+筐体热耦合:COMSOL的RF+Heat Transfer模块
  • 学习、导线天线概念:NEC(免费)足够

辐射方向图的先进研究

大规模MIMO与波束成形

🧑‍🎓

5G的Massive MIMO跟方向图分析怎么关联?

🎓

Massive MIMO是64~256素子阵列,动态成形波束。方向图分析面临3大挑战:

  • 计算规模:256素×数GHz自由度达数千万。需用区域分解或等效原理子域法
  • 互耦合:素子间相互阻抗改变各素子的方向图(有源素子方向图)。单体方向图×AF不够
  • 波束权重量化:数字移相器是6~8比特量化,零点方向会偏离设计值。仿真常用连续值,与实机脱节
🧑‍🎓

256素子全波全场仿真能做吗?

🎓

实际不行。业界用的是混合法。中心几素子用全波求解(HFSS/CST)得有源素子方向图,再复制到阵列各位置,与阵列因子合成。周边素子的边缘效应用子阵补正。这样计算时间削减100倍以上。

机器学习融合

🧑‍🎓

最近听说用机器学习预测方向图…

🎓

有两个热点方向:

  • 代理模型:从参数扫描数据训练DNN快速预测方向图。优化中成本函数计算快1000倍。但外推能力弱
  • PINN(物理信息神经网络):将麦克斯韦方程编入损失函数保证物理一致。无需网格,大规模问题有优势,但精度还达不到FEM

实际应用中,两者都是"用FEM/MoM替代",而是"快速缩小初期设计探索空间"的工具。

Coffee Break 小知识

5G基站"绕过人体"的零点控制实际

5G NR的大规模MIMO已实现让方向图的零点指向人体方向来降低SAR的技术。基站通过GPS和惯导确定用户端末位置,实时更新波束权重。但实际问题是波束权重的量子化误差让零点容易"堵上"。精确的仿真与测量对比是必须的。这是个对"仿真≠现实"的深刻洞察。

辐射方向图的故障排除

方向图计算的典型问题

🧑‍🎓

方向图仿真常见的"坑"有哪些?

🎓
现象原因对策
增益比理论值高3dB以上对称面设置错(全空间功率未正确除)检查PEC/PMC,去掉对称重新算
旁瓣左右不对称网格非对称或馈电点偏移网格对称化、确认几何精度
反向瓣异常大辐射边界距离不足或PML参数不当离天线λ/4以上,检查PML设置
共振频率偏差介质网格不足或材料参数错λ_d/10细分,确认ε_r、tanδ
交叉偏波理论上为零却有输出网格非对称或数值误差对称条件严格设置,细分网格
🧑‍🎓

对称面造成3dB差别,那用对称计算真的能省计算量吗?

🎓

能省一半。但前提是$P_{\text{rad}}$的积分要自动扩展到全空间。HFSS有自动补正,但COMSOL需手工乘2。这是新手最常掉的坑。

工具别错误对策

🧑‍🎓

HFSS和CST特有的方向图错误有哪些?

🎓

常遇到的:

🎓

Ansys HFSS:

  • "自适应网格未收敛":初始网格太粗或辐射边界太近。手动设λseeding
  • "波导口模式计算失败":波导口尺寸不当(过大或过小)。调整口宽到基板的5~10倍
  • "负辐射功率":辐射边界(Radiation Boundary)与PML干扰。检查距离
🎓

CST Studio Suite:

  • "远场监视无效数据":远场频率超出计算范围或时间信号衰减不足(能量衰减<0.01%)
  • "馈电点网格太小":自动生成网格在口处坍缩。Local Mesh指定最小单元大小
🧑‍🎓

归根结底"仿真合不上"的万能法则是什么?

🎓
  1. 最小模型验证:用半波振子这种有解析解的结构,验证求解器配置
  2. 一次改一个参数:网格、边界、口设置逐个调。不要同时改多个
  3. 多工具交叉验证:把HFSS结果跟NEC或MATLAB天线工具箱对比
  4. 回到物理直觉:50dBi增益不可能存在。全向辐射0dBi、半波振子2.15dBi这个感觉要有
  5. 相关仿真器

    用这个领域的交互式仿真器来体验理论

    仿真器列表

    相关领域

    电磁场解析耦合解析结构分析
    这篇文章的评价
    谢谢您的反馈!
    有帮助
    更详细
    报告
    错误
    有帮助
    0
    更详细
    0
    报告错误
    0
    Written by NovaSolver Contributors
    Anonymous Engineers & AI — 网站地图
    查看资料