On a new modeling strategy: The logarithmically-exponential class of distributions

Abdulhakim A. Al-Babtain; Ibrahim Elbatal; Christophe Chesneau; Mohammed Elgarhy; Abdulhakim A. Al-Babtain; Ibrahim Elbatal; Christophe Chesneau; Mohammed Elgarhy

doi:10.3934/math.2021456

AIMS Mathematics

2021, Volume 6, Issue 7: 7845-7871. doi: 10.3934/math.2021456

Previous Article Next Article

Research article

On a new modeling strategy: The logarithmically-exponential class of distributions

1.
Department of Statistics and Operations Research, King Saud University, Riyadh 11362, Saudi Arabia
2.
Department of Mathematics and Statistics, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
3.
LMNO, Université de Caen Normandie, Campus II, Science 3, Caen 14032, France
4.
The Higher Institute of Commercial Sciences, Al mahalla Al kubra, Algarbia 31951, Egypt

Received: 28 January 2021 Accepted: 07 May 2021 Published: 18 May 2021
MSC : 60E05, 62E15, 62F30

In this paper, a promising modeling strategy for data fitting is derived from a new general class of univariate continuous distributions. This class is governed by an original logarithmically-exponential one-parameter transformation which has the ability to enhance some modeling capabilities of any parental distribution. In relation to the current literature, it appears to be a "limit case" of the well-established truncated generalized Fréchet generated class. In addition, it offers a natural alternative to the famous odd inverse exponential generated class. Some special distributions are presented, with particular interest in a novel heavy-tailed three-parameter distribution based on the Lomax distribution. Functional equivalences, modes analysis, stochastic ordering, functional expansions, moment measures, information measures and reliability measures are derived. From generic or real data, our modeling strategy is based on the new class combined with the maximum likelihood approach. We apply this strategy to the introduced modified Lomax model. The efficiency of the three parameter estimates is validated by a simulation study. Subsequently, two referenced real data sets are adjusted according to the rules of the art; the first one containing environmental data and the second one, financial data. In particular, we show that the proposed model is preferable to four concurrents also derived from the Lomax model, including the odd inverse exponential Lomax model.

Keywords:

Citation: Abdulhakim A. Al-Babtain, Ibrahim Elbatal, Christophe Chesneau, Mohammed Elgarhy. On a new modeling strategy: The logarithmically-exponential class of distributions[J]. AIMS Mathematics, 2021, 6(7): 7845-7871. doi: 10.3934/math.2021456

Related Papers:

[1]	C. Kavitha, A. Gowrisankar . Fractional integral approach on nonlinear fractal function and its application. Mathematical Modelling and Control, 2024, 4(3): 230-245. doi: 10.3934/mmc.2024019
[2]	Anil Chavada, Nimisha Pathak, Sagar R. Khirsariya . A fractional mathematical model for assessing cancer risk due to smoking habits. Mathematical Modelling and Control, 2024, 4(3): 246-259. doi: 10.3934/mmc.2024020
[3]	Yangtao Wang, Kelin Li . Exponential synchronization of fractional order fuzzy memristor neural networks with time-varying delays and impulses. Mathematical Modelling and Control, 2025, 5(2): 164-179. doi: 10.3934/mmc.2025012
[4]	Xue Zhang, Bo Sang, Bingxue Li, Jie Liu, Lihua Fan, Ning Wang . Hidden chaotic mechanisms for a family of chameleon systems. Mathematical Modelling and Control, 2023, 3(4): 400-415. doi: 10.3934/mmc.2023032
[5]	Shuli Zhang, Yansheng Liu . Existence of solutions for a class of fractional dynamical systems with two damping terms in Banach space. Mathematical Modelling and Control, 2023, 3(3): 168-180. doi: 10.3934/mmc.2023015
[6]	Abduljawad Anwar, Shayma Adil Murad . On the Ulam stability and existence of $L^p$ -solutions for fractional differential and integro-differential equations with Caputo-Hadamard derivative. Mathematical Modelling and Control, 2024, 4(4): 439-458. doi: 10.3934/mmc.2024035
[7]	Zhen Xin, Yuhe Yang, Qiaoxia Li . Controllability of nonlinear ordinary differential equations with non-instantaneous impulses. Mathematical Modelling and Control, 2022, 2(1): 1-6. doi: 10.3934/mmc.2022001
[8]	Rashid Jan, Normy Norfiza Abdul Razak, Sania Qureshi, Imtiaz Ahmad, Salma Bahramand . Modeling Rift Valley fever transmission: insights from fractal-fractional dynamics with the Caputo derivative. Mathematical Modelling and Control, 2024, 4(2): 163-177. doi: 10.3934/mmc.2024015
[9]	K. Venkatachalam, M. Sathish Kumar, P. Jayakumar . Results on non local impulsive implicit Caputo-Hadamard fractional differential equations. Mathematical Modelling and Control, 2024, 4(3): 286-296. doi: 10.3934/mmc.2024023
[10]	Saravanan Shanmugam, R. Vadivel, S. Sabarathinam, P. Hammachukiattikul, Nallappan Gunasekaran . Enhancing synchronization criteria for fractional-order chaotic neural networks via intermittent control: an extended dissipativity approach. Mathematical Modelling and Control, 2025, 5(1): 31-47. doi: 10.3934/mmc.2025003

Abstract

1. Introduction

The role of mathematical modeling is significant in different fields of science. Mathematical modeling has been identified as a helpful research area for investigating many aspects of circuit problems ^[1,2,3]. Many researchers have used integer-order derivatives, but integer order fails to describe the desired behavior of circuit problems. This way, most researchers have been using fractional derivative operators to obtain more realistic results of physical problems ^[4,5,6]. It is seen that models involving non-integer order differential equations and integrals are more natural than to the ordinary derivative model ^[7,8,9].

Applications of fractional calculus to many problems have been studied extensively in previously published papers ^[10,11]. All of these mentioned papers analyzed engineering problems and infectious diseases to decrease the transmission of infectious diseases. Researchers have proved that fractional derivatives can describe such a circuit problem efficiently and provide better results as compared to ordinary derivatives ^{[12,13,14,15,16]}. Therefore, we want to study a memristor-based chaotic system in fractional type derivatives. This paper considers two fractional derivatives: Atangana-Baleanu–Caputo and Caputo-Fabrizio derivatives. In ^[17], the author proposed a circuit integrator model as

$\begin{equation} \begin{cases} _{0}^{C}\mathscr{D}_{\zeta}^{\alpha}[\mathscr{X}(\zeta)] = \mathscr{Y}(\zeta), \\ _{0}^{C}\mathscr{D}_{\zeta}^{\alpha}[\mathscr{Y}(\zeta)] = -\frac{1}{3}\left(\mathscr{X}(\zeta)+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta)-1\right)\mathscr{Y}(\zeta)\right), \\ _{0}^{C}\mathscr{D}_{\zeta}^{\alpha}[\mathscr{Z}(\zeta)] = -\frac{3}{5}\mathscr{Z}(\zeta)-\mathscr{Y}(\zeta)+\mathscr{Z}(\zeta)\mathscr{Y}(\zeta). \end{cases} \end{equation}$

