4. 星载雷达回波建模与仿真
星载雷达的接收信号不是简单的“发射脉冲延迟后返回”,而是发射波形、目标散射、星地几何、天线方向图与接收机噪声共同作用的结果。回波建模的目标,是在给定轨道、姿态、波束指向和目标场景的条件下,生成可用于脉冲压缩、成像或检测算法验证的复基带信号。本节围绕点目标、分布式面目标和星地几何与多普勒历史三条主线展开,并给出可直接在仿真中复用的建模流程与参数示例。
一、点目标回波:从斜距到复基带
点目标模型把目标视为一个位于特定斜距、具有特定雷达截面积(RCS)的散射中心。对线性调频(LFM)脉冲雷达,第 n 个发射脉冲的复基带回波可写为
s_r(t, n) = A_n · exp(j·π·K_r·(t - τ_n)²) · exp(-j·4π·R_n/λ) · rect((t - τ_n)/T_p)
A_n = sqrt(P_t · G_t · G_r · λ² · σ / ((4π)³ · R_n⁴ · L))
其中 τ_n = 2R_n/c 为双程时延,K_r = B/T_p 为调频斜率,R_n 为第 n 个脉冲时刻的斜距,λ 为波长,σ 为目标 RCS,L 为系统损耗。相位项 exp(-j4πR_n/λ) 是回波建模中最关键的部分:它把斜距的微小变化转化为相位变化,进而形成多普勒历史。
仿真实现时,通常先按快时间采样率 f_s 生成距离向时间轴,再对每个脉冲计算 R_n 和 τ_n,将 LFM 波形延迟并乘以相位项。需要注意两点:一是延迟通常不是采样间隔的整数倍,需用频域相位斜坡或插值实现分数延迟;二是 R_n 随脉冲变化,若直接按每个脉冲重新生成波形,计算量会随脉冲数线性增长,实际仿真中常采用“距离向 FFT—相位补偿—距离向 IFFT”的频域方法提高效率。
二、分布式面目标回波:从单点到面元叠加
分布式面目标(如海面、地面、云雨区)不能再用单个散射中心描述,而应把照射区域划分为若干面元,每个面元视为一个具有独立 RCS 和相位的点目标,总回波为所有面元回波的相干叠加:
s_face(t, n) = Σ_k A_k · exp(j·π·K_r·(t - τ_k(n))²) · exp(-j·4π·R_k(n)/λ)
σ_k = σ⁰ · ΔA_k · G_k² / R_k⁴
其中 σ⁰ 为归一化后向散射系数,ΔA_k 为面元面积,G_k 为该面元方向上的天线增益。与点目标相比,面目标回波的核心难点在于:面元数量大、相位随机性强、回波在距离和方位上都会展宽。若面元划分过粗,会丢失散射统计特性;划分过细,则计算量急剧上升。工程上常按“一个距离分辨单元内至少 4~8 个面元”的原则确定网格,并利用距离向和方位向的独立采样降低冗余。
面目标仿真的另一个关键是随机相位与幅度统计。对均匀分布的地面或海面,面元复散射系数可建模为复高斯随机变量,其幅度服从瑞利分布;对存在强散射体的场景,可叠加一个确定性分量,使幅度服从莱斯分布。这样生成的回波经过脉冲压缩后,距离向旁瓣和起伏特性才接近真实观测。
三、星地几何与多普勒历史建模
星载雷达的多普勒历史由卫星与目标之间的相对运动决定。设卫星位置为 P_s(t),目标位置为 P_t,则瞬时斜距 R(t) = |P_s(t) − P_t|,多普勒频率为
f_d(t) = -(2/λ) · dR(t)/dt
R(t) ≈ R_0 + v_r·t + (v_t²/(2R_0))·t²
其中 v_r 为径向速度,v_t 为切向速度。一次近似下多普勒频率近似线性变化,二次项对应距离弯曲。对低轨卫星,轨道高度约 500~800 km,地速约 7.5 km/s,若波束斜视,多普勒中心频率可达数十千赫,多普勒调频率约几千赫每秒。仿真中需要先由轨道根数或两行根数(TLE)计算卫星在地固系或惯性系中的位置,再转换到天线相位中心,最后投影到波束指向方向得到斜距历史。
图中给出了从卫星位置、斜距历史到面元叠加和距离—慢时间矩阵的完整关系。实际仿真时,建议把几何计算与信号生成解耦:先由轨道模型生成每个脉冲的 R_n 和 f_d,n,再统一送入回波合成模块。这样既便于替换轨道参数,也便于对点目标和面目标复用同一套相位计算。
四、仿真参数与结果判读
下面给出一组可直接用于教学仿真的参数配置,用于生成点目标与面目标混合场景的回波数据:
| 参数 | 取值 | 说明 |
|---|---|---|
| 载频 f_c | 9.6 GHz | X 波段 |
| 带宽 B | 20 MHz | 距离分辨率约 7.5 m |
| 脉宽 T_p | 20 μs | 调频斜率 1×10¹² Hz/s |
| PRF | 2 kHz | 脉冲数 256 |
| 采样率 f_s | 25 MHz | 过采样系数 1.25 |
| 轨道高度 | 700 km | 地速约 7.5 km/s |
| 点目标 RCS | 10 m² | 初始斜距 850 km |
| 面目标 σ⁰ | −20 dB | 面元 4 m × 4 m |
生成回波后,可按以下顺序判读:先对单脉冲做距离向脉冲压缩,检查点目标是否出现在预期距离门,主瓣宽度是否接近 c/(2B);再沿慢时间观察峰值位置是否随 R_n 移动,移动斜率应与径向速度一致;最后对面目标区域做方位向 FFT,观察多普勒频谱是否展宽到理论值 B_d ≈ 2v_t·θ_az/λ。若点目标压缩后旁瓣异常抬高,通常说明分数延迟实现有误或相位项符号不一致;若面目标回波在距离向出现周期性条纹,则可能是面元网格与采样率发生了周期性耦合。
把上述三部分组合起来,就得到星载雷达回波仿真的基本闭环:几何模型提供斜距和多普勒历史,点目标模型验证相位与延迟的正确性,面目标模型提供统计起伏和展宽特性,最终输出可用于脉冲压缩、成像和检测算法验证的复基带数据矩阵。