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?c2?=?c2（la2）;
c24=c24（la2）;

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