(1.1)

For solving nonlinear fractional differential equations, numerical techniques have been considered effective mathematical tools ^[16,18]. An Adam-Bashforth numerical method has been used in many fields of applied science such as engineering and epidemiology problems ^[19,20,21]. However, some numerical scheme are not fully applied to differential equations with local and non-local operators, as they were designed only for integer-order derivatives ^[19]. In this paper, we study the system (1.1) under ABC and CF types of fractal-fractional derivatives with order $\alpha$ .

2. Analysis of the fractal-fractional model (1.1) in ABC sense

In this section, we study the chaotic system (1.1) under ABC type fractal-fractional derivatives with order $\alpha$ as follows:

$\begin{equation} \begin{cases} _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \mathscr{Y}(\zeta), \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Y}(\zeta)] = -\frac{1}{3}\left(\mathscr{X}(\zeta)+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta)-1\right)\mathscr{Y}(\zeta)\right), \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Z}(\zeta)] = -\frac{3}{5}\mathscr{Z}(\zeta)-\mathscr{Y}(\zeta)+\mathscr{Z}(\zeta)\mathscr{Y}(\zeta). \end{cases} \end{equation}$

(2.1)

The following are the circuital equations related to the system (2.1):

$\begin{equation} \begin{cases} _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[-\left(\frac{\mathcal{R}_{2}}{\mathcal{R}_{1}\mathcal{R}_{3}}\right)\mathscr{Y}\right], \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Y}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[\left(\frac{\mathcal{R}_{7}}{\mathcal{R}_{8}\mathcal{R}_{5}}\right)\mathscr{Y}-\left(\frac{1}{\mathcal{R}_{4}}\right)\mathscr{X}-\left(\frac{1}{\mathcal{R}_{6}}\right)\mathscr{Z}^{2}\mathscr{Y}\right], \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[\left(\frac{\mathcal{R}_{12}}{\mathcal{R}_{13}\mathcal{R}_{9}}\right)\mathscr{YZ}-\left(\frac{1}{\mathcal{R}_{10}}\right)\mathscr{Y}-\left(\frac{1}{\mathcal{R}_{11}}\right)\mathscr{Z}\right], \end{cases} \end{equation}$

(2.2)

where parameters are given as follows: $\mathcal{R}_{1} = \mathcal{R}_{3} = \mathcal{R}_{2} = \mathcal{R}_{5} = \mathcal{R}_{9} =$ $\mathcal{R}_{7} = \mathcal{R}_{10} = \mathcal{R}_{13} = \mathcal{R}_{12} = 10K\Omega, \mathcal{R}_{6} = \mathcal{R}_{8} = 10.62K\Omega, \mathcal{R}_{4} = 11.62K\Omega, \mathcal{R}_{11} = 11.49K\Omega.$

Let $\mathscr{P}(\zeta)$ be a continuous function in $(a, b)$ ^[22], then the ABC fractal-fractional derivative of $\mathscr{P}(\zeta)$ with order $0 < \alpha \leq 1$ is given by

$\begin{equation} _{a}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)] = \frac{\mathcal{B}(\alpha)}{1-\alpha}\frac{d}{d\zeta^{\gamma}}\int_{a}^{\zeta}\mathscr{P}(\xi)E_{\alpha}\left[-\alpha\frac{(\zeta-\xi)^{\alpha}}{1-\alpha}\right]d\xi, \end{equation}$

(2.3)

where $\mathcal{B}(\alpha)$ is normalization function such that $\mathcal{B}(0) = \mathcal{B}(1) = 1$ and $\gamma \le1$ .

Let $0 < \alpha \leq 1$ ; $\gamma\leq 1$ ; then, the fractal-fractional integral is given by

$\begin{equation} \begin{array}{lll} &_{0}^{ABC}\mathscr{I}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)]\\ = &\frac{\gamma(1-\alpha)\zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\mathscr{P}(\zeta)+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}\int_{0}^{\zeta}\xi^{\gamma-1}\mathscr{P}(\xi)\left(\zeta-\xi\right)^{\alpha-1}d\xi. \end{array} \end{equation}$

(2.4)

3. Existence of the solution

In this section, we prove that the model (2.1) has at least one solution.

Consider the model (2.1) as

$\begin{equation} \begin{cases} _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \gamma \zeta^{\gamma-1}\mathfrak{X}, \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Y}(\zeta)] = \gamma \zeta^{\gamma-1}\mathfrak{Y}, \\ _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Z}(\zeta)] = \gamma \zeta^{\gamma-1}\mathfrak{Z}, \end{cases} \end{equation}$

(3.1)

where

$\begin{eqnarray*} \begin{cases} \mathfrak{X}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = \mathscr{Y}(\zeta)\\ \mathfrak{Y}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = -\frac{1}{3}\left(\mathscr{X}(\zeta)+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta)-1\right)\mathscr{Y}(\zeta)\right)\\ \mathfrak{Z}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = -\frac{3}{5}\mathscr{Z}(\zeta)-\mathscr{Y}(\zeta)+\mathscr{Z}(\zeta)\mathscr{Y}(\zeta). \end{cases} \end{eqnarray*}$

We can write system (3.1) as:

$\begin{equation} \begin{cases} _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)] = \gamma \zeta^{\gamma-1}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right), & \zeta\geq0, \\ \mathscr{P}(0) = \mathscr{P}_{0}. \end{cases} \end{equation}$

(3.2)

By replacing $_{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha}$ with $_{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}$ and applying a fractional integral, we get

$\begin{equation} \begin{alignedat}{1} \mathscr{P}(\zeta)-\mathscr{P}(0) = & \frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)+\frac{\alpha\gamma}{\mathcal{B}(\alpha) \varGamma(\alpha)}\\ &\int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta-\xi\right)^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi.\end{alignedat} \end{equation}$

(3.3)

We define a Banach space $\Omega = \mathscr{F}\times\mathscr{F}\times\mathscr{F},$ for the existence theory, where $\mathscr{F} = \mathbb{C}[0, T]$ with the norm

$\left\Vert \mathscr{P}\right\Vert = \max_{\zeta\in[0, T]}\left|\mathscr{X}(\zeta)+\mathscr{Y}(\zeta)+\mathscr{Z}(\zeta)\right|.$

We define an operator $\mathscr{T}:\Omega\rightarrow \Omega$ as follows:

$\begin{equation} \begin{alignedat}{1}\mathscr{T}\mathscr{P}(\zeta) = & \mathscr{P}(0)+ \frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)\\&+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta-\xi\right)^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi.\end{alignedat} \end{equation}$

