fit multiple positive and negative gaussian surface peaks that changes in time with shared kinetic parameters in matlab

main.m
function y = gaussianSurfaceModel(x, t, parameters)
    % Initialize output
    y = zeros(size(t));
    
    % Extract parameters
    amplitudePositive1 = parameters(1);
    amplitudeNegative1 = parameters(2);
    sigmaPositive1 = parameters(3);
    sigmaNegative1 = parameters(4);
    % ... add more parameters as needed
    
    % Compute contributions from each peak
    y = y + amplitudePositive1 * exp(-((t-x(1)).^2) / (2*sigmaPositive1^2));
    y = y + amplitudeNegative1 * exp(-((t-x(1)).^2) / (2*sigmaNegative1^2));
    % ... add more peaks as needed
end
563 chars
17 lines

main.m
% Define initial parameter values
initialParameters = [initialAmplitudePositive1, initialAmplitudeNegative1, initialSigmaPositive1, initialSigmaNegative1, ...];

% Perform the fit
fitParameters = lsqcurvefit(@gaussianSurfaceModel, initialParameters, t, y);
257 chars
6 lines

main.m
% Extract the fitted parameters
amplitudePositive1_fit = fitParameters(1);
amplitudeNegative1_fit = fitParameters(2);
sigmaPositive1_fit = fitParameters(3);
sigmaNegative1_fit = fitParameters(4);
% ... extract more parameters as needed

% Evaluate the model at desired time points
y_fit = gaussianSurfaceModel(fitParameters, t);
329 chars
10 lines