fmincon maximise ODE Output

Hi

Im busy with a problem and I think its closely related to https://www.mathworks.com/matlabcentral/profile/authors/733182-jason-nicholson

I have a series reaction (A->B->C ) occurring in a CSTR and this can be described by the following differential equations (Component balance):

dCa/dt=(F/V)*(Cao)-(F/V)*Ca-kao*exp(-Ea/(R*T))*Ca

dCc/dt=-(F/V)*Cc+kco*exp(-Ec/(R*T))*Cb

dCb/dt= -(F/V)*Cb-kco*exp(-Ec/(R*T))*Cb+kao*exp(-Ea/(R*T))*Ca

I need to obtain the T and V parameters (Temperature and Volume) that maximise the production of Cb.

Can you help?

Kind regards

[Star Strider]

Here is one method of using a differential equation solver with an optimisation function (here lsqcurvefit): Monod kinetics and curve fitting. You can adapt it to use your differential equation and constraints with fmincon. Remember that to maximise a function you have to give the optimisation functions the negative of the function output. Depending on your constants, you may want to use a stiff solver such as ode15s.

I’m not listing this as an Answer because there is not enough information to provide a specific reply.

Answer 1

[Alan Weiss]

It is possible that you can adapt this example to your problem. You might be able use fmincon as a solver, too, keeping in mind these considerations, especially on taking larger-than-default finite differences.

Good luck,

Alan Weiss

MATLAB mathematical toolbox documentation