(3.4)

The Lipschitz condition is applied to non-linear function $\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)$ as:

● For each $\mathscr{P}\in\Omega,$ there exist constants $\mathcal{M} > 0,$ and $\mathcal{M}_{a}$ such that

$\begin{equation} \left|\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)\right|\leq\mathcal{M}\left|\mathscr{P}(\zeta)\right|+\mathcal{M}_{a}. \end{equation}$

(3.5)

● For each $\mathscr{P}, \mathscr{P}_{a}\in\Omega,$ there exists a constant $\mathcal{M}_{a} > 0$ such that

$\begin{equation} \left|\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)-\mathcal{G}\left(\zeta, \mathscr{P}_{a}(\zeta)\right)\right|\leq\mathcal{M}_{a}\left|\mathscr{P}(\zeta)-\mathscr{P}_{a}(\zeta)\right|. \end{equation}$

(3.6)

Theorem 3.1. Assume that the condition (3.5) is satisfied. Let $\Xi:[0, T]\times\varOmega\rightarrow \mathcal {R}$ be a function that is continuous (2.1) must have at least one solution.

Proof. First, we prove that the operator $\mathscr{T}$ described by Eq. (3.4) is fully continuous. Because $\Xi$ is continuous, $\mathscr{T}$ is also continuous. Let $\mathscr{A} = \left\{ \mathscr{P}\in\varOmega:\left\Vert \mathscr{P}\right\Vert \leq M, M > 0\right\}$ . Now for any $\mathscr{P}\in\varOmega,$ we have

$\begin{align*} &\left\Vert \mathscr{T}\mathscr{P}\right\Vert \\ = & \max\limits_{\zeta\in[0, T]}\Big{|}\mathscr{P}(0)+\frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)\\ &+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta-\xi\right)^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi\Big{|}, \\ \leq & \mathscr{P}(0)+\frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert +\mathcal{M}_{a}\right)\\ &+\max\limits_{\zeta\in[0, T]}\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta-\xi\right)^{\gamma-1}\left|\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)\right|d\xi, \\ \leq & \mathscr{P}(0)+\frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert +\mathcal{M}_{a}\right)\\ &+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert +\mathcal{M}_{a}\right)\zeta_{2}^{\alpha+\gamma-1}\mathscr{B}(\alpha, \gamma)\leq\mathcal{R}. \end{align*}$

where $\mathscr{B}(\alpha, \gamma)$ represents the beta function. Thus the operator $\mathscr{T}$ is uniformly period,

For eqicontinuity of $\mathscr{T},$ let us take $\mathscr{K}_{1} < \mathscr{K}_{2}\leq$ T.

Then consider

$\begin{align*} &\left\Vert\mathscr{T}\mathscr{P}(\zeta_{2})-\mathscr{T}\mathscr{P}_{a}(\zeta_{1})\right\Vert\\ = &\max\limits_{\zeta\in[0, T]}\Big{|}\frac{\gamma \zeta_{2}^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta_{2}, \mathscr{P}(\zeta_{2})\right)\\ &+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta_{2}-\xi\right)^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi\\ &-\frac{\gamma \zeta_{1}^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta_{1}, \mathscr{P}(\zeta_{1})\right)\\ &-\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta_{1}-\xi\right)^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi\Big{|}\\ \leq & \frac{\gamma \zeta_{1}^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert+\mathcal{M}_{a}\right)+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert\right.\\ &\left.+\mathcal{M}_{a}\right)\zeta_{1}^{\alpha+\gamma-1}\mathscr{B}(\alpha, \gamma)-\frac{\gamma \zeta_{2}^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert+\mathcal{M}_{a}\right)\\ &-\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}\left(\mathcal{M}\left\Vert \mathscr{P}\right\Vert +\mathcal{M}_{a}\right)\zeta_{2}^{\alpha+\gamma-1}\mathscr{B}(\alpha, \gamma). \end{align*}$

when $\zeta_{1}\rightarrow \zeta_{2},$ then $\left|\mathscr{T}\mathscr{P}(t_{2})-\mathscr{T}\mathscr{P}_{a}(\zeta_{1})\right|\rightarrow0,$ Hence $\mathscr{T}$ is equicontinuous. So, by the Arzela-Ascoli theorem is completely continuous. Thus, by Schauder's fixed point result the proposed model has at least one solution.

4. Uniqueness of the solution

Theorem 4.1. Let (3.6) hold. If $C < 1,$ where

$C = \left(\frac{\gamma T^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}T^{\alpha+\gamma-1}\mathscr{B}(\alpha, \gamma)\right)\mathcal{M}_{a},$ then the considered model has a unique solution

Proof. For $\mathscr{P}, \mathscr{P}_{a}\in\Omega$ , we have

$\begin{align*} & \left\Vert \mathscr{T}\mathscr{P}-\mathscr{T}\mathscr{P}_{a}\right\Vert \\ = & \max\limits_{\zeta\in[0, T]}\Big{|}\frac{\gamma \zeta^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}\left(\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)-\mathcal{G}\left(\zeta, \mathscr{P}_{a}(\zeta)\right)\right)\\ &+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)} \int_{0}^{\zeta}\xi^{\gamma-1}\left(t-\xi\right)^{\gamma-1}\left(\mathcal{G}\left.(\xi, \mathscr{P}(\xi)\right)\right.\\ & -\mathcal{G}\left(\xi, \mathscr{P}_{a}(\xi)\right)) d\xi\Big{|}\\ \leq & \frac{\gamma T^{\gamma-1}\left(1-\alpha\right)}{\mathcal{B}(\alpha)}+\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}T^{\alpha+\gamma-1}\mathscr{B}(\alpha, \gamma)\left\Vert \mathscr{P}-\mathscr{P}_{a}\right\Vert\\ \leq & C\left\Vert \mathscr{P}-\mathscr{P}_{a}\right\Vert. \end{align*}$

As a result, $\mathscr{T}$ is a contraction. So, the model has a unique solution according to the Banach contraction principle.

5. Numerical procedure

Toufik and Atangana ^[23,24] presented a numerical scheme. In this section, we illustrate the numerical scheme in details:

$\begin{equation} \begin{cases} _{0}^{ABC}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)] = \mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right), & \zeta\geq0, \\ \mathscr{P}(0) = \mathscr{P}_{0}. \end{cases} \end{equation}$

(5.1)

Using the fundamental theorem of fractional calculus Eq. (5.1) can be written as:

