一类Caputo-Fabrizio脉冲分数阶微分方程解的存在性
Existence of Solutions for a Class of Caputo-Fabrizio Impulsive Fractional Differential Equations
摘要: 在Banach空间中研究一类具有Caputo-Fabrizio分数阶导数且在非局部条件下的脉冲分数阶微分方程解的存在性。利用Schaefer不动点定理,压缩映射原理,Arzela-Ascoli定理,得到了该脉冲分数阶问题至少一个解和唯一解,并用一个例子验证其中一个结论。
Abstract: The existence of a class of impulsive fractional differential equations with Caputo-Fabrizio fractional derivatives under non-local conditions is studied in Banach space. Based on Schaefer’s fixed point theorem, compression mapping principle and Arzela-Ascoli theorem, at least one solution and the only solution of the pulse fractional order problem are obtained, and one of the conclusions is verified by an example.
文章引用:张德霞, 王仲平. 一类Caputo-Fabrizio脉冲分数阶微分方程解的存在性[J]. 应用数学进展, 2024, 13(12): 5120-5128. https://doi.org/10.12677/aam.2024.1312494

1. 引言

分数阶微分方程在描述物理、生物、经济等领域的非整数阶动力系统中有着广泛的应用[1] [2]。近年来,随着分数阶微分方程理论的发展,越来越多的学者开始关注分数阶微分方程。这类方程在描述具有脉冲效应的系统时,如金融市场、人口动态等,具有显著优势。

2015年,Caputo和Fabrizio [3]定义了一种没有奇异核的新的分数阶导数,在经典的Caputo分数阶导数

D C α f( t )= 1 Γ( 1α ) a t ( ts ) α f ( s )ds

基础上,用 M( α ) 1α 替换 1 Γ( 1α ) exp( α ts 1α ) 替换 ( ts ) α ,得到Caputo-Fabrizio分数阶导数。

同年,Losada和Nieto [4]提出了与Caputo-Fabrizio分数阶导数相对应的分数阶积分,并证明了Caputo-Fabrizio线性分数阶微分方程

