Loading [MathJax]/jax/output/HTML-CSS/jax.js
  数学物理学报  2015, Vol. 35 Issue (6): 1059-1070   PDF (360 KB)    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
冯育强
王蔚敏
李寿贵
带有奇异非线性项的分数微分方程周期解的存在性与唯一性
冯育强1, 王蔚敏2, 李寿贵2    
1 武汉科技大学理学院 武汉 430065;
2 冶金工业过程系统科学湖北省重点实验室 武汉 430081
摘要: 该文目的在于给出如下分数阶微分方程解的存在唯一性结论
{Dαx(t)=f(t,x(t)),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1), (PBVP)
其中ft=0可以是奇异的.主要的工具是上下解方法、最大值原理和单调迭代技术.最后举例说明所获结论的应用.
关键词: 分数微分方程     周期边值问题     存在性     唯一性     奇异性    
Existence and Uniqueness Results for the Periodic Boundary Value Problems of Fractional Differential Equations with Singular Nonlinearities
Feng Yuqiang1, Wang Weimin2, Li Shougui2    
1 School of Science, Wuhan University of Science and Technology, Wuhan 430065;
2 Hubei Province Key Laboratory of Systems Science in Metallurgical Process, Wuhan 430081
Abstract: The purpose of this paper is to give some sufficient conditions for the existence and uniqueness of solutions to the fractional differential equation as follows
{Dαx(t)=f(t,x(t)),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1), (PBVP)
where Dα denotes the Riemann-Liouville fractional derivative, f may be singular at t=0. Lower and upper solutions method, maximum principle together with iterative technique are employed. An example is presented to illustrate the application of results obtained.
Key words: Fractional differential equation     Periodic boundary value problem     Existence     Uniqueness     Singularity    
1 引 言

由于分数微积分具有刻画记忆和遗传性质的特性,所以在描述具有记忆和遗传性质的各种物质、 流变学、材料和力学系统、信号处理和系统辨识、ANN (神经网络)、分形和混沌时, 分数微积分比整数阶微积分更加准确. 分数阶模型涵盖了流体力学,流变学,粘弹性力学, 分数控制系统与分数控制器,各种电子回路,电分析化学,生物系统的电传导, 神经的分数模型以及分数回归模型等领域(见文献[1]及其参考文献). 正是由于其广泛而重要的应用,分数微分方程成为了目前理论研究的热点, 涌现出大量的研究成果,如文献[2,3,4,5,6,7]及其参考文献.

最近,分数阶微分方程的周期边值问题引起了诸多关注. 文献[8]首次考察了如 下分数阶微分方程边值问题的可解性

{Dαx(t)λx(t)=f(t,x(t)),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1). 作者计算了该问题的Green函数,并利用Schaeffer不动点定理和Bananch不动点定理证 明了一些解的存在性和唯一性结论. 随后,在文献[9]和[11]中, 作者给出了以下周期边值问题的一些比较原理

{Dαx(t)λx(t)=σ(t),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1). 文献[10]利用单调迭代方法研究了分数微分方程周期解的存在性与唯一性. 文献[12]则将该项研究推广到时滞微分方程的情形. 在这些文献中,大都要求非线 性项是Lipschitz 连续的或者单增且连续.

受以上工作以及减算子不动点研究[13]的启发,本文探讨了当非线性项为奇异的且满足 一定的递减性质时,如下的分数阶微分方程接的存在性与唯一性

{Dαx(t)=f(t,x(t)),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1). (PBVP)
主要工具是上下解方法、极大值原理和单调迭代技术.

本文的结构为: 第2节给出一些必要的记号与预备知识; 解的存在性与唯一性结论在第3节给出; 第4节举例说明所获结论的应用.

2 预备知识

设C[0, 1]表示[0,1]所有连续函数的集合. 若在该线性空间赋予范数

u=max0t1|u(t)|, 则(C[0,1],||.||)是一个Banach空间. 对于 r>0,u C(0,1],定义

tru(t)|t=0=limt0+tru(t), Cr[0,1]={uC(0,1]|tru(t)C[0,1]}.Cr[0,1]在范数 ur=max0t1|tru(t)| 下也成为一个Banach空间.

