【自适应天线与相控阵技术】适应天线系统中的 DPCA 雷达杂波抑制性能分析

在自适应天线系统中，阵列互耦效应对位移相位中心天线（Displaced Phase Center Antenna, DPCA）雷达杂波抑制性能具有显著影响。本文系统阐述 DPCA 杂波对消的基本原理，建立考虑阵列互耦的双通道信号模型，推导杂波协方差矩阵与对消因子的解析表达式，并基于矩量法仿真与实测数据，定量分析阵列规模、辐射单元类型、栅格形式、边缘保护带、相位中心位移及扫描角度等关键设计参数对杂波对消性能的影响规律。

1. 引言

自适应天线与相控阵技术自二十世纪五十年代以来已被广泛研究。在雷达应用中，自适应天线系统的核心功能包括：抑制干扰与噪声、抑制雷达杂波，以及检测运动目标。对于机载或星载雷达平台，地杂波（Ground Clutter）是限制动目标检测性能的主要因素。

当雷达平台相对于地面运动时，静止杂波的多普勒频谱会被平台运动展宽，从而淹没慢速运动目标。位移相位中心天线（DPCA）技术通过在同一雷达脉冲重复间隔（PRI）内利用空间分离的两个相位中心进行收发，形成等效的静止单基地雷达观测点，从而补偿平台运动引起的杂波相位变化，使静止杂波在两次观测中保持高度相关而被对消，同时运动目标因额外的径向速度产生相位差得以保留。

然而，实际相控阵天线存在单元间的互耦效应。当阵列被划分为两个接收子阵形成 DPCA 相位中心时，前导子阵与后随子阵所处的电磁环境具有镜像不对称性，导致两个相位中心的方向图、相位中心位置及杂波统计特性存在差异。这种由互耦引起的失配会直接限制 DPCA 的杂波对消能力。本文将建立完整的理论模型，系统分析互耦效应对 DPCA 性能的影响。

2. DPCA 杂波对消原理

2.1 运动平台与杂波几何

考虑一个运动的雷达平台（机载或星载），以速度 $v$ 沿航迹方向飞行。设雷达天线为相控阵，在脉冲重复周期 $T_r$ （PRI）内，天线相位中心从位置 $A$ 移动到位置 $B$，位移量为 $d = vT_r$。

在第一个发射时刻，天线以全孔径从相位中心 $A$ 发射脉冲；经过一个 PRI 后，天线移动到相位中心 $B$ 并发射第二个脉冲。对于接收，需要在两个时刻分别形成接收相位中心 $A'$ 和 $B'$，使得等效的单基地雷达位置重合。

具体而言，当接收第一个脉冲回波时，天线已运动到 $B$ 附近，此时通过子阵加权将接收相位中心后移（向 $A$ 的方向调整），使等效接收相位中心与发射相位中心 $A$ 重合；同理处理第二个脉冲。这样，两次观测对应于空间中的同一等效相位中心，静止杂波在两次回波中的相位完全一致，而运动目标因自身径向速度 $v_t$ 在间隔 $T_r$ 内产生额外的多普勒相移 $\Delta \phi = 4\pi v_t T_r / \lambda$，从而可被检测。

2.2 杂波反射特性

对于 L 波段星载或机载雷达，地杂波的平均雷达截面积（RCS）随擦地角（Grazing Angle）变化。当擦地角接近 $90^\circ$（即正下视）时，杂波反射最强；随着擦地角减小（向地平线方向），杂波强度逐渐降低。因此，天线单元在法向（天底方向）的辐射特性对杂波抑制具有重要影响。

3. 阵列架构与互耦效应

3.1 双相位中心接收架构

考虑一个由 $N$ 个单元组成的相控阵。通过 corporate 馈电网络或 T/R 组件的双路输出，阵列被划分为两个独立的接收相位中心：前导相位中心 A（Leading Phase Center）和后随相位中心 B（Trailing Phase Center）。

