1. Introduction
The fractional differential equation (FDE) and the fractional calculus (FC) are now the most growing research areas and receiving good feedback from scientists and authors working in modeling real-life problems in the scientific and engineering fields. It is well-known that the fractional differential equation has many applications in scientific and engineering sectors. In particular, the dynamical systems and the chaos theory are also getting much attention from the researchers due to their wide applications in many fields such as autonomous mobile robots, chaotic communication, image watermarking, and circuits, etcetera [1-3]. Many researchers studied the chaotic and other models in fractional operators see for example [4-15]. These operators in power law, exponential decay, and Mittag-Leffler kernel were utilized successfully to present the various models of practical purpose. For example, the authors considered in Ref. [4] the chaotic financial model and its control strategy. The time delay chaotic model for the circuit and its synchronization has been investigated in Ref. [5]. The fractal-fractional (FF) operators introduced recently has been getting an impressive response from researchers. Using the idea of FF and its application to the existing chaotic models in the literature are studied in Ref. [6]. The new chaotic system in fractional derivative and its comparative study is suggested in Ref. [7], while the application of fractional calculus to African tortoise and image processing has been discussed in Ref. [8,9]. The application of fractional derivative to zoonotic disease [10], application to dengue epidemic, see [11], tuberculosis, see [12-14], and the Chua model [15] has been investigated in arbitrary derivative.
Mostly the problems formulated in scientific and engineering areas are of high nonlinear nature, and its exact solution is not possible. The solution of a fractional differentiable equation is also very complicated due to its order and difficult to solve analytically, and thus the researchers always looking to have their solution numerically. Fractional differential equation and its application to nonlinear problems of different fields with novel numerical investigations have been studied in Refs. [10,16,17] and the references therein. Nowadays the researchers from different areas developing mathematical models in fractional calculus due to its applications. One of the interesting applications is the best fitting to the real cases for an epidemic disease and another real data setting that come from different areas of engineering and sciences. The memory and the heredity properties are attached to the fractional operators and can best determine the dynamics of dynamical systems and other related areas problems efficiently than the classical orders. It is the advantage of fractional calculus that can determine well the analysis between any two suggested points. Some recent work related to fractional operators and their applications can be seen in Refs. [18-20]. In this work, the authors focused on studying the circuit system, fractional-order capacitor, and sound encryption system using arbitrary order derivatives. The fractal order derivative and the fractional-order, as well as their combinations for the solution of physical problems have been suggested in Refs. [21-26]. The chaotic Shinriki’s oscillator fractional model has been investigated by the authors in Ref. [21]. In Ref. [22], the authors studied the periodicity and multiple attractors for Villis model. The idea of FF has been considered to study influenza and malaria model in Refs. [23,24]. The fractal-fractional operators have been applied effectively for the related problems, see [25-27], where the authors obtain results for new chaotic attractors, the study of the Hepatitis C model, and for system identification. Some more interesting work, where the authors proposed the fractional derivatives as an application to the scientific problems, see [28-35]. For example, the author in Ref. [28] used the Laplace transform for getting the solution of a fractional telegraph equation. The formulation of the tumor fractional model using nonsingular kernel has been studied in Ref. [29]. In Ref. [30], the authors studied the solution of the fifth-order evolution equation in a fractional derivative. Using the Haar wavelet method, the authors in Ref. [31] obtained the solution of the Lotka-Volterra model related to population in a fractional derivative. Using Yang-Abdel-Aty-Cattani fractional operator, the authors discussed the solution of heat equations in Ref. [32]. The classical optical solution for the propagation model has been derived in Ref. [33]. The two interacting species model using numerical method is discussed in Ref. [34]. A fractional-order SEIR model using Genocchi has been studied in Ref. [35].
In this investigation, we considered the FF designed operators suggested in the sense of power law, exponential decay, and the Mittag-Leffler kernel. These suggested operators were found interesting from a researchers point of view and have been found reliable for chaotic order systems. Using these operators, many hidden attractors have been found which were not possible for ordinary fractional operators. It is clear that these operators attract various attractors by considering their different orders. Recently, fractal-fractional orders problems are formulated in Ref. [6,36], etc. It is wellknown that the results shown in Ref. [6,36] are the applications for a physical systems. Further, these operators were found best for data fitting; one can see [36]. Moreover, the results suggested in Refs. [6,36-38] demonstrate the combinations of these FF operators provide many insights about the dynamics, which can only be observed for the FF orders. With the above many useful applications of the FF operators, we give a chaotic model and present its analysis and provide a useful and effective scheme for its solution. We will find many useful numerical results by considering the combinations of the FF orders.
The work shown in this paper has been arranged sectionwise is as follows. We provide the details of the FF operators and model description in Secs. 2 and 3. Section 4 presents the novel solution technique for the FF model in Caputo, CF, and AB operators. We discussed the simulation results in Secs. 5 and 6, we summarized the finding.
2. Fundamental of FF operators
The present section determines the necessary results associated with fractal-fractional differential and their corresponding integral [6,36].
Definition 1. A function g(t) with fractal order β and fractional order α is continuous and differentiable over (m 1 ,m 2), then, in a power-law sense, the following is presented:
with
Definition 2. A function g(t) with fractal order β and fractional order α is continuous and differentiable over (m 1 ,m 2), then, in exponential decay kernel, the following is shown:
with
Definition 3. A function g(t) with fractal and fractional order β and α respectively, is continuous and differentiable over (m 1 ,m 2), then, the Mittag-Leffler kernel follows the formulate given by:
with α> 0, β≤ 1 ∈ ℕ and AB(α) = 1 −α+ (α/Γ(α)).
Definition 4. A function g(t) continuous over (m 1 ,m 2), with fractal order β and fractional order α then, for the power-law kernel, we have the following integral:
Definition 5. A function g(t) continuous over (m 1 ,m 2), with fractal order β and fractional order α then, the FF integral exponentially decaying kernel for f(t) is described by the following:
Definition 6. A function z(t) continuous on (a,b), with fractal order β and fractional order α then, for the Mittag-Leffler kernel, the following integral is suggested:
Lemma 1. The following result holds for fractional and fractal order α, and β, respectively,
Proof. Consider
Riemann-Liouville integral gives the following,
Therefore,
Lemma 2. The following result holds for fractional and fractal order α and β, respectively,
Proof. In order to show this result given by
Since
is differentiable, so, we can use
So,
3. The chaotic model
We consider a model which chaotic and has been presented in Ref. [39]. The model has three equations, and its differential equations form considered in Ref. [39] is shown below:
where x(t), y(t) and z(t) represent the state variables and become chaotic for the initial values of the variables x(0) = 0, y(0) = 1 and z(0) = 0. The fractal-fractional representation of the model (7) is given by:
where α represents the fractional order while β is fractal order. No equilibria exist for the model (8). Further, in the following, we are going to formulate the problem above using FF operators and present the numerical scheme for each case.
4. Numerical procedure for the chaotic model
This section presents briefly the procedure to handle the model formulated in FF with fractional operators sense, Caputo, CF, and AB. One by one, we apply the fractal operator to the model described above (8).
4.1. Numerical scheme for FF-Caputo model
We start the model to provide a numerical algorithm by considering the model in power-law case. To start the algorithm, we need first the write the model (8) in Volterra representations, because the fractional integral is differentiable, so in the Riemann-Liouville (RL) sense, we have,
we present the results in the following,
Consider now replacing the operator shown by Riemann-Liouville(RL) in the above system with the Caputo operator to get the initial conditions in the integer case. Further, we get the result by using RL fractional integral:
where
For t = t n +1, we have the following,
Simplifying the above results further, we get
where t(0) = 0. The functions δ β−1 f i (x,y,z,δ) for i = 1,2,3 shown in Eq. (13) are approximated in the interval [t j ,t j+1 ] using the Lagrangian interpolation concept, and getting the following,
Using the results shown in (14), we have the below system,
Further solution of system (15) leads to the following,
where
where n = 0 to N and j = 1 to n.
4.2. Numerical solution of FF-CF model
Now, we consider the model in the FF-CF case and present in detail the numerical scheme. For this, we need to write the model in FF exponential decay form using the following structure:
We obtain the following by using the Caputo-Fabrizio integral,
The following is presented at t n +1:
Further, we have the following:
The idea of Lagrange polynomial is used to obtain the following result:
Further, we obtain the following,
Further in next subsection, we use the FF model and present a numerical procedure in order to obtain their solution numerically. We provide the details of the scheme in the below subsection.
4.3. Numerical solution of FF-AB model
In order to have a numerical procedure to handle the FF model in Atangana-Baleanu (AB) operator sense, the model shown in (8) can be re-written using the following way:
The following is presented based on Atangana-Baleanu (AB) integral,
Further, at t n +1, we get,
We obtain the following after simplifying the integral in the above equations:
Now, approximating the expressions in (26) given by δ β−1 f 1(x,y,z,δ), δ β−1 f 2(x,y,z,δ) and δ β−1 f 3(x,y,z,δ) in the given interval [t j ,t j +1], the following numerical scheme is presented,
where
where n = 0 to N and j = 1 to n.
5. Discussion
We studied above the detailed algorithms for the solution of FF models numerically using the power law, exponential decay, and the Mittag-Leffler cases. For each operator, the scheme has been presented in a detailed manner. We use the numerical algorithm shown above in order to obtain the solution of the chaotic model in graphical form. Using the scheme presented in (16) for the FF system in Caputo sense, we used the scheme, and the corresponding graphical results are shown in Figs. 1-4 with different fractal and fractional orders. We see, for different orders of fractal and fractional, the model behaves chaotically. We also utilized a different set of fractal and fractional orders and obtained interesting chaotic results. For the model formulated in FF-CF derivative, we considered the suggested scheme (22), and the solution of the model has been shown graphically in Figs. 5-10, with arbitrary orders of fractal and fractional orders. One can observe that different combinations of the suggested operators provide useful numerical results. Further, the procedure given in (27) is used for the numerical results of the fractalfractional Atangana-Baleanu model with arbitrary order of fractal and fractional order parameters values and obtain the Figs. 11-14. From all these results, it can be concluded that the chaotic model provides different behaviors for the suggested values of the fractal and fractional orders. Using these fractal-fractional operators, many chaotic results and hidden attractors can be studied effectively. These operators can be applied to other nonlinear problems arising from physical sciences, social sciences, and engineering areas where one can find the best possible results for their proposed problems.
6. Conclusion
We presented a chaotic model in FF operators where the fractional operators were considered in Caputo, CF, and AB sense. We presented briefly the simulation of the chaotic model by using novel numerical approaches. We obtained the numerical results presented for the solution of the fractalfractional model in the sense of three different fractalfractional operators. We used the initial value for t = 0.1 and the step size h = 0.01. We presented the graphical results for each operator by selecting some suitable values for the fractal and fractional orders. We see that for each value of the fractal and fractional order, different chaotic behaviors of the model are observed. It is observed that by changing the values of the fractal and fractional orders, some new chaotic results are obtained graphically. This new analysis of the FF model in different operators bring new doors of investigations for scientists and researchers working in modeling electrical circuits theory. The results reveal that the fractalfractional operators and their application to a chaotic model provide useful results that cannot be seen for an ordinary operator. These operators can be beneficial in modeling other engineering fields and sciences related areas.