2、阵列信号基础:阵列几何构型、阵列流形向量与窄带信号模型

阵列信号处理的核心任务,是从多个空间分布的传感器(阵元)同步采集到的信号中,提取信源的方向、幅度、频率等信息。要完成这一任务,首先必须建立一套严谨的数学模型,将物理世界的“空间位置”与数学世界的“相位差”对应起来。本章作为全课程的基石,将系统梳理三种典型的阵列几何构型(均匀线阵ULA、均匀面阵URA、任意阵列),推导描述阵列空间响应的核心工具——阵列流形向量(Array Manifold Vector),并建立窄带信号模型。理解这三者之间的内在联系,是后续学习波束形成权值计算、DOA估计算法(如MUSIC、ESPRIT)的前提。

阵列信号基础:几何构型 → 流形向量 → 窄带模型 均匀线阵 (ULA) d=λ/2 平面波前 (θ) 相位差: 2πd·sinθ/λ 均匀面阵 (URA) Y轴 X轴 dx=dy=λ/2 任意阵列 X Y 位置矢量: p = [x,y,z]ᵀ 阵列流形向量 a(θ) 与窄带信号模型 a(θ) = [1, e^{-j2πd·sinθ/λ}, ..., e^{-j2π(M-1)d·sinθ/λ}]ᵀ M: 阵元数 | d: 阵元间距 | θ: 入射角 相位差 = 2π·(空间程差)/λ 空间程差 = d·sinθ (ULA) x(t) = a(θ)·s(t) + n(t) x(t): M×1 接收快拍 | s(t): 信号复包络 n(t): 噪声向量 | 窄带: 包络变化远慢于载波 带宽 B << 载频 f₀, 且 B·(阵列孔径)/c << 1

上图从整体上勾勒了本章的知识框架:左侧的几何构型决定了阵元间的空间相位关系,中间的阵列流形向量是该关系的数学封装,右侧的窄带模型则将这种封装嵌入到接收信号的时域表达中。下面我们逐一展开。

2.1 阵列几何构型:从一维到三维的空间采样

核心原理:阵列几何构型决定了空间采样方式,直接影响阵列的孔径、分辨能力和模糊度。不同构型在相同阵元数下,其方向图主瓣宽度、旁瓣电平以及DOA估计的CRB(克拉美-罗界)均有显著差异。

  • 均匀线阵(ULA):所有阵元等间距排列在一条直线上。设阵元间距为d,阵元数为M,则第m个阵元的位置为pm = (m-1)·d。ULA只能区分一维角度(方位角θ),且存在前后向模糊(θ与π-θ不可区分)。工程中通常取d = λ/2,以避免栅瓣(Grating Lobe)。
  • 均匀面阵(URA):阵元在二维平面上按矩形或三角形网格排列。URA可同时估计方位角θ和俯仰角φ,消除了ULA的前后向模糊。设x方向间距为dx,y方向间距为dy,则第(m,n)个阵元位置为pm,n = ((m-1)dx, (n-1)dy)。
  • 任意阵列:阵元位置无规则约束,可以是共形阵(贴附于载体表面)、稀疏阵或随机阵。任意阵列的流形向量无法用闭式表达式统一描述,通常需要实测或数值计算。其优势在于可灵活设计波束形状,但代价是算法复杂度显著上升。
构型类型维度可估计角度流形向量表达式典型应用
ULA1D方位角 θ闭式解析式雷达测向、通信基站
URA2D方位角 θ + 俯仰角 φ闭式解析式(Kronecker积)相控阵雷达、卫星通信
任意阵列2D/3Dθ, φ(或更多)数值表或插值声纳、共形天线、无人机集群
经验提醒:在系统设计初期,不要盲目追求复杂构型。若目标只测方位角,ULA是最优选择;若需要俯仰角且安装空间允许,URA是标准方案;任意阵列通常用于特殊载体(如导弹导引头、机翼蒙皮),此时必须配合高精度标定。

2.2 阵列流形向量:空间响应的数学指纹

核心原理:阵列流形向量a(θ)描述了来自方向θ的单位功率窄带信号在M个阵元上引起的复振幅响应。其本质是各阵元相对于参考阵元的相位差向量。对于ULA,以第0个阵元为参考,第m个阵元接收到的信号相对于参考阵元的时延为τm = (m·d·sinθ)/c,对应的相位差为 -2πf0τm = -2πm·d·sinθ/λ。

