✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
在现代脉冲雷达系统中,相位编码信号以其较好的抗干扰性能,越来越被重视和使用.MATLAB作为一种仿真工具,经常被用于雷达信号处理方案设计中.本文用MATLAB对相位编码信号的信号处理过程进行仿真,对信号处理过程中各节点信号进行分析,为雷达系统的总体设计提供了参考依据.
⛄ 完整代码
clear, clc
%% 绘制雷达信号处理结果图
% 定义文件名称
filenames = {'fall_aside_Sxx_1602.csv', ...
'fall_backward_Sxx_1602.csv', ...
'fall_forward_Sxx_1601.csv', ...
'jump_Sxx_1601.csv', ...
'sit_Sxx_1602.csv', ...
'walk_Sxx_1602.csv'};
% 定义标题名称
titles = {'(a)', '(b)', '(c)', '(d)', '(e)', '(f)'};
% 读取时间向量和频率刻度
t = load('./radar_time.csv');
f = load('./radar_freq.csv');
% 初始化画布
figure(1)
set(gcf, 'position', [100, 100, 1200, 600])
i = 0;
% 定义均值滤波器
filter = fspecial('average', 3);
% 读取并绘制时频谱图
for filename = filenames
Sxx = load(strcat('./data/', filename{1}));
Sxx = reshape(Sxx, [128, 128, 4]);
subplot(231+i)
Stf = Sxx(:, :, 2);
Stf_shifted = fftshift(Stf, 2)';
%Stf_filted = imfilter(Stf_shifted, filter, 'replicate');
imagesc(t, -f/1e3, 2*Stf_shifted)
xlabel('Time/(s)'), ylabel('Frequency/(kHz)')
title(titles{i+1})
set(gca, 'YDir', 'normal')
set(gca, 'XTick', 0:0.5:2.5)
set(gca, 'XTicklabel', {'0', '0.5', '1.0', '1.5' '2.0', '2.5'})
colorbar
colormap jet
i = i + 1;
end
⛄ 运行结果
⛄ 参考文献
[1] 宋光然, 侯庆文, 倪梓明. 基于Faster R-CNN的雷达时序信号频谱图的料线检测[J]. 2019.
[2] 郭利荣, 何明浩, 郁春来,等. 基于频谱图像中心矩特征的雷达信号分选方法[J]. 探测与控制学报, 2015, 37(2):5.
[3] 殷俊丽, 丁康利, 郝鹏飞. 基于MATLAB的雷达信号处理仿真[J]. 电子技术与软件工程, 2017(18):1.
[4] 成鸿飞, 王莹, 段柳云,等. 基于Matlab的雷达系统仿真软件包设计与实现[J]. 现代电子技术, 2010, 33(16):4.