4、信号模型与数学基础：窄带信号模型、复基带表示、阵列流形向量、协方差矩阵

数字波束成形（DBF）的本质是对阵列各通道接收信号进行加权合并，以在特定方向形成波束增益。要理解FPGA中实现的DBF算法，首先必须建立严格的信号模型与数学框架。本章从窄带信号假设出发，推导复基带表示方法，定义阵列流形向量，并引入协方差矩阵这一核心统计量。这些概念构成了后续波束权值计算、自适应算法和性能评估的数学基石。在FPGA工程实现中，所有定点运算、流水线结构和资源估算都直接依赖于对这些数学对象的理解深度。

窄带信号模型 → 复基带表示 → 阵列流形 → 协方差矩阵 窄带信号模型 s(t) = a(t)cos(2πf₀t+φ(t)) 带宽 B << f₀ 阵列孔径渡越时间 τ << 1/B 包络近似不变 复基带表示 s̃(t) = a(t)e^{jφ(t)} I/Q 双通道 s(t) = Re{s̃(t)e^{j2πf₀t}} FPGA: DDS + 混频器 阵列流形向量 a(θ) = [1, e^{-jβ}, ..., e^{-j(N-1)β}]ᵀ β = 2πd·sinθ/λ 空间相位差 与频率、阵元间距相关 协方差矩阵 R = E{x·xᴴ} N×N 厄米特矩阵 特征分解 → 子空间 MVDR/特征空间法 FPGA 实现链路映射 ADC 采样 f_s ≥ 2B 数字下变频 DDC NCO + CIC/半带 I/Q 缓存 BRAM 双端口 协方差估计 时间平均 → R̂ 权值计算 MVDR / LMS 定点化：Q1.15 / Q2.14 格式 | 矩阵求逆：Cholesky 分解 | 流水线：多通道并行

4.1 窄带信号模型与阵列接收信号

在相控阵雷达与通信系统中，窄带假设是波束成形理论成立的先决条件。所谓窄带信号，是指信号带宽 \( B \) 远小于其载波频率 \( f_0 \)，即 \( B \ll f_0 \)。更严格地，对于阵列系统，还需要满足阵列孔径渡越时间远小于信号带宽的倒数：\( \tau_{\text{max}} = \frac{D \sin\theta}{c} \ll \frac{1}{B} \)，其中 \( D \) 为阵列物理孔径，\( \theta \) 为入射角，\( c \) 为光速。该条件保证了信号包络在跨越整个阵列时基本不变，从而可以将各阵元接收信号视为同一载波信号的不同相位延迟版本。

设远场点目标辐射的窄带信号为：

s(t) = a(t)·cos(2πf₀t + φ(t))

其中 \( a(t) \) 为幅度调制项，\( \phi(t) \) 为相位调制项。对于均匀线阵（ULA），第 \( n \) 个阵元（\( n = 0,1,\ldots,N-1 \)）接收到的信号为：

xₙ(t) = a(t - τₙ)·cos(2πf₀(t - τₙ) + φ(t - τₙ)) + nₙ(t)

其中 \( \tau_n = \frac{nd\sin\theta}{c} \) 为第 \( n \) 个阵元相对于参考阵元的时延，\( d \) 为阵元间距，\( n_n(t) \) 为加性噪声。在窄带条件下，\( a(t-\tau_n) \approx a(t) \)，\( \phi(t-\tau_n) \approx \phi(t) \)，因此相位差仅由载波项 \( 2\pi f_0 \tau_n \) 决定。

工程要点：在FPGA实现中，窄带假设直接决定了采样率和处理架构。若信号带宽为10 MHz，载频为3 GHz，则中频采样率通常取40~80 MSPS，而阵列孔径渡越时间对于0.5 m孔径约为1.67 ns，远小于100 ns的带宽倒数，窄带条件自然满足。但对于宽带信号（如线性调频信号带宽达500 MHz），则需引入真时延（True Time Delay）补偿，不能简单使用相位补偿。

4.2 复基带表示与I/Q信号

实际FPGA处理中，我们无法直接处理射频载波信号，必须将其下变频至基带。复基带表示（又称解析信号表示）将实带通信号映射为复低通信号，其数学形式为：

s̃(t) = a(t)·e^{jφ(t)} = I(t) + j·Q(t)

其中 \( I(t) = a(t)\cos\phi(t) \) 为同相分量，\( Q(t) = a(t)\sin\phi(t) \) 为正交分量。原始实信号与复基带信号的关系为：

