Caputo fractional derivative

Last updated September 27, 2024

In mathematics, the Caputo fractional derivative, also called Caputo-type fractional derivative, is a generalization of derivatives for non-integer orders named after Michele Caputo. Caputo first defined this form of fractional derivative in 1967.^[1]

Motivation

The Caputo fractional derivative is motivated from the Riemann–Liouville fractional integral. Let ${\textstyle f}$ be continuous on $\left(0,\,\infty \right)$ , then the Riemann–Liouville fractional integral ${\textstyle {^{\text{RL}}\operatorname {I} }}$ states that

${_{0}^{\text{RL}}\operatorname {I} _{x}^{\alpha }}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(-\alpha \right)}}\cdot \int \limits _{0}^{x}{\frac {f\left(t\right)}{\left(x-t\right)^{1-\alpha }}}\,\operatorname {d} t$

where ${\textstyle \Gamma \left(\cdot \right)}$ is the Gamma function.

Let's define ${\textstyle \operatorname {D} _{x}^{\alpha }:={\frac {\operatorname {d} ^{\alpha }}{\operatorname {d} x^{\alpha }}}}$ , say that ${\textstyle \operatorname {D} _{x}^{\alpha }\operatorname {D} _{x}^{\beta }=\operatorname {D} _{x}^{\alpha +\beta }}$ and that ${\textstyle \operatorname {D} _{x}^{\alpha }={^{\text{RL}}\operatorname {I} _{x}^{-\alpha }}}$ applies. If ${\textstyle \alpha =m+z\in \mathbb {R} \wedge m\in \mathbb {N} _{0}\wedge 0<z<1}$ then we could say ${\textstyle \operatorname {D} _{x}^{\alpha }=\operatorname {D} _{x}^{m+z}=\operatorname {D} _{x}^{z+m}=\operatorname {D} _{x}^{z-1+1+m}=\operatorname {D} _{x}^{z-1}\operatorname {D} _{x}^{1+m}={^{\text{RL}}\operatorname {I} }_{x}^{1-z}\operatorname {D} _{x}^{1+m}}$ . So if $f$ is also $C^{m}\left(0,\,\infty \right)$ , then

${\operatorname {D} _{x}^{m+z}}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(1-z\right)}}\cdot \int \limits _{0}^{x}{\frac {f^{\left(1+m\right)}\left(t\right)}{\left(x-t\right)^{z}}}\,\operatorname {d} t.$

This is known as the Caputo-type fractional derivative, often written as ${\textstyle {^{\text{C}}\operatorname {D} }_{x}^{\alpha }}$ .

Definition

The first definition of the Caputo-type fractional derivative was given by Caputo as:

${^{\text{C}}\operatorname {D} _{x}^{m+z}}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(1-z\right)}}\cdot \int \limits _{0}^{x}{\frac {f^{\left(m+1\right)}\left(t\right)}{\left(x-t\right)^{z}}}\,\operatorname {d} t$

where $C^{m}\left(0,\,\infty \right)$ and ${\textstyle m\in \mathbb {N} _{0}\wedge 0<z<1}$ .^[2]

A popular equivalent definition is:

${^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(\left\lceil \alpha \right\rceil -\alpha \right)}}\cdot \int \limits _{0}^{x}{\frac {f^{\left(\left\lceil \alpha \right\rceil \right)}\left(t\right)}{\left(x-t\right)^{\alpha +1-\left\lceil \alpha \right\rceil }}}\,\operatorname {d} t$

where ${\textstyle \alpha \in \mathbb {R} _{>0}\setminus \mathbb {N} }$ and ${\textstyle \left\lceil \cdot \right\rceil }$ is the ceiling function. This can be derived by substituting ${\textstyle \alpha =m+z}$ so that ${\textstyle \left\lceil \alpha \right\rceil =m+1}$ would apply and ${\textstyle \left\lceil \alpha \right\rceil +z=\alpha +1}$ follows.^[3]

Another popular equivalent definition is given by:

${^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(n-\alpha \right)}}\cdot \int \limits _{0}^{x}{\frac {f^{\left(n\right)}\left(t\right)}{\left(x-t\right)^{\alpha +1-n}}}\,\operatorname {d} t$

where ${\textstyle n-1<\alpha <n\in \mathbb {N} .}$ .