操作流程(构建流形向量):

  1. 确定阵列几何:列出所有阵元的空间坐标pm = [xm, ym, zm]ᵀ。
  2. 定义入射方向单位矢量:u(θ,φ) = [cosθ·cosφ, sinθ·cosφ, sinφ]ᵀ(球坐标定义)。
  3. 计算空间程差:Δm = pmᵀ·u(θ,φ)。
  4. 生成流形向量:a(θ,φ) = [e^{-j2πΔ0/λ}, e^{-j2πΔ1/λ}, ..., e^{-j2πΔM-1/λ}]ᵀ。

对于URA,流形向量可表示为x方向与y方向流形向量的Kronecker积:aURA(θ,φ) = ay(θ,φ) ⊗ ax(θ,φ),这一性质可大幅简化计算。

案例拆解:设某ULA阵列,M=4,d=λ/2,信号从θ=30°入射。则相邻阵元相位差为 -2π·(λ/2)·sin30°/λ = -π/2。流形向量为 a(30°) = [1, e^{-jπ/2}, e^{-jπ}, e^{-j3π/2}]ᵀ = [1, -j, -1, j]ᵀ。若信号复包络s(t)=1,则四路接收信号分别为1, -j, -1, j。这一组复数恰好构成了一个空间旋转的相位序列,后续波束形成正是通过加权求和将这些相位对齐。

2.3 窄带信号模型:从物理信号到复基带表示

核心原理:窄带假设是阵列信号处理中最常用的简化条件。所谓“窄带”,是指信号带宽B远小于载频f0,且信号包络在跨越阵列孔径所需的时间τmax内近似不变。数学上,接收信号可表示为:

x(t) = a(θ)·s(t) + n(t)
其中:
  x(t) ∈ C^{M×1}:M个阵元的接收快拍
  s(t) ∈ C^{1×1}:信号复包络(含幅度和相位)
  n(t) ∈ C^{M×1}:加性噪声,通常假设为空间白噪声
  a(θ) ∈ C^{M×1}:阵列流形向量

该模型将载波项e^{j2πf₀t}隐去,所有信息都包含在复包络s(t)中。对于多个信源(K个),模型扩展为x(t) = A(Θ)·s(t) + n(t),其中A(Θ) = [a(θ₁), a(θ₂), ..., a(θ_K)]为M×K维方向矩阵。

风险提示:窄带假设失效的典型场景包括:超宽带雷达(如UWB)、大孔径阵列接收宽带信号、或信号带宽与阵列孔径的乘积不满足B·(D/c) << 1。此时若仍使用窄带模型,会导致波束形成主瓣偏移和DOA估计偏差。工程上需改用宽带聚焦算法(如CSSM、WAVES)或时延补偿方法。

2.4 常见误区与设计复盘

误区一:认为d越小越好。实际上,d < λ/2会增大阵元间互耦,且降低阵列有效孔径,导致波束宽度变宽。d > λ/2则会产生栅瓣,使DOA估计出现模糊。标准做法是d = λ/2,若需扩展孔径,可采用稀疏阵列配合解模糊算法。

误区二:忽略阵元位置误差。实际阵列中,阵元安装位置偏差、通道幅相不一致都会破坏流形向量的理想结构。仿真中性能优异的算法,在实测中往往因流形失配而性能骤降。因此,高精度阵列标定(如使用旋转源或近场校正)是工程落地的必要环节。

误区三:将窄带模型直接套用于宽带信号。若信号带宽较宽,不同频率分量对应的流形向量不同,简单使用中心频率的a(θ)会引入模型误差。此时应使用频域子带处理或聚焦变换。

设计参数推荐值/做法影响验证方法
阵元间距 dλ/2(最高频点)避免栅瓣,平衡互耦仿真方向图检查栅瓣
阵元数 M按分辨率需求M越大,波束越窄,DOA精度越高计算CRB或仿真MUSIC谱
参考阵元选几何中心或第0号阵元影响流形向量相位基准一致性检查
通道校准幅度误差<0.5dB,相位误差<5°保证流形向量准确性近场校准或注入已知信号
复盘总结:本章的核心逻辑链为“几何构型 → 空间相位差 → 流形向量 → 窄带接收模型”。在后续章节中,无论是计算最优权值(如MVDR),还是进行子空间分解(如MUSIC),其数学推导的起点都是本章建立的模型。建议读者在进入下一章前,务必亲手推导一遍ULA的流形向量,并理解当d=λ/2时,流形向量恰好构成一个空间DFT矩阵的列,这一性质将直接关联到波束形成与空间谱估计的统一视角。