处理雷达数据，噪声过大

Question

使用TI公司的IWR1443采集了跌倒数据，得到的谱图是这样的，出现这种情况的原因是什么，我个人认为是背景噪音过大，如何用matlab处理这种情况？

目前做的操作是读取雷达数据后做一维FFT变换，静态滤波，快速傅里叶变换得到的数据

raw_data = ReadDataFile(filepath);% 读取数据文件 raw_data：原始数据，按列排放
lvds_data = Real2Complex(raw_data).';      
radar_data = reshape(lvds_data,NumofSamples,NumofChirps*NumofFrames,n_RX);  % 维度信息
radar_data_avg = sum(radar_data,3)/n_RX;                                       % 通道非相干积累
radar_data_avg = radar_data_avg-mean(radar_data_avg,1);                     % 直流偏置
data3D_avg = reshape(radar_data_avg,NumofSamples,NumofChirps,NumofFrames);  % 维度信息

%-----时频处理-----
%--------------------------------------------------
range_win = hamming(NumofSamples);     % 加海明窗
range_profile = zeros(NofPoints_1DFFT,NumofChirps*NumofFrames);
for i=1:NumofChirps*NumofFrames
  temp=radar_data_avg(:,i).*range_win;    %加窗函数
  temp_fft=fft(temp,NofPoints_1DFFT);                 %对每个chirp做N点FFT
  range_profile(:,i)=temp_fft;
% range_profile_cfar(:,m)=CFAR((temp_fft));
end

rangecell_avg = sum(range_profile(1:NofPoints_1DFFT/2,:));
rangecell_avg = rangecell_avg-mean(rangecell_avg,2);%静态杂波滤除，平均相消法

TF = FastTFanls(rangecell_avg); %快速时频分析，使用什么方法在函数内调整

figure
colormap(jet(256))
imagesc([0,Ttotal],[-Fds/2,Fds/2],20*log10(fftshift(abs(TF),1)+eps))
xlabel('Time (s)')
ylabel('Doppler (Hz)')
title('Doppler Signature of Human Walk')
axis xy
colorbar
drawnow

应该在何处去噪，最好可以提供相关代码

Answer 1

根据您提供的信息，出现这种谱图的原因可能是背景噪音过大。要处理这种情况，您可以使用MATLAB中的滤波器对数据进行降噪处理。常用的滤波器包括均值滤波器、中值滤波器、高斯滤波器等。您可以根据实际情况选择合适的滤波器进行处理。另外，您还可以尝试使用时域和频域的方法对数据进行处理，以达到更好的降噪效果。