2. 星载雷达几何与运动模型
星载雷达与机载、地基雷达最本质的差别,在于平台本身处在高速绕地运动中,而目标(地面或海面散射体)又随地球自转一起运动。多普勒频率不再只是“平台相对目标径向速度”的简单函数,而是卫星轨道运动、地球自转、波束指向和斜距历程共同作用的结果。要把多普勒处理做对,第一步就是把这几者的解析关系理清楚。
一、坐标系与参照约定
本节统一采用地心惯性坐标系(ECI)描述卫星轨道运动,采用地心地固坐标系(ECEF)描述地面目标位置。卫星在 ECI 中的位置记为 Rs(t),地面目标在 ECEF 中的位置记为 Rt,转换到 ECI 时需乘以地球自转矩阵。斜距定义为
R(t) = | R_s(t) - M_eci←ecef(t) · R_t |
其中 Meci←ecef(t) 是随时间变化的地球自转旋转矩阵,自转角速度 ωe ≈ 7.292115×10⁻⁵ rad/s。约定:卫星速度方向为航向,波束指向与航向的夹角为方位角,与星下点方向的夹角为俯仰角(下视角)。全文正负号以“目标接近卫星时多普勒为正”为约定。
二、卫星轨道与星下点运动
近圆轨道下,卫星在 ECI 中的位置可近似用轨道六根数描述:半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 ν。对星载 SAR 常用的近圆太阳同步轨道,典型参数为:轨道高度 H ≈ 700 km,倾角 i ≈ 98.2°,偏心率 e ≈ 0.001,轨道周期 T ≈ 99 min,卫星地速约 7.6 km/s。
卫星在 ECI 中的速度 Vs 由轨道动力学给出,其星下点在地球表面的移动速度约为
V_g ≈ V_s · R_e / (R_e + H)
代入 Re = 6371 km、H = 700 km,可得 Vg ≈ 6.7 km/s。但星下点在地固系中的实际移动还要叠加地球自转,因此地面轨迹相对地球的等效速度会随纬度变化,这正是后续多普勒中心频率随纬度漂移的根源之一。
三、地球自转的贡献
地球自转使地面目标在 ECI 中具有切向速度
V_e = ω_e × R_t
其大小为 ωe·Re·cos(纬度) ≈ 465·cos(纬度) m/s。在赤道处约 465 m/s,在中纬度约 300 m/s。这个速度虽然远小于卫星速度,但它与波束指向的耦合会显著改变多普勒中心频率,尤其在斜视或大下视角时不可忽略。工程上常把地球自转效应折算为“等效斜视角偏移”,用于修正多普勒中心估计。
四、斜距历程与多普勒频率
斜距历程 R(t) 是星载雷达多普勒处理的核心量。对近圆轨道,R(t) 在波束驻留时间内可展开为关于慢时间 t 的二次多项式:
R(t) ≈ R_0 + V_r·t + (1/2)·a_r·t²
其中 R0 为波束中心时刻斜距,Vr 为径向速度,ar 为径向加速度。对多普勒频率求导可得
f_D(t) = -2/λ · dR/dt = -2/λ · (V_r + a_r·t)
因此多普勒频率在慢时间上呈线性调频(LFM)形式,调频率为
K_a = -2·a_r/λ
这就是方位向匹配滤波的调频率来源。以 C 波段 λ = 5.6 cm、ar ≈ 0.5 m/s² 为例,Ka ≈ -17.9 Hz/s,对应方位向合成孔径时间约 1 s 时多普勒带宽约 18 Hz。
五、波束指向与多普勒中心
多普勒中心频率 fDC 由波束中心视线方向与平台速度的夹角决定。对正侧视,fDC ≈ 0;对斜视,fDC 随斜视角线性增大。地球自转会使地面目标相对卫星产生额外切向速度,等效于引入一个附加斜视角,从而把 fDC 推离零频。工程上常用“多普勒中心估计—去斜—方位向 FFT”流程处理,其中去斜量即由 fDC 决定。
下面给出一个可直接用于教学示例的参数表,用于估算某近圆轨道下的多普勒中心频率。表中数值为合理假设的示例参数,非实测结果。
| 参数 | 符号 | 示例取值 | 说明 |
|---|---|---|---|
| 轨道高度 | H | 700 km | 近圆太阳同步轨道 |
| 卫星速度 | V_s | 7.6 km/s | ECI 中速度 |
| 波长 | λ | 5.6 cm | C 波段 |
| 下视角 | θ | 35° | 波束中心与星下点夹角 |
| 斜视角 | θ_sq | 3° | 波束中心偏离正侧视 |
| 目标纬度 | φ | 30° N | 用于计算 V_e |
| 地球自转速度 | V_e | ≈ 403 m/s | 465·cos30° |
| 多普勒中心 | f_DC | ≈ 1.2 kHz | 由 V_s、V_e、θ、θ_sq 合成 |
六、从几何到处理的映射
把上述关系串起来,可以得到一条清晰的映射链:轨道六根数 → 卫星位置与速度 → 星下点与波束指向 → 斜距历程 R(t) → 多普勒频率 fD(t) → 多普勒中心 fDC 与调频率 Ka → 方位向匹配滤波参数。实际实现时,通常先用轨道预报得到 Rs(t) 与 Vs(t),再对每个距离门计算 R(t) 并求导,最后把 fDC 和 Ka 送入方位向处理链。
需要特别注意的是,地球自转不仅改变 fDC,还会使等多普勒锥发生旋转,导致同一距离门内不同方位位置的多普勒历史不再严格对称。对高分辨率或大斜视场景,必须把地球自转项完整保留在斜距模型中,否则方位向聚焦会出现明显展宽或旁瓣抬升。