相位编码脉内线性调频混合调制雷达系统仿真
项目简介
本项目实现了一个基于MATLAB的深度雷达信号处理仿真平台。核心重点在于构建和评估一种混合调制波形——结合了相位编码(13位巴克码)与脉内线性调频(LFM)的信号体制。这种混合波形旨在兼顾相位编码的低截获概率特性与LFM的大时宽带宽积优势。
仿真系统完整覆盖了从发射波形设计、复杂的环境回波模拟、接收端数字化处理到最终目标检测(CFAR)的全链路流程,能够直观展示复杂波形在现代雷达系统中的应用效果。
功能特性
- 混合调制波形设计:实现了13位巴克码对LFM信号的相位调制,生成大带宽、复杂结构的发射信号。
- 高保真回波模拟:基于物理模型模拟回波,包含时延、多普勒频移(针对每个脉冲精确计算相位)以及距离衰减。
- 数字化接收机前端:模拟中频(IF)采样,通过数字下变频(DDC)和低通滤波获取基带I/Q信号。
- 全链路信号处理:包含脉冲压缩、MTI杂波抑制、MTD相干积累和CFAR恒虚警检测。
系统要求
- MATLAB R2016a 或更高版本
- Signal Processing Toolbox(用于滤波器设计和FFT处理)
使用方法
- 确保MATLAB路径中包含本项目脚本。
- 直接运行主函数(
main),脚本将自动执行以下步骤:
* 初始化系统参数与环境配置。
* 生成混合调制发射信号。
* 逐脉冲模拟回波并进行DDC处理。
* 依次执行脉冲压缩、MTI、MTD和CFAR处理。
* 最后弹出图形窗口展示仿真结果(部分可视化)。
详细功能与实现逻辑分析
本项目代码逻辑清晰划分为九个主要模块,具体实现细节如下:
1. 系统参数设置
构建了典型脉冲多普勒雷达的参数体系:
- 射频与时序:设定了30MHz带宽、20us脉宽、50MHz中频载波以及2000Hz的脉冲重复频率(PRF)。
- 混合调制参数:定义了13位巴克码(Barker Code),计算码元宽度,将脉宽分割为对应的时间片。
- 时域采样:基于奈奎斯特采样定理设定采样率(4倍中频),并构建了快时间(距离维)和慢时间(脉冲维)坐标轴。
2. 目标与环境配置
设置了三个具有代表性的点目标用于验证算法性能:
- 目标1:高速接近的运动目标(15km处,150m/s)。
- 目标2:远离的弱小运动目标(20km处,-80m/s,RCS较小)。
- 目标3:强静止杂波/地物(10km处,0m/s,高RCS),用于测试MTI的对消性能。
同时定义了信噪比(SNR),并换算为高斯白噪声的功率水平。
3. 发射信号生成 (混合调制)
实现了LFM与相位编码的深度融合:
- 首先生成基础的线性调频信号(LFM chirp)。
- 根据时间轴映射,将13位巴克码的相位跳变(0或180度)强加于LFM信号之上。
- 上变频:将生成的复基带混合信号调制到50MHz中频上,并取实部以模拟真实的发射射频/中频信号。
4. 回波模拟与数字化正交解调 (DDC)
这是仿真的核心物理层模型:
- 回波生成:在每个PRI(脉冲重复周期)内,根据目标的距离和速度计算瞬时时延和多普勒相移。代码通过构建复基带回波并叠加多普勒项,再上变频的方式,精确保留了相位信息。
- 噪声注入:在接收信号中叠加复高斯白噪声。
- DDC实现:
1.
混频:接收到的中频实信号与数字本振(NCO)相乘。
2.
低通滤波:使用6阶巴特沃斯低通滤波器滤除倍频分量,提取基带I/Q正交信号。
5. 脉冲压缩 (Matched Filtering)
- 采用频域匹配滤波算法。
- 将接收到的回波数据与发射基带信号的共轭频谱在频域相乘,然后进行IFFT。
- 此步骤显著提高了距离分辨率和信噪比,实现了大时宽带宽积的压缩增益。
6. 动目标显示 (MTI)
- 设计了一个三脉冲对消器(系数为
[1, -2, 1])。 - 该滤波器本质上是一个沿慢时间维度的高通FIR滤波器,旨在滤除零多普勒频率的分量(即滤除代码中设置的目标3——静止强杂波)。
- 代码处理了滤波器带来的瞬态效应,截去了前两个无效脉冲。
7. 动目标检测 (MTD)
- 对MTI处理后的数据进行慢时间维度的加窗(Hamming窗)FFT处理。
- 通过相干积累,将不同速度的目标聚焦在不同的多普勒通道上,进一步提升了运动目标的信杂比(SCR)。
8. 恒虚警检测 (CFAR)
- 实现了单元平均恒虚警检测器 (CA-CFAR)。
- 算法逻辑:
1. 选取多普勒维度的最大能量切片作为检测剖面。
2. 定义保护单元和参考单元,滑动窗口计算背景杂波/噪声的平均功率。
3. 根据预设的虚警概率(Pfa=1e-4)计算自适应门限。
4. 对比信号幅度与门限,输出检测结果(目标距离索引)。
9. 结果可视化
代码包含绘图模块,用于通过图表展示系统的关键中间状态,包括:
- 发射混合波形的时域图与频谱图。
- (注:代码可视化部分在发射波形展示后截断,但逻辑上支持完整的处理结果展示)。