步进频 ISAR 成像仿真项目
项目介绍
本项目致力于利用 MATLAB 平台实现基于步进频信号(Stepped Frequency Waveform, SFW)的逆合成孔径雷达(ISAR)成像仿真。系统首先通过建立空间旋转目标的散射中心模型,设置目标在笛卡尔坐标系中的几何分布及散射强度。核心成像流程采用步进频信号波形设计,通过逐个脉冲改变载频来实现合成大带宽,从而在硬件处理瞬时带宽有限的情况下获得极高的距离向分辨率。本项目实现了从回波生成到二维成像的完整流程,是研究雷达信号处理和目标探测特性的基础实验环境。
功能特性
- 步进频信号体制:通过合成带宽技术克服瞬时带宽限制,实现高分辨率距离向探测。
- 多点散射建模:预设典型的飞机形状散射模型,支持自定义目标的几何分布。
- 动态回波仿真:模拟目标在旋转过程中的相干回波,包含相位延迟计算及复高斯白噪声注入。
- 二维图像重构:采用距离-多普勒(Range-Doppler)算法,结合加窗抑制旁瓣技术,生成清晰的目标图像。
- 定量指标评估:系统能够自动估算成像后的目标物理尺寸,并对比理论分辨率。
使用方法
- 启动 MATLAB 软件(建议 2018b 及以上版本)。
- 将该仿真代码复制并保存为 .m 文件。
- 直接点击“运行”按钮或在命令行窗口输入函数名称。
- 程序将自动执行仿真并在窗口中显示四个维度的对比图(建模图、一维距离像、距离走动图、二维重构图)。
- 运行结束后,可在 MATLAB 控制台查看理论分辨率及目标尺寸估计报告。
系统要求
- 软件要求:MATLAB 基础版(无需特殊工具箱,主要依赖矩阵运算与 FFT 函数)。
- 硬件保障:标准个人电脑即可,由于采用 128x128 矩阵运算,内存占用极低。
实现逻辑与算法细节说明
1. 系统参数初始化
系统设定起始频率为 10GHz,频率步进为 10MHz,通过 128 个子脉冲合成 1.28GHz 的总带宽。仿真设置了目标的旋转角速度以及脉冲串的重复周期,为后续的方位向多普勒处理提供时间基准。
2. 散射模型与动态回波生成
- 几何建模:代码定义了一个包含 9 个散射点的飞机模型,每个点拥有独立的 [x, y] 坐标和散射强度。
- 旋转模拟:在慢时间维度(M 个脉冲串),程序利用旋转矩阵更新目标在每个采样时刻的瞬时空间位置。
- 回波计算:基于远场模型,通过计算各散射点到雷达的射影距离,叠加各个频率点上的相位响应。其核心数学表达为:
exp(-j * 4 * pi * fk * range_p / c)。 - 噪声处理:通过功率计算在回波中加入指定的复高斯白噪声,以增强仿真的真实性。
3. 一维脉冲压缩(距离向处理)
针对步进频信号,系统对每个脉冲串的频率采样点执行逆快速傅里叶变换(IFFT)。在变换前,施加了 Hamming 窗函数以降低距离向的旁瓣电平。这一步骤成功将频域回波转化为时域的一维高分辨距离像(HRRP),反映了目标沿雷达视线方向的散射分布。
4. 二维成像处理(方位向处理)
- 相干累积:系统对处理后的距离像矩阵进行方位向(慢时间轴)的 FFT 操作。
- 多普勒映射:由于目标旋转引起的各散射点多普勒频率差异,通过 FFT 可以将目标在方位向的分散能量进行聚焦。
- 中心偏移:使用
fftshift 函数将零频分量移至频谱中心,确保图像符合物理几何视角。
5. 坐标变换与结果可视化
- 物理坐标计算:根据光速、合成带宽、波长和旋转角速度,将 IFFT 和 FFT 后的采样点索引转化为真实的物理距离(米)。
- 多维度展示:代码利用
subplot 功能同步展示散射中心模型、单次采样的一维距离像、展现距离走动特性的 Sinogram 图,以及最终生成的二维 ISAR 彩色重构图。
6. 尺寸估算逻辑
程序通过图像分割算法,预设灵敏度阈值来提取二维图像中的目标轮廓。通过查找二值化图像中目标占据的像素行列范围,乘以对应的距离向和方位向分辨率,得出目标的估计长度和宽度,实现了从成像到特征提取的完整闭环。