Not enough input arguments (2024)

tic

%% Preparation

clc; clear

data = importdata("FCPIB-293K-2.5mW-400nm-Jan072021 -ibg -bg -chirp.csv"); % insert file path within parenthesis

load Steady_State_Parameter_Values.mat

z = load("Steady_State_Parameter_Values.mat");

%% Preamble

% Fundamental constants

h = 4.1356677*10^-15; % units: eV/ Hz

c = 3*10^8; % SI units

kB = 8.617333268*10^-5; % units: eV/ K

% Parameters from steady state fitting

A1 = p(1,1);

A2 = p(1,2);

Eg = p(1,3);

Eb = p(1,4);

R = p(1,5);

g = p(1,6);

% Data

Wavelength = data(:, 1);% units: nm

E = (h*c)./(Wavelength*10^-9);

delay_t = data(1, :);

carrier_T = [1198.8, 816.7, 446.8, 328.7];

col1 = 56;

col2 = 63;

col3 = 74;

col4 = 87;

% Data for fitting

Range_E = E >= 1.5 & E <= 2.0;

Range_W = Wavelength >= (h*c)/(2.0*10^-9) & Wavelength <= (h*c)/(1.5*10^-9);

E_p = E(Range_E); % selected probe energies

data_new2 = data(Range_W, [col1,col2,col3,col4]);

delta_Abs1 = data(Range_W,col1);

delta_Abs2 = data(Range_W,col2);

delta_Abs3 = data(Range_W,col3);

delta_Abs4 = data(Range_W,col4);

% Fitting function: Elliott's Model for Transient Absorption

function F = EM_TA_wR1(x, e_p, z)

P = z.p;

kB = 8.617333268*10^-5; % units: eV/ K

A1 = P(1,1);

A2 = P(1,2);

Eg = P(1,3);

Eb = P(1,4);

R = P(1,5);

g = P(1,6);

carrier_T = [1198.8, 816.7, 446.8, 328.7];

for i = numel(e_p)

E_p = e_p(i);

F(i) = x(5)*A1.*((2.*pi.*sqrt(x(2)))./E_p).*1/g.*(integral(@(E)sech(((E_p - E)./g)).*(1 + 10.*R.*(E - x(1)) + ...

126.*(R).^2.*(E - x(1)).^2)./(1 - exp(-2.*pi.*sqrt(x(2)./(E - x(1))))), x(1), Inf, 'ArrayValued', 1)).*(1 - 1./(1 + exp((E_p - x(4))./(kB.*carrier_T(1,1))))).^2 - ...

A1.*(2.*pi.*sqrt(Eb)/E_p).*(1/g).*...

(integral(@(e)sech(((E_p - e)./g)).*(1 + 10.*R.*(e - Eg) + ...

126.*(R).^2.*(e - Eg).^2)./(1 - exp(-2.*pi.*sqrt(Eb./(e - Eg)))), Eg, Inf, 'ArrayValued', 1)) + ...

x(6).*A2.*(4.*pi.*(x(2)).^3/2).*1/x(3).*(...

(1/1^3).*sech((E_p - x(1) + x(2)./1^2)./x(3)) + ...

(1/2^3).*sech((E_p - x(1) + x(2)./2^2)./x(3)) + ...

(1/3^3).*sech((E_p - x(1) + x(2)./3^2)./x(3)) + ...

(1/4^3).*sech((E_p - x(1) + x(2)./4^2)./x(3)) + ...

(1/5^3).*sech((E_p - x(1) + x(2)./5^2)./x(3)) + ...

(1/6^3).*sech((E_p - x(1) + x(2)./6^2)./x(3)) + ...

(1/7^3).*sech((E_p - x(1) + x(2)./7^2)./x(3))) - ...

A2.*(4.*pi.*Eb.^3/2).*1/g.*(...

(1./1.^3).*sech((E_p - Eg + Eb./1.^2)./g) + ...

(1./2.^3).*sech((E_p - Eg + Eb./2.^2)./g) + ...

(1./3.^3).*sech((E_p - Eg + Eb./3.^2)./g) + ...

(1./4.^3).*sech((E_p - Eg + Eb./4.^2)./g) + ...

(1./5.^3).*sech((E_p - Eg + Eb./5.^2)./g) + ...

(1./6.^3).*sech((E_p - Eg + Eb./6.^2)./g) + ...

(1./7.^3).*sech((E_p - Eg + Eb./7.^2)./g));

end

F = F(:);

end

% Solver

lb = [Eg, Eb, g, 0.3, 0.5, 0.2]; ub = [55, 0.05, 0.05, 20, 2, 1];

x0 = [1.65, 0.03, 0.03, 1.3, 1, 0.3];

optim_lsq = optimoptions('lsqcurvefit', 'Algorithm', 'levenberg-marquardt', 'MaxFunctionEvaluations',10^5, 'MaxIterations', 10^5, 'FunctionTolerance',10^-10, 'StepTolerance', 10^-10);

delta_Abs(:, 1) = data_new2(:, 1);

carrierT = carrier_T(1, 1);

x(1, :) = lsqcurvefit(@(x, e_p) EM_TA_wR1, x0, E_p, delta_Abs(:, 1), lb, ub, optim_lsq);

plot(E_p, delta_Abs(:, 1), 'o')

hold on

plot(E_p, EM_TA_wR1(x(1, :), E_p, Z), 'LineWidth', 4.0)

hold off

for n = 2:4

delta_Abs(:, n) = data_new2(:,n);

carrierT = carrier_T(1, n);

[x(n, :), residualnorm, residual, exitflag, output, lambda, jacobian] = lsqcurvefit(@(x, E_p) EM_TA_wR(x, E_p, Z), x(n-1, :), E_p, delta_Abs(:, n), lb, ub, optim_lsq);

plot(E_p, delta_Abs(:, n), 'o')

hold on

plot(E_p, EM_TA_wRn(x(n, :), E_p, Z), 'LineWidth', 4.0)

xlabel('Probe Photon Energy (eV)')

ylabel('\Delta A (O.D.)')

legend('0.5 ps', '1.0 ps', '2.0 ps', '4.0 ps', 'Location', 'southeast')

end

Bearbeitet: Stephen23 am 11 Jul. 2024 um 11:54