4. 海杂波特性分析与抑制方法

海杂波是雷达照射海面时由大量随机散射体回波叠加形成的干扰，它和目标回波在时域、频域、统计域上高度交叠，是海面目标检测中最主要的虚警来源。理解海杂波，核心不是记住某个公式，而是建立“统计分布—检测门限—自适应策略—实测效果”这条因果链：海况和雷达参数决定杂波幅度服从什么分布，分布决定 CFAR 采用哪种参考单元估计和哪种门限系数，门限策略又直接决定检测概率与虚警率能否同时满足要求。

一、海杂波的统计特性与模型选择

海杂波幅度随海况、擦地角、极化方式、雷达分辨率和驻留时间变化显著。低分辨率、低海况、大擦地角条件下，回波由大量独立小散射体叠加，幅度接近瑞利分布；当分辨率提高或海况增强时，海尖峰出现，尾部变厚，瑞利分布会严重低估大幅度杂波，导致虚警激增。工程上常用韦布尔分布、对数正态分布和 K 分布描述这类非均匀杂波，其中 K 分布由慢变的调制分量（反映大尺度海浪结构）和快变的散斑分量共同构成，对海尖峰的拟合较好。

选择模型时可用实测杂波数据做拟合检验：对幅度样本取累积分布，与候选分布的理论曲线对比，或用二阶、高阶矩估计形状参数。例如 K 分布的形状参数 ν 越小，尾部越厚、海尖峰越明显；ν 增大时逐渐退化为瑞利型。实际系统中不必追求最优拟合，而应关注“模型尾部是否足以覆盖强杂波”，因为 CFAR 的虚警率主要由分布尾部决定。

二、CFAR 检测原理

CFAR（恒虚警率）检测的基本思想是：不用固定门限，而是用被检测单元（CUT）周围参考单元的杂波功率估计背景电平，再乘以由虚警率要求确定的门限系数得到自适应门限。设参考单元平均功率估计为 Z，门限为 T = α·Z，则判决规则为：CUT 功率超过 T 判为目标，否则判为杂波。门限系数 α 由参考单元数 N 和目标虚警率 Pfa 共同决定，对指数型杂波，单元平均 CFAR 有 α = N·(Pfa^(−1/N) − 1)。

参考窗的排布方式决定了 CFAR 的抗干扰能力。单元平均 CFAR（CA-CFAR）在均匀杂波中最优，但参考窗内若混入目标或强杂波边缘，背景估计会被抬高，造成遮蔽或虚警。为此有取大（GO-CFAR）、取小（SO-CFAR）、有序统计（OS-CFAR）等变体：GO-CFAR 取两侧均值较大者，适合杂波边缘，虚警控制好但边缘目标易漏检；SO-CFAR 取较小者，适合多目标环境但虚警偏高；OS-CFAR 对参考单元排序后取第 k 个值，抗多目标能力强，代价是排序运算量。

三、自适应阈值选择的关键参数

门限系数 α 的取值直接体现“虚警率优先还是检测率优先”。以 N=16、Pfa=10⁻⁴ 为例，CA-CFAR 的 α ≈ 16×(10⁻⁴^(−1/16) − 1) ≈ 16×(1.778 − 1) ≈ 12.4，即门限约为背景平均功率的 12.4 倍。若把 Pfa 放宽到 10⁻³，α 降至约 8.9，检测率提高但虚警增多。实际系统还要考虑参考单元数：N 太小，估计方差大，门限抖动；N 太大，跨越不同杂波区域，背景估计失真。工程上常按距离分辨率与海况选取 N 在 16～32 之间，并在杂波边缘处切换到 GO-CFAR 或加保护单元。

另一个关键参数是参考窗与保护单元的布局。保护单元用于隔离目标能量泄漏，避免目标自身抬高背景估计造成自遮蔽。若目标展宽为 3 个距离单元，保护单元至少取 2～3 个。此外，可引入杂波图（clutter map）做慢变背景记忆，对海尖峰区域提高门限，对平静海面降低门限，实现空域和时域的双重自适应。

四、海杂波抑制实战对比

下面给出一个教学示例（参数为合理假设，非实测数据）：某 X 波段雷达距离分辨率 6 m，擦地角约 1°，海况 3 级，在 200 个距离单元上注入目标（信杂比约 12 dB）与海尖峰。分别用固定门限、CA-CFAR（N=16，Pfa=10⁻⁴）、OS-CFAR（N=16，k=12）处理同一组回波，观察检测结果。

判读结果时不能只看虚警数，还要看目标是否在杂波边缘被“吃掉”。若发现目标附近门限异常抬高，应检查保护单元是否足够、参考窗是否跨越了强杂波区；若虚警集中在海尖峰处，说明分布尾部建模不足，可改用 OS-CFAR 或引入杂波图提高该区域门限。最终以“给定虚警率下检测概率是否达标、边缘目标是否丢失”作为抑制方法是否有效的判断依据，而不是单纯追求画面干净。