Numerical Optimal Control of HIV Transmission in Octave/MATLAB

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.

[1]  Delfim F. M. Torres,et al.  Modeling and optimal control of HIV/AIDS prevention through PrEP , 2017, 1703.06446.

[2]  John W. Eaton,et al.  GNU octave : a high-level interactive language for numerical computations : Octave version 2.0.17 (stable) , 1997 .

[3]  Radu Burlacu On a class of optimal control problems in Mathematical Biology , 1999 .

[4]  Delfim F. M. Torres,et al.  Direct transcription methods based on fractional integral approximation formulas for solving nonlinear fractional optimal control problems , 2018, Commun. Nonlinear Sci. Numer. Simul..

[5]  Gary An,et al.  The Crisis of Reproducibility, the Denominator Problem and the Scientific Role of Multi-scale Modeling , 2018, Bulletin of Mathematical Biology.

[6]  Singiresu S. Rao,et al.  Optimization Theory and Applications , 1980, IEEE Transactions on Systems, Man, and Cybernetics.

[7]  Delfim F. M. Torres,et al.  Mathematical Modelling, Simulation, and Optimal Control of the 2014 Ebola Outbreak in West Africa , 2015, 1503.07396.

[8]  Suzanne Lenhart,et al.  Optimal control of treatments in a two-strain tuberculosis model , 2002 .

[9]  Delfim F. M. Torres,et al.  Ebola model and optimal control with vaccination constraints , 2017, 1703.01368.

[10]  John T. Workman,et al.  Optimal Control Applied to Biological Models , 2007 .

[11]  F. Brauer,et al.  Mathematical Models in Population Biology and Epidemiology , 2001 .

[12]  Delfim F. M. Torres,et al.  Optimal Control and Sensitivity Analysis of a Fractional Order TB Model , 2018, Statistics, Optimization & Information Computing.

[13]  Delfim F. M. Torres,et al.  Optimal control for a tuberculosis model with reinfection and post-exposure interventions. , 2013, Mathematical biosciences.

[14]  Delfim F. M. Torres,et al.  A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde , 2016, 1612.00732.

[15]  L. S. Pontryagin,et al.  Mathematical Theory of Optimal Processes , 1962 .

[16]  Delfim F. M. Torres,et al.  Enhancement of chemotherapy using oncolytic virotherapy: Mathematical and optimal control analysis. , 2018, Mathematical biosciences and engineering : MBE.

[17]  Delfim F. M. Torres,et al.  Analysis of a SIRI Epidemic Model with Distributed Delay and Relapse , 2018, Statistics, Optimization & Information Computing.

[18]  John W. Eaton,et al.  GNU Octave manual version 3: a high-level interactive language for numerical computations , 2008 .

[19]  Somayeh Nemati,et al.  A numerical approach for solving fractional optimal control problems using modified hat functions , 2019, Commun. Nonlinear Sci. Numer. Simul..

[20]  Sunil Deshpande,et al.  Optimal Input Signal Design for Data-Centric Identification and Control with Applications to Behavioral Health and Medicine , 2014 .

[21]  E. Blum,et al.  The Mathematical Theory of Optimal Processes. , 1963 .

[22]  Delfim F. M. Torres,et al.  Analysis and Optimal Control of an Intracellular Delayed HIV Model with CTL Immune Response , 2018, Math. Comput. Sci..