fit_abs_single.m

function [fit_parameters,igout,plotdata] = call_fit_abs() 

root='C:\Users\Manfred Beutel\Dropbox\Masterarbeit\Messung\NbSTO(0.7)_4.2.2016\measurement2\'
[filename,pathname]=uigetfile({'*.dat'},'title',root,'MultiSelect','off');
lst=importdata(strcat(pathname,filename));
lst=lst.data;
lst(:,1)=lst(:,1)/10^9;

figure(1);
plot(lst(:,1),lst(:,4).^2,'o','MarkerSize',2);
[xpts,ypts]=getpts(figure(1));
lb=find(lst(:,1)>xpts(1),1,'first');
ub=find(lst(:,1)>xpts(2),1,'first');
lst=lst(lb:ub,:);
figure(1)
plot(lst(:,1),lst(:,4).^2,'o','MarkerSize',2);drawnow;
pause(2);


colabs=4;

IL=max(lst(:,4)).^2;
f0ind=find(lst(:,4)==max(lst(:,4)));
f0=lst(f0ind);
[indrowl,indcoll]=find(lst(:,colabs)>max(lst(:,colabs))/2,1,'first');
[indrowr,indcolr]=find(lst(:,colabs)>max(lst(:,colabs))/2,1,'last');
fB=lst(indrowr,1)-lst(indrowl,1);
if fB==0
    fB=100/1000000000;
end;
y0=lst(10,colabs);
ig=[IL*fB,fB,f0,y0]

lb=[0, 0, 0, 0];
ub=[100,60,4,10];
%fit

objfcn=@(v,f)(2/pi)*v(1)*v(2)./(4*(f-v(3)).^2+v(2)^2)+v(4); %L0=v(1);fB=v(2);f0=v(3);y0=v(4)
%right fit function: objfcn=@(v,f)(2/3.14159)*v(1)*v(2)./(v(2)^2+4*(f-v(3)).^2)+v(4);
%to find init guess for amplitude use a0/fB, actually it should be a0*fB
opt=optimoptions(@lsqcurvefit,'Algorithm','levenberg-marquardt','MaxFunEvals',100000);
[vfit,resnorm] = lsqcurvefit(objfcn,ig,lst(:,1),lst(:,colabs).^2,[],[],opt);
objfcn_fit=objfcn(vfit,lst(:,1));
igout=[vfit(1),vfit(2),vfit(3),vfit(4),1,1,1,1,1,1];    

f0fit=vfit(3);fBfit=vfit(2);Qfit=f0fit/fBfit;
fit_parameters=[f0fit,fBfit,Qfit];

plotdata=cat(2,[],lst(:,1),objfcn_fit);


figure(3);
plot(lst(:,1),lst(:,4).^2,'o','MarkerSize',2);
hold on;
plot(lst(:,1),objfcn_fit,'LineWidth',1.5);
l=legend('abs data_i=1','fit');
set(l,'Interpreter','none');
hold off;