Parameter estimation in kinetic study of reforming reaction involving ODEs

1 次查看（过去 30 天）

Ghansyam Vadodaria 2019-6-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/466866-parameter-estimation-in-kinetic-study-of-reforming-reaction-involving-odes

评论： Ghansyam Vadodaria 2019-6-14

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 个评论
显示 3更早的评论隐藏 3更早的评论

Are Mjaavatten 2019-6-13

在 MATLAB Online 中打开

I am afraid this problem is quite a bit more challenging than you imagine. You should read the Matlab tutorials on differential equation solving and curve fitting.

You should always test and debug the parts before combing them in a large system:

First, you should make Difeq a separate file and have P as an extra input. Next you must debug it: You refer to V(5,:), but V should be a 1 by 4 vector. V(5) seems to be time, so maybe you should use t instead? Make sure that the results of Difeq(0,V0,P0) look reasonable.

Next try to solve the ODE for P = P0:

P = P0;
rhs = @(t,V) Difeq(0, V, P);
[t,s] = ode45(rhs,[0;[Time]],V0);

It may be wise to first try with a much smaller time interval, though.

Once this works, you may in principle be ready to try to fit the output. However, trying to fit 13 parameters with values that vary by a factor of 10^11 is not likely to succeed. It is a good idea to normalise your parameters, for instance by dividing them with your initial guess P0. Even so, fitting all 13 at once is going to be quite difficult. I suggest you try to find good guesses for most parameters and fit only one or two at first. Then increase the number of unknowns by one parameter at a time. Still, 13 seems a bit much.

Good luck! (You will need it.)

Torsten 2019-6-14

V has 4 elements, namely [YH, YH2,YCO, YCO2].

What are you trying to do when you refer to V(5,:) ?

Ghansyam Vadodaria 2019-6-14

@ Torsten, In my system I have four dependent variables (YB, YH2, YCO, YCO2) and two independent variables (Time and Temperature) V(5,:) is temperature.

请先登录，再进行评论。

请先登录，再回答此问题。