3、波束赋形核心理论：加权向量、幅度锥削、相位补偿、波束形成公式与方向图计算

相控阵天线的本质，是通过控制阵列中每个单元的激励幅度与相位，使各单元辐射场在空间指定方向同相叠加，在其他方向相互抵消，从而形成具有指向性的波束。这一过程在数学上可抽象为加权向量对阵列流形向量的线性组合。本章从加权向量的物理含义出发，依次推导幅度锥削、相位补偿与波束形成公式，并给出方向图计算与可视化方法，为后续版图参数对齐提供理论基准。

3.1 加权向量的数学定义与物理意义

相控阵天线由 N 个辐射单元组成，第 n 个单元在远场观察方向 θ 处的辐射贡献可表示为：

Eₙ(θ) = wₙ · e^{j·k·dₙ·sinθ}

其中，wₙ 为第 n 个单元的复加权系数，k = 2π/λ 为自由空间波数，dₙ 为第 n 个单元相对于参考点的位置坐标。将所有单元的贡献求和，即得到阵列因子（Array Factor, AF）：

AF(θ) = Σₙ wₙ · e^{j·k·dₙ·sinθ}   (n = 1, 2, ..., N)

加权向量 w = [w₁, w₂, ..., wₙ]ᵀ 完整描述了阵列的激励状态。每个加权系数可分解为幅度项与相位项：

wₙ = aₙ · e^{jφₙ}

其中 aₙ 为幅度激励（通常取正值），φₙ 为相位激励。幅度项决定各单元对波束贡献的“大小”，相位项决定各单元辐射场的“时间对齐”程度。加权向量的设计直接决定了波束的主瓣指向、主瓣宽度、旁瓣电平以及零陷位置等核心指标。

3.2 幅度锥削：窗函数的选择与旁瓣抑制

均匀激励（所有 aₙ = 1）时，阵列方向图的第一旁瓣电平约为 -13.3 dB，这在许多雷达和通信系统中无法满足要求。幅度锥削通过调整各单元的幅度激励，使阵列口径场分布从中心向边缘逐渐衰减，从而降低旁瓣电平。

常用的幅度锥削窗函数包括：

幅度锥削的本质是用主瓣展宽和增益损失换取旁瓣抑制。设计时需要根据系统指标权衡：例如，某相控阵要求旁瓣低于 -35 dB，则 Taylor 窗（n̄=4）通常比 Hamming 窗更合适，因为其主瓣展宽更小、增益损失更低。

3.3 相位补偿：波束指向控制的核心

波束指向控制的核心思想是：通过调整各单元的相位，使各单元在期望方向 θ₀ 上的辐射场同相叠加。第 n 个单元在 θ₀ 方向的空间相位延迟为 k·dₙ·sinθ₀，因此补偿相位应为：

φₙ = -k·dₙ·sinθ₀

代入波束形成公式后，阵列因子变为：

AF(θ) = Σₙ aₙ · e^{j·k·dₙ·(sinθ - sinθ₀)}

当 θ = θ₀ 时，所有单元的相位项均为零，辐射场同相叠加，形成主瓣峰值。当 θ ≠ θ₀ 时，各单元相位不一致，部分抵消，形成旁瓣或零陷。

对于均匀线阵（单元间距 d），波束指向 θ₀ 时的相位补偿量呈等差数列分布：

φₙ = -n · (2πd/λ) · sinθ₀   (n = 0, 1, ..., N-1)

这一线性相位分布是相控阵波束扫描的基础。实际工程中，移相器通常为数字移相器（如 5 位、6 位），相位量化误差会导致波束指向偏差和旁瓣抬升，需要在版图参数对齐时予以考虑。

3.4 波束形成公式的完整推导与方向图计算流程

综合幅度锥削与相位补偿，完整的波束形成公式为：

AF(θ) = Σₙ aₙ · e^{j(φₙ + k·dₙ·sinθ)}
       = Σₙ aₙ · e^{j·k·dₙ·(sinθ - sinθ₀)}

方向图计算通常采用以下流程：

1. 确定阵列参数：单元数 N、单元间距 d、工作频率 f（对应波长 λ）；
2. 选择幅度锥削窗函数：根据旁瓣指标确定 aₙ 序列；
3. 计算相位补偿：根据期望波束指向 θ₀ 计算 φₙ = -k·dₙ·sinθ₀；
4. 构建加权向量：wₙ = aₙ · e^{jφₙ}；
5. 扫描观察角 θ：在 [-90°, 90°] 范围内以 0.1° 步进计算 AF(θ)；
6. 归一化与 dB 转换：以最大值归一化，计算 20·log₁₀|AF(θ)|；
7. 提取关键指标：主瓣峰值位置、3dB 波束宽度、第一旁瓣电平、零陷位置。

