4. 海杂波多普勒谱特征提取与 Bragg 峰识别
为什么先看多普勒谱，而不是先看幅度
海面回波在慢时间维上并非白噪声：当雷达照射到满足 Bragg 条件的海浪分量时，回波会在多普勒谱上形成一对关于零频近似对称的尖峰，即 Bragg 峰。对 X 波段、擦地角较小的对海雷达，Bragg 峰的多普勒频率可由深水重力波色散关系近似给出：
fB ≈ ±(2vr/λ)·cosψ，vr = √(gλB/(2π))，λB = λ/(2sinθg)
其中 λ 为雷达波长，θg 为擦地角，ψ 为 Bragg 波传播方向与雷达视线夹角，g≈9.81 m/s²。可见：雷达频率越高（λ 越小），Bragg 峰偏移越大；擦地角越小，Bragg 波长越长、相速度越高，偏移也越大。这正是判断“谱上某个尖峰是海杂波还是目标”的第一层依据——目标回波的多普勒位置由径向速度决定，与海况、擦地角无固定色散关系。
谱估计的实现要点
对某一距离单元取 N 个脉冲的复回波 x[n]，先去除直流分量（海杂波零频附近能量极强，不去直流会淹没 Bragg 峰），再加窗抑制旁瓣，然后做 FFT 并取功率谱：
import numpy as np
def doppler_spectrum(x, prf, win="hann"):
x = np.asarray(x, dtype=complex)
x = x - x.mean() # 去直流，避免零频泄漏
N = x.size
w = np.hanning(N) if win == "hann" else np.ones(N)
X = np.fft.fftshift(np.fft.fft(x * w))
P = np.abs(X) ** 2 / (N * np.sum(w ** 2)) # 功率谱密度归一化
f = np.fft.fftshift(np.fft.fftfreq(N, d=1.0 / prf))
return f, P
def bragg_freq(lam, grazing_deg, psi_deg=0.0):
g = 9.81
th = np.deg2rad(grazing_deg)
lam_b = lam / (2.0 * np.sin(th))
vr = np.sqrt(g * lam_b / (2.0 * np.pi))
return 2.0 * vr / lam * np.cos(np.deg2rad(psi_deg))关键参数：脉冲数 N 决定多普勒分辨率 Δf = PRF/N；PRF 决定无模糊多普勒范围 ±PRF/2。若目标径向速度对应的多普勒接近 PRF/2，会发生多普勒模糊，此时不能仅凭谱峰位置下结论。
Bragg 峰识别与误判判断
识别流程不是简单找最大值，而是三步：①在理论 fB 附近（考虑 ±20% 容差与海流引起的整体平移）搜索局部极大；②检查正负频是否成对出现且幅度量级相近（Bragg 散射的对称性）；③检查峰宽与谱对称性——Bragg 峰通常较窄且成对，目标回波常为单侧、可能更窄或更宽，且不满足色散关系。
海况与雷达频率对 Bragg 峰偏移的影响
用上面的 bragg_freq 做一组教学示例计算（假设值，非实测）：取 X 波段 λ=0.03 m、擦地角 5°，则 λB≈0.172 m，vr≈0.52 m/s，fB≈69 Hz；若改为 S 波段 λ=0.10 m、同擦地角，λB≈0.574 m，vr≈0.95 m/s，fB≈19 Hz。可见频率升高约 3.3 倍，Bragg 偏移也近似同比例增大。海况的影响体现在谱宽与峰强：高海况下连续谱底抬升、Bragg 峰相对减弱甚至被淹没，此时单看峰高容易漏判。
| 观察项 | Bragg 峰特征 | 目标回波特征 | 判读结论 |
|---|---|---|---|
| 正负频对称性 | 成对出现，幅度量级相近 | 通常单侧 | 成对→倾向海杂波 |
| 与理论 fB 偏差 | 在 ±20% 容差内（扣除海流平移） | 无固定关系 | 超容差→需排查目标或干扰 |
| 峰宽 | 较窄，随海况展宽 | 由目标起伏与驻留时间决定 | 结合展宽趋势判断 |
| 随擦地角变化 | θg 减小，fB 增大 | 不随 θg 规律变化 | 多距离单元对比可区分 |
实操判断路径
- 选一个明显只有海面、无已知目标的距离单元，取 N=256 脉冲、记录 PRF，调用
doppler_spectrum得到 f 与 P。 - 用当前 λ 与擦地角算理论 fB，在 ±20% 区间内找局部极大，记录正负峰位置，计算峰对间距并与 2fB 比较。
- 若峰对间距与 2fB 吻合且正负峰幅度比在合理范围（如 0.5~2），判为 Bragg 峰；若某距离单元出现单侧强峰且不满足色散关系，标记为疑似目标，转入检测流程复核。
- 改变擦地角（不同距离单元）重复上述步骤，观察 fB 是否按理论趋势变化，作为二次确认。