Fractal derivative

Last updated January 23, 2024

In applied mathematics and mathematical analysis, the fractal derivative or Hausdorff derivative is a non-Newtonian generalization of the derivative dealing with the measurement of fractals, defined in fractal geometry. Fractal derivatives were created for the study of anomalous diffusion, by which traditional approaches fail to factor in the fractal nature of the media. A fractal measure t is scaled according to t^α. Such a derivative is local, in contrast to the similarly applied fractional derivative. Fractal calculus is formulated as a generalization of standard calculus.^[1]

Physical background

Porous media, aquifers, turbulence, and other media usually exhibit fractal properties. Classical diffusion or dispersion laws based on random walks in free space (essentially the same result variously known as Fick's laws of diffusion, Darcy's law, and Fourier's law) are not applicable to fractal media. To address this, concepts such as distance and velocity must be redefined for fractal media; in particular, scales for space and time are to be transformed according to (x^β, t^α). Elementary physical concepts such as velocity are redefined as follows for fractal spacetime (x^β, t^α):^[2]

v'={\frac {dx'}{dt'}}={\frac {dx^{\beta }}{dt^{\alpha }}}\,,\quad \alpha ,\beta >0

,

where S^α,β represents the fractal spacetime with scaling indices α and β. The traditional definition of velocity makes no sense in the non-differentiable fractal spacetime.^[2]

Definition

Based on above discussion, the concept of the fractal derivative of a function u(t) with respect to a fractal measure t has been introduced as follows:^[3]

{\frac {\partial f(t)}{\partial t^{\alpha }}}=\lim _{t_{1}\rightarrow t}{\frac {f(t_{1})-f(t)}{t_{1}^{\alpha }-t^{\alpha }}}\,,\quad \alpha >0

,

A more general definition is given by^[3]

{\frac {\partial ^{\beta }f(t)}{\partial t^{\alpha }}}=\lim _{t_{1}\rightarrow t}{\frac {f^{\beta }(t_{1})-f^{\beta }(t)}{t_{1}^{\alpha }-t^{\alpha }}}\,,\quad \alpha >0,\beta >0

.

For a function y(t) on $F^{\alpha }$ -perfect fractal set F the fractal derivative or $F^{\alpha }$ -derivative of y(t) at t is defined by

D_{F}^{\alpha }y(t)=\left\{{\begin{array}{ll}{\underset {x\rightarrow t}{F_{-}\lim }}~{\frac {y(x)-y(t)}{S_{F}^{\alpha }(x)-S_{F}^{\alpha }(t)}},&{\text{if}}~t\in F;\\0,&{\text{otherwise}}.\end{array}}\right.

.

Motivation

The derivatives of a function f can be defined in terms of the coefficients a_k in the Taylor series expansion:

$f(x)=\sum _{k=1}^{\infty }a_{k}\cdot (x-x_{0})^{k}=\sum _{k=1}^{\infty }{1 \over k!}{d^{k}f \over dx^{k}}(x_{0})\cdot (x-x_{0})^{k}=f(x_{0})+f'(x_{0})\cdot (x-x_{0})+o(x-x_{0})$

From this approach one can directly obtain:

$f'(x_{0})={f(x)-f(x_{0})-o(x-x_{0}) \over x-x_{0}}=\lim _{x\to x_{0}}{f(x)-f(x_{0}) \over x-x_{0}}$

This can be generalized approximating f with functions (x^α-(x₀)^α)^k:

$f(x)=\sum _{k=1}^{\infty }b_{k}\cdot (x^{\alpha }-x_{0}^{\alpha })^{k}=f(x_{0})+b_{1}\cdot (x^{\alpha }-x_{0}^{\alpha })+o(x^{\alpha }-x_{0}^{\alpha })$

Note that the lowest order coefficient still has to be b₀=f(x₀), since it's still the constant approximation of the function f at x₀.

Again one can directly obtain:

$b_{1}=\lim _{x\to x_{0}}{f(x)-f(x_{0}) \over x^{\alpha }-x_{0}^{\alpha }}{\overset {\underset {\mathrm {def} }{}}{=}}{df \over dx^{\alpha }}(x_{0})$

The Fractal Maclaurin series of f(t) with fractal support F is as follows:

$f(t)=\sum _{m=0}^{\infty }{\frac {(D_{F}^{\alpha })^{m}f(t)|_{t=0}}{m!}}(S_{F}^{\alpha }(t))^{m}$

Properties

Expansion coefficients

Just like in the Taylor series expansion, the coefficients b_k can be expressed in terms of the fractal derivatives of order k of f:

$b_{k}={1 \over k!}{\biggl (}{d \over dx^{\alpha }}{\biggr )}^{k}f(x=x_{0})$

Proof idea: Assuming

${\textstyle ({d \over dx^{\alpha }})^{k}f(x=x_{0})}$

exists, b_k can be written as

${\textstyle b_{k}=a_{k}\cdot ({d \over dx^{\alpha }})^{k}f(x=x_{0})}$

One can now use

${\textstyle f(x)=(x^{\alpha }-x_{0}^{\alpha })^{n}\Rightarrow ({d \over dx^{\alpha }})^{k}f(x=x_{0})=n!\delta _{n}^{k}}$

and since

${\textstyle b_{n}{\overset {\underset {\mathrm {!} }{}}{=}}1\Rightarrow a_{n}={1 \over n!}}$

Chain rule

If for a given function f both the derivative Df and the fractal derivative D_αf exist, one can find an analog to the chain rule:

${df \over dx^{\alpha }}={df \over dx}{dx \over dx^{\alpha }}={1 \over \alpha }x^{1-\alpha }{df \over dx}$

The last step is motivated by the implicit function theorem which, under appropriate conditions, gives us

${\frac {dx}{dx^{\alpha }}}=({\frac {dx^{\alpha }}{dx}})^{-1}$

Similarly for the more general definition:

${d^{\beta }f \over d^{\alpha }x}={d(f^{\beta }) \over d^{\alpha }x}={1 \over \alpha }x^{1-\alpha }\beta f^{\beta -1}(x)f'(x)$

Application in anomalous diffusion

As an alternative modeling approach to the classical Fick's second law, the fractal derivative is used to derive a linear anomalous transport-diffusion equation underlying anomalous diffusion process,^[3]

{\frac {du(x,t)}{dt^{\alpha }}}=D{\frac {\partial }{\partial x^{\beta }}}\left({\frac {\partial u(x,t)}{\partial x^{\beta }}}\right),\quad (1)

u(x,0)=\delta (x)

where 0 < α < 2, 0 < β < 1, $x\in \mathbb {R}$ , and δ(x) is the Dirac delta function.

To obtain the fundamental solution, we apply the transformation of variables

t'=t^{\alpha }\,,\quad x'=x^{\beta }.

then the equation (1) becomes the normal diffusion form equation, the solution of (1) has the stretched Gaussian kernel:^[3]

u(x,t)={\frac {1}{2{\sqrt {\pi t^{\alpha }}}}}e^{-{\frac {x^{2\beta }}{4t^{\alpha }}}}

The mean squared displacement of above fractal derivative diffusion equation has the asymptote:^[3]

\left\langle x^{2}(t)\right\rangle \propto t^{(3\alpha -\alpha \beta )/2\beta }.

Fractal-fractional calculus

The fractal derivative is connected to the classical derivative if the first derivative exists. In this case,

{\frac {\partial f(t)}{\partial t^{\alpha }}}=\lim _{t_{1}\rightarrow t}{\frac {f(t_{1})-f(t)}{t_{1}^{\alpha }-t^{\alpha }}}\ ={\frac {df(t)}{dt}}{\frac {1}{\alpha t^{\alpha -1}}},\quad \alpha >0

.

However, due to the differentiability property of an integral, fractional derivatives are differentiable, thus the following new concept was introduced by Prof Abdon Atangana from South Africa.

The following differential operators were introduced and applied very recently.^[4] Supposing that y(t) be continuous and fractal differentiable on (a, b) with order β, several definitions of a fractal–fractional derivative of y(t) hold with order α in the Riemann–Liouville sense:^[4]

Having power law type kernel:

$^{FFP}D_{0,t}^{\alpha ,\beta }{\Big (}y(t){\Big )}={\dfrac {1}{\Gamma (m-\alpha )}}{\dfrac {d}{dt^{\beta }}}\int _{0}^{t}(t-s)^{m-\alpha -1}y(s)ds$

Having exponentially decaying type kernel:

$^{FFE}D_{0,t}^{\alpha ,\beta }{\Big (}y(t){\Big )}={\dfrac {M(\alpha )}{1-\alpha }}{\dfrac {d}{dt^{\beta }}}\int _{0}^{t}\exp {\Big (}-{\dfrac {\alpha }{1-\alpha }}(t-s){\Big )}y(s)ds$ ,

Having generalized Mittag-Leffler type kernel:

${}_{a}^{FFM}D_{t}^{\alpha }f(t)={\frac {AB(\alpha )}{1-\alpha }}{\frac {d}{dt^{\beta }}}\int _{a}^{t}f(\tau )E_{\alpha }\left(-\alpha {\frac {\left(t-\tau \right)^{\alpha }}{1-\alpha }}\right)\,d\tau \,.$

The above differential operators each have an associated fractal-fractional integral operator, as follows:^[4]

Power law type kernel:

$^{FFP}J_{0,t}^{\alpha ,\beta }{\Big (}y(t){\Big )}={\dfrac {\beta }{\Gamma (\alpha )}}\int _{0}^{t}(t-s)^{\alpha -1}s^{\beta -1}y(s)ds$

Exponentially decaying type kernel:

$^{FFE}J_{0,t}^{\alpha ,\beta }{\Big (}y(t){\Big )}={\dfrac {\alpha \beta }{M(\alpha )}}\int _{0}^{t}s^{\beta -1}y(s)ds+{\dfrac {\beta (1-\alpha )t^{\beta -1}y(t)}{M(\alpha )}}$ .

Generalized Mittag-Leffler type kernel:

$^{FFM}J_{0,t}^{\alpha ,\beta }{\Big (}y(t){\Big )}={\dfrac {\alpha \beta }{AB(\alpha )}}\int _{0}^{t}s^{\beta -1}y(s)(t-s)^{\alpha -1}ds+{\dfrac {\beta (1-\alpha )t^{\beta -1}y(t)}{AB(\alpha )}}$ .

FFM refers to fractal-fractional with the generalized Mittag-Leffler kernel.

Fractal non-local calculus

Fractal analogue of the right-sided Riemann-Liouville fractional integral of order $\beta \in \mathbb {R}$ of f is defined by:^[1]

${x}{\mathcal {I}}_{b}^{\beta }f(x)={\frac {1}{\Gamma (\beta )}}\int _{x}^{b}{\frac {f(t)}{(S_{F}^{\alpha }(t)-S_{F}^{\alpha }(x))^{1-\beta }}}d_{F}^{\alpha }t$ .

Fractal analogue of the left-sided Riemann-Liouville fractional integral of order $\beta \in \mathbb {R}$ of f is defined by:

${a}{\mathcal {I}}_{x}^{\beta }f(x)={\frac {1}{\Gamma (\beta )}}\int _{a}^{x}{\frac {f(t)}{(S_{F}^{\alpha }(x)-S_{F}^{\alpha }(t))^{1-\beta }}}d_{F}^{\alpha }t.$

Fractal analogue of the right-sided Riemann-Liouville fractional derivative of order $\beta \in \mathbb {R}$ of f is defined by:

${x}{\mathcal {D}}_{b}^{\beta }f(x)={\frac {1}{\Gamma (n-\beta )}}(-D_{F}^{\alpha })^{n}\int _{x}^{b}{\frac {f(t)}{(S_{F}^{\alpha }(t)-S_{F}^{\alpha }(x))^{-n+\beta +1}}}d_{F}^{\alpha }t$

Fractal analogue of the left-sided Riemann-Liouville fractional derivative of order $\beta \in \mathbb {R}$ of f is defined by:

${a}{\mathcal {D}}_{x}^{\beta }f(x)={\frac {1}{\Gamma (n-\beta )}}(D_{F}^{\alpha })^{n}\int _{a}^{x}{\frac {f(t)}{(S_{F}^{\alpha }(x)-S_{F}^{\alpha }(t))^{-n+\beta +1}}}d_{F}^{\alpha }t$

Fractal analogue of the right-sided Caputo fractional derivative of order $\beta \in \mathbb {R}$ of f is defined by:

${x}^{C}{\mathcal {D}}_{b}^{\beta }f(x)={\frac {1}{\Gamma (n-\beta )}}\int _{x}^{b}(S_{F}^{\alpha }(t)-S_{F}^{\alpha }(x))^{n-\beta -1}(-D_{F}^{\alpha })^{n}f(t)d_{F}^{\alpha }t$

Fractal analogue of the left-sided Caputo fractional derivative of order $\beta \in \mathbb {R}$ of f is defined by:

${a}^{C}{\mathcal {D}}_{x}^{\beta }f(x)={\frac {1}{\Gamma (n-\beta )}}\int _{a}^{x}(S_{F}^{\alpha }(x)-S_{F}^{\alpha }(t))^{n-\beta -1}(D_{F}^{\alpha })^{n}f(t)d_{F}^{\alpha }t$

Related Research Articles

Bessel functions, first defined by the mathematician Daniel Bernoulli and then generalized by Friedrich Bessel, are canonical solutions $y (x)$ of Bessel's differential equation

In probability theory and statistics, the exponential distribution or negative exponential distribution is the probability distribution of the distance between events in a Poisson point process, i.e., a process in which events occur continuously and independently at a constant average rate; the distance parameter could be any meaningful mono-dimensional measure of the process, such as time between production errors, or length along a roll of fabric in the weaving manufacturing process. It is a particular case of the gamma distribution. It is the continuous analogue of the geometric distribution, and it has the key property of being memoryless. In addition to being used for the analysis of Poisson point processes it is found in various other contexts.

<span class="mw-page-title-main">Taylor's theorem</span> Approximation of a function by a truncated power series

In calculus, Taylor's theorem gives an approximation of a $-times differentiable function around a given point by a polynomial of degree, called the -th-order Taylor polynomial . For a smooth function, the Taylor polynomial is the truncation at the order of the Taylor series of the function. The first-order Taylor polynomial is the linear approximation of the function, and the second-order Taylor polynomial is often referred to as the quadratic approximation . There are several versions of Taylor's theorem, some giving explicit estimates of the approximation error of the function by its Taylor polynomial.$

In probability theory and statistics, the beta distribution is a family of continuous probability distributions defined on the interval [0, 1] or in terms of two positive parameters, denoted by alpha (α) and beta (β), that appear as exponents of the variable and its complement to 1, respectively, and control the shape of the distribution.

Flight dynamics is the science of air vehicle orientation and control in three dimensions. The three critical flight dynamics parameters are the angles of rotation in three dimensions about the vehicle's center of gravity (cg), known as pitch, roll and yaw. These are collectively known as aircraft attitude, often principally relative to the atmospheric frame in normal flight, but also relative to terrain during takeoff or landing, or when operating at low elevation. The concept of attitude is not specific to fixed-wing aircraft, but also extends to rotary aircraft such as helicopters, and dirigibles, where the flight dynamics involved in establishing and controlling attitude are entirely different.

In mathematics, differential forms provide a unified approach to define integrands over curves, surfaces, solids, and higher-dimensional manifolds. The modern notion of differential forms was pioneered by Élie Cartan. It has many applications, especially in geometry, topology and physics.

Fractional calculus is a branch of mathematical analysis that studies the several different possibilities of defining real number powers or complex number powers of the differentiation operator

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

In mathematics, the Hodge star operator or Hodge star is a linear map defined on the exterior algebra of a finite-dimensional oriented vector space endowed with a nondegenerate symmetric bilinear form. Applying the operator to an element of the algebra produces the Hodge dual of the element. This map was introduced by W. V. D. Hodge.

In mathematics, the classical orthogonal polynomials are the most widely used orthogonal polynomials: the Hermite polynomials, Laguerre polynomials, Jacobi polynomials.

In mathematics, the Riemann–Liouville integral associates with a real function $another function Iαf of the same kind for each value of the parameter α > 0. The integral is a manner of generalization of the repeated antiderivative of f in the sense that for positive integer values of α, Iαf is an iterated antiderivative of f of order α. The Riemann–Liouville integral is named for Bernhard Riemann and Joseph Liouville, the latter of whom was the first to consider the possibility of fractional calculus in 1832. The operator agrees with the Euler transform, after Leonhard Euler, when applied to analytic functions. It was generalized to arbitrary dimensions by Marcel Riesz, who introduced the Riesz potential.$

Multi-index notation is a mathematical notation that simplifies formulas used in multivariable calculus, partial differential equations and the theory of distributions, by generalising the concept of an integer index to an ordered tuple of indices.

In quantum field theory, the Lehmann–Symanzik–Zimmermann (LSZ) reduction formula is a method to calculate S-matrix elements from the time-ordered correlation functions of a quantum field theory. It is a step of the path that starts from the Lagrangian of some quantum field theory and leads to prediction of measurable quantities. It is named after the three German physicists Harry Lehmann, Kurt Symanzik and Wolfhart Zimmermann.

In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.

In calculus, the Leibniz integral rule for differentiation under the integral sign states that for an integral of the form

The Engel expansion of a positive real number x is the unique non-decreasing sequence of positive integers $such that$

An $- superprocess,, within mathematics probability theory is a stochastic process on that is usually constructed as a special limit of near-critical branching diffusions.$

In the following we solve the second-order differential equation called the hypergeometric differential equation using Frobenius method, named after Ferdinand Georg Frobenius. This is a method that uses the series solution for a differential equation, where we assume the solution takes the form of a series. This is usually the method we use for complicated ordinary differential equations.

In the fields of dynamical systems and control theory, a fractional-order system is a dynamical system that can be modeled by a fractional differential equation containing derivatives of non-integer order. Such systems are said to have fractional dynamics. Derivatives and integrals of fractional orders are used to describe objects that can be characterized by power-law nonlocality, power-law long-range dependence or fractal properties. Fractional-order systems are useful in studying the anomalous behavior of dynamical systems in physics, electrochemistry, biology, viscoelasticity and chaotic systems.

In mathematics, Katugampola fractional operators are integral operators that generalize the Riemann–Liouville and the Hadamard fractional operators into a unique form. The Katugampola fractional integral generalizes both the Riemann–Liouville fractional integral and the Hadamard fractional integral into a single form and It is also closely related to the Erdelyi–Kober operator that generalizes the Riemann–Liouville fractional integral. Katugampola fractional derivative has been defined using the Katugampola fractional integral and as with any other fractional differential operator, it also extends the possibility of taking real number powers or complex number powers of the integral and differential operators.

References

1 2 Khalili Golmankhaneh, Alireza (2022). Fractal Calculus and its Applications. Singapore: World Scientific Pub Co Inc. p. 328. doi:10.1142/12988. ISBN 978-981-126-110-7. S2CID 248575991.
1 2 Chen, Wen (May 2006). "Time-space fabric underlying anomalous diffusion". Chaos, Solitons & Fractals. 28 (4): 923–929. arXiv: math-ph/0505023 . doi:10.1016/j.chaos.2005.08.199.
1 2 3 4 5 Chen, Wen; Sun, Hongguang; Zhang, Xiaodi; Korošak, Dean (March 2010). "Anomalous diffusion modeling by fractal and fractional derivatives". Computers & Mathematics with Applications. 59 (5): 1754–1758. doi: 10.1016/j.camwa.2009.08.020 .
1 2 3 Atangana, Abdon; Sania, Qureshi (2019). "Modeling attractors of chaotic dynamical systems with fractal–fractional operators". Chaos, Solitons & Fractals. 123: 320–337. Bibcode:2019CSF...123..320A. doi:10.1016/j.chaos.2019.04.020. S2CID 145861887.

Bibliography

Chen, Wen (2006). "Time–space fabric underlying anomalous diffusion". Chaos, Solitons and Fractals. 28 (4): 923–929. arXiv: math-ph/0505023 . Bibcode:2006CSF....28..923C. doi:10.1016/j.chaos.2005.08.199. S2CID 18369880.
Kanno, R. (1998). "Representation of random walk in fractal space-time". Physica A. 248 (1–2): 165–175. Bibcode:1998PhyA..248..165K. doi:10.1016/S0378-4371(97)00422-6.
Chen, W.; Sun, H.G.; Zhang, X.; Korosak, D. (2010). "Anomalous diffusion modeling by fractal and fractional derivatives". Computers and Mathematics with Applications. 59 (5): 1754–8. doi: 10.1016/j.camwa.2009.08.020 .
Sun, H.G.; Meerschaert, M.M.; Zhang, Y.; Zhu, J.; Chen, W. (2013). "A fractal Richards' equation to capture the non-Boltzmann scaling of water transport in unsaturated media". Advances in Water Resources. 52 (52): 292–5. Bibcode:2013AdWR...52..292S. doi:10.1016/j.advwatres.2012.11.005. PMC 3686513 . PMID 23794783.
Cushman, J.H.; O'Malley, D.; Park, M. (2009). "Anomalous diffusion as modeled by a nonstationary extension of Brownian motion". Phys. Rev. E. 79 (3): 032101. Bibcode:2009PhRvE..79c2101C. doi:10.1103/PhysRevE.79.032101. PMID 19391995.
Mainardi, F.; Mura, A.; Pagnini, G. (2010). "The M-Wright Function in Time-Fractional Diffusion Processes: A Tutorial Survey". International Journal of Differential Equations. 2010: 104505. arXiv: 1004.2950 . Bibcode:2010arXiv1004.2950M. doi: 10.1155/2010/104505 . S2CID 37271918.

External links

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[Ali-1] 1 2 Khalili Golmankhaneh, Alireza (2022). Fractal Calculus and its Applications. Singapore: World Scientific Pub Co Inc. p. 328. doi:10.1142/12988. ISBN 978-981-126-110-7. S2CID 248575991.

[:1-2] 1 2 Chen, Wen (May 2006). "Time-space fabric underlying anomalous diffusion". Chaos, Solitons & Fractals. 28 (4): 923–929. arXiv: math-ph/0505023 . doi:10.1016/j.chaos.2005.08.199.

[:0-3] 1 2 3 4 5 Chen, Wen; Sun, Hongguang; Zhang, Xiaodi; Korošak, Dean (March 2010). "Anomalous diffusion modeling by fractal and fractional derivatives". Computers & Mathematics with Applications. 59 (5): 1754–1758. doi: 10.1016/j.camwa.2009.08.020 .

[CDC2019Out-4] 1 2 3 Atangana, Abdon; Sania, Qureshi (2019). "Modeling attractors of chaotic dynamical systems with fractal–fractional operators". Chaos, Solitons & Fractals. 123: 320–337. Bibcode:2019CSF...123..320A. doi:10.1016/j.chaos.2019.04.020. S2CID 145861887.

[1]

[2]

[3]

[4]