Parameter estimation in kinetic study of reforming reaction involving ODEs
3 views (last 30 days)
Show older comments
Hi!
I am trying to estimate 13 unknown parameters involved in 4 ordinary differential equations (ODEs) by fitting it to experimental data. There are six variables (4 dependent variables and 2 independent variables). Eight set of data points are available. Equations are as follows.
- dYB/dt = -(r1) - (r2)
- dYH2/dt = (8*r1) + (12*r2) + (r3) - (r4)
- dYCO/dt = (4*r1) - (r3) + (r4)
- dYCO2/dt = (4*r2) + (r3) - (r4)
Where r1 = k10*e^(-E1/(RT))*(PB^a), r2 = k20*e^(-E2/(RT))*(PB^b), r3 = k30*e^(-E3/(RT))*(PCO^c), r4 = k40*e^(-E4/(RT))*(PCO2^d)*(PH2^e).
I tried to estimate parameters using lsqcurvefit and ode15s, but there are errors which I am not able to rectify.
Thank you in advance.
5 Comments
Torsten
on 14 Jun 2019
V has 4 elements, namely [YH, YH2,YCO, YCO2].
What are you trying to do when you refer to V(5,:) ?
Answers (0)
See Also
Categories
Find more on Event Functions in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!