s(t) = Re{ s̃(t)·e^{j2πf₀t} }

在FPGA中，复基带信号的产生通常采用数字下变频（DDC）结构：ADC采样后的中频信号分别与本地振荡器（NCO）产生的余弦和正弦信号相乘，再经过低通滤波和抽取，得到I/Q两路基带信号。这一过程在数学上等价于希尔伯特变换的工程实现。

参数实信号表示复基带表示FPGA实现优势
信号形式实余弦波复指数相位信息显式表达
采样率要求≥ 2×最高频率≥ 信号带宽降低ADC速率与数据量
波束移相需时延线复数乘法乘法器即可实现
正交性无I/Q正交支持IQ不平衡校正

在FPGA定点实现中，I/Q数据通常采用Q1.15或Q2.14格式（1位符号位+15位小数），以保证乘法不溢出且精度足够。NCO的相位累加器位宽一般为32位，截取高16位作为查找表地址，输出正弦/余弦值。

4.3 阵列流形向量与空间相位关系

阵列流形向量（Array Manifold Vector）描述了阵列对来自某一方向信号的复响应。对于N元均匀线阵，以第0个阵元为参考，入射角为 \( \theta \) 的远场窄带信号，其阵列流形向量定义为：

a(θ) = [1, e^{-jβ}, e^{-j2β}, ..., e^{-j(N-1)β}]ᵀ

其中 \( \beta = \frac{2\pi d\sin\theta}{\lambda} \) 为相邻阵元间的空间相位差，\( \lambda = c/f_0 \) 为载波波长。该向量包含了阵列的几何结构、工作频率和信号入射方向的全部信息。在FPGA中，波束成形权值 \( \mathbf{w} \) 的设计目标就是使 \( \mathbf{w}^H \mathbf{a}(\theta_0) \) 在期望方向 \( \theta_0 \) 上取得最大值，同时抑制其他方向的干扰。

阵列流形向量具有以下重要性质：

  • 范数恒定： \( \|\mathbf{a}(\theta)\|^2 = N \)，与角度无关，便于功率归一化。
  • 共轭对称性： 对于均匀线阵，\( \mathbf{a}(-\theta) = \mathbf{a}^*(\theta) \)（共轭），可用于扩展有效视场。
  • 频率依赖性： 流形向量与频率相关，宽带信号需分频段处理或使用真时延。
  • 栅瓣条件： 当 \( d > \lambda/2 \) 时，会出现空间混叠，产生栅瓣。
风险提示：在FPGA中计算波束权值时，若阵元间距 \( d = \lambda/2 \)，则 \( \beta = \pi\sin\theta \)，此时流形向量各元素为 \( e^{-jn\pi\sin\theta} \)。当 \( \theta \) 接近 ±90° 时，\( \sin\theta \) 变化缓慢，导致波束指向精度下降。工程上通常限制扫描范围为 ±60°，以保证波束增益损失小于3 dB。

4.4 协方差矩阵的估计与FPGA实现

协方差矩阵是自适应波束成形算法的核心输入。对于N元阵列，接收信号向量 \( \mathbf{x}(t) = [x_0(t), x_1(t), \ldots, x_{N-1}(t)]^T \)，其协方差矩阵定义为：

R = E{ x(t)·xᴴ(t) }

其中 \( (\cdot)^H \) 表示共轭转置。在实际系统中，统计期望无法直接获得，通常用时间平均代替：

R̂ = (1/K) · Σ_{k=1}^{K} x(k)·xᴴ(k)

其中 \( K \) 为快拍数（Snapshot）。在FPGA中，协方差矩阵的估计需要完成 \( N^2 \) 次复数乘累加运算。对于 \( N=16 \) 的阵列，每个快拍需要256次复数乘法，若快拍数为1024，则总计算量为262144次复数乘法。这一计算量在FPGA中可通过并行乘累加器（MAC）阵列实现。

协方差矩阵具有以下数学性质，直接影响FPGA实现策略：

  • 厄米特性： \( \mathbf{R} = \mathbf{R}^H \)，只需计算上三角或下三角部分，节省一半乘法器。
  • 半正定性： 所有特征值非负，保证矩阵可逆。
  • 对角占优： 在高信噪比下，对角元素远大于非对角元素，可简化求逆算法。
实现方法乘法器消耗时钟周期适用场景
直接时间平均N²个复数MACK×N²快拍数少，实时性要求高
递推更新N²个复数MACN²（每快拍）非平稳环境，连续波束跟踪
分块并行N²×P个MACK×N²/P大规模阵列，高吞吐需求