0<r<1,Cr[0,1]C(0,1]L1[0,1],这里 L1[0,1] 表示[0,1] 上所有Lebesgue可积的函数集合.

P={uCr[0,1]u(t)0,t(0,1]}.

2.1 P 是一个正规锥.

(1) PCr[0,1]中的一个锥.

这可以用P的定义来验证.

(2) P 是正规锥.

事实上,设 u,vP 满足 0u(t)v(t),t(0,1],则有 0tru(t)trv(t),t[0,1],即 urvr.

因此,P 是一个正规锥.

一个连续函数u:(0,1]Rα(0<α<1) 阶 Riemann-Liouville导数定义为

Iαu(t)=1Γ(α)t0(tτ)α1u(τ)dτ, 其中 Γ(.) 表示Γ函数.

引理2.20<α<1,如果 uCα[0,1],IαuC[0,1].

uCα[0,1],则

Iαu(t)=1Γ(α)t0(tτ)α1u(τ)dτ=1Γ(α)t0(tτ)α1(τ)αU(τ)dτ  (U(τ)=ταu(τ))=1Γ(α)10(1s)α1(s)αU(st)ds  (s=τt).

注意到 U(t)[0,1] 上一致连续,即对 ϵ>0, 存在 δ>0,当 t1,t2[0,1]|t1t2|<δ就有

|U(t1)U(t2)|<ϵ. 所以

|Iαu(t1)Iαu(t2)|=|1Γ(α)10(1s)α1(s)αU(st1)ds1Γ(α)10(1s)α1(s)αU(st2)ds|=|1Γ(α)10(1s)α1(s)α[U(st1)U(st2)]ds|1Γ(α)10(1s)α1(s)αU(st1)U(st2)dsϵ1Γ(α)10(1s)α1(s)αds  ( st1st2∣≤|t1t2|<δ)=ϵB(α,1α)Γ(α)=Γ(1α)ϵ, 其中 B(.,.)表示Beta函数.

因此有 IαuC[0,1].

一个连续函数uC(0,1]α(0<α<1)阶Riemann-Liouville导数定义为

Dαu(t)=1Γ(1α)ddtt0(tτ)αu(τ)dτ=ddtI1αu(t). 作为该定义的直接结论,分数阶导数具有以下性质

(1)当 β>1时,

Dα(tβ)=1Γ(1α)ddtt0(tτ)ατβdτ=Γ(β+1)Γ(β+2α)ddt(t1+βα).

(2)对于任意的 uC(0,1]L1[0,1],DαIαu=u.

(3)当 uC(0,1]L1[0,1]时,分数微分方程 Dαu(t)=0 有解为 u(t)=ctα1,cR.

(4)当 uC(0,1]L1(0,1]DαuC(0,1]L1(0,1](0<α<1)时,下式成立

IαDαu(t)=u(t)+ctα1, cR.

对于 α>0,β>0,二元的 Mittag-Leffler 函数定义如下

Eα,β(z)=k=0zkΓ(αk+β),  zC. 据文献[10]可知,当 x,yR,  x<0<y

0<Eα,α(x)<Eα,α(0)=1Γ(α)<Eα,α(y), limx+Eα,α(x)=+,  limxEα,α(x)=0.

为证明主要结论,需要以下的引理.

2.3[9]hC[0,1],α(0,1),λREα,α(λ)<1Γ(α). 如果 u(t)C1α[0,1] 满足

(i) Dαu(t)λu(t)=h(t)(0<t<1);

(ii) limt0+t1αu(t)=u(1).

则有

u(t)=10Gλ,α(t,s)h(s)ds, 其中 Gλ,α(t,s) 的定义为: 当 0st1,

Γ(α)Eα,α(λtα)Eα,α(λ(1s)α)tα1(1s)α11Γ(α)Eα,α(λ)+(ts)α1Eα,α(λ(ts)α).

0ts1,

Γ(α)Eα,α(λtα)Eα,α(λ(1s)α)tα1(1s)α11Γ(α)Eα,α(λ).

2.1 据文献[10,引理2.2],当 α(0,1), λR, λ<0时, 0<Eα,α(λ)<1Γ(α). 如果Eα,α(λ)<1Γ(α), 容易验证 Gλ,α(t,s)>0 for (t,s)(0,1)×(0,1).

