2. 阵列信号处理基础
阵列信号处理是数字波束形成与混合波束系统的数学地基。无论是毫米波相控阵、大规模MIMO基站，还是车载雷达，其核心都归结为同一件事：用多个空间上分离的传感器采样同一波场，通过相位差反推信号来向，并通过加权合成把能量聚焦到期望方向。本章围绕均匀线阵（ULA）、均匀平面阵（UPA）、导向矢量、阵列流形、空间频率与波达方向（DOA）六个基本概念展开，建立"几何—相位—数学表示"的完整链条。
一、核心原理：几何决定相位，相位编码方向
假设远场窄带信号以平面波形式入射，波长为 λ。当平面波到达相邻两个阵元时，由于传播路径不同，会产生一个固定的路程差 Δd。这个路程差在载波上表现为相位差：
Δφ = 2π · Δd / λ = 2π · d · sinθ / λ其中 d 为阵元间距，θ 为入射角（以阵列法线为 0°）。这个式子是整个阵列处理的"第一性原理"：相位差是方向的函数。反过来，只要测出阵元间的相位差，就能估计 θ。这就是DOA估计的物理基础。
把 N 个阵元的相位差按顺序排列，就得到导向矢量（Steering Vector）。对均匀线阵：
a(θ) = [1, e^(j·2πd·sinθ/λ), ..., e^(j·2πd·(N-1)·sinθ/λ)]^T导向矢量是"某个方向在阵列上的签名"。不同方向对应不同签名，因此阵列可以像指纹识别一样区分来向。所有可能方向的导向矢量集合，称为阵列流形（Array Manifold），它是 N 维复空间中的一条曲线（线阵）或一张曲面（面阵）。
二、均匀线阵与均匀平面阵的结构对比
均匀线阵（ULA）是最简单的阵列：N 个阵元等间距排成一条直线。它只能分辨一维角度（方位或俯仰之一），且存在前后模糊和栅瓣问题。均匀平面阵（UPA）在平面内按 M×N 网格排列，可同时估计方位角和俯仰角，是二维波束形成和混合波束系统的常用结构。
下面这张图展示了 ULA 与 UPA 的几何结构、入射平面波、阵元间相位差以及导向矢量的对应关系，是本节最核心的可视化对象。
三、空间频率：把角度问题变成频率问题
观察导向矢量 a(θ) 的指数项，可以定义空间频率：
u = d · sinθ / λ于是 a(θ) = [1, e^(j2πu), e^(j4πu), …, e^(j2π(N-1)u)]ᵀ。这与时间域离散信号的复指数形式完全一致，只是"频率"变成了"空间频率"。这一类比极其重要：所有时间域频谱分析工具（FFT、加窗、谱估计）都可以平移到空间域。对 ULA 做空间 FFT，就能得到波束扫描谱，峰值位置对应 DOA。
空间频率的取值范围是 [-0.5, 0.5]（对应 sinθ ∈ [-1, 1] 且 d = λ/2）。当 d > λ/2 时，u 会超出该范围，产生栅瓣（grating lobe），即多个角度对应同一空间频率，导致方向模糊。这就是为什么工程上通常取 d = λ/2。
四、阵列流形与DOA估计流程
阵列流形 A = [a(θ₁), a(θ₂), …, a(θ_K)] 把所有方向的导向矢量堆叠成矩阵。接收信号模型为：
X(t) = A · S(t) + N(t)其中 S(t) 为信源信号，N(t) 为噪声。DOA 估计就是从 X(t) 中恢复 θ。典型流程如下：
- 采集：N 个阵元同步采样，得到 N×T 数据矩阵。
- 协方差估计：R = E[XXᴴ]，实际用有限快拍平均。
- 谱构造：MUSIC 用噪声子空间，Capon 用 R⁻¹，构造空间谱 P(θ)。
- 谱峰搜索：在 θ ∈ [-90°, 90°] 上找 P(θ) 的峰值。
- 角度输出：峰值位置即 DOA 估计值。
五、关键参数与常见误区
| 参数/概念 | 数学表示 | 物理含义 | 工程取值/注意 |
|---|---|---|---|
| 阵元间距 d | d | 相邻阵元空间距离 | 通常 d = λ/2，避免栅瓣 |
| 空间频率 u | u = d·sinθ/λ | 角度在阵列上的归一化相位斜率 | |u| ≤ 0.5 无模糊 |
| 导向矢量 a(θ) | N×1 复向量 | 阵列对方向 θ 的响应 | 与阵列几何强相关 |
| 阵列流形 | {a(θ) | θ∈Θ} | 所有方向响应的集合 | N 维复空间中的曲线/曲面 |
| 波达方向 θ | θ | 信号入射角（法线为 0°） | ULA 仅一维，UPA 二维 |
| 栅瓣 | d > λ/2 时出现 | 虚假谱峰，方向模糊 | 必须通过 d ≤ λ/2 抑制 |
六、案例：8 元 ULA 的 DOA 估计
设 8 元 ULA，d = λ/2，信源从 θ = 30° 入射。则空间频率 u = 0.5·sin30° = 0.25，相邻阵元相位差 Δφ = 2π×0.25 = π/2。导向矢量为：
a(30°) = [1, e^(jπ/2), e^(jπ), e^(j3π/2), e^(j2π), e^(j5π/2), e^(j3π), e^(j7π/2)]ᵀ
= [1, j, -1, -j, 1, j, -1, -j]ᵀ可以看到相位以 π/2 为步进循环，这正是空间频率 0.25 的离散复指数特征。若用 MUSIC 算法，噪声子空间与 a(30°) 正交，空间谱在 30° 处出现尖峰。若把 d 改为 λ，则 u = 0.5，相位差变为 π，此时 θ = 30° 与 θ = -30° 产生相同空间频率，出现方向模糊——这就是栅瓣的直观体现。
掌握 ULA、UPA、导向矢量、阵列流形、空间频率与 DOA 这六个概念及其数学关系，是后续理解波束形成权重设计、混合波束架构和自适应算法的前提。它们共同构成了阵列信号处理的"坐标系"：几何定义相位，相位定义空间频率，空间频率反解方向，而阵列流形则把所有可能方向组织成一个可搜索的数学空间。