Dynamical study of fractional order differential equations of predator‐pest models

To explore the impact of pest‐control strategy through a fractional derivative, we consider three predator‐prey systems by simple modification of Rosenzweig‐MacArthur model. First, we consider fractional‐order Rosenzweig‐MacArthur model. Allee threshold phenomena into pest population is considered for the second case. Finally, we consider additional food to the predator and harvesting in prey population. The main objective of the present investigation is to observe which model is most suitable for the pest control. To achieve this goal, we perform the local stability analysis of the equilibrium points and observe the basic dynamical properties of all the systems. We observe fractional‐order system has the ability to stabilize Rosenzweig‐MacArthur model with low pest density from oscillatory state. In the numerical simulations, we focus on the bistable regions of the second and third model, and we also observe the effect of the fractional order α throughout the stability region of the system. For the third model, we observe a saddle‐node bifurcation due to the additional food and Allee effect to the pest densities. Also, we numerically plot two parameter bifurcation diagram with respect to the harvesting parameter and fractional order of the system. We finally conclude that fractional‐order Rosenzweig‐MacArthur model and the modified Rosenzweig‐MacArthur model with additional food for the predator and harvested pest population are more suitable models for the pest management.     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Analytical and numerical solutions of a two‐dimensional multi‐term time‐fractional Oldroyd‐B model

In this paper, we consider a two‐dimensional multi‐term time‐fractional Oldroyd‐B equation on a rectangular domain. Its analytical solution is obtained by the method of separation of variables. We employ the finite difference method with a discretization of the Caputo time‐fractional derivative to obtain an implicit difference approximation for the equation. Stability and convergence of the approximation scheme are established in the L_∞ ‐norm. Two examples are given to illustrate the theoretical analysis and analytical solution. The results indicate that the present numerical method is effective for this general two‐dimensional multi‐term time‐fractional Oldroyd‐B model.     

Analytical and approximate solutions of fractional‐order susceptible‐infected‐recovered epidemic model of childhood disease

In this work, we make use of the conformable fractional differential transform method (CFDTM) in order to compute an approximate solution of the fractional‐order susceptible‐infected‐recovered (SIR) epidemic model of childhood disease. The method provides the solution in the form of a rapidly convergent series. Two explanatory and illustrative examples are given to represent the efficacy of the obtained results.     

Synchronization of fractional‐order chaotic systems with disturbances via novel fractional‐integer integral sliding mode control and application to neuron models

In this paper, a novel fractional‐integer integral type sliding mode technique for control and generalized function projective synchronization of different fractional‐order chaotic systems with different dimensions in the presence of disturbances is presented. When the upper bounds of the disturbances are known, a sliding mode control rule is proposed to insure the existence of the sliding motion in finite time. Furthermore, an adaptive sliding mode control is designed when the upper bounds of the disturbances are unknown. The stability analysis of sliding mode surface is given using the Lyapunov stability theory. Finally, the results performed for synchronization of three‐dimensional fractional‐order chaotic Hindmarsh‐Rose (HR) neuron model and two‐dimensional fractional‐order chaotic FitzHugh‐Nagumo (FHN) neuron model.     

Hopf bifurcation of a diffusive Gause‐type predator‐prey model induced by time fractional‐order derivatives

Since population behaviors possess the characteristic of history memory, we, in this paper, introduce time fractional‐order derivatives into a diffusive Gause‐type predator‐prey model, which is time fractional‐order reaction‐diffusion equations and a generalized form of its corresponding first‐derivative model. For this kind of model, we prove the existence and uniqueness of a global positive solution by using the theory of evolution equations and the comparison principle of time fractional‐order partial differential equations. Besides, we obtain the stability and Hopf bifurcation of the Gause‐type predator‐prey model in the forms of the time fractional‐order ordinary equations and of the time fractional‐order reaction‐diffusion equations, respectively. Our results show that the stable region of the parameters in these 2 models can be enlarged by the time fractional‐order derivatives. Some numerical simulations are made to verify our results.     

Application of the modified operational matrices in multiterm variable‐order time‐fractional partial differential equations

In this paper, we present a novel discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variable‐order time‐fractional partial differential equations (M‐V‐TFPDEs) in the large interval. In this purpose, the accurate modified operational matrices are constructed to reduce the problems into a system of algebraic equations. Also, the computational algorithm based on the method and modified operational matrices in the large interval is easily implemented. Furthermore, we discuss the error estimation of the proposed method. Ultimately, to confirm our theoretical analysis and accuracy of numerical approach, several examples are presented.     

