fitngaussian(xy, parameters)

fit n number of Gaussian peaks / sum of Gaussian peaks to data with automated peak detection

fits n number of Gaussian peaks to xy data with automated peak detection (tweak minpeakheightm, minpeakwidth for better detection)
see also https://www.mathworks.com/matlabcentral/fileexchange/63748-fitnlorentzian-xydata--parameters-
how to use:
1. Download fitngaussian.m plus example data in 4peaksdata.dat (https://github.com/ebauch/matlab/blob/master/4peaksdata.dat)
2. In matlab run
% import xy data
xydata = load('4peaksdata.dat');

% fit data with n = 1,...,4 peaks
n = 3;
fitngaussian(xydata, n)

code is maintained at

https://github.com/ebauch/matlab