【雷达信号处理】基于matlab雷达信号采集数据处理（时频谱图）【含Matlab源码 3410期】.zip

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('时间/(s)'), ylabel('频率/(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