每个单元信号经过功分器后分别馈入相位中心 $A$ 和 $B$ 的接收通道。通过调节 T/R 模块的幅度与相位权重，可以控制两个相位中心的位置。设两个相位中心沿航迹方向的位移为 $\Delta$。

存在两种极端的子阵划分方式：

重叠孔径（Overlapped Apertures）： 两个相位中心共享大部分阵元，仅存在较小的位移 $\Delta$；
分离孔径（Split Apertures）： 两个相位中心由完全不重叠的子阵构成，位移 $\Delta$ 达到最大，约为阵列孔径的一半。

实际系统中，通过 T/R 模块的幅度加权可在硬件上实现相位中心的连续调节。

3.2 互耦不对称性

由于阵列边缘效应，前导子阵 $A$ 与后随子阵 $B$ 所处的互耦环境存在镜像不对称性。具体而言，若将整个孔径视为具有几何对称性，子阵 $A$ 在其右侧（航迹前方）看到的邻近单元数量与子阵 $B$ 在其左侧（航迹后方）看到的单元数量不同。这种不对称的互耦环境会轻微改变两个子阵的波束形状，使得它们的辐射方向图、相位中心位置及阻抗特性产生差异，进而降低杂波对消性能。

为缓解这种不对称性，通常在有效阵列边缘设置若干行被动端接单元（Passively Terminated Elements），称为保护带（Guard Bands）。保护带通过模拟无限大阵列的局部环境，使两个相位中心的边缘条件趋于一致，从而改善互耦对称性。

3.3 辐射单元类型与仿真方法

本文主要分析两种典型的线天线单元：

单极子（Monopole）： 长度约为 $\lambda/4$，垂直于地面安装。其理论单元方向图在法向（$\theta = 0^\circ$，天底方向）存在零点，有利于抑制正下视强杂波；但宽角扫描能力受限，通常用于 $30^\circ$ – $60^\circ$ 的扫描扇区。
偶极子（Dipole）： 长度约为 $\lambda/2$，通常平行于地面或按特定方式架设。其覆盖范围较宽，可扫描至法向附近，但法向无零点，会引入更多天底杂波。

在仿真中，阵列置于无限大导电地面上方。利用镜像法（Method of Images），将地面上方的阵列等效为自由空间中的原阵列与其镜像阵列的组合，从而可方便地应用矩量法（Method of Moments, MoM）计算单元间的互耦系数与有源单元方向图。

4. 信号模型与杂波对消理论推导

4.1 双通道接收信号模型

设阵列在存在互耦时的接收信号向量为 $\mathbf{x} \in \mathbb{C}^{N \times 1}$。由于互耦效应，实际接收电压与理想各向同性阵列不同。互耦矩阵 $\mathbf{C} \in \mathbb{C}^{N \times N}$ 描述了单元间的耦合关系，使得理想入射信号 $\mathbf{x}_0$ 与实际信号满足：

$$\mathbf{x} = \mathbf{C} \mathbf{x}_0$$

对于 DPCA 的双通道接收，通过子阵选择矩阵 $\mathbf{S}_A, \mathbf{S}_B \in \mathbb{C}^{N_A \times N}, \mathbb{C}^{N_B \times N}$（通常为对角加权矩阵），提取两个子阵的输出：

$$x_A = \mathbf{w}_A^H \mathbf{S}_A \mathbf{x}, \quad x_B = \mathbf{w}_B^H \mathbf{S}_B \mathbf{x}$$

其中 $\mathbf{w}_A, \mathbf{w}_B$ 为子阵 $A$ 和 $B$ 的波束形成权重。

4.2 杂波协方差矩阵

假设杂波为广义平稳随机过程，两个通道接收的杂波信号构成随机向量 $\mathbf{y} = [y_A, y_B]^T$。其协方差矩阵（即杂波对消矩阵或方向图相关矩阵）为：

$$\mathbf{R}_c = E[\mathbf{y} \mathbf{y}^H] = \begin{bmatrix} R_{11} & R_{12} \\ R_{21} & R_{22} \end{bmatrix}$$

