fitting routines(1)

4
FO consecutive intermediate, 4 parameters Equation f =A20*exp(-k2*t)+ k1*A10/(k1-k2)*(exp(-k2*t) - exp(-k1*t)) fit f to y ''fit f to y with weight reciprocal_y ''fit f to y with weight reciprocal_ysquare Variables t = col(26) y = col(27) reciprocal_y = 1/abs(y) reciprocal_ysquare = 1/y^2 'Automatic Initial Parameter Estimate Functions first(q)=if(size(q)<10,size(q)-1,int(0.9*size(q))) ylast(q)=mean(q[data(first(q),size(q))]) xnear0(q)=max(abs(q))-abs(q) yatxnear0(q,r)=xatymax(q,xnear0(r)) Initial parameters A20 = 0.0424 ' {{previous: 0.0439974}} A10 = 1.8 ' {{previous: 6.00835}} k1 = 0.001 ' {{previous: 0.00141284}} k2 = 0.01 ' {{previous: 0.0086247}} Constraints k1 > 0 k2 > 0 ;k1 > k2 ;A10=1.7912 ;k2=0.0121 ;k1=3.8949e-3 Iterations 200 Step size 1 Tolerance 1e-10

Upload: fredrick-mutunga

Post on 29-Sep-2015

217 views

Category:

Documents


0 download

DESCRIPTION

sigmaplot

TRANSCRIPT

FO consecutive intermediate, 4 parameters

Equationf =A20*exp(-k2*t)+ k1*A10/(k1-k2)*(exp(-k2*t) - exp(-k1*t))fit f to y''fit f to y with weight reciprocal_y''fit f to y with weight reciprocal_ysquare

Variablest = col(26)y = col(27)reciprocal_y = 1/abs(y)reciprocal_ysquare = 1/y^2'Automatic Initial Parameter Estimate Functionsfirst(q)=if(size(q) 0k2 > 0;k1 > k2;A10=1.7912;k2=0.0121;k1=3.8949e-3

Iterations200Step size1Tolerance1e-10

FO consecutive product, 3 parameters

Equationf = A10*(1 - k2/(k2-k1)*exp(-k1*t) + k1/(k2-k1)*exp(-k2*t))fit f to y''fit f to y with weight reciprocal_y''fit f to y with weight reciprocal_ysquare

Variablest = col(26)y = col(28)reciprocal_y = 1/abs(y)reciprocal_ysquare = 1/y^2'Automatic Initial Parameter Estimate Functionsfirst(q)=if(size(q) k1;A20=0.0389;A10=1.791;k1=3.8949e-3;k2=8.5472e-3;A10=1.8336

Iterations200Step size1Tolerance1e-10