光谱拟合matlab代码

function [II_nI_nmI_rRSD1RSD2RSD3] = calculate1（aySsepstart）


tic
for n=1:520
value=lsqcurvefit（@lorentzpstarta（s:e）y（s:en）[0 0 pstart（3）-5][10^8 pstart（2）*5 pstart（3）+5] optimset（‘Display‘‘off‘））;
c24（n）=value（1）;
end
toc

c2 = c24./S;
cc = sort（c2）;
la2 = find（（c2>cc（60）） & （c2 < cc（461）））;
c2 = c2（la2）;
c24=c24（la2）;

S = S（la2）;
mc=mean（c24）;
mS=mean（S）;