3. 多普勒效应与回波信号建模

星载雷达与地面静止雷达最本质的差别，在于雷达本身以约 7 km/s 的速度在轨道上运动。这一运动叠加在电磁波的传播之上，使回波信号在时间延迟之外，还携带了一个由平台运动引入的相位调制项。理解这一相位调制的来源、形式与频谱后果，是后续脉冲压缩、多普勒滤波和运动补偿的物理基础。本节从发射信号出发，逐步推导到回波频谱结构，并给出可用于仿真的建模方法。

一、发射信号与传播时延

设雷达发射线性调频（LFM）脉冲串，第 m 个脉冲的基带形式可写为

s_tx(t) = rect(t/T_p) · exp(jπK_r t²) · exp(j2πf_c t)
其中：T_p 为脉宽，K_r = B/T_p 为调频斜率，B 为带宽，f_c 为载频

脉冲以重复周期 PRI 依次发射，慢时间 t_m = m·PRI（m = 0,1,…,M−1）标记脉冲序号。电磁波经目标反射后返回，单程传播时延为

τ(t) = R(t) / c
R(t) = |R_sat(t) − R_tgt|  （斜距，随慢时间变化）

其中 R_sat(t) 为卫星在惯性系中的位置矢量，R_tgt 为目标位置。由于平台运动，R(t) 在合成孔径时间内并非常数，而是随慢时间缓变。这一缓变正是多普勒效应的几何来源。

二、平台运动引起的相位调制

接收回波可表示为发射信号的延迟副本：

s_rx(t) = A · s_tx(t − τ(t))
        = A · rect((t−τ)/T_p) · exp(jπK_r (t−τ)²) · exp(j2πf_c (t−τ))

关键在最后一项 exp(−j2πf_c τ(t))。将 τ(t) = R(t)/c 代入，得到由斜距变化引起的相位项：

φ(t) = −2πf_c · R(t) / c = −(2π/λ) · R(t)

这就是平台运动引起的相位调制。对慢时间求导即得多普勒频率：

f_d(t) = (1/2π) · dφ/dt = −(2/λ) · dR(t)/dt = −(2/λ) · v_r(t)

其中 v_r = dR/dt 为径向速度。符号约定：目标靠近雷达时 R 减小，v_r < 0，f_d > 0；目标远离时 f_d < 0。全文与图示均采用此约定。

在合成孔径时间 T_sar 内，若目标相对雷达的径向速度近似线性变化，可将 R(t) 在波束中心时刻 t_0 处作泰勒展开：

R(t) ≈ R_0 + v_r·(t−t_0) + ½·a_r·(t−t_0)²
φ(t) ≈ −(2π/λ)[R_0 + v_r·(t−t_0) + ½·a_r·(t−t_0)²]

一次项对应多普勒中心频率 f_dc = −2v_r/λ，二次项对应多普勒调频率 K_a = −2a_r/λ。对于正侧视条带模式，v_r 近似为零，f_dc ≈ 0；对于斜视或前视，f_dc 可达数千赫兹，必须精确估计并补偿。

三、回波频谱结构

回波在快时间维和慢时间维上具有不同的频谱特征。快时间维上，LFM 信号经匹配滤波后形成 sinc 型距离压缩峰，主瓣宽度约 c/(2B)；慢时间维上，相位调制 φ(t) 使同一距离单元内的回波在频域展开为多普勒谱。

以典型星载 SAR 参数为例：λ = 0.056 m（C 波段），v_sat = 7.6 km/s，斜距 R_0 = 850 km，方位波束宽度 θ_az = 0.3°。则合成孔径时间 T_sar ≈ λR_0/(v_sat·θ_az·π/180) ≈ 2.1 s，多普勒带宽 B_d ≈ 2v_sat·θ_az/λ ≈ 1.4 kHz。这意味着慢时间采样率（PRF）必须大于 B_d，否则多普勒谱将发生混叠。

回波频谱的完整结构可分解为：

• 距离频率：由 LFM 带宽 B 决定，范围 [−B/2, B/2]；
• 方位频率：由多普勒调频率 K_a 和合成孔径时间决定，范围 [f_dc − B_d/2, f_dc + B_d/2]；
• 耦合项：当存在斜视或平台加速度时，距离频率与方位频率产生线性耦合，表现为距离走动（range walk）和距离弯曲（range curvature）。

在二维频域中，回波信号可近似表示为：

S(f_r, f_a) ≈ A · rect(f_r/B) · rect((f_a−f_dc)/B_d)
              · exp(−jπ f_r²/K_r)
              · exp(−jπ (f_a−f_dc)²/K_a)
              · exp(−j4πR_0/λ)
              · exp(j2π f_r · f_a / K_coupling)

最后一项即距离-方位耦合项，K_coupling 与斜视角和平台速度有关。忽略该项会导致距离压缩后目标位置偏移，在宽测绘带或高分辨率场景下不可接受。

四、建模与仿真实现

在工程实践中，回波信号建模通常采用“几何-相位-采样”三步法。以下给出一个可复现的仿真流程示例（参数为合理假设，非实测数据）：

% 参数设置（示例）
c = 3e8; fc = 5.4e9; lambda = c/fc;
B = 30e6; Tp = 10e-6; Kr = B/Tp;
PRI = 1e-3; M = 512; PRF = 1/PRI;
vs = 7600; R0 = 850e3; theta_az = 0.3*pi/180;
Tsar = lambda*R0/(vs*theta_az);  % ≈ 2.1 s
Ka = -2*vs^2/(lambda*R0);        % 多普勒调频率

% 慢时间序列
tm = (0:M-1)*PRI - Tsar/2;
R = sqrt(R0^2 + (vs*tm).^2);     % 斜距历史
tau = 2*R/c;                     % 双程时延

% 快时间采样
fs = 2*B; dt = 1/fs;
Nf = round(Tp*fs);
tf = (0:Nf-1)*dt;

% 生成回波（单目标）
for m = 1:M
    s_tx = exp(1j*pi*Kr*tf.^2);
    s_rx(m,:) = exp(-1j*4*pi*R(m)/lambda) .* ...
                [zeros(1,round(tau(m)*fs)), s_tx(1:end-round(tau(m)*fs))];
end

仿真结果可通过二维 FFT 观察频谱结构：距离维应呈现 LFM 的矩形谱，方位维应呈现由 K_a 决定的线性调频谱。若在方位谱中观察到谱峰偏离零频，则说明存在多普勒中心偏移，需在后续处理中估计并补偿。

五、关键参数关系与影响

从表中可见，星载条件下 B_d 通常在 1～3 kHz 量级，而 PRF 受距离模糊和方位模糊的双重约束，往往在 1～2 kHz 之间。这意味着多普勒谱经常处于临界采样状态，对 PRF 的选择和后续多普勒中心估计精度提出了很高要求。若 f_dc 估计误差超过 B_d/2，方位压缩将出现信噪比损失和主瓣偏移。

上图将发射信号、传播时延、相位调制与二维频谱结构串联起来，并标注了各环节的关键参数与符号约定。实际建模时，应首先确认几何关系（斜距历史 R(t)），再计算相位 φ(t)，最后通过快时间采样和慢时间采样生成回波矩阵。频谱分析时，距离维关注带宽与耦合项，方位维关注多普勒中心、带宽和调频率。只有这三者同时正确，后续的脉冲压缩和运动补偿才有可靠输入。