2.4[9]α(0,1),λREα,α(λ)<1Γ(α). 如果 u(t)C1α[0,1] 满足

(i) Dαu(t)λu(t)0(0<t<1);

(ii) limt0+t1αu(t)=u(1).

u(t)0,t(0,1].

3 主要结论

在本节中,始终假定以下条件成立.

(H1) f:(0,1]×RR 是连续函数;

(H2) 对于xC1α[0,1], f(t,x(t))C1α[0,1], 其中

t1αf(t,x(t))|t=0:=limt0+t1αf(t,x(t)). 值得注意的是: 当 (H1)(H2) 成立时,f(t,x) 可以在t=0 是奇异的.

易见 (PBVP) 与如下问题等价

{Dαx(t)λx(t)=f(t,x(t))λx(t),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1), 其中参数 λR 且满足 Eα,α(λ)<1Γ(α).

对于 xC1α[0,1],定义映射 T 如下

(Tx)(t)=10Gλ, α(t,s)[f(s,x(s))λx(s)]ds, 其中 Gλ, α(t,s) 的定义见引理2.3. 于是 xC1α(PBVP) 的一个解当且仅当 xC1αT的一个不动点.

3.1 当条件 (H1)(H2)成立时, 算子 TC1α[0,1]C1α[0,1]T 是连续的.

结论分两步证明.

(1) 若 xC1α[0,1],则TxC1α[0,1].

显然当 xC1α[0,1] 时,[f(t,x(t))λx(t)]C1α[0,1].

σ(t)=f(t,x(t))λx(t),u(t)=(Tx)(t),则有

