
The two-dimensional (2D) cine cardiovascular magnetic resonance (CMR) technique is the reference standard for assessing cardiac function. However, one challenge with 2D cine is that the acquisition time for the whole cine stack is long and requires multiple breath holds, which may not be feasible for pediatric or ill patients. Though single breath-hold multi-slice cine may address the issue, it can only acquire low-resolution images, and hence, affect the accuracy of cardiac function assessment. To address these challenges, a Ferumoxytol-enhanced, free breathing, isotropic high-resolution 3D cine technique was developed. The method produces high-contrast cine images with short acquisition times by using compressed sensing together with a manifold-based method for image denoising. This study included fifteen patients (9.1 ± 5.6 yrs.) who were referred for clinical cardiovascular magnetic resonance imaging (MRI) with Ferumoxytol contrast and were prescribed the 3D cine sequence. The data was acquired on a 1.5T scanner. Statistical analysis shows that the manifold-based denoised 3D cine can accurately measure ventricular function with no significant differences when compared to the conventional 2D breath-hold (BH) cine. The multiplanar reconstructed images of the proposed 3D cine method are visually comparable to the golden standard 2D BH cine method in terms of clarity, contrast, and anatomical precision. The proposed method eliminated the need for breath holds, reduced scan times, enabled multiplanar reconstruction within an isotropic data set, and has the potential to be used as an effective tool to access cardiovascular conditions.
Citation: Anna Andrews, Pezad Doctor, Lasya Gaur, F. Gerald Greil, Tarique Hussain, Qing Zou. Manifold-based denoising for Ferumoxytol-enhanced 3D cardiac cine MRI[J]. Mathematical Biosciences and Engineering, 2024, 21(3): 3695-3712. doi: 10.3934/mbe.2024163
[1] | Anastacia Dlamini, Emile F. Doungmo Goufo, Melusi Khumalo . On the Caputo-Fabrizio fractal fractional representation for the Lorenz chaotic system. AIMS Mathematics, 2021, 6(11): 12395-12421. doi: 10.3934/math.2021717 |
[2] | A. E. Matouk . Chaos and hidden chaos in a 4D dynamical system using the fractal-fractional operators. AIMS Mathematics, 2025, 10(3): 6233-6257. doi: 10.3934/math.2025284 |
[3] | Ihtisham Ul Haq, Shabir Ahmad, Sayed Saifullah, Kamsing Nonlaopon, Ali Akgül . Analysis of fractal fractional Lorenz type and financial chaotic systems with exponential decay kernels. AIMS Mathematics, 2022, 7(10): 18809-18823. doi: 10.3934/math.20221035 |
[4] | Xiaolei Dong, Yuming Qin . Finite fractal dimension of pullback attractors for a nonclassical diffusion equation. AIMS Mathematics, 2022, 7(5): 8064-8079. doi: 10.3934/math.2022449 |
[5] | Muhammad Farman, Aqeel Ahmad, Ali Akgül, Muhammad Umer Saleem, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Dynamical behavior of tumor-immune system with fractal-fractional operator. AIMS Mathematics, 2022, 7(5): 8751-8773. doi: 10.3934/math.2022489 |
[6] | Aziz Khan, Thabet Abdeljawad, Manar A. Alqudah . Neural networking study of worms in a wireless sensor model in the sense of fractal fractional. AIMS Mathematics, 2023, 8(11): 26406-26424. doi: 10.3934/math.20231348 |
[7] | Muhammad Sarwar, Aiman Mukheimer, Syed Khayyam Shah, Arshad Khan . Existence of solutions of fractal fractional partial differential equations through different contractions. AIMS Mathematics, 2024, 9(5): 12399-12411. doi: 10.3934/math.2024606 |
[8] | Narcisse Batangouna . A robust family of exponential attractors for a time semi-discretization of the Ginzburg-Landau equation. AIMS Mathematics, 2022, 7(1): 1399-1415. doi: 10.3934/math.2022082 |
[9] | Amir Ali, Abid Ullah Khan, Obaid Algahtani, Sayed Saifullah . Semi-analytical and numerical computation of fractal-fractional sine-Gordon equation with non-singular kernels. AIMS Mathematics, 2022, 7(8): 14975-14990. doi: 10.3934/math.2022820 |
[10] | Daniel Clemente-López, Esteban Tlelo-Cuautle, Luis-Gerardo de la Fraga, José de Jesús Rangel-Magdaleno, Jesus Manuel Munoz-Pacheco . Poincaré maps for detecting chaos in fractional-order systems with hidden attractors for its Kaplan-Yorke dimension optimization. AIMS Mathematics, 2022, 7(4): 5871-5894. doi: 10.3934/math.2022326 |
The two-dimensional (2D) cine cardiovascular magnetic resonance (CMR) technique is the reference standard for assessing cardiac function. However, one challenge with 2D cine is that the acquisition time for the whole cine stack is long and requires multiple breath holds, which may not be feasible for pediatric or ill patients. Though single breath-hold multi-slice cine may address the issue, it can only acquire low-resolution images, and hence, affect the accuracy of cardiac function assessment. To address these challenges, a Ferumoxytol-enhanced, free breathing, isotropic high-resolution 3D cine technique was developed. The method produces high-contrast cine images with short acquisition times by using compressed sensing together with a manifold-based method for image denoising. This study included fifteen patients (9.1 ± 5.6 yrs.) who were referred for clinical cardiovascular magnetic resonance imaging (MRI) with Ferumoxytol contrast and were prescribed the 3D cine sequence. The data was acquired on a 1.5T scanner. Statistical analysis shows that the manifold-based denoised 3D cine can accurately measure ventricular function with no significant differences when compared to the conventional 2D breath-hold (BH) cine. The multiplanar reconstructed images of the proposed 3D cine method are visually comparable to the golden standard 2D BH cine method in terms of clarity, contrast, and anatomical precision. The proposed method eliminated the need for breath holds, reduced scan times, enabled multiplanar reconstruction within an isotropic data set, and has the potential to be used as an effective tool to access cardiovascular conditions.
Very recently the concept of fractal differentiation and fractional differentiation have been combined to produce new differentiation operators [1,2,3,4]. The new operators were constructed using three different kernels namely, power law, exponential decay and the generalised Mittag-Leffler function. The new operators have two parameters, the first is considered as fractional order and the second as fractal dimension.
The operators were tested to model some real world problems [9,10,11] surprisingly, the new differential and integral operators were found to be powerful mathematical tools able to capture even hidden complexities of nature [5,6,16]. Very detailed attractors could be captured when modelling with such differential and integral operators [8]. One of the great advantage of these new operators is that they can at the same time capture processes following power law, exponential decay law, fading memory, crossover and self-similar processes. This unique and outstanding capacity makes these new operators suitable tools to model many complex real world problems. Additionally when the fractal dimension is one, we recover all the existing fractional differential and integral operators, if the fractional order is one, we recover the so-call fractal differential and integral operators. Finally if the fractal and fractional orders tend to 1 we recover classical differential and integral operators. Therefore, it is believed that these new operators are the future for modeling.
On the other hand, modeling real-life problems have been always of interest of researchers. Mathematicians especially have introduced and developed a variety of well-known and sophisticated tools to model problems arising from other sciences such as biology, finance, epidemiology, and geology. The list is not exhaustive. Moreover, they used these tools to investigate solutions to many important and crucial problems that humans are challenging in their daily life. Among these tools are differential equations and fractal differential equations for instance. Nevertheless, the continuous improvement leads constantly their researches in order to obtain better solutions that aim at reaching the reality. Fractal-fractional differential operators can be seen as an enhanced instrument for modeling important practical problems. One of the main motivations for this article was to introduce methods for solving fractal-fractional problems that have not been considered before. By using the proposed methods, approximate solutions to these problems can be determined more easily and with higher accuracy. In addition, the methods can be applied to real-world problems. In this work, we present applications of such numerical schemes in solving chaotic models that involve the new class of differential operators. We consider some well-known chaotic models with fractal-fractional differential operators, so that we can compare our results with those in the literature. The article is organized as follows: In section 2, we give a brief overview of some basic definitions of fractal-fractional differential calculus. Two efficient and effective numerical methods for determining approximate solutions to fractal-fractional problems are presented in section 3. The first is for the Caputo fractal derivative second is for Caputo-Fabrizio-Caputo fractal-fractional derivative and the thierd is for the Atangana-Baleanu-Caputo fractal-fractional derivative. The kernels used in these two definitions are singular and non-singular, respectively. Several numerical simulations for chaotic systems are described in section 4 we find the numerical analysis for Duffing attractor. Numerical analysis for EL-Niˆno Southern Oscillation is find in section 5. In section 6–9, we analysis numerical analysis for Ikeda system, numerical analysis for Dadras attractor and numerical analysis for Aizawa attractor, numerical analysis for Thomas attractor and numerical analysis for 4 Wings attractor. In section 9 we introduce a new Atangana Sonal attractor and find the solution of this. The results obtained are accurate, interesting, and meaningful. Finally, we present our overall conclusions.
Definition 2.1. [3] The fractal fractional derivative of f(t) with order ϖ−κ in the Riemann-Liouville sense is defined as follows:
FF−RLDκ,ϖ0,t{f(t)}=1Γ(m−κ)ddtϖ∫t0(t−s)m−κ−1f(s)ds, | (2.1) |
where m−1<κ,ϖ≤m∈N and df(s)dsϖ=limt→sf(t)−f(s)tϖ−sϖ.
Definition 2.2. [3] The Caputo-Fabrizio fractal-fractional derivative of f(t) with order ϖ−κ in the Riemann-Liouville sense is defined as follows:
FF−CFRDκ,ϖ0,t{f(t)}=M(κ)(1−κ)ddtϖ∫t0exp(−κ1−κ(t−s))f(s)ds, | (2.2) |
where κ>0,ϖ≥m∈N and M(0)=M(1)=1.
Definition 2.3. [2] The Atangan-Baleanu fractal-fractional derivative of f(t) with order ϖ−κ in the Riemann-Liouville sense is defined as follows:
FF−ABRDκ,ϖ0,t{f(t)}=AB(κ)(1−κ)ddtϖ∫t0Eκ(−κ1−κ(t−s)κ)f(s)ds, | (2.3) |
where 0<κ,ϖ≤1 and AB(κ)=1−κ+κΓ(κ).
Definition 2.4. [3] The fractal-fractional derivative of with order ϖ−κ in the Liouville-Caputo sense is defined as follows:
FF−CDκ,ϖ0,t{f(t)}=1Γ(m−κ)ddtϖ∫t0(t−s)m−κ−1(ddsτf(s))ds, | (2.4) |
where m−1<κ,ϖ≤m∈N and df(s)dsϖ=limt→sf(t)−f(s)tϖ−sϖ.
Definition 2.5. [2] The Caputo-Fabrizio fractal-fractional derivative of f(t) with order ϖ−κ in the Liouville-Caputo sense is defined as follows:
FF−CFCDκ,ϖ0,t{f(t)}=M(κ)(1−κ)∫t0exp(−κ1−κ(t−s))(ddsϖf(s))ds, | (2.5) |
where κ>0,ϖ≥m∈N and M(0)=M(1)=1.
Definition 2.6. [4] The Atangana-Baleanu fractal-fractional derivative of f(t) with order ϖ−κ in the Liouville-Caputo sense is defined as follows:
FF−ABCDκ,ϖ0,t{f(t)}=AB(κ)(1−κ)∫t0Eκ(−κ1−κ(t−s)κ)(ddsϖf(s))ds, | (2.6) |
where 0<κ,ϖ≤1 and AB(κ)=1−κ+κΓ(κ).
Definition 2.7. [2] The Liouville-Caputo fractal-fractional integral of f(t) with order κ is defined as follows:
FF−CIκ0,t{f(t)}=ϖΓ(κ)∫t0(t−s)κ−1sϖ−1f(s)ds. | (2.7) |
Definition 2.8. [3] The Caputo-Fabrizio fractal-fractional integral of f(t) with order κ is defined as follows:
FF−CFIκ,ϖ0,t{f(t)}=κϖM(κ)∫t0sκ−1f(s)ds+ϖ(1−κ)tϖ−1M(κ)f(t). | (2.8) |
Definition 2.9. [3] The Atangana-Baleanu fractal-fractional integral of f(t) with order κ is defined as follows:
FF−ABIκ,ϖ0,t{f(t)}=κϖAB(κ)∫t0sκ−1(t−s)κ−1f(s)ds+ϖ(1−κ)tϖ−1AB(κ)f(t). | (2.9) |
In this section, we have given three numerical schemes for Caputo-fractal-fractional, Caputo-Fabrizio-fractal fractional and the Atangana-Baleanu fractal-fractional derivative operators [7].
Consider the following differential equations in the fractal-fractional Liouville-Caputo sense
FF−CDκ,ϖ0,t{u(t)}=f(t,u(t),v(t)). | (3.1) |
Equation (3.1) can be converted to the Volterra case and the numerical scheme of this system using a Caputo-fractal-fractional approach at tn+1 is given by
u(tn+1)=u0+ϖΓ(κ)∫tn+10sϖ−1(tn+1−s)κ−1f(s,u,v)ds. | (3.2) |
We can approximate the above integral to
u(tn+1)=u0+ϖΓ(κ)n∑j=0∫tj+1tjsϖ−1(tn+1−s)κ−1f(s,u,v)ds. | (3.3) |
With in the finite interval [tj,tj+1], we approximate the function sϖ−1,f(t,u,v) using the Lagrangian piecewise interpolation such that
E(s)=s−tj−1tj−tj−1tϖ−1jf(tj,uj,vj)−s−tjtj−tj−1tϖ−1j−1f(tj−1,uj−1,vj−1). | (3.4) |
So we obtain
u(tn+1)=u0+ϖΓ(κ)n∑j=0∫tj+1tjsϖ−1(tn+1−s)κ−1E(s)ds. | (3.5) |
Solving the integral of the right hand side, we obtain the following numerical scheme
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1jf(tj,uj,vj)((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1f(tj−1,uj−1,vj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (3.6) |
Consider the following differential equations in the Caputo-Fabrizio-fractal-fractional derivative
FF−CFCDκ,ϖ0,t{u(t)}=f(t,u(t),v(t)). | (3.7) |
Applying the Caputo-Fabrizio integral, we obtain
u(t)=u0+ϖtϖ−1(1−κ)M(κ)f(t,u,v)+κϖM(κ)∫t0sϖ−1f(s,u,v)ds. | (3.8) |
Here we present the detailed derivation of the numerical scheme. Thus, at tn+1 we have
u(tn+1)=u0+ϖtϖ−1n(1−κ)M(κ)f(tn,un,vn)+κϖM(κ)∫tn+10sϖ−1f(s,u,v)ds. | (3.9) |
Taking the difference between the consecutive terms, we obtain
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)f(tn,un,vn)−ϖtϖ−1n−1(1−κ)M(κ)f(tn−1,un−1,vn−1)+κϖM(κ)∫tn+1tnsϖ−1f(s,u,v)ds. | (3.10) |
Now using the Lagrange polynomial piece-wise interpolation and integrating, we obtain
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)f(tn,un,vn)−ϖtϖ−1n−1(1−κ)M(κ)f(tn−1,un−1,vn−1)+κϖM(κ)×[32(Δt)tϖ−1nf(tn,un,vn)−Δt2tϖ−1n−1f(tn−1,un−1,vn−1)]. | (3.11) |
Consider the following differential equations in the Atangana-Baleanu fractal-fractional derivative in the Liouville-Caputo sense
FF−ABCDκ,ϖ0,t{u(t)}=f(t,u(t),v(t)). | (3.12) |
Applying the Atangana-Baleanu integral, we have
u(t)=u(0)+ϖtϖ−1n(1−κ)AB(κ)f(t,u,v)+κϖAB(κ)Γ(κ)∫t0sϖ−1(t−s)κ−1f(s,u,v)ds. | (3.13) |
At tn+1, we have the following
u(tn+1)=u(0)+κϖAB(κ)Γ(κ)∫tn+10sϖ−1(tn+1−s)κ−1f(s,u,v)ds+ϖtϖ−1n(1−κ)AB(κ)f(tn,un,vn), | (3.14) |
the above system can be expressed as
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)f(tn,un,vn)+κϖAB(κ)Γ(κ)n∑j=0∫tj+1tjsϖ−1(tn+1−s)κ−1f(s,u,v)ds. | (3.15) |
Approximating sϖ−1f(s,u,v) in [tj,tj+1], we have the following numerical scheme
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)f(tn,un,vn)+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1jf(tj,uj,vj)[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1f(tj−1,uj−1,vj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (3.16) |
The Duffing Pendulum is a kind of a forced oscillator with damping, the mathematical model represented in state variable as when we apply Caputo-fractal-fractional derivative is given by
FF−CDκ,ϖ0,t{u1(t)}=Au2(t), | (4.1) |
FF−CDκ,ϖ0,t{u2(t)}=u1(t)−u31(t)−au2(t)+bcos(ωt), | (4.2) |
where initial conditions are u1(0)=0.01, u2(0)=0.5 and the constant are κ=0.25, b=0.3, ω=1, h=0.01, and t=100. The numerical scheme is given by
u1(tn+1)=u01+Aϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1ju2(t)j((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1Au2(t)j−1((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],u2(tn+1)=u01+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(u1(t)j−u31(t)j−au2(tj−τ)+bcosωtj)×((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−(tϖ−1j−1u1(t)j−1−u31(t)j−1−au2(tj−1−τ)+bcosωtj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (4.3) |
Numerical simulation of the Caputo-power law case are depicted in Figures 1–3 for different values of fractional order κ and fractal dimension ϖ.
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u1(t)}=Au2(t), | (4.4) |
FF−CFCDκ,ϖ0,t{u2(t)}=u1(t)−u31(t)−au2(t)+bcosωt, | (4.5) |
The numerical scheme is given by
u1(tn+1)=un1+ϖtϖ−1n(1−κ)M(κ)u2(t)n−ϖtϖ−1n−1(1−κ)M(κ)u2(t)n−1+κϖM(κ)×[32(Δt)tϖ−1nu2(t)n−Δt2tϖ−1n−1u2(t)n−1],u2(tn+1)=un1+ϖtϖ−1n(1−κ)M(κ)(u1(t)n−u31(t)n−au2(tn−τ)+bcosωtn)−ϖtϖ−1n−1(1−κ)M(κ)(u1(t)n−1−u31(t)n−1−au2(tn−1−τ)+bcosωtn−1)+κϖM(κ)×[32(Δt)tϖ−1n(u1(t)n−u31(t)n−au2(tn−τ)+bcosωtn)−Δt2tϖ−1n−1(u1(t)n−1−u31(t)n−1−au2(tn−1−τ)+bcosωtn−1)]. | (4.6) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u1(t)}=Au2(t), | (4.7) |
FF−ABCDκ,ϖ0,t{u2(t)}=u1(t)−u31(t)−au2(t)+bcosωt, | (4.8) |
The differential and integral operators used in the case of Caputo and Caputo-Fabrizio have kernel with no crossover in waiting time distribution. Therefore to include the effect of crossover in waiting time distribution, we make use of the differential operator whose kernel has crossover behavior, thus the Atangana-Baleanu operators are used here. Also the numerical simulation is presented in Figures 4–8.
Additionally adapting the previously used numerical scheme, we obtain the following:
u1(tn+1)=u1(0)+ϖtϖ−1n(1−κ)AB(κ)u2(t)n+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1ju2(t)j[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1u2(t)j−1((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],u2(tn+1)=u2(0)+ϖtϖ−1n(1−κ)AB(κ)(u1(t)n−u31(t)n−au2(tn−τ)+bcosωtn)+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(u1(t)j−u31(t)j−au2(tj−τ)+bcosωtj)[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(u1(t)j−1−u31(t)j−1−au2(tj−1−τ)+bcosωtj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (4.9) |
The EL-Niˆno Southern Oscillation known as ENSO is an unconventional repetitive variation in sea area temperatures over the tropical eastern Pacific Ocean and winds that affect the climate change of much of the tropics and sub-tropics The natural behavior of this dynamical system can be classified in two phases. The first phase is the warming phase of the sea temperature, this phase was named EL-Niˆno. The second phase known as EL-Niˆno is the cooling phase. The results of some data collection showed that, the Southern Oscillation is following atmospheric component, more importantly the phenomena is coupled with the sea temperature change. On the other hand, the EL-Niˆno is followed by high air surface pressure in the side of tropical western Pacific and La Nina is accompanied with low pressure. In this section, we consider the mathematical model able to replicate such natural occurrence, nevertheless, here we consider the model with different non-local differential and integral operators.
The behavior of this phenomena is described by the following equation as in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=−tanh(κu(t−τ))+bcos(2πt), | (5.1) |
where κ=100 and b=1.
The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−tanh(κu(tj−τ))+bcos(2πtj))×((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−tanh(κu(tj−1−τ))+bcos(2πtj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (5.2) |
The above numerical solution is depicted in Figures 9–13, for different value of fractional order and fractal dimension. In this exercise, we chose to simulate the case with Caputo derivative. The numerical simulation depend on the fractional order and the fractal dimension. More precisely for the values 0.8 and 0.5, we observed a very strange behavior where the oscillation are vanishing.
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u(t)}=−tanh(κu(t−τ))+bcos(2πt). | (5.3) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)(−tanh(κu(tn−τ))+bcos(2πtn))−ϖtϖ−1n−1(1−κ)M(κ)(−tanh(κu(tn−1−τ))+bcos(2πtn−1))+κϖM(κ)×[32(Δt)tϖ−1n(−tanh(κu(tn−τ))+bcos(2πtn))−Δt2tϖ−1n−1−tanh(κu(tn−1−τ))+bcos(2πtn−1)]. | (5.4) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=−tanh(κu(t−τ))+bcos(2πt). | (5.5) |
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)(−tanh(κu(tn−τ))+bcos(2πtn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−tanh(κu(tj−τ))+bcos(2πtj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−tanh(κu(tj−1−τ))+bcos(2πtj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (5.6) |
Consider the ikeda delay system in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=−bu(t)+asin(cx(t−τ)). | (6.1) |
In [12], the author proposed a second delay parameter. The below given equation is resultant solution
FF−CDκ,ϖ0,t{u(t)}=−bu(t−τ1)+asin(cu(t−τ2)), | (6.2) |
here we consider a = 24, b = 3 and c = 1 and the delays τ1=0.01 and τ2=0.1.
The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−bu(tj−τ1)+asin(cu(tj−τ2)))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−bu(tj−1−τ1)+asin(cu(tj−1−τ2)))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (6.3) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u1(t)}=−bu(t−τ1)+asin(cu(t−τ2)). | (6.4) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)(−bu(tn−τ1)+asin(cu(tn−τ2)))−ϖtϖ−1n−1(1−κ)M(κ)(−bu(tn−1−τ1)+asin(cu(tn−1−τ2)))+κϖM(κ)×[32(Δt)tϖ−1n(−bu(tn−τ1)+asin(cu(tn−τ2)))−Δt2tϖ−1n−1(−bu(tn−1−τ1)+asin(cu(tn−1−τ2)))]. | (6.5) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=−bu(t−τ1)+asin(cu(t−τ2)). | (6.6) |
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)(−bu(tn−τ1)+asin(cu(tn−τ2)))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−bu(tj−τ1)+asin(cu(tj−τ2)))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−bu(tj−1−τ1)+asin(cu(tj−1−τ2)))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (6.7) |
Consider the Dadras system described [13] by the following equation as in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=w(t)−Au(t)+Bv(t)w(t),FF−CDκ,ϖ0,t{v(t)}=Cv(t)−u(t)w(t)+w(t),FF−CDκ,ϖ0,t{w(t)}=Du(t)v(t)−Ew(t), | (7.1) |
where A = 3, B = 2.7, C = 1.7, D = 2 and E = 9.
The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(w(tj)−Au(tj)+Bv(tj)w(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(w(tj−1)−Au(tj−1)+Bv(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],v(tn+1)=v0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(Cv(tj)−u(tj)w(tj)+w(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Cv(tj−1)−u(tj−1)w(tj−1)+w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],w(tn+1)=w0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(Du(tj)v(tj)−Ew(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Du(tj−1)v(tj−1)−Ew(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (7.2) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u(t)}=w(t)−Au(t)+Bv(t)w(t),FF−CFCDκ,ϖ0,t{v(t)}=Cv(t)−u(t)w(t)+w(t),FF−CFCDκ,ϖ0,t{w(t)}=Du(t)v(t)−Ew(t). | (7.3) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)(w(tn)−Au(tn)+Bv(tn)w(tn))−ϖtϖ−1n−1(1−κ)M(κ)(w(tn−1)−Au(tn−1)+Bv(tn−1)w(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(w(tn)−Au(tn)+Bv(tn)w(tn))−Δt2tϖ−1n−1(w(tn−1)−Au(tn−1)+Bv(tn−1)w(tn−1))], | (7.4) |
v(tn+1)=vn+ϖtϖ−1n(1−κ)M(κ)(Cv(tn)−u(tn)w(tn)+w(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Cv(tn−1)−u(tn−1)w(tn−1)+w(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Cv(tn)−u(tn)w(tn)+w(tn))−Δt2tϖ−1n−1(Cv(tn−1)−u(tn−1)w(tn−1)+w(tn−1))], | (7.5) |
w(tn+1)=wn+ϖtϖ−1n(1−κ)M(κ)(Du(tn)v(tn)−Ew(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Du(tn−1)v(tn−1)−Ew(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Du(tn)v(tn)−Ew(tn))−Δt2tϖ−1n−1(Du(tn−1)v(tn−1)−Ew(tn−1))]. | (7.6) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=w(t)−Au(t)+Bv(t)w(t),FF−ABCDκ,ϖ0,t{v(t)}=Cv(t)−u(t)w(t)+w(t),FF−ABCDκ,ϖ0,t{w(t)}=Du(t)v(t)−Ew(t). | (7.7) |
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)(w(tn)−Au(tn)+Bv(tn)w(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(w(tj)−Au(tj)+Bv(tj)w(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(w(tj−1)−Au(tj−1)+Bv(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (7.8) |
v(tn+1)=v(0)+ϖtϖ−1n(1−κ)AB(κ)(Cv(tn)−u(tn)w(tn)+w(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Cv(tj)−u(tj)w(tj)+w(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Cv(tj−1)−u(tj−1)w(tj−1)+w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (7.9) |
w(tn+1)=w(0)+ϖtϖ−1n(1−κ)AB(κ)(Du(tn)v(tn)−Ew(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Du(tj)v(tj)−Ew(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Du(tj−1)v(tj−1)−Ew(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (7.10) |
Numerical simulation are depicted in Figures 14–17 for the value of κ=0.8 and value of ϖ=0.9. The used numerical model is that with the fractal-fractional differential operator. The model with this new differential operator display some important behavior.
A very strange attractor have been studied in the last past years, although such system of equations have not being attracting attention of many researchers, but the attractor is very strange as the system is able to display very interesting attractor in form of sphere. The model under investigation is called Aizawa attractor, this system when applied iteratively on three-dimensional coordinates, it is important to point out that, evolving in such a way as to have the consequential synchronizes map out a three dimensional shape, in this case a sphere with a tube-like structure powerful one of it's axis. In this section, we consider the model using the fractal-fractional with power law, exponential decay law and the generalized Mittag-Leffler function.
Consider the Aizawa system described by the following equation as in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=(w(t)−B)u(t)−Dv(t), | (8.1) |
FF−CDκ,ϖ0,t{v(t)}=Du(t)+(w(t)−B)v(t), | (8.2) |
FF−CDκ,ϖ0,t{w(t)}=C+Aw(t)−13w3(t)−(u2(t)+v2(t))(1+Ew(t))+Fw(t)u3(t), | (8.3) |
where A=0.95, B=0.7, C=0.6, D=3.5, E=0.25, and F=0.1. The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j((w(tj)−B)u(tj)−Dv(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1((w(tj−1)−B)u(tj−1)−Dv(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (8.4) |
v(tn+1)=v0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(Du(tj)+(w(tj)−B)v(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Du(tj−1)+(w(tj−1)−B)v(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (8.5) |
w(tn+1)=w0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(C+Aw(tj)−13w3(tj)−(u2(tj)+v2(tj))(1+Ew(tj))+Fw(tj)u3(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(C+Aw(tj−1)−13w3(tj−1)−(u2(tj−1)+v2(tj−1))(1+Ew(tj−1))+Fw(tj−1)u3(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (8.6) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u(t)}=(w(t)−B)u(t)−Dv(t), | (8.7) |
FF−CFCDκ,ϖ0,t{v(t)}=Du(t)+(w(t)−B)v(t), | (8.8) |
FF−CFCDκ,ϖ0,t{w(t)}=C+Aw(t)−13w3(t)−(u2(t)+v2(t))(1+Ew(t))+Fw(t)u3(t). | (8.9) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)((w(tn)−B)u(tn)−Dv(tn))−ϖtϖ−1n−1(1−κ)M(κ)((w(tn−1)−B)u(tn−1)−Dv(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n((w(tn)−B)u(tn)−Dv(tn))−Δt2tϖ−1n−1((w(tn−1)−B)u(tn−1)−Dv(tn−1))], | (8.10) |
v(tn+1)=vn+ϖtϖ−1n(1−κ)M(κ)(Du(tn)+(w(tn)−B)v(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Du(tn−1)+(w(tn−1)−B)v(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Du(tn)+(w(tn)−B)v(tn))−Δt2tϖ−1n−1(Du(tn−1)+(w(tn−1)−B)v(tn−1))], | (8.11) |
w(tn+1)=wn+ϖtϖ−1n(1−κ)M(κ)(C+Aw(tn)−13w3(tn)−(u2(tn)+v2(tn))(1+Ew(tn))+Fw(tn)u3(tn))−ϖtϖ−1n−1(1−κ)M(κ)(C+Aw(tn−1)−13w3(tn−1)−(u2(tn−1)+v2(tn−1))×(1+Ew(tn−1))+Fw(tn−1)u3(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(C+Aw(tn)−13w3(tn)−(u2(tn)+v2(tn))×(1+Ew(tn))+Fw(tn)u3(tn))−Δt2tϖ−1n−1(C+Aw(tn−1)−13w3(tn−1)−(u2(tn−1)+v2(tn−1))×(1+Ew(tn−1))+Fw(tn−1)u3(tn−1))]. | (8.12) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=(w(t)−B)u(t)−Dv(t), | (8.13) |
FF−ABCDκ,ϖ0,t{v(t)}=Du(t)+(w(t)−B)v(t), | (8.14) |
FF−ABCDκ,ϖ0,t{w(t)}=C+Aw(t)−13w3(t)−(u2(t)+v2(t))(1+Ew(t))+Fw(t)u3(t). | (8.15) |
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)((w(tn)−B)u(tn)−Dv(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j((w(tj)−B)u(tj)−Dv(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(w(tj−1)−B)u(tj−1)−Dv(tj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (8.16) |
v(tn+1)=v(0)+ϖtϖ−1n(1−κ)AB(κ)(Du(tn)+(w(tn)−B)v(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Du(tj)+(w(tj)−B)v(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Du(tj−1)+(w(tj−1)−B)v(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (8.17) |
w(tn+1)=w(0)+ϖtϖ−1n(1−κ)AB(κ)(C+Aw(tn)−13w3(tn)−(u2(tn)+v2(tn))(1+Ew(tn))+Fw(tn)u3(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(C+Aw(tj)−13w3(tj)−(u2(tj)+v2(tj))(1+Ew(tj))+Fw(tj)u3(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(C+Aw(tj−1)−13w3(tj−1)−(u2(tj−1)+v2(tj−1))(1+Ew(tj−1))+Fw(tj−1)u3(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (8.18) |
Numerical simulation are depicted in Figures 18–21 for the value of κ=0.8 and value of ϖ=0.9. The used numerical model is that with the fractal-fractional differential operator. The model with this new differential operator display some important behavior.
Consider the Thomas system described by the following equation as in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=−Au(t)+sinv(t), | (9.1) |
FF−CDκ,ϖ0,t{v(t)}=−Av(t)+sinw(t), | (9.2) |
FF−CDκ,ϖ0,t{w(t)}=−Aw(t)+sinu(t). | (9.3) |
The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−Au(tj)+sinv(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−Au(tj−1)+sinv(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (9.4) |
v(tn+1)=v0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−Av(tj)+sinw(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−Av(tj−1)+sinw(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (9.5) |
w(tn+1)=w0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−Aw(tj)+sinu(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−Aw(tj−1)+sinu(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (9.6) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u(t)}=−Au(t)+sinv(t), | (9.7) |
FF−CFCDκ,ϖ0,t{v(t)}=−Av(t)+sinw(t), | (9.8) |
FF−CFCDκ,ϖ0,t{w(t)}=−Aw(t)+sinu(t). | (9.9) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)(−Au(tn)+sinv(tn))−ϖtϖ−1n−1(1−κ)M(κ)(−Au(tn−1)+sinv(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(−Au(tn)+sinv(tn))−Δt2tϖ−1n−1(−Au(tn−1)+sinv(tn−1))], | (9.10) |
v(tn+1)=vn+ϖtϖ−1n(1−κ)M(κ)(−Av(tn)+sinw(tn))−ϖtϖ−1n−1(1−κ)M(κ)(−Av(tn−1)+sinw(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(−Av(tn)+sinw(tn))−Δt2tϖ−1n−1(−Av(tn−1)+sinw(tn−1))], | (9.11) |
w(tn+1)=wn+ϖtϖ−1n(1−κ)M(κ)(−Aw(tn)+sinu(tn))−ϖtϖ−1n−1(1−κ)M(κ)(−Aw(tn−1)+sinu(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(−Aw(tn)+sinu(tn))−Δt2tϖ−1n−1(−Aw(tn−1)+sinu(tn−1))]. | (9.12) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=−Au(t)+sinv(t), | (9.13) |
FF−ABCDκ,ϖ0,t{v(t)}=−Av(t)+sinw(t), | (9.14) |
FF−ABCDκ,ϖ0,t{w(t)}=−Aw(t)+sinu(t), | (9.15) |
where A=0.19.
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)(−Au(tn)+sinv(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−Au(tj)+sinv(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−Au(tj−1)+sinv(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (9.16) |
v(tn+1)=v(0)+ϖtϖ−1n(1−κ)AB(κ)(−Av(tn)+sinw(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−Av(tj)+sinw(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−Av(tj−1)+sinw(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (9.17) |
w(tn+1)=w(0)+ϖtϖ−1n(1−κ)AB(κ)(−Aw(tn)+sinu(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−Aw(tj)+sinu(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−Aw(tj−1)+sinu(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (9.18) |
Numerical simulation are depicted in Figures 22–25 for the value of κ=0.9 and value of ϖ=0.9. The used numerical model is that with the fractal-fractional differential operator. The model with this new differential operator display some important behavior.
Consider the 4 Wings system proposed in [15] and modifying in [14] is described by the following equation as in Caputo-fractal-fractional derivative sense
FF−CDκ,ϖ0,t{u(t)}=Au(t)−Bv(t)w(t), | (10.1) |
FF−CDκ,ϖ0,t{v(t)}=−Cv(t)+u(t)w(t), | (10.2) |
FF−CDκ,ϖ0,t{w(t)}=Du(t)−Ew(t)+u(t)v(t), | (10.3) |
where A=4, B=6, C=10, D=5, and E=1.
The numerical scheme is given by
u(tn+1)=u0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(Au(tj)−Bv(tj)w(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Au(tj−1)−Bv(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (10.4) |
v(tn+1)=v0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(−Cv(tj)+u(tj)w(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(−Cv(tj−1)+u(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (10.5) |
w(tn+1)=w0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(Du(tj)−Ew(tj)+u(tj)v(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Du(tj−1)−Ew(tj−1)+u(tj−1)v(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (10.6) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{u(t)}=Au(t)−Bv(t)w(t), | (10.7) |
FF−CFCDκ,ϖ0,t{v(t)}=−Cv(t)+u(t)w(t), | (10.8) |
FF−CFCDκ,ϖ0,t{w(t)}=Du(t)−Ew(t)+u(t)v(t). | (10.9) |
The numerical scheme is given by
u(tn+1)=un+ϖtϖ−1n(1−κ)M(κ)(Au(tn)−Bv(tn)w(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Au(tn−1)−Bv(tn−1)w(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Au(tn)−Bv(tn)w(tn))−Δt2tϖ−1n−1(Au(tn−1)−Bv(tn−1)w(tn−1))], | (10.10) |
v(tn+1)=vn+ϖtϖ−1n(1−κ)M(κ)(−Cv(tn)+u(tn)w(tn))−ϖtϖ−1n−1(1−κ)M(κ)(−Cv(tn−1)+u(tn−1)w(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(−Cv(tn)+u(tn)w(tn))−Δt2tϖ−1n−1(−Cv(tn−1)+u(tn−1)w(tn−1))], | (10.11) |
w(tn+1)=wn+ϖtϖ−1n(1−κ)M(κ)(Du(tn)−Ew(tn)+u(tn)v(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Du(tn−1)−Ew(tn−1)+u(tn−1)v(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Du(tn)−Ew(tn)+u(tn)v(tn))−Δt2tϖ−1n−1(Du(tn−1)−Ew(tn−1)+u(tn−1)v(tn−1))], | (10.12) |
w(tn+1)=wn+ϖtϖ−1n(1−κ)M(κ)(Du(tn)−Ew(tn)+u(tn)v(tn))−ϖtϖ−1n−1(1−κ)M(κ)(Du(tn−1)−Ew(tn−1)+u(tn−1)v(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(Du(tn)−Ew(tn)+u(tn)v(tn))−Δt2tϖ−1n−1(Du(tn−1)−Ew(tn−1)+u(tn−1)v(tn−1))]. | (10.13) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{u(t)}=Au(t)−Bv(t)w(t), | (10.14) |
FF−ABCDκ,ϖ0,t{v(t)}=−Cv(t)+u(t)w(t), | (10.15) |
FF−ABCDκ,ϖ0,t{w(t)}=Du(t)−Ew(t)+u(t)v(t). | (10.16) |
The numerical scheme is given by
u(tn+1)=u(0)+ϖtϖ−1n(1−κ)AB(κ)(Au(tn)−Bv(tn)w(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Au(tj)−Bv(tj)w(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Au(tj−1)−Bv(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (10.17) |
v(tn+1)=v(0)+ϖtϖ−1n(1−κ)AB(κ)(−Cv(tn)+u(tn)w(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(−Cv(tj)+u(tj)w(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(−Cv(tj−1)+u(tj−1)w(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (10.18) |
w(tn+1)=w(0)+ϖtϖ−1n(1−κ)AB(κ)(Du(tn)−Ew(tn)+u(tn)v(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Du(tj)−Ew(tj)+u(tj)v(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Du(tj−1)−Ew(tj−1)+u(tj−1)v(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (10.19) |
Numerical simulation are depicted in Figures 26–29 for the value of κ=0.9 and value of ϖ=0.8. The used numerical model is that with the fractal-fractional differential operator. The model with this new differential operator display some important behavior.
We introduce the following attractor called AS attractor
dx(t)dt=Ay(t),dy(t)dt=x(t)−x3(t)−ay(t)+bcos(ct),dz(t)dt=By(t)+sin(y(t))+x(t), | (11.1) |
where the initial conditions
x(0)=0.01,y(0)=0.5,z(0)=−0.6 |
a=0.25,b=0.3,c=1. |
The mathematical model represented in state variable as when we apply Caputo-fractal-fractional derivative is given by
FF−CDκ,ϖ0,t{x(t)}=Ay(t),FF−CDκ,ϖ0,t{y(t)}=x(t)−x3(t)−ay(t)+bcos(ct),FF−CDκ,ϖ0,t{z(t)}=By(t)+sin(y(t))+x(t). | (11.2) |
The numerical scheme is given by
x(tn+1)=x0+Aϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(y(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(Ay(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],y(tn+1)=y0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(x(tj)−x3(tj)−ay(tj)+bcos(ctj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(x(tj−1)−x3(tj−1)−ay(tj−1)+bcos(ctj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))],z(tn+1)=z0+ϖ(Δt)κΓ(κ+2)n∑j=0[tϖ−1j(By(tj)+sin(y(tj))+x(tj))((n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ))−tϖ−1j−1(By(tj−1)+sin(y(tj−1))+x(tj−1)((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (11.3) |
Involving the fractal-fractional derivative in the Caputo-Fabrizio-Caputo, we have
FF−CFCDκ,ϖ0,t{x(t)}=Ay(t),FF−CFCDκ,ϖ0,t{y(t)}=x(t)−x3(t)−ay(t)+bcos(ct),FF−CFCDκ,ϖ0,t{z(t)}=By(t)+sin(y(t))+x(t). | (11.4) |
x(tn+1)=xn+ϖtϖ−1n(1−κ)M(κ)(Ay(tn)))−ϖtϖ−1n−1(1−κ)M(κ)(Ay(tn−1)))+κϖM(κ)×[32(Δt)tϖ−1n((Ay(tn)))−Δt2tϖ−1n−1(Ay(tn−1))], | (11.5) |
y(tn+1)=yn+ϖtϖ−1n(1−κ)M(κ)(x(tn)−x3(tn)−ay(tn)+bcos(ctn))−ϖtϖ−1n−1(1−κ)M(κ)(x(tn−1)−x3(tn−1)−ay(tn−1)+bcos(ctn−1))+κϖM(κ)×[32(Δt)tϖ−1n(x(tn)−x3(tn)−ay(tn)+bcos(ctn))−Δt2tϖ−1n−1(x(tn−1)−x3(tn−1)−ay(tn−1)+bcos(ctn−1))], | (11.6) |
z(tn+1)=zn+ϖtϖ−1n(1−κ)M(κ)(By(tn)+sin(y(tn))+x(tn))−ϖtϖ−1n−1(1−κ)M(κ)(By(tn−1)+sin(y(tn−1))+x(tn−1))+κϖM(κ)×[32(Δt)tϖ−1n(By(tn)+sin(y(tn))+x(tn))−Δt2tϖ−1n−1(By(tn−1)+sin(y(tn−1))+x(tn−1))]. | (11.7) |
Involving the fractal-fractional derivative in the Atangana-Baleanu-Caputo sense, we have
FF−ABCDκ,ϖ0,t{x(t)}=Ay(t),FF−ABCDκ,ϖ0,t{y(t)}=x(t)−x3(t)−ay(t)+bcos(ct),FF−ABCDκ,ϖ0,t{z(t)}=By(t)+sin(y(t))+x(t). | (11.8) |
The numerical scheme is given by
x(tn+1)=x(0)+ϖtϖ−1n(1−κ)AB(κ)(Ay(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(Ay(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(Ay(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (11.9) |
y(tn+1)=y(0)+ϖtϖ−1n(1−κ)AB(κ)(x(tn)−x3(tn)−ay(tn)+bcos(ctn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(x(tj)−x3(tj)−ay(tj)+bcos(ctj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(x(tj−1)−x3(tj−1)−ay(tj−1)+bcos(ctj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))], | (11.10) |
z(tn+1)=z(0)+ϖtϖ−1n(1−κ)AB(κ)(By(tn)+sin(y(tn))+x(tn))+ϖΔtκAB(κ)Γ(κ+2)×n∑j=0[tϖ−1j(By(tj)+sin(y(tj))+x(tj))[(n+1−j)κ(n−j+2+κ)−(n−j)κ(n−j+2+2κ)]−tϖ−1j−1(By(tj−1)+sin(y(tj−1))+x(tj−1))((n+1−j)κ+1−(n−j)κ(n−j+1+κ))]. | (11.11) |
We present the numerical simulation in the following Figures 30–42 for different values of fractional order and fractal dimension.
Fractal-fractional differential operators have been introduced very recently; however, the new concept has not yet attracted attention of many scholars. In fact, few works have being done where such differential and integral operators are used. The advantages if this work is that it considers operators capturing new complexities of nature with great success. Such differential and integral operators are sophisticated tools to model complex real world problems. Moreover, in order to evaluate the efficiency and the capabilities of the new differential and integral operators, we investigate the behavior of some well-known chaotic attractors and see if one will capture more complexities compared to existing differential and integral operators. Additionally, we introduced a new chaotic model with alternative attractors and we showed by numerical simulation that these new differential and integral operators are powerful mathematical operators able to capture heterogeneity. Since any new numerical method should be validated in terms of convergence, stability and consistency of solutions, these are important research directions left to future work.
The authors would like to express their sincere appreciation to the United Arab Emirates University for the financial support through UPAR Grant No. 31S369.
The authors have declared no conflict of interest.
[1] |
S. V. Raman, M. Markl, A. R. Patel, J. Bryant, B. D. Allen, S. Plein, et al., 30-minute cmr for common clinical indications: a society for cardiovascular magnetic resonance white paper, J. Cardiovasc. Magn. Reson., 24 (2022), 13. https://doi.org/10.1186/s12968-022-00844-6 doi: 10.1186/s12968-022-00844-6
![]() |
[2] | M. A. Syed, S. V. Raman, O. P. Simonetti, Basic principles of cardiovascular MRI: physics and imaging techniques, Springer, 2015. |
[3] |
R. Menchón-Lara, F. Simmross-Wattenberg, P. Higuera, M. Martín-Fernández, C. Alberola-López, Reconstruction techniques for cardiac cine mri, Insights Imag., 10 (2019), 1–16. https://doi.org/10.1186/s13244-019-0754-2 doi: 10.1186/s13244-019-0754-2
![]() |
[4] |
C. M. Kramer, J. Barkhausen, C. Bucciarelli-Ducci, S. D. Flamm, R. J. Kim, E. Nagel, Standardized cardiovascular magnetic resonance imaging (cmr) protocols: 2020 update, J. Cardiovasc. Magn. Reson., 22 (2020), 1–18. https://doi.org/10.1186/s12968-020-00607-1 doi: 10.1186/s12968-020-00607-1
![]() |
[5] |
J. Liu, P. Spincemaille, N. C. F. Codella, T. D. Nguyen, M. R. Prince, Y. Wang, Respiratory and cardiac self-gated free-breathing cardiac cine imaging with multiecho 3d hybrid radial ssfp acquisition, Magn. Reson. Med., 63 (2010), 1230–1237. https://doi.org/10.1002/mrm.22306 doi: 10.1002/mrm.22306
![]() |
[6] |
Q. Lyu, H. Shan, Y. Xie, A. C. Kwan, Y. Otaki, K. Kuronuma, et al., Cine cardiac mri motion artifact reduction using a recurrent neural network, IEEE Trans. Med. Imaging, 40 (2021), 2170–2181. https://doi.org/10.1109/TMI.2021.3073381 doi: 10.1109/TMI.2021.3073381
![]() |
[7] |
Q. Zou, A. H. Ahmed, P. Nagpal, S. Kruger, M. Jacob, Dynamic imaging using a deep generative storm (gen-storm) model, IEEE Trans. Med. Imaging, 40 (2021), 3102–3112. https://doi.org/10.1109/TMI.2021.3065948 doi: 10.1109/TMI.2021.3065948
![]() |
[8] |
Q. Zou, A. H. Ahmed, S. Dzelebdzic, T. Hussain, Free-breathing and ungated cardiac mri reconstruction using a deep kernel representation, Appl. Sci., 13 (2023), 2281. https://doi.org/10.3390/app13042281 doi: 10.3390/app13042281
![]() |
[9] |
D. C. Peters, R. Nezafat, H. Eggers, C. Stehning, W. J. Manning, 2d free-breathing dual navigator-gated cardiac function validated against the 2d breath-hold acquisition, J. Magn. Reson. Imaging, 28 (2008), 773–777. https://doi.org/10.1002/jmri.21417 doi: 10.1002/jmri.21417
![]() |
[10] |
M. H. Moghari, A. Barthur, M. E. Amaral, T. Geva, A. J. Powell, Free-breathing whole-heart 3d cine magnetic resonance imaging with prospective respiratory motion compensation, Magn. Reson. Med., 80 (2018), 181–189. https://doi.org/10.1002/mrm.27021 doi: 10.1002/mrm.27021
![]() |
[11] |
M. Usman, B. Ruijsink, M. S. Nazir, G. Cruz, C. Prieto, Free breathing whole-heart 3d cine mri with self-gated cartesian trajectory, Magn. Reson. Med., 38 (2017), 129–137. https://doi.org/10.1016/j.mri.2016.12.021 doi: 10.1016/j.mri.2016.12.021
![]() |
[12] |
J. Liu, L. Feng, H. -W. Shen, C. Zhu, Y. Wang, K. Mukai, et al., Highly-accelerated self-gated free-breathing 3d cardiac cine mri: Validation in assessment of left ventricular function, Magn. Reson. Mater. Phys. Biol. Med., 30 (2017), 337–346. https://doi.org/10.1007/s10334-017-0607-2 doi: 10.1007/s10334-017-0607-2
![]() |
[13] |
T. Küstner, A. Bustin, O. Jaubert, R. Hajhosseiny, P. G. Masci, R. Neji, et al., Isotropic 3d cartesian single breath-hold cine mri with multi-bin patch-based low-rank reconstruction, Magn. Reson. Med., 84 (2020), 2018–2033. https://doi.org/10.1002/mrm.28267 doi: 10.1002/mrm.28267
![]() |
[14] |
H. Y. Carr, Steady-state free precession in nuclear magnetic resonance, Phys. Rev., 112 (1958), 1693. https://doi.org/10.1103/PhysRev.112.1693 doi: 10.1103/PhysRev.112.1693
![]() |
[15] |
J. Hamilton, D. Franson, N. Seiberlich, Recent advances in parallel imaging for mri, Prog. Nucl. Magn. Reson. Spectrosc., 101 (2017), 71–95. https://doi.org/10.1016/j.pnmrs.2017.04.002 doi: 10.1016/j.pnmrs.2017.04.002
![]() |
[16] |
M. H. Moghari, T. Geva, A. J. Powell, Prospective heart tracking for whole-heart magnetic resonance angiography, Magn. Reson. Med, 77 (2017), 759–765. https://doi.org/10.1002/mrm.26117 doi: 10.1002/mrm.26117
![]() |
[17] |
M. Lustig, D. L. Donoho, J. M. Santos, J. M. Pauly, Compressed sensing mri, IEEE Signal Process. Mag., 25 (2008), 72–82. https://doi.org/10.1109/MSP.2007.914728 doi: 10.1109/MSP.2007.914728
![]() |
[18] |
M. R. Bashir, L. Bhatti, D. Marin, R. C. Nelson, Emerging applications for ferumoxytol as a contrast agent in mri, J. Magn. Reson. Imaging, 41 (2015), 884–898. https://doi.org/10.1002/jmri.24691 doi: 10.1002/jmri.24691
![]() |
[19] |
Y. Zhang, H. Lin, Y. Li, H. Ma, A patch based denoising method using deep convolutional neural network for seismic image, IEEE Access, 7 (2019), 156883–156894. https://doi.org/10.1109/ACCESS.2019.2949774 doi: 10.1109/ACCESS.2019.2949774
![]() |
[20] |
J. He, L. Ding, L. Jiang, Z. Li, Q. Hu, Intrinsic dimensionality estimation based on manifold assumption, J. Vis. Commun. Image Represent., 25 (2014), 740–747. https://doi.org/10.1016/j.jvcir.2014.01.006 doi: 10.1016/j.jvcir.2014.01.006
![]() |
[21] | T. Kessler, G. Dorian, J. H. Mack, Application of a rectified linear unit (relu) based artificial neural network to cetane number predictions, in Internal Combustion Engine Division Fall Technical Conference, American Society of Mechanical Engineers, 58318 (2017), V001T02A006. https://doi.org/10.1115/ICEF2017-3614 |
[22] |
T. K. Kim, T test as a parametric statistic, Korean J. Anesthesiol., 68 (2015), 540–546. https://doi.org/10.4097/kjae.2015.68.6.540 doi: 10.4097/kjae.2015.68.6.540
![]() |
[23] | G. B. Barrett, The coefficient of determination: Understanding r squared and r squared, Math. Teacher, 93 (2000), 230–234. |
[24] |
D. Giavarina, Understanding bland altman analysis, Biochem. Med., 25 (2015), 141–151. https://doi.org/10.11613/BM.2015.015 doi: 10.11613/BM.2015.015
![]() |
1. | Manal Alqhtani, Kolade M. Owolabi, Khaled M. Saad, Edson Pindza, Spatiotemporal chaos in spatially extended fractional dynamical systems, 2023, 119, 10075704, 107118, 10.1016/j.cnsns.2023.107118 | |
2. | Nadiyah Hussain Alharthi, Abdon Atangana, Badr S. Alkahtani, Some results on the existence of Fractal-fractional ODE with power law, 2023, 20, 25900374, 100398, 10.1016/j.rinam.2023.100398 | |
3. | Tianxian Zhang, Yongqi Zhao, Xiangliang Xu, Si Wu, Yujuan Gu, Solution and dynamics analysis of fractal-fractional multi-scroll Chen chaotic system based on Adomain decomposition method, 2024, 178, 09600779, 114268, 10.1016/j.chaos.2023.114268 | |
4. | Sonal Jain, Salah Boulaaras, A new approach of modeling crossover behaviors in heat conduction is a piecewise heat equation with variable coefficients and constant fractional order, 2024, 0, 1937-1632, 0, 10.3934/dcdss.2024155 | |
5. | Yanli Ma, Xuewu Zuo, Kolade M. Owolabi, Stability Analysis of SIRS Model considering Pulse Vaccination and Elimination Disturbance, 2024, 2024, 2314-4785, 1, 10.1155/2024/6617911 | |
6. | V. S. Verma, Harshita Kaushik, Ram Singh, Sonal Jain, Ali Akgül, Stability and Hopf bifurcation of TB-COVID-19 coinfection model with impact of time delay, 2024, 2254-3902, 10.1007/s40324-024-00366-4 | |
7. | Ahlem Yahiaoui, Med-Salem Rezaoui, Omar Djidel, Rafik Guefaifia, Salah Boulaaras, Kirchhoff-type problems with the non-local fractional d(z,.)-Laplacian operator, 2025, 2025, 1029-242X, 10.1186/s13660-025-03291-w |