0
$\begingroup$

I want to do the following,

I have a set of 20 first order differential equations and I want to estimate some of the parameters. I've got the following initial and boundary conditions. The initial conditions of 20 variables (at t=0) Boundary values of 20 variables (at time , t = 2000s)

Using the shooting method, I want to estimate the values of 10 parameters using least squares minimization, using fmincon. The objective function will minimize the difference between the target(Boundary values of 20 variables (at time , t = 2000s)) and the values obtained at the end of shoot using ode15s.

I'm not sure how to implement this. Could someone help?

I have a function which solves the 10 differential equation using ode15s, how do I proceed after this step?

Edit: I have a defined afunction that computes the residual,

function residual = LQminimization(initialConditions,tSpan,y_SS)
y0 = initialConditions;
[t,y] = ode15s(@myModel,tSpan,y0);
residual= (y[:,end] - y_SS).^2;

I have reached till this point.I'm not sure how to proceed with estimating parameters from here. Any suggestions?

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.