Multi‐switching combination–combination synchronization of non‐identical fractional‐order chaotic systems

In this paper, multi‐switching combination–combination synchronization scheme has been investigated between a class of four non‐identical fractional‐order chaotic systems. The fractional‐order Lorenz and Chen's systems are taken as drive systems. The combination–combination of multi drive systems is then synchronized with the combination of fractional‐order Lü and Rössler chaotic systems. In multi‐switching combination–combination synchronization, the state variables of two drive systems synchronize with different state variables of two response systems simultaneously. Based on the stability of fractional‐order chaotic systems, the multi‐switching combination–combination synchronization of four fractional‐order non‐identical systems has been investigated. For the synchronization of four non‐identical fractional‐order chaotic systems, suitable controllers have been designed. Theoretical analysis and numerical results are presented to demonstrate the validity and feasibility of the applied method. Copyright © 2017 John Wiley & Sons, Ltd.     

Chaotic fractional‐order model for muscular blood vessel and its control via fractional control scheme

This article studies the chaotic and complex behavior in a fractional‐order biomathematical model of a muscular blood vessel (MBV). It is shown that the fractional‐order MBV (FOMBV) model exhibits very complex and rich dynamics such as chaos. We show that the corresponding maximal Lyapunov exponent of the FOMBV system is positive which implies the existence of chaos. Strange attractors of the FOMBV model are depicted to validate the chaotic behavior of the system. We change the fractional order of the model and investigate the dynamics of the system. To suppress the chaotic behavior of the model, we propose a single input fractional finite‐time controller and prove its stability using the fractional Lyapunov theory. In addition, the effects of the model uncertainties and external disturbances are taken into account and a robust fractional finite‐time controller is constructed. The upper bound of the chaos suppression time is also given. Some computer simulations are presented to illustrate the findings of this article. © 2014 Wiley Periodicals, Inc. Complexity 20: 37–46, 2014     

An efficient numerical treatment of fourth‐order fractional diffusion‐wave problems

In this paper, we consider the numerical treatment of a fourth‐order fractional diffusion‐wave problem. Our proposed method includes the use of parametric quintic spline in the spatial dimension and the weighted shifted Grünwald‐Letnikov approximation of fractional integral. The solvability, stability, and convergence of the numerical scheme are rigorously proved. It is shown that the theoretical convergence order improves those of earlier work. Simulation is further carried out to demonstrate the numerical efficiency of the proposed scheme and to compare with other methods.     

A study of the fractional‐order mathematical model of diabetes and its resulting complications

Diabetes is a worldwide problem that affects one of every 11 persons nowadays. The IDF Diabetes Atlas (Eighth edition, 2017) states that approximately 415 million people in the world are living with the disease and that this number will rise to 629 million by the year 2045. It is a very serious problem of the world. A major part of the world population is affected by this disease and its resulting complications. In this paper, we propose to investigate a fractional‐order model of diabetes and its resulting complications. The mathematical model's parameters define the population of diabetic patients and those who are diabetic with complications at a given time t. We have also discussed the existence, uniqueness, and stability of the fractional‐order model, which we consider here. We make use of the homotopy decomposition method (HDM) in order to solve the problem.     

A new approach by two‐dimensional wavelets operational matrix method for solving variable‐order fractional partial integro‐differential equations

In this paper, a new computational scheme based on operational matrices (OMs) of two‐dimensional wavelets is proposed for the solution of variable‐order (VO) fractional partial integro‐differential equations (PIDEs). To accomplish this method, first OMs of integration and VO fractional derivative (FD) have been derived using two‐dimensional Legendre wavelets. By implementing two‐dimensional wavelets approximations and the OMs of integration and variable‐order fractional derivative (VO‐FD) along with collocation points, the VO fractional partial PIDEs are reduced into the system of algebraic equations. In addition to this, some useful theorems are discussed to establish the convergence analysis and error estimate of the proposed numerical technique. Furthermore, computational efficiency and applicability are examined through some illustrative examples.     

A gray model for increasing sequences with nonhomogeneous index trends based on fractional‐order accumulation

By comparing the class ratio deviation and restoring error of first‐order accumulation with that of fractional‐order accumulation, a gray model for monotonically increasing sequences can obtain optimal simulation accuracy via selecting a proper cumulative order. In this study, a gray model for increasing sequences with nonhomogeneous index trends based on fractional‐order accumulation is proposed. To reduce the modeling error caused by the background value and to improve the prediction accuracy of the model, an optimized model using the 3/8 Simpson formula is constructed. Finally, the 2 proposed models are used to predict the total energy consumption in China and the monthly sales of new products in an enterprise. Compared with the GM(1,1) model based on fractional‐order accumulation, the proposed model exhibits better simulation and prediction accuracy.     

