DOI:10.32604/cmc.2021.015161
| Computers, Materials & Continua DOI:10.32604/cmc.2021.015161 | |
| Article |
Optimal and Memristor-Based Control of A Nonlinear Fractional Tumor-Immune Model
1Department of Mathematics, College of Science, Taif University, Taif, 21944, Saudi Arabia
2Department of Mathematics, Faculty of Science, Zagazig University, Zagazig, 44519, Egypt
3Department of Physics and Engineering Mathematics, Faculty of Electronic engineering, Menoufia University, Menouf, 32952, Egypt
4Department of Mathematics, Faculty of Science, Al-Azher University, Nasr City, 11884, Egypt
*Corresponding Author: Amr M. S. Mahdy. Email: amattaya@tu.edu.sa
Received: 08 November 2020; Accepted: 05 January 2021
Abstract: In this article, the reduced differential transform method is introduced to solve the nonlinear fractional model of Tumor-Immune. The fractional derivatives are described in the Caputo sense. The solutions derived using this method are easy and very accurate. The model is given by its signal flow diagram. Moreover, a simulation of the system by the Simulink of MATLAB is given. The disease-free equilibrium and stability of the equilibrium point are calculated. Formulation of a fractional optimal control for the cancer model is calculated. In addition, to control the system, we propose a novel modification of its model. This modification is based on converting the model to a memristive one, which is a first time in the literature that such idea is used to control this type of diseases. Also, we study the system’s stability via the Lyapunov exponents and Poincare maps before and after control. Fractional order differential equations (FDEs) are commonly utilized to model systems that have memory, and exist in several physical phenomena, models in thermoelasticity field, and biological paradigms. FDEs have been utilized to model the realistic biphasic decline manner of elastic systems and infection of diseases with a slower rate of change. FDEs are more useful than integer-order in modeling sophisticated models that contain physical phenomena.
Keywords: RDTM; tumor-immune; optimal control; caputo derivative; signal flow; simulink; disease-free equilibrium; stability; memristive; lyapunov exponents; poincare map
Delayed ordinary differential equations have been utilized in modeling cancer diseases [1–10]. Fractional order differential equations (FDEs) are commonly utilized to model systems that have memory that exists in several physical phenomena, models in thermoelasticity field and biological paradigms. In [1], a system of FDEs was applied for modeling the interactions in the cancer-immune system. The model comprises double immune effectors: 
(as an example for natural killer cells and cytotoxic T cells), interactive ?? against the cancer cells, 
, with a type III’s Holling function. A type III’s Holling function represents a case where the number of victims wasted per each predator at start increases slowly as the density of victim increases, but then levels decrease with further increase in victim density. That means the response of predators to the victim is decreased at low victim density, then levels decrease with further increase in victim density.
The form of the model is (see [11]):
where 
is Caputo fractional derivative operator ( [12–15]) with 
, 
, 
and 
are positive constants.
Using a signal flow graph for representing the dynamical systems is very useful (see for example [16,17]).
The signal flow graph is a graph tool that can be used to show the interrelation between the system states and enable us to use the graph theoretic tools to discover new features of the system.
The above system Eq. (1) can be represented by the signal flow graph 
as shown in Fig. 1. The signal flow diagram of the calculated system Eq. (1) has the following adjacency matrix 
. Where each state variable is modeled by a vertex and there is an edge about two states (a, b) if the state a affects directly the state b according to Eq. (1).
Figure 1: The signal flow graph of the studied system (Eq. (1))
For more details about signal flow graph and theoretical graph theory, see for example ([16,18–20]) and the references therein.
In this work, the reduced differential transform method is introduced for solving the nonlinear fractional model of Tumor-Immune. The fractional derivatives are described in the Caputo sense. The solutions derived using this method are easy and very accurate. The model is given by its signal flow diagram. Moreover, a simulation of the system using the Simulink of MATLAB is given. The disease-free equilibrium and stability of the equilibrium point are calculated. Formulation of a fractional optimal control for the cancer model is calculated. In addition, To control the system, we proposed a novel modification of this model utilizing the method of converting it to be a memristive system which is the first time in the literature to use such idea to control this type of diseases. Also, we study the system’s stability via the Lyapunov exponent and Poincare map before and after control. FDEs have been utilized to model the realistic biphasic decline manner of elastic systems and infection of diseases but at a slower rate of change. FDEs are more useful than integer-order in modeling sophisticated models that contain physical phenomena.
The paper is structured in eight sections. In Section 2, the basic definitions of fractional calculus are presented. In Section 3, the formulation of a fractional optimal control of the cancer model is studied. In Section 4, we describe the Reduced differential transform method (RDTM) with illustration examples. In Section 5, fixed points and the stability of the system are investigated. In Section 6, the Cancer model is simulated using Simulink
Matlab and campared with the RDTM. In Section 7, we propose a novel method to control the system based on the idea of the memristor. Also, the equilibrium point and stability of nonlinear fractional memristor-based cancer model are studied. In addition, we study the system’s stability via the Lyapunov exponent and Poincare map before and after control. Conclusions are given in Section 8.
2 Basic Definitions of Fractional Calculus
In this section, we present the basic definitions and properties of the fractional calculus theory, which are used in this paper.
Definition 1 Areal function 
, z > 0, is said to be in the space 
, 
, if there exists a real number 
such that 
where 
and it is said to be in the space 
if 
.
Definition 2 The Riemann–Liouville integral operator of order 
with 
is defined as ([21–29]):
Definition 3 The Caputo fractional derivative operator 
of order 
is defined in the following form ([11–14,23–30]):
3 Formulation of Fractional Optimal Control of Cancer Model
Consider the state system given in Eq. (1), in R3, with the set of admissible control functions:
where Tfis the final time, 
are the control functions.
The objective function is defined as follows (quadratic is the control variable).
Wherever both immune effectors coexist, the non biological inner solution is measured by A, B and C.
Minimizing the following objective function is the main aim in FOCPs by finding the optimal controls 
:
subjected to the constraints
The following inital conditions are satisfied:
In order to give a definition of the FOCP, define a modified objective (cost) function as:
where we can define the Hamiltonian of the objective function (8) and the cancer immune system (1) as follows:
From (8) and (10), the necessary and sufficient conditions of FOPC can be derived as follows.
where 
are the Lagrange multipliers. Eqs. (12) and (13) produce the necessary conditions in terms of a Hamiltonian of the FOPC. We arrive at the following theorem.
Theorem 1. 
are optimal controls with corresponding state 
then there are adjoint elements, 
satisfy the following conditions.
(i) Co-state equations (adjoint equations)
(ii) Transversal cases
(iii) Optimality conditions
Moreover, the controlling functions 
are presented by
For more details about problem optimal control, see for example [31–34].
4 Applications of Reduced Differential Transform Method (RDTM)
To clarify the efficiency of our proposed method [34–41], we shall apply it to a special case of fractional-order biological systems that presented in [11]. The symbolic calculus software MatLab is used to calculate all of the results given here. Recently, there is some growth in the area of numerical study as well as their applications ([42–55]).
Example 4.1 Take into account the cancer fractional model (1) (c.f [11]):
Applying RDTM Eq. (22), we get:
By substituting Eq. (1) in Eq. (23).
The series solution of the system (1.1) can be calculated via applying the differential inverse transform, where a = r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and different 
and get it as:
We get the solution as a series (for illustration see Figs. 2–4):
Figure 2: The numerical simulations T, E1 and E2 of model cancer at 
Figure 3: The numerical simulations of T, E1 and E2 model cancer at 
Figure 4: The numerical simulations T, E1 and E2 of model cancer at 
From model (1) the numerical results are given in Figs. 2–4. Where in Figs. 2–4, we present the behavior of numerical simulations of T, E1 and E2 model cancer with different values of 
.
5 Fixed Point and the Stability of The System
To compute the fixed points, equate all right hand sides of (1) to zero.
We calculate the Jacobian matrix as:
Stability of 
:
To calculate the Eigen-values, we write:
Then we have: 
.
Stability of all fixed points using Lyapunov exponents and poincare map satisfy the behavior of the uncontrolled system, see Figs. 5–8.
Figure 5: The dependency of Lyapunov exponents on the system parameters
Figure 6: The behavior of Lyapunov exponenets vs. time of the uncontrolled system where 
and 
Figure 7: Poincere map of 
vs. 
of the uncontrolled system where 
and 
Figure 8: Poincere map of 
vs. 
of the uncontrolled system where 
d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and 
In this section, in Fig. 9, we ready a simulation of the system by Simulink of MATLAB. From which, as shown in Figs. 10–12, it is clear that the solutions of the system by the proposed method are the same as it from the Simulink. In addition, the diagram of the Simulink shows the dependency of the system components on each other.
Figure 9: System simulation by MATLAB
Simulink
Figure 10: T(t) of cancer model at 
Figure 11: E1(t) of cancer model at 
Figure 12: E2(t) of cancer model at 
The following Figs. 10–12, show the simulation responses that completely agree with the analytical solution. Figs. 13–15, represents the phase spaces.
Figure 13: E1(t) vs. T(t) space from Simulink simulation of the cancer model
Figure 14: E2(t) vs. T(t) space from Simulink simulation of the cancer model
Figure 15: E2(t) vs. E1(t) space from Simulink simulation of the cancer model
7 The Memristor-Based Conrelled Cancer Model
The memristive system is the system that has a memristor. The memristor is a variable that store the history of a selected state variable. Converting systems to be memristive means using the history of the system to improve and control the dynamics of the system.
Considering the symmetry of the cancer model (1), we modified this model to become a memristive system via the method of adding a memristor effect to the first equation. In this work, we use the memristor that has a quadratic nonlinearity for flux 
and the electric charge q .
So, the memductance is given by 
where c1 and c2 are two memristor parameters with positive values (for more details about types of memristors and its applications see for example [35]). We suggest to add a memristor effect as the feedback term on the first equation of the original system T(t) and use the second equation as the internal state equation of memristor, then a 3-dimensional memristive controlled cancer system can be constructed:
In the following, we study the fixed point and the stability of the proposed memristive controlled cancer system.
I Fixed Points of the Proposed Memristive Controlled Cancer System
To compute the fixed points solve the following system:
II Studying the Stability:
We calculate the Jacobian matrix as: 
Stability of E0 : 
That has the following Eigen-values:
The above system has solved numerically using Rung–Kuta method where a = r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7, c1 = 0.0005, c2 = 0.005 and the initial conditions are 
The solution of the memristive cancer model is shown in Figs. 16–18. The proposed memristive cancer model is well controlled as shown in its time response see Figs. 19–21. The stability of the equilibrium points are shown by the behavior of Lyapunov exponents and Poincere maps see Figs. 22–27.
Figure 16: T(t) of the proposed memristive cancer model is well controlled
Figure 17: E1(t) of the proposed memristive cancer model is well controlled
Figure 18: E2(t) of the proposed memristive cancer model is well controlled
Figure 19: The phase plane E1(t) vs. T(t) of the memristive controlled cancer model
The phase planes of the memristive canser model are shown in Figs. 19–21.
Figure 20: The phase plane E2(t) vs. T(t) of the memristive controlled cancer model
Figure 21: The phase plane E2(t) vs. E1(t) of the memristive controlled cancer model
Figure 22: The behavior of Lyapunov exponenets vs. time of the controlled system where 
Figure 23: Poincere map of 
vs. 
of the controlled system where 
d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7
Figure 24: Poincere map of 
vs. 
of the controlled system where a = 3 , r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and 
Figure 25: Poincere map of 
vs. 
of the controlled system where a = 2 , r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and 
Figure 26: Poincere map of 
vs. 
of the controlled system where a = 1 , r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and 
Figure 27: Poincere map of 
vs. 
of the controlled system where a = 3 , r1 = r2 = 1, d1 = 0.3, d2 = 0.7, k1 = 0.3, k2 = 0.7 and 
In this paper, a nonlinear cancer fractional model of the Tumor-Immune problem has been calculated which plays a necessary role in applied sciences. The fractional derivatives have been described in the Caputo sense. Also, RDTM has been applied to get the approximate solutions of this model. A signal flow diagram of the calculated system has been proposed and disscussed. The free disease equilibrium (FDE) and stability of equilibrium point have been studied. A simulation of the system using Simulink of MATLAB has been presesnted. The phase space has been displayed. To control the nonlinear fractional model of the Tumor-Immune, we have proposed a novel modification of this model via converting it to be a memristive system. It is the first time in the literature to convert such models to be memristive. Also, we have studied the system’s stability via Lyapunov exponents and Poincare maps before and after control. The numerical simulations are very consistent with the analytical ones.
Acknowledgement: The authors are thankful of the Taif University. Taif University researchers supporting project number (TURSP-2020/160), Taif University, Taif, Saudi Arabia.
Funding Statement: This paper was funded by “Taif University Researchers Supporting Project number (TURSP-2020/160), Taif University, Taif, Saudi Arabia”.
Conflicts of interest: The authors declare that they have no conflicts of interest to report regarding the present study.
 | This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. |