近年来, 分数阶微分方程在物理, 化学, 工程中得到了广泛应用, 分数阶微分方程的理论研究获得了快速的发展.科研工作者利用泛函分析的技巧, 如Banach压缩映像原理, Krasnoselskii's不动点定理, Schauder不动点定理获得了很多分数阶初边值问题解的存在性结果, 如文献[1-3, 18, 23-24].
2008年, Lakshmikantham[19]等, 研究了下列分数阶微分方程初值问题
其中$0 < q\leq1$, $D_{0+}^{q}$表示Riemann-Liouville分数阶导数.作者利用经典的分析方法获得了上述微分方程解的存在性.
2011年, 寇春海[20]研究了下列在右半轴上的分数阶初值问题
其中$0 < \alpha\leq1$, $D_{0+}^{\alpha}$表示Riemann-Liouville分数阶导数, $f:(0, +\infty)\times {\Bbb R} \rightarrow {\Bbb R}$.作者构造了一个新的Banach空间并利用不动点定理获得了上述方程解的存在性结果.
2012年, Baleanu[21]研究了如下分数阶初值问题
其中$0 < \alpha < \beta < 1$, $D_{0+}^{\alpha}, D_{0+}^{\beta}$表示Riemann-Liouville分数阶导数, $f:[0, 1]\times [0, +\infty)$是连续的.作者利用锥上的不动点定理得到了解与多个正解的存在性结果.
近年来, 变量阶分数次导数在工程, 物理, 金融, 水文等领域中发挥了越来越重要的作用, 试验证明它能够更好地模拟现实问题.与此同时, 变量阶的分数次微分方程受到了诸多学者的关注, 其中对于含有变量阶导数的初值问题解的存在性研究也取得了一些成果, 如文献[4-16].
2015年, 李秀英[22]研究了如下分数阶泛函边值问题的数值解
其中$1 < \alpha (x)\leq 2, {\mu }_{1}, {\mu }_{2}$是常数, $D_{0+}^{\alpha (x)}$表示变量阶Caputo导数, $\tau (x)\in C^1[0, 1], a(x), b(x)\in W^1[0, 1], u(x)\in W^4[0, 1], \tau (x)$大于或小于$x$. $D_{0+}^{\alpha (x)}$定义为
虽然分数阶微分方程解的存在性的研究结果已经很多, 但鲜有文章涉及变量阶分数次微分方程解的存在性的研究.受前面文献的启发, 此文研究如下带有变量阶导数的分数次微分方程初值问题
这里$0 < \beta < \alpha < 1$, $0 < p(t) < \alpha-\beta$, $0\leq t\leq T < +\infty$, $\lambda\in {\Bbb R}\backslash\{0\}$, 变量阶导数$D_{0+}^{p(t)}$定义为
且
上式表示指标为$1-p(t)$的变量阶积分, $D_{0+}^\alpha, D_{0+}^\beta$是Riemann-Liouville分数次导数.
很显然常微分方程或分数阶微分方程是变量阶微分方程的特殊情形.由(1.2)和(1.3)式知, 当$p(t)$是一个常值函数时, 即$p(t)\equiv q$ ($q$是一个正常数), $I_{0+}^{p(t)}, D_{0+}^{p(t)}$便是通常意义下的Riemann-Liouville分数阶积分和导数[17].
该文安排如下.第二章给出变量阶分数次导数和积分的某些性质.在第三章中, 利用Schauder不动点定理, 我们得到了初值问题(1.1)解的存在性结果.
首先我们给出分数阶算子的一些基本性质.
命题2.1[17] 如果$f\in L(0, b), 0 < b < +\infty$, 则$I_{0+}^\gamma I_{0+}^\delta f(t)=I_{0+}^{\gamma+\delta}f(t)$, $\gamma>0, \delta>0$.
命题2.2[17] 如果$f\in L(0, b), 0 < b < +\infty$, 则$D_{0+}^\gamma I_{0+}^\gamma f(t)=f(t)$, $\gamma>0$.
命题2.3[17] 设$0 < \alpha\leq 1$, 则微分方程
有解
命题2.4[17] 若$0 < \alpha\leq 1$, $u(t)\in L(0, b)$, $D_{0+}^\alpha u\in L(0, b)$, 则
命题2.5[17] Riemann-Liouville分数阶积分算子$I_{0+}^\delta : C[0, b]\rightarrow C[0, b], \delta>0$是有界的.
上述分数阶算子的性质是否仍适用于变量阶分数次算子是值得我们研究的.以命题2.1为例, 我们考虑一种特殊的情形: $f(t)\equiv 1, t\in [0, T]$.接下来, 我们将验证$I_{0+}^{p(t)}I_{0+}^{q(t)}f(t)$和$I_{0+}^{p(t)+q(t)}f(t)$是否相等.
首先, 根据(1.3)式, 有
另一方面
从而, 对$f(t)\equiv 1$, 仅当
才有下式成立
然而, 我们不能确定上式的正确性.我们仅知
因此, 对于$p(t), q(t), f(t), t\in[0, T]$, 对于$[0, T]$中的某些点, 我们有
特别地, 对于一般的函数$0 < p(t) < 1, f(t)$, 对于$[0, T]$中的某些点, 有
下面, 我们通过下例来说明(2.1)式是正确的.
例2.1 设$p(t)=t, f(t)=1, ~0\leq t\leq 3, $
下面计算$I_{0+}^{p(t)}I_{0+}^{q(t)}f(t)|_{t=1}$和$I_{0+}^{p(t)+q(t)}f(t)|_{t=1}$.
首先, 有
另一方面, 有
因此
通过上述分析, 可知命题2.1-2.4对于变量阶分数次算子而言是无效的, 从而我们不能将一个变量阶分数次微分方程转化成与之等价的积分方程.因此, 研究变量阶分数次微分方程及其初边值问题解的存在性, 虽然困难但很有意义.
接下来, 我们将给出关于变量阶分数次积分算子的一个重要结论.
引理2.1 设$x(t), p(t)$是实函数且$I_{0+}^{p(t)}x(t)$存在.若$q>0$, 则
证 由(1.3)式, 得
上式意味着(2.2)式是成立的.
注2.1 设$x(t), p(t)$是实函数且$I_{0+}^{p(t)}x(t)$存在.对于$q>0$, 则我们无法断定下式是否成立
引理2.2 设$q:[0, T]\rightarrow (0, q^{*}]$ ($0 < q^{*} < 1$)是连续的, 则对于$x(t)\in C[0, T]$, 有$I_{0+}^{1-q(t)}x(t)\in C[0, T]$.
证 由函数$q(t)$的连续性及Gamma函数, 可知$\Gamma(1-q(t))\in C[0, T]$.从而
对于$0 < T < +\infty$, 有
首先, 我们将证明
对$x\in C[0, T]$, 有
这意味着(2.6)式是成立的.
同样, 函数
是递减的.事实上, 由于$\ln a < \ln b < 0$, $0 < b^t < a^t$, 我们有
这表明$k(t)$是递减的.类似的
(这里$0 < \frac{t_1-s}T < \frac{t_2-s}T < 1$)也是递减的.
对于$t_0=0$时, 设$t>t_0=0$, $t-t_0=t\leq T$, 记
对于$\forall\varepsilon>0$, 取$\delta=\big(\frac{(1-q^{*})\varepsilon}{LMT^{*}T^{q^{*}}}\big)^{\frac 1{1-q^{*}}}$.对$t\in (0, T]$及$0 < t < \delta$, 根据(2.4)-(2.7)式, 可得
故$I_{0+}^{1-q(t)}x(t)$在$t=0$点连续.
对于$t_0\in (0, T]$, 取$t>t_0$, $t-t_0 < T$.因此, 当$x\in C[0, T]$, 对$\forall\varepsilon>0$, 取$\delta=(\frac{(1-q^{*})\varepsilon}{2LMT^{*}T^{q^{*}}})^{\frac 1{1-q^{*}}}$, $t, t_0\in [0, T]$.不失一般性, 设$t_0 < t$, 当$t-t_0 < \delta$时, 由(2.4)-(2.7)式, 可得
即$I_{0+}^{1-q(t)}x(t)$在$t_0$点连续.由$t_0$的任意性知, $I_{0+}^{1-q(t)}x(t)\in C[0, T], x\in C[0, T]$.证毕.
本章将给出本文的主要结果.现假设如下条件成立
$(H_1)$ 设$0 < \beta < \alpha < 1$, $p:[0, T]\rightarrow (0, \alpha-\beta]$为连续函数, 对于$h\in C[0, T]$, 都有$I_{0+}^{1-p(t)}I_{0+}^\beta h(t)\in AC[0, T]$;
$(H_2)$ $t^rf:[0, T]\times {\Bbb R}\times {\Bbb R}\rightarrow {\Bbb R}$ ($0\leq r < \alpha-\beta$)是连续的;
$(H_3)$ 存在正常数$c_1, c_2, c_3, 0 < u_1 < 1, 0 < u_2 < 1$, 使得
引理3.1 设
其范数定义为
则$(E, \|\cdot\|)$是一个Banach空间.
证 此引理的证明类似于文献[2]中的引理2.2.该证明我们此处省略了.
为了得到本文结果, 我们首先给出如下定理.
定理3.1 假设$(H_1), (H_2)$成立, 则
(A) 如果初值问题$(1.1)$有一个解为$x\in C[0, T]$且$D_{0+}^\beta x\in C[0, T]$, 则$z(t)=I_{0+}^{1-\alpha}x(t)\in C[0, T]$是如下积分方程的一个解
(B) 如果方程$(3.1)$有一个解为$z\in E$, 则$x(t)=I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z\in C[0, T]$是初值问题$(1.1)$的一个解.
证 假设$x\in C[0, T]$且$D_{0+}^\beta x\in C[0, T]$是初值问题(1.1)的一个解.根据(1.2)式及Riemann-Liouville分数阶导数的定义, 对方程(1.1)两边应用积分算子$I_{0+}^1$, 得
由命题2.1, 2.2, 2.4及(2.6)式, 有
设$D_{0+}^\beta x(t)=y(t)$, 由命题2.4, 得$x(t)=I_{0+}^\beta y(t)$.由$(3.2)$和命题2.1, 我们有
同样的, 设$I_{0+}^{1-\alpha+\beta}y(t)=z(t)$, 由命题2.2, 得$y(t)=D_{0+}^{1-\alpha+\beta}z(t)$.根据$(3.3)$式, 我们有
即$z(t)=I_{0+}^{1-\alpha+\beta}y(t)=I_{0+}^{1-\alpha+\beta}D_{0+}^\beta x(t)=I_{0+}^{1-\alpha}x(t)$是方程$(3.1)$的一个解.
另一方面, 如果积分方程(3.1)有一个解为$z\in E$, 则我们可以得到问题(1.1)的一个解.
事实上, 若对于方程(3.1)的一个解$z\in E$, 设
对(3.4)式两边作用算子$I_{0+}^{1-\alpha+\beta}$, 由命题2.4知
由$\beta < \alpha$及$y(t)=D_{0+}^{1-\alpha+\beta}z(t)\in C[0, T]$得$c=0$.所以
由此, 将$z(t)=I_{0+}^{1-\alpha+\beta}y(t)$代入$(3.1)$式, 得
同理, 对于$(3.6)$式, 可设
根据命题2.2, 对$(3.7)$式两边同时应用$D_{0+}^\beta$, 可得
将$y(t)=D_{0+}^\beta x(t)$带入$(3.6)$式, 由命题2.1, 可得
根据命题2.5, 知
同样的, 根据命题2.5, $D_{0+}^{1-\alpha+\beta}z\in C[0, T]$及$(H_1)$, 有
因此, 根据(3.9)及(3.10)式, 我们有$I_{0+}^{1-\alpha}x\in AC[0, T]$.对$(3.9)$式两边同时求导, 可得
根据$D_{0+}^{1-\alpha+\beta}z\in C[0, T]$我们有$x(0)=I_{0+}^\beta y(t)|_{t=0}=I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t)|_{t=0}=0$, 这意味着$x(t)=I_{0+}^\beta y(t)=I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t)$是初值问题(1.1)的一个解.
下面给出方程$(1.1)$解的存在性结论.
定理3.2 假设$(H_1)$-$(H_3)$成立且有下面的不等式成立, 则初值问题$(1.1)$至少有一个解, 这里
其中
证 由定理$3.2$可知, 我们将讨论方程$(3.1)$解的存在性.定义算子$T:E\rightarrow E$:
易知$T:E\rightarrow E$是有意义的.
首先, 我们将证明$Tz(t), D_{0+}^{1-\alpha+\beta}Tz(t)\in C[0, T], z(t)\in E$.
事实上, 设
由$(H_2)$, 可知$g:[0, T]\times {\Bbb R}\times {\Bbb R}\rightarrow{\Bbb R}$是连续的.同样, 由命题2.5, 若$D_{0+}^{1-\alpha+\beta}z\in C[0, T]$, 有$I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z\in C[0, T]$.因为$D_{0+}^{1-\alpha+\beta}z(s)$在任意点$t_0\tau\in [0, T]$($\tau\in [0, 1]$)都连续, 从而对于任意的$\varepsilon>0$, 存在$\delta_1>0$, 当$t_0\tau, t\tau\in [0, T], \tau\in [0, 1]$, $|t\tau-t_0\tau| < \delta_1$, 下式恒成立
类似的, $I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(s)$也在任意点$t_0\tau\in [0, T]$($\tau\in [0, 1]$)连续.对任意的$\varepsilon>0$, 存在$\delta_2>0$, 当$t_0\tau, t\tau\in [0, T], \tau\in [0, 1]$, $|t\tau-t_0\tau| < \delta_2$, 恒有
由于$g\big(s, I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z, D_{0+}^{1-\alpha+\beta}z\big)$在$(t_0\tau, I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t_0\tau), D_{0+}^{1-\alpha+\beta}z(t_0\tau))$ ($t_0\in [0, T], \tau\in [0, 1]$)处连续.从而, 对于任意的$\varepsilon>0$, 存在$\delta>0$, 当$t_0\tau, t\tau\in [0, T], \tau\in [0, 1]$, $|t\tau-t_0\tau| < \delta$, 结合公式(3.12)和(3.13), 我们有
另一方面, 由Riemann-Liouville分数阶导数和积分的定义, 我们有
及
因此, 通过前面的讨论, 对于$t, t_0\in [0, T]$, 使得$t>t_0$, $t-t_0 < T$, 我们有
可知$I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t\tau)$, $D_{0+}^{1-\alpha+\beta}z(t\tau)$均在$t_0\tau\in [0, T]$ $(\tau\in [0, 1])$处连续, $g\big(t\tau, I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t\tau), $ $D_{0+}^{1-\alpha+\beta}z(t\tau)\big)$在$(t_0\tau, I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t_0\tau), D_{0+}^{1-\alpha+\beta}z(t_0\tau))$ ($t_0\in [0, T]$, $\tau\in [0, 1]$)处连续.易知$t^{\alpha-\beta-r}$在$t_0\in [0, T]$处连续.继而我们有
从而, $I_{0+}^{\alpha-\beta}f\big(t, I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t), D_{0+}^{1-\alpha+\beta}z(t)\big)$在$t_0$处连续.由于$t_0$是任意的, 故有
同样的分析可知
根据前面的分析和引理2.2得$I_{0+}^{1-p(t)}I_{0+}^\beta D_{0+}^{1-\alpha+\beta}z(t)\in C[0, T]$及
从而$Tz(t)\in C[0, T], ~D_{0+}^{1-\alpha+\beta}Tz(t)\in C[0, T], z\in C[0, T]$.这意味着$T: E\rightarrow E$是有意义的.
另一方面, 易知$T: E\rightarrow E$是全连续的.
设$\Omega=\{z\in E; \|z\|\leq R\}$, 其中
对$z\in \Omega$, 有
类似的, 对$z\in E$, 有
这意味着$T(\Omega)\subseteq \Omega$.根据Schauder不动点定理知$T$在$z\in E$中至少有一个不动点, 即方程$(3.1)$至少有一个解.证毕.
注3.1 在$(H_3)$中, 若$\mu_1, \mu_2\geq 1$, 则适当的增加对$c_1, c_2, c_3$的约束条件, 利用类似于定理2.1的方法, 我们同样可以得到初值问题$(1.1)$的解.