The problem with these definitions is that they only allow arguments in ${\textstyle \left(0,\,\infty \right)}$ . This can be fixed by replacing the lower integral limit with ${\textstyle a}$ : ${\textstyle {_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]={\frac {1}{\Gamma \left(\left\lceil \alpha \right\rceil -\alpha \right)}}\cdot \int \limits _{a}^{x}{\frac {f^{\left(\left\lceil \alpha \right\rceil \right)}\left(t\right)}{\left(x-t\right)^{\alpha +1-\left\lceil \alpha \right\rceil }}}\,\operatorname {d} t}$ . The new domain is ${\textstyle \left(a,\,\infty \right)}$ .^[4]

Properties and theorems

Basic properties and theorems

A few basic properties are:^[5]

A table of basic properties and theorems
Properties	$f\left(x\right)$	${_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]$	Condition
Definition	$f\left(x\right)$	$f^{\left(\alpha \right)}\left(x\right)-f^{\left(\alpha \right)}\left(a\right)$
Linearity	$b\cdot g\left(x\right)+c\cdot h\left(x\right)$	$b\cdot {_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[g\left(x\right)\right]+c\cdot {_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[h\left(x\right)\right]$
Index law	$\operatorname {D} _{x}^{\beta }$	${_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha +\beta }}$	$\beta \in \mathbb {Z}$
Semigroup property	${_{a}^{\text{C}}\operatorname {D} _{x}^{\beta }}$	${_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha +\beta }}$	$\left\lceil \alpha \right\rceil =\left\lceil \beta \right\rceil$

Non-commutation

The index law does not always fulfill the property of commutation:

$\operatorname {_{a}^{\text{C}}D} _{x}^{\alpha }\operatorname {_{a}^{\text{C}}D} _{x}^{\beta }=\operatorname {_{a}^{\text{C}}D} _{x}^{\alpha +\beta }\neq \operatorname {_{a}^{\text{C}}D} _{x}^{\beta }\operatorname {_{a}^{\text{C}}D} _{x}^{\alpha }$

where $\alpha \in \mathbb {R} _{>0}\setminus \mathbb {N} \wedge \beta \in \mathbb {N}$ .

Fractional Leibniz rule

The Leibniz rule for the Caputo fractional derivative is given by:

$\operatorname {_{a}^{\text{C}}D} _{x}^{\alpha }\left[g\left(x\right)\cdot h\left(x\right)\right]=\sum \limits _{k=0}^{\infty }\left[{\binom {a}{k}}\cdot g^{\left(k\right)}\left(x\right)\cdot \operatorname {_{a}^{\text{RL}}D} _{x}^{\alpha -k}\left[h\left(x\right)\right]\right]-{\frac {\left(x-a\right)^{-\alpha }}{\Gamma \left(1-\alpha \right)}}\cdot g\left(a\right)\cdot h\left(a\right)$

where ${\textstyle {\binom {a}{b}}={\frac {\Gamma \left(a+1\right)}{\Gamma \left(b+1\right)\cdot \Gamma \left(a-b+1\right)}}}$ is the binomial coefficient.^[6]^[7]

Relation to other fractional differential operators

Caputo-type fractional derivative is closely related to the Riemann–Liouville fractional integral via its definition:

${_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]={_{a}^{\text{RL}}\operatorname {I} _{x}^{\left\lceil \alpha \right\rceil -\alpha }}\left[\operatorname {D} _{x}^{\left\lceil \alpha \right\rceil }\left[f\left(x\right)\right]\right]$

Furthermore, the following relation applies:

