oxy_fit_poly2

0001 function [new_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2]=...
0002     oxy_fit_poly2(volt,OXSAT,temp,press,sta,bot_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2);
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)=p1;
0010 beta(1,7)=p2;
0011 beta0=non_lin_fit(x,bot_oxy,'oxygen_fit_poly2',beta);
0012 Soc=beta0(1);
0013 Voffset=beta0(2);
0014 Boc=beta0(3);
0015 tcor=beta0(4);
0016 pcor=beta0(5);
0017 p1=beta0(6);
0018 p2=beta0(7);
0019 new_oxy=((Soc.*(volt+Voffset))+Boc.*exp(-0.03.*temp)+(p1.*sta+p2.*sta.^2)).*exp(tcor.*temp).*OXSAT.*exp(pcor.*press);
0020 return