function that outputs a gaussian fitted value-use with nlinfit- also subtracts a baseline Z=A.*exp^(-(x-x0).^2/(2*sigma.^2)) e.g [betanew]=nlinfit(X,Y,@gaussfit,beta0)
0001 function Z = gaussfit_sig(beta,X) 0002 %function that outputs a gaussian fitted value-use with nlinfit- also 0003 %subtracts a baseline 0004 %Z=A.*exp^(-(x-x0).^2/(2*sigma.^2)) 0005 %e.g [betanew]=nlinfit(X,Y,@gaussfit,beta0) 0006 0007 xin=X;%(:,1); 0008 A= beta(1); 0009 x0 = beta(2); 0010 sigma=beta(3); 0011 m = beta(4); 0012 c = beta(5); 0013 0014 Z = A*exp(-((xin-x0).^2)./(2.*sigma^2)) - (m.*xin+c);