4. 干扰信号建模：噪声调频干扰建模、距离欺骗干扰建模、速度欺骗干扰建模

各位好，我是老张。今天咱们聊聊干扰信号建模。说实话，这是整个火控抗干扰课程里最“实在”的一章。你想想看，连对手的干扰长什么样都不知道，你怎么去对抗它？

我个人习惯是，先搞清楚干扰的数学模型，再谈怎么对付它。这就像打仗前先摸清敌人的装备，对吧？

4.1 噪声调频干扰建模

噪声调频干扰，说白了就是让干扰信号的频率跟着噪声随机变化。它是最常见的压制式干扰之一。我在项目中遇到过好几次，对方一开这个干扰，雷达屏幕上直接一片雪花，啥也看不见。

它的数学模型长这样：

s_j(t) = A_j * cos(2π * f_c * t + 2π * K_f * ∫_0^t n(τ) dτ + φ_0)

这里：

• A_j：干扰信号幅度，决定了干扰功率
• f_c：载波频率，一般对准雷达工作频率
• K_f：调频斜率，控制频率偏移范围
• n(τ)：零均值高斯白噪声，方差为σ²
• φ_0：初始相位

为什么要用噪声调频？因为它的瞬时频率一直在变，雷达的窄带滤波器很难锁定它。嗯，这里要注意：调频斜率K_f越大，干扰带宽越宽，压制效果越好，但功率谱密度会下降。

下面给一段Python代码，方便大家直接跑仿真：

import numpy as np
import matplotlib.pyplot as plt

def noise_fm_jamming(fs, T, fc, Kf, sigma, Aj=1.0):
    """
    生成噪声调频干扰信号
    fs: 采样率 (Hz)
    T: 信号时长 (s)
    fc: 载波频率 (Hz)
    Kf: 调频斜率 (Hz/V)
    sigma: 噪声标准差
    Aj: 幅度
    """
    t = np.arange(0, T, 1/fs)
    n = sigma * np.random.randn(len(t))  # 高斯白噪声
    phase = 2 * np.pi * Kf * np.cumsum(n) / fs  # 积分求相位
    s_j = Aj * np.cos(2 * np.pi * fc * t + phase)
    return t, s_j

# 示例：采样率10MHz，时长1ms，载波1MHz
t, s_j = noise_fm_jamming(fs=10e6, T=1e-3, fc=1e6, Kf=1e5, sigma=1.0)
plt.plot(t[:500], s_j[:500])
plt.title('噪声调频干扰信号（时域片段）')
plt.xlabel('时间 (s)')
plt.ylabel('幅度')
plt.show()

4.2 距离欺骗干扰建模

距离欺骗干扰，说白了就是让雷达以为目标在假的距离上。它属于欺骗式干扰，比压制式干扰更阴险。你想想看，雷达明明看到目标在10公里处，结果实际在5公里，这仗还怎么打？

它的核心思路是：干扰机收到雷达脉冲后，延迟一段时间再转发回去。延迟量Δτ决定了假目标的距离偏移：

R_false = R_true + c * Δτ / 2

其中c是光速。如果Δτ为正，假目标比真目标远；Δτ为负，假目标更近（需要干扰机提前转发，技术上更难）。

我建议把距离欺骗干扰分成两种类型：

• 固定延迟欺骗：延迟量固定，产生一个固定的假目标
• 拖引式欺骗：延迟量逐渐增大，让假目标慢慢远离真目标，把雷达的跟踪波门拖走

拖引式欺骗的数学模型更复杂一些：

Δτ(t) = Δτ_0 + a * t  (0 ≤ t ≤ T_drag)

这里a是拖引加速度，T_drag是拖引时间。我在项目中遇到过最狠的一次，对方用0.5μs/s的加速度拖引，我们的跟踪波门硬生生被拖走了2公里才反应过来。

代码实现也不复杂：

def range_deception_jamming(pulse_timing, delay, fs):
    """
    生成距离欺骗干扰信号
    pulse_timing: 原始脉冲到达时间 (s)
    delay: 延迟量 (s)
    fs: 采样率
    """
    # 假设原始脉冲是矩形脉冲
    pulse_width = 1e-6  # 1微秒脉宽
    t_j = pulse_timing + delay  # 干扰脉冲时间
    # 生成干扰脉冲（这里简化处理）
    n_samples = int(pulse_width * fs)
    jamming_pulse = np.ones(n_samples)
    return t_j, jamming_pulse

# 拖引式：延迟随时间线性增加
t_drag = np.linspace(0, 5, 100)  # 5秒拖引
delay_drag = 0.1e-6 + 0.2e-6 * t_drag  # 加速度0.2μs/s
plt.plot(t_drag, delay_drag * 1e6)
plt.title('拖引式距离欺骗：延迟随时间变化')
plt.xlabel('时间 (s)')
plt.ylabel('延迟 (μs)')
plt.show()

4.3 速度欺骗干扰建模

速度欺骗干扰，说白了就是让雷达以为目标在错误的速度上。它利用的是多普勒频移原理。你想想看，如果雷达以为目标以300m/s的速度接近，实际只有100m/s，那拦截弹的引导头就全偏了。

速度欺骗的核心是改变干扰信号的多普勒频率。真实目标的多普勒频移是：

f_d = 2 * v_r / λ

其中v_r是径向速度，λ是雷达波长。干扰机通过调制信号频率，让雷达解算出错误的速度：

f_j = f_d + Δf

Δf就是速度欺骗的频率偏移量。对应的假速度是：

v_false = v_true + λ * Δf / 2

我建议把速度欺骗也分成两种：

• 固定速度欺骗：Δf固定，产生一个固定速度的假目标
• 速度拖引：Δf逐渐变化，让假目标的速度慢慢偏离真目标

速度拖引的数学模型：

Δf(t) = Δf_0 + b * t  (0 ≤ t ≤ T_drag_v)

b是频率变化率，单位Hz/s。我在项目中遇到过最头疼的一次，对方用500Hz/s的速度拖引，我们的速度跟踪环路花了3秒才重新锁定，那3秒里导弹已经飞偏了。

代码实现：

def velocity_deception_jamming(t, fc, v_true, delta_f, Aj=1.0):
    """
    生成速度欺骗干扰信号
    t: 时间向量
    fc: 雷达载频
    v_true: 真实径向速度
    delta_f: 速度欺骗频偏
    """
    lambda_c = 3e8 / fc  # 波长
    f_d_true = 2 * v_true / lambda_c  # 真实多普勒
    f_j = f_d_true + delta_f  # 干扰多普勒
    s_j = Aj * np.cos(2 * np.pi * (fc + f_j) * t)
    return s_j

# 示例：速度拖引，频偏从0线性增加到2kHz
t = np.arange(0, 5, 1/10e6)  # 5秒，10MHz采样
delta_f_t = 0 + 400 * t  # 400Hz/s变化率
s_j = velocity_deception_jamming(t, fc=10e9, v_true=200, delta_f=delta_f_t)
# 这里只展示原理，实际应用中需要和雷达脉冲同步

好了，三种干扰建模就讲到这里。噪声调频是压制，距离欺骗和速度欺骗是诱骗。你想想看，如果对手把这三种组合起来用，那才叫头疼。下一章咱们就聊聊怎么对付这些干扰。