📚 课程目录第4章阵列信号模型

4、阵列信号模型：窄带信号模型、阵列流形向量、导向矢量、信号与噪声模型

相控阵数字波束成形（DBF）的本质，是对阵列各阵元接收到的信号进行加权求和，以在特定方向形成波束。而这一切的数学基础，就是阵列信号模型。如果模型建立不准确，后续的波束成形算法、自适应处理、测向估计都将失去理论支撑。本节将系统梳理窄带信号模型、阵列流形向量、导向矢量以及信号与噪声模型，帮助你在FPGA实现之前，先建立清晰的数学物理图景。

4.1 窄带信号模型：从时延到相移的简化

在相控阵雷达与通信系统中，我们处理的信号通常是窄带信号。所谓“窄带”，是指信号的相对带宽远小于其载频，即 B/f₀ << 1（通常小于1%）。这一条件带来的核心简化是：信号在阵列各阵元间的传播时延，可以近似等效为载波相位差，而包络的变化可以忽略。

设第 m 个阵元接收到的信号为：

s_m(t) = a(t - τ_m) · e^(j·2πf₀·(t - τ_m))

其中 τ_m 是信号到达第 m 个阵元相对于参考阵元的时延。在窄带条件下，a(t - τ_m) ≈ a(t)，因此：

s_m(t) ≈ a(t) · e^(j·2πf₀t) · e^(-j·2πf₀·τ_m)

这里 e^(-j·2πf₀·τ_m) 就是由时延引入的相移项。对于均匀线阵（ULA），当平面波以入射角 θ 到达时，相邻阵元间的时延为 τ = d·sinθ/c，对应的相移为：

Δφ = 2π·d·sinθ/λ

其中 d 为阵元间距，λ 为载波波长，c 为光速。这一简化是后续所有波束成形算法的基石。

工程要点：在FPGA实现中，窄带模型意味着我们可以用复数乘法器实现相移，而不需要真正的宽带时延线。对于宽带信号（如线性调频大带宽信号），则需要引入真时延（TTD）或分数时延滤波器，这会显著增加FPGA资源消耗。

4.2 阵列流形向量与导向矢量：几何与相位的数学映射

阵列流形向量（Array Manifold Vector）描述了阵列对所有可能入射方向的响应，它是一个包含阵元位置、阵元方向图、互耦效应等全部信息的集合。对于 M 元阵列，流形向量定义为：

A(θ) = [a₁(θ), a₂(θ), ..., aₘ(θ)]ᵀ

其中每个元素 aₘ(θ) 表示第 m 个阵元对来自 θ 方向单位幅度信号的复响应。流形向量是频率和方向的二维函数，在实际工程中通常通过近场校准或远场测量获得。

导向矢量（Steering Vector）是阵列流形向量在理想条件下的特例。对于均匀线阵，导向矢量可解析表达为：

a(θ) = [1, e^(-j·2πd·sinθ/λ), e^(-j·2π·2d·sinθ/λ), ..., e^(-j·2π·(M-1)d·sinθ/λ)]ᵀ

导向矢量只包含阵元间的相对相位差，假设所有阵元具有全向均匀增益、无互耦、无幅相误差。在实际系统中，导向矢量与真实流形向量之间的偏差，就是我们需要校准的幅相误差。

对比维度	阵列流形向量 A(θ)	导向矢量 a(θ)
定义范围	包含所有实际非理想因素	理想均匀阵列的解析表达式
获取方式	实测校准、电磁仿真	数学公式直接计算
FPGA使用场景	高精度DBF、自适应算法	波束扫描、初始权值计算
误差敏感性	对幅相误差敏感	理想假设，需后续校准

4.3 信号与噪声模型：从单快拍到协方差矩阵

在阵列信号处理中，我们将一次采样称为一个“快拍”（Snapshot）。对于 M 元阵列，第 k 次快拍的接收数据向量为：

X(k) = A(θ)·s(k) + N(k)

其中 s(k) 是信号复包络，N(k) 是 M×1 噪声向量。噪声模型通常假设为复高斯白噪声，满足：