{Dαu(t)λu(t)=σ(t),tJ:=(0,1], 0<α<1,limt0+t1αu(t)=u(1).

据文献[9],当 σC1α[0,1]Eα,α(λ)<1Γ(α) 时,下述线性分数微分方程边值问题

{Dαu(t)λu(t)=σ(t),tJ:=(0,1], 0<α<1,limt0+t1αx(t)=x(1) 的解是

u(t)=u0Γ(α)t1αEα,α(λtα)+t0(ts)α1Eα,α(λ(ts)α)σ(s)ds, 其中 u0为一常数.

以下来证明 uC1α[0,1].

事实上,

t1αu(t)=u0Γ(α)t2(1α)Eα,α(λtα)+t1αt0(ts)α1Eα,α(λ(ts)α)σ(s)ds

只需说明 t1αt0(ts)α1Eα,α(λ(ts)α)σ(s)dst的连续函数. 注意到

t1αt0(ts)α1Eα,α(λ(ts)α)σ(s)ds=k=0t1αλkαΓ(αk+α)t0(ts)(k+1)α1σ(s)ds=k=0t1αλkαΓ(αk+α)t0(ts)(k+1)α1sα1ϕ(s)ds=k=0t(k+1)αλkαΓ(αk+α)10(1u)(k+1)α1uα1ϕ(ut)du, 其中ϕ(t)=t1ασ(t)C[0,1]. 由于 ϕ(t)[0,1] 上一致连续,即 ϵ>0,存在δ>0,当 t1,t2[0,1]|t1t2|<δ时,

|ϕ(t1)ϕ(t2)|<ϵ. 此时,

|10(1u)(k+1)α1uα1ϕ(ut1)du10(1u)(k+1)α1uα1ϕ(ut2)du|=|10(1u)(k+1)α1uα1(ϕ(ut1)ϕ(ut2))du|10(1u)(k+1)α1uα1|ϕ(ut1)ϕ(ut2)|duϵ10(1u)(k+1)α1uα1du=ϵB((k+1)α,α). 这说明对于任意的kt(k+1)αΓ(αk+α)10(1u)(k+1)α1uα1ϕ(ut)dut 的连续函数.

因此,t1αu(t)C[0,1].

(2) 在假设 (H1)(H2)下,利用Lebesgue控制收敛定理可以证明 T 是连续的.

u,vC1α[0,1] 满足 u(t)v(t),t[0,1],记

[u,v]={xC1α[0,1]|u(t)x(t)v(t),t(0,1]}.

3.1 如果存在 u,vC1α[0,1],uv 和常数 λ (λ 满足Eα,α(λ)<1Γ(α)), 使得以下条件成立

(A1) 对于任意的t[0,1],x,y[u,v],xy,

f(t,y(t))f(t,x(t))λ(y(t)x(t));

(A2) 对于任意的t[0,1],x,y[u,v],

f(t,lx(t)+(1l)y(t)lf(t,x(t))+(1l)f(t,y(t));

(A3) 对于任意的 t(0,1),

Dα(v+u)(t)+λ(vu)(t)2f(t,v(t)), f(t,u(t))Dαv(t)λ(vu)(t);

(A4) limt0+t1αu(t)=u(1), limt0+t1αv(t)=v(1).

那么(PBVP) 存在唯一解 xC1α[0,1][u,v].

分数阶微分方程周期边值问题 (PBVP) 等价于如下边值问题

{Dαx(t)λx(t)=f(t,x(t))λx(t),0<t<1,limt0+t1αx(t)=x(1).

x0 为下述边值问题的解

{Dαx(t)λx(t)=f(t,v(t))λv(t),limt0+t1αx(t)=x(1),

x0=10Gλ,α(t,s)[f(s,v(s))λv(s)]ds, 这里 Gλ,α(t,s) 的定义见引理 2.2. 构造迭代列 {xn} 如下

{Dαxn+1(t)λxn+1(t)=f(t,xn(t))λxn(t),limt0+t1αxn+1(t)=xn+1(1),n=0,1,2,.

由于 f 连续,所以点列 {xn} 有意义且有 {xn}C1α[0,1].

以下来证明 {xn} 收敛到问题 (PBVP) 在序区间 [u,v] 中的唯一解.

证明过程分为5步.

第1步) v(t)x0(t)(u+v2)(t)u(t).

由定理条件 (A1),(A3),(A4),可知

Dα(u+v2)(t)λ(u+v2)(t)f(t,v(t))λv(t)f(t,u(t))λu(t)Dαv(t)λv(t). 于是,根据 x0 的定义易得

Dαv(t)λv(t)Dαx0(t)λx0(t)Dα(u+v2)(t)λ(u+v2)(t). 注意到

limt0+t1αu(t)=u(1),limt0+t1αv(t)=v(1),limt0+t1αx0(t)=x0(1). 据引理 2.5,有

v(t)x0(t)(u+v2)(t)u(t).

第2步)对于n=0,1,2,, 成立着

v(t)x2n+1(t)x2n(t)(u+v2)(t)u(t).

利用数学归纳法来证明结论.

事实上,由于

Dαx1(t)λx1(t)=f(t,x0(t))λx0(t).

根据第1步证明及条件(A1)

f(t,u(t))λu(t)f(t,x0(t))λx0(t)f(t,v(t))λv(t). 注意到

f(t,u(t))λu(t)Dαv(t)λv(t), 所以有

Dαv(t)λv(t)Dαx1(t)λx1(t)Dαx0(t)λx0(t).

根据引理 2.4,

v(t)x1(t)x0(t). 所以当 n=0结论成立.

设若 n=k 时结论也成立,即

v(t)x2k+1(t)x2k(t)(u+v2)(t)u(t).

则由假设 (A1),对于任意的 t[0,1]

f(t,u(t))λu(t)f(t,x2k(t))λx2k(t)f(t,x2k+1(t))λx2k+1(t)f(t,v(t))λv(t). 这说明

Dαv(t)λv(t)Dαx2k+1(t)λx2k+1(t)Dαx2k+2(t)λx2k+2(t)Dαx0(t)λx0(t). 再根据引理 2.3,可知

v(t)x2k+1(t)x2k+2(t)x0(t)u(t).

重复此过程,可得

{f(t,u(t))λu(t)f(t,x2k+2(t))λx2k+2(t)f(t,x2k+1(t))λx2k+1(t)f(t,v(t))λv(t),Dαv(t)λv(t)Dαx2k+3(t)λx2k+3(t)Dαx2k+2(t)λx2k+2(t)Dαx0(t)λx0(t),v(t)x2k+3(t)x2k+2(t)x0(t)u(t), 即当 n=k+1 结论也成立.

所以对所有的 n{0}N,结论都成立.

第3步) {x2n(t)} 递增,{x2n+1(t)} 递减.

由第2步证明过程可知

v(t)x1(t)u(t). 于是根据条件 (A1)

f(t,u(t))λu(t)f(t,x1(t))λx1(t)f(t,v(t))λv(t). 这表明

Dαv(t)λv(t)Dαx2(t)λx2(t)Dαx0(t)λx0(t). v(t)x2(t)x0(t).

类似于步骤2,利用数学归纳法可以证明 {x2n(t)} 是递增的. 同理可证{x2n+1(t)} 是递减的.

第4步) {xn(t)} 收敛到问题 (PBVP)的一个解.

由步骤1--3,对于 n=0,1,2,,

u(t)(u+v2)(t)x2n(t)x2n+2(t)x2n+3(t)x2n+1(t)v(t).

yn(t)=xn(t)u(t),则对 n=0,1,2,,

0v(t)u(t)2y2n(t)y2n+2(t)y2n+3(t)y2n+1(t)v(t)u(t).

rn=sup{rR|y2n(t)ry2n+1(t)}, 则数列 {rn}是良定的. 易验证 12rn1{rn} 递增,所以 {rn} 收敛.

由假设 (A1),(A2),可得

y2n+3(t)y2n+1(t)=x2n+1(t)u(t)=10Gλ,α(t,s)[f(s,x2n(s))λx2n(s)]dsu(t)=10Gλ,α(t,s)[f(s,(y2n+u)(s))λ(y2n+u)(s)]dsu(t)10Gλ,α(t,s)[f(s,(rny2n+1+u)(s))λ(rny2n+1+u)(s)]dsu(t)=10Gλ,α(t,s)[f(s,rnx2n+1+(1rn)u(s))λ(rnx2n+1(s)+(1rn)u(s))]dsu(t)10Gλ,α(t,s)[rnf(s,x2n+1(s))+(1rn)f(s,u(s))λ(rnx2n+1(s)+(1rn)u(s))]dsu(t)=rn10Gλ,α(t,s)[f(s,x2n+1(s))λx2n+1(s)]ds+(1rn)10Gλ,α(t,s)[f(s,u(s))λu(s)]dsu(t)rnx2n+2(t)+(1rn)v(t)u(t)=rny2n+2(t)+(1rn)(v(t)u(t))=rny2n+2(t)+2(1rn)(v(t)u(t))2rny2n+2(t)+2(1rn)y2n+2(t)(2rn)y2n+2(t), 从而

rn+1=sup{rR|y2n+3(t)ry2n+2(t)}12rn.

limnrn=r,上述不等式表明 r12r,所以 r=1.

对任意的偶数 p>0,

0y2n+py2ny2n+1y2n(1rn)y2n+1(1rn)(vu)(t).rn1 及锥P 的正规性可知,{y2n}C1α[0,1] 中的Cauchy列. 类似地,可以证明 {y2n+1} 也是C1α[0,1] 中的Cauchy列. 于是

limny2n=limny2n+1, 所以 {yn} 收敛.

y=limnyn,x=y(t)+u(t), 则有

limnxn(t)=x(t). 以及

x(t)=limn10Gλ,α(t,s)[f(s,xn(s))λxn(s)]ds=10Gλ,α(t,s)[f(s,x(s))λx(s)]ds, 这表明 x 是问题 (PBVP) 的一个解.

第5步) x 是问题 (PBVP) 在序区间 [u,v]中的唯一解.

¯x[u,v](PBVP) 的一个解,则

f(t,v(t))λv(t)f(t,¯x(t))λ¯x(t)f(t,u(t))λu(t). 从而

Dαx0(t)λx0(t)Dα¯x(t)λ¯x(t)Dαv(t)λv(t).

据引理 2.4,有 x0(t)¯x(t)v(t).

因而

f(t,v(t))λv(t)f(t,¯x(t))λ¯x(t)f(t,x0(t))λx0(t). 这表明

Dαx0(t)λx0(t)Dα¯x(t)λ¯x(t)Dαx1(t)λx1(t). 由引理2.4可知 x0(t)¯x(t)x1(t).

继续这一过程,容易验证

x2n(t)¯x(t)x2n+1(t),n=0,1,2,. 因此

¯x(t)=limnxn(t)=x(t). 证毕.

3.1 定理3.1是一个局部性结论, 即在一个局部其解是存在且唯一的; 但从全局来看,边值问题可能会有多个解.

3.1 如果存在常数 c>0,λ<0 满足

1) λc2t1αf(t,ct1α)t1αf(t,0)λc;

2) 对于 t[0,1],f(t,.)为凸函数;

