Jump to content

LV Shared Resource: Difference between revisions

From mintOC
Line 15: Line 15:
<math>
<math>
\begin{array}{llclr}
\begin{array}{llclr}
  \displaystyle \min_{u} & \int_0^{t_f} (x_0(t) - 1.5)^2 + (x_1(t) - 1)^2 + (x_2(t) - 1)^2 dt \\[1.5ex]
  \displaystyle \min_{u} & \int_0^{t_f} && (x_0(t) - 1.5)^2 + (x_1(t) - 1)^2 + (x_2(t) - 1)^2 \ dt \\[1.5ex]
  \mbox{s.t.}  
  \mbox{s.t.}  
  & \dot{x}_0 & = &  x_0(t) - x_0(t) x_1(t) - x_0(t) x_2(t), \\
  & \dot{x}_0 & = &  x_0(t) - x_0(t) x_1(t) - x_0(t) x_2(t), \\
Line 21: Line 21:
  & \dot{x}_2 & = & -x_2(t) + \alpha x_0(t) x_2(t) - c_2 x_2(t) u(t),  \\[1.5ex]
  & \dot{x}_2 & = & -x_2(t) + \alpha x_0(t) x_2(t) - c_2 x_2(t) u(t),  \\[1.5ex]
  & x(0) &=& x_0, \\
  & x(0) &=& x_0, \\
  & u(t) &\in& [0,1],
  & u(t) &\in& [0,1], \\
  & \alpha &>&  1.
  & \alpha &>&  1.
\end{array}  
\end{array}  

Revision as of 06:58, 25 August 2025

LV Shared Resource
State dimension: 1
Differential states: 3
Discrete control functions: 5
Interior point equalities: 3


This Lotka Volterra problem with explicit inclusion of a shared resource is a variant of the Lotka Volterra fishing problem. Its dynamics are given via a three-dimensional ODE model.

Mathematical formulation

The optimal control problem is given by

minu0tf(x0(t)1.5)2+(x1(t)1)2+(x2(t)1)2 dts.t.x˙0=x0(t)x0(t)x1(t)x0(t)x2(t),x˙1=x1(t)+x0(t)x1(t)c1x1(t)u(t),x˙2=x2(t)+αx0(t)x2(t)c2x2(t)u(t),x(0)=x0,u(t)[0,1],α>1.

Parameters

These fixed values are used within the model.

[t0,tf]=[0,12],(c0,1,c1,1)=(0.2,0.1),(c0,2,c1,2)=(0.4,0.2),(c0,3,c1,3)=(0.01,0.1),(c0,4,c1,4)=(0,0),(c0,5,c1,5)=(0.1,0.2).

Reference Solutions

If the problem is relaxed, i.e., we demand that w(t) is in the continuous interval [0,1] rather than being binary, the optimal solution can be determined by means of direct optimal control.

The optimal objective value of the relaxed problem with nt=12000,nu=150 is x2(tf)=0.345768563. The objective value of the solution with binary controls obtained by Combinatorial Integral Approximation (CIA) is x2(tf)=0.348617982.