數學 上,黎曼-劉維爾積分 (Riemann–Liouville integral )將定義域和值域均為實數集的實函數
f
:
R
→
R
{\displaystyle f:\mathbb {R} \rightarrow \mathbb {R} }
與另一個同類型的函數 I α f 相關聯,其中參數 α > 0 。該積分是f 的重複反導數 (原函數)的一種推廣形式,當α 為正整數值時, I α f 就是f 的α 階反導數。黎曼-劉維爾積分以伯恩哈德·黎曼 和約瑟夫·劉維爾 的名字命名,後者於1832年首次考慮分數微積分 的可能性。[ 1] [ 2] [ 3] [ 4] 當應用於解析函數 時,該算符與萊昂哈德·歐拉 (Leonhard Euler) 提出的歐拉變換 一致。 [ 5] 它被Marcel Riesz 推廣到任意維度,並引入了Riesz勢 。該算符也被稱為黎曼-劉維爾微分積分 (Riemann–Liouville differintegral ),微分積分 (differintegral)一詞體現它是一種微分和積分的聯合運算。
黎曼-劉維爾積分源自柯西重複積分公式。對於在區間 [a , x ] 上連續的函數f ,柯西n 重重複積分公式如下
I
n
f
(
x
)
=
f
(
−
n
)
(
x
)
=
1
(
n
−
1
)
!
∫
a
x
(
x
−
t
)
n
−
1
f
(
t
)
d
t
.
{\displaystyle I^{n}f(x)=f^{(-n)}(x)={\frac {1}{(n-1)!}}\int _{a}^{x}\left(x-t\right)^{n-1}f(t)\,\mathrm {d} t.}
現在,通過用α 代替正整數n ,這個公式可以推廣到任何正實數,因此我們得到黎曼-劉維爾分數階積分的定義
I
α
f
(
x
)
=
1
Γ
(
α
)
∫
a
x
(
x
−
t
)
α
−
1
f
(
t
)
d
t
{\displaystyle I^{\alpha }f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}\left(x-t\right)^{\alpha -1}f(t)\,\mathrm {d} t}
黎曼-劉維爾積分定義為
I
α
f
(
x
)
=
1
Γ
(
α
)
∫
a
x
f
(
t
)
(
x
−
t
)
α
−
1
d
t
{\displaystyle I^{\alpha }f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}f(t)(x-t)^{\alpha -1}\,dt}
其中Γ 是伽馬函數 , a 是任意但固定的基點。如果f 是局部可積函數 ,且α 是半平面 Re(α ) > 0 中的複數 ,則積分定義明確。其中對基點a 的依賴往往被隱去,體現了積分常數 的選取自由。顯然, I 1 f 是f 的(一階)反導數,當α 取正整數值時,根據柯西重複積分公式 , I α f 是其α 階反導數。另一種強調基點的記法如下:[ 6]
a
D
x
−
α
f
(
x
)
=
1
Γ
(
α
)
∫
a
x
f
(
t
)
(
x
−
t
)
α
−
1
d
t
.
{\displaystyle {}_{a}D_{x}^{-\alpha }f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}f(t)(x-t)^{\alpha -1}\,dt.}
如果a = −∞ ,且對f 有適當的限制,上述式子也是適用的。
有以下基本關係成立:
d
d
x
I
α
+
1
f
(
x
)
=
I
α
f
(
x
)
,
I
α
(
I
β
f
)
=
I
α
+
β
f
,
{\displaystyle {\frac {d}{dx}}I^{\alpha +1}f(x)=I^{\alpha }f(x),\quad I^{\alpha }(I^{\beta }f)=I^{\alpha +\beta }f,}
後式表明其具有半群 性質。 [ 1] 這些性質不僅使得分數階積分的定義成為可能,而且通過對I α f 取足夠多的導數,還可以定義分數階微分。
對一個確定有界區間(a, b) 。算子I α 將每個可積函數 f 在(a, b)上與函數I α f 在(a, b)上相關聯。根據富比尼定理 ,函數I α f 也是可積的。因此I α 在L 1 (a ,b ) 上定義了一個線性算子 :
I
α
:
L
1
(
a
,
b
)
→
L
1
(
a
,
b
)
.
{\displaystyle I^{\alpha }:L^{1}(a,b)\to L^{1}(a,b).}
富比尼定理還表明,該算子關於L 1 上的Banach 空間 結構是連續的 ,並且以下不等式成立:
‖
I
α
f
‖
1
≤
|
b
−
a
|
ℜ
(
α
)
ℜ
(
α
)
|
Γ
(
α
)
|
‖
f
‖
1
.
{\displaystyle \left\|I^{\alpha }f\right\|_{1}\leq {\frac {|b-a|^{\Re (\alpha )}}{\Re (\alpha )|\Gamma (\alpha )|}}\|f\|_{1}.}
這裡‖ · ‖1 表示L 1 (a ,b ) 上的範數 。
更一般地,根據赫爾德不等式 ,如果f ∈ L p (a , b ) ,則I α f ∈ L p (a , b ) 同樣成立。下列類似的不等式
‖
I
α
f
‖
p
≤
|
b
−
a
|
ℜ
(
α
)
/
p
ℜ
(
α
)
|
Γ
(
α
)
|
‖
f
‖
p
{\displaystyle \left\|I^{\alpha }f\right\|_{p}\leq {\frac {|b-a|^{\Re (\alpha )/p}}{\Re (\alpha )|\Gamma (\alpha )|}}\|f\|_{p}}
成立。其中‖ · ‖p 是區間(a, b)上的L p 範數 。
因此我們有一個有界線性算子I α : L p (a , b ) → L p (a , b ) 。此外,當參數α沿實軸趨於0,在L p 意義上, I α f → f 。即對於所有p ≥ 1
lim
α
→
0
+
‖
I
α
f
−
f
‖
p
=
0
{\displaystyle \lim _{\alpha \to 0^{+}}\|I^{\alpha }f-f\|_{p}=0}
此外,通過估計I 的最大函數,可以證明極限I α f → f 幾乎處處 逐點成立。
算子I α 在整個實數軸上的局部可積函數集上有明確定義
R
{\displaystyle \mathbb {R} }
。它定義了任意Banach 空間 中指數型 函數的有界變換
X
σ
=
L
1
(
e
−
σ
|
t
|
d
t
)
{\displaystyle X_{\sigma }=L^{1}(e^{-\sigma |t|}dt)}
,由局部可積函數組成,其範數
‖
f
‖
=
∫
−
∞
∞
|
f
(
t
)
|
e
−
σ
|
t
|
d
t
{\displaystyle \|f\|=\int _{-\infty }^{\infty }|f(t)|e^{-\sigma |t|}\,dt}
是有限的。對於f ∈ X σ , I α f 的拉普拉斯變換 形式特別簡單
(
L
I
α
f
)
(
s
)
=
s
−
α
F
(
s
)
{\displaystyle ({\mathcal {L}}I^{\alpha }f)(s)=s^{-\alpha }F(s)}
其中Re(s ) > σ 。這裡F (s ) 表示f 的拉普拉斯變換,該性質表示I α 是傅立葉乘子 。
f 的分數階導數也可以定義為
d
α
d
x
α
f
=
def
d
⌈
α
⌉
d
x
⌈
α
⌉
I
⌈
α
⌉
−
α
f
{\displaystyle {\frac {d^{\alpha }}{dx^{\alpha }}}f\,{\overset {\text{def}}{=}}{\frac {d^{\lceil \alpha \rceil }}{dx^{\lceil \alpha \rceil }}}I^{\lceil \alpha \rceil -\alpha }f}
其中⌈ · ⌉ 表示上限函數 。通過定義
D
x
α
f
(
x
)
=
{
d
⌈
α
⌉
d
x
⌈
α
⌉
I
⌈
α
⌉
−
α
f
(
x
)
α
>
0
f
(
x
)
α
=
0
I
−
α
f
(
x
)
α
<
0.
{\displaystyle D_{x}^{\alpha }f(x)={\begin{cases}{\frac {d^{\lceil \alpha \rceil }}{dx^{\lceil \alpha \rceil }}}I^{\lceil \alpha \rceil -\alpha }f(x)&\alpha >0\\f(x)&\alpha =0\\I^{-\alpha }f(x)&\alpha <0.\end{cases}}}
可以得到微分和積分之間的微分積分插值。
1967 年,卡普托 (Caputo) 提出了一種分數階導數, [ 7] 由此產生的導數具有與前者不同的性質:這種導數對常數函數求導得到零,更重要的是,拉普拉斯變換 的初值項是用該函數的值及其整數階導數來表示的,而不是像黎曼-劉維爾導數那樣用分數階導數來表示。 [ 8] 以x 為基點的 Caputo 分數階導數為:
D
x
α
f
(
y
)
=
1
Γ
(
1
−
α
)
∫
x
y
f
′
(
y
−
u
)
(
u
−
x
)
−
α
d
u
.
{\displaystyle D_{x}^{\alpha }f(y)={\frac {1}{\Gamma (1-\alpha )}}\int _{x}^{y}f'(y-u)(u-x)^{-\alpha }du.}
另一種表示是:
a
D
~
x
α
f
(
x
)
=
I
⌈
α
⌉
−
α
(
d
⌈
α
⌉
f
d
x
⌈
α
⌉
)
.
{\displaystyle {}_{a}{\tilde {D}}_{x}^{\alpha }f(x)=I^{\lceil \alpha \rceil -\alpha }\left({\frac {d^{\lceil \alpha \rceil }f}{dx^{\lceil \alpha \rceil }}}\right).}
函數f (x ) = x (藍色曲線)的半導數(紫色曲線)與一階導數(紅色曲線)。
動畫顯示導數算子在簡單冪函數y = x 的反導數( α = −1 : y = 1 / 2 x 2 )和導數(α = +1 : y = 1 )之間連續來回震盪。
假設f (x ) 是單項式 ,形式如下:
f
(
x
)
=
x
k
.
{\displaystyle f(x)=x^{k}\,.}
一階導數與通常情形一致:
f
′
(
x
)
=
d
d
x
f
(
x
)
=
k
x
k
−
1
.
{\displaystyle f'(x)={\frac {d}{dx}}f(x)=kx^{k-1}\,.}
重複此過程可得到更一般的結果
d
a
d
x
a
x
k
=
k
!
(
k
−
a
)
!
x
k
−
a
,
{\displaystyle {\frac {d^{a}}{dx^{a}}}x^{k}={\dfrac {k!}{(k-a)!}}x^{k-a}\,,}
用伽馬函數 代替階乘 後,得到
d
a
d
x
a
x
k
=
Γ
(
k
+
1
)
Γ
(
k
−
a
+
1
)
x
k
−
a
,
k
>
0.
{\displaystyle {\frac {d^{a}}{dx^{a}}}x^{k}={\dfrac {\Gamma (k+1)}{\Gamma (k-a+1)}}x^{k-a},\quad \ k>0.}
對於k = 1 和a = 1 / 2
d
1
2
d
x
1
2
x
=
Γ
(
1
+
1
)
Γ
(
1
−
1
2
+
1
)
x
1
−
1
2
=
Γ
(
2
)
Γ
(
3
2
)
x
1
2
=
1
π
2
x
1
2
.
{\displaystyle {\frac {d^{\frac {1}{2}}}{dx^{\frac {1}{2}}}}x={\frac {\Gamma (1+1)}{\Gamma \left(1-{\frac {1}{2}}+1\right)}}x^{1-{\frac {1}{2}}}={\frac {\Gamma (2)}{\Gamma \left({\frac {3}{2}}\right)}}x^{\frac {1}{2}}={\frac {1}{\frac {\sqrt {\pi }}{2}}}x^{\frac {1}{2}}.}
為了證明這就是真正的「半導數」(其中H 滿足H 2 f (x ) = Df (x ) ),我們重複該過程得到:
d
1
2
d
x
1
2
2
x
1
2
π
=
2
π
Γ
(
1
+
1
2
)
Γ
(
1
2
−
1
2
+
1
)
x
1
2
−
1
2
=
2
π
Γ
(
3
2
)
Γ
(
1
)
x
0
=
2
π
2
x
0
π
=
1
,
{\displaystyle {\dfrac {d^{\frac {1}{2}}}{dx^{\frac {1}{2}}}}{\dfrac {2x^{\frac {1}{2}}}{\sqrt {\pi }}}={\frac {2}{\sqrt {\pi }}}{\dfrac {\Gamma (1+{\frac {1}{2}})}{\Gamma ({\frac {1}{2}}-{\frac {1}{2}}+1)}}x^{{\frac {1}{2}}-{\frac {1}{2}}}={\frac {2}{\sqrt {\pi }}}{\frac {\Gamma \left({\frac {3}{2}}\right)}{\Gamma (1)}}x^{0}={\frac {2{\frac {\sqrt {\pi }}{2}}x^{0}}{\sqrt {\pi }}}=1\,,}
(因為
Γ
(
3
2
)
=
π
2
{\textstyle \Gamma \!\left({\frac {3}{2}}\right)={\frac {\sqrt {\pi }}{2}}}
且Γ(1) = 1 ),這確實是預期的結果
(
d
1
2
d
x
1
2
d
1
2
d
x
1
2
)
x
=
d
d
x
x
=
1
.
{\displaystyle \left({\frac {d^{\frac {1}{2}}}{dx^{\frac {1}{2}}}}{\frac {d^{\frac {1}{2}}}{dx^{\frac {1}{2}}}}\right)\!x={\frac {d}{dx}}x=1\,.}
對於負整數冪k ,1/
Γ
{\textstyle \Gamma }
為 0,因此使用以下關係很方便: [ 9]
d
a
d
x
a
x
−
k
=
(
−
1
)
a
Γ
(
k
+
a
)
Γ
(
k
)
x
−
(
k
+
a
)
for
k
≥
0.
{\displaystyle {\frac {d^{a}}{dx^{a}}}x^{-k}=\left(-1\right)^{a}{\dfrac {\Gamma (k+a)}{\Gamma (k)}}x^{-(k+a)}\quad {\text{ for }}k\geq 0.}
上述微分算子的擴展不僅限於實數冪;它也適用於複數冪。例如, (1 + i ) 次導數的(1 − i ) 次導數得出二階導數。將a 設為負值會得到積分。
對於一般函數f (x ) 且0 < α < 1 ,完全分數階導數為
D
α
f
(
x
)
=
1
Γ
(
1
−
α
)
d
d
x
∫
0
x
f
(
t
)
(
x
−
t
)
α
d
t
.
{\displaystyle D^{\alpha }f(x)={\frac {1}{\Gamma (1-\alpha )}}{\frac {d}{dx}}\int _{0}^{x}{\frac {f(t)}{\left(x-t\right)^{\alpha }}}\,dt.}
對於任意α ,由於伽馬函數對於負(實)整數是無限的,因此有必要在進行整數導數之後應用分數階導數。例如,
D
3
2
f
(
x
)
=
D
1
2
D
1
f
(
x
)
=
D
1
2
d
d
x
f
(
x
)
.
{\displaystyle D^{\frac {3}{2}}f(x)=D^{\frac {1}{2}}D^{1}f(x)=D^{\frac {1}{2}}{\frac {d}{dx}}f(x).}
我們也可以通過拉普拉斯變換 來解決這個問題。已知
L
{
J
f
}
(
s
)
=
L
{
∫
0
t
f
(
τ
)
d
τ
}
(
s
)
=
1
s
(
L
{
f
}
)
(
s
)
{\displaystyle {\mathcal {L}}\left\{Jf\right\}(s)={\mathcal {L}}\left\{\int _{0}^{t}f(\tau )\,d\tau \right\}(s)={\frac {1}{s}}{\bigl (}{\mathcal {L}}\left\{f\right\}{\bigr )}(s)}
和
L
{
J
2
f
}
=
1
s
(
L
{
J
f
}
)
(
s
)
=
1
s
2
(
L
{
f
}
)
(
s
)
{\displaystyle {\mathcal {L}}\left\{J^{2}f\right\}={\frac {1}{s}}{\bigl (}{\mathcal {L}}\left\{Jf\right\}{\bigr )}(s)={\frac {1}{s^{2}}}{\bigl (}{\mathcal {L}}\left\{f\right\}{\bigr )}(s)}
等等,我們斷言
J
α
f
=
L
−
1
{
s
−
α
(
L
{
f
}
)
(
s
)
}
{\displaystyle J^{\alpha }f={\mathcal {L}}^{-1}\left\{s^{-\alpha }{\bigl (}{\mathcal {L}}\{f\}{\bigr )}(s)\right\}}
。
例如,
J
α
(
t
k
)
=
L
−
1
{
Γ
(
k
+
1
)
s
α
+
k
+
1
}
=
Γ
(
k
+
1
)
Γ
(
α
+
k
+
1
)
t
α
+
k
{\displaystyle J^{\alpha }(t^{k})={\mathcal {L}}^{-1}\left\{{\frac {\Gamma (k+1)}{s^{\alpha +k+1}}}\right\}={\frac {\Gamma (k+1)}{\Gamma (\alpha +k+1)}}t^{\alpha +k}}
正如預期的那樣。事實上,考慮到卷積 規則
L
{
f
∗
g
}
=
(
L
{
f
}
)
(
L
{
g
}
)
{\displaystyle {\mathcal {L}}\{f*g\}={\bigl (}{\mathcal {L}}\{f\}{\bigr )}{\bigl (}{\mathcal {L}}\{g\}{\bigr )}}
並且為了清楚起見,我們將p (x ) 簡寫為p (x ) = x α − 1 ,由此可得
(
J
α
f
)
(
t
)
=
1
Γ
(
α
)
L
−
1
{
(
L
{
p
}
)
(
L
{
f
}
)
}
=
1
Γ
(
α
)
(
p
∗
f
)
=
1
Γ
(
α
)
∫
0
t
p
(
t
−
τ
)
f
(
τ
)
d
τ
=
1
Γ
(
α
)
∫
0
t
(
t
−
τ
)
α
−
1
f
(
τ
)
d
τ
{\displaystyle {\begin{aligned}\left(J^{\alpha }f\right)(t)&={\frac {1}{\Gamma (\alpha )}}{\mathcal {L}}^{-1}\left\{{\bigl (}{\mathcal {L}}\{p\}{\bigr )}{\bigl (}{\mathcal {L}}\{f\}{\bigr )}\right\}\\&={\frac {1}{\Gamma (\alpha )}}(p*f)\\&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{t}p(t-\tau )f(\tau )\,d\tau \\&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{t}\left(t-\tau \right)^{\alpha -1}f(\tau )\,d\tau \\\end{aligned}}}
這就是上面柯西推導的式子。
拉普拉斯變換對相對較少的函數起作用,但它們通常對求解分數階微分方程很有用 。
^ 1.0 1.1 Lizorkin 2001 harvnb模板錯誤: 無指向目標: CITEREFLizorkin2001 (幫助 )
^ Liouville, Joseph , Mémoire sur quelques questions de géométrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces questions , Journal de l'École Polytechnique (Paris), 1832, 13 : 1–69 .
^ Liouville, Joseph , Mémoire sur le calcul des différentielles à indices quelconques , Journal de l'École Polytechnique (Paris), 1832, 13 : 71–162 .
^ Riemann, Georg Friedrich Bernhard , Versuch einer allgemeinen Auffassung der integration und differentiation , Weber, H. (編), Gesammelte Mathematische Werke, Leipzig, 1896 [1847] .
^ Brychkov & Prudnikov 2001 harvnb模板錯誤: 無指向目標: CITEREFBrychkovPrudnikov2001 (幫助 )
^ Miller & Ross 1993
^ Caputo 1967
^ Loverro 2004
^ Bologna, Mauro, Short Introduction to Fractional Calculus (PDF) , Universidad de Tarapaca, Arica, Chile, [2014-04-06 ] , (原始內容 (PDF) 存檔於2016-10-17)
Brychkov, Yu.A.; Prudnikov, A.P., Euler transformation , Hazewinkel, Michiel (編), 数学百科全书 , Springer , 2001, ISBN 978-1-55608-010-4 .
Caputo, Michele, Linear model of dissipation whose Q is almost frequency independent. II, Geophysical Journal International, 1967, 13 (5): 529–539, Bibcode:1967GeoJ...13..529C , doi:10.1111/j.1365-246x.1967.tb02303.x .
Hille, Einar; Phillips, Ralph S., Functional analysis and semi-groups, Providence, R.I.: American Mathematical Society , 1974, MR 0423094 .
Hazewinkel, Michiel (編), Fractional integration and differentiation , 数学百科全书 , Springer , 2001, ISBN 978-1-55608-010-4 .
Loverro, Adam, Fractional Calculus: History, Definitions and Applications for the Engineer (PDF) , Notre Dame, IN: University of Notre Dame, 2004-05-08, (原始內容 (PDF) 存檔於2005-10-29)
Miller, Kenneth S.; Ross, Bertram, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, 1993, ISBN 0-471-58884-9 .
Riesz, Marcel , L'intégrale de Riemann-Liouville et le problème de Cauchy, Acta Mathematica , 1949, 81 (1): 1–223, ISSN 0001-5962 , MR 0030102 , doi:10.1007/BF02395016 .