fitnlorentzian(xyda​ta, parameters)

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

fits n number of lorentizian peaks to xy data with automated peak detection (tweak minpeakheight and minpeakwidth for better detection)
see also https://www.mathworks.com/matlabcentral/fileexchange/63771-fitngaussian-xy--parameters-

how to use:

1. Download fitnlorentzian.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;
fitnlorentzian(xydata, n)

code is maintained at

https://github.com/ebauch/matlab