分数阶微分方程在工程、化学、物理、生物等领域有着广泛应用,例如热传导领域和流体学领域[1-3], 而且分数阶导数模型克服了经典整数阶微分模型理论与实验结果不吻合的缺点[4], 因此研究分数阶微分方程边值问题有着重要的意义.近年来,大量文献报道微分方程[5-6]和分数阶微分方程[4, 7-10]边值问题解的存在性.2005年,当1 < α≤2时,Bai等[7]推导了分数阶微分方程边值问题
$ \left\{ {\begin{array}{*{20}{l}} {D_{0 + }^\alpha u\left( t \right) + f\left( {t, u\left( t \right)} \right) = 0, t \in \left( {0, 1} \right), }\\ {u\left( 0 \right) = u\left( 1 \right) = 0, } \end{array}} \right. $ |
的格林函数及重要的性质(f:[0, +∞)×
$ \left\{ {\begin{array}{*{20}{l}} {D_{0 + }^\alpha u\left( t \right) + f\left( {t, u\left( t \right)} \right) = 0, t \in \left( {0, 1} \right), }\\ {u\left( 0 \right) = u'\left( 1 \right) = u''\left( 0 \right) = 0, } \end{array}} \right. $ |
至少存在一个正解,其中f:[0, +∞)×
$ \left\{ {\begin{array}{*{20}{l}} {u''\left( t \right) + f\left( {u\left( t \right)} \right) = 0, t \in \left[{0, 1} \right], }\\ {u'\left( 0 \right) = u\left( 1 \right) = 0, } \end{array}} \right. $ |
存在正解,其中f:
目前研究分数阶微分方程边值问题的主要工具有锥拉伸与锥压缩不动点原理、Krasnoselskii不动点原理、Schauder不动点原理上下解等.本文利用Leggett-Williams不动点定理, 参照文献[9]中的方法研究Caputo分数阶微分方程边值问题
$ \left\{ {\begin{array}{*{20}{l}} {D_{0 + }^\alpha u\left( t \right) + f\left( {t,u\left( t \right)} \right) = 0,t \in \left( {0,1} \right),}\\ {u'\left( 0 \right) = u\left( 1 \right) = 0,} \end{array}} \right. $ | (0.1) |
正解的存在性,其中1 < α≤2,f:[0, +∞)×
定义1.1[11]一个连续函数u:(0, +∞)→
$ {}^CD_{0 + }^\alpha u\left( t \right) = \frac{1}{{\mathit{\Gamma} \left( {n-\alpha } \right)}}\int_0^t {\frac{{{u^{\left( n \right)}}\left( s \right)}}{{{{\left( {t-s} \right)}^{\alpha-n + 1}}}}} {\rm{d}}s, $ |
其中α>0,n=[α]+1,[α]代表实数α的整数部分。上式右边在(0, +∞)内逐点有定义.
引理1.1[11]令α>0,若u∈ACn[0,1]或u∈Cn[0,1],则
$ I_{0 + }^\alpha {}^CD_{0 + }^\alpha u\left( t \right) = u\left( t \right)-\sum\limits_{k = 0}^{n-1} {\frac{{{u^{\left( k \right)}}\left( 0 \right)}}{{k!}}{t^k}}, $ |
其中n=[α]+1,I0+α代表α阶Riemann-Liouville型积分.
引理1.2令α∈(1, 2],给定h∈C[0,1],则
$ {}^CD_{0 + }^\alpha u\left( t \right) + h\left( t \right) = 0,t \in \left( {0,1} \right), $ | (1.1) |
$ 2u'\left( 0 \right) = u\left( 1 \right) = 0, $ | (1.2) |
的唯一解为
证明由引理1.1可得,(1.1)式等价于方程u(t)=-I0+αh(t)+C1+C2t,其中C1, C2∈
因此,(1.1)~(1.2)式的唯一解是
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\;\;u\left( t \right) = \frac{{- 1}}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^t {{{\left( {t- s} \right)}^{\alpha- 1}}h\left( s \right){\rm{d}}} s + \frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^1 {\left( {1 - } \right.} }\\ {{{\left. s \right)}^{\alpha - 1}}h\left( s \right){\rm{d}}s = \frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^t {\left[{{{\left( {1-s} \right)}^{\alpha-1}}-{{\left( {t - s} \right)}^{\alpha - 1}}} \right]h\left( s \right){\rm{d}}} s + }\\ {\frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_t^1 {{{\left( {t -s} \right)}^{\alpha -1}}h\left( s \right){\rm{d}}} s = \int_0^1 {G\left( {t, s} \right)h\left( s \right)} {\rm{d}}s.} \end{array} $ |
引理1.3引理1.2中的G(t, s)有下列性质:
(ⅰ)G(t, s)∈C([0,1]×[0,1],
(ⅱ)
(ⅲ)
(ⅳ)
证明(i)~(iii)显然可得,只需证明(iv).
令
$ \mathop {\max }\limits_{t \in \left[{\frac{1}{r}, 1-\frac{1}{r}} \right]} G\left( {t, s} \right) = \left\{ {\begin{array}{*{20}{l}} {{g_1}\left( {\frac{1}{r}, s} \right), s \in \left( 0 \right., \left. {\frac{1}{r}} \right], }\\ {{g_2}\left( {t, s} \right), s \in \left[{\frac{1}{r}, 1-\frac{1}{r}} \right], }\\ {{g_2}\left( {t, s} \right), s \in \left[{1-\frac{1}{r}} \right., \left. 1 \right), } \end{array}} \right. $ |
$ \mathop {\min }\limits_{t \in \left[{\frac{1}{r}, 1-\frac{1}{r}} \right]} G\left( {t, s} \right) = \left\{ {\begin{array}{*{20}{l}} {{g_1}\left( {1 -\frac{1}{r}, s} \right), s \in \left( 0 \right., \left. {\frac{1}{r}} \right], }\\ {{g_1}\left( {1 - \frac{1}{r}, s} \right), s \in \left[{\frac{1}{r}, 1-\frac{1}{r}} \right], }\\ {{g_2}\left( {t, s} \right), s \in \left[{1-\frac{1}{r}} \right., \left. 1 \right), } \end{array}} \right. $ |
从而
$ M = \frac{{{{\left( {1-\frac{1}{r}} \right)}^{\alpha-1}}}}{{{{\left( {1-\frac{1}{r}} \right)}^{\alpha - 1}} - {{\left( {1 - \frac{1}{{2r}}} \right)}^{\alpha - 1}}}}. $ |
又
$ \begin{array}{*{20}{l}} {\mathop {\max }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} G\left( {t, s} \right) =\\ \left\{ {\begin{array}{*{20}{l}} {{g_1}\left( {{t_1}, s} \right), s \in \left( 0 \right., \left. {{t_1}} \right]}\\ {{g_2}\left( {t, s} \right), s \in \left[{1- {t_2}, \left. 1 \right)} \right.} \end{array}, } \right.}\\ {\mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} G\left( {t, s} \right) =\\ \left\{ {\begin{array}{*{20}{l}} {{g_1}\left( {1 -{t_1}, s} \right), s \in \left( 0 \right., \left. {1 -{t_2}} \right]}\\ {{g_2}\left( {t, s} \right), s \in \left[{1-{t_2}, \left. 1 \right)} \right.} \end{array}, } \right.} \end{array} $ |
从而
令γ, β, θ是锥P上的非负连续凸函数,α, ψ是锥P上的非负连续凹函数,那么对非负实数h, a, b, d和c,定义下列凸集:
P(γ, c)={u∈P:γ(u) < c},P(γ, α, a, c)={u∈P:a≤α(u), γ(u)≤c},
Q(γ, β, d, c)={u∈P:β(u)≤d, γ(u)≤c},
P(γ, θ, α, a, b, c)={u∈P:a≤α(u), θ(u)≤b, γ(u)≤c},
Q(γ, β, ψ, h, d, c)={u∈P:h≤ψ(u), β(u)≤d, γ(u)≤c}.
定理1.1[12]令E是一个实Banach空间,且P⊂E是一个锥.假设存在正数c和M,使锥P上的非负连续凹函数α, ψ及非负连续凸函数γ, β, θ满足
α(u)≤β(u), ‖u‖≤γ(u), u∈P(γ, c).
若F:P(γ, c)→P(γ, c)是一个全连续算子且存在非负实数h, d, a, b, 0 < d < a使得
(B1){u∈P(γ, θ, α, a, b, c):α(u)>a}≠∅且α(F(u)) > a, u∈P(γ, θ, α, a, b, c);
(B2){u∈Q(γ, β, ψ, h, d, c):β(u) < d}≠∅且β(F(u)) < d, u∈Q(γ, β, ψ, h, d, c);
(B3)若u∈P(γ, α, a, c)且θ(F(u))>b,则α(F(u))>a;
(B4)若u∈Q(γ, β, d, c)且ψ(F(u)) < h,则β(F(u)) < d.
那么F至少有3个不动点u1, u2, u3∈P(γ, c), 使得
β(u1) < d, a < α(u2), d < β(u3), α(u3) < a.
2 主要结果令E=C[0,1],其范数为
$ \begin{array}{*{20}{l}} {\alpha \left( u \right) = \mathop {\min }\limits_{t \in \left[ {{t_1},{t_2}} \right] \cup \left[ {1 - {t_2},1 - {t_1}} \right]} u\left( t \right),\beta \left( u \right) = }\\ {\mathop {\max }\limits_{t \in \left[ {1/r,1 - 1/r} \right]} u\left( t \right)\gamma \left( u \right) = \mathop {\max }\limits_{t \in \left[ {0,{t_3}} \right] \cup \left[ {1 - {t_3},1} \right]} u\left( t \right),\theta \left( u \right) = }\\ {\mathop {\max }\limits_{t \in \left[ {{t_1},{t_2}} \right] \cup \left[ {1 - {t_2},1 - {t_1}} \right]} u\left( t \right),\psi \left( u \right) = \mathop {\min }\limits_{t \in \left[ {1/r,1 - 1/r} \right]} u\left( t \right).} \end{array} $ |
其中t1, t2和r是非负的实数满足
$ \alpha \left( u \right) \le \beta \left( u \right), $ | (2.1) |
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\gamma \left( u \right) = \mathop {\max }\limits_{t \in \left[ {0,{t_3}} \right] \cup \left[ {1 - {t_3},1} \right]} u\left( t \right) \ge \mathop {\max }\limits_{t \in \left[ {0,{t_3}} \right]} u\left( t \right) = }\\ {\mathop {\max }\limits_{t \in \left[ {0,{t_3}} \right]} \left( {\int_0^{{t_3}} {G\left( {t,s} \right)f\left( {s,u\left( s \right)} \right){\rm{d}}s + } \int_{{t_3}}^{1 - {t_3}} {G\left( {t,s} \right)f\left( {s,} \right.} } \right.}\\ {\left. {\left. {u\left( s \right)} \right){\rm{d}}s + \int_{1 - {t_3}}^1 {G\left( {t,s} \right)f\left( {s,u\left( s \right)} \right){\rm{d}}s} } \right) = }\\ {\frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^1 {{{\left( {1 - s} \right)}^{\alpha - 1}}} f\left( {s,u\left( s \right)} \right){\rm{d}}s = \left\| u \right\|.} \end{array} $ | (2.2) |
定理2.1假设存在非负实数a, b和c使得
(H1)
(H2)
(H3)
那么,边值问题(1)至少有3个正解u1, u2和u3, 满足
$ \begin{array}{*{20}{l}} {\mathop {\max }\limits_{t \in \left[{0, {t_3}} \right] \cup \left[{1-{t_3}, 1} \right]} {u_i} \le c, i = 1, 2, 3, }\\ {\mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} {u_1}\left( t \right) > b, \mathop {\max }\limits_{t \in \left[{1/r, 1-1/r} \right]} {u_2}\left( t \right) < a, }\\ {\mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} {u_3}\left( t \right) < b, \mathop {\max }\limits_{t \in \left[{1/r, 1-1/r} \right]} {u_3}\left( t \right) > a.} \end{array} $ |
证明在锥P上定义算子A为
$ Au\left( t \right) = \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s. $ |
因为
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\;Au\left( t \right) = \frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^t {\left( {{{\left( {1-s} \right)}^{\alpha-1}}-{{\left( {t - s} \right)}^{\alpha - 1}}} \right)} f\left( {s, } \right.}\\ {\left. {u\left( s \right)} \right){\rm{d}}s + \frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_t^1 {{{\left( {1 - s} \right)}^{\alpha - 1}}f\left( {s, u\left( s \right)} \right){\rm{d}}s} \ge }\\ {\frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_0^t {{{\left( {1 - s} \right)}^{\alpha - 1}}f\left( {s, u\left( s \right)} \right){\rm{d}}s} + \frac{1}{{\mathit{\Gamma} \left( \alpha \right)}}\int_t^1 {{{\left( {1 - s} \right)}^{\alpha - 1}}f\left( {s, } \right.} }\\ {\left. {u\left( s \right)} \right){\rm{d}}s = \left\| {Au} \right\|, } \end{array} $ |
所以A:P→P连续.由Arzela-Ascoli定理易证A:P→P是全连续的.
首先,对任意u∈P,由(2.1)式和(2.2)式可知α(u)≤β(u), ‖u‖≤γ(u).若u∈P(γ, c),则‖u‖≤c.又由(H3)得
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\;\gamma \left( {Au} \right) = \mathop {\max }\limits_{t \in \left[{0, {t_3}} \right] \cup \left[{1-{t_3}, 1} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s \le }\\ {c\alpha \mathit{\Gamma} \left( \alpha \right)\int_0^1 {G\left( {s, s} \right){\rm{d}}s} = c, } \end{array} $ |
因此,A:P(γ, c)→P(γ, c).说明{u∈P(γ, θ, α, a,
(1)若u∈Q(γ, β, a, c)且
$ \begin{array}{*{20}{l}} {\;\;\;\;\beta \left( {Au} \right) = \mathop {\max }\limits_{t \in \left[{1/r, 1-1/r} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s \le }\\ {M\mathop {\min }\limits_{t \in \left[{1/r, 1-1/r} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right){\rm{d}}s = M\psi \left( {Au} \right)} < a;} \end{array} $ |
(2)若u∈Q(γ, β, ψ, $ \frac{a}{M} $, a, c),由(H1)得
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\;\beta \left( {Au} \right) = \mathop {\max }\limits_{t \in \left[{1/r, 1-1/r} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s < }\\ {\frac{a}{{\mathit{\Gamma} \left( \alpha \right)\left( {1 - \frac{1}{r}} \right)\alpha - 1}}\mathop {\max }\limits_{t \in \left[{1/r, 1-1/r} \right]} \int_0^1 {G\left( {t, s} \right)} {\rm{d}}s = a;} \end{array} $ |
(3)若u∈Q(γ, α, b, c),且θ(A(u))>
$ \begin{array}{*{20}{l}} {\;\;\;\;\alpha \left( {Au} \right) = \mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s \ge m}\\ {\mathop {\max }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, u\left( s \right)} \right)} {\rm{d}}s = m\theta \left( {Au} \right) > b;} \end{array} $ |
(4)若u∈Q(γ, θ, b,
$ \begin{array}{*{20}{l}} {\;\;\;\;\;\alpha \left( {Au} \right) = \mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} \int_0^1 {G\left( {t, s} \right)f\left( {s, } \right.} }\\ {\left. {u\left( s \right)} \right){\rm{d}}s > \frac{{\alpha b\mathit{\Gamma} \left( \alpha \right)}}{{1 - {{\left( {1 - {t_1}} \right)}^\alpha }}}\mathop {\min }\limits_{t \in \left[{{t_1}, {t_2}} \right] \cup \left[{1-{t_2}, 1-{t_1}} \right]} \int_0^1 {G\left( {t, s} \right) \cdot } }\\ {{\rm{d}}s = b;} \end{array} $ |
最后,由定理1.1可得,边值问题(0.1)至少有3个正解u1, u2, u3∈P(γ, c), 满足
α(u1) > b, β(u2) < a, α(u3) < b, β(u3) > a.
3 结论本文研究了一类Caputo分数阶微分方程边值问题多解的存在性.证明时,将微分方程边值问题转化为积分方程,进一步转化为讨论积分算子不动点的问题,然后通过运用Leggett-Williams不动点定理该分数阶微分方程边值问题至少有3个正解存在的结果,其中格林函数的性质和非线性项的条件至关重要.
[1] |
PODLUBNY I. Fractional Differential Equations, Mathematics in Science and Engineering[M]. New York: Academic Press, 1999.
|
[2] |
ADOMIAN G, ELROD M, RACH R. A new approach to boundary value equations and application to a generalization of Airy's equation[J]. J Math Anal Appl, 1989, 140(2): 554-568. DOI:10.1016/0022-247X(89)90083-8 |
[3] |
AGARWAL R P, MEEHAN M, O'REGAN D. Fixed Point Theory and Applications[M]. Cambridge: Cambridge University Press, 2001.
|
[4] |
ABDELJAWAD T, BALEANU D. Fractional differen-ces and integration by parts[J]. Journal of Computational Analysis and Applications, 2011, 13(3): 574-582. |
[5] |
王勇, 韦煜明. 二阶非线性时滞微分方程边值问题正解的存在性[J]. 广西科学, 2012, 19(1): 40-43. WANG Y, WEI Y M. Existence of positive solutions for boundary value problems of nonlinear second-order delay differential equations[J]. Guangxi Sciences, 2012, 19(1): 40-43. |
[6] |
严建明. 中立型微分方程的正解存在性及非振动解的渐近性[J]. 广西科学, 2008, 15(1): 7-9. YAN J M. Existence of asymptotic behavious of positive solution of neutral differential equation[J]. Guangxi Sciences, 2008, 15(1): 7-9. |
[7] |
BAI Z B, LU H S. Positive solutions for boundary value problem of nonlinear fractional differential equation[J]. Journal of Mathematical Analysis and Applications, 2005, 311(2): 495-505. DOI:10.1016/j.jmaa.2005.02.052 |
[8] |
BAI Z B, QIU T T. Existence of positive solution for singular fractional differential equation[J]. Applied Mathematics and Computation, 2009, 215(7): 2761-2767. DOI:10.1016/j.amc.2009.09.017 |
[9] |
DOGAN A. On the existence of positive solutions for the second-order boundary value problem[J]. Applied Mathematics Letters, 2015, 49: 107-112. DOI:10.1016/j.aml.2015.05.004 |
[10] |
XIE W Z, XIAO J, LUO Z G. Existence of extremal solutions for nonlinear fractional differential equation with nonlinear boundary conditions[J]. Applied Mathematics Letters, 2015, 41: 46-51. DOI:10.1016/j.aml.2014.10.014 |
[11] |
OLDHAM K B, SPANIER J. The Fractional Calculus[M]. New York: Academic Press, 1974.
|
[12] |
AVERY R I. A generalization of the Leggett-Williams fixed point theorem[J]. Math Sci Res Hot-Line, 1999, 3(7): 9-14. |