$\begin{equation} \begin{alignedat}{1}\mathscr{P}(\zeta)-\mathscr{P}(0) = &\frac{\gamma \zeta^{\gamma-1}(1-\alpha)}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)+\frac{\gamma\alpha}{\mathcal{B}(\alpha)\varGamma(\alpha)}\\&\int_{0}^{\zeta}\xi^{\gamma-1}\left(\zeta-\xi\right)^{\alpha-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi.\end{alignedat} \end{equation}$

(5.2)

At a point $\zeta = \zeta_{\mathscr{K}+1}$ , for $\mathscr{K} = 0, 1, 2, ...,$ equation Eq. (5.2) becomes

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}(\zeta_{\mathscr{K}+1})-\mathscr{P}(0)\\ = & \frac{\gamma(1-\alpha)\zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\alpha\gamma}{\mathcal{B}(\alpha)\varGamma(\alpha)}\int_{0}^{\zeta_{\mathscr{K}+1}}\xi^{\gamma-1}\left(\zeta_{\mathscr{K}+1}-\xi\right)^{\alpha-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi, \end{alignedat} \end{equation}$

(5.3)

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}_{\mathscr{K}+1} = \mathscr{P}(\zeta_{\mathscr{K}+1}) \\ = &\mathscr{P}(0)+\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathcal{G}\left(t_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)+\frac{\alpha}{\mathcal{B}(\gamma\alpha)\varGamma(\alpha)}\\&\sum\limits_{i = 0}^{\mathscr{K}}\int_{\zeta_{i}}^{\zeta_{i+1}}\xi^{\gamma-1}\left(\zeta_{\mathscr{K}+1}-\xi\right)^{\alpha-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi. \end{alignedat} \end{equation}$

(5.4)

We can approximated the function $\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)$ in the interval $[\zeta_{i}, \zeta_{i+1}].$

$\begin{equation} \begin{alignedat}{1}q_{\mathscr{K}} = &\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right) \\ = &\frac{\xi-\zeta_{i-1}}{\zeta_{i}-\zeta_{i-1}}\mathcal{G}\left(\zeta_{i}, \mathscr{P}_{i}\right)+\frac{\xi-\zeta_{i}}{\zeta_{i-1}-\zeta_{i}}\mathcal{G}\left(\zeta_{i-1}, \mathscr{P}_{i-1}\right).\end{alignedat} \end{equation}$

(5.5)

Substituting the above value in Eq. (5.4), we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}_{\mathscr{K}+1}\\ = &\mathscr{P}(0)+\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)\varGamma(\alpha)}\sum\limits_{i = 0}^{\mathscr{K}}\left(\frac{\mathcal{G}\left(\zeta_{i}, \mathscr{P}(\zeta_{i})\right)}{n}\int_{\zeta_{i}}^{\zeta_{i+1}}\xi^{\gamma-1}\left(\zeta-\zeta_{i-1}\right)\right.\\&\left.\left(\zeta_{\mathscr{K}+1}-\xi\right)^{\alpha-1}d\xi\right. -\frac{\mathcal{G}\left(\zeta_{i-1}, \mathscr{P}(\zeta_{i-1})\right)}{n}\\&\left.\int_{\zeta_{i}}^{\zeta_{i+1}}\xi^{\gamma-1}\left(\zeta-\zeta_{i-1}\right)\left(\zeta_{\mathscr{K}+1}-\xi\right)^{\alpha-1}d\xi\right). \end{alignedat} \end{equation}$

(5.6)

After substituting $n = \zeta_{i}-\zeta_{i-1}$ and simplification, we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}_{\mathscr{K}+1}\\ & = \mathscr{P}_{0}+\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)+\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\\&\sum\limits_{i = 0}^{\mathscr{K}}\left[\frac{n^{\alpha}\mathcal{G}\left(\zeta_{i}, \mathscr{P}(\zeta_{i})\right)}{\varGamma(\alpha+2)}\left((\mathscr{K}-i+1)^{\alpha}(\mathscr{K}+2-i+\alpha)\right)\right.\\ & \left(\mathscr{K}-i\right)^{\alpha}\left(\mathscr{K}-i+2+2\alpha\right)-\frac{n^{\alpha}\mathcal{G}\left(\zeta_{i-1}, \mathscr{P}(\zeta_{i-1})\right)}{\varGamma(\alpha+2)}\\ & \left.\left((\mathscr{K}-i+1)^{\alpha+1}-(\mathscr{K}-i)^{\alpha}(\mathscr{K}-i+1+\alpha\right)\right]. \end{alignedat} \end{equation}$

(5.7)

For simplicity, we define

$\begin{eqnarray*} \begin{cases} \mathfrak{X}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = \mathscr{Y}(\zeta)\\ \mathfrak{Y}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = -\frac{1}{3}\left(\mathscr{X}(\zeta)+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta)-1\right)\mathscr{Y}(\zeta)\right)\\ \mathfrak{Z}(\mathscr{X}, \mathscr{Y}, \mathscr{Z}, \zeta) = -\mathscr{Y}(t)-\frac{3}{5}\mathscr{Z}(\zeta)+\mathscr{Z}(\zeta)\mathscr{Y}(\zeta). \end{cases} \end{eqnarray*}$

The numerical scheme for the non-linear integrator circuit under the ABC fractal-fractional derivative is given below:

