oxy_fit_exp

0001 function [new_oxy,Soc,Voffset,Boc,tcor,pcor,cal,gama]=...
0002     oxy_fit_exp(volt,OXSAT,temp,press,sta,bot_oxy,Soc,Voffset,Boc,tcor,pcor,cal,gama);
0003 x=[volt OXSAT temp press sta];
0004 beta(1,1)=Soc;
0005 beta(1,2)=Voffset;
0006 beta(1,3)=Boc;
0007 beta(1,4)=tcor;
0008 beta(1,5)=pcor;
0009 beta(1,6)=cal;
0010 beta(1,7)=gama;
0011 beta0=non_lin_fit(x,bot_oxy,'oxygen_fit_exp',beta);
0012 Soc=beta0(1);
0013 Voffset=beta0(2);
0014 Boc=beta0(3);
0015 tcor=beta0(4);
0016 pcor=beta0(5);
0017 cal=beta0(6);
0018 gama=beta0(7);
0019 new_oxy=((Soc.*(volt+Voffset))+Boc.*exp(-0.03.*temp)+cal.*exp(gama*sta)).*exp(tcor.*temp).*OXSAT.*exp(pcor.*press);
0020 return