Tempered分数阶微分方程Nagumo型唯一解
Nagumo-Type Uniqueness Result for Tempered Fractional Differential Equations
DOI: 10.12677/PM.2023.1311339, PDF, HTML, XML, 下载: 98  浏览: 137 
作者: 黄洪鸿:福建师范大学数学与统计学院,福建 福州
关键词: Tempered分数阶算子Nagumo型唯一性迭代法Tempered Fractional Operators Nagumo-Type Uniqueness Iterative Approach
摘要: 本文重点讨论Tempered分数阶微分方程柯西问题Nagumo迭代近似值的唯一性和收敛性。首先把柯西问题转化成等价Volterra积分方程,证出解的唯一性。然后,运用迭代方法,我们将Nagumo型唯一性和迭代近似值序列扩展到Tempered分数阶微分方程。
Abstract: This work is concerned with Nagumo-type uniqueness and convergence of successive approxima-tions to Cauchy problem for Tempered fractional differential equations. Firstly, to prove the uniqueness of the solution, the Cauchy problem is transformed into an equivalent Volterra integral equation. Then, using the iterative method, we extend Nagumo-type uniqueness and convergence of successive approximations to Tempered fractional differential equations.
文章引用:黄洪鸿. Tempered分数阶微分方程Nagumo型唯一解[J]. 理论数学, 2023, 13(11): 3254-3261. https://doi.org/10.12677/PM.2023.1311339

1. 引言

分数阶导数是一个重要的数学工具,它具有一些非常有用的性质。分数阶导数具有非局域性,也就是说,它能够反映物理量在整个空间范围内的变化情况。相比于整数阶导数,分数阶导数可以更好地描述非局部的物理现象,如扩散、传输和耗散等。分数阶导数发展出几种不同类型的分数阶导数:如Riemann-Liouville分数阶、Caputo型分数阶导数( [1] )、Riesz型分数阶导数( [2] )以及Hilfer分数阶导数( [3] )。现在广泛使用的一种分数阶导数的定义是Riemann-Liouville分数阶积分与导数。函数在区间 [ 0 , T ] 0 < α < 1 阶Riemann-Liouville分数阶积分与导数的定义为( [4] [5] ):

I 0 + α x ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 x ( s ) d s

D x α ( t ) = D I 0 + 1 α x ( t )

其中 D = d d t 。此后,分数阶微积分开始在各个领域得到应用。例如,在物理学中,分数阶微积分被用于

描述扩散过程、黏弹性材料、地震等。在工程学中,它被用于控制系统的稳定性分析、信号处理等。这类方程能更准确地描述某些现象,因为它们考虑了记忆效应和遗传效应。

分数阶导数在物理背景与实验数据进一步吻合的过程中修正和发展,其中就拓展出Tempered分数阶导数与积分。据我们所知,Buschman在文献 [6] 中首次提出具有弱奇异核和指数核的分数阶积分的定义。目前为止,Tempered型分数阶导数在基本性质、初值问题、边值问题和中值问题解的存在性、解的数值计算等等方面取得部分成果,其中Tempered分数阶微分方程的柯西问题有助于我们更深入地理解分数阶微分方程的性质和行为,这包括了解的唯一性、稳定性、收敛性等。Tempered分数阶微分方程可以描述许多具有记忆和遗传特性的复杂现象,它在处理某些具有特定衰减特性的问题时非常有用。例如材料科学、生物医学、工程控制等领域。研究Nagumo迭代近似值的唯一性和收敛性可以为实际问题的求解提供有效方法和理论支持,从而提高求解分数阶微分方程的效率和精度。因此,本文将Nagumo型迭代近似值的唯一性和收敛性从分数阶微分方程推广到Tempered分数阶微分方程。

例如,在 [7] 中,J. Chu和Z. W证明了半无限区间上非线性二阶微分方程