其中各元素的物理含义为：

$$R_{ij} = \int_{\Omega} P_t^2(\theta) \cdot P_{r,i}(\theta) P_{r,j}^*(\theta) \cdot \sigma_c(\theta) \cdot e^{j \phi_{ij}(\theta)} d\Omega$$

式中：

$P_t(\theta)$：全孔径发射方向图；
$P_{r,i}(\theta), P_{r,j}(\theta)$：第 $i,j$ 个接收相位中心的方向图（包含互耦效应）；
$\sigma_c(\theta)$：杂波散射系数（与擦地角相关）；
$\phi_{ij}(\theta)$：由相位中心位移引入的相位因子，理想情况下用于补偿平台运动。

对于窄带远场条件，若平台速度为 $v$，PRI 为 $T_r$，理想相位中心位移应为 $d = v T_r / 2$（单基地等效），则相位因子近似为：

$$\phi_{AB} = \frac{2\pi}{\lambda} \cdot 2d \sin\theta = \frac{4\pi d}{\lambda} \sin\theta$$

当实际相位中心位置偏离理想值时，该相位因子无法完全补偿杂波的多普勒相位，导致 $R_{12}$ 的模值下降。

4.3 杂波对消因子的严格推导

DPCA 杂波对消的本质是利用两个通道的杂波相关性进行自适应或固定加权对消。考虑最优两通道杂波抑制：寻找复加权向量 $\mathbf{w} = [w_1, w_2]^T$，使得输出杂波功率最小：

$$\min_{\mathbf{w}} \mathbf{w}^H \mathbf{R}_c \mathbf{w}, \quad \text{s.t.} \quad \mathbf{w}^H \mathbf{u} = 1$$

其中 $\mathbf{u}$ 为约束向量（例如 $\mathbf{u} = [1, 0]^T$ 保证目标信号无失真）。利用拉格朗日乘子法，最优权重为：

$$\mathbf{w}_{opt} = \frac{\mathbf{R}_c^{-1} \mathbf{u}}{\mathbf{u}^H \mathbf{R}_c^{-1} \mathbf{u}}$$

对应的输出杂波功率（残余杂波）为：

$$P_{res} = \frac{1}{\mathbf{u}^H \mathbf{R}_c^{-1} \mathbf{u}}$$

对于简单的两通道对消（如 DPCA 相减），更常用的性能度量是杂波对消因子（Clutter Cancellation Factor） $C$，定义为最优对消后残余杂波功率与输入杂波功率之比。通过直接计算 $2 \times 2$ 矩阵的逆，可得：

$$C = \frac{\det(\mathbf{R}_c)}{R_{11} R_{22}} = 1 - \frac{|R_{12}|^2}{R_{11} R_{22}}$$

推导过程：

设 $\mathbf{R}_c = \begin{bmatrix} R_{11} & R_{12}^* \\ R_{12} & R_{22} \end{bmatrix}$，其逆矩阵为：

$$\mathbf{R}_c^{-1} = \frac{1}{\det(\mathbf{R}_c)} \begin{bmatrix} R_{22} & -R_{12}^* \\ -R_{12} & R_{11} \end{bmatrix}$$

取 $\mathbf{u} = [1, 1]^T$（对应两通道相减对消，即约束 $w_1 + w_2 = 1$ 的简化情形，或直接考虑相干系数），则：

$$\mathbf{u}^H \mathbf{R}_c^{-1} \mathbf{u} = \frac{R_{11} + R_{22} - 2\text{Re}(R_{12})}{\det(\mathbf{R}_c)}$$

若进一步假设两通道增益平衡（$R_{11} = R_{22} = P_c$），且 $R_{12} = \rho P_c$（$\rho$ 为复相关系数），则：

$$C = 1 - |\rho|^2 = 1 - \frac{|R_{12}|^2}{R_{11} R_{22}}$$

