Double monod maximum rate estimation
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Maria Teresa Aguirrezabala Campano
am 26 Jan. 2023
Kommentiert: Torsten
am 26 Jan. 2023
I want to model a rate as a double monod, and I have the experimental data of two subtrates in time. I only need to estimate one parameter (maximum rate) and have been checking my notes from 10 years ago, browsing answers here, but none are similar to my case and I don't know how to proceed (I also forgot how these types of fittings are done). I need to use the data of both substrates in time for my rate to make sense, but I cannot understand how to implement that into ODE. Perhaps it's because I shouldn't be using ODE, but I'm just lost. Below are my equations.
2 Kommentare
Akzeptierte Antwort
Torsten
am 26 Jan. 2023
Bearbeitet: Torsten
am 26 Jan. 2023
% Initial guess for rmax
rmax0 = ...;
% Time is a (nx1) vector,
% Sdata is a (nx2) matrix with measurements of c_CH4 (1st column) and c_O2 (2nd column)
Time = ...;
Sdata = ...;
% Call fitting function
[B,Rsdnrm,Rsd,ExFlg,OptmInfo,Lmda,Jmat] = lsqcurvefit(@(x,xdata)MonodKinetics(x,xdata,Sdata),rmax0,Time,Sdata);
% Define kinetic equations
function S = MonodKinetics(rmax, Time, Sdata)
x0 = [Sdata(1,1);Sdata(1,2)]; % (2x1) vector of initial conditions for c_CH4 and c_O2
[T,S] = ode45(@DifEq, Time, x0);
function xdot = DifEq(t, x)
F = ...;
MG = ...;
K_CH4 = ...;
K_O2 = ...;
MTT = rmax*x(1)/(x(1)+K_CH4)*x(2)/(x(2)+K_O2);
R = 0.2*MTT;
xdot = zeros(2,1);
xdot(1) = MG - MTT - F;
xdot(2) = R + MTT;
end
end
2 Kommentare
Torsten
am 26 Jan. 2023
You should plot the solution curves for several values of the parameter rmax in order to get a feeling in how far this parameter can influence the overall behaviour of C_CH4 and C_O2 if all other parameters remain constant.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Interpolation finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!