${_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]={_{a}^{\text{RL}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]-\sum \limits _{k=0}^{\left\lceil \alpha \right\rceil }\left[{\frac {x^{k-\alpha }}{\Gamma \left(k-\alpha +1\right)}}\cdot f^{\left(k\right)}\left(0\right)\right]$

where ${_{a}^{\text{RL}}\operatorname {D} _{x}^{\alpha }}$ is the Riemann–Liouville fractional derivative.

Laplace transform

The Laplace transform of the Caputo-type fractional derivative is given by:

${\mathcal {L}}_{x}\left\{{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[f\left(x\right)\right]\right\}\left(s\right)=s^{\alpha }\cdot F\left(s\right)-\sum \limits _{k=0}^{\left\lceil \alpha \right\rceil }\left[s^{\alpha -k-1}\cdot f^{\left(k\right)}\left(0\right)\right]$

where ${\textstyle {\mathcal {L}}_{x}\left\{f\left(x\right)\right\}\left(s\right)=F\left(s\right)}$ .^[8]

Caputo fractional derivative of some functions

The Caputo fractional derivative of a constant $c$ is given by:

${\begin{aligned}{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[c\right]&={\frac {1}{\Gamma \left(\left\lceil \alpha \right\rceil -\alpha \right)}}\cdot \int \limits _{a}^{x}{\frac {\operatorname {D} _{t}^{\left\lceil \alpha \right\rceil }\left[c\right]}{\left(x-t\right)^{\alpha +1-\left\lceil \alpha \right\rceil }}}\,\operatorname {d} t={\frac {1}{\Gamma \left(\left\lceil \alpha \right\rceil -\alpha \right)}}\cdot \int \limits _{a}^{x}{\frac {0}{\left(x-t\right)^{\alpha +1-\left\lceil \alpha \right\rceil }}}\,\operatorname {d} t\\{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[c\right]&=0\end{aligned}}$

The Caputo fractional derivative of a power function $x^{b}$ is given by:^[9]

${\begin{aligned}{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[x^{b}\right]&={_{a}^{\text{RL}}\operatorname {I} _{x}^{\left\lceil \alpha \right\rceil -\alpha }}\left[\operatorname {D} _{x}^{\left\lceil \alpha \right\rceil }\left[x^{b}\right]\right]={\frac {\Gamma \left(b+1\right)}{\Gamma \left(b-\left\lceil \alpha \right\rceil +1\right)}}\cdot {_{a}^{\text{RL}}\operatorname {I} _{x}^{\left\lceil \alpha \right\rceil -\alpha }}\left[x^{b-\left\lceil \alpha \right\rceil }\right]\\{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[x^{b}\right]&={\begin{cases}{\frac {\Gamma \left(b+1\right)}{\Gamma \left(b-\alpha +1\right)}}\left(x^{b-\alpha }-a^{b-\alpha }\right),\,&{\text{for }}\left\lceil \alpha \right\rceil -1<b\wedge b\in \mathbb {R} \\0,\,&{\text{for }}\left\lceil \alpha \right\rceil -1\geq b\wedge b\in \mathbb {N} \\\end{cases}}\end{aligned}}$

The Caputo fractional derivative of a exponential function $e^{a\cdot x}$ is given by:

${\begin{aligned}{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[e^{b\cdot x}\right]&={_{a}^{\text{RL}}\operatorname {I} _{x}^{\left\lceil \alpha \right\rceil -\alpha }}\left[\operatorname {D} _{x}^{\left\lceil \alpha \right\rceil }\left[e^{b\cdot x}\right]\right]=b^{\left\lceil \alpha \right\rceil }\cdot {_{a}^{\text{RL}}\operatorname {I} _{x}^{\left\lceil \alpha \right\rceil -\alpha }}\left[e^{b\cdot x}\right]\\{_{a}^{\text{C}}\operatorname {D} _{x}^{\alpha }}\left[e^{b\cdot x}\right]&=b^{\alpha }\cdot \left(E_{x}\left(\left\lceil \alpha \right\rceil -\alpha ,\,b\right)-E_{a}\left(\left\lceil \alpha \right\rceil -\alpha ,\,b\right)\right)\\\end{aligned}}$

where ${\textstyle E_{x}\left(\nu ,\,a\right)={\frac {a^{-\nu }\cdot e^{a\cdot x}\cdot \gamma \left(\nu ,\,a\cdot x\right)}{\Gamma \left(\nu \right)}}}$ is the ${\textstyle \operatorname {E} _{t}}$ -function and ${\textstyle \gamma \left(a,\,b\right)}$ is the lower incomplete gamma function.^[10]

Related Research Articles

In mathematical physics and mathematics, the Pauli matrices are a set of three $2 \times 2$ complex matrices that are traceless, Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

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.

In mathematics, the floor function is the function that takes as input a real number $x$ , and gives as output the greatest integer less than or equal to $x$ , denoted $⌊ x ⌋$ or $floor(x)$ . Similarly, the ceiling function maps $x$ to the smallest integer greater than or equal to $x$ , denoted $⌈ x ⌉$ or $ceil(x)$ .

In trigonometry, the law of sines, sine law, sine formula, or sine rule is an equation relating the lengths of the sides of any triangle to the sines of its angles. According to the law, $where a, b, and c are the lengths of the sides of a triangle, and α, β, and γ are the opposite angles, while R is the radius of the triangle's circumcircle. When the last part of the equation is not used, the law is sometimes stated using the reciprocals; The law of sines can be used to compute the remaining sides of a triangle when two angles and a side are known—a technique known as triangulation. It can also be used when two sides and one of the non-enclosed angles are known. In some such cases, the triangle is not uniquely determined by this data and the technique gives two possible values for the enclosed angle.$

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 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.$

<span class="mw-page-title-main">Dirichlet distribution</span> Probability distribution

In probability and statistics, the Dirichlet distribution, often denoted $, is a family of continuous multivariate probability distributions parameterized by a vector of positive reals. It is a multivariate generalization of the beta distribution, hence its alternative name of multivariate beta distribution (MBD). Dirichlet distributions are commonly used as prior distributions in Bayesian statistics, and in fact, the Dirichlet distribution is the conjugate prior of the categorical distribution and multinomial distribution.$

In mathematics, a Lambert series, named for Johann Heinrich Lambert, is a series taking the form

In the field of mathematics, norms are defined for elements within a vector space. Specifically, when the vector space comprises matrices, such norms are referred to as matrix norms. Matrix norms differ from vector norms in that they must also interact with matrix multiplication.

In mathematics, the Mittag-Leffler function $is a special function, a complex function which depends on two complex parameters and . It may be defined by the following series when the real part of is strictly positive:$

An osculating circle is a circle that best approximates the curvature of a curve at a specific point. It is tangent to the curve at that point and has the same curvature as the curve at that point. The osculating circle provides a way to understand the local behavior of a curve and is commonly used in differential geometry and calculus.

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

Expected shortfall (ES) is a risk measure—a concept used in the field of financial risk measurement to evaluate the market risk or credit risk of a portfolio. The "expected shortfall at q% level" is the expected return on the portfolio in the worst $of cases. ES is an alternative to value at risk that is more sensitive to the shape of the tail of the loss distribution.$

A ratio distribution is a probability distribution constructed as the distribution of the ratio of random variables having two other known distributions. Given two random variables X and Y, the distribution of the random variable Z that is formed as the ratio Z = X/Y is a ratio distribution.

<span class="mw-page-title-main">Normal-inverse-gamma distribution</span>

In probability theory and statistics, the normal-inverse-gamma distribution is a four-parameter family of multivariate continuous probability distributions. It is the conjugate prior of a normal distribution with unknown mean and variance.

A product distribution is a probability distribution constructed as the distribution of the product of random variables having two other known distributions. Given two statistically independent random variables X and Y, the distribution of the random variable Z that is formed as the product $is a product distribution .$

In mathematics, the Fox–Wright function (also known as Fox–Wright Psi function, not to be confused with Wright Omega function) is a generalisation of the generalised hypergeometric function _pF_q(z) based on ideas of Charles Fox (1928) and E. Maitland Wright (1935):

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 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

↑ Diethelm, Kai (2019). "General theory of Caputo-type fractional differential equations". Fractional Differential Equations. pp. 1–20. doi:10.1515/9783110571660-001. ISBN 978-3-11-057166-0 . Retrieved 2023-08-10.
↑ Caputo, Michele (1967). "Linear Models of Dissipation whose Q is almost Frequency Independent-II". ResearchGate. 13 (5): 530. Bibcode:1967GeoJ...13..529C. doi: 10.1111/j.1365-246X.1967.tb02303.x .
↑ Lazarević, Mihailo; Rapaić, Milan Rade; Šekara, Tomislav (2014). "Introduction to Fractional Calculus with Brief Historical Background". ResearchGate: 8.
↑ Dimitrov, Yuri; Georgiev, Slavi; Todorov, Venelin (2023). "Approximation of Caputo Fractional Derivative and Numerical Solutions of Fractional Differential Equations". Fractal and Fractional. 7 (10): 750. doi: 10.3390/fractalfract7100750 .
↑ Sikora, Beata (2023). "Remarks on the Caputo fractional derivative" (PDF). Matematyka I Informatyka Na Uczelniach Technicznych (5): 78–79.
↑ Huseynov, Ismail; Ahmadova, Arzu; Mahmudov, Nazim (2020). "Fractional Leibniz integral rules for Riemann-Liouville and Caputo fractional derivatives and their applications". ResearchGate: 1. arXiv: 2012.11360 .
↑ Weisstein, Eric W. (2024). "Binomial Coefficient". mathworld.wolfram.com. Retrieved 2024-05-20.
↑ Sontakke, Bhausaheb Rajba; Shaikh, Amjad (2015). "Properties of Caputo Operator and Its Applications to Linear Fractional Differential Equations" (PDF). Journal of Engineering Research and Applications. 5 (5): 23–24. ISSN 2248-9622.
↑ Weisstein, Eric W. "Fractional Derivative". mathworld.wolfram.com. Retrieved 2024-05-20.
↑ Weisstein, Eric W. (2024). "E_t-Function". mathworld.wolfram.com. Retrieved 2024-05-20.