单基三通道SAR-GMTI原理
本文系统推导机载单基三通道合成孔径雷达地面动目标指示SAR‑GMTI的信号处理流程。处理步骤依次为1. 系统几何与等效相位中心模型2. LFM 回波信号生成与雷达方程3. 距离压缩匹配滤波4. 方位向傅里叶变换——进入距离多普勒域5. 距离多普勒域通道配准6. DPCA 杂波对消7. ATI 沿航迹干涉测速8. 二维 CFAR 检测与多普勒凹口9. 距离多普勒算法RDA背景成像10. 性能指标评估1. 系统几何与等效相位中心1.1 坐标定义飞机以恒定速度沿轴正方向平飞高度为。定义慢时间方位时间零时刻飞机位于。天线采用“一发三收”体制- 中间天线Ch2同时用作发射和接收其相位中心瞬时坐标- 前向天线Ch3仅接收坐标- 后向天线Ch1仅接收坐标。为沿航迹基线长度。仿真典型值雷达工作于正侧视条带模式斜视角。1.2 点目标斜距设地面点目标坐标为。目标到发射/接收天线的瞬时距离决定了回波时延和相位。发射路径均由 Ch2 完成定义零多普勒斜距零多普勒时刻上式可写为接收路径- Ch2收发共用接收距离同发射距离故双程斜距为- Ch3前向接收接收天线坐标为- Ch1后向接收1.3 等效相位中心EPC近似当基线远小于斜距时可将 Ch1 和 Ch3 分别等效为一个收发共用的单基雷达其相位中心沿航迹偏移。以 Ch3 为例考虑其斜距并作泰勒展开在合成孔径时间内故一阶项可近似为进一步将该一阶项吸收到根号内的时间偏移中类似地物理解释Ch3 的等效相位中心比 Ch2 超前秒Ch1 则滞后秒。三通道可视为三个沿航迹排列的虚拟单基雷达它们接收到的静止目标回波仅在慢时间上存在固定偏移。1.4 DPCA 条件偏置相位中心天线DPCA技术要求前后两副天线在慢时间上实现“自补偿”相邻脉冲间前一副天线移动到后一副天线的位置。若脉冲重复周期则平台在一个 PRT 内飞过的距离为。当成立时Ch1 和 Ch3 的等效相位中心恰好相差两个于是静止目标的斜距满足此时只要将 Ch1 和 Ch3 的信号在方位时间上平移一个 PRT 即可与 Ch2 完全对齐使静止杂波在配准后完全相同。2. LFM 回波信号模型2.1 发射波形雷达发射线性调频LFM脉冲串。单个脉冲的基带复数形式为其中为脉宽为调频率单位 Hz/s为信号带宽。载频对应的真实射频信号为。2.2 点目标基带接收信号目标散射后回波被各通道接收。经下变频至基带第通道的信号为式中-为包含天线方向图、传播衰减和目标雷达截面积的幅度因子-为随机起始相位对各通道相同-为光速。2.3 天线方向图与雷达方程方向图模型采用均匀口径的 sinc² 近似。方位向天线长度、俯仰向天线高度最大方向性系数。对于偏离波束中心的角度方向图因子为收发增益假设 Ch2 收发共用。单基雷达方程接收功率其中为峰值功率为系统损耗包括馈线、大气衰减等。复基带信号幅度取以使信号平均功率等于。2.4 运动目标回波若目标具有恒定方位向速度和地距向速度正值为远离雷达在慢时间其坐标更新为其斜距历程不再满足静止目标的简单双曲线需逐脉冲代入精确斜距公式不采用 EPC 近似以准确保留运动引起的相位变化。2.5 接收机噪声噪声建模为加性复高斯白噪声每个通道的噪声功率由下式决定其中J/K 为玻尔兹曼常数K 为标准温度为噪声带宽取采样频率为噪声系数线性值。噪声独立添加至各通道。3. 距离压缩匹配滤波距离向进行脉冲压缩以提高距离分辨率。匹配滤波器的冲激响应为发射 LFM 信号的共轭反转在频域实现匹配滤波更高效对每个脉冲的回波快时间信号作 FFT乘以传递函数再经 IFFT 得到压缩输出。为压低旁瓣通常对加窗如汉明窗。物理效果压缩后点目标回波被压缩成宽度约为的 sinc 状脉冲峰值位于延迟处距离分辨率。三通道独立进行距离压缩得到。4. 方位向傅里叶变换与距离多普勒RD域对距离压缩后的信号沿方位向慢时间进行 FFT变换到距离多普勒RD域为多普勒频率轴范围。在 RD 域静止目标的方位信号近似为线性调频LFM信号其频谱在方位带宽内较平坦运动目标因径向速度会产生多普勒中心偏移使得其能量在多普勒谱上偏离静止杂波的主瓣。这为后续通过频率选择性如多普勒凹口和对消处理提供了基础。定义与径向速度的关系正侧视条件下5. RD 域通道配准由于等效相位中心存在沿航迹偏移需将 Ch1 和 Ch3 的信号在方位时间上平移使三通道的静止目标回波对齐。时间平移对应于频域线性相位。由 EPC 近似Ch1 的斜距历程相当于 Ch2 滞后秒。故。傅里叶变换的时移性质给出因此要将 Ch1 与 Ch2 对齐需乘以相反的相位因子同理Ch3 的斜距相当于 Ch2 超前所以配准相位为实现时对所有距离门乘以同一个相位因子因为基线偏移与距离无关至此理想情况下静止杂波在三个通道的 RD 谱中将完全重合。6. DPCA 杂波对消配准后静止目标在三个通道的 RD 谱相同因此可通过相减消除。对消方式为综合对消结果取两者之和增大动目标残差物理意义满足 DPCA 条件且通道完全一致时静止杂波被完全对消残差仅由噪声、通道误差和系统非理想因素构成运动目标由于自身运动破坏了时移关系对消后仍有显著残留从而突显目标。杂波对消比 CRR 定义为对消前 Ch2 平均功率与对消后残差中值功率之比采用中值是为了避免运动目标能量拉高残差功率的均值。CRR 是衡量 DPCA 性能的关键指标通常要求 10 dB。7. ATI 沿航迹干涉测速利用配准后的 Ch3 与 Ch1 进行干涉测量运动目标的径向速度。干涉图干涉相位7.1 干涉相位与径向速度的关系考虑一个径向速度为的运动目标。其多普勒频率为正侧视。对于这样的目标Ch1 和 Ch3 的原始回波之间存在时间差因为等效相位中心相距对应时间差。频域表现为线性相位差。配准过程是对齐静止目标的运动目标的配准相位与静止目标相同因此配准后 RD3r 和 RD1r 仍保留有因而产生的相对相位差所以代入得到于是径向速度估计式为无模糊速度干涉相位只能在之间无模糊测量。当时对应最大可测速度对于仿真参数, 得。7.2 相位缠绕处理实际速度超出无模糊范围时干涉相位会发生缠绕。本仿真暂未进行解缠但无模糊速度范围足以覆盖典型车辆目标的速度。8. 二维 CFAR 检测与多普勒凹口8.1 单元平均 CFAR 检测器在 DPCA 对消后的幅度图或功率图上进行二维恒虚警CFAR检测。对每个被检单元设置保护窗半宽度和训练窗半宽度构成一个空心矩形窗。用窗内但不包括保护区的单元估算背景杂波功率其中为训练单元总数。检测门限为门限因子与设计虚警概率和训练单元数的关系为若被检单元功率大于则判决为目标。8.2 多普勒凹口DC Notch即使经 DPCA 对消零多普勒DC附近的杂波残余仍然很强容易产生虚警。为此在检测前设置一个凹口排除多普勒中心附近的若干单元代码中取 ±4 个 bin。物理上这些区域对应径向速度极低的目标或静止地物残余不是 GMTI 关注的对象。9. 距离多普勒算法RDA背景成像为提供视觉背景对配准后的距离压缩数据通常取 Ch2进行完整的 RDA 成像处理得到聚焦的 SAR 图像。主要步骤9.1 距离徙动校正RCMC目标在合成孔径时间内由于平台运动其斜距随时间变化导致回波包络在快时间上发生弯曲即距离徙动。徙动量与多普勒频率有关其中为徙动因子。校正时在 RD 域进行插值如 sinc 插值将每个距离门的数据按平移。9.2 方位压缩校正后方位向信号变为理想的线性调频信号其匹配滤波器为乘以该相位后作 IFFT即完成方位向压缩得到二维 SAR 图像。图像中静止目标聚焦良好运动目标则因多普勒参数失配可能散焦或位移但其信息并不直接用于检测仅作背景。10. 性能评估指标代码最后计算并输出了多项性能指标用以评估 GMTI 系统的能力。10.1 分辨率- 距离分辨率(m)。- 方位分辨率(m)其中为方位天线长度。10.2 杂波对消比 CRR如前所述CRR 反映 DPCA 对静止杂波的抑制效果。10.3 信杂噪比改善 SCNR其中为检测到的目标在 DPCA 图中的平均功率。10.4 最小可检测速度 MDVMDV 粗略估计为无模糊速度与 SCNR 开方的比值它表征系统在给定 SCNR 下能可靠检测的最小径向速度。10.5 目标检测统计统计检测点个数、估计速度的均值、标准差并与真实速度比较得到测速误差。