$\begin{equation} \begin{alignedat}{1}&\mathscr{X}_{\mathscr{K}+1} = \mathscr{X}_{0}+\frac{(1-\alpha)\gamma t^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathfrak{X}\left(t_{\mathscr{K}}, \mathscr{X}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\sum\limits_{i = 0}^{\mathscr{K}}\left[\frac{n^{\alpha}\mathfrak{X}\left(\zeta_{i}, \mathscr{X}(\zeta_{i})\right)}{\varGamma(\alpha+2)}\left((\mathscr{K}+1-i)^{\alpha}(\mathscr{K}-i+2+\alpha)\right)\right.\\ & \left(\mathscr{K}-i\right)^{\alpha}\left(\mathscr{K}+2-i+2\alpha\right)-\frac{n^{\alpha}\mathfrak{X}\left(\zeta_{i-1}, \mathscr{X}(\zeta_{i-1})\right)}{\varGamma(\alpha+2)}\\ & \left.\left((\mathscr{K}-i+1)^{\alpha+1}-(\mathscr{K}-i)^{\alpha}(\mathscr{K}+1-i+\alpha\right)\right]. \end{alignedat} \end{equation}$

(5.8)

$\begin{equation} \begin{alignedat}{1}&\mathscr{Y}_{\mathscr{K}+1} = \mathscr{Y}_{0}+\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathfrak{Y}\left(\zeta_{\mathscr{K}}, \mathscr{Y}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\sum\limits_{i = 0}^{\mathscr{K}}\left[\frac{n^{\alpha}\mathfrak{Y}\left(\zeta_{i}, \mathscr{Y}(\zeta_{i})\right)}{\varGamma(\alpha+2)}\left((\mathscr{K}+1-i)^{\alpha}(\mathscr{K}-i+2+\alpha)\right)\right.\\ & \left(\mathscr{K}-i\right)^{\alpha}\left(\mathscr{K}+2-i+2\alpha\right)-\frac{n^{\alpha}\mathfrak{Y}\left(\zeta_{i-1}, \mathscr{Y}(\zeta_{i-1})\right)}{\varGamma(\alpha+2)}\\ & \left.\left((\mathscr{K}-i+1)^{\alpha+1}-(\mathscr{K}-i)^{\alpha}(\mathscr{K}-i+1+\alpha\right)\right]. \end{alignedat} \end{equation}$

(5.9)

$\begin{equation} \begin{alignedat}{1}&\mathscr{Z}_{\mathscr{K}+1} = \mathscr{P}_{0}+\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\gamma\alpha)}\mathfrak{Z}\left(\zeta_{\mathscr{K}}, \mathscr{Z}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\sum\limits_{i = 0}^{\mathscr{K}}\left[\frac{n^{\alpha}\mathfrak{Z}\left(\zeta_{i}, \mathscr{Z}(\zeta_{i})\right)}{\varGamma(\alpha+2)}\left((\mathscr{K}-i+1)^{\alpha}(\mathscr{K}+2-i+\alpha)\right)\right.\\ & \left(\mathscr{K}-i\right)^{\alpha}\left(\mathscr{K}-i+2+2\alpha\right)-\frac{n^{\alpha}\mathfrak{Z}\left(\zeta_{i-1}, \mathscr{Z}(\zeta_{i-1})\right)}{\varGamma(\alpha+2)}\\ & \left.\left((\mathscr{K}-i+1)^{\alpha+1}-(\mathscr{K}-i)^{\alpha}(\mathscr{K}+1-i+\alpha\right)\right]. \end{alignedat} \end{equation}$

(5.10)

Now, we present some numerical results displaying complex behaviors of the novel fractal-fractional system (2.1). to present the graphical results for the ABC fractal-fractional model (2.1) by considered the above novel numerical scheme with the initial conditions $\mathscr{X}(0) = 0.2$ , $\mathscr{Y}(0) = 0.2$ , and $\mathscr{Z}(0) = 0.3$ . The simulation results in both cases include the time-domain responses of the state variables and the two- and three- dimensional phase portraits. The plots in and describe the dynamical behavior of $\mathscr{X}(\zeta)$ , $\mathscr{Y}(\zeta)$ , and $\mathscr{Z}(\zeta)$ for different fractal-fractional order in the sense of the ABC fractal operator with $\mathcal{R}_{1} = \mathcal{R}_{3} = \mathcal{R}_{2} = \mathcal{R}_{5} = \mathcal{R}_{9} =$ $\mathcal{R}_{7} = \mathcal{R}_{10} = \mathcal{R}_{12} = \mathcal{R}_{13} = 20K\Omega, \mathcal{R}_{6} = \mathcal{R}_{8} = 11.62K\Omega, \mathcal{R}_{4} = 12.62K\Omega, \mathcal{R}_{11} = 15.49K\Omega.$ The simulation results in for the two- and three-dimensional phase portraits are displayed in the plots of Figures 1 and 3. In Figure (2)(a-c), we can observe that when the order of the ABC fractal-derivative decreases, the dynamics of the ABC fractal-fractional system (2.1) decrease as well. This means that the ABC has an acceleration effect in the system process when the fractional order decreases. The figures mentioned above confirm the effectiveness of the proposed numerical scheme in exhibiting both chaotic and non-chaotic behaviors of the novel fractal fractional dynamical system (2.1).

Figure 1. Dynamical behavior of the chaotic attractor for

$\mathscr{X}(0) = 0.2, \mathscr{Y}(0) = 0.2$ and

$\mathscr{Z}(0) = 0.3$ , as obtained for

$\gamma = 1$ and

$\alpha = 0.95$ by using the ABC fractal-fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 2. Numerical simulation of (2.1) for

$\gamma = 1$ and different values of

$\alpha = 0.95$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. Dynamical behavior of the chaotic attractor for

$\mathscr{X}(0) = 0.2, \mathscr{Y}(0) = 0.2$ and

$\mathscr{Z}(0) = 0.3$ , as obtained for

$\gamma = 1$ and

$\alpha = 0.85$ by using the ABC fractal-fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 4. Numerical simulation of (2.1) for

$\gamma = 0.85$ and different values of

$\alpha = 0.85$ .

DownLoad: Full-Size Img PowerPoint

6. Analysis of the chaotic system (1.1) under the Caputo-Fabrizo fractal-fractional derivative

In this section, we study the fractal-fractional chaotic system (1.1) under the CF derivative with the fractional order $\alpha$ as follows:

$\begin{equation} \begin{cases} _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \mathscr{Y}(\zeta), \\ _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Y}(\zeta)] = -\frac{1}{3}\left(\mathscr{X}(\zeta)+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta)-1\right)\mathscr{Y}(\zeta)\right), \\ _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Z}(\zeta)] = -\frac{3}{5}\mathscr{Z}(\zeta)-\mathscr{Y}(\zeta)+\mathscr{Z}(\zeta)\mathscr{Y}(\zeta). \end{cases} \end{equation}$

(6.1)

The circuital equations associated with the system equations are given below:

$\begin{equation} \begin{cases} _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{X}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[-\left(\frac{\mathcal{R}_{2}}{\mathcal{R}_{1}\mathcal{R}_{3}}\right)\mathscr{Y}\right], \\ _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Y}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[\left(\frac{\mathcal{R}_{7}}{\mathcal{R}_{8}\mathcal{R}_{5}}\right)\mathscr{Y}-\left(\frac{1}{\mathcal{R}_{4}}\right)\mathscr{X}-\left(\frac{1}{\mathcal{R}_{6}}\right)\mathscr{Z}^{2}\mathscr{Y}\right], \\ _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{Z}(\zeta)] = \frac{1}{\mathcal{C}_{0}}\left[\left(\frac{\mathcal{R}_{12}}{\mathcal{R}_{13}\mathcal{R}_{9}}\right)\mathscr{YZ}-\left(\frac{1}{\mathcal{R}_{10}}\right)\mathscr{Y}-\left(\frac{1}{\mathcal{R}_{11}}\right)\mathscr{Z}\right], \end{cases} \end{equation}$

(6.2)

where the parameters values are can be considered as follows: $\mathcal{R}_{1} = \mathcal{R}_{3} = \mathcal{R}_{2} = \mathcal{R}_{5} = \mathcal{R}_{9} = \mathcal{R}_{7} =$ $\mathcal{R}_{12} = \mathcal{R}_{10} = \mathcal{R}_{13} = 10K\Omega, \mathcal{R}_{6} = \mathcal{R}_{8} = 10.62K\Omega, \mathcal{R}_{4} = 11.62K\Omega, \mathcal{R}_{11} = 11.49K\Omega.$

