source: MML/trunk/machine/SOLEIL/Booster/applications/fitprofil.m @ 4

% profil de perte

function [estimates, model] = fitprofil(xdata, ydata)
% Call fminsearch with a random starting point.
start_point = [2 0 90]; % sig et zb et amplitude
model = @expfun;
estimates = fminsearch(model, start_point);
% expfun accepts curve parameters as inputs, and outputs sse,
% the sum of squares error for A * exp(-lambda * xdata) - ydata, 
% and the FittedCurve. FMINSEARCH only needs sse, but we want to 
% plot the FittedCurve at the end.
    function [sse, FittedCurve] = expfun(params)
        sig = params(1);
        zb  = params(2);
        A   = params(3);
        chambre=10;
        z   =xdata;
        
        ch=min(abs(chambre-z-zb),abs(-chambre-z-zb));
        for i=1:length(z)
            x=z(i);
            if (x+zb)>chambre
                ch(i)=0;
            elseif (x+zb)<-chambre
                ch(i)=0;
            end
        end
        profil=erf(ch/sqrt(2)/sig);
        FittedCurve = A*(profil);
        ErrorVector = FittedCurve - ydata;
        sse = sum(ErrorVector .^ 2);
    end
end