Numerical Fractional Optimal Control of Respiratory Syncytial Virus Infection in Octave/MATLAB

Abstract

In this article, we develop a simple mathematical GNU Octave/MATLAB code that is easy to modify for the simulation of mathematical models governed by fractional-order differential equations, and for the resolution of fractional-order optimal control problems through Pontryagin’s maximum principle (indirect approach to optimal control). For this purpose, a fractional-order model for the respiratory syncytial virus (RSV) infection is considered. The model is an improvement of one first proposed by the authors in 2018. The initial value problem associated with the RSV infection fractional model is numerically solved using Garrapa’s fde12 solver and two simple methods coded here in Octave/MATLAB: the fractional forward Euler’s method and the predict-evaluate-correct-evaluate (PECE) method of Adams–Bashforth–Moulton. A fractional optimal control problem is then formulated having treatment as the control. The fractional Pontryagin maximum principle is used to characterize the fractional optimal control and the extremals of the problem are determined numerically through the implementation of the forward-backward PECE method. The implemented algorithms are available on GitHub and, at the end of the paper, in appendixes, both for the uncontrolled initial value problem as well as for the fractional optimal control problem, using the free GNU Octave computing software and assuring compatibility with MATLAB.

Optimal Control for Fractional Order Dynamics of Tumor Growth

“…A modified Adams-Bashforth-Moulton method's Predictor-Evaluate-Correct-Evaluate (PECE) will be used to solve our FOCPs, as shown in (46 -49) following [8]. We start by solving equation (53) forward in time given initial conditions and the control, and then we use the Adams-type predictor-corrector approach, see also [17,26,8]. For a general problem…”

Fractional Optimal Control Model and Bifurcation Analysis of Human Syncytial Respiratory Virus Transmission Dynamics

“…In this study, we modified the research conducted in [34]; the total population in the study is classified as follows. (i) Susceptible (S): This group consists of population free of virus, but with chance to become infected with possible exposure to the virus in the future.…”