Let $\mathscr{P}(\zeta)\in\mathcal{H}\left(a, b\right)$ ^[25]; then, the Caputo-Fabrizo fractal fractional derivative can be defined as

$\begin{equation} \begin{alignedat}{1}_{a}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)] = & \frac{\mathcal{B}(\alpha)}{1-\alpha}\int_{a}^{\zeta}\frac {d\mathscr{P}(\xi)}{d\xi^{\gamma}}\exp\left[-\alpha\frac{(\zeta-\xi)}{1-\alpha}\right]d\xi, \end{alignedat} \end{equation}$

(6.3)

where, $\gamma \leq 1$ and $0 < \alpha \leq 1$ $\mathcal{B}(0) = \mathcal{B}(1) = 1.$

Here, we develop the numerical scheme for the fractional model (6.2) by using the methods described in ^[23,26]. We contemplate the following equation:

$\begin{equation} \begin{cases} _{0}^{CF}\mathscr{D}_{\zeta}^{\alpha, \gamma}[\mathscr{P}(\zeta)] = \mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right), & \zeta\geq0, \\ \mathscr{P}(0) = \mathscr{P}_{0}, \end{cases} \end{equation}$

(6.4)

where $\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)$ represents the right hand sides of the equations in the model (1.1). $\mathscr{P}(\zeta)$ stands for $\mathscr{X}(\zeta), \mathscr{Y}(\zeta)$ and $\mathscr{Z}(\zeta)$ .

Equation (6.4) can be expressed by using the fundamental theorem of fractional calculus:

$\begin{equation} \begin{alignedat}{1}\mathscr{P}(\zeta) = & \mathscr{P}(0)+ \frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta, \mathscr{P}(\zeta)\right)\\ &+\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\int_{a}^{\zeta}\xi^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi. \end{alignedat} \end{equation}$

(6.5)

For $\zeta = \zeta_{\mathscr{K}+1}, \mathscr{K} = 0, 1, 2, ...,$ we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}(\zeta_{\mathscr{K}+1})-\mathscr{P}(0)\\ = & \frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)\varGamma(\alpha)}\int_{0}^{\zeta_{\mathscr{K}+1}}\xi^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi. \end{alignedat} \end{equation}$

(6.6)

Now at $\zeta = \zeta_{\mathscr{K}}, \mathscr{K} = 0, 1, 2, ...,$ we have:

$\begin{equation} \begin{alignedat}{1}\mathscr{P}(\zeta_{\mathscr{K}})-\mathscr{P}_{0} = & \frac{(1-\alpha)\gamma t^{\gamma-1}}{\mathcal{B}(\alpha)}\mathcal{G}\left(\zeta_{\mathscr{K}-1}, \mathscr{P}(\zeta_{\mathscr{K}-1})\right)\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\int_{0}^{\zeta_{\mathscr{K}}}\xi^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi. \end{alignedat} \end{equation}$

(6.7)

Using the successive terms with difference from the above equations, we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}(\zeta_{\mathscr{K}+1})-\mathscr{P}(\zeta_{\mathscr{K}})\\ = & \frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\left[\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)-\mathcal{G}\left(\zeta_{\mathscr{K}-1}, \mathscr{P}(\zeta_{\mathscr{K}-1})\right)\right]\\ & +\frac{\gamma\alpha}{\mathcal{B}(\alpha)}\int_{\zeta_{\mathscr{K}}}^{\zeta_{\mathscr{K}+1}}\xi^{\gamma-1}\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)d\xi. \end{alignedat} \end{equation}$

(6.8)

Let the function $\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right)$ be closed in the interval $\left[\zeta_{\mathscr{K}}, \zeta_{\mathscr{K}+1}\right]:$ then, the function can be approximated as:

$\begin{equation} \begin{alignedat}{1}q_{\mathscr{K}}\cong\mathcal{G}\left(\xi, \mathscr{P}(\xi)\right) & = \frac{\xi-\zeta_{j-1}}{\zeta_{j}-\zeta_{j-1}}\mathcal{G}\left(\zeta_{j}, \mathscr{P}_{j}\right)\\ &+\frac{\xi-\zeta_{j}}{\zeta_{j-1}-\zeta_{j}}\mathcal{G}\left(\zeta_{j-1}, \mathscr{P}_{j-1}\right).\end{alignedat} \end{equation}$

(6.9)

Substituting the above approximation in Eq. (6.8), we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}_{\mathscr{K}+1}-\mathscr{P}_{\mathscr{K}} \\ = &\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}\left[\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)-\mathcal{G}\left(\zeta_{\mathscr{K}-1}, \mathscr{P}(\zeta_{\mathscr{K}-1})\right)\right]\\& +\frac{\alpha}{\mathcal{B}(\alpha)} \int_{\zeta_{i}}^{\zeta_{i+1}}\xi^{\gamma-1}\left(\frac{\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}_{\mathscr{K}}\right)}{\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1}}\left(\xi-\zeta_{\mathscr{K}-1}\right)\right.\\ &\left.-\frac{\mathcal{G}\left(\zeta_{\mathscr{K}-1}, \mathscr{P}_{\mathscr{K}-1}\right)}{\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1}}\left(\xi-\zeta_{\mathscr{K}}\right)\right)d\xi. \end{alignedat} \end{equation}$

(6.10)

After simplification, we get