{ D CF α x( t )=λx( t )+u( t ),t0,λ, x( 0 )= x 0 .

解的存在性。

2021年,Abbas [5]利用压缩映射的拓扑度不动点定理,研究了Caputo-Fabrizio型分数阶导数脉冲分数阶微分方程的初值问题

{ D CF α x( t )=f( t,x( t ) ),tJ=[ 0,T ],t t i , x( 0 )= x 0 , Δx( t i )= I i ( x( t i ) ),i=1,2,,n.

解的存在唯一性。

在某些应用场景中,系统的历史状态对当前状态的影响不能简单地用一个局部条件来刻画。因此,需要引入非局部条件来更准确地描述这种动态行为。文献[6]研究了任意Hilbert空间中一类具有非局部条件的非线性脉冲中立型含Caputo分数阶导数的微分方程,利用解析半群理论和不动点方法证明了一类近似解的存在唯一性;文献[7]利用解算子理论、近似技术、半群理论和不动点理论,研究了非局部和脉冲中立型积分微分方程的近似可控性。文献[8]研究了包含Caputo分数阶导数的分数阶微分方程柯西问题

{ D C α x( t )=f( t,x( t ) ),tJ:=[ 0,T ], x( 0 )+g( x )= x 0 .

解的存在唯一性和数据依赖性。

基于上述工作,本文研究了在Banach空间中包含Caputo-Fabrizio分数阶导数且具有非局部条件的下列脉冲分数阶微分方程

{ D CF α x( t )=f( t,x( t ) ),tJ=[ 0,T ],t t i , Δx( t i )= I i ( x( t i ) ),i=1,2,,n, x( 0 )+g( x )= x 0 . (1)

解的存在性,其中 D CF t α x( t ) 0<α<1 阶的Caputo-Fabrizio分数阶导数,函数 f:J×XX 连续, Δx( t i )=x( t i + )x( t i ) x( t i ),x( t i + ) 表示 x= t i 处的左、右极限, I i 是脉冲函数,非局部项 g:PC( J,X )X 是连续映射, PC( J,X ) 是所有从JX的连续函数的Banach空间。对于 xPC( J,X ) ,其范数为 x :=sup{ | x( t ) |:tJ }

2. 预备知识

定义1 [4] 0<α<1 f C 1 [ 0,b ) b>0 ,Caputo-Fabrizio分数阶导数的定义为

D CF α f( t )= ( 2α )M( α ) 2( 1α ) 0 t exp( α( ts ) 1α ) f ( s )ds .

其中 M( α ) 是一个依赖于 α 的归一化常数,当且仅当f是常数函数时, D CF α f( t )=0

定义2 [4] 0<α<1 ,函数f α 阶Caputo-Fabrizio分数阶积分的定义为

I CF α f( t )= 2( 1α ) ( 2α )M( α ) f( t )+ 2α ( 2α )M( α ) 0 t f( s )ds ,t0.

根据前面的定义,函数f 0<α<1 的Caputo-Fabrizio型分数阶积分是函数f与其一阶积分的平均值。

2( 1α ) ( 2α )M( α ) + 2α ( 2α )M( α ) =1.

由此可得出结论

M( α )= 2 2α ,0α1.

定理1 [9] (Banach压缩映射原理)设X是完备的度量空间, T:XX 是一个压缩映射,那么T有且只有一个不动点。

定理2 [10] (Schaefer不动点定理)设X是Banach空间, Q:XX 是一个完全连续算子,如果集合

S( Q )={ xX:x=λQx,0<λ<1 }

有界,则Q有一个不动点。

定理3 [11] (Arzela-Ascoli定理)设X紧, HC( X ) ,则H相对紧,当且仅当H有界且等度连续。

3. 主要结论

引理1 [4]函数 xPC( J,X ) 是积分方程

x( t )= x 0 + c α ( σ( t )σ( 0 ) )+ d α 0 t σ( s )ds ,t0,

的解,当且仅当x是以下方程

{ D CF α x( t )=σ( t ),0<α<1,t0, x( 0 )= x 0 X.

的解,其中 c α = 2( 1α ) ( 2α )M( α ) =1α d α = 2α ( 2α )M( α ) =α

引理2函数 xPC( J,X ) 是积分方程

x( t )={ x 0 g( x )+ c α [ f( t,x( t ) )f( 0,x( 0 ) ) ]+ d α 0 t f( s,x( s ) )ds ,t[ 0, t 1 ], x 0 g( x )+ 0< t i <t I i ( x( t i ) )+ c α [ f( t,x( t ) )f( 0,x( 0 ) ) ] + d α 0< t i <t 0 t i f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds ,t( t i , t i+1 ],i=1,,n. (2)

的解,当且仅当x是脉冲分数阶微分方程(1)的一个解。

证明 首先设 xPC( J,X ) 是脉冲分数阶微分方程(1)的一个解,那么,对于 t[ 0, t 1 ] D CF α x( t )=f( t,x( t ) ) x( 0 )+g( x )= x 0 ,由引理1,得

x( t )= x 0 g( x )+ c α [ f( t,x( t ) )f( 0,x( 0 ) ) ]+ d α 0 t f( s,x( s ) )ds

对于 t( t 1 , t 2 ] D CF α x( t )=f( t,x( t ) ) x( t 1 + )=x( t 1 )+ I 1 ( x( t 1 ) ) ,由引理1,有

x( t )=x( t 1 + )+ c α [ f( t,x( t ) )f( t 1 ,x( t 1 ) ) ]+ d α 0 t f( s,x( s ) )ds =x( t 1 )+ I 1 ( x( t 1 ) )+ c α [ f( t,x( t ) )f( t 1 ,x( t 1 ) ) ]+ d α 0 t f( s,x( s ) )ds = x 0 g( x )+ I 1 ( x( t 1 ) )+ c α [ f( t 1 ,x( t 1 ) )f( 0,x( 0 ) ) ] + c α [ f( t,x( t ) )f( t 1 ,x( t 1 ) ) ]+ d α 0 t 1 f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds = x 0 g( x )+ I 1 ( x( t 1 ) )+ c α [ f( t,x( t ) )f( 0,x( 0 ) ) ] + d α 0 t 1 f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds

对于 t( t 2 , t 3 ] D CF α x( t )=f( t,x( t ) ) x( t 2 + )=x( t 2 )+ I 2 ( x( t 2 ) ) ,有

x( t )=x( t 2 + )+ c α [ f( t,x( t ) )f( t 2 ,x( t 2 ) ) ]+ d α 0 t f( s,x( s ) )ds =x( t 2 )+ I 2 ( x( t 2 ) )+ c α [ f( t,x( t ) )f( t 2 ,x( t 2 ) ) ]+ d α 0 t f( s,x( s ) )ds = x 0 g( x )+ I 1 ( x( t 1 ) )+ I 2 ( x( t 2 ) )+ c α [ f( t 2 ,x( t 2 ) )f( 0,x( 0 ) ) ] + c α [ f( t,x( t ) )f( t 2 ,x( t 2 ) ) ]+ d α 0 t 1 f( s,x( s ) )ds + d α 0 t 2 f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds = x 0 g( x )+ I 1 ( x( t 1 ) )+ I 2 ( x( t 2 ) )+ c α [ f( t,x( t ) )f( 0,x( 0 ) ) ] + d α 0 t 1 f( s,x( s ) )ds + d α 0 t 2 f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds

对于 t( t i , t i+1 ] ,有

x( t )= x 0 g( x )+ 0< t i <t I i ( x( t i ) ) + c α [ f( t,x( t ) )f( 0,x( 0 ) ) ]          + d α 0< t i <t 0 t i f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds ,i=1,,n.

反之,如果x满足(2),可以证明x是(1)的解。证毕。

定义算子: T:PC( J,X )PC( J,X )

( Tx )( t )= b α g( x )+ 0< t i <t I i ( x( t i ) )+ c α f( t,x( t ) )              + d α 0< t i <t 0 t i f( s,x( s ) )ds + d α 0 t f( s,x( s ) )ds ,t( t i , t i+1 ],i=1,,n.

其中 b α = x 0 c α f( 0,x( 0 ) )

因此,脉冲分数阶微分方程(1)的解的存在性等价于算子T不动点的存在性。

为了证明方便,我们做了以下假设:

(H1) 对 x,yX tJ ,存在 δ f >0 ,使得

f( t,x )f( t,y ) δ f xy ,

( t,x )( J,X ) ,存在 δ 1 >0, α 1 [ 0,1 ) δ 2 >0 使得

f( t,x ) δ 1 x α 1 + δ 2 ,

(H2) 脉冲函数 I i :XX ,存在 γ I [ 0, 1 n ) ,使得

I i ( x ) I i ( y ) γ I xy ,

xX tJ ,存在 γ 1 >0 α 2 [ 0,1 ) γ 2 >0 使得

I i ( x ) γ 1 x α 2 + γ 2 ,

(H3) 对 x,yX tJ ,存在 δ f >0 C g >0 使得

g( x )g( y ) C g xy ,

xPC( J,X ) ,存在 C 1 >0 α 3 [ 0,1 ) C 2 >0 使得

g( x ) C 1 x α 3 + C 2 ,

引理3算子T PC( J,X ) 上连续。

证明 设 { x n } PC( J,X ) 中的一个序列,且 x n x ,对 tJ ,有

( T x n )( t )( Tx )( t ) g( x n )g( x ) + 0< t i <t I i ( x n ( t i ) ) I i ( x( t i ) ) + c α f( t, x n ( t ) )f( t,x( t ) ) + d α 0< t i <t 0 t i f( s, x n ( s ) )f( s,x( s ) ) ds + d α 0 t f( s, x n ( s ) )f( s,x( s ) ) ds C g x n x + 0< t i <t γ I x n x + c α δ f x n x + d α 0< t i <t 0 t i δ f x n x ds + d α 0 t δ f x n x ds ( C g +n γ I + c α δ f + d α nT δ f + d α T δ f ) x n x

n ( T x n )( t )( Tx )( t ) 0 。因此,T是连续的。

引理4算子T PC( J,X ) 上有界。

证明 需要证明,对 μ>0 ,存在一个 K>0 的常数,使得对于每个 x B μ ={ x PC μ:xPC( J,X ) } ,使得 Tx PC K 。现在,设 { x n } 是有界子集 M B μ 上的一个序列,对于每个 x n M ,通过假设(H1)~(H3),有

( T x n )( t ) | b α |+ g( x n ) + 0< t i <t I i ( x n ( t i ) ) + c α f( t, x n ( t ) ) + d α 0< t i <t 0 t i f( s, x n ( s ) ) ds + d α 0 t f( s, x n ( s ) ) ds | b α |+ C 1 x α 3 + C 2 +n( γ 1 x α 2 + γ 2 ) + c α ( δ 1 x α 1 + δ 2 )+ d α nT( δ 1 x α 1 + δ 2 )+ d α T( δ 1 x α 1 + δ 2 ) | b α |+ C 1 μ α 3 + C 2 +n( γ 1 μ α 2 + γ 2 ) +( c α + d α nT+ d α T ) δ 1 μ α 1 +( c α + d α nT+ d α T ) δ 2 :=K.

因此, ( Tx )( t ) K ( T x n ) M上一致有界。

引理5算子T在 PC( J,X ) 上是等度连续的。

证明 设 { x n } 是在引理4中定义的有界子集M的一个序列,对于 t 1 , t 2 J, t 1 < t 2 ,有

( T x n )( t 2 )( T x n )( t 1 ) 0< t i < t 2 t 1 I i ( x n ( t i ) ) + c α [ f( t 2 , x n ( t 2 ) )f( t 1 , x n ( t 1 ) ) ] + d α 0< t i < t 2 t 1 0 t i f( s, x n ( s ) )ds + d α t 1 t 2 f( s,x( s ) )ds 0< t i < t 2 t 1 I i ( x n ( t i ) ) + c α f( t 2 , x n ( t 2 ) )f( t 1 , x n ( t 1 ) ) + d α 0< t i < t 2 t 1 0 t i f( s, x n ( s ) ) ds + d α t 1 t 2 f( s,x( s ) ) ds 0< t i < t 2 t 1 I i ( x n ( t i ) ) + c α f( t 2 , x n ( t 2 ) )f( t 1 , x n ( t 1 ) ) + d α 0< t i < t 2 t 1 0 t i f( s, x n ( s ) ) ds + d α ( δ 1 x α 1 + δ 2 ) ( t 2 t 1 ).

t 2 t 1 时,上述不等式右端趋于0。因此, ( T x n ) 是等度连续的。

引理6算子T PC( J,X ) 上是紧的。

证明 设一个闭子集 KPC( J,X ) ,由引理4,5知, T:PC( J,X )PC( J,X ) 是有界且等度连续的。根据Arzela-Ascoli定理,可以得出T(K)是 PC( J,X ) 的一个相对紧的子集。因此, T:PC( J,X )PC( J,X ) 是紧算子。

定理4假设(H1)-(H3)成立,脉冲分数阶微分方程(1)至少有一个解。

证明 易知,算子T的不动点是(1)的解。显然,算子 T:PC( J,X )PC( J,X ) 是连续且完全连续的,则证明 S( T )={ xPC( J,X ):x=λTx,λ[ 0,1 ] } 是有界的,设 xS( T ) ,则对某个 λ[ 0,1 ] x=λTx ,有

x( t ) PC =λ ( Tx )( t ) | b α |+ g( x ) + 0< t i <t I i ( x( t i ) ) + c α f( t,x( t ) ) + d α 0< t i <t 0 t i f( s,x( s ) ) ds + d α 0 t f( s,x( s ) ) ds          | b α |+ C 1 x α 3 + C 2 +n( γ 1 x α 2 + γ 2 ) + c α ( δ 1 x α 1 + δ 2 )+ d α nT( δ 1 x α 1 + δ 2 )+ d α T( δ 1 x α 1 + δ 2 )          | b α |+ C 1 μ α 3 + C 2 +n( γ 1 μ α 2 + γ 2 ) +( c α + d α nT+ d α T ) δ 1 μ α 1 +( c α + d α nT+ d α T ) δ 2          K.

其中 α 1 , α 2 , α 3 [ 0,1 ) ,由引理4的结果表明,S PC( J,X ) 中有界。根据Schaefer不动点定理,T至少有一个不动点,该不动点是(1)的解。

定理5假设(H1)~(H3)成立,如果

C g +n γ I + c α δ f + d α nT δ f + d α T δ f <1

那么分数阶脉冲问题(1)在 PC( J,X ) 上有唯一解。

证明 T PC( J,X ) 上的压缩映射。对 x,yPC( J,X ) ,有

( Tx )( t )( Ty )( t ) g( x )g( y ) + 0< t i <t I i ( x( t i ) ) I i ( y( t i ) ) + c α f( t,x( t ) )f( t,y( t ) ) + d α 0< t i <t 0 t i f( s,x( s ) )f( s,y( s ) ) ds + d α 0 t f( s,x( s ) )f( s,y( s ) ) ds C g xy + 0< t i <t γ I xy + c α δ f xy + d α 0< t i <t 0 t i δ f xy ds + d α 0 t δ f xy ds ( C g +n γ I + c α δ f + d α nT δ f + d α T δ f ) xy

因此,T PC( J,X ) 上的压缩映射。由压缩映射原理,得出算子T PC( J,X ) 上有唯一不动点,即问题(1)有唯一解。

4. 例子

考虑下面的Caputo-Fabrizio型脉冲分数阶微分方程:

{ D CF 1 5 x( t )= x( t ) ( 1+ e 2t )( 1+ x( t ) ) ,tJ=[ 0,1 ],t 1 3 , Δx( 1 3 )= 1 8 x( 1 3 ) , x( 0 )+ 1 9 x( t )=0. (3)

解 由(3)和定义2可知, α= 1 5 n=1 T=1 M( 1 5 )= 10 9 c 1 5 = 4 5 d 1 5 = 1 5 f( t,x )= x ( 1+ e 2t )( 1+x ) I 1 ( x )= 1 8 x

t[ 0,1 ] ,有

f( t,x )f( t,y ) = 1 1+ e 2t x 1+x y 1+y 1 2 xy ( 1+x )( 1+y ) 1 2 xy ,

因此, δ f = 1 2 δ 1 = 1 2 α 1 =1 。同样,有

I( x )I( y ) = 1 8 xy ,

因此, γ I = 1 8 α 2 =1

g( x )g( y ) = 1 9 xy ,

因此, C g = 1 9 α 3 =1 。所以,定理4的条件成立,

C g +n γ I + c α δ f + d α nT δ f + d α T δ f 0.836<1.

因此,根据定理4,给定分数阶脉冲问题(3)在 [ 0,1 ] 上有唯一解。

5. 小结

本文主要讨论了Banach空间中含有脉冲点的Caputo-Fabrizio型分数阶微分方程解的存在性问题,在方程同时具有脉冲项及非局部条件的情况下,利用Schaefer不动点定理,压缩映射原理和Arzela-Ascoli定理,讨论了 0<α<1 阶的分数阶脉冲微分方程,得到了解的存在性定理,最后用一个例子验证结论的正确性。这一结果不仅丰富和发展了分数阶微分方程的理论体系,而且为解决实际问题提供了理论依据。同时,这也为今后进一步探讨此类方程的解的性质、稳定性和数值计算方法奠定了基础。然而,本文的研究仍具有一定的局限性,未来研究可以在此基础上,考虑更广泛的脉冲条件和分数阶导数,以拓展方程的应用范围。

NOTES

*第一作者。

#通讯作者。

参考文献

[1] Baleanu, D., Jajarmi, A., Mohammadi, H. and Rezapour, S. (2020) A New Study on the Mathematical Modelling of Human Liver with Caputo-Fabrizio Fractional Derivative. Chaos, Solitons & Fractals, 134, Article 109705.
https://doi.org/10.1016/j.chaos.2020.109705
[2] Atangana, A. and Alqahtani, R.T. (2016) Numerical Approximation of the Space-Time Caputo-Fabrizio Fractional Derivative and Application to Groundwater Pollution Equation. Advances in Difference Equations, 2016, Article No. 156.
https://doi.org/10.1186/s13662-016-0871-x
[3] Caputo, M. and Fabrizio, M. (2015) A New Definition of Fractional Derivative without Singular Kernel. Progress in Fractional Differentiation & Applications, 1, 73-85.
[4] Losada, J. and Nieto, J.J. (2015) Properties of a New Fractional Derivative without Singular Kernel. Progress in Fractional Differentiation & Applications, 1, 87-92.
[5] Abbas, M.I. (2020) On the Initial Value Problems for the Caputo-Fabrizio Impulsive Fractional Differential Equations. Asian-European Journal of Mathematics, 14, Article 2150073.
https://doi.org/10.1142/s179355712150073x
[6] Kaliraj, K., Manjula, M. and Ravichandran, C. (2022) New Existence Results on Nonlocal Neutral Fractional Differential Equation in Concepts of Caputo Derivative with Impulsive Conditions. Chaos, Solitons & Fractals, 161, Article 112284.
https://doi.org/10.1016/j.chaos.2022.112284
[7] Jeet, K. and Sukavanam, N. (2020) Approximate Controllability of Nonlocal and Impulsive Neutral Integro-Differential Equations Using the Resolvent Operator Theory and an Approximating Technique. Applied Mathematics and Computation, 364, Article 124690.
https://doi.org/10.1016/j.amc.2019.124690
[8] Wang, J., Zhou, Y. and Wei, W. (2012) Study in Fractional Differential Equations by Means of Topological Degree Methods. Numerical Functional Analysis and Optimization, 33, 216-238.
https://doi.org/10.1080/01630563.2011.631069
[9] 许天周. 应用泛函分析[M]. 北京: 科学出版社, 2002.
[10] Zhou, Y. (2022). Basic Theory of Fractional Differential Equations. 3rd Edition, World Scientific.
https://doi.org/10.1142/13289
[11] 程民德, 邓东皋, 龙瑞麟. 实分析[M]. 北京: 高等教育出版社, 2008.