毫米波FMCW雷达目标探测与全链路仿真系统

资 源 简 介

本项目旨在构建一个高精度的毫米波FMCW（调频连续波）雷达全链路仿真平台，用于模拟和验证雷达在复杂环境下的探测性能。系统首先实现发射端信号建模，生成具有可配置带宽、斜率和时宽的线性调频脉冲（Chirp）信号，并支持MIMO（多输入多输出）天线阵列的发射模式配置。在信道传播模拟阶段，系统能够计算多目标的电磁波传播损耗、由于距离产生的时延以及目标运动产生的多普勒频移，同时叠加高斯白噪声和瑞利/对数正态分布的杂波干扰，以模拟真实的物理环境。在接收信号处理环节，项目包含混频解调、低通滤波及离散化采样模块，生成原始中频（IF）数据。核心算法模块依次执行距离维FFT（快时间处理）和速度维FFT（慢时间处理），生成距离-多普勒热力图（Range-Doppler Map）。随后，集成单元平均CFAR（CA-CFAR）或有序统计CFAR（OS-CFAR）算法，在动态噪声背景下精准提取目标点。最后，利用数字波束形成（DBF）或超分辨MUSIC算法进行到达角（DOA）估计，从而输出所有探测目标的精确距离、径向速度和方位角信息，为雷达波形优化和感知算法开发提供数据支撑。

详 情 说 明

毫米波雷达目标探测与信号处理仿真系统

项目介绍

本项目构建了一个高精度的毫米波FMCW（调频连续波）雷达全链路仿真平台。该系统基于MATLAB开发，能够模拟从回波信号生成、信道传播、接收信号处理到最终目标检测与角度估计的完整雷达信号处理流程。

仿真系统工作在77GHz汽车雷达频段，采用标准的线性调频脉冲（Chirp）波形，支持多接收天线（MIMO/SIMO）配置。通过该平台，用户可以深入理解FMCW雷达的测距、测速和测角原理，并验证CFAR检测与DBF波束形成算法在复杂多目标场景下的性能。

功能特性

• 雷达参数化建模：支持自定义载频、带宽、Chirp周期、采样率、天线数量及间距等核心参数。
• 复杂目标场景模拟：支持设置多个具有不同距离、速度、角度和雷达散射截面（RCS）的目标。
• 中频信号仿真：直接生成经过混频后的原始中频（IF）数据，包含精确的距离时延相移、多普勒频移及阵列相位差。
• 环境噪声模拟：在原始数据中叠加指定信噪比（SNR）的高斯白噪声。
• 全流程信号处理：

* 2D-FFT处理：包含距离维FFT（加窗）和多普勒维FFT（加窗及移频）。 * 非相参积累：对多天线数据进行功率叠加，提高信噪比。

• 恒虚警检测（CFAR）：实现了二维单元平均CFAR（2D CA-CFAR）算法，自适应提取目标点。
• 角度估计（DOA）：基于数字波束形成（DBF）技术，对检测到的目标进行高精度方位角估计。
• 多维可视化：提供时域波形、距离谱、距离-多普勒（RD）热力图及检测结果点云图。

系统要求

• MATLAB R2018a 或更高版本
• Signal Processing Toolbox（信号处理工具箱）
• Phased Array System Toolbox（相控阵系统工具箱，可选，本项目核心算法为手写实现，不强依赖此工具箱）

核心算法与实现逻辑

本项目代码主要分为八个处理阶段，具体逻辑如下：

1. 系统参数配置

初始化雷达的物理层参数。设定载波频率为77GHz，带宽为150MHz，Chirp周期60us。根据这些参数系统自动计算调频斜率、距离分辨率、速度分辨率以及最大探测范围。天线阵列被配置为4阵元的均匀线阵（ULA），间距为半波长。

2. 目标场景定义

定义仿真场景中的目标列表。每个目标包含四个属性：距离（m）、速度（m/s，远离为正）、角度（deg）和RCS。通过矩阵形式配置多个目标，用于验证多目标分辨能力。

3. 回波信号生成（IF建模）

这是仿真的物理核心。代码不分别生成发射和接收射频信号，而是根据FMCW原理直接合成去斜后的中频（IF）信号。

• 多维数据生成：创建 [采样点数 x Chirp数 x 接收天线数] 的三维数据立方体。
• 相位模型：对每个采样点，计算由于距离产生的差拍频率（$f_r$）和由于速度产生的多普勒频率（$f_d$）。
• 阵列流形：根据由于入射角导致的天线波程差，计算每个接收天线的附加相位偏移。
• 动态更新：在慢时间维度（Chirp间）更新目标距离，模拟目标运动对回波的影响。
• 幅度模型：基于雷达方程简化模型，幅度与RCS成正比，与距离的平方成反比。

4. 噪声注入

使用MATLAB内置函数生成高斯白噪声，并按照设定的SNR（如15dB）叠加到生成的中频信号上，模拟真实链路中的热噪声。

5. 距离-多普勒处理（2D-FFT）

对原始数据立方体进行频域变换：

• 距离维FFT：在快时间维度（采样点方向）加汉明窗后进行FFT，提取距离信息。
• 多普勒维FFT：在慢时间维度（Chirp方向）对距离FFT的结果加窗后进行FFT，并使用 fftshift 将零频移至中心，提取速度信息。
• 非相参积累：将4个接收天线的RD图模值平方后求和，生成最终的功率谱（RD Map），并转换为对数（dB）刻度。

6. 目标检测（CFAR）

在RD图上执行二维单元平均恒虚警检测（2D CA-CFAR）：

• 滑动窗口：遍历RD图的有效区域。
• 参考单元与保护单元：定义训练窗口（Tr, Td）用于估计噪声底，定义保护窗口（Gr, Gd）防止目标能量泄露到训练单元。
• 自适应阈值：根据虚警概率（$P_{fa}$）计算门限因子，结合局部噪声估值动态设定检测阈值。
• 点提取：记录所有超过阈值的单元坐标（距离索引，多普勒索引）。

7. 到达角估计（DOA）

对CFAR检测到的目标点进行角度解算：

• 数据提取：针对每个检测到的 (Range, Doppler) 单元，提取原始数据立方体中对应的多天线阵列向量（1x4复数向量）。
• 数字波束形成（DBF）：构建标准导向矢量，在-90度到+90度范围内进行空间扫描。
• 峰值搜索：计算空间谱，寻找谱峰对应的角度作为目标的方位角估计值。
• 信息聚合：将解算出的距离、速度、角度和功率信息整合成最终的目标列表。

8. 结果可视化

利用MATLAB绘图功能展示仿真全过程的关键数据：

• 时域图：展示第一个Chirp在第一个天线上的中频信号实部波形。
• 距离谱：展示经过一维FFT后的频谱，呈现目标的距离分布。
• RD热力图：通过颜色映射展示距离-多普勒平面的能量分布，直观反映目标的速度和距离关系。
• 检测结果：在图中标记出CFAR检测到的目标点及其解算出的角度信息（注：代码逻辑包含此部分，用于后续的点云展示）。

使用方法

1. 打开MATLAB软件。
2. 将包含上述逻辑的脚本文件设置为当前工作目录或添加到路径中。
3. 直接运行主函数。
4. 程序将在控制台输出雷达设计指标及检测到的目标数量，并弹出一个综合图形窗口显示仿真结果。