Saturation curve simulation (CW T1, T2 measurement)

[dominikk]

Dear Authors,

Is there a way to read in a set of EPR spectra measured as a function of MW power in order to calculate the saturation curve and fit the T1 and T2?

Thank you!

All best,
Dominik

[Stefan Stoll]

Use the eprload function. If your data comes as a single file, one call is enough.

Otherwise, use a for loop to read in all the files. Assuming you named them systematically, you could use something like

Code: Select all

P_dB = [0 3 9 12 15 18 21 24 27 30 33 36 39 42]; % mw power in dB

for iP = 1:numel(P_dB)
  filename = sprintf('mydata_%d_dB.DTA',P_dB(iP));  % make filename
  [x,y_] = eprload(filename);   % read in 1D data
  spc(:,iP) = y_(:);            % store data as column in a matrix
end