Record Raw ADC Data for Offline Processing from TI mmWave Radar Board Using DCA1000EVM Capture Card

dca = dca1000("IWR6843ISK");

clear dca
% Create connection to the TI Radar board and DCA1000EVM Capture card
dca = dca1000("IWR6843ISK");
% Specify the duration to record ADC data
dca.RecordDuration = 100;
% Specify the location at which you want to store the recorded data along
% with the recording parameters
dca.RecordLocation = "C:\TIRadarADCData\dca1000Data";
% Start recording.
% The function startRecording opens a window. Ensure that you do not
% close this window. It will automatically close when the recording
% finishes.
startRecording(dca);
% The startRecording function captures data in background.
% While this is happening, you are free to utilize MATLAB for any other
% tasks you might need to perform. The following code is designed to
% prevent MATLAB from proceeding until the recording has finished.
% The isRecording() function will return true if the recording is still
% in progress. Once the recording has concluded or if it has not started,
% the function will return false.
while isRecording(dca)
end
% Remember the record location for post-processing.
% In this example, we will save the recording location in a variable.
recordLocation = dca.RecordLocation;
% Clear the dca1000 object and remove the hardware connections if required
clear dca

% Load the recording parameters and Radar Configurations from the record
% location. This data is stored as a .mat file in the record location
temp = load(fullfile(recordLocation,'iqData_RecordingParameters.mat'));
dcaRecordingParams = temp.RecordingParameters;
% Define a variable to set the sampling rate in Hz for the
% phased.RangeResponse object. Because the dca1000 object provides the
% sampling rate in kHz, convert this rate to Hz.
fs = dcaRecordingParams.ADCSampleRate*1e3;
% Define a variable to set the FMCW sweep slope in Hz/s for the
% phased.RangeResponse object. The dca1000 object provides the
% sweep slope in GHz/us; convert this sweep slope to Hz/s.
sweepSlope = dcaRecordingParams.SweepSlope * 1e15;
% Define a variable to set the number of range samples
nr = dcaRecordingParams.SamplesPerChirp;
% Define a variable to set the center frequency in Hz for the
% phased.RangeDopplerScope object. The dca1000 object provides the
% center frequency in GHz; convert this rate to Hz.
fc = dcaRecordingParams.CenterFrequency*1e9;
% Define a variable to set the pulse repetition period in seconds for the
% phased.RangeDopplerScope object. Because the dca1000 object provides the
% chirp cycle time in us, convert this rate to seconds. The value is
% multiplied by 2 since in the example we are plotting the data only for first
% transmit channel out of 2 channels
tpulse = 2*dcaRecordingParams.ChirpCycleTime*1e-6;
% Pulse repetition interval
prf = 1/tpulse;
% Number of active receivers
nrx = dcaRecordingParams.NumReceivers;
% Number of chirps
nchirp = dcaRecordingParams.NumChirps;
% Create phased.RangeResponse System object that performs range filtering
% on fast-time (range) data, using an FFT-based algorithm
rangeresp = phased.RangeResponse(RangeMethod = 'FFT', ...
    RangeFFTLengthSource = 'Property', ...
    RangeFFTLength = nr, ...
    SampleRate = fs, ...
    SweepSlope = sweepSlope, ...
    ReferenceRangeCentered = false);
% Create range doppler scope to compute and display the response map.
rdscope = phased.RangeDopplerScope(IQDataInput=true, ...
    SweepSlope = sweepSlope,SampleRate = fs, ...
    DopplerOutput="Speed",OperatingFrequency=fc, ...
    PRFSource="Property",PRF=prf, ...
    RangeMethod="FFT",RangeFFTLength=nr, ...
    ReferenceRangeCentered = false);
% Create a dca1000FileReader object that enables you to read ADC data (IQ data) cubes
% from the binary files recorded using the object DCA1000 in MATLAB
fr = dca1000FileReader(recordLocation = recordLocation);
% Execute the loop till all the IQ data cubes are read from the record location
while (fr.CurrentPosition <= fr.NumDataCubes)
% Read a single IQ data cube from the binary file,
% starting from the location previously accessed.
iqData = read(fr,1);
% The function read() return cell array of IQ data cube,
% extract the radar data cube from it
iqData = iqData{1};
% Get the data from first receiver antenna
iqDataRx1 = squeeze(iqData(:,1,:));
% Plots the range response corresponding to the input signal, iqData.
plotResponse(rangeresp, iqDataRx1);
% Get the data from first receiver antenna and first transmitter antenna
% (Every alternate pulse is from one transmitter antenna for this
% configuration)
iqDataRx1Tx1 = squeeze(iqData(:,1,1:2:end));
% Plots the range doppler response corresponding to the input signal,
% iqData.
rdscope(iqDataRx1Tx1);
end