DeepAI AI Chat
Log In Sign Up

Numerical Optimal Control of HIV Transmission in Octave/MATLAB

12/19/2019
by   Carlos Campos, et al.
Universidade de Aveiro
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

11/01/2019

Review: Ordinary Differential Equations For Deep Learning

To better understand and improve the behavior of neural networks, a rece...
03/22/2020

Generalization of partitioned Runge–Kutta methods for adjoint systems

This study computes the gradient of a function of numerical solutions of...
03/04/2021

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 ...
04/06/2021

Extraction of a computer-certified ODE solver

Reliably determining system trajectories is essential in many analysis a...
02/20/2023

Neural Optimal Control using Learned System Dynamics

We study the problem of generating control laws for systems with unknown...
11/17/2019

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

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