3) 对于 t[0,1],x,yR,xy,成立着f(t,y)f(t,x)λ(yx);

4) Eα,α(λ)<1Γ(α).

(PBVP) 有唯一解 u 满足 0u(t)ctα1,t[0,1]. 进一步, 如果存在 t0(0,1) 使得 f(t0,0)0,那么 u(t)>0,t(0,1).

在推论3.1的假设下,设 u(t)0,v(t)=ctα1,可以验证定理3.1的条件(A1)--(A5)都满足. 于是根据定理3.1可知周期边值问题存在唯一解 u 满足 0u(t)ctα1,t[0,1].

进一步,如果存在 t0(0,1] 使得 f(t0,0)0,那么 u(t)(0,1)内不会恒为0,从而根据引理2.3得 u(t)>0, t(0,1).

4 例子

考察如下的边值问题

{Dαx(t)=ex(t)+λx(t)+btα1,0<t<1,limt0+t1αx(t)=x(1), 其中 f(t,x)=ex+λx+btα1,显然ft=0 是奇异的. 如果 Eα,α(λ)<1Γ(α), λ2b(λ+1),那么据推论3.1可知该问题有唯一解 u满足

0u(t)ctα1,  t(0,1]. 进一步,由f(t,0)0 可知 u(t)>0,t(0,1).

参考文献
[1] Podlubny I. Fractional Differential Equations. San Diego:Academic Press, 1999
[2] 胡雷, 张淑琴, 侍爱玲.分数阶微分方程耦合系统共振边值问题解的存在性. 数学物理学报, 2014, 34(5):1313-1326
[3] Chen P, Li Y, Li Q. Existence of mild solutions for fractional evolution equations with nonlocal initial conditions. Ann Polon Math, 2014, 110(1):13-24
[4] Lakshmikantham V. Theory of fractional functional differential equations. Nonlinear Anal, 2008, 69(10):3337-3343
[5] Bai Z, Qiu T. Existence of positive solution for singular fractional differential equation. Appl Math Comput, 2009, 215(1):2761-2767
[6] Zhang S. Monotone iterative nethod for initial value problem involving Riemann-Liouville fractional derivatives. Nonlinear Anal, 2009, 71(5/6):2087-2093
[7] Benchohra M, Hellar M. Global uniqueness results for fractional partial hyperbolic differential equations with state-dependent delay. Ann Polon Math, 2014, 110(3):259-281
[8] Belmekki M, Nieto Juan J, Rodriguez-Lopez R. Existence of periodic solution for a nonlinear fractional differential equations. Boundary Value Problem, 2009, Article ID 324561, 18 pages
[9] Nieto Juan J. Maximum principles for fractional differential equations derived from Mittag-Leffler functions. Appl Math Lett, 2010, 23(10):1248-1251
[10] Wei Z, Dong W, Che J. Periodic boundary value problems for fractional differential equations involving a Riemann-Liouville fractional derivative. Nonlinear Anal, 2010, 73(10):3232-3238
[11] Nieto Juan J. Comparision results for periodic boundary value problem of fractional differential equations. Fractional Differential Calculus, 2011, 1(1):99-104
[12] Karami H, Babakhani A, Baleanu D. Existence results for a class of fractional differential equations with periodic boundary conditions and with delay. Abstract and Appllied Analysis, 2013, Article ID 176180, 8 pages
[13] 李福义, 冯锦峰, 沈沛龙. 一类减算子的不动点定理及其应用. 数学学报, 1999, 42(2):193-196
带有奇异非线性项的分数微分方程周期解的存在性与唯一性
冯育强, 王蔚敏, 李寿贵