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
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
radar.zip (27个子文件)
radar
data
jump_Sxx_2.csv 515KB
sit_Sxx_1601.csv 515KB
fall_backward_Sxx_2.csv 514KB
sit_Sxx_2.csv 516KB
fall_backward_Sxx_1.csv 516KB
jump_Sxx_1602.csv 513KB
fall_forward_Sxx_1602.csv 513KB
fall_aside_Sxx_1.csv 514KB
sit_Sxx_1602.csv 514KB
sit_Sxx_1.csv 515KB
fall_aside_Sxx_2.csv 513KB
fall_aside_Sxx_1602.csv 513KB
walk_Sxx_1602.csv 513KB
fall_backward_Sxx_1601.csv 515KB
fall_forward_Sxx_2.csv 513KB
fall_forward_Sxx_1.csv 513KB
jump_Sxx_1.csv 513KB
fall_forward_Sxx_1601.csv 513KB
fall_backward_Sxx_1602.csv 515KB
jump_Sxx_1601.csv 513KB
walk_Sxx_1601.csv 513KB
walk_Sxx_2.csv 513KB
fall_aside_Sxx_1601.csv 512KB
walk_Sxx_1.csv 513KB
main.m 1KB
radar_time.csv 768B
radar_freq.csv 1KB
共 27 条
- 1
资源评论
洋洋Young
- 粉丝: 4534
- 资源: 25
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功