Integrable couplings of fractional L‐hierarchy and its Hamiltonian structures

Based on a general isospectral problem of fractional order, a fractional bilinear form variational identity, the new integrable coupling of fractional L‐hierarchy and the Hamiltonian structures of the integrable coupling of fractional L‐hierarchy are obtained. Copyright © 2016 John Wiley & Sons, Ltd.     

Fractional‐order orthogonal Bernstein polynomials for numerical solution of nonlinear fractional partial Volterra integro‐differential equations

In this paper, a new two‐dimensional fractional polynomials based on the orthonormal Bernstein polynomials has been introduced to provide an approximate solution of nonlinear fractional partial Volterra integro‐differential equations. For this aim, the fractional‐order orthogonal Bernstein polynomials (FOBPs) are constructed, and its operational matrices of integration, fractional‐order integration, and derivative in the Caputo sense and product operational matrix are derived. These operational matrices are utilized to reduce the under study problem to a nonlinear system of algebraic equations. Using the approximation of FOBPs, the convergence analysis and error estimate associated to the proposed problem have been investigated. Finally, several examples are included to clarify the validity, efficiency, and applicability of the proposed technique via FOBPs approximation.     

Two classes of conformable fractional Sturm‐Liouville problems: Theory and applications

In this paper, we investigate in more detail some useful theorems related to conformable fractional derivative (CFD) and integral and introduce two classes of conformable fractional Sturm‐Liouville problems (CFSLPs): namely, regular and singular CFSLPs. For both classes, we study some of the basic properties of the Sturm‐Liouville theory. In the class of r‐CFSLPs, we discuss two types of CFSLPs which include left‐ and right‐sided CFDs, each of order α∈(n,n+1], and prove properties of the eigenvalues and the eigenfunctions in a certain Hilbert space. Also, we apply a fixed‐point theorem for proving the existence and uniqueness of the eigenfunctions. As an operator for the class of s‐CFSLPs, we first derive two fractional types of the hypergeometric differential equations of order α∈(0,1] and obtain their analytical eigensolutions as Gauss hypergeometric functions. Afterwards, we define the conformable fractional Legendre polynomial/functions (CFLP/Fs) as Jacobi polynomial and investigate their basic properties. Moreover, the conformable fractional integral Legendre transforms (CFILTs) based on CFLP/Fs‐I and ‐II are introduced, and using these new transforms, an effective procedure for solving explicitly certain ordinary and partial conformable fractional differential equations (CFDEs) are given. Finally, as a theoretical application, some fractional diffusion equations are solved.     

Optimal synchronization of two different in‐commensurate fractional‐order chaotic systems with fractional cost function

This article investigates the optimal synchronization of two different fractional‐order chaotic systems with two kinds of cost function. We use calculus of variations for minimizing cost function subject to synchronization error dynamics. We introduce optimal control problem to solve fractional Euler–Lagrange equations. Optimal control signal and minimum time of synchronization are obtained by proposed method. Examples show the optimal synchronization of two different systems with two different cost functions. First, we use an ordinary integer cost function then we use a fractional‐order cost function and comparing the results. Finally, we suggest a cost function which has the optimal solution of this problem, and we can extend this solution to solve other synchronization problems. © 2016 Wiley Periodicals, Inc. Complexity 21: 401–416, 2016     

Application of the Caputo‐Fabrizio and Atangana‐Baleanu fractional derivatives to mathematical model of cancer chemotherapy effect

In this paper, we obtain approximate‐analytical solutions of a cancer chemotherapy effect model involving fractional derivatives with exponential kernel and with general Mittag‐Leffler function. Laplace homotopy perturbation method and the modified homotopy analysis transform method were applied. The first method is based on a combination of the Laplace transform and homotopy methods, while the second method is an analytical technique based on homotopy polynomial. The cancer chemotherapy effect equations are solved numerically and analytically using the aforesaid methods. Illustrative examples are included to demonstrate the validity and applicability of the presented technique with new fractional‐order derivatives with exponential decay law and with general Mittag‐Leffler law.     

High‐order compact schemes for fractional differential equations with mixed derivatives

In this article, we consider two‐dimensional fractional subdiffusion equations with mixed derivatives. A high‐order compact scheme is proposed to solve the problem. We establish a sufficient condition and show that the scheme converges with fourth order in space and second order in time under this condition.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2141–2158, 2017