波束成形天线阵列的电磁分析
理论与物理
波束成形是什么
波束成形在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波束成形的样子。
阵列因子的推导与物理意义
- 相位项 $k \cdot d \cdot n \cdot \sin\theta$：从第$n$单元向角度$\theta$方向辐射时的空间相位延迟。$k = 2\pi/\lambda$ 是自由空间波数，$d \cdot \sin\theta$ 是相邻单元间的路径差。
- 相移 $\beta_n$：通过移相器人为施加给各单元的相位。用于波束转向或零陷转向。模拟移相器通常进行4~6比特的量化，相位量化误差会劣化旁瓣。
- 幅度权重 $w_n$：各单元的幅度控制。等幅度时旁瓣为-13.2dB，使用泰勒窗或切比雪夫窗可改善至-30~-40dB。但存在波束宽度变宽的权衡。
- SINR：$\mathrm{SINR} = \frac{|\mathbf{w}^H \mathbf{a}(\theta_s)|^2}{\mathbf{w}^H \mathbf{R}_{i+n} \mathbf{w}}$。权重向量$\mathbf{w}$、转向向量$\mathbf{a}(\theta_s)$、干扰+噪声的协方差矩阵$\mathbf{R}_{i+n}$。MVDR波束成形器就是求解使此SINR最大化的权重。
假设条件与适用范围
- 远场假设：在 $R > 2D^2/\lambda$（$D$为阵列口径）条件下AF有效。近距离需要菲涅耳区域的相位校正
- 忽略互耦（Pattern Multiplication）：假设所有单元方向图相同。$d > \lambda$ 时大致有效，但 $d < 0.5\lambda$ 时失效
- 窄带假设：带宽小于中心频率的5%时，相位控制足够。宽带时需要真时延（TTD）控制
- 平面波入射：假设信号源足够远。室内或NFC应用需要近距离模型
波束成形主要参数一览
| 参数 | 符号 | 单位 | 典型值（28GHz, 64单元ULA） |
|---|---|---|---|
| 波数 | $k$ | rad/m | $\approx 586$ |
| 单元间距 | $d$ | mm | 5.35（= λ/2） |
| HPBW | — | 度 | $\approx 1.6°$ |
| 阵列增益 | $G_A$ | dBi | $\approx 10\log_{10}(64) = 18\,\mathrm{dBi}$ |
| 第一旁瓣 | SLL | dB | -13.2（等幅度）/ -30（泰勒窗） |
数值解法与实现
全波解析方法
分析天线阵列的电磁特性时，听说有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$ 的所有单元重复此操作
- 叠加整体方向图：
なった
詳しく
報告