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”的频域方法提高效率。

教学示例:取 f_c = 9.6 GHz、B = 20 MHz、T_p = 20 μs、PRF = 2 kHz、f_s = 25 MHz,点目标初始斜距 850 km、径向速度 7.2 km/s。按上述参数生成 256 个脉冲的回波,可在距离—慢时间二维矩阵中观察到一条近似直线的距离走动轨迹,其斜率对应径向速度。

二、分布式面目标回波:从单点到面元叠加

分布式面目标(如海面、地面、云雨区)不能再用单个散射中心描述,而应把照射区域划分为若干面元,每个面元视为一个具有独立 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)计算卫星在地固系或惯性系中的位置,再转换到天线相位中心,最后投影到波束指向方向得到斜距历史。

星载雷达回波建模:几何、面元与多普勒历史 地球表面(分布式面目标区域) 卫星平台 P_s(t) 天线波束 目标点 P_t R(t) = |P_s(t) − P_t| v_s ≈ 7.5 km/s 面元划分:σ_k = σ⁰·ΔA_k·G_k²/R_k⁴ 多普勒历史 f_d(t) 慢时间 t f_d 近似线性调频 回波生成链 ① 轨道/姿态 → P_s(t) ② 斜距 R(t) → τ(t) ③ LFM 延迟 + 相位 ④ 面元相干叠加 距离—慢时间矩阵 快时间 × 脉冲序号

图中给出了从卫星位置、斜距历史到面元叠加和距离—慢时间矩阵的完整关系。实际仿真时,建议把几何计算与信号生成解耦:先由轨道模型生成每个脉冲的 R_n 和 f_d,n,再统一送入回波合成模块。这样既便于替换轨道参数,也便于对点目标和面目标复用同一套相位计算。

四、仿真参数与结果判读

下面给出一组可直接用于教学仿真的参数配置,用于生成点目标与面目标混合场景的回波数据:

参数取值说明
载频 f_c9.6 GHzX 波段
带宽 B20 MHz距离分辨率约 7.5 m
脉宽 T_p20 μs调频斜率 1×10¹² Hz/s
PRF2 kHz脉冲数 256
采样率 f_s25 MHz过采样系数 1.25
轨道高度700 km地速约 7.5 km/s
点目标 RCS10 m²初始斜距 850 km
面目标 σ⁰−20 dB面元 4 m × 4 m

生成回波后,可按以下顺序判读:先对单脉冲做距离向脉冲压缩,检查点目标是否出现在预期距离门,主瓣宽度是否接近 c/(2B);再沿慢时间观察峰值位置是否随 R_n 移动,移动斜率应与径向速度一致;最后对面目标区域做方位向 FFT,观察多普勒频谱是否展宽到理论值 B_d ≈ 2v_t·θ_az/λ。若点目标压缩后旁瓣异常抬高,通常说明分数延迟实现有误或相位项符号不一致;若面目标回波在距离向出现周期性条纹,则可能是面元网格与采样率发生了周期性耦合。

注意:多普勒频率的符号取决于径向速度的定义方向。本文约定目标靠近雷达时 v_r > 0、f_d > 0。若采用相反约定,相位项 exp(-j4πR/λ) 中的符号需同步调整,否则脉冲压缩后的多普勒方向会反转。

把上述三部分组合起来,就得到星载雷达回波仿真的基本闭环:几何模型提供斜距和多普勒历史,点目标模型验证相位与延迟的正确性,面目标模型提供统计起伏和展宽特性,最终输出可用于脉冲压缩、成像和检测算法验证的复基带数据矩阵。