Research Article  Open Access
Optimal Control of a Delayed HIV Infection Model via Fourier Series
Abstract
We present a delayed optimal control which describes the interaction of the immune system with the human immunodeficiency virus (HIV) and CD4^{+} Tcells. In order to improve the therapies, treatment and the intracellular delays are incorporated into the model. The optimal control in this model represents the efficiency of drug treatment in preventing viral production and new infections. The optimal pair of control and trajectories of this nonlinear delay system with quadratic cost functional is obtained by Fourier series approximation. The method is based on expanding time varying functions in the nonlinear delay system into their Fourier series with unknown coefficients. Using operational matrices for integration, product, and delay, the problem is reduced to a set of nonlinear algebraic equations.
1. Introduction
Delays occur frequently in biological, chemical, electronic, and transportation systems [1]. Many mathematical models have been developed in order to understand the dynamics of HIV infection [2–7]. Moreover, optimal control methods have been applied to the derivation of optimal therapies for this HIV infection [8–14]. All these methods are based on HIV models which ignore the intracellular delay by assuming that the infectious process is instantaneous; that is, as soon as the virus enters an uninfected cell, it starts to produce virus particles; however, this is not reasonable biologically. In this paper, we consider the mathematical model of HIV infection with intracellular delay presented in [15] in order to make the model more tangible and closer to what happens in reality.
Orthogonal functions (OFs) have received considerable attention in dealing with various problems of dynamic systems. Using operational matrices, the technique is based on reduction of these problems to systems of algebraic equations. Special attention has been given to applications of Walsh functions [16], blockpulse functions [17], Laguerre polynomials [18], Legendre polynomials [19], Chebyshev polynomials [20], and Fourier series [21].
In this paper, we apply Fourier series approximation to find the optimal pair of control and trajectories of the nonlinear delayed optimal control system governed by ordinary delay differential equations which describe the interaction of the human immunodeficiency virus (HIV). Operational matrices of integration, product, and delay have the most important role in our method. The paper is organized as follows. Section 2 consists of an introduction to Fourier series approximation and operational and other matrices, being used in Section 4. Section 3 proposes the delayed optimal control model of HIV infection. In Section 4, we utilize the Fourier series approximation to solve our model and results are demonstrated by some figures. Finally, the conclusions are summarized in Section 5.
2. Properties of Fourier Series
A measurable function defined over the interval to may be expanded into a Fourier series as follows: where the Fourier coefficients and are given by The series in (1) has an infinite number of terms. To obtain an approximate expression for , one can truncate the series up to the th term as follows: where the Fourier series coefficient vector and the Fourier series vector are defined as where
The elements of are orthogonal in the interval . By integrating both sides of (6a) and (6b) with respect to , we obtain
Now, we approximate the function in the (7) by a truncated Fourier series. Consequently, the forward integral of the Fourier series vector can be represented by where is the Fourier series operational matrix of forward integration and is given as As a result, we obtain Moreover, where The delay function is the shift of the function defined in (5) along the time axis by . The general expression is given by where is the delay operational matrix of Fourier series, which is as follows:
Also, we have (see [21]) Using (3) and (10) leads to We can get where is the product operational matrix for the vector (see [22]).
As a result,
3. HIV Delayed Optimal Control Problem
Using time delay dynamic systems, many physical systems are best modeled as follows: where the state is a vector function, is an vector control function, and (, ) are nonnegative constant time delays, and the vector functions and are defined appropriately and given (see [23–26]).
In some systems of this type, it is desirable to select the optimal pair which satisfies (22a), (22b), and (22c) and minimizes performance criterion modeled by a cost function of the form
The following model is a delayed control system of HIV infection of CD Tcells (for more details, see [15]): where control function is showing the effect of drugs on HIV virus production. For the delayed control model (24a), (24b), (24c), (24d), (24e), (24f), (24g), and (24h), we consider the objective functional to be defined as where , the weight of current costs of treatment, is assumed to be . Our goal is to maximize the objective functional (25) subject to the delayed control system (24a), (24b), (24c), (24d), (24e), (24f), (24g), and (24h), that is, to maximize the total count of Tcells and to minimize the costs of treatment. The treatment interval is assumed to be , which shows the duration on treatment in terms of days.
4. Solving the Delayed Optimal Control Model of the HIV Infection Using Fourier Series
According to Section 2, each function of , and in the interval can be approximated by a Fourier series. So, where , and are, respectively, the Fourier series coefficient vectors (all are unknown) of , and . For time delay functions and , we have
Therefore, for the interval , we have
By using (17), (18), and (26), one can rewrite (25) as follows: where is Fourier series of , while .
By integrating both sides of (24a) from to and using (26), (20), and (10), we obtain where and are the Fourier series coefficient vectors of and , respectively; that is, , .
So, Similarly, for (24b), we have by using (28a) and (28b) for the righthand side of (24b), one can obtain where , , and .
To compute , we consider , which is approximated by Fourier series as where is the Fourier series coefficient vector of in the interval and .
By using (13) and (34), we have
To find , one can use (12), (21), and (15). Hence, Similarly, Finally, by considering , and , we have Now, we equalize (32) and (38); from the two sides of (24b),
Similarly, by integrating (24c) from to , we obtain
Eliminating from (31), (39), and (40) results in the following equations:
To apply the constraint , in (24h), one can first discretize the time interval into grid points as follows: where . Thus, can be rewritten as Similarly, for , and , we have Finally, the optimal control problem now is reduced to
Now, one can use various toolboxes like Matlab software to solve the above problem. The parameters of viral spread in model (24a), (24b), (24c), (24d), (24e), (24f), (24g), and (24h) are presented in Table 1.

