孙文超, 苏有慧, 孙 爱

(1. 沈阳工业大学 理学院, 沈阳 110870; 2. 徐州工程学院 数学与统计学院, 江苏 徐州 221018)

0 引 言

微分方程问题在物理、 化学、 生物和工程等领域应用广泛[1-5], 目前, 关于分数阶微分方程解的存在性研究已有很多结果[6-12]. 但实际问题中抽象出的方程通常可能是较复杂的且含有积分的分数阶微分方程, 如扩散过程、 封闭系统中种群增长的Volterra模型及控制理论等[13-16]. 因此, 研究分数阶积分微分方程尤为重要.

文献[14]研究了如下分数阶积分微分方程正解的存在性:

得到了其解的存在性和唯一性结果, 其中:n-1<α≤n, 0<μ

使用Krasnoselskii不动点定理, 得到了其至少存在一个解的结果, 其中: 1<α≤2, 0<β≤1;f: [0,1]×3→是连续的; 定义映射但上述结果都只从理论方面研究了分数阶积分微分解的存在性, 关于其数值解方面的研究目前报道较少[17-18].基于此, 本文考虑如下非线性分数阶积分微分方程:

(1)

且满足

通过使用Banach压缩映射原理和Krasnoselskii不动点定理, 给出分数阶积分微分方程(1)至少存在一个解和存在唯一解的结论.特别地, 通过模拟仿真实例说明本文主要结果的正确性.

1 预备知识

定义1[19]令α>0, 函数u: (0,+∞)→的Riemann-Liouville分数阶积分定义为

其中等式右端在(0,+∞)上有定义.

定义2[19]令α>0, 函数u: (0,+∞)→的Riemann-Liouville分数阶导数定义为

其中n-1≤α

引理1[20]令α>0, 则有分数阶微分方程

这里Ci∈,i=1,2,…,n,n=[α]+1.

引理2假设y∈C[0,1]且n-1<α1,α-β≥1, 0<η<1,k0>k1>0, 则边值问题

(2)

有唯一解

(3)

这里m=k0-k1ηα-β-1>0.

证明: 根据引理1, 有

由u(0)=u′(0)=…=u(n-2)(0)=0, 可得

C2=C3=…=Cn-2=0,

因此

(4)

进一步可得

(5)

将式(5)代入式(4)即可推出式(3).

引理3(Krasnoselskii不动点定理)[21]设X是Banach空间E中的非空闭子集, 映射A,B满足如下条件：

1) 当x,y∈X时,Ax+By∈X;

2)A是连续且紧的;

3)B是一个压缩映射.

则存在z∈X, 使得z=Az+Bz.

2 主要结果

假设:

(H3) 存在一个常数b∈(0,α-n+1)和非负函数g(t)∈L1/b([0,1],(0,+∞)), 使得对任意的t∈[0,1],u,x,y,z,u′,x′,y′,z′∈, 均有

|f(t,u,x,y,z)-f(t,u′,x′,y′,z′)|≤g(t)(|u-u′|+|x-x′|+|y-y′|+|z-z′|).

定理1假设条件(H1)～(H3)成立, 且存在一个常数ρ, 如果

(6)

则边值问题(1)有唯一解.

证明: 定义算子A:X→X为

令ΩR={u∈X: ‖u‖

1) 对任意的u∈ΩR……