View source for Van der Pol Oscillator (JModelica)

This page contains the source code to solve the [[Van der Pol Oscillator]] problem with JModelica. The automatic differentiation tool CasADI and the solver IPOPT were used to solve the problem.

Model file (VDP_Opt.mop)
<source lang="Optimica">
//-------------------------------------------------------------------------
//Van der Pol Oscillator with direct collocation using JModelica
//(c) Madeleine Schroter
//--------------------------------------------------------------------------

optimization VDP_Opt(objectiveIntegrand = y^2+x^2+u^2, startTime = 0, finalTime = 20)

//The states
Real x(start=1, fixed=true); //position coordinate
Real y(start=0, fixed=true); 


//The control signal
input Real u; //damping of the oscillation

equation
der(y) = (1-x^2) * y - x +u;
der(x) = y;

constraint
u<=0.75;

end VDP_Opt;

</source>

Run file
<source lang ="Python">
#Import the function for transfering a model to CasADiInterface
from pyjmi import transfer_optimization_problem

op=transfer_optimization_problem("VDP_Opt", "VDP_Opt.mop")

res=op.optimize()

</source>

[[Category:JModelica]]