Part of a series of articles about |
Calculus |
---|
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
and of the integration operator [Note 1]
and developing a calculus for such operators generalizing the classical one.
In this context, the term powers refers to iterative application of a linear operator to a function , that is, repeatedly composing with itself, as in
For example, one may ask for a meaningful interpretation of
as an analogue of the functional square root for the differentiation operator, that is, an expression for some linear operator that, when applied twice to any function, will have the same effect as differentiation. More generally, one can look at the question of defining a linear operator
for every real number in such a way that, when takes an integer value , it coincides with the usual -fold differentiation if , and with the -th power of when .
One of the motivations behind the introduction and study of these sorts of extensions of the differentiation operator is that the sets of operator powers defined in this way are continuous semigroups with parameter , of which the original discrete semigroup of for integer is a denumerable subgroup: since continuous semigroups have a well developed mathematical theory, they can be applied to other branches of mathematics.
Fractional differential equations, also known as extraordinary differential equations, [1] are a generalization of differential equations through the application of fractional calculus.
In applied mathematics and mathematical analysis, a fractional derivative is a derivative of any arbitrary order, real or complex. Its first appearance is in a letter written to Guillaume de l'Hôpital by Gottfried Wilhelm Leibniz in 1695. [2] Around the same time, Leibniz wrote to one of the Bernoulli brothers describing the similarity between the binomial theorem and the Leibniz rule for the fractional derivative of a product of two functions.[ citation needed ]
Fractional calculus was introduced in one of Niels Henrik Abel's early papers [3] where all the elements can be found: the idea of fractional-order integration and differentiation, the mutually inverse relationship between them, the understanding that fractional-order differentiation and integration can be considered as the same generalized operation, and the unified notation for differentiation and integration of arbitrary real order. [4] Independently, the foundations of the subject were laid by Liouville in a paper from 1832. [5] [6] [7] Oliver Heaviside introduced the practical use of fractional differential operators in electrical transmission line analysis circa 1890. [8] The theory and applications of fractional calculus expanded greatly over the 19th and 20th centuries, and numerous contributors have given different definitions for fractional derivatives and integrals. [9]
Let f(x) be a function defined for x > 0. Form the definite integral from 0 to x. Call this
Repeating this process gives
and this can be extended arbitrarily.
The Cauchy formula for repeated integration, namely leads in a straightforward way to a generalization for real n: using the gamma function to remove the discrete nature of the factorial function gives us a natural candidate for applications of the fractional integral operator as
This is in fact a well-defined operator.
It is straightforward to show that the J operator satisfies
Proof of this identity |
---|
where in the last step we exchanged the order of integration and pulled out the f(s) factor from the t integration. Changing variables to r defined by t = s + (x − s)r, The inner integral is the beta function which satisfies the following property: Substituting back into the equation: Interchanging α and β shows that the order in which the J operator is applied is irrelevant and completes the proof. |
This relationship is called the semigroup property of fractional differintegral operators.
The classical form of fractional calculus is given by the Riemann–Liouville integral, which is essentially what has been described above. The theory of fractional integration for periodic functions (therefore including the "boundary condition" of repeating after a period) is given by the Weyl integral. It is defined on Fourier series, and requires the constant Fourier coefficient to vanish (thus, it applies to functions on the unit circle whose integrals evaluate to zero). The Riemann–Liouville integral exists in two forms, upper and lower. Considering the interval [a,b], the integrals are defined as
Where the former is valid for t > a and the latter is valid for t < b. [10]
It has been suggested [11] that the integral on the positive real axis (i.e. ) would be more appropriately named the Abel–Riemann integral, on the basis of history of discovery and use, and in the same vein the integral over the entire real line be named Liouville–Weyl integral.
By contrast the Grünwald–Letnikov derivative starts with the derivative instead of the integral.
The Hadamard fractional integral was introduced by Jacques Hadamard [12] and is given by the following formula,
The Atangana–Baleanu fractional integral of a continuous function is defined as:
Unfortunately, the comparable process for the derivative operator D is significantly more complex, but it can be shown that D is neither commutative nor additive in general. [13]
Unlike classical Newtonian derivatives, fractional derivatives can be defined in a variety of different ways that often do not all lead to the same result even for smooth functions. Some of these are defined via a fractional integral. Because of the incompatibility of definitions, it is frequently necessary to be explicit about which definition is used.
The corresponding derivative is calculated using Lagrange's rule for differential operators. To find the αth order derivative, the nth order derivative of the integral of order (n − α) is computed, where n is the smallest integer greater than α (that is, n = ⌈α⌉). The Riemann–Liouville fractional derivative and integral has multiple applications such as in case of solutions to the equation in the case of multiple systems such as the tokamak systems, and Variable order fractional parameter. [14] [15] Similar to the definitions for the Riemann–Liouville integral, the derivative has upper and lower variants. [16]
Another option for computing fractional derivatives is the Caputo fractional derivative. It was introduced by Michele Caputo in his 1967 paper. [17] In contrast to the Riemann–Liouville fractional derivative, when solving differential equations using Caputo's definition, it is not necessary to define the fractional order initial conditions. Caputo's definition is illustrated as follows, where again n = ⌈α⌉:
There is the Caputo fractional derivative defined as: which has the advantage that is zero when f(t) is constant and its Laplace Transform is expressed by means of the initial values of the function and its derivative. Moreover, there is the Caputo fractional derivative of distributed order defined as
where ϕ(ν) is a weight function and which is used to represent mathematically the presence of multiple memory formalisms.
In a paper of 2015, M. Caputo and M. Fabrizio presented a definition of fractional derivative with a non singular kernel, for a function of given by:
where . [18]
In 2016, Atangana and Baleanu suggested differential operators based on the generalized Mittag-Leffler function . The aim was to introduce fractional differential operators with non-singular nonlocal kernel. Their fractional differential operators are given below in Riemann–Liouville sense and Caputo sense respectively. For a function of given by [19] [20]
If the function is continuous, the Atangana–Baleanu derivative in Riemann–Liouville sense is given by:
The kernel used in Atangana–Baleanu fractional derivative has some properties of a cumulative distribution function. For example, for all , the function is increasing on the real line, converges to in , and . Therefore, we have that, the function is the cumulative distribution function of a probability measure on the positive real numbers. The distribution is therefore defined, and any of its multiples is called a Mittag-Leffler distribution of order . It is also very well-known that, all these probability distributions are absolutely continuous. In particular, the function Mittag-Leffler has a particular case , which is the exponential function, the Mittag-Leffler distribution of order is therefore an exponential distribution. However, for , the Mittag-Leffler distributions are heavy-tailed. Their Laplace transform is given by:
This directly implies that, for , the expectation is infinite. In addition, these distributions are geometric stable distributions.
The Riesz derivative is defined as
where denotes the Fourier transform. [21] [22]
Classical fractional derivatives include:
New fractional derivatives include:
The Coimbra derivative is used for physical modeling: [31] A number of applications in both mechanics and optics can be found in the works by Coimbra and collaborators, [32] [33] [34] [35] [36] [37] [38] as well as additional applications to physical problems and numerical implementations studied in a number of works by other authors [39] [40] [41] [42]
For
where the lower limit can be taken as either or as long as is identically zero from or to . Note that this operator returns the correct fractional derivatives for all values of and can be applied to either the dependent function itself with a variable order of the form or to the independent variable with a variable order of the form .
The Coimbra derivative can be generalized to any order, [43] leading to the Coimbra Generalized Order Differintegration Operator (GODO) [44]
For
where is an integer larger than the larger value of for all values of . Note that the second (summation) term on the right side of the definition above can be expressed as
so to keep the denominator on the positive branch of the Gamma () function and for ease of numerical calculation.
The -th derivative of a function at a point is a local property only when is an integer; this is not the case for non-integer power derivatives. In other words, a non-integer fractional derivative of at depends on all values of , even those far away from . Therefore, it is expected that the fractional derivative operation involves some sort of boundary conditions, involving information on the function further out. [45]
The fractional derivative of a function of order is nowadays often defined by means of the Fourier or Mellin integral transforms.[ citation needed ]
The Erdélyi–Kober operator is an integral operator introduced by Arthur Erdélyi (1940). [46] and Hermann Kober (1940) [47] and is given by
which generalizes the Riemann–Liouville fractional integral and the Weyl integral.
In the context of functional analysis, functions f(D) more general than powers are studied in the functional calculus of spectral theory. The theory of pseudo-differential operators also allows one to consider powers of D. The operators arising are examples of singular integral operators; and the generalisation of the classical theory to higher dimensions is called the theory of Riesz potentials. So there are a number of contemporary theories available, within which fractional calculus can be discussed. See also Erdélyi–Kober operator, important in special function theory ( Kober 1940 ), ( Erdélyi 1950–1951 ).
As described by Wheatcraft and Meerschaert (2008), [48] a fractional conservation of mass equation is needed to model fluid flow when the control volume is not large enough compared to the scale of heterogeneity and when the flux within the control volume is non-linear. In the referenced paper, the fractional conservation of mass equation for fluid flow is:
When studying the redox behavior of a substrate in solution, a voltage is applied at an electrode surface to force electron transfer between electrode and substrate. The resulting electron transfer is measured as a current. The current depends upon the concentration of substrate at the electrode surface. As substrate is consumed, fresh substrate diffuses to the electrode as described by Fick's laws of diffusion. Taking the Laplace transform of Fick's second law yields an ordinary second-order differential equation (here in dimensionless form):
whose solution C(x,s) contains a one-half power dependence on s. Taking the derivative of C(x,s) and then the inverse Laplace transform yields the following relationship:
which relates the concentration of substrate at the electrode surface to the current. [49] This relationship is applied in electrochemical kinetics to elucidate mechanistic behavior. For example, it has been used to study the rate of dimerization of substrates upon electrochemical reduction. [50]
In 2013–2014 Atangana et al. described some groundwater flow problems using the concept of a derivative with fractional order. [51] [52] In these works, the classical Darcy law is generalized by regarding the water flow as a function of a non-integer order derivative of the piezometric head. This generalized law and the law of conservation of mass are then used to derive a new equation for groundwater flow.
This equation[ clarification needed ] has been shown useful for modeling contaminant flow in heterogenous porous media. [53] [54] [55]
Atangana and Kilicman extended the fractional advection dispersion equation to a variable order equation. In their work, the hydrodynamic dispersion equation was generalized using the concept of a variational order derivative. The modified equation was numerically solved via the Crank–Nicolson method. The stability and convergence in numerical simulations showed that the modified equation is more reliable in predicting the movement of pollution in deformable aquifers than equations with constant fractional and integer derivatives [56]
Anomalous diffusion processes in complex media can be well characterized by using fractional-order diffusion equation models. [57] [58] The time derivative term corresponds to long-time heavy tail decay and the spatial derivative for diffusion nonlocality. The time-space fractional diffusion governing equation can be written as
A simple extension of the fractional derivative is the variable-order fractional derivative, α and β are changed into α(x, t) and β(x, t). Its applications in anomalous diffusion modeling can be found in the reference. [56] [59] [60]
Fractional derivatives are used to model viscoelastic damping in certain types of materials like polymers. [11]
Generalizing PID controllers to use fractional orders can increase their degree of freedom. The new equation relating the control variableu(t) in terms of a measured error valuee(t) can be written as
where α and β are positive fractional orders and Kp, Ki, and Kd, all non-negative, denote the coefficients for the proportional, integral, and derivative terms, respectively (sometimes denoted P, I, and D). [61]
The propagation of acoustical waves in complex media, such as in biological tissue, commonly implies attenuation obeying a frequency power-law. This kind of phenomenon may be described using a causal wave equation which incorporates fractional time derivatives:
See also Holm & Näsholm (2011) [62] and the references therein. Such models are linked to the commonly recognized hypothesis that multiple relaxation phenomena give rise to the attenuation measured in complex media. This link is further described in Näsholm & Holm (2011b) [63] and in the survey paper, [64] as well as the Acoustic attenuation article. See Holm & Nasholm (2013) [65] for a paper which compares fractional wave equations which model power-law attenuation. This book on power-law attenuation also covers the topic in more detail. [66]
Pandey and Holm gave a physical meaning to fractional differential equations by deriving them from physical principles and interpreting the fractional-order in terms of the parameters of the acoustical media, example in fluid-saturated granular unconsolidated marine sediments. [67] Interestingly, Pandey and Holm derived Lomnitz's law in seismology and Nutting's law in non-Newtonian rheology using the framework of fractional calculus. [68] Nutting's law was used to model the wave propagation in marine sediments using fractional derivatives. [67]
The fractional Schrödinger equation, a fundamental equation of fractional quantum mechanics, has the following form: [69] [70]
where the solution of the equation is the wavefunction ψ(r, t) – the quantum mechanical probability amplitude for the particle to have a given position vector r at any given time t, and ħ is the reduced Planck constant. The potential energy function V(r, t) depends on the system.
Further, is the Laplace operator, and Dα is a scale constant with physical dimension [Dα] = J1 − α·mα·s−α = kg1 − α·m2 − α·sα − 2, (at α = 2, for a particle of mass m), and the operator (−ħ2Δ)α/2 is the 3-dimensional fractional quantum Riesz derivative defined by
The index α in the fractional Schrödinger equation is the Lévy index, 1 < α ≤ 2.
As a natural generalization of the fractional Schrödinger equation, the variable-order fractional Schrödinger equation has been exploited to study fractional quantum phenomena: [71]
where is the Laplace operator and the operator (−ħ2Δ)β(t)/2 is the variable-order fractional quantum Riesz derivative.
In mechanics, the virial theorem provides a general equation that relates the average over time of the total kinetic energy of a stable system of discrete particles, bound by a conservative force, with that of the total potential energy of the system. Mathematically, the theorem states that where T is the total kinetic energy of the N particles, Fk represents the force on the kth particle, which is located at position rk, and angle brackets represent the average over time of the enclosed quantity. The word virial for the right-hand side of the equation derives from vis, the Latin word for "force" or "energy", and was given its technical definition by Rudolf Clausius in 1870.
In mathematical analysis, the Dirac delta function, also known as the unit impulse, is a generalized function on the real numbers, whose value is zero everywhere except at zero, and whose integral over the entire real line is equal to one. Thus it can be represented heuristically as
Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems published by mathematician Emmy Noether in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.
In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).
In fractional calculus, an area of mathematical analysis, the differintegral is a combined differentiation/integration operator. Applied to a function ƒ, the q-differintegral of f, here denoted by
In mathematical analysis, initialization of the differintegrals is a topic in fractional calculus, a branch of mathematics dealing with derivatives of non-integer order.
In applied mathematics, discretization is the process of transferring continuous functions, models, variables, and equations into discrete counterparts. This process is usually carried out as a first step toward making them suitable for numerical evaluation and implementation on digital computers. Dichotomization is the special case of discretization in which the number of discrete classes is 2, which can approximate a continuous variable as a binary variable.
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.
In probability theory and related fields, Malliavin calculus is a set of mathematical techniques and ideas that extend the mathematical field of calculus of variations from deterministic functions to stochastic processes. In particular, it allows the computation of derivatives of random variables. Malliavin calculus is also called the stochastic calculus of variations. P. Malliavin first initiated the calculus on infinite dimensional space. Then, the significant contributors such as S. Kusuoka, D. Stroock, J-M. Bismut, Shinzo Watanabe, I. Shigekawa, and so on finally completed the foundations.
In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.
In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.
In calculus, the Leibniz integral rule for differentiation under the integral sign, named after Gottfried Wilhelm Leibniz, states that for an integral of the form where and the integrands are functions dependent on the derivative of this integral is expressible as where the partial derivative indicates that inside the integral, only the variation of with is considered in taking the derivative.
In many-body theory, the term Green's function is sometimes used interchangeably with correlation function, but refers specifically to correlators of field operators or creation and annihilation operators.
In mathematics – specifically, in stochastic analysis – an Itô diffusion is a solution to a specific type of stochastic differential equation. That equation is similar to the Langevin equation used in physics to describe the Brownian motion of a particle subjected to a potential in a viscous fluid. Itô diffusions are named after the Japanese mathematician Kiyosi Itô.
In mathematics, the Struve functionsHα(x), are solutions y(x) of the non-homogeneous Bessel's differential equation:
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.
In theoretical physics, relativistic Lagrangian mechanics is Lagrangian mechanics applied in the context of special relativity and general relativity.
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.
Prabhakar function is a certain special function in mathematics introduced by the Indian mathematician Tilak Raj Prabhakar in a paper published in 1971. The function is a three-parameter generalization of the well known two-parameter Mittag-Leffler function in mathematics. The function was originally introduced to solve certain classes of integral equations. Later the function was found to have applications in the theory of fractional calculus and also in certain areas of physics.
{{cite book}}
: |journal=
ignored (help){{cite book}}
: |work=
ignored (help)