Numerical Optimal Control of HIV Transmission in Octave/MATLAB

by   Carlos Campos, et al.

We provide easy and readable GNU Octave/MATLAB code for the simulation of mathematical models described by ordinary differential equations and for the solution of optimal control problems through Pontryagin's maximum principle. For that, we consider a normalized HIV/AIDS transmission dynamics model based on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde. Ecol. Complex. 2017, 30, 70–75), given by a system of four ordinary differential equations. An HIV initial value problem is solved numerically using the ode45 GNU Octave function and three standard methods implemented by us in Octave/MATLAB: Euler method and second-order and fourth-order Runge-Kutta methods. Afterwards, a control function is introduced into the normalized HIV model and an optimal control problem is formulated, where the goal is to find the optimal HIV prevention strategy that maximizes the fraction of uninfected HIV individuals with the least HIV new infections and cost associated with the control measures. The optimal control problem is characterized analytically using the Pontryagin Maximum Principle, and the extremals are computed numerically by implementing a forward-backward fourth-order Runge-Kutta method. Complete algorithms, for both uncontrolled initial value and optimal control problems, developed under the free GNU Octave software and compatible with MATLAB are provided along the article.


page 1

page 2

page 3

page 4


Numerical solution of optimal control of atherosclerosis using direct and indirect methods with shooting/collocation approach

We present a direct numerical method for the solution of an optimal cont...

Review: Ordinary Differential Equations For Deep Learning

To better understand and improve the behavior of neural networks, a rece...

Generalization of partitioned Runge–Kutta methods for adjoint systems

This study computes the gradient of a function of numerical solutions of...

Differential Flatness as a Sufficient Condition to Generate Optimal Trajectories in Real Time

As robotic systems increase in autonomy, there is a strong need to plan ...

Extraction of a computer-certified ODE solver

Reliably determining system trajectories is essential in many analysis a...

Neural Optimal Control using Learned System Dynamics

We study the problem of generating control laws for systems with unknown...

Optimal Control of a Differentially Flat 2D Spring-Loaded Inverted Pendulum Model

This paper considers the optimal control problem of an extended spring-l...

Please sign up or login with your details

Forgot password? Click here to reset