Jump to content

LV Shared Resource: Difference between revisions

From mintOC
Line 17: Line 17:
  \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(t) & = &  x_0(t) - x_0(t) x_1(t) - x_0(t) x_2(t), \\
  & \dot{x}_1 & = & - x_1(t) + x_0(t) x_1(t) - c_1 x_1(t) u(t),  \\
  & \dot{x}_1(t) & = & - x_1(t) + x_0(t) x_1(t) - c_1 x_1(t) u(t),  \\
  & \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(t) & = & -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], \\

Revision as of 07:01, 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(t)=x0(t)x0(t)x1(t)x0(t)x2(t),x˙1(t)=x1(t)+x0(t)x1(t)c1x1(t)u(t),x˙2(t)=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,40],(c1,c2)=(0.1,0.4),x0=(1.5,0.5,1) or (1.5,1,0.5),α=1.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.