雷达有源噪声干扰效能仿真系统
项目简介
本项目构建了一个高保真的雷达对抗仿真环境,基于MATLAB平台开发。系统专注于模拟线性调频(LFM)脉冲多普勒雷达在面临不同类型有源噪声干扰时的探测性能。通过从信号产生、环境传输、干扰叠加到接收机信号处理的全链路仿真,系统能够定性和定量地评估射频噪声、调幅噪声、调频噪声对雷达脉冲压缩和CFAR检测的具体影响。
主要功能特性
- 高逼真信号建模:支持线性调频(LFM)信号发射与具备多普勒频移和时延的点目标回波生成。
- 多样式干扰生成:内置三种经典压制性干扰模型(射频噪声、调幅噪声、调频噪声),支持自定义干信比(JSR)和干扰带宽。
- 全链路信号处理:实现了数字接收机的核心算法,包括匹配滤波(脉冲压缩)和恒虚警检测(CFAR)。
- 量化效能评估:提供无干扰与有干扰环境下的对比分析,并具备自动化扫描干信比(JSR)以分析信噪比损失(PCSNR)的功能。
- 多维可视化展示:集成时域波形、频域谱图、距离维脉压输出以及JSR趋势曲线的综合显示。
系统要求
- MATLAB R2016b 或更高版本
- Signal Processing Toolbox(信号处理工具箱)
使用方法
直接运行
main.m 脚本即可启动仿真。程序将自动执行参数初始化、信号生成、处理与绘图,最终输出包含6个子图的综合分析窗口。
用户可以通过修改代码开头的“系统参数设置”区域来调整雷达载频、带宽、目标距离、干扰类型及干信比等关键参数。
---
代码实现逻辑详解
main.m 脚本包含了完整的仿真控制流,其内部处理逻辑严格遵循以下步骤:
1. 系统参数初始化
定义了光速、雷达参数(载频10GHz、带宽30MHz、脉宽10us)、目标参数(距离5km、速度100m/s、RCS 1平米)以及干扰机参数。特别地,干扰类型通过变量
Jamming_Type 切换(1=RF, 2=AM, 3=FM),干信比由
JSR_dB 控制。
2. LFM信号与目标回波建模
- 发射信号:构建复包络LFM信号
S_ref,根据带宽和脉宽计算调频斜率 $K$。 - 目标回波:计算目标往返时延
tau 和多普勒频移 fd。在接收序列中对应的位置插入衰减后的参考信号,并叠加多普勒相位项 exp(j*2*pi*fd*t),真实模拟动目标回波特性。
3. 有源噪声干扰生成
调用子函数
generate_jamming 生成干扰信号 $J(t)$。该过程首先计算干扰所需的功率水平,并结合设定的干扰带宽和中心频率偏差生成干扰波形。
4. 环境混合
仿真真实电磁环境,生成接收信号
Rx_signal,其组成为:目标回波 + 干扰信号 + 热噪声(根据设定SNR生成的高斯白噪声)。同时生成一个不含干扰的参考信号
Rx_no_jamming 用于对比。
5. 脉冲压缩处理(匹配滤波)
在频域执行脉冲压缩。
- 对参考信号添加 Hamming 窗以抑制副瓣。
- 利用FFT/IFFT实现快速卷积:$PC = IFFT(FFT(Rx) times conj(FFT(S_{ref})))$。
- 结果进行归一化并转换为对数刻度(dB)以便于观察。
6. CFAR 恒虚警检测
调用
cf_ar_algo 子函数执行单元平均恒虚警检测(CA-CFAR)。算法使用滑动窗口机制,根据参考单元的平均功率和预设的虚警概率(Pfa=1e-4)动态计算检测门限。
7. JSR 参数扫描分析
代码包含一段自动化循环,在 -10dB 到 40dB 的范围内扫描干信比(JSR)。针对每一个 JSR 值,重新生成干扰并在脉压后计算“峰值平均底噪比”(PCSNR)的衰减情况,从而得到干扰压制曲线。
8. 可视化分析
代码最后生成一个包含六个子图的窗口:
- 无干扰时域图:展示微弱的目标回波淹没在噪声中。
- 有干扰时域图:展示大功率噪声干扰对时域信号的覆盖。
- 频谱对比:对比目标回波与干扰信号的频谱重叠情况。
- 无干扰脉压输出:清晰的目标尖峰与自适应CFAR门限。
- 有干扰脉压输出:目标尖峰被干扰抬升的底噪淹没或CFAR门限被抬高的情况。
- JSR影响曲线:定量展示随着干扰增强,雷达输出信噪比的下降趋势。
---
关键算法与函数说明
1. 干扰生成算法 (generate_jamming)
该函数实现了频域滤波法来生成带限噪声:
- 首先生成宽带复高斯白噪声。
- 通过FFT变换到频域,利用矩形窗函数截取指定的干扰带宽(
Jamming_Bj)并进行频移。 - 通过IFFT变换回时域得到基带噪声 $u(t)$。
- 根据干扰类型进行调制:
*
RF Noise (射频噪声):直接输出带限噪声。
*
AM Noise (调幅噪声):构造包络 $A(1 + m cdot u(t))$ 乘以载波项。
*
FM Noise (调频噪声):对噪声实部积分作为相位调制项 $exp(j(2pi f_j t + 2pi K_{fm} int u(tau)dtau))$。
2. CA-CFAR 检测算法 (cf_ar_algo)
实现了经典的单元平均CFAR:
- 输入:脉压后的各距离单元幅度。
- 结构:包含检测单元、保护单元(防止目标能量泄漏到参考单元)和参考单元。
- 门限计算:$T = alpha times Mean(Reference_Cells)$。其中乘积因子 $alpha$ 由虚警概率 $P_{fa}$ 和参考单元数 $N$ 决定,公式为 $alpha = N(P_{fa}^{-1/N} - 1)$。
3. 多普勒处理
在生成目标回波时,显式引入了多普勒相位项
phase_doppler。这保证了仿真不仅能验证距离维度的压缩,在扩展到多脉冲处理时也能验证速度维度的特性(虽然本代码主要展示单脉冲结果)。