基于MATLAB的SAR成像RD与CS算法全流程仿真系统

资 源 简 介

本项目旨在通过MATLAB平台，深入实现并对比分析合成孔径雷达（SAR）领域最经典的两种成像算法：距离多普勒（RD）算法和Chirp Scaling（CS）算法。项目首先建立雷达回波仿真模型，能够根据用户自定义的雷达系统参数（如载波频率、信号带宽、脉冲宽度、脉冲重复频率、平台速度、飞行高度等）以及目标场景几何分布，生成多点目标的原始基带回波数据。在成像处理环节，RD算法模块完整实现了距离向脉冲压缩、距离徙动曲线计算、基于Sinc插值的距离徙动校正（RCMC）、方位向脉冲压缩以及相位残差补偿；CS算法模块则利用变标方程，通过随距离变化的频率调制，在二维频域内通过纯相位操作完成距离徙动校正，有效避免了插值带来的计算负担和相位误差。此外，项目还集成了成像质量评估功能，能够自动提取点目标的冲激响应函数（IRF），计算并输出空间分辨率、峰值旁瓣比（PSLR）和积分旁瓣比（ISLR）等关键指标，帮助用户直观理解算法物理机制及性能差异，适用于雷达信号处理教学、算法研究及工程原理验证。

详 情 说 明

基于MATLAB的SAR成像算法(RD与CS)全流程仿真系统

项目简介

本项目是一个基于MATLAB平台开发的合成孔径雷达（SAR）成像仿真系统。项目旨在对比分析SAR信号处理领域中最经典的 距离多普勒（RD）算法 与 Chirp Scaling（CS）算法。

系统实现了从雷达系统参数配置、原始回波信号生成、成像算法处理到最终图像质量评估的全流程闭环。通过该项目，用户可以直观地理解点目标在SAR成像过程中的二维压缩机制，观测距离徙动校正（RCMC）的效果，并定量评估不同算法的聚焦性能。

功能特性

• 全参数化系统设计：支持自定义载波频率、信号带宽、脉冲宽度、PRF、平台高度与速度等关键雷达参数（默认配置为C波段机载SAR）。
• 原始回波仿真：基于“Stop-and-Go”模型，生成包含多个点目标（中心、近端、远端）的基带回波数据，模拟真实的相位历程。
• RD算法实现：包含完整的距离压缩、距离徙动校正（基于插值）、方位压缩流程，并应用了窗函数以抑制旁瓣。
• CS算法实现：利用Chirp Scaling原理，通过纯相位操作在二维频域完成距离徙动校正，避免了时域插值操作。
• 可视化与评估：提供原始回波的实部与幅度显示，对比RD与CS算法的最终成像结果，并包含针对点目标的质量评估模块（分辨率、PSLR、ISLR）。

系统与环境要求

• MATLAB R2016b 或更高版本
• Signal Processing Toolbox（用于FFT、窗函数等操作）

快速开始

1. 下载本项目的所有文件。
2. 在MATLAB中打开 main.m 文件。
3. 直接运行 main.m。
4. 程序将自动依次执行：参数初始化 -> 回波生成 -> RD成像 -> CS成像 -> 结果绘图与指标输出。

---

核心算法与代码实现逻辑

本节详细描述 main.m 中实际实现的功能逻辑。

1. 系统参数配置与场景构建

代码首先定义了名为 P 的结构体，用于封装所有系统参数。

• 雷达参数：设定了5.3GHz载频、30MHz带宽、2.5us脉宽等典型机载SAR参数。
• 几何参数：计算了场景中心斜距、合成孔径长度及对应的时间窗。
• 采样空间：根据奈奎斯特采样定理及过采样率，自动计算距离向（Nr）和方位向（Na）的离散采样点数，并建立了快时间（tr）和慢时间（ta）轴。

2. 原始回波生成 (generate_echo)

该模块负责生成雷达的原始基带信号：

• 多目标支持：定义了由中心点、左近点、右远点组成的点目标阵列。
• 信号模型：采用时域逐点累加方式。对于每一个目标，计算其随慢时间变化的瞬时斜距 $R(eta)$。
• 相位历程：精确模拟了由距离变化引起的双程相位延迟 $-4pi f_c R/c$ 以及线性调频信号的二次相位 $pi K_r (t - 2R/c)^2$。
• 包络控制：加入了距离向的矩形窗函数，方位向则通过仿真时间长度隐式截断。

3. 距离多普勒算法 (run_rd_algorithm)

实现了标准的RD成像流程，核心步骤如下：

• 距离压缩：将数据变换到距离频域，生成匹配滤波器（包含Hamming窗），相乘后变换回距离时域。此步骤完成后，目标在距离向聚焦，但存在距离徙动。
• 变换域：对数据进行方位向FFT，进入“距离时域-方位频域”（Range-Doppler Domain）。
• 距离徙动校正 (RCMC)：

* 计算徙动因子 $D(f_a)$。 * 插值实现：代码利用 interp1 函数，逐行（即针对每个方位频率）对距离向数据进行线性插值，将弯曲的轨迹校正为直线。

• 方位压缩：

* 在RD域计算方位向匹配滤波器。 * 构造包含 $D(f_a)$ 的精确相位补偿项。 * 应用Hamming窗进行加权，最后通过IFFT完成方位聚焦。

4. Chirp Scaling算法 (run_cs_algorithm)

实现了基于变标原理的CS成像流程，核心步骤如下：

• 变换域：首先进行方位向FFT，进入RD域。
• Chirp Scaling操作 (H1)：计算Scaling因子，构造二次相位函数，对信号进行相位调制。这一步通过调整调频斜率，使不同距离门的徙动曲线趋于一致。
• 距离压缩与一致RCMC (H2)：

* 变换到二维频域。 * 构造综合相位滤波器 $H_2$，同时完成距离压缩、二次距离压缩（SRC）以及一致（Bulk）距离徙动校正。

• 方位压缩 (H3/Azimuth MF)：

* 将数据变换回RD域。 * 由于代码实现了简化的CS流程，此处采用逐个距离门处理的方式。 * 针对每一个距离门，计算对应的方位匹配滤波器（去调频），并应用相位校正，最后通过方位IFFT完成成像。

5. 质量评估 (analyze_quality)

*注：根据代码逻辑，该模块负责对成像后的点目标进行指标分析。*

• 峰值提取：自动搜索图像中的最大值位置（即点目标中心）。
• 切片分析：分别提取通过峰值的距离向和方位向剖面。
• 指标计算：计算并输出空间分辨率（主瓣宽度）、峰值旁瓣比（PSLR）和积分旁瓣比（ISLR），用于量化评估算法生成的图像质量。

结果展示

运行程序后，系统将弹出两个主要窗口：

1. 原始回波数据：展示回波实部和幅度的二维分布，可见明显的距离徙动曲线。
2. 成像结果对比：左右分栏显示RD算法和CS算法处理后的最终聚焦图像。

控制台（Command Window）将实时输出各阶段的执行状态、算法耗时以及针对中心目标的各项成像质量指标。