阵列天线原理：阵列因子、方向图乘积原理、均匀线性阵列
单个辐射单元的方向图是有限的，但工程中我们往往需要高增益、窄波束或可控波束扫描的天线。阵列天线通过将多个单元按一定规律排列，利用电磁波在空间中的相干叠加，实现对辐射方向图的“重塑”。理解这一过程的核心，在于掌握阵列因子与方向图乘积原理。本节以最基础的均匀线性阵列（ULA）为对象，拆解其数学本质与物理意义。
3.1 方向图乘积原理：单元因子与阵列因子
设阵列由 \( N \) 个相同的辐射单元组成，第 \( n \) 个单元在远区观察点 \( (\theta, \phi) \) 处产生的电场可表示为单元方向图 \( f(\theta, \phi) \) 与激励复振幅 \( I_n e^{j\psi_n} \) 的乘积，再乘上由单元空间位置引起的相位延迟 \( e^{j k \hat{r} \cdot \vec{r}_n } \)。总场为所有单元贡献的叠加：
\[ \vec{E}_{\text{total}}(\theta,\phi) = f(\theta,\phi) \cdot \sum_{n=0}^{N-1} I_n e^{j\psi_n} e^{j k \hat{r} \cdot \vec{r}_n} \]观察上式，总方向图被分解为两个独立部分的乘积：
- 单元因子（Element Factor）：\( f(\theta,\phi) \)，仅取决于单个辐射单元的类型与尺寸，与阵列布局无关。
- 阵列因子（Array Factor, AF）：\( \sum_{n=0}^{N-1} I_n e^{j\psi_n} e^{j k \hat{r} \cdot \vec{r}_n} \)，仅取决于单元间的相对位置、激励幅度与相位，与单元本身无关。
下图展示了方向图乘积原理的直观含义：一个半波偶极子（单元因子为面包圈形）沿 x 轴排列成 8 元均匀线性阵列（阵列因子为多波瓣形），两者相乘后得到实际总方向图——在阵列因子主瓣方向保留单元辐射，在阵列因子零点方向总场被抑制。
3.2 均匀线性阵列的阵列因子推导
均匀线性阵列指 \( N \) 个单元等间距 \( d \) 排列在一条直线上，各单元激励幅度相同（\( I_n = I_0 \)），激励相位按等差级数递增（\( \psi_n = n\alpha \)）。设阵列沿 z 轴放置，第一个单元位于原点，则第 \( n \) 个单元的位置为 \( z_n = nd \)。
在远区，观察方向与阵列轴（z 轴）夹角为 \( \theta \)，相邻单元到观察点的波程差为 \( d\cos\theta \)，对应的空间相位差为 \( k d \cos\theta \)。阵列因子为等比级数求和：
\[ AF(\theta) = \sum_{n=0}^{N-1} e^{jn(\alpha + k d \cos\theta)} = \frac{1 - e^{jN\psi}}{1 - e^{j\psi}}, \quad \psi = \alpha + k d \cos\theta \]利用欧拉公式化简，取幅度绝对值，得到常用的阵因子方向函数：
\[ |AF(\theta)| = \left| \frac{\sin\left( \frac{N\psi}{2} \right)}{\sin\left( \frac{\psi}{2} \right)} \right| \]当 \( \psi = 0 \) 时，分子分母同时趋于零，由洛必达法则可得峰值 \( |AF|_{\max} = N \)。归一化后，阵列因子可写为：
\[ AF_{\text{norm}}(\theta) = \frac{\sin\left( \frac{N}{2} (k d \cos\theta + \alpha) \right)}{N \sin\left( \frac{1}{2} (k d \cos\theta + \alpha) \right)} \]- ψ —— 相邻单元之间的总相位差（激励相位差 + 空间波程差），是分析阵列因子的核心变量。
- α —— 馈电相位递增量（度/单元），控制波束扫描方向。
- d —— 单元间距，决定栅瓣出现条件。
- N —— 单元数量，决定波束宽度与副瓣电平。
3.3 均匀线性阵列的方向图特性
基于上述公式，我们可以系统分析 ULA 方向图的几个核心指标。下表总结了各参数对方向图的影响规律：
| 参数 | 主瓣方向 θ₀ | 半功率波束宽度 HPBW | 第一副瓣电平 SLL | 栅瓣条件 |
|---|---|---|---|---|
| 单元数 N 增大 | 不变 | 变窄（≈ 0.886λ/(Nd)） | 基本不变（约 -13.5 dB） | 无影响 |
| 间距 d 增大 | 不变 | 变窄 | 不变 | d ≥ λ 时出现栅瓣 |
| 相位递增量 α 增大 | 向端射方向偏移 | 略有展宽 | 可能恶化 | 扫描角增大时更易出现 |
主瓣方向由 \( \psi = 0 \) 条件决定：
\[ \cos\theta_0 = -\frac{\alpha}{k d} = -\frac{\alpha \lambda}{2\pi d} \]当 \( \alpha = 0 \) 时，主瓣指向阵列法线方向（\( \theta_0 = 90^\circ \)），称为侧射阵；当 \( \alpha = -k d \) 时，主瓣指向阵列轴线方向（\( \theta_0 = 0^\circ \)），称为端射阵。通过连续调节 \( \alpha \)，可实现波束在空间的电扫描。
3.4 仿真案例：8 元 ULA 方向图计算
下面通过一个 Python 脚本，使用 NumPy 计算并绘制 8 元 ULA 在侧射与扫描两种状态下的归一化方向图。该脚本可直接运行，用于验证本节公式。
import numpy as np
import matplotlib.pyplot as plt
N = 8 # 单元数
d = 0.5 # 间距（波长单位）
theta = np.linspace(0, np.pi, 1801) # 观察角 0~180°
def array_factor(N, d, alpha, theta):
"""计算归一化阵列因子（线性刻度）"""
kd = 2 * np.pi * d
psi = alpha + kd * np.cos(theta)
# 避免除零：psi 接近 0 时取极限值 N
AF = np.sin(N * psi / 2) / (N * np.sin(psi / 2 + 1e-12))
return np.abs(AF)
# 侧射阵：alpha = 0
AF_broadside = array_factor(N, d, 0.0, theta)
# 扫描阵：主瓣指向 60°（相对阵列轴）
theta0 = np.deg2rad(60)
alpha_scan = -2 * np.pi * d * np.cos(theta0)
AF_scan = array_factor(N, d, alpha_scan, theta)
# 转换为 dB
AF_db1 = 20 * np.log10(AF_broadside + 1e-6)
AF_db2 = 20 * np.log10(AF_scan + 1e-6)
plt.figure(figsize=(8, 4))
plt.plot(np.rad2deg(theta), AF_db1, label='Broadside (α=0)')
plt.plot(np.rad2deg(theta), AF_db2, '--', label='Scan to 60°')
plt.xlabel('θ (deg)'); plt.ylabel('Normalized AF (dB)')
plt.ylim(-40, 5); plt.grid(True); plt.legend()
plt.title(f'{N}-element ULA, d=λ/2')
plt.tight_layout(); plt.show()运行上述代码，可以观察到：侧射时主瓣在 90° 方向，第一副瓣约为 -13.5 dB；当 α 设置为扫描到 60° 时，主瓣偏移至 60° 方向，但波束有所展宽，且一侧副瓣电平略有升高。这验证了表 3.1 中的规律。
3.5 常见误区与设计复盘
- 先明确工作频率与扫描范围，确定最大允许间距 \( d_{\max} = \lambda / (1 + |\cos\theta_{\max}|) \)。
- 根据波束宽度要求确定单元数 \( N \approx 0.886\lambda / (\text{HPBW} \cdot d) \)。
- 在仿真软件中分别绘制单元方向图与阵列因子，验证乘积原理。
- 检查扫描到最大角度时的副瓣抬升与阻抗失配情况，必要时引入幅度加权（如泰勒加权）压低副瓣。
均匀线性阵列是阵列天线中最基础也最重要的构型。通过方向图乘积原理，我们将复杂的阵列辐射问题拆解为“单元设计”与“阵列综合”两个相对独立的子问题；通过阵列因子的闭式表达式，我们能够快速预测波束指向、波束宽度与栅瓣位置。这些概念将直接推广到平面阵列、共形阵列以及数字波束形成等更复杂的阵列系统中。