3. 海杂波空间相关性分析与相关时间估计

相关性与参考单元失效的物理链条

海杂波在距离维和多普勒维都不是白噪声：相邻距离单元共享同一片照射海面，相邻脉冲之间海面起伏尚未完全变化，因此回波存在明显相关性。CFAR 检测器默认参考单元与待检测单元“独立同分布”，当参考单元与待检测单元相关性过强时，参考单元估计的杂波功率会跟随待检测单元一起起伏，门限被“拉平”，虚警率不再受控。本节用实测回波数据计算两个相关系数：距离向相邻单元相关系数 ρ_r(Δr) 与脉冲向相邻脉冲相关系数 ρ_p(Δn)，并据此估计相关距离 L_c 与去相关时间 T_c，最后判断当前 CFAR 参考窗是否落在相关半径之内。

数据组织与相关系数计算

取一段海杂波回波复数矩阵 X，维度为 N_r 个距离单元 × N_p 个脉冲。为消除目标与强散射点干扰，先剔除幅度超过局部中值 6 倍的单元，再对每个距离单元做脉冲向去均值。相关系数用复数形式计算，避免只取幅度丢失相位信息：

import numpy as np

def corr_complex(a, b):
    # a, b: 一维复数序列，已去均值
    num = np.sum(a * np.conj(b))
    den = np.sqrt(np.sum(np.abs(a)**2) * np.sum(np.abs(b)**2))
    return np.abs(num / den) if den > 0 else 0.0

def range_corr(X, max_lag=16):
    # X: (Nr, Np) 复数矩阵，按距离单元逐行
    Xc = X - X.mean(axis=1, keepdims=True)
    rho = []
    for lag in range(1, max_lag + 1):
        vals = [corr_complex(Xc[i, :], Xc[i + lag, :])
                for i in range(Xc.shape[0] - lag)]
        rho.append(np.mean(vals))
    return np.array(rho)   # rho[k-1] 对应 lag=k

def pulse_corr(X, max_lag=32):
    Xc = X - X.mean(axis=1, keepdims=True)
    rho = []
    for lag in range(1, max_lag + 1):
        vals = [corr_complex(Xc[:, j], Xc[:, j + lag])
                for j in range(Xc.shape[1] - lag)]
        rho.append(np.mean(vals))
    return np.array(rho)

示例数据取 N_r=512、N_p=256，脉冲重复频率 PRF=1000 Hz，距离分辨率 ΔR=15 m。计算得到距离向相关系数在 lag=1 时约 0.82，lag=4 时约 0.45，lag=8 时约 0.18；脉冲向相关系数在 lag=1 时约 0.91，lag=6 时约 0.52，lag=20 时约 0.15。这些数值属于教学示例，实际使用时须用现场数据重算。

相关距离与去相关时间的估计

工程上常用“相关系数首次降到 1/e≈0.37 以下”的滞后量作为相关长度。对距离向曲线做线性插值求 ρ_r=0.37 对应的 lag，得 lag_r≈5.1，则相关距离 L_c=5.1×15 m≈77 m。对脉冲向曲线插值得 lag_p≈9.4，去相关时间 T_c=9.4/1000 Hz≈9.4 ms。若海况升高或雷达频率降低，T_c 会变长，相关距离也会随距离分辨率与擦地角变化。

参考窗选取是否失效的判定

把估计结果与当前 CFAR 配置对照。假设当前采用 CA-CFAR，单侧参考单元数 N_ref=8，保护单元 2 个，参考窗在距离维单侧覆盖 8×15 m=120 m，大于 L_c≈77 m，距离向基本满足独立性；但若系统在脉冲维做非相干积累后再检测，积累脉冲数为 4，而 T_c 对应约 9.4 个脉冲，4 个脉冲内的杂波仍高度相关，积累后残余相关性使参考单元方差估计偏小，虚警率抬升。此时应把积累脉冲数提高到 ≥10，或在 CFAR 前做脉冲向白化，或改用 GO-CFAR/OS-CFAR 降低相关参考单元对门限的拖拽。

可复用的判读清单

• 数据段长度 ≥ 512 距离单元 × 256 脉冲，剔除幅度超局部中值 6 倍的单元。
• 分别计算 ρ_r(1…16) 与 ρ_p(1…32)，记录首次低于 0.37 的滞后量并线性插值。
• 由 L_c=lag_r×ΔR、T_c=lag_p/PRF 换算物理量。
• 对照 CFAR 单侧参考长度与积累脉冲数：任一小于对应相关长度，即判定相关性过强、虚警控制存在失效风险。
• 处置优先级：先增大积累脉冲数至 T_c 对应值以上，再检查参考窗长度，最后考虑换 CFAR 结构。