$\begin{equation} \begin{alignedat}{1}&\mathscr{P}_{\mathscr{K}+1} = \mathscr{P}_{0}\\&+\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{3(t_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\mathcal{G}\left(\zeta_{\mathscr{K}}, \mathscr{P}(\zeta_{\mathscr{K}})\right)\\ & -\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{\gamma\alpha((\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\mathcal{G}\left(\zeta_{\mathscr{K}-1}, \mathscr{P}(\zeta_{\mathscr{K}-1})\right). \end{alignedat} \end{equation}$

(6.11)

Now, applying the above numerical scheme for the fractional non-linear integrator circuit gives

$\begin{equation} \begin{alignedat}{1}&\mathscr{X}(\zeta_{\mathscr{K}+1})\\ = & \mathscr{X}_{0}+\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{3(\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\mathscr{Y}(\zeta_{\mathscr{K}})\\ & -\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{\gamma\alpha((\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\mathscr{Y}(\zeta_{\mathscr{K}-1}), \end{alignedat} \end{equation}$

(6.12)

$\begin{equation} \begin{alignedat}{1}&\mathscr{Y}(\zeta_{\mathscr{K}+1})\\ = & \mathscr{Y}_{0}+\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{3(\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\\ &\left(-\frac{1}{3}\left(\mathscr{X}(\zeta_{\mathscr{K}})+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta_{\mathscr{K}})-1\right)\mathscr{Y}(\zeta_{\mathscr{K}})\right)\right)\\ & -\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{\gamma\alpha((\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\\ &\left(-\frac{1}{3}\left(\mathscr{X}(\zeta_{\mathscr{K}-1})+\frac{3}{2}\left(\mathscr{Z}^{2}(\zeta_{\mathscr{K}-1})-1\right)\mathscr{Y}(\zeta_{\mathscr{K}-1})\right)\right), \end{alignedat} \end{equation}$

(6.13)

$\begin{equation} \begin{alignedat}{1}&\mathscr{Z}(\zeta_{\mathscr{K}+1})\\ = & \mathscr{Z}_{0}+\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{3(\zeta_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\\&\left(-\mathscr{Y}(\zeta_{\mathscr{K}})-\frac{3}{5}\mathscr{Z}(\zeta_{\mathscr{K}})+\mathscr{Z}(\zeta_{\mathscr{K}})\mathscr{Y}(\zeta_{\mathscr{K}})\right)\\ & -\left(\frac{(1-\alpha)\gamma \zeta^{\gamma-1}}{\mathcal{B}(\alpha)}+\frac{\gamma\alpha((t_{\mathscr{K}}-\zeta_{\mathscr{K}-1})}{2\mathcal{B}(\alpha)}\right)\\&\left(-\mathscr{Y}(\zeta_{\mathscr{K}-1})-\frac{3}{5}\mathscr{Z}(\zeta_{\mathscr{K}-1})+\mathscr{Z}(\zeta_{\mathscr{K}-1})\mathscr{Y}(\zeta_{\mathscr{K}-1})\right). \end{alignedat} \end{equation}$

(6.14)

We present some numerical simulations to study the complex behaviors of the fractal-fractional system (6.1). to present the graphical results for the CF fractal-fractional model (2.1) that were obtained by considering the above novel numerical scheme with the initial conditions $\mathscr{X}(0) = 0.2$ , $\mathscr{Y}(0) = 0.2$ , and $\mathscr{Z}(0) = 0.3$ . For different fractal-fractional orders, we give numerical results for the considered model (6.1) that was developed by using two-phase and three-phase simulations. The plots in and shown the dynamical behavior of the $\mathscr{X}(\zeta)$ , $\mathscr{Y}(\zeta)$ , and $\mathscr{Z}(\zeta)$ for different fractional orders in the sense of the CF fractal operator $\mathcal{R}_{1} = \mathcal{R}_{3} = \mathcal{R}_{2} = \mathcal{R}_{5} = \mathcal{R}_{9} = \mathcal{R}_{7} =$ $\mathcal{R}_{10} = \mathcal{R}_{12} = \mathcal{R}_{13} = 20K\Omega, \mathcal{R}_{6} = \mathcal{R}_{8} = 12.62K\Omega, \mathcal{R}_{4} = 10.62K\Omega, \mathcal{R}_{11} = 16.49K \Omega$ . The simulation results for the two- and three-dimensional phase portraits are displayed in the plots of Figures 5 and 7. By decreasing the fractional orders, we can see how the dynamics of the model variables change. We can observe that when the order of the CF fractal-fractional derivative decreases, the dynamics of the system (6.1) decrease as well.

Figure 5. Dynamical behavior of the chaotic attractor for

$\mathscr{X}(1) = 0.2, \mathscr{Y}(1) = 0.2$ and

$\mathscr{Z}(1) = 0.3$ , as obtained for

$\gamma = 1$ and

$\alpha = 0.95$ by using the CF fractal-fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 6. Numerical simulation of (6.2) for

$\gamma = 1$ and different values of

$\alpha = 0.95$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. Dynamical behavior of the chaotic attractor for

$\mathscr{X}(1) = 0.2, \mathscr{Y}(1) = 0.2$ and

$\mathscr{Z}(1) = 0.3$ , as obtained for

$\gamma = 1$ and

$\alpha = 0.85$ by using the CF fractal-fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 8. Numerical simulation of (6.2) for

$\gamma = 0.95$ and different values of

$\alpha = 0.85$ .

DownLoad: Full-Size Img PowerPoint

Two newly proposed fractal-fractional operators based on the Mittag-Leffler and exponential functions were applied to extend the memristor-based chaotic system. The one based on the exponential function was developed by Caputo and Fabrizio, while Atangana and Baleanu proposed the one based on the Mittag-Leffler function. We evaluated the model's more complex behavior in the form of a fractal-fractional operator, which is often challenging to achieve by using the integer-order operator. The fractal-fractional operator is thus a better technique for examining the more complicated behavior of the proposed system.

7. Conclusion and discussion

This article has studied the fractal-fractional memristor-based chaotic model by using different types of fractional derivatives. We established results for the uniqueness and existence of the solution of the model. In each case, a numerical scheme has been established to find approximate solutions of the model, which accurately assess the numerical behavior of the system. Finally, numerical simulations were demonstrated. We defined the sensitivity of these fractal-fractional systems by using the results obtained by modifying some parameters in the model. We have seen that as the values of $\alpha$ decreases, the value of $\mathscr{X}(\zeta)$ decreases significantly for the ABC fractal operator as compared to the CF fractal operator, while the value of $\mathscr{Y}(\zeta)$ and $\mathscr{Z}(\zeta)$ increases significantly for the ABC fractal operator as compared to the CF fractal operator. We can find approximate solutions to some other fractal-fractional systems by using these developed numerical schemes.

Conflict of interest

All authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	M. V. Aarset, How to identify bathtub hazard rate, IEEE Trans. Reliab., 36 (1987), 106–108.
[2]	M. A. Abdullah, H. A. Abdullah, Estimation of Lomax parameters based on generalized probability weighted moment, J. King Saud Univ. Sci., 22 (2010), 171–184.
[3]	A. Afaq, S. P. Ahmad, A. Ahmed, Bayesian Analysis of shape parameter of Lomax distribution under different loss functions, Int. J. Stat. Math., 2 (2010), 55–65.
[4]	M. Ahsanullah, Record values of Lomax distribution, Stat. Neerl., 41 (1991), 21–29.
[5]	A. A. Al-Babtain, I. Elbatal, C. Chesneau, F. Jamal, The transmuted Muth generated class of distributions with applications, Symmetry, 12 (2020), 1677. doi: 10.3390/sym12101677
[6]	A. A. Al-Babtain, I. Elbatal, C. Chesneau, F. Jamal, Box-Cox gamma-G family of distributions: Theory and applications, Mathematics, 8 (2020), 1801. doi: 10.3390/math8101793
[7]	S. Al-Marzouki, F. Jamal, C. Chesneau, M. Elgarhy, Topp-Leone odd Fréchet generated family of distributions with applications to Covid-19 data sets, Comput. Model Eng. Sci., 125 (2020), 437–458.
[8]	M. A. Aldahlan, F. Jamal, C. Chesneau, I. Elbatal, M. Elgarhy, Exponentiated power generalized Weibull power series family of distributions: Properties, estimation and applications, PloS One, 15 (2020), 1–25.
[9]	M. A. Aldahlan, F. Jamal, C. Chesneau, M. Elgarhy, I. Elbatal, The truncated Cauchy power family of distributions with inference and applications, Entropy, 22 (2020), 1–25.
[10]	A. M. Almarashi, M. Elgarhy, F. Jamal, C. Chesneau, The exponentiated truncated inverse Weibull generated family of distributions with applications, Symmetry, 12 (2020), 1–21.
[11]	J. M. Amigo, S. G. Balogh, S. Hernandez, A brief review of generalized entropies, Entropy, 20 (2018), 813. doi: 10.3390/e20110813
[12]	M. Aslam, Z. Asghar, Z. Hussain, S. F. Shah, A modified T-X family of distributions: classical and Bayesian analysis, J. Taibah Univ. Sci., 14 (2020), 254–264. doi: 10.1080/16583655.2020.1732642
[13]	M. A. Badr, I. Elbatal, F. Jamal, C. Chesneau, M. Elgarhy, The transmuted odd Fréchet-G family of distributions: Theory and applications, Mathematics, 8 (2020), 1–20.
[14]	N. Balakrishnan, M. Ahsanullah, Relations for single and product moments of record values from Lomax distribution, Sankhya B, 56 (1994), 140–146.
[15]	A. Balkema, L. de Haan, Residual life time at great age, Ann. Probab., 2 (1974), 792–804.
[16]	R. A. R. Bantan, C. Chesneau, F. Jamal, M. Elgarhy, On the analysis of new Covid-19 cases in Pakistan using an exponentiated version of the M family of distributions, Mathematics, 8 (2020), 1–20.
[17]	R. A. R. Bantan, F. Jamal, C. Chesneau, M. Elgarhy, Type II Power Topp-Leone generated family of distributions with statistical inference and applications, Symmetry, 12 (2020), 1–24.
[18]	M. A. Beg, On the estimation of $P(Y <X)$ for the two parameter exponential distribution, Metrika, 27 (1980), 29–34. doi: 10.1007/BF01893574
[19]	C. C. R. Brito, L. C. Rêgo, W. R. Oliveira, F. Gomes-Silva, Method for generating distributions and classes of probability distributions: The univariate case. Hacettepe J. Math. Stat., 48 (2019), 897–930.
[20]	M. C. Bryson, Heavy-tailed distributions: Properties and tests, Technometrics, 16 (1974), 61–68. doi: 10.1080/00401706.1974.10489150
[21]	G. Casella, R. L. Berger, Statistical Inference, Brooks/Cole Publishing Company: Bel Air, CA, USA, 1990.
[22]	C. Chesneau, T. El Achi, Modified odd Weibull family of distributions: Properties and applications, J. Indian Soc. Prob. Stat., 21 (2020), 259–286. doi: 10.1007/s41096-020-00075-x
[23]	G. M. Cordeiro, R. B. Silva, A. D. C. Nascimento, Recent Advances in Lifetime and Reliability Models, Bentham Sciences Publishers, Sharjah, UAE, 2020.
[24]	A. El-Bassiouny, N. Abdo, H. Shahen, Exponential Lomax distribution, Int. J. Comput. Appl., 121 (2015), 24–29.
[25]	W. Gilchrist, Statistical Modelling with Quantile Functions, CRC Press, Abingdon, 2000.
[26]	M. A. Haq, M. Elgarhy, The odd Fréchet-G family of probability distributions, J. Stat. Appl. Prob., 7 (2018), 189–203. doi: 10.18576/jsap/070117
[27]	F. Jamal, C. Chesneau, M. Elgarhy, Type II general inverse exponential family of distributions, J. Stat. Manage. Syst., 23 (2020), 617–641.
[28]	D. Kumar, U. Singh, S. K. Singh, Life time distribution: Derived from some minimum guarantee distribution, Sohag J. Math., 4 (2017), 7–11. doi: 10.18576/sjm/040102
[29]	A. J. Lemonte, G. M. Cordeiro, An extended Lomax distribution, Statistics, 47 (2013), 800–816. doi: 10.1080/02331888.2011.568119
[30]	M. Mead, A new generalization of Burr XII distribution, J. Stat.: Adv. Theo. Appl., 12 (2014), 53–71.
[31]	N. U. Nair, P. G. Sankaran, Quantile based reliability analysis, Commun. Stat. Theory Methods, 38 (2009), 222–232. doi: 10.1080/03610920802187430
[32]	M. A. Nasir, M. H. Tahir, C. Chesneau, F. Jamal, M. A. A. Shah, The odd generalized gamma-G family of distributions: Properties, regressions and applications, Statistica, 80 (2020), 3–38.
[33]	S. Nasiru, Extended odd Fréchet-G family of distributions, J. Probab. Stat., 1 (2018), 1–12.
[34]	D. Bates, J. Chambers, P. Dalgaard, S. Falcon, R Development Core Team, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0, Available from: http://www.R-project.org, 2005.
[35]	E. A. Rady, W. A. Hassanein, T. A. Elhaddad, The power Lomax distribution with an application to bladder cancer data, SpringerPlus, 5 (2016), 1838. doi: 10.1186/s40064-016-3464-y
[36]	A. Rényi, On measures of entropy and information, Proc. 4th Berkeley symposium on mathematical statistics and probability, 1 (1960), 47–561.
[37]	M. Shaked, J. G. Shanthikumar, Stochastic Orders, Wiley, New York, N.Y., USA, 2007.
[38]	M. H. Tahir, G. M. Cordeiro, Compounding of distributions: A survey and new generalized classes, J. Stat. Distrib. Appl., 3 (2016), 1–35.
[39]	C. Tsallis, Possible generalization of Boltzmann-Gibbs statistics, J. Stat. Phys., 52 (1988), 479–487. doi: 10.1007/BF01016429
[40]	M. A. J. Van Montfort, On testing that the distribution of extremes is of type I when type II is the alternative, J. Hydrol., 11 (1970), 421–427. doi: 10.1016/0022-1694(70)90006-5
[41]	H. M. Yousof, M. Rasekhi, M. Alizadeh, G. G. Hamedani, The Marshall-Olkin exponentiated generalized G family of distributions: properties, applications and characterizations, J. Nonlinear Sci. Appl., 13 (2020), 34–52.
[42]	R. A. ZeinEldin, C. Chesneau, F. Jamal, M. Elgarhy, A. M. Almarashi, S. Al-Marzouki, Truncated generalized Fréchet generated family of distributions with applications, Comput. Model Eng. Sci., 126 (2021), 791–819.

This article has been cited by:

A. E. Matouk, Chaos and hidden chaos in a 4D dynamical system using the fractal-fractional operators, 2025, 10, 2473-6988, 6233, 10.3934/math.2025284

Reader Comments

Your name:*

Email:*
© 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)