add calibration script to be run before running labs

Author: George Mencoff

Diff for: demos/radar/helpers/arrangePulseData.m

@@ -3,6 +3,10 @@
 %
 % Copyright 2023 The MathWorks, Inc.
 
+% Combine data from channels with calibration weights
+calweights = loadCalibrationWeights().DigitalWeights;
+indata = indata * conj(calweights);
+
 % Extract timing from pluto and phaser setup
 fs = rx.SamplingRate;
 tsweep = double(bf.FrequencyDeviationTime) / 1e6;

Diff for: shared/calibration/CalibrationValueFormat.m

@@ -5,8 +5,8 @@
     % Copyright 2023 The MathWorks, Inc.
 
     properties
-        AnalogWeights
-        DigitalWeights
+        AnalogWeights = ones(4,2);
+        DigitalWeights = ones(2,1);
     end
 
     methods
@@ -15,5 +15,4 @@
             s.DigitalWeights = this.DigitalWeights;
         end
     end
-end
-
+end

Diff for: shared/calibration/generateCalibrationWeights.m

@@ -18,9 +18,13 @@ function generateCalibrationWeights()
 % Find the hb100 center frequency
 fc_hb100 = findTxFrequency();
 
+pause(5);
+
 % Save the gain profile for the antenna.
 saveGainProfile(fc_hb100);
 
+pause(5);
+
 % Get the calibration data
 CalibrationData = calibrationRoutine(fc_hb100);
 

Diff for: shared/calibration/loadCalibrationWeights.m

@@ -3,6 +3,11 @@
 
 filepath = fileparts(which('generateCalibrationWeights'));
 calweights_filename = [filepath,'\','CalibrationWeights.mat'];
-load(calweights_filename,'calibrationweights');
+if isfile(calweights_filename)
+    load(calweights_filename,'calibrationweights');
+else
+    calibrationweights = CalibrationValueFormat;
+end
+    
 
 end

Diff for: shared/phasercontrol/helperGainCodes.m

@@ -1,7 +1,7 @@
 function calibGainCode = helperGainCodes(analogWeights)
     % Get gain code for each subarray
 
-    % Copyright 2023 The MathWorks, Inc.
+    % Copyright 2023-2025 The MathWorks, Inc.
 
     sub1weights = analogWeights(:,1);
     sub2weights = analogWeights(:,2);
@@ -10,7 +10,11 @@
 
     % Loading the measured gain profiles. These are captured using the
     % script saveGainProfile.m.
-    load('GainProfile.mat','subArray1_NormalizedGainProfile','subArray2_NormalizedGainProfile','gaincode'); 
+    try
+        load('GainProfile.mat','subArray1_NormalizedGainProfile','subArray2_NormalizedGainProfile','gaincode'); 
+    catch
+        error('No GainProfile.mat file was found. Please run calibration');
+    end
 
     calibGainCode = zeros(1,8);
     for nch = 1 : 4

Diff for: workshop/Calibrate Antenna.mlx

@@ -29,10 +29,12 @@
 tx.DataSource = "DMA";
 
 % Setup beamformers all to max gain with no phase shifts
-% TODO: this is where we should include calibration coefficients
+calibrationweights = loadCalibrationWeights();
 bf = setupPhaser(rx,fc);
 bf.RxPowerDown(:) = 0;
-bf.RxGain(:) = 127;
+
+% Load calibration data into beamformer
+setAnalogBfWeights(bf,calibrationweights.AnalogWeights);
 
 % Setup Phase Locked Loop (PLL)
 bf.Frequency = (fc+rx.CenterFrequency)/4;