八种典型雷达辐射源信号仿真系统
项目介绍
本项目是一个基于MATLAB开发的综合性雷达辐射源信号仿真平台。系统能够生成电子战和信号处理领域中最具代表性的八种雷达信号,模拟真实的电磁传播环境(添加高斯白噪声),并结合时频分析工具对生成信号的脉内调制特征进行可视化验证。该平台旨在为雷达信号截获、特征提取及调制识别算法的研究提供标准化的数据源。
功能特性
- 多类型信号生成:支持生成8种特定调制类型的雷达信号,包括:
* 常规脉冲信号 (Normal Pulse/CW)
* 线性调频信号 (LFM)
* 非线性调频信号 (NLFM)
* 二相编码信号 (BPSK - Barker码)
* 四相编码信号 (QPSK - P4码)
* 频率编码信号 (Costas码)
* 多相编码信号 (Frank码)
* 频移键控信号 (FSK - 2FSK)
- 参数配置灵活:支持自定义采样频率、载波频率、脉冲宽度、信噪比、信号带宽等核心参数。
- 环境模拟:内置加性高斯白噪声(AWGN)模块,可根据设定的SNR值生成含噪信号。
- 信号分析与可视化:
* 时域分析:展示信号实部波形。
* 时频分析:利用短时傅里叶变换(STFT)生成时频图(Spectrogram),直观展示频率随时间的变化规律。
- 自动化流程:一键执行仿真,自动完成所有信号的生成、噪化处理及图表绘制(4x4网格布局)。
系统要求
- MATLAB R2016a及以上版本
- Signal Processing Toolbox(信号处理工具箱,用于
spectrogram 函数)
使用方法
- 打开MATLAB软件,定位到项目所在目录。
- 直接运行主程序函数。
- 程序将自动清除工作区,加载默认参数,并弹出一个绘图窗口展示8种信号的时域和时频特性。
- 如需修改仿真参数,可在主程序顶部的
SysParam 结构体中调整数值(如 Fs, fc, PW, SNR, B 等)。
---
代码实现逻辑详解
该系统主要由主控流程、信号生成核心、噪声处理模块和可视化分析模块四部分组成。
1. 全局参数配置与主控流程
主程序入口首先定义了全局系统参数结构体
SysParam,默认配置如下:
- 采样频率:100MHz
- 载波频率:10MHz
- 脉冲宽度:20us
- 信噪比:15dB
- 信号带宽:5MHz
主循环通过遍历预定义的信号类型列表({'Normal', 'LFM', 'NLFM', 'BPSK', 'QPSK', 'Costas', 'Frank', 'FSK'}),依次执行“生成 -> 加噪 -> 分析绘图”三个步骤。
2. 核心信号生成逻辑
系统通过
generate_signal 函数实现不同调制方式的数学模型:
- Normal (常规脉冲):生成单频连续波脉冲,表现为复指数信号 $exp(j cdot 2pi f_c t)$。
- LFM (线性调频):根据带宽和脉宽计算调频斜率 $k$,相位包含二次项 $t^2$,实现频率随时间线性变化。
- NLFM (非线性调频):代码实现了一种组合调制模型。在LFM的基础上叠加了一个余弦调制项 $4cos(2pi t/PW)$,形成S型频率变化曲线以改善旁瓣特性。
- BPSK (二相编码):采用 Barker-13 码序列
[+1 +1 +1 +1 +1 -1 -1 +1 +1 -1 +1 -1 +1]。将脉冲划分为13个子码元,根据码值翻转相位(0或$pi$)。 - QPSK (四相编码):采用 P4码 实现。码长设为64,相位遵循具体的二次方规律生成,本质上是LFM信号的离散相位近似。
- Costas (频率编码):使用特定的跳芬序列
[4 7 1 6 5 2 3]。将脉冲划分为7个跳频周期,每个周期内保持固定频率,实现脉内频率跳变。 - Frank (多相编码):基于 $4 times 4$ 矩阵(码长16)生成Frank码相位系数。相位步进由矩阵行列索引决定,近似于LFM特性的多相编码。
- FSK (频移键控):实现 2FSK 调制。生成长度为20的随机二进制序列,根据比特值(0或1)在 $f_c-1MHz$ 和 $f_c+1MHz$ 两个频率点之间切换。
3. 噪声添加模块
通过
add_awgn 函数实现环境模拟:
- 计算输入信号的平均功率。
- 根据用户设定的
SNR(dB值)反算噪声功率。 - 生成服从标准正态分布的复高斯噪声,并按噪声功率系数缩放。
- 将噪声叠加至原始信号上。
4. 信号分析与可视化模块
通过
analyze_and_plot 函数完成数据处理与绘图:
* 使用
spectrogram 函数进行短时傅里叶变换(STFT)。
* 使用汉宁窗(Hann window),窗口长度约为数据长度的1/64,重叠率为75%,FFT点数为512。
* 创建一个画布,通过
subplot 将其划分为 4行 x 4列 的网格。
* 每个信号占用两个子图位置:
1.
时域图:显示含噪信号实部的波形变化。
2.
时频图:显示STFT幅度谱,以颜色映射表征能量分布,横轴为时间,纵轴为频率。
* 8种信号依次排列,直观对比不同调制样式在时域和时频域的特征差异。