{ x = f ( t , x ) , t t 0 0 x ( t 0 ) = x 0 lim t x ( t ) = ψ 0 , lim t x ( t ) exp ( t ) = 0 , ψ 0 (1.1)

Nagumo型唯一性解。

后来,V. Lakshmikantham和S. Leela把Nagumo型迭代近似值的唯一性和收敛性推广到分数阶微分方程:

{ D α x = f ( t , x ) , 0 < α < 1 x ( t 0 ) = x 0 = x ( t ) ( t t 0 ) 1 α | t = t 0 . (1.2)

受此启发,本文证明了Tempered分数阶微分方程柯西问题Nagumo迭代近似解的唯一性和收敛性。首先,通过将柯西问题转化为等价的积分方程,证明了该方程解的唯一性。然后,利用迭代方法,证明了迭代Nagumo近似值序列的唯一性和收敛性,并收敛于该积分方程的唯一解。这些结果表明,对于Tempered分数阶微分方程柯西问题,可以通过迭代方法得到唯一且收敛的近似解。

本文结构如下:第二节介绍了Tempered分数阶导数的定义(定义2.1)以及Tempered分数阶微分方程的等价Volterra积分方程(引理2.2),第三节介绍了本文的主要结果Tempered分数阶微分方程Nagumo型迭代近似值的唯一性和收敛性。

2. 预备知识

本节给出Tempered分数阶导数定义和引理。以下均假设 0 < q < 1 ,参数 λ 0 为正数,设可积空间 L ( [ a , b ] ) ,其中包括有限区间 [ a , b ] 上的勒贝格可测函数,

L ( [ a , b ] ) = { u : u L ( [ a , b ] ) = a b | x ( t ) | d t < } .

定义2.1 ( [8] [9] ):设 x ( t ) 是定义在 [ 0 , T ] 上的实值函数,其q阶Riemann-Liouville Tempered分数阶导数定义为

a D t q , λ x ( t ) = e λ t D a t q ( e λ t x ( t ) ) = e λ t Γ ( 1 q ) d d t t a e λ s x ( s ) ( t s ) q d s , (2.1)

其中 a D t q ( e λ t x ( t ) ) 为Riemann-Liouville分数阶导数定义 [5]

a D t q ( e λ t x ( t ) ) = d d t ( a I t 1 q ( e λ t x ( t ) ) ) = 1 Γ ( 1 q ) d d t a t ( e λ s x ( s ) ) ( t s ) α d s . (2.2)

为了证明定理1.1,还需以下引理,有兴趣的读者可参考文献 [10] 中引理1的证明。

引理2.2

{ a D t q , λ x ( t ) = f ( t , x ( t ) ) , 0 < q < 1 , λ 0 [ a D t q 1 ( e λ t x ( t ) ) ] | t = a = g 0 (2.3)

若函数 f ( t , x ( t ) ) x ( t ) 属于 L ( [ a , b ] ) x ( t ) 是柯西问题(2.3)的解,问题(2.3)等价Volterra积分方程:

x ( t ) = g 0 e λ t ( t a ) q 1 Γ ( q ) + 1 Γ ( q ) a t e λ ( t s ) ( t s ) q 1 f ( s , x ( s ) ) d s .

3. 主要结果及证明

在上述定义和引理的基础上,本文考虑分数阶微分方程(2.1)迭代近似值的唯一性和收敛性。

{ D q , λ x ( t ) = f ( t , x ( t ) ) , f C ( R 0 , R ) x ( t 0 ) = x ( t ) ( t t 0 ) 1 q | t = t 0 = x 0 (3.1)

其中 R 0 = { ( t , x ) , t 0 t t 0 + a , | x x 0 ( t ) | b } D q , λ 表示为q阶Riemann-Liouville Tempered分数阶导数(见定义2.1)。由引理2.1可知柯西问题(3.1)等价于第二类Volterra积分方程:

x ( t ) = x 0 ( t ) + 1 Γ ( q ) t 0 t e λ ( t s ) ( t s ) q 1 f ( s , x ( s ) ) d s (3.2)

其中 x 0 ( t ) = x 0 e λ t ( t t 0 ) q 1 Γ ( q ) 。因此(3.1)初值问题的每个解等价于(3.2) Volterra积分方程的解,反之亦然。

我们将推广到Tempered分数阶微分方程(3.1) Nagumo型解的存在性与唯一性。

定理3.1假设函数 x ( t ) C p ( [ t 0 , t 0 + a ] , R ) = { u : u C ( [ t 0 , t 0 + a ] , R ) } 是柯西问题(3.1)的一个解, ( t t 0 ) p u ( t ) C ( [ t 0 , t 0 + a ] , R ) , p = 1 q D q , λ 在区间 [ t 0 , t 0 + a ] 内有意义且连续。那么假设(3.1)中f满足Nagumo条件:

| f ( t , x ) f ( t , y ) | | x y | Γ ( q ) ( t t 0 ) q , t t 0 , ( t , x ) R 0 . (3.3)

那么迭代近似值

x n + 1 = x 0 ( t ) + 1 Γ ( q ) t 0 t e λ ( t s ) ( t s ) q 1 f ( s , x n ( s ) ) d s (3.4)

收敛于(3.1)在 [ t 0 , t 0 + η ] 的唯一解 x ( t ) ,其中 η = min ( a , ( b Γ ( 1 + q ) M ) ) 1 q ,在 R 0 | f ( t , x ) | M

证明:第一步,要证(3.1)解 x ( t ) 的唯一性。首先假设 x ( t ) , y ( t ) 是(3.1)的解。那么 ϕ ( t ) = | x ( t ) y ( t ) | 。因为 ϕ ( t 0 ) = 0 ,可得

lim t t 0 + ϕ ( t ) t t 0 = lim t t 0 + ( 1 t t 0 1 Γ ( q ) t 0 t e λ ( t s ) ( t s ) q 1 [ f ( s , x ( s ) ) f ( s , y ( s ) ) ] d s ) = 1 Γ ( q ) [ ( t t 0 ) q 1 e λ ( t t 0 ) [ f ( t 0 , x ( t 0 ) ) f ( t 0 , y ( t 0 ) ) ] ] = 0. (3.5)

故而

lim t t 0 + ϕ ( t ) ( t t 0 ) q e λ ( t t 0 ) = lim t t 0 + ϕ ( t ) ( t t 0 ) ( t t 0 ) 1 q e λ ( t t 0 ) = 0. (3.6)

接着使用Nagumo (3.3)条件得

ϕ ( t ) 1 Γ ( q ) t 0 t ( t s ) q 1 e λ ( t s ) | ϕ ( s ) | ( s t 0 ) q d s Γ ( q ) q . (3.7)

ϕ ( t ) e λ ( t t 0 ) ( t t 0 ) q = ψ ( t ) ,同样地 ψ ( t 0 ) = 0 ,(3.7)可转化为

ψ ( t ) q ( t t 0 ) q t 0 t ( t s ) q ψ ( s ) d s .

为了证(3.1)解的唯一性,等价于证明

ψ ( t ) 0. (3.8)

令(3.8)不成立,则存在一点 t 1 ( t 0 , t 0 + η ) ,使得 m = max [ t 0 , t 0 + η ] ψ ( t ) = ψ ( t 1 )

那么

m = ψ ( t 1 ) q e λ ( t 1 t 0 ) ( t 1 t 0 ) q t 0 t 1 ( t 1 s ) q ψ ( s ) d s < q m ( t 1 t 0 ) q t 0 t 1 ( t 1 s ) q d s = q m ( t 1 t 0 ) q ( t 1 s ) q q | t 0 t 1 = q m ( t 1 t 0 ) q ( t 1 t 0 ) q q = m . (3.9)

这证得了 ψ ( t ) 0 ,意味着 ϕ ( t ) 0 ,即(3.1)柯西问题有唯一解。

现在我们要证逐次近似值 x n + 1 ( t ) n = 0 , 1 , 2 , 满足:

(i) 是良定义和连续的;

(ii) 是一致有界的;

(iii) 在 [ t 0 , t 0 + η ] 上等度连续的。

其中 η = min ( a , ( b Γ ( 1 + q ) M ) ) 1 q ,以及在 R 0 | f ( t , x ) | M

现在由(3.4)的定义可知,

| x n + 1 ( t ) x 0 ( t ) | t 0 t ( t s ) q 1 | f ( x , x n ( s ) ) | e λ ( t s ) d s .

n = 0 时,结合 | f ( t , x ) | M 上式可化为

| x 1 ( t ) x 0 ( t ) | C M Γ ( q ) ( t s ) q q | t 0 t = C M ( t t 0 ) q q Γ ( q ) = C M a q q Γ ( q ) b

其中C是常数。

通过数学归纳法可得 { x n + 1 ( t ) } 符合(ii)条件,即函数族 { x n + 1 ( t ) } [ t 0 , t 0 + η ] 是一致有界的。

对于 t 0 t 1 t 2 ,结合 x n ( t ) 的定义和 | f ( t , x ) | M ,我们有

| x n ( t 1 ) x n ( t 2 ) | 1 Γ ( q ) [ | t 0 t 1 ( t 1 s ) q 1 e λ ( t 1 s ) f ( s , x n ( s ) ) d s t 0 t 2 ( t 2 s ) q 1 e λ ( t 2 s ) f ( s , x n ( s ) ) d s | ] = 1 Γ ( q ) [ | t 0 t 1 ( ( t 1 s ) q 1 e λ ( t 1 s ) ( t 2 s ) q 1 e λ ( t 2 s ) ) f ( s , x n ( s ) ) d s t 1 t 2 ( t 2 s ) q 1 e λ ( t 2 s ) f ( s , x n ( s ) ) d s | ]

M Γ ( q ) [ t 0 t 1 ( ( t 1 s ) q 1 ( t 2 s ) q 1 ) d s + t 1 t 2 ( t 2 s ) q 1 d s ] = M Γ ( q + 1 ) [ ( t 1 t 0 ) q ( t 2 t 0 ) q + 2 ( t 2 t 1 ) q ] 2 M Γ ( q + 1 ) ( t 2 t 1 ) q . (3.10)

| t 2 t 1 | < δ = ( ϵ Γ ( q + 1 ) 2 M ) 1 q ,对于n足够大,则(3.10)满足小于 ϵ ,即

| x n ( t 1 ) x n ( t 2 ) | < ϵ .

故而容易得函数族 { x n + 1 ( t ) } [ t 0 , t 0 + η ] 是等度连续的,因此存在一致收敛的子列 { x n k ( t ) } 。于是,假设当 n 时, { x n ( t ) x n + 1 ( t ) } 0 ,则(3.4)说明任何这样的子列的极限都是(3.1)的唯一解 x ( t ) 。在此基础上,由文献 [4] 可知选择一个子列是不必要的,并且完整的序列 x n + 1 ( t ) 一致收敛于 x ( t ) 。因此,我们可以得到

ϕ ( t ) = lim sup n ( | x n + 1 ( t ) x n ( t ) | ) 0. (3.11)

对于 t 0 t 3 t 4 ,同理可得:

| x n ( t 3 ) x n 1 ( t 3 ) | | x n ( t 4 ) x n 1 ( t 4 ) | | x n ( t 3 ) x n 1 ( t 3 ) x n ( t 4 ) + x n 1 ( t 4 ) | 1 Γ ( q ) [ | t 0 t 3 ( t 3 s ) q 1 e λ ( t 3 s ) D ( S ) d s t 0 t 4 ( t 4 s ) q 1 e λ ( t 4 s ) D ( S ) d s | ] = 1 Γ ( q ) [ | t 0 t 3 ( ( t 3 s ) q 1 e λ ( t 3 s ) ( t 4 s ) q 1 e λ ( t 4 s ) ) D ( s ) d s t 3 t 4 ( t 4 s ) q 1 e λ ( t 4 s ) D ( s ) d s | ]

2 M Γ ( q ) [ t 0 t 3 ( ( t 3 s ) q 1 ( t 4 s ) q 1 ) d s + t 3 t 4 ( t 4 s ) q 1 d s ] = 2 M Γ ( q + 1 ) [ ( t 3 t 0 ) q ( t 4 t 0 ) q + 2 ( t 4 t 3 ) q ] 4 M Γ ( q + 1 ) ( t 4 t 3 ) q . (3.12)

其中在 R 0 | D ( s ) | = | f ( s , x n 1 ( s ) ) f ( s , x n 2 ( t ) ) | 2 M ,以及 x n ( t ) 与(3.4)中定义相同。

| t 4 t 3 | < δ = ( ϵ Γ ( q + 1 ) 4 M ) 1 q ,则(3.12)满足小于 ϵ 。对于n足够大, ϵ > 0 ,(3.12)可化为

因为 t 3 , t 4 可任意替换,以及 ϵ > 0 是任意,故而

| ϕ ( t 3 ) ϕ ( t 4 ) | 4 M Γ ( q + 1 ) ( t 4 t 3 ) q .

由此证得 ϕ ( t ) 的连续性。

接下来,结合Nagumo条件和(3.4)中 x n + 1 ( t ) 的定义,可得

| x n + 1 ( t ) x n ( t ) | q Γ ( q ) Γ ( q ) t 0 t ( t s ) q 1 e λ ( t s ) | x n ( s ) x n 1 ( s ) | ( s t 0 ) q d s (3.13)

任意固定 t [ t 0 , t 0 + η ] ,则存在一个整数序列 n 1 < n 2 < ,使得当 n = n k

| x n + 1 ( t ) x n ( t ) | ϕ ( t )

ϕ * ( t ) = lim n = n k | x n ( t ) x n 1 ( t ) | 一致存在于 t [ t 0 , t 0 + η ]

因此,

ϕ ( t ) q t 0 t ( t s ) q 1 e λ ( t s ) ϕ * ( t ) ( s t 0 ) q d s . (3.14)

由于 ϕ * ( t ) ϕ ( t ) ,(3.14)放缩为

ϕ ( t ) q t 0 t ( t s ) q 1 e λ ( t s ) ϕ ( t ) ( s t 0 ) q d s . (3.15)

与上文证明相似,我们可以得出

lim t t 0 + ϕ ( t ) t t 0 = lim t t 0 + ( 1 t t 0 1 Γ ( q ) t 0 t e λ ( t s ) ( t s ) q 1 [ f ( s , x n ( s ) ) f ( s , x n 1 ( s ) ) ] d s ) = 1 Γ ( q ) [ ( t t 0 ) q 1 e λ ( t t 0 ) [ f ( t 0 , x n ( t 0 ) ) f ( t 0 , x n 1 ( t 0 ) ) ] ] = 0.

更进一步可得,

lim t t 0 + ϕ ( t ) ( t t 0 ) q e λ ( t t 0 ) = lim t t 0 + ϕ ( t ) ( t t 0 ) ( t t 0 ) 1 q e λ ( t t 0 ) = 0.

假设 ψ ( t ) = ϕ ( t ) ( t t 0 ) q e λ ( t t 0 ) ,从上式可得 ψ ( t 0 ) = 0 。接下来,为了证 x n + 1 ( t ) 的收敛性,等价于证

ψ ( t ) 0. (3.16)

假设(3.16)不成立,则存在 ψ ( t 1 ) = max [ t 0 , t 0 + η ] ψ ( t ) = m 。那么

m = ψ ( t 1 ) = ϕ ( t 1 ) ( t 1 t 0 ) q e λ ( t 1 t 0 ) q e λ ( t 1 t 0 ) ( t 1 t 0 ) q t 0 t 1 ( t 1 s ) q ψ ( s ) d s < q m ( t 1 t 0 ) q t 0 t 1 ( t 1 s ) q d s = q m ( t 1 t 0 ) q ( t 1 s ) q q | t 0 t 1 = q m ( t 1 t 0 ) q ( t 1 t 0 ) q q = m . (3.17)

显然矛盾,(3.16)成立,即逐次近似值 x n + 1 ( t ) (见(3.4))收敛于(3.1)的唯一解。定理3.1成立。

参考文献

[1] Podlubny, I. (1999) Fractional Differential Equations. Academic Press, San Diego.
[2] Samko, S., Kilbas, A. and Marichev, O. (1993) Fractional Integrals and Derivatives: Theory and Applications. Gordon and Breach, Lon-don.
[3] Hilfer, R. (2000) Applications of Fractional Calculus in Physics. World Scientific, Singapore.
https://doi.org/10.1142/3779
[4] Samko, S.G., Kilbas, A.A. and Marichev, O.I. (1993) Fractional Integrals and Derivatives. Gordon and Breach Science Publishers, Switzerland, Philadelphia.
[5] Diethelm, K. (2004) The Analysis of Fractional Differential Equations. Journal of Mathematical Analysis and Applications, 256, 229-248.
[6] Buschman, R.G. (1972) Decomposition of an Integral Operator by Use of Mikusinski Calculus. SIAM Journal on Mathematical Analysis, 3, 83-85.
https://doi.org/10.1137/0503010
[7] Chu, J. and Wang, Z. (2023) Nagumo-Type Uniqueness and Stability for Nonlinear Differential Equations on Semi-Infinite Intervals. Journal of Differential Equations, 367, 229-245.
https://doi.org/10.1016/j.jde.2023.05.001
[8] Baeumera, B. and Meerschaert, M. (2010) Tempered Stable Levy Motion and Transient Super-Diffusion. Journal of Computational and Applied Mathematics, 233, 2438-2448.
https://doi.org/10.1016/j.cam.2009.10.027
[9] Cartea, A. and del-Castillo-Negrete, D. (2007) Fluid Limit of the Continuous-Time Random Walk with General Levy Jump Distribution Functions. Physical Review E, 76, Article ID: 041105.
https://doi.org/10.1103/PhysRevE.76.041105
[10] Li, C., Deng, W. and Zhao, L. (2015) Well-Posedness and Numerical Algorithm for the Tempered Fractional Ordinary Differential Equations. arXiv:1501.00376.