MIMO雷达系统全流程仿真工具包 (MIMO Radar System Simulation Toolkit)
项目简介
本项目是一个基于 MATLAB 开发的毫米波 MIMO(多输入多输出)雷达系统仿真平台。该工具包旨在演示 MIMO 技术相对于传统相控阵雷达在角度分辨率和阵列孔径方面的优势。项目完整模拟了从发射波形生成、信道传播、接收回波模拟,到信号处理(脉冲压缩、虚拟阵列合成)及高级参数估计(RD 图、DBF、Capon、MUSIC)的全流程。
该仿真主要针对 77GHz 车载毫米波雷达场景,采用 TDM(时分多址)MIMO 体制。
功能特性
- TDM-MIMO 发射机制模拟:模拟 4 发 4 收(4Tx, 4Rx)的天线配置,通过时分复用方式实现正交性,降低了有效脉冲重复频率(PRF)。
- 虚拟孔径合成:利用稀疏发射阵列和接收阵列的卷积效应,将 4x4 的物理通道合成为 16 个等效虚拟阵元的均匀线阵(ULA),显著提高了角度分辨率。
- 高保真回波信号生成:支持多目标环境配置,能够模拟目标的距离、角度、速度(多普勒效应)和雷达散射截面(RCS),并叠加高斯白噪声。
- 基础信号处理:实现了基于 FFT 的匹配滤波(脉冲压缩)和维度重塑,构建“快时间-虚拟阵元-慢时间”的数据立方体。
- 距离-多普勒分析:生成 RD 图(Range-Doppler Map),用于同时估计目标的距离和相对速度。
- 高级 DOA 估计算法:集成了三种不同的波达方向估计算法(DBF、MIMO-Capon、MIMO-MUSIC),用于对比不同算法在空间谱估计上的性能差异。
系统要求
- MATLAB R2018a 或更高版本
- Signal Processing Toolbox(建议安装,用于部分信号处理函数)
使用方法
- 确保 MATLAB 工作路径已包含本项目的脚本文件。
- 直接运行主函数
main。 - 程序将自动执行参数初始化、回波生成、信号处理算法,并弹出包含 5 个子图的综合结果窗口。
代码实现逻辑与算法详解
本项目采用单文件 main.m 实现全流程,具体实现逻辑如下:
1. 系统与阵列参数配置
- 物理参数:设定载频为 77GHz,带宽 150MHz,脉宽 10us。
- MIMO 阵列设计:
* 接收阵列(Rx):4 阵元,间距为半波长(lambda/2)。
* 发射阵列(Tx):4 阵元,采用稀疏布置,间距为 Rx 阵列的总长度(4 * lambda/2),以此最大化虚拟孔径。
- 仿真控制:设定信噪比(SNR)为 10dB,慢时间维包含 64 个脉搏用于多普勒处理。
2. 目标环境构建
代码中预设了 4 个典型目标,用于测试系统的各项分辨能力:
- 目标1 & 2:距离相近(80m 和 85m),角度相同,测试距离分辨率。
- 目标3:远距离、负速度目标。
- 目标4:近距离静止目标。
3. TDM-MIMO 回波生成 (核心逻辑)
代码未直接调用工具箱函数,而是通过原理公式手动生成回波数据
raw_data:
- 信号模型:采用 LFM(线性调频)基带信号模型。
- 相位计算:精确计算了由几何位置(Tx 和 Rx 坐标)引起的空间相位差,以及由目标运动引起的多普勒相位差。
- TDM 机制:在生成回波时,慢时间变量
slow_time 随发射天线序号 tx 偏移,模拟了不同天线在不同时间片发射的过程。 - 噪声叠加:根据设定的 SNR 计算信号功率,并叠加复高斯白噪声。
4. 接收端信号处理
- 脉冲压缩 (Matched Filtering):对原始回波数据在快时间维进行 FFT,与发射波形的频谱共轭相乘,再进行 IFFT,以获得高信噪比的距离像。
- 虚拟阵列合成:将数据维度从
[采样点, Rx, Tx, 脉冲] 重塑为 [采样点, Rx*Tx, 脉冲]。由此构建出 16 个虚拟阵元的等效数据,且阵元间距满足半波长条件。
5. 结果可视化与高级算法
代码最终绘制了一个包含 5 个子图的 Figure,展示处理结果:
- 虚拟阵列示意图:直观展示了 Tx、Rx 物理位置以及合成后的虚拟阵元位置,验证了虚拟孔径的扩展效果。
- 一维距离像:取第一个虚拟通道的数据展示脉冲压缩结果,可清晰观察到目标的距离峰值。
- 距离-多普勒 (RD) 图:
* 在慢时间维加 Hanning 窗以抑制旁瓣。
* 进行 2D-FFT 处理。
* 根据 TDM 体制下的有效 PRF 计算速度轴刻度。
* 针对特定目标距离门的数据,计算协方差矩阵。
* 生成标准导向矢量,进行加权求和($w^H R w$),绘制空间谱。
* 计算协方差矩阵的逆矩阵(引入对角加载
1e-6 * eye 防止矩阵奇异)。
* 基于 MVDR 准则(最小方差无畸变响应)计算空间谱,相比 DBF 具有更好的分辨率和旁瓣抑制能力。
* 对协方差矩阵进行特征分解。
* 区分信号子空间和噪声子空间。
* 利用信号子空间与噪声子空间的正交性,构建极窄的谱峰以实现超高分辨率的角度估计。