单基三通道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 最小可检测速度 MDV
MDV 粗略估计为无模糊速度与 SCNR 开方的比值:
它表征系统在给定 SCNR 下能可靠检测的最小径向速度。
10.5 目标检测统计
统计检测点个数、估计速度的均值、标准差,并与真实速度比较,得到测速误差。
