波束成形天线阵列的电磁分析
波束成形天线阵列的电磁的理论基础
波束成形是什么
波束成形在5G中经常听到，具体是做什么的呢？
简单来说，就是通过控制多个天线单元的相位，将无线电波束集中到特定方向的技术。就像把几十个扬声器排成一列，错开时间发出声音，只有某个方向的声音会变大吧？波束成形就是用无线电波实现同样的原理。
原来如此！那5G基站是怎么使用的呢？
5G基站使用128单元等大规模MIMO阵列，为每个用户单独定向波束。这相当于从传统的"照亮整个房间的荧光灯"方式，进化到了"用聚光灯照亮每个人"的方式。这样SINR（信号与干扰加噪声比）会急剧上升。
用电磁分析处理波束成形时，特别困难的地方在哪里？
是准确建模单元间的互耦（mutual coupling）。如果忽略这一点，旁瓣会比预期更大，零陷方向也会偏移。特别是在 d < 0.5λ 的密集配置下，互耦非常强烈，仅靠单个单元方向图的乘积（Pattern Multiplication）完全无法保证精度。
阵列因子（Array Factor）
阵列因子的公式在教科书上见过，但总感觉不太明白...
对于N单元的等间距线性阵列（ULA: Uniform Linear Array），阵列因子是这样的：
各参数的含义如下：
- $N$ — 阵列单元数
- $w_n$ — 第$n$单元的幅度权重（用于锥削）
- $k = 2\pi/\lambda$ — 波数（自由空间）
- $d$ — 单元间距
- $\theta$ — 与侧射方向的夹角
- $\beta_n$ — 施加给第$n$单元的相移
直观地说，每个单元发出的电波，其"波前到达时间"会相差 $d \sin\theta$。这个差值与人为添加的相位 $\beta_n$ 相互抵消的方向，所有单元的电波会相互增强——那就是主波束。
例如，对于N=8、半波长间距 d=λ/2 的阵列，主波束的3dB波束宽度大概是多少呢？
等幅度的情况下，侧射方向的HPBW（半功率波束宽度）近似为：
代入 $N=8$, $d=\lambda/2$ 得到 $\mathrm{HPBW} \approx 0.886/(8 \times 0.5) = 0.221\,\mathrm{rad} \approx 12.7°$。单元数越多，波束越尖锐。5G的128单元阵列可以聚焦到约 $1°$。
波束转向与相位控制
要怎么把波束转向任意方向呢？
将每个单元的相移 $\beta_n$ 设置如下：
$\theta_s$ 就是转向角度。设置这个相位后，在 $\theta = \theta_s$ 方向上，AF中的指数项变为零，所有单元同相叠加。例如，在28GHz的5G基站中，想要转向 $\theta_s = 30°$ 时，相邻单元间的相位差为 $\Delta\beta = -k \cdot d \cdot \sin 30° = -\pi/2$（当d=λ/2时）。
我听说转向后波束宽度会变化...
你注意到了关键点。波束转向时，阵列的有效口径会缩小为 $\cos\theta_s$ 倍。因此波束宽度会变宽为 $\mathrm{HPBW}(\theta_s) \approx \mathrm{HPBW}_0 / \cos\theta_s$。$\theta_s = 60°$ 时是侧射方向的2倍。实际工作中，在 $\pm 60°$ 以上的转向时增益下降严重，所以通常用扇区划分来应对。
栅瓣条件
栅瓣是什么？和旁瓣不一样吗？
栅瓣是与主波束强度相同的伪波束。旁瓣比主波束弱，但栅瓣是在非预期方向出现的、与主波束同等级的"另一个主波束"。其物理原理与衍射光栅（diffraction grating）相同，因此得名。
栅瓣不进入可见区域（$|\sin\theta| \le 1$）的条件是：
| 转向角 $\theta_s$ | 单元间距上限 | 28GHz频段具体值 |
|---|---|---|
| 0° (侧射) | $d < \lambda$ | $d < 10.7\,\mathrm{mm}$ |
| 30° | $d < 0.67\lambda$ | $d < 7.1\,\mathrm{mm}$ |
| 60° | $d < 0.54\lambda$ | $d < 5.7\,\mathrm{mm}$ |
| 全方向对应 | $d < 0.5\lambda$（安全侧） | $d < 5.35\,\mathrm{mm}$ |
实际工作中 $d = 0.5\lambda$ 是标准。不过在Sub-6GHz频段，由于单元尺寸的物理限制，容易维持 $d = 0.5\lambda$，而在毫米波频段，单元间的制造公差会成为问题。
单元间互耦
老师一开始说"互耦很重要"，具体影响有多大呢？
用数字来说明比较容易理解。例如，对于 $d = 0.5\lambda$ 的贴片天线阵列，比较忽略互耦和考虑互耦的情况：
- 旁瓣电平(SLL)：忽略 → -26dB / 考虑 → -18dB（恶化8dB）
- 零陷深度：忽略 → -60dB / 考虑 → -30dB（零陷变浅）
- 主波束方向：最大1.5°的指向误差
- 输入阻抗：因单元位置最大有20%的变动（边缘单元 vs 中央单元）
居然差8dB！考虑互耦的分析要怎么做呢？
主要有三种方法：
- 全单元全波分析：用FEM或MoM一次性求解所有单元。最精确，但对于128单元，计算成本巨大
- 嵌入单元方向图（Embedded Element Pattern）：只激励阵列中的一个单元，其他单元接匹配终端，获取其方向图。对N个单元重复此操作。计算成本和精度平衡较好
- 耦合矩阵校正：测量/计算S参数的互耦分量 $S_{ij}$，用于校正权重向量的方法。$\mathbf{w}_{\mathrm{corrected}} = \mathbf{C}^{-1}\mathbf{w}_{\mathrm{ideal}}$
波束成形与"鸡尾酒会效应"
在嘈杂的派对上能只听清特定人声音的"鸡尾酒会效应"，是大脑利用双耳信号的到达时间差（ITD）来估计方向。波束成形的数学本质与此完全相同，给每个单元的信号施加适当的相位延迟后叠加，特定方向的信号就会被增强。人类的听觉只有2个单元（左右耳）的阵列，而5G基站有128个单元。拥有128只耳朵的超人在派对上能分别听清所有人的对话——这就是大规模MIMO波束成形的样子。
波束成形天线阵列的电磁的数值方法
全波分析方法
分析天线阵列的电磁特性时，听说有FEM、FDTD等各种方法，波束成形用哪种呢？
主要有四种方法，根据阵列规模和所需精度来选用：
| 方法 | 原理 | 擅长的问题 | 对大规模阵列的适用性 |
|---|---|---|---|
| FEM（有限元法） | 弱形式 + 边单元（Nedelec单元） | 复杂形状、介质基板、多层结构 | 自适应网格，适用于中等规模 |
| FDTD（时域有限差分法） | 麦克斯韦方程组的Yee网格离散化 | 宽带特性、瞬态响应 | 内存需求大，大规模困难 |
| MoM（矩量法） | 积分方程 + 基函数展开 | 仅金属结构、开放区域问题 | 稠密矩阵$O(N^2)$，可用FMM加速 |
| MLFMM（多层快速多极子法） | 将MoM加速至$O(N\log N)$ | 大规模金属阵列 | 实用至数千单元 |
对于5G基站的128单元贴片阵列，现实中选择哪种方法？
现实中的主流是混合方法："用FEM求解单个单元+周边少数单元，获取嵌入单元方向图 → 用阵列因子合成"。如果尝试用FEM求解全部128个单元，在28GHz下网格会达到数亿单元，即使是128GB内存的工作站也很困难。
区域分解与混合法
混合法是什么机制呢？
思路很简单。近场区域（单元及其邻近环境）用FEM精确求解，远场区域用MoM或FMM高效处理。具体方法有：
- FEM-BI（边界积分）法：在FEM区域的外边界上放置等效电磁流，用积分方程求解外部。Ansys HFSS的Domain Decomposition法就是这种
- FE-BI + MLFMM：用MLFMM加速BI部分。FEKO（Altair）擅长此方法
- 特征基函数法：以子阵列为单位生成基函数，减少整体问题的自由度
计算成本会变化多少呢？
以64单元贴片阵列（28GHz）为例进行比较：
| 方法 | 自由度数 | 内存 | 计算时间（16核） |
|---|---|---|---|
| 全单元FEM | ~2亿 | ~200GB | ~48小时 |
| FEM-BI（区域分解） | ~3000万/子域 | ~32GB | ~6小时 |
| 嵌入单元方向图 + AF合成 | ~500万（单单元+邻近） | ~8GB | ~30分钟 |
有两位数的差距。所以实际工作中嵌入单元方向图法占绝大多数。
嵌入式单元方向图法
请告诉我嵌入单元方向图的具体步骤！
步骤如下：
- 创建包含全部N个单元的阵列模型
- 仅以1W激励第$m$单元，其他所有单元接匹配终端（50Ω负载）
- 在此状态下获取远场方向图 $\mathbf{E}_m(\theta, \phi)$——这就是第$m$单元的"嵌入单元方向图"
- 对 $m = 0, 1, \ldots, N-1$ 的所有单元重复此操作
- 叠加整体方向图：
其中 $w_m$ 是第$m$单元的幅度权重，$\beta_m$ 是相位控制。这样就把互耦的影响自动包括在内了。计算效率比全元素求解高得多。
为什么其他单元要接匹配终端而不是开路？
好问题。开路（高阻）会产生反射，导致阵列中的驻波和失配。匹配终端（50Ω）表示这些单元吸收入射信号，模拟实际馈电网络中的情况。这样获取的方向图才能准确反映带宽工作时的互耦效应。