阵列因子推导：从各向同性单元推导阵列因子的一般表达式

在阵列天线分析中，我们面临的核心问题是：当多个辐射单元按一定规律排列时，远场总辐射方向图如何由单元方向图与阵列布局共同决定？本节将从最简单的各向同性单元出发，严格推导阵列因子（Array Factor, AF）的一般表达式。这一推导不仅是后续分析均匀直线阵、平面阵的基础，更是理解波束扫描、栅瓣抑制等工程问题的数学根基。

核心原理：方向图乘积定理与相位叠加

阵列天线的远场总场强可以表示为单元方向图与阵列因子的乘积，即方向图乘积定理：总方向图 = 单元方向图 × 阵列因子。其中阵列因子仅取决于单元的排列方式、激励幅度与相位，与单元本身无关。因此，采用各向同性单元（方向图为球面均匀辐射）作为推导起点，可以剥离单元影响，纯粹研究阵列几何与激励对辐射方向图的作用。

考虑N个各向同性点源沿任意空间位置分布，第n个单元位于坐标矢量 rn，其激励复电流为 In = Anejαn（A_n为幅度，α_n为相位）。在远场观察点方向 r̂ 处，第n个单元相对于坐标原点的相位延迟为：

Δψn = k · (r̂ · rn) = (2π/λ) · (xncosφsinθ + ynsinφsinθ + zncosθ)

其中k为波数，λ为波长，θ为俯仰角（从z轴起算），φ为方位角。总远场电场为各单元贡献的叠加：

Etotal(θ,φ) = Σn=1N In · ejk(r̂·rn) / (4πr)

提取公因子后，阵列因子定义为：

AF(θ,φ) = Σn=1N Anej(αn + k·r̂·rn)

这就是阵列因子的一般表达式。它本质上是各单元激励的相位加权叠加，其中包含两项相位：一项是激励自身相位α_n，另一项是空间传播路径差引入的相位 k·r̂·r_n。

方法流程：从一般到特殊的推导路径

为了得到工程中常用的显式表达式，我们按以下步骤推进推导：

1. 建立坐标系与位置矢量：将阵列几何参数化。对于直线阵，令单元沿x轴等间距d排列，则 r_n = (n-1)d · x̂（n=1,2,…,N）。
2. 计算空间相位差：对于沿x轴的直线阵，观察方向与x轴的夹角余弦为 cosψ = sinθcosφ。因此空间相位差为 k·r̂·r_n = k(n-1)d·sinθcosφ。
3. 引入均匀激励条件：设各单元幅度相等 A_n=A₀，相位按等差递增 α_n = (n-1)β，其中β为相邻单元间的激励相位差（由移相器提供）。
4. 化简为等比级数：将上述条件代入一般表达式，得到等比级数求和形式。
5. 提取幅度归一化形式：利用欧拉公式将复数求和化为正弦比形式，得到经典的均匀直线阵阵列因子。

对于均匀直线阵，令 u = kd·sinθcosφ + β，则阵列因子可写为：

AF(u) = A₀ · Σn=0N-1 ejnu = A₀ · [sin(Nu/2) / sin(u/2)] · ej(N-1)u/2

若以阵列中心为相位参考点，则相位因子 e^j(N-1)u/2 消失，得到实函数：

AF(θ) = A₀ · sin[N·(kd·sinθ + β)/2] / sin[(kd·sinθ + β)/2]

这一表达式是后续分析波束指向、主瓣宽度、副瓣电平的基础。

案例/场景：8单元均匀直线阵的阵列因子计算

设有一个8单元均匀直线阵，工作频率 f = 10 GHz（λ = 30 mm），单元间距 d = 15 mm（即0.5λ），各单元等幅激励，相邻相位差 β = 0（边射阵）。我们计算其阵列因子并分析方向图特征。

第一步：确定参数

N = 8, d = 0.5λ, β = 0
u = kd·sinθ = (2π/λ)·(0.5λ)·sinθ = π·sinθ

第二步：代入阵列因子公式

AF(θ) = sin(8·π·sinθ/2) / sin(π·sinθ/2)
      = sin(4π·sinθ) / sin(π·sinθ/2)

第三步：求主瓣方向：当 u=0 时，即 θ=0°，AF取最大值N=8。主瓣指向阵列法线方向。

第四步：求零陷位置：令分子 sin(4π·sinθ)=0 且分母不为0，得 4π·sinθ = mπ（m=1,2,3,4），即 sinθ = m/4。对应 θ = 14.48°, 30°, 48.59°, 90°。

第五步：求副瓣峰值：副瓣出现在分子接近±1处，第一副瓣电平约为 -13.3 dB（相对于主瓣）。

下表总结了关键参数与计算结果：

下图展示了阵列因子幅度随角度变化的极坐标方向图，可以直观看到主瓣、副瓣和零陷的分布规律：