由于分数微积分具有刻画记忆和遗传性质的特性,所以在描述具有记忆和遗传性质的各种物质、 流变学、材料和力学系统、信号处理和系统辨识、ANN (神经网络)、分形和混沌时, 分数微积分比整数阶微积分更加准确. 分数阶模型涵盖了流体力学,流变学,粘弹性力学, 分数控制系统与分数控制器,各种电子回路,电分析化学,生物系统的电传导, 神经的分数模型以及分数回归模型等领域(见文献[1]及其参考文献). 正是由于其广泛而重要的应用,分数微分方程成为了目前理论研究的热点, 涌现出大量的研究成果,如文献[2,3,4,5,6,7]及其参考文献.
最近,分数阶微分方程的周期边值问题引起了诸多关注. 文献[8]首次考察了如 下分数阶微分方程边值问题的可解性
受以上工作以及减算子不动点研究[13]的启发,本文探讨了当非线性项为奇异的且满足 一定的递减性质时,如下的分数阶微分方程接的存在性与唯一性
本文的结构为: 第2节给出一些必要的记号与预备知识; 解的存在性与唯一性结论在第3节给出; 第4节举例说明所获结论的应用.
设C[0, 1]表示[0,1]所有连续函数的集合. 若在该线性空间赋予范数
当 0<r<1,Cr[0,1]⊂C(0,1]⋂L1[0,1],这里 L1[0,1] 表示[0,1] 上所有Lebesgue可积的函数集合.
设 P={u∈Cr[0,1]∣u(t)≥0,∀t∈(0,1]}.
引理2.1 P 是一个正规锥.
证 (1) P 是Cr[0,1]中的一个锥.
这可以用P的定义来验证.
(2) P 是正规锥.
事实上,设 u,v∈P 满足 0≤u(t)≤v(t),∀t∈(0,1],则有 0≤tru(t)≤trv(t),∀t∈[0,1],即 ‖u‖r≤‖v‖r.
因此,P 是一个正规锥.
一个连续函数u:(0,1]→R的 α(0<α<1) 阶 Riemann-Liouville导数定义为
引理2.2 设 0<α<1,如果 u∈Cα[0,1], 则 Iαu∈C[0,1].
证 设 u∈Cα[0,1],则
注意到 U(t) 在 [0,1] 上一致连续,即对 ∀ϵ>0, 存在 δ>0,当 t1,t2∈[0,1] 且|t1−t2|<δ就有
因此有 Iαu∈C[0,1].
一个连续函数u∈C(0,1]的α(0<α<1)阶Riemann-Liouville导数定义为
(1)当 β>−1时,
Dα(tβ)=1Γ(1−α)ddt∫t0(t−τ)−ατβdτ=Γ(β+1)Γ(β+2−α)ddt(t1+β−α).
(2)对于任意的 u∈C(0,1]⋂L1[0,1],DαIαu=u.
(3)当 u∈C(0,1]⋂L1[0,1]时,分数微分方程 Dαu(t)=0 有解为 u(t)=ctα−1,c∈R.
(4)当 u∈C(0,1]⋂L1(0,1] 且 Dαu∈C(0,1]⋂L1(0,1](0<α<1)时,下式成立
对于 α>0,β>0,二元的 Mittag-Leffler 函数定义如下
为证明主要结论,需要以下的引理.
引理2.3[9] 设 h∈C[0,1],α∈(0,1),λ∈R 且 Eα,α(λ)<1Γ(α). 如果 u(t)∈C1−α[0,1] 满足
(i) Dαu(t)−λu(t)=h(t)(0<t<1);
(ii) limt→0+t1−αu(t)=u(1).
当 0≤t≤s≤1,
注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),λ∈R 且 Eα,α(λ)<1Γ(α). 如果 u(t)∈C1−α[0,1] 满足
(i) Dαu(t)−λu(t)≥0(0<t<1);
在本节中,始终假定以下条件成立.
(H1) f:(0,1]×R⟶R 是连续函数;
(H2) 对于∀x∈C1−α[0,1], f(t,x(t))∈C1−α[0,1], 其中
易见 (PBVP) 与如下问题等价
对于 x∈C1−α[0,1],定义映射 T 如下
引理3.1 当条件 (H1)−(H2)成立时, 算子 T 映 C1−α[0,1] 到 C1−α[0,1] 且 T 是连续的.
证 结论分两步证明.
(1) 若 x∈C1−α[0,1],则Tx∈C1−α[0,1].
显然当 x∈C1−α[0,1] 时,[f(t,x(t))−λx(t)]∈C1−α[0,1].
设 σ(t)=f(t,x(t))−λx(t),u(t)=(Tx)(t),则有
据文献[9],当 σ∈C1−α[0,1] 且 Eα,α(λ)<1Γ(α) 时,下述线性分数微分方程边值问题
以下来证明 u∈C1−α[0,1].
事实上,
只需说明 t1−α∫t0(t−s)α−1Eα,α(λ(t−s)α)σ(s)ds 是 t的连续函数. 注意到
因此,t1−αu(t)∈C[0,1].
(2) 在假设 (H1)−(H2)下,利用Lebesgue控制收敛定理可以证明 T 是连续的.
设 u,v∈C1−α[0,1] 满足 u(t)≤v(t),t∈[0,1],记
定理3.1 如果存在 u,v∈C1−α[0,1],u≤v 和常数 λ (λ 满足Eα,α(λ)<1Γ(α)), 使得以下条件成立
(A1) 对于任意的t∈[0,1],x,y∈[u,v],x≤y,
(A2) 对于任意的t∈[0,1],x,y∈[u,v],
(A3) 对于任意的 t∈(0,1),
(A4) limt→0+t1−αu(t)=u(1), limt→0+t1−αv(t)=v(1).
那么(PBVP) 存在唯一解 x∗∈C1−α[0,1]∩[u,v].
证 分数阶微分方程周期边值问题 (PBVP) 等价于如下边值问题
设 x0 为下述边值问题的解
由于 f 连续,所以点列 {xn} 有意义且有 {xn}⊂C1−α[0,1].
以下来证明 {xn} 收敛到问题 (PBVP) 在序区间 [u,v] 中的唯一解.
证明过程分为5步.
第1步) v(t)≥x0(t)≥(u+v2)(t)≥u(t).
由定理条件 (A1),(A3),(A4),可知
第2步)对于n=0,1,2,⋯, 成立着
利用数学归纳法来证明结论.
事实上,由于
根据第1步证明及条件(A1) 有
根据引理 2.4,
设若 n=k 时结论也成立,即
则由假设 (A1),对于任意的 t∈[0,1]
重复此过程,可得
所以对所有的 n∈{0}⋃N,结论都成立.
第3步) {x2n(t)} 递增,{x2n+1(t)} 递减.
由第2步证明过程可知
类似于步骤2,利用数学归纳法可以证明 {x2n(t)} 是递增的. 同理可证{x2n+1(t)} 是递减的.
第4步) {xn(t)} 收敛到问题 (PBVP)的一个解.
由步骤1--3,对于 n=0,1,2,⋯,
设 yn(t)=xn(t)−u(t),则对 n=0,1,2,⋯,
令
由假设 (A1),(A2),可得
设 limn→∞rn=r,上述不等式表明 r≥12−r,所以 r=1.
对任意的偶数 p>0,
设 y∗=limn→∞yn,x∗=y∗(t)+u(t), 则有
第5步) x∗ 是问题 (PBVP) 在序区间 [u,v]中的唯一解.
设 ¯x∈[u,v] 是 (PBVP) 的一个解,则
据引理 2.4,有 x0(t)≤¯x(t)≤v(t).
因而
继续这一过程,容易验证
注3.1 定理3.1是一个局部性结论, 即在一个局部其解是存在且唯一的; 但从全局来看,边值问题可能会有多个解.
推论3.1 如果存在常数 c>0,λ<0 满足
1) λc2≤t1−αf(t,ct1−α)≤t1−αf(t,0)≤−λc;
2) 对于 ∀t∈[0,1],f(t,.)为凸函数;
3) 对于 ∀t∈[0,1],x,y∈R,x≤y,成立着f(t,y)−f(t,x)≤λ(y−x);
4) Eα,α(λ)<1Γ(α).
证 在推论3.1的假设下,设 u(t)≡0,v(t)=ctα−1,可以验证定理3.1的条件(A1)--(A5)都满足. 于是根据定理3.1可知周期边值问题存在唯一解 u∗ 满足 0≥u∗(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).
考察如下的边值问题