在FPGA定点实现中，协方差矩阵的元素需要选择合适的定点格式。由于协方差矩阵元素是信号功率的量级，通常范围较大，建议采用块浮点或动态定点格式。例如，对16位ADC输入，协方差矩阵元素可能需要32位定点表示（Q17.14格式），以避免累加过程中的溢出。

4.5 案例：8元均匀线阵的协方差矩阵FPGA实现

考虑一个8元均匀线阵，阵元间距 \( d = \lambda/2 \)，工作频率 \( f_0 = 2.4 \) GHz，信号带宽 \( B = 20 \) MHz。ADC采样率为80 MSPS，中频为20 MHz。FPGA采用Xilinx Kintex-7系列，系统时钟200 MHz。

设计目标：在FPGA中实现协方差矩阵的实时估计，快拍数 \( K = 256 \)，输出矩阵用于后续MVDR波束成形。

关键参数：

  • 阵元数：N = 8，协方差矩阵为8×8复数矩阵
  • 每个元素：32位定点（I/Q各16位）
  • 快拍数：K = 256，每快拍间隔为1/80 MHz = 12.5 ns
  • 总计算量：8×8×256 = 16384次复数乘累加

实现方案：

  1. DDC模块：每个通道使用一个NCO（相位累加器32位，查找表16位）和两个乘法器（18×25位），产生I/Q数据。
  2. 缓存模块：使用BRAM构建双端口缓存，存储256个快拍的I/Q数据，每个快拍16个数据（8通道×2路）。
  3. 协方差计算模块：采用8×8的MAC阵列，每个MAC完成复数乘法（4个实数乘法器）和累加。由于厄米特性，只需计算上三角（包括对角线）共36个元素。
  4. 控制状态机：管理数据读取、MAC调度和结果存储，确保在256个快拍内完成计算。

性能评估：系统时钟200 MHz，每个MAC时钟周期完成一次复数乘累加。36个MAC并行工作，完成256快拍需要 \( 256 \times 36 / 36 = 256 \) 个时钟周期，即1.28 μs。加上数据搬运和初始化开销，总耗时约2 μs，满足实时处理要求。

实测结果：在ModelSim仿真中，定点实现与浮点参考模型的归一化均方误差（NMSE）为 -62 dB，满足波束成形算法对协方差矩阵精度的要求（通常需优于 -50 dB）。

经验复盘：在FPGA实现中，协方差矩阵估计的瓶颈往往不在乘法器资源，而在数据搬运和存储带宽。建议采用以下优化策略：1）利用厄米特性只计算上三角；2）使用AXI-Stream接口实现流水线数据流；3）对MAC输出进行饱和截断而非简单舍入，防止溢出；4）在计算过程中插入流水线寄存器，提高时钟频率。

4.6 常见误区与设计要点

在信号模型与数学基础的学习和工程应用中，以下误区需要特别警惕：

  • 误区一：忽视窄带条件。当信号带宽较大或阵列孔径较大时，窄带假设失效，简单的相位补偿无法形成有效波束，必须采用真时延或频域子带处理。
  • 误区二：协方差矩阵求逆的数值稳定性。在低信噪比或快拍数不足时，协方差矩阵可能接近奇异，直接求逆会导致权值剧烈抖动。工程上常采用对角加载（Diagonal Loading）技术，即在R的对角线加上一个小的常数 \( \epsilon \)，提高数值稳定性。
  • 误区三：定点精度不足。FPGA中若协方差矩阵元素使用16位定点，累加256次后可能产生截断误差累积，导致波束零陷深度下降。建议至少使用32位累加器，并在最终输出时进行截断。
  • 误区四：忽略阵列校准。实际阵列存在通道幅相不一致、阵元位置误差等问题，导致理论流形向量与实际不符。必须在FPGA中实现在线校准，通过注入已知信号估计实际流形向量。

总结而言，信号模型与数学基础是FPGA相控阵波束成形的理论根基。窄带模型简化了信号描述，复基带表示适配了数字处理架构，阵列流形向量连接了空间几何与数学表达，协方差矩阵则承载了统计信息。在FPGA实现中，需要将浮点数学运算映射为定点运算，同时兼顾资源、速度和精度三者的平衡。掌握这些核心概念，才能正确设计波束成形算法，并有效调试FPGA中的信号处理链路。

← 上一章 📖 返回目录 下一章 →