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资源评论
洋洋Young
- 粉丝: 1954
- 资源: 13