以下为 Python 计算示例（使用 NumPy）：

import numpy as np

def array_factor(theta_deg, theta0_deg, N, d_over_lambda, window=None):
    theta = np.deg2rad(theta_deg)
    theta0 = np.deg2rad(theta0_deg)
    n = np.arange(N)
    # 相位补偿
    phi = -2 * np.pi * d_over_lambda * n * np.sin(theta0)
    # 幅度锥削
    if window is None:
        a = np.ones(N)
    elif window == 'hamming':
        a = np.hamming(N)
    elif window == 'taylor':
        # 简化 Taylor 窗（n_bar=4, SLL=-35dB）
        a = np.sinc(n / N) * np.hamming(N)  # 近似
    # 加权向量
    w = a * np.exp(1j * phi)
    # 阵列因子
    AF = np.sum(w[:, None] * np.exp(1j * 2 * np.pi * d_over_lambda * n[:, None] * np.sin(theta)), axis=0)
    return 20 * np.log10(np.abs(AF) / np.max(np.abs(AF)) + 1e-12)

3.5 案例拆解：16 单元线阵的波束赋形设计与方向图分析

背景与约束：某车载雷达相控阵天线，工作频率 24 GHz（λ = 12.5 mm），采用 16 单元微带贴片线阵，单元间距 d = 7.5 mm（0.6λ）。系统要求：波束扫描范围 ±30°，旁瓣电平 ≤ -30 dB，3dB 波束宽度 ≤ 8°。

初始设计：采用均匀激励，波束指向 θ₀ = 0°（法向）。计算得到方向图指标：

分析判断：均匀激励旁瓣电平 -13.3 dB 远高于 -30 dB 指标，必须采用幅度锥削。Hamming 窗虽然旁瓣最低，但主瓣展宽至 8.5°，超过 8° 的指标限制。Taylor 窗（n̄=4）在旁瓣 -35 dB 与主瓣 7.5° 之间取得最佳平衡，满足全部指标。

实际操作步骤：

1. 在 MATLAB/Python 中生成 Taylor 窗系数（n̄=4，SLL=-35 dB），得到 16 个幅度值 aₙ；
2. 计算波束指向 θ₀ = 30° 时的相位补偿 φₙ = -n·(2π·0.6)·sin30° = -n·1.885 rad；
3. 构建加权向量 wₙ = aₙ·e^{jφₙ}，代入阵列因子公式计算方向图；
4. 验证扫描至 30° 时主瓣峰值位置、旁瓣电平及波束宽度变化；
5. 将加权向量转换为版图参数：幅度 → 衰减器 dB 值，相位 → 移相器数字码。

结果对比：采用 Taylor 窗后，法向方向图第一旁瓣为 -35.2 dB，3dB 波束宽度 7.5°，满足指标。扫描至 30° 时，由于阵列投影口径减小，波束宽度展宽至 8.7°，旁瓣略有抬升至 -33.5 dB，仍在可接受范围内。

风险与代价：Taylor 窗导致增益损失 0.7 dB，意味着雷达探测距离缩短约 3.5%。此外，边缘单元的幅度激励仅为中心单元的 0.45 倍，对衰减器精度要求较高（±0.2 dB 以内），否则实际旁瓣电平可能恶化 2~3 dB。

复盘经验：幅度锥削设计必须在旁瓣、波束宽度、增益三者之间权衡。对于车载雷达这类对探测距离敏感的系统，可考虑“部分锥削”方案——仅对边缘 4 个单元进行锥削，中心 12 个单元保持均匀激励，可在旁瓣 -28 dB 与增益损失 0.3 dB 之间取得折中。

3.6 常见误区与工程复盘

误区一：认为幅度锥削越强越好。实际上，过度锥削会导致主瓣过度展宽、增益严重下降，且对幅度误差更敏感。当旁瓣指标为 -30 dB 时，Taylor 窗（n̄=4）通常比 Hamming 窗更合适，因为后者主瓣展宽过大。

误区二：忽略相位量化误差的影响。数字移相器的量化误差（如 5 位移相器步进 11.25°）会导致波束指向偏差和旁瓣抬升。对于 16 单元阵列，5 位量化引起的最大指向偏差约为 0.5°，旁瓣抬升约 1~2 dB。在版图参数对齐时，必须将量化误差纳入预算。

误区三：将加权向量直接映射到版图参数，未考虑单元互耦。实际阵列中，单元间的互耦会改变有效激励幅度和相位，导致方向图畸变。工程上应在版图仿真（如 HFSS、CST）中提取实际单元方向图，再与理论方向图对比，通过迭代优化补偿互耦影响。