caiyu/PKU-Millimeter-Wave-Radar-Tutorial
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Add files via upload

Browse Source
This commit is contained in:
St. Zhang 2024-03-11 21:06:00 +08:00 committed by GitHub
parent 2abac9828b
commit ced7e84e7e
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
2 changed files with 92 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

55
基于热力图的人体定位demo/DCA1000_Read_Data.m Normal file
View File

@ -0,0 +1,55 @@
function data_radar=DCA1000_Read_Data(fname)
global parameter;
global frame;
%% 读取回波数据
fid = fopen(fname,'rb');
%16bits，复数形式(I/Q两路)，4RX,3TX,有符号16bit
sdata = fread(fid,frame*parameter.Samples*parameter.Chirps*4*3*2,'int16');
sdata = sdata((frame-1)*parameter.Samples*parameter.Chirps*4*3*2+1:frame*parameter.Samples*parameter.Chirps*4*3*2);
%% 1843+DCA1000
fileSize = size(sdata, 1);
lvds_data = zeros(1, fileSize/2);
count = 1;
for i=1:4:fileSize-5
lvds_data(1,count) = sdata(i) + 1i*sdata(i+2);
lvds_data(1,count+1) = sdata(i+1)+1i*sdata(i+3);
count = count + 2;
end
lvds_data = reshape(lvds_data, parameter.Samples*parameter.txNum*parameter.rxNum, parameter.Chirps);
lvds_data = lvds_data.';
cdata = zeros(parameter.txNum*parameter.rxNum,parameter.Chirps*parameter.Samples);
for row = 1:parameter.txNum*parameter.rxNum
for i = 1: parameter.Chirps
cdata(row,(i-1)*parameter.Samples+1:i*parameter.Samples) = lvds_data(i,(row-1)*parameter.Samples+1:row*parameter.Samples);
end
end
fclose(fid);
data_radar_1 = reshape(cdata(1,:),parameter.Samples,parameter.Chirps); %TX1 RX1
data_radar_2 = reshape(cdata(2,:),parameter.Samples,parameter.Chirps); %TX1 RX2
data_radar_3 = reshape(cdata(3,:),parameter.Samples,parameter.Chirps); %TX1 RX3
data_radar_4 = reshape(cdata(4,:),parameter.Samples,parameter.Chirps); %TX1 RX4
data_radar_5 = reshape(cdata(5,:),parameter.Samples,parameter.Chirps); %TX2 RX1
data_radar_6 = reshape(cdata(6,:),parameter.Samples,parameter.Chirps); %TX2 RX2
data_radar_7 = reshape(cdata(7,:),parameter.Samples,parameter.Chirps); %TX2 RX3
data_radar_8 = reshape(cdata(8,:),parameter.Samples,parameter.Chirps); %TX2 RX4
data_radar_9 = reshape(cdata(9,:),parameter.Samples,parameter.Chirps); %TX3 RX1
data_radar_10 = reshape(cdata(10,:),parameter.Samples,parameter.Chirps); %TX3 RX2
data_radar_11 = reshape(cdata(11,:),parameter.Samples,parameter.Chirps); %TX3 RX3
data_radar_12 = reshape(cdata(12,:),parameter.Samples,parameter.Chirps); %TX3 RX4
data_radar=[];
data_radar(:,:,1)=data_radar_1; %三维雷达回波数据
data_radar(:,:,2)=data_radar_2;
data_radar(:,:,3)=data_radar_3;
data_radar(:,:,4)=data_radar_4;
data_radar(:,:,5)=data_radar_5;
data_radar(:,:,6)=data_radar_6;
data_radar(:,:,7)=data_radar_7;
data_radar(:,:,8)=data_radar_8;
data_radar(:,:,9)=data_radar_9;
data_radar(:,:,10)=data_radar_10;
data_radar(:,:,11)=data_radar_11;
data_radar(:,:,12)=data_radar_12;
end

37
基于热力图的人体定位demo/generateParameter.m Normal file
View File

@ -0,0 +1,37 @@
%% 雷达参数设置
function parameter = generateParameter()
parameter.c = 3e8; %光速
parameter.stratFreq = 77e9; %起始频率
parameter.Tr = 60e-6; %扫频时间
parameter.Idle_time = 100e-6; %空闲时间
parameter.Tc = parameter.Tr+parameter.Idle_time; %Chirp之间的间隔
parameter.Samples = 256; %采样点
parameter.Fs = 10e6; %采样率
parameter.Tframe_set = 80e-3; % 帧周期
parameter.rangeBin = parameter.Samples ; %rangebin
parameter.Chirps = 128; %chirp数
parameter.dopplerBin = parameter.Chirps; %dopplerbin
parameter.Slope = 29.982e12; %chirp斜率
parameter.Bandwidth = parameter.Slope * parameter.Tr ; %发射信号有效带宽
parameter.BandwidthValid = parameter.Samples/parameter.Fs*parameter.Slope; %发射信号带宽
parameter.centerFreq = parameter.stratFreq + parameter.Bandwidth / 2; %中心频率
parameter.lambda = parameter.c / parameter.centerFreq; %波长
parameter.txAntenna = ones(1,3); %发射天线个数
parameter.rxAntenna = ones(1,4); %接收天线个数
parameter.txNum = length(parameter.txAntenna);
parameter.rxNum = length(parameter.rxAntenna);
parameter.virtualAntenna = length(parameter.txAntenna) * length(parameter.rxAntenna);
parameter.angleBin = 180; %anglebin
parameter.dz = parameter.lambda / 2; %接收天线俯仰间距
parameter.dx = parameter.lambda / 2; %接收天线水平间距
parameter.numCPI = 50; % 帧数
end