该式即为 DPCA 杂波对消因子的标准表达式。当两通道杂波完全相关时（$|\rho| = 1$），$C = 0$，残余杂波为零，对消比无限大；当存在方向图失配、互耦不对称或相位中心位置误差时，$|\rho| < 1$，$C > 0$，残余杂波功率增大。

杂波对消比（Cancellation Ratio） 定义为：

$$CR = 10 \log_{10} \left( \frac{1}{C} \right) = -10 \log_{10} \left( 1 - \frac{|R_{12}|^2}{R_{11} R_{22}} \right)$$

在后续分析中，该对消比将作为评估 DPCA 性能的核心指标。

5. 仿真与实测验证

5.1 实验阵列参数

为验证理论模型，构建了一个 L 波段实验阵列，主要参数如下：

单元数： 96 个有源单元；
孔径尺寸： 约 1.5m；
栅格： 六边形（Hexagonal Lattice）；
地面： 直径 2.1m 的导电地面（仿真中近似为无限大地面）；
边缘终端： 2 圈被动端接保护带单元。

在仿真中，分别采用单极子与偶极子单元，保持相同的六边形栅格与阵列布局，通过矩量法计算互耦系数及有源方向图。

5.2 方向图与对消性能对比

5.2.1 单极子阵列仿真与实测

对于 96 单元单极子相控阵，发射采用全孔径，并施加 10dB 余弦幅度加权。发射波束偏离法向（Broadside）扫描。接收时形成两个相位中心 $A$ 和 $B$。

当两相位中心间隔为 3 列单元时，两个接收方向图的幅度与相位在扫描角（如 $40^\circ$）附近表现出良好的一致性，但仍存在细微差异。这种差异导致杂波对消比限制在约 42dB（仿真值，无 T/R 模块误差）。

实测数据中，由于 T/R 模块存在幅度/相位误差，对消比进一步受限至约 40dB。仿真结果（实线）与实测结果（点）随扫描角的变化趋势一致：随着扫描角从 $30^\circ$ 增大到 $55^\circ$，对消性能逐渐下降，表明互耦效应在宽角扫描时更为显著。

5.2.2 偶极子阵列仿真

对于相同布局的偶极子阵列，扫描可从法向附近（$10^\circ$）开始至 $55^\circ$。仿真显示，偶极子阵列的杂波对消曲线相对平坦，对消比在 40dB 至 60dB 之间变化。在小尺寸阵列（96 单元）条件下，偶极子与单极子的互耦行为存在明显差异：单极子阵列的对消性能随相位中心位移增加而下降更快。

5.3 相位中心位移的影响

固定阵列规模（96 单元），改变两相位中心沿航迹方向的列间距（位移 $\Delta$）：

单极子阵列： 当位移从 3 列增加到 6 列时，每个子阵的有效单元数减少，边缘效应增强，互耦不对称性加剧，导致对消性能显著下降。在 $40^\circ$ – $55^\circ$ 扫描范围内，对消比从约 50dB 降至 30dB 以下。
偶极子阵列： 同样呈现随位移增加而性能下降的趋势，但下降幅度较单极子缓和。

物理机制： 位移增大意味着每个接收子阵的孔径减小，子阵边缘单元占比增加，互耦对方向图的影响更为突出，使得两个相位中心的方向图失配加剧。

6. 参数敏感性分析

6.1 保护带（Guard Bands）的作用

对于 96 单元阵列（3 列位移，$40^\circ$ 与 $55^\circ$ 扫描角），对比无保护带与 2 圈保护带的仿真结果：

扫描角	无保护带对消比	2圈保护带对消比
$40^\circ$	$\approx 37\text{dB}$	$> 50\text{dB}$
$55^\circ$	$\approx 35\text{dB}$	$\approx 43\text{dB}$

保护带的引入使对消性能提升约 10 – 15dB。其物理机制在于：

改善互耦对称性： 保护带使两个相位中心的边缘环境趋于一致，减少了因镜像不对称引起的方向图差异；
校准相位中心位置： 无保护带时，互耦效应使实际相位中心偏离理论设计位置；保护带的存在使实际相位中心更接近预期值，从而改善了 DPCA 的相位补偿精度。