The results of this optimization problem for and are depicted in Figure 1. As can be seen, utilizing the treatment increases the Tcells (Figure 1(a)) and decreases the infected Tcells (Figure 1(b)) and viral load (Figure 1(c)).
(a) Tcells ()
(b) Infected Tcells ()
(c) Viral load ()
(d) Control ()
5. Conclusion
In this paper, we consider a delayed mathematical model describing HIV infection of Tcells during therapy.
Our objective is to minimize the cost of treatment and to maximize the uninfected Tcells. We use Fourier series approximation in order to solve the socalled mathematical model. By applying operational matrices for integration, product, and delay, the nonlinear optimal control model is reduced to a set of nonlinear algebraic equations. Since the set of sines and cosines is orthogonal, the operational matrices contain many zero elements. Hence, it makes the method computationally easy and attractive. The results show the efficiency of the method.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research was in part supported by a grant from the Center of Excellence on Modelling and Control Systems (CEMCS).
References
 M. Jamshidi and C. M. Wang, “A computational algorithm for largescale nonlinear timedelay systems,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 14, no. 1, pp. 2–9, 1984. View at: Publisher Site  Google Scholar  MathSciNet
 K. Hattaf and N. Yousfi, “Dynamics of HIV infection model with therapy and cure rate,” International Journal of Tomography and Statistics, vol. 16, no. 11, pp. 74–80, 2011. View at: Google Scholar
 K. Hattaf, N. Yousfi, and A. Tridane, “Mathematical analysis of a virus dynamics model with general incidence rate and cure rate,” Nonlinear Analysis: Real World Applications, vol. 13, no. 4, pp. 1866–1872, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. S. Perelson and P. W. Nelson, “Mathematical analysis of HIV1 dynamics in vivo,” SIAM Review, vol. 41, no. 1, pp. 3–44, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 A. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, and D. D. Ho, “HIV1 dynamics in vivo: virion clearance rate, infected cell lifespan, and viral generation time,” Science, vol. 271, no. 5255, pp. 1582–1586, 1996. View at: Publisher Site  Google Scholar
 X. Zhou, X. Song, and X. Shi, “A differential equation model of HIV infection of CD4+ Tcells with cure rate,” Journal of Mathematical Analysis and Applications, vol. 342, no. 2, pp. 1342–1355, 2008. View at: Google Scholar
 H. Zhu and X. Zou, “Dynamics of a HIV1 infection model with cellmediated immune response and intracellular delay,” Discrete and Continuous Dynamical Systems B, vol. 12, no. 2, pp. 511–524, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 K. Hattaf and N. Yousfi, “Two optimal treatments of HIV infection model,” World Journal of Modelling and Simulation, vol. 8, no. 1, pp. 27–36, 2012. View at: Google Scholar
 K. R. Fister, S. Lenhart, and J. S. McNally, “Optimizing chemotherapy in an HIV model,” Electronic Journal of Differential Equations, vol. 1998, 12 pages, 1998. View at: Google Scholar  MathSciNet
 H. R. Joshi, “Optimal control of an HIV immunology model,” Optimal Control Applications & Methods, vol. 23, no. 4, pp. 199–213, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 J. Karrakchou, M. Rachik, and S. Gourari, “Optimal control and infectiology: application to an HIV/AIDS model,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 807–818, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 D. Kirschner, S. Lenhart, and S. Serbin, “Optimal control of the chemotherapy of HIV,” Journal of Mathematical Biology, vol. 35, no. 7, pp. 775–792, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 H.D. Kwon, “Optimal treatment strategies derived from a HIV model with drugresistant mutants,” Applied Mathematics and Computation, vol. 188, no. 2, pp. 1193–1204, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 R. Culshaw, S. Ruan, and R. J. Spiteri, “Optimal HIV treatment by maximising immune response,” Journal of Mathematical Biology, vol. 48, no. 5, pp. 545–562, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 R. V. Culshaw and S. Ruan, “A delaydifferential equation model of HIV infection of CD4^{+} Tcells,” Mathematical Biosciences, vol. 165, no. 1, pp. 27–39, 2000. View at: Publisher Site  Google Scholar
 W. Chen and Y. Shih, “Shift Walsh matrix and delaydifferential equations,” IEEE Transactions on Automatic Control, vol. 23, no. 6, pp. 1023–1028, 1978. View at: Publisher Site  Google Scholar
 G. P. Rao and T. Srinivasan, “Analysis and synthesis of dynamic systems containing time delays via blockpulse functions,” Proceedings of the Institution of Electrical Engineers, vol. 125, no. 10, pp. 1064–1068, 1978. View at: Google Scholar  MathSciNet
 F. Kung and H. Lee, “Solution and parameter estimation in linear tim einvariant delay systems using Laguerre polynomial expansion,” Journal of Dynamic Systems, Measurement and Control, vol. 105, no. 4, pp. 297–301, 1983. View at: Publisher Site  Google Scholar
 L. Lee and F. Kung, “Shifted Legendre series solution and paramet er estimation of linear delayed systems,” International Journal of Systems Science, vol. 16, no. 10, pp. 1249–1256, 1985. View at: Publisher Site  Google Scholar
 I. R. Horng and J. H. Chou, “Analysis, parameter estimation and optimal control of timedelay systems via Chebyshev series,” International Journal of Control, vol. 41, no. 5, pp. 1221–1234, 1985. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Razzaghi and M. Razzaghi, “Fourier series direct method for variational problems,” International Journal of Control, vol. 48, no. 3, pp. 887–895, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Gh. Ghanbari and M. H. Farahi, “Optimal control strategy for a HIV infection model via Fourier series,” Journal of Applied Mathematics, Statistics and Informatics, vol. 9, no. 2, pp. 87–94, 2014. View at: Publisher Site  Google Scholar
 D. H. Eller, J. K. Aggarwal, and H. T. Banks, “Optimal control of linear timedelay systems,” IEEE Transactions on Automatic Control, vol. 14, pp. 678–687, 1969. View at: Google Scholar  MathSciNet
 N. N. Krasovskii, Stability of Motion: applications of Lyapunov's second Method to Differential Systems and Equations with Delay, 1963. View at: MathSciNet
 R. Loxton, K. L. Teo, and V. Rehbock, “An optimization approach to statedelay identification,” IEEE Transactions on Automatic Control, vol. 55, no. 9, pp. 2113–2119, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 L. Göllmann, D. Kern, and H. Maurer, “Optimal control problems with delays in state and control variables subject to mixed controlstate constraints,” Optimal Control Applications and Methods, vol. 30, no. 4, pp. 341–365, 2009. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2014 Gh. Ghanbari and M. H. Farahi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.