1.前言
分数阶微积分具有历史依赖性,在描述复杂系统方面具有明显优势n分数阶微积分是将整数阶微积分算子的阶数推广到实数或者复数。 目前在理论研究和实际工程中,分数阶微积分有三种主要定义方式,分别是Grünwald-Letnikov、 Riemann-Liouville 和 Caputo 定义。 Grünwald-Letnikov 定义是将整数阶微积分的差分极限直接推广而来,适用于数值求解。 Riemann-Liouville 定义则采用微分-积分形式,避免了复杂的极限运算,更适用于数学统计分析。 Caputo定义则采用积分-微分形式,其 Laplace 变换简洁明了,在实际工程中被广泛应用。
2.Grünwald-Letnikov 定义
Grünwald-Letnikov 分数阶导数是整数阶导数差分定义的推广。考虑n 阶可
导的函数f(t),它的 1 阶导数定义为:
f
(
1
)
(
t
)
=
lim
h
→
0
f
(
t
)
−
f
(
t
−
h
)
h
f^{(1)}(t)=\lim_{h\to0}\frac{f(t)-f(t-h)}h
f(1)(t)=h→0limhf(t)−f(t−h)
因此,n阶导数为:
f
(
n
)
(
t
)
=
lim
h
→
0
∑
k
=
0
n
(
−
1
)
k
(
n
k
)
f
(
t
−
k
h
)
h
n
f^{(n)}(t)=\underset{h\to0}{\operatorname*{\lim}}\frac{\sum_{k=0}^n\left(-1\right)^k\binom nkf(t-kh)}{h^n}
f(n)(t)=h→0limhn∑k=0n(−1)k(kn)f(t−kh)
其中,
(
n
k
)
=
n
(
n
−
1
)
(
n
−
2
)
⋯
(
n
−
k
+
1
)
k
!
\binom nk=\frac{n(n-1)(n-2)\cdots(n-k+1)}{k!}
(kn)=k!n(n−1)(n−2)⋯(n−k+1)
Grünwald-Letnikov 分数阶导数是将正整数推广为正实数
α
\alpha
α。 对有限的区域[a,b]内的函数f(t),函数f(t)的
α
\alpha
α 阶左、右分数阶导数的定义分别为:
α
G
D
t
α
f
(
t
)
=
lim
h
→
0
+
1
h
α
∑
k
=
0
[
t
−
a
h
]
(
−
1
k
)
k
(
α
k
)
f
(
t
−
k
h
)
t
G
D
b
α
f
(
t
)
=
lim
h
→
0
+
1
h
α
∑
k
=
0
[
b
−
t
h
]
(
−
1
k
)
k
(
α
k
)
f
(
t
+
k
h
)
\begin{aligned}&_\alpha^GD_t^\alpha f(t)=\lim_{h\to0^+}\frac1{h^\alpha}\sum_{k=0}^{\left[\frac{t-a}h\right]}\binom{-1}k^k\binom\alpha kf(t-kh)\\&_t^GD_b^\alpha f(t)=\lim_{h\to0^+}\frac1{h^\alpha}\sum_{k=0}^{\left[\frac{b-t}h\right]}\binom{-1}k^k\binom\alpha kf(t+kh)\end{aligned}
αGDtαf(t)=h→0+limhα1k=0∑[ht−a](k−1)k(kα)f(t−kh)tGDbαf(t)=h→0+limhα1k=0∑[hb−t](k−1)k(kα)f(t+kh)
当
α
∈
[
0
,
1
)
\alpha \in [0,1)
α∈[0,1) 时, Grünwald-Letnikov 分数阶导数的 Laplace变换为:
L
[
G
0
D
t
α
f
(
t
)
]
=
s
α
F
(
s
)
L\bigg[\begin{smallmatrix}G\\0\end{smallmatrix}D_t^\alpha f(t)\bigg]=s^\alpha F(s)
L[G0Dtαf(t)]=sαF(s)
当
α
>
1
\alpha>1
α>1时,Grünwald-Letnikov 分数阶导数的 Laplace 变换不存在。
3.Riemann-Liouville 定义
Riemann-Liouville 分数阶导数是分数阶导数的一种最早被提出的形式。相
比于 Grünwald-Letnikov 分数阶导数, Riemann-Liouville 分数阶导数在计算时更为简便。 Riemann-Liouville 分数阶导数理论发展较为完备, 是目前最常用的分数阶微积分定义。
设f(t)在[a,b]逐段连续,则函数f(t)的
α
\alpha
α阶 Riemann-Liouville 型分数阶左、右导数定义分别为:
a
R
L
D
t
α
f
(
t
)
=
1
Γ
(
α
)
d
n
d
t
n
∫
a
t
(
t
−
u
)
n
−
α
−
1
f
(
u
)
d
u
,
n
−
1
≤
α
≤
n
,
t
R
L
D
b
α
f
(
t
)
=
1
Γ
(
α
)
d
n
d
t
n
∫
t
b
(
t
−
u
)
n
−
α
−
1
f
(
u
)
d
u
,
n
−
1
≤
α
≤
n
\begin{gathered} _{a}^{RL}D_{t}^{\alpha}f(t)=\frac1{\Gamma\left(\alpha\right)}\frac{d^{n}}{dt^{n}}\int_{a}^{t}\left(t-u\right)^{n-\alpha-1}f(u)du,n-1\leq\alpha\leq n, \\ _{t}^{RL}D_{b}^{\alpha}f(t)=\frac1{\Gamma\left(\alpha\right)}\frac{d^{n}}{dt^{n}}\int_{t}^{b}\left(t-u\right)^{n-\alpha-1}f(u)du,n-1\leq\alpha\leq n \end{gathered}
aRLDtαf(t)=Γ(α)1dtndn∫at(t−u)n−α−1f(u)du,n−1≤α≤n,tRLDbαf(t)=Γ(α)1dtndn∫tb(t−u)n−α−1f(u)du,n−1≤α≤n
其中,
Γ
(
α
)
\Gamma(\alpha)
Γ(α)是Gamma 函数:
Γ
(
α
)
=
∫
0
∞
e
−
t
t
α
−
1
d
t
,
R
e
(
α
)
>
0
\Gamma(\alpha)=\int_0^\infty e^{-t}t^{\alpha-1}dt,\mathrm{Re}(\alpha)>0
Γ(α)=∫0∞e−ttα−1dt,Re(α)>0
Riemann-Liouville 分数阶导数是先对函数做
n
−
α
n-\alpha
n−α阶分数阶积分,然后进行n 次导数。这样定义使得该算子的定义域扩大,不需要函数的连续性就可以对其求分数阶导数。
设f(t)定义在区间
(
0
,
∞
)
\left(0,\infty\right)
(0,∞)上, Riemann-Liouville 分数阶导数的 Laplace 变换:
L
[
R
L
0
D
t
α
f
(
t
)
]
=
s
α
F
(
s
)
−
∑
j
=
0
n
−
1
s
j
[
R
L
0
D
t
α
−
j
−
1
f
(
t
)
]
t
=
0
,
n
−
1
≤
α
<
n
L\bigg[\begin{array}{c}RL\\0\end{array}D_t^\alpha f(t)\bigg]=s^\alpha F(s)-\sum_{j=0}^{n-1}s^j\bigg[\begin{array}{c}RL\\0\end{array}D_t^{\alpha-j-1}f(t)\bigg]_{t=0},n-1\leq\alpha<n
L[RL0Dtαf(t)]=sαF(s)−j=0∑n−1sj[RL0Dtα−j−1f(t)]t=0,n−1≤α<n
4.Caputo 定义
Caputo 分数阶导数的定义与 Riemann-Liouville 分数阶导数定义类似,只是
交换了积分和微分的顺序。对于阶数
α
>
0
\alpha>0
α>0 , Caputo 分数阶导数是先求n 阶导数,再求
α
−
n
\alpha-n
α−n阶积分。
设函数f(t)在区间[a,b]上n 阶可导, 则函数f(t)的
α
\alpha
α阶Caputo型分数阶导数:
a
c
D
t
α
f
(
t
)
=
1
Γ
(
n
−
α
)
d
n
d
t
n
∫
a
t
(
t
−
u
)
n
−
α
−
1
f
(
n
)
(
u
)
d
u
,
n
−
1
<
α
<
n
_a^cD_t^\alpha f(t)=\frac1{\Gamma(n-\alpha)}\frac{d^n}{dt^n}\int_a^t\left(t-u\right)^{n-\alpha-1}f^{(n)}(u)du,n-1<\alpha<n
acDtαf(t)=Γ(n−α)1dtndn∫at(t−u)n−α−1f(n)(u)du,n−1<α<n
Caputo 分数阶导数对函数f(t)要求较高,需要函数f(t)满足n 阶可导条件,其经常用在初值问题中对时间求导数,在工程和物理中应用较为广泛。
设f(t)定义在区间
(
0
,
∞
)
\left(0,\infty\right)
(0,∞)上, Caputo 分数阶导数的 Laplace 变换:
L
[
C
0
D
t
α
f
(
t
)
]
=
s
α
F
(
s
)
−
∑
j
=
0
n
−
1
s
α
−
j
−
1
f
(
j
)
(
0
)
,
n
−
1
<
α
≤
n
L\bigg[\begin{smallmatrix}C\\0\end{smallmatrix}D_t^\alpha f(t)\bigg]=s^\alpha F(s)-\sum_{j=0}^{n-1}s^{\alpha-j-1}f^{(j)}(0),n-1<\alpha\leq n
L[C0Dtαf(t)]=sαF(s)−j=0∑n−1sα−j−1f(j)(0),n−1<α≤n
5.参考文献
[1] Debnath L. A brief historical introduction to fractional calculus[J]. International Journal of Mathematical Education in Science and Technology, 2004, 35(4): 487-501.
[2] Diethelm K. The analysis of fractional differential equations: an application-oriented exposition using differential operators of Caputo type[M]. Berlin: Springer, 2010.