6.2 阵列行数的影响

对于大型 DPCA 相控阵（固定 128 列，孔径约 16m，采用分离孔径即最大相位中心位移，2 圈保护带，六边形栅格），研究有源行数从 4 行变化到 16 行时的对消性能。

结果表明：在 $40^\circ$ 与 $55^\circ$ 扫描角下，杂波对消比基本保持平坦，不随行数显著变化。这说明对于大型阵列，行数（即阵列在俯仰向的尺寸）并非影响 DPCA 对消性能的主要因素。核心影响因素在于航迹向的相位中心位移与扫描角度。

6.3 阵列列数（长度）的影响

固定行数为 8 行，改变有源列数（即阵列在航迹向的长度），分别研究方形栅格与六边形栅格下的单极子与偶极子阵列：

方形栅格： 随着列数增加（阵列增长），单元总数增加，子阵孔径相对边缘效应减弱，对消性能整体提升。在 $55^\circ$ 扫描角下，小阵列（少列数）中偶极子的对消性能略优于单极子；但在大型阵列中，单极子的对消性能反超偶极子。这种交叉现象源于不同尺寸下互耦环境的改变。
六边形栅格： 单极子与偶极子的性能差异较方形栅格有所减小，但列数增加带来的性能提升趋势保持一致。

总体而言，阵列长度（列数）的增加有助于提升 DPCA 杂波对消性能，因为更大的孔径使子阵方向图更稳定，互耦引起的相对失配减小。

6.4 大阵列相位中心位移的不敏感性

对于 $8 \times 128$（1024 个有源单元）的大型阵列（2 圈保护带，六边形栅格），研究相位中心位移从较小值变化到半孔径分离时的对消性能：

仿真结果表明，当阵列规模足够大时，杂波对消比对相位中心位移的变化相对不敏感。在固定扫描角下，不同位移对应的对消曲线趋于平坦，对消比可维持在 65dB 以上。这是因为大阵列的子阵即使经过位移分割，仍包含足够多的单元，边缘互耦效应对整体方向图的影响被平均化。

需要注意的是，上述高对消比是在无 T/R 模块误差的理想仿真条件下获得的；实际硬件中，T/R 通道的幅度/相位误差通常会将性能限制在 40dB 左右。

7. 结论

本文系统分析了阵列互耦效应对 DPCA 雷达杂波抑制性能的影响，主要结论如下：

DPCA 原理： 通过双相位中心形成等效静止单基地雷达观测点，利用静止杂波的高相关性与运动目标的多普勒相移实现杂波对消与目标检测。对消性能由杂波协方差矩阵的非对角相关性决定，对消因子 $C = 1 - |R_{12}|^2 / (R_{11} R_{22})$ 定量描述了残余杂波功率。
互耦不对称性： 前导与后随相位中心因处于阵列的不同边缘环境，其互耦状态存在镜像不对称，导致两个子阵的方向图与相位中心位置发生失配，降低杂波相关性。
保护带效应： 在有效阵列边缘设置被动端接的保护带单元，可显著改善互耦对称性，使对消性能提升 10 – 20dB，同时使实际相位中心更接近理论设计位置。
阵列规模影响：
- 行数： 对于大型阵列，俯仰向行数变化对对消性能影响甚微；
- 列数/长度： 航迹向长度增加可提升对消性能，因大孔径削弱了边缘互耦的相对影响；
- 相位中心位移： 小型阵列中，位移增大会显著降低对消性能（子阵孔径减小，边缘效应增强）；大型阵列对该参数相对不敏感。
单元类型差异： 单极子单元在法向具有天然零点，有利于抑制天底强杂波，但宽角扫描时互耦引起的性能下降较偶极子更为显著；偶极子覆盖范围宽，但在小阵列条件下与单极子的优劣关系会随阵列尺寸发生交叉变化。