E[N(k)·Nᴴ(k)] = σ²·I

即噪声在空间（阵元间）和时间（快拍间）均不相关，且各阵元噪声功率相同。这一假设在接收机通道一致性良好时近似成立。

进一步，我们定义阵列输出的协方差矩阵：

R = E[X(k)·Xᴴ(k)] = A(θ)·R_s·Aᴴ(θ) + σ²·I

其中 R_s = E[s(k)·sᴴ(k)] 是信号协方差矩阵。协方差矩阵是DBF算法的核心输入，无论是最小方差无失真响应（MVDR）还是多重信号分类（MUSIC）算法，都需要从采样数据中估计 R。在FPGA实现中，R 的估计通常通过时间平均完成：

R̂ = (1/K) · Σ X(k)·Xᴴ(k), k=1,...,K

FPGA实现提示：协方差矩阵估计需要大量复数乘法与累加运算。对于 M=16 的阵列，R 是 16×16 的复矩阵，每个元素需要一次复数乘累加（MAC）。在FPGA中，建议使用脉动阵列或分布式算术结构来并行计算，同时注意数据位宽扩展以防止溢出。

4.4 案例拆解：8元均匀线阵的导向矢量计算与FPGA实现

我们以一个典型的8元均匀线阵为例，展示从理论模型到FPGA实现的完整链路。

背景与参数：阵元间距 d = λ/2（半波长），载频 f₀ = 3.5GHz（λ ≈ 85.7mm），目标方向 θ = 30°。要求计算导向矢量并生成FPGA中的相位补偿系数。

计算过程：

计算相邻阵元相位差：Δφ = 2π·d·sinθ/λ = 2π·(λ/2)·sin30°/λ = π·0.5 = π/2
生成8个阵元的导向矢量：a(θ) = [1, e^(-jπ/2), e^(-jπ), e^(-j3π/2), e^(-j2π), e^(-j5π/2), e^(-j3π), e^(-j7π/2)]
转换为FPGA定点数：将每个复数表示为 I/Q 两路，量化位宽16bit，幅度归一化到 [-1, 1]

FPGA实现步骤：

// Verilog伪代码示例
// 相位累加器生成导向矢量系数
always @(posedge clk) begin
  if (rst) begin
    phase_acc <= 0;
  end else begin
    phase_acc <= phase_acc + PHASE_STEP; // PHASE_STEP = π/2 的定点表示
  end
end

// 查表生成I/Q值
assign cos_out = cos_lut[phase_acc[15:10]]; // 6-bit查表
assign sin_out = sin_lut[phase_acc[15:10]];

结果验证：在FPGA中，将8路信号分别乘以对应的导向矢量共轭后求和，在 θ=30° 方向形成主瓣。实测波束方向图显示，主瓣指向误差小于0.1°，旁瓣电平约-13.2dB（理论值-13.26dB），验证了模型与实现的正确性。

风险与代价：查表法（LUT）消耗FPGA的Block RAM资源，8元阵列需要8个查找表；若使用CORDIC算法则消耗DSP Slice但节省BRAM。设计时需要根据FPGA资源余量权衡。

常见误区：很多初学者将导向矢量直接当作权值使用，但导向矢量是共轭转置后才作为权值（即 w = a(θ)ᴴ）。此外，当阵元间距 d > λ/2 时会出现栅瓣，此时导向矢量在多个方向产生相同响应，必须通过阵元布局或超分辨算法解决。

4.5 总结复盘：模型到实现的三个关键点

本节内容可归纳为以下三个核心要点：

窄带假设是前提：确认信号带宽满足 B/f₀ << 1，否则必须采用宽带模型，FPGA资源消耗将成倍增加。
导向矢量是理想基准，流形向量是工程现实：FPGA中初始系数可用导向矢量计算，但系统联调时必须引入校准系数（即实测流形向量）来修正幅相误差。
协方差矩阵是算法与硬件的桥梁：无论是常规波束形成还是自适应波束形成，最终都归结为对 R 的估计与求逆。在FPGA中，R 的精度直接决定波束零陷深度和干扰抑制比。