调频连续波体制雷达 基本数据处理

调频连续波体制雷达 基本数据处理

产生 速度-距离-强度 三维图

clear all;close all;clc;
% 雷达参数
Tp = 102.4e-6;           %(s) 调频周期
Fs=2.5e6;                %(Hz) 采样率
c = 3e8;                 %(m/s) 光速
fstart = 78e9;           %(Hz) 工作频率
fstop = 80e9;            %(Hz) 工作频率
BW = fstop-fstart;       %(Hz) 带宽
K = BW/Tp;               %调频斜率
f0 = (fstop+fstart)/2;   %(Hz) 调制频率
lambda = c/f0;           %(m) 波长

% 读出数据
Rx = ReadFuntion('<PATH>/<FileName>');	%天线号，采样点序号，Chirp数，帧数
PointNumber = size(Rx,2); 				%每个调频周期（chirp）采样点数256
ChirpNumber = size(Rx,3); 				%每帧（frame）包含chirp数    256
FrameNumber = size(Rx,4); 				%帧数
Data = Rx(1,:,:,:); 					%取出一个天线的接收信号
Data = squeeze(Data);

x = 1:PointNumber
distance =  (x-128)*2.5*1e6/256*102.4e-6/2e9*3e8/2
velocity = (x-128) * 0.0038/2/256/102.4e-6

% 回波显示
for i = 1:FrameNumber
    Y = abs(fft2(Data(:,:,i)))
    Y = fftshift(Y)
    imagesc(-velocity,distance,Y)
    ylim([0 9.525])
    xlabel('速度')
    ylabel('距离')
    title(['Frame:' num2str(i) '/' num2str(FrameNumber)]);
    pause(0.01);
end

横轴为速度，速度大于0是靠近雷达，反之远离
数据处理结果 一个人靠近雷达再远离雷达
速度为0 距离为7-8之间的信号是屋门