由于分数微积分具有刻画记忆和遗传性质的特性,所以在描述具有记忆和遗传性质的各种物质、 流变学、材料和力学系统、信号处理和系统辨识、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^{\alpha}u\in C[0,1].
一个连续函数u\in C(0,1]的\alpha(0<\alpha<1)阶Riemann-Liouville导数定义为
(1)当 \beta>-1时,
D^{\alpha}(t^{\beta})=\frac{1}{\Gamma(1-\alpha)}\frac{\rm d}{{\rm d}t}\int_0^t(t-\tau)^{-\alpha}{\tau}^{\beta}{\rm d}\tau=\frac{\Gamma(\beta+1)}{\Gamma(\beta+2-\alpha)}\frac{\rm d}{{\rm d}t}(t^{1+\beta-\alpha}).
(2)对于任意的 u\in C(0,1]\bigcap L^1[0,1],D^{\alpha}I^{\alpha}u=u.
(3)当 u\in C(0,1]\bigcap L^1[0,1]时,分数微分方程 D^{\alpha}u(t)=0 有解为 u(t)=ct^{\alpha-1},c\in R.
(4)当 u\in C(0,1]\bigcap L^1(0,1] 且 D^{\alpha}u\in C(0,1]\bigcap L^1(0,1] (0<\alpha<1)时,下式成立
对于 \alpha>0,\beta>0,二元的 Mittag-Leffler 函数定义如下
为证明主要结论,需要以下的引理.
{\bf 引理2.3}[9] 设 h\in C[0,1],\alpha\in (0,1),\lambda\in {\Bbb R} 且 E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)}. 如果 u(t)\in C_{1-\alpha}[0,1] 满足
(i) D^{\alpha}u(t)-\lambda u(t)=h(t)(0<t<1);
(ii) \lim\limits_{t\rightarrow 0^+}t^{1-\alpha}u(t)=u(1).
当 0\leq t\leq s\leq1,
{\bf 注2.1}\quad 据文献[10,引理2.2],当 \alpha\in (0,1),\ \lambda\in R,\ \lambda<0时, 0<E_{\alpha,\alpha}(\lambda)< \frac{1}{\Gamma(\alpha)}. 如果E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)}, 容易验证 G_{\lambda,\alpha}(t,s)>0 for (t,s)\in (0,1)\times (0,1).
{\bf 引理2.4}[9]\quad 设 \alpha\in (0,1),\lambda\in {\Bbb R} 且 E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)}. 如果 u(t)\in C_{1-\alpha}[0,1] 满足
(i) D^{\alpha}u(t)-\lambda u(t)\geq 0(0<t<1);
在本节中,始终假定以下条件成立.
(H_1) f:(0,1]\times R\longrightarrow R 是连续函数;
(H_2) 对于 \forall x\in C_{1-\alpha}[0,1],\ f(t,x(t))\in C_{1-\alpha}[0,1], 其中
易见 (PBVP) 与如下问题等价
对于 x\in C_{1-\alpha}[0,1],定义映射 T 如下
{\bf 引理3.1}\quad 当条件 (H_1)-(H_2)成立时, 算子 T 映 C_{1-\alpha}[0,1] 到 C_{1-\alpha}[0,1] 且 T 是连续的.
{\bf 证}\quad 结论分两步证明.
(1) 若 x\in C_{1-\alpha}[0,1],则Tx\in C_{1-\alpha}[0,1].
显然当 x\in C_{1-\alpha}[0,1] 时,[f(t,x(t))-\lambda x(t)]\in C_{1-\alpha}[0,1].
设 \sigma(t)=f(t,x(t))-\lambda x(t),u(t)=(Tx)(t),则有
据文献[9],当 \sigma\in C_{1-\alpha}[0,1] 且 E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)} 时,下述线性分数微分方程边值问题
以下来证明 u\in C_{1-\alpha}[0,1].
事实上,
只需说明 t^{1-\alpha}\int_0^t {(t-s)}^{\alpha-1}E_{\alpha,\alpha}(\lambda {(t-s)^{\alpha})}\sigma(s){\rm d}s 是 t的连续函数. 注意到
因此,t^{1-\alpha}u(t)\in C[0,1].
(2) 在假设 (H_1)-(H_2)下,利用Lebesgue控制收敛定理可以证明 T 是连续的.
设 u,v\in C_{1-\alpha}[0,1] 满足 u(t)\leq v(t),t\in [0,1],记
{\bf 定理3.1}\quad 如果存在 u,v \in C_{1-\alpha}[0,1],u\leq v 和常数 \lambda (\lambda 满足E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)}), 使得以下条件成立
(A_1) 对于任意的t\in [0,1],x,y\in [u,v],x\leq y,
(A_2) 对于任意的t\in [0,1],x,y\in [u,v],
(A_3) 对于任意的 t\in (0,1),
(A_4) \lim\limits_{t\rightarrow 0^+}t^{1-\alpha}u(t)=u(1),\ \lim\limits_{t\rightarrow 0^+}t^{1-\alpha}v(t)=v(1).
那么(PBVP) 存在唯一解 x^*\in C_{1-\alpha}[0,1]\cap [u,v].
{\bf证}\quad 分数阶微分方程周期边值问题 (PBVP) 等价于如下边值问题
设 x_0 为下述边值问题的解
由于 f 连续,所以点列 \{x_n\} 有意义且有 \{x_n\}\subset C_{1-\alpha}[0,1].
以下来证明 \{x_n\} 收敛到问题 (PBVP) 在序区间 [u,v] 中的唯一解.
证明过程分为5步.
第1步) v(t)\geq x_0(t)\geq (\frac{u+v}{2})(t)\geq u(t).
由定理条件 (A1),(A3),(A4),可知
第2步)对于n=0,1,2,\cdots , 成立着
利用数学归纳法来证明结论.
事实上,由于
根据第1步证明及条件(A1) 有
根据引理 2.4,
设若 n=k 时结论也成立,即
则由假设 (A1),对于任意的 t\in [0,1]
重复此过程,可得
所以对所有的 n\in \{0\}\bigcup N,结论都成立.
第3步) \{x_{2n}(t)\} 递增,\{x_{2n+1}(t)\} 递减.
由第2步证明过程可知
类似于步骤2,利用数学归纳法可以证明 \{x_{2n}(t)\} 是递增的. 同理可证\{x_{2n+1}(t)\} 是递减的.
第4步) \{x_{n}(t)\} 收敛到问题 (PBVP)的一个解.
由步骤1--3,对于 n=0,1,2,\cdots,
设 y_n(t)=x_n(t)-u(t),则对 n=0,1,2,\cdots,
令
由假设 (A1),(A2),可得
设 \lim\limits_{n\rightarrow \infty}r_n=r,上述不等式表明 r\geq \frac{1}{2-r},所以 r=1.
对任意的偶数 p>0,
设 y^*=\lim\limits_{n\rightarrow \infty}y_{n},x^*=y^*(t)+u(t), 则有
第5步) x^* 是问题 (PBVP) 在序区间 [u,v]中的唯一解.
设 \overline{x}\in [u,v] 是 (PBVP) 的一个解,则
据引理 2.4,有 x_0(t)\leq \overline{x}(t)\leq v(t).
因而
继续这一过程,容易验证
{\bf 注3.1}\quad 定理3.1是一个局部性结论, 即在一个局部其解是存在且唯一的; 但从全局来看,边值问题可能会有多个解.
{\bf 推论3.1}\quad 如果存在常数 c>0,\lambda<0 满足
1) \frac{\lambda c}{2}\leq t^{1-\alpha}f(t,c t^{1-\alpha})\leq t^{1-\alpha}f(t,0)\leq -\lambda c;
2) 对于 \forall t\in [0,1] ,f(t,.)为凸函数;
3) 对于 \forall t\in [0,1],x,y\in R,x\leq y,成立着f(t,y)-f(t,x)\leq \lambda (y-x);
4) E_{\alpha,\alpha}(\lambda)<\frac{1}{\Gamma(\alpha)}.
{\bf 证}\quad 在推论3.1的假设下,设 u(t)\equiv 0,v(t)=ct^{\alpha-1},可以验证定理3.1的条件(A1)--(A5)都满足. 于是根据定理3.1可知周期边值问题存在唯一解 u^* 满足 0 \geq {u^*}(t)\geq ct^{\alpha-1},t\in [0,1].
进一步,如果存在 t_0\in(0,1] 使得 f(t_0,0)\neq 0,那么 u^*(t) 在(0,1)内不会恒为0,从而根据引理2.3得 u^*(t)>0,\ t\in(0,1).
考察如下的边值问题