一类非线性Caputo型分数阶微分方程耦合系统的正解

合集下载

一类非线性分数阶微分方程边值问题解的存在性

一类非线性分数阶微分方程边值问题解的存在性

一类非线性分数阶微分方程边值问题解的存在性张海丽【摘要】利用锥拉伸和压缩不动点定理研究了非线性分数阶微分方程边值问题:{cDα0+u(t)=λf(t,u(t),u'(t)),0<t<1 u(0)=u'(0)=u(1)=0得到解的存在性,其中0<λ<1,2<α≤3,f∈C([0,1]×R2,R),cDDα0+为Caputo型分数阶微分.【期刊名称】《宜宾学院学报》【年(卷),期】2017(017)012【总页数】5页(P78-81,91)【关键词】分数阶微分方程;边值问题;不动点定理【作者】张海丽【作者单位】山西交通职业技术学院基础部,山西太原030031【正文语种】中文【中图分类】O175.8近年来,越来越多的学者开始研究分数阶微分方程,他们发现用分数阶微分方程能更准确地描述事物的变化规律,在众多领域都有着广泛的应用.许多学者研究了分数阶微分方程边值问题解的存在性与唯一性,获得了很多成果[1-6].文献[1]运用Banach压缩映象原理证明了分数阶微分方程边值问题解的存在唯一性.文献[2]运用上下解方法研究了分数阶微分方程边值问题得到正解的存在性.文献[3]运用Schauder不动点定理及Krasnoselskii’s不动点定理研究了非线性分数阶微分方程边值问题得到正解的存在性,其中是Caputo型分数阶微分.本文讨论非线性分数阶微分方程边值问题:用锥拉伸和压缩不动点定理得到解的存在性.其中是Caputo型分数阶微分,f∶[0,1]×[0,+∞)×R→[0,+∞)连续.定义1若函数u∈ACn([0,1],R),且α>0,则u的α阶Caputo型分数阶导数定义为:其中:当α为非正整数时,n=[α]+1(其中[α]为α的整数部分),当α为正整数时,n=α.引理1[7]假设E是Banach空间,P为E中的锥,Ω1,Ω2为E中两个开集,且满足θ∈Ω1,,算子全连续,且下列条之一成立,则算子A在中至少有一个不动点. 引理2[3]令2<α≤3,y∈C[0,1],则分数阶微分方程存在唯一解:其中引理3由式(4)'定义的函数满足:证明:①由定义可知②当0≤s≤t≤1时当0≤t≤s≤1时③当0<ξ1≤s≤t≤ξ2<1时当0<ξ1≤t≤s≤ξ2<1时记则X的范数为方便,记且定义算子T∶X→X如下:这里G(t,s)由式(4)定义,则边值问题(1)在Ya中的解等价于Tu=u在X中的不动点.引理4若u∈Ya是边值问题(1)的一个正解,则其满足:其中证明:若u∈Ya是边值问题(1)的一个正解,则由注1和引理3中②有即另一方面,利用引理3中③和式(6)有于是利用(6)和(7)有引理5算子T∶X→X是全连续.证明:先证算子T∶X→X.易知Tu(t)≥0,∀t∈[0,1],对∀u∈X,由引理3,类似于式(6)和(7)可知即因此算子T∶X→X.下证T是全连续算子.由于函数G,f都是连续的,故T也是连续的.令是X中的有界集,记首先,对任意的u∈Br,由式(5)和引理3知:因此算子T在Br上一致有界.其次,设0≤t1<t2≤1,对任意的u∈Br,有当t1→t2时:因此,算子T在Br上等度连续,由Ascoli-Arzela定理知,算子T∶X→X是全连续的.为方便,引入如下记号:令下列定理成立.定理1若f0∈[0,r)且f∞∈(R,+∞],则边值问题(1)在X中至少有一个正解.证明:只需证明算子T在X中至少有一个不动点,由假设f0∈[0,r),则存在μ1>0,及充分小的ε1>0,使得其中f0+ε1≤r.令对∀u∈∂Ω1⋂X,利用(5)和(8)可知:即另一方面,由f∞∈(R,+∞],则存在l>μ1>0及充分小的ε2>0,使得其中f∞-ε2≥R.令则对∀u∈∂Ω2⋂X,由引理4可知由式(5)和(10)可知:从而由式(9)和(11)及引理1可知,算子T在上至少有一个不动点u,即边值问题(1)在X中至少有一个正解u.定理2若f∞∈[0,r)且f0∈(R,+∞],则边值问题(1)在X中至少有一个正解.同理可证,略.参考文献:[1]BAI Z B.On positive solutions of a nonlocal fractional boundary value problem[J].Nonlinear Analysis,2010,72(2):916-924.[2]ZHAO Y G,SUN S R,HAN Z L et al.The existence of multiple positive solutions for boundary value problems of nonlinear fractional diffusion equations[J].Communications in Nonlinear Science and Numerical Simulation,2011,16(4):2086-2097.[3]XIONG Y,WEI Z L,WEI D.Existence of positive solutions for the boundary value problem of nonlinear fractional differentialequations[J].Communications in Nonlinear Science and Numerical Simulation,2012,17(1):85-92.[4]SUN Y,ZHAO M.Positive solutions for a class of fractional differential equations with integral boundary conditions[J].Applied Mathematics Letters,2014(34):17-21.[5]WANG J R,ZHANG Y R.A class of nonlinear differential equations with fractional integrable impulses[J].Communications in Nonlinear Science and Numerical Simulation,2014,19(9):3001-3010.[6]LIU Z H,LI X W.Existence and uniqueness of solutions for the nonlinear impulsive fractional differential equation[J].Communications in NonlinearScience and Numerical Simulation,2013,18(6):1362-1373. [7]郭大钧.非线性泛函分析[M].第二版.济南:山东科学技术出版社,2004.。

Caputo分数阶微分方程解的唯一性

Caputo分数阶微分方程解的唯一性

Caputo分数阶微分方程解的唯一性
徐超宇;王颖;訾玉梅
【期刊名称】《应用数学进展》
【年(卷),期】2024(13)4
【摘要】本文主要研究一类具有Riemann-Stieltjes边值条件的Caputo分数阶微分方程。

利用Green函数的性质,Banach收缩原理,证明了方程解的唯一性。

【总页数】7页(P1210-1216)
【作者】徐超宇;王颖;訾玉梅
【作者单位】临沂大学数学与统计学院临沂
【正文语种】中文
【中图分类】O17
【相关文献】
1.状态依赖脉冲Caputo分数阶微分方程解的存在唯一性
2.带积分边界条件的非线性Caputo型分数阶微分方程边值问题解的存在唯一性
3.Caputo型分数阶微分系统正解的唯一性
因版权原因,仅展示原文概要,查看原文内容请购买。

《两类非线性时间分数阶耦合方程组的Galerkin混合元方法研究》范文

《两类非线性时间分数阶耦合方程组的Galerkin混合元方法研究》范文

《两类非线性时间分数阶耦合方程组的Galerkin混合元方法研究》篇一一、引言近年来,随着科学计算与数值分析的飞速发展,分数阶微分方程因其能够描述多种复杂的物理现象而备受关注。

特别地,时间分数阶微分方程在描述流体流动、扩散、波传播等过程中表现出独特的优势。

本文将重点研究两类非线性时间分数阶耦合方程组,并采用Galerkin混合元方法进行求解。

二、问题描述1. 第一类非线性时间分数阶耦合方程组主要涉及流体动力学、多孔介质中的物质扩散等问题;2. 第二类非线性时间分数阶耦合方程组则用于描述其他物理过程,如化学反应过程等。

这两类方程均具有较强的非线性和时变特性,因此给求解带来了极大的挑战。

三、Galerkin混合元方法简介Galerkin混合元方法是一种常用的数值求解方法,其基本思想是利用有限维空间中的函数集来逼近真实解。

该方法具有较高的求解精度和灵活性,适用于求解各类复杂的微分方程问题。

在处理时间分数阶耦合方程组时,Galerkin混合元方法可以有效地处理方程组中的非线性和耦合特性。

四、两类非线性时间分数阶耦合方程组的Galerkin混合元方法1. 第一类非线性时间分数阶耦合方程组的Galerkin混合元方法针对第一类非线性时间分数阶耦合方程组,我们首先构建合适的有限元空间和混合元空间。

然后,利用Galerkin混合元方法,将原问题转化为一系列的线性或非线性代数方程组。

接着,采用迭代法或直接法求解这些代数方程组,从而得到原问题的近似解。

2. 第二类非线性时间分数阶耦合方程组的Galerkin混合元方法对于第二类非线性时间分数阶耦合方程组,我们同样采用Galerkin混合元方法进行求解。

首先,根据问题的特点,选择合适的有限元空间和混合元空间。

然后,利用Galerkin混合元方法将原问题转化为一系列的代数方程组。

在求解过程中,我们需要注意处理方程组中的非线性和耦合特性,以确保求解的准确性和效率。

五、数值实验与结果分析为验证Galerkin混合元方法在求解两类非线性时间分数阶耦合方程组中的有效性,我们进行了大量的数值实验。

无穷区间上分数阶耦合微分系统积分边值问题正解的存在性

无穷区间上分数阶耦合微分系统积分边值问题正解的存在性

无穷区间上分数阶耦合微分系统积分边值问题 正解的存在性鄢
许文序袁 周宗福鄢鄢
渊 安徽大学 数学科学学院袁合肥 圆猿园远园员冤
摘摇 要院分数阶微积分被广泛应用于流体力学尧电化学分析尧生物系统的电传导等领域袁分数阶微分方 程的边值问题已成为研究热点袁无限区间上的边值问题是其中比较困难的部分袁针对这种边值问题袁提出了 一类无穷区间上具有积分边界条件的分数阶耦合微分方程曰应用格林函数及分数阶微积分的有关结论袁将 这类无穷区间上具积分边界条件的分数阶耦合微分方程边值问题转化为等价的积分系统曰引入函数乘积空 间和二维积分算子袁借助锥上 运则葬泽灶燥泽藻造泽噪蚤蚤 不动点定理袁并利用一些分析技巧袁得到此边值问题至少存在一 个正解的充分条件袁建立了无限区间上分数阶耦合边值问题正解存在性的新结果遥
第灾燥猿造郾远猿远卷摇 第晕韵员 郾期员 摇








重庆工商大学学报渊 自然科学版冤
悦澡燥灶早择蚤灶早 栽藻糟澡灶燥造 月怎泽蚤灶藻泽泽 哉灶蚤增郾 渊 晕葬贼 杂糟蚤
耘凿冤 摇
摇摇摇摇摇摇

圆园员怨 年 圆 月 云藻遭郾 圆园员怨
凿燥蚤院员园郾 员远园缘缘 辕 躁郾 蚤泽泽灶郾 员远苑圆原园缘愿载郾 圆园员怨郾 园园园员郾 园园圆
关键词院无穷区间曰分数阶耦合微分系统曰锥上不动点定理曰积分边界条件 中图分类号院韵员苑缘郾 愿摇 摇 摇 文献标志码院粤摇 摇 摇 文章编号院员远苑圆原园缘愿载渊圆园员怨冤园员原园园园远原园苑
园摇 引摇 言
分数微积分是整数微积分的一种推广援 在实际 问题中袁分 数 阶 模 型 比 整 数 阶 模 型 更 有 应 用 价 值援 近年来袁分 数 阶 微 分 方 程 引 起 了 人 们 极 大 的 关 注袁 除了其自身 理 论 的 深 入 发 展袁 它 在 分 子 物 理 学袁 流 体力学袁黏弹性力学袁电化学分析袁生物系统的电传 导等领域有广泛的应用咱员原苑暂 援 大部分学者致力于有 限区间上的分数阶边值问题咱愿原员员暂 袁对无穷区间上的

关于一类四阶非线性微分方程的正解

关于一类四阶非线性微分方程的正解
= s u p



, s ) ,
厂 :l i m i n f m i n丛
J 。 G , ) G ( , r ) d r 。
为 了方便 , 定义 :

广 :l i m s u p m a ) 【

K r = { u E K : I l u l l < r ) , = { ∈ K : I l u l = r ) 。
摘 要:通过使 用不动 点指数 定 理 ,在 适 当的 条件 下 ,给 出一 类四 阶非 线性微 分方程 的 一 个正解 的存 在
性 结果 。
关 键 词 :正解 ; 锥 ;不 动 点指 数
中 图分类号 : O1 7 5 . 8
文献标 识码: A
两端 简单支撑 的弯 曲弹性梁 的平衡 状态可用 四阶微 分方程 的两点边 值 问题 来 描述 , 关 于该 问题 解 的存 在 性 已有 很 多 作 者 研 究 过 并 获 得 一 些 存 在 性结果 】 。 下面讨论非线性四阶边值问题的一个正 解 的存 在性 :
m r I i n 3 1 u ( t )  ̄ o ' l l u l l > H , 于 是 V £ ∈ [ L 叶 1 叶 , } ] 』 ,
M J 。 G ( s , s ) f ( s , ) ) ,
r 1
引理 1 对于任意 t , s ∈ [ 0 , 1 】 , G , 以 下性质:
( 1 ) , S ) >O ;
( 2 ) , s )
, s ) ;
( 3 ) , s ) , f ) , s ) 。 方程( 1 ) 相 当于方程
J 。 J 。 c ( t , £ ) ( G ( , r ) ) c ( s , s , u ) ) d T :

分数阶微分方程数值解的一种逼近方法讲解

分数阶微分方程数值解的一种逼近方法讲解

分数阶微分方程数值解的一种逼近方法By:Pankaj Kumar, Om Prakash Agrawal摘要本文提出了一类分数阶微分方程(FDEs)的数值解方案.在这种方法中,FDEs 被Caputo型分数阶导数所表现. Caputo型分数阶导数的属性可以让一个分数阶微分方程减弱为一个Volterra型积分方程. 这样做了之后,许多研究Volterra 型积分方程的数值方法也可以应用于寻找FDEs的数值解. 本文总时间被划分为一组小区间,在两个连续区间中,用二次多项式逼近未知函数. 这些近似被替换成转化的Volterra型积分方程由此获得一组方程. 这些方程的解提供了FDE的解. 这种方法被应用于解决两种类型的FDEs,线性和非线性. 用这里给出的方法得到的解能与解析解和其他方法的数值解较好的吻合. 同时结果说明这种数值方法是稳定的.1.引言本文讨论分数阶微分方程的数值解. 分数阶导数和分数阶积分近年来收到了广泛的关注. 在许多实际应用中,分数阶导数和分数阶积分为考虑的系统提供了更加精确地模型. 比如,分数阶导数已经被成功地运用到模拟许多粘性材料的依赖频率的阻尼行为.1980年之前,Bagley 和Torvik提出了这个领域已经被研究的工作的一个回顾,并且说明了半阶导数模型可以非常好地描述阻尼材料的频率以来. 另一些学者说明了分数阶导数和分数阶积分在电化学过程,电解质极化,有色噪声,粘性材料和混沌领域的应用. Mainardi,Rossikhin和Shitikova 提出了分数阶导数和分数阶积分在一般固体力学,特定粘弹性阻尼模型中的应用的调查. Magin提出了分数阶微积分在生物工程的三个关键部分的回顾. 分数阶导数和分数阶积分在其他领域的应用以及相关的数学工具和技巧还可以在许多其他文献上找到.系统模型中分数阶导数的引进大多会导致分数阶微分方程的出现. 对某些特定的分数阶微分方程在通常系统条件下的解,已经有几种方法被找到. 这些方法包括,拉普拉斯变换,傅里叶变换,模态综合法和特征向量展开法,数值法以及基于Laguerre积分公式的方法. 然而,这些方法中大多数不能被应用到非线性分数阶微分方程. 更进一步的,正如Diethelm等人指出的,这些方法很多只能应用到特定类型的分数阶微分方程,并且人们并不知道他们能否被推广. 并且,在很多作者的研究成果中,并没有出现系统性的收敛性分析.最近,对于能被应用到线性和非线性分数阶微分方程的数值稳定数值逼近技巧,人们的兴趣愈发浓厚. 这些方法技巧大多利用了分数阶微分方程可以被减弱为Volterra型积分方程的特性. 因此,Volterra型积分方程的数值解法也可以应用到分数阶微分方程的解当中. Diethelm等人提出了分数阶微分方程数值解的一种PECE方法,其中P,C,E分别代表预测,校正和估计. 这样一来很多学者又推广了应用于常微分方程和分数阶微分方程的Adams–Bashforth–Moulton 型预测-校正格式. 这种方法的提出也是利用分数阶微分方程可以被转化为Volterra型积分方程的特点. 这些作者同时提出了误差分析和用Richardson外推法改善数值精度的延伸. Ford和Simpson提出了一种阶数大于1的分数阶微分方程的数值解法. 在该公式中,阶大于1的分数阶微分方程被减弱为阶小于1的分数阶微分方程,然后用相应的数值解法解由此导出的系统. 在所有这些方法当中,节点之间的未知函数用线性函数逼近. Kumar and Agrawal提出了阶数大于1的分数阶微分方程的数值解法. 这种方法要求就y(t)和它的导数在时间节点上连续.本文基于古典分数阶微分方程可以转化为Volterra型积分方程的特点也提出了一种数值方法来逼近分数阶微分方程的解. 特别地,我们用二次逼近函数来建立这种算法,结果说明这种方法可以被应用到寻求分数阶微分方程的数值解. 我们还通过两个例子,线性和非线性问题的解决,说明了这种方法的高效和准确,并且这种数值方法是稳定的.2.数值算法关于分数阶导数的定义已经出现有好几种,它们包括Riemann–Liouville, Grun-wald–Letnikov, Weyl, Caputo, Marchaud,和Riesz分数阶导数. 这里,我们规定使用Caputo导数.其中,Caputo导数的定义是, (n-1<α<n),(1)其中,α是导数的阶数,n是比α大的最小的整数.式(1)早在19世纪就在Liouville的论文中被提出,在Caputo的论文发表前一年它被Rabotnov所用. 然而,在文献中,被(1)式所定义的分数阶导数作为Caputo导数被广泛认知.在接下来的讨论中,我们考虑含有Caputo导数的初值问题:(2)在初始条件:, k=0,1,...,n-1,(3)下的解,其中,f是任意函数,是y的k阶导数,,k=0,1,…,n-1是指定初值条件. 假设这个函数关于参数和积分区间都是连续的,并且对于它的第二个参数满足Lipschitz条件.在纯数学中,Riemann–Liouville导数比Caputo导数应用更加广泛. 然而,这里考虑Caputo导数是因为以Riemann–Liouville导数为基础的分数阶微分方程要求y(t)在t=0点的导数和积分为0.一般来说,这些条件的物理意义不是已知的,并且在实际应用中,他们是不可用的. Lorenzo and Hartley讨论了寻找在更一般的情况满足下初始条件的正确格式的问题. 在Diethelm and Ford的文章中,方程(2)和(3)被证明可以等价描述为:,(4)其中g(t)为:. (5)为了解释以二次多项式为基础的数值方法,我们假设我要求的是由(2)式定义的分数阶微分方程从0到T的积分. 为了达到这个目的,我们把时间T等分成N 份,令h=T/N,作为时间区间的每一个部分. 时间在网格点上被表示为. 同时假设y(t)的数值逼近值被网格点所决定. 该方法的基本思想是在相邻的两个时间节点和上数值地获取函数y(t)的值,然后重复这个过程来接近所求积分直到取到终点T.为了便于接下来的讨论,我们规定如下记号:, ,这里的方法需要对方程(4)每一步求两次积分值. 这里有两种方法来达成目的.第一种,用一些近似函数逼近y(t),然后用一种数值方法确定式(4)的积分值.这里需要在未知积分的情况下对和作初始的估计. 第二种, 都用近似函数来显式地逼近y(t)和f(t,y(t))以及确定式(4)的积分. 注意在这种情况下,和会作为参数出现在函数f当中. 本文利用的是第二种逼近方法.现在,我们给出算法的详细思路. 首先我们需要确定y(t),,的值. 用二次插值函数可以在区间[0,]上逼近y(t)和f(t,y(t)):(6)以及(7)其中,是函数在第k个时间节点的值,,k=0,1,2是二次插值函数,其中下标(j,k)代表在第j+1,j=0,…,m步的第k个近似函数.我们首先确定y(t)在和处的值. 把(7)式带入(4)式,并积分得到:(8)其中,,(9)可以精确计算得到. 注意式(8)需要知道f在和的值(或者间接地说,y的值). 为了得到,在[,]上把f(t,y(t))近似为:,(10)其中,,是另一个二次插值函数. 函数,k=1,2,3由下式给出,函数由相似的办法定义.把(10)代到(4)中,积分得到:, (11)其中,, (12)可以(9)中一样被精确计算出来. 由(7)可以得到的值为:(13) 这里,我们充分利用了二次多项式的性质. 在非二次多项式的情况下,将会有不同的参数.把(13)代到(11)得到,+, (14)注意到(8)和(14)是关于两个未知量和的方程,可以用Newton–Raphson法,不动点迭代或者其他非线性方法求解. 这里,我们用Newton–Raphson法求解这些非线性方程. 这个方法需要对和作一个初始的估计. 当α大于1, 由t=0处的斜率可以得到关于和的更好的估计. 然而,在这里,对于α>1和α<1我们对把这些变量的初值估计为. 注意在每一次迭代式,时间步长取2h.现在我们假定在处,y的值是已知的,我们要求的是和处的值. 根据以上的逼近方法,和可以被表示为:+, (15)++(16)其中,,k=2m,2m+2,2m+2,,k=2m,2m+1/2,2m+1是和,用同样方法确定的系数. 注意(15),(16)的积分可以被数值确定因为y(t)在,处的值是已知的. 这些方程含有两个未知量和,而他们可以通过Newton–Raphson法得到. 本文中,我们把作为和的初值估计. 这样一来,方程(1)就可以在需要的区间上求积分.作为特殊情况,考虑如下非线性系统:.这种条件下,,式(16)和(15)减弱为:, (17) 其中=, (18)=, (19)=-, (20)=1+(21)=-, (22)=- . (23)(17)是一组线性联立方程,可以用线性方法求解.请注意以下两个补充说明. 第一,方程(1)只含有一个y(t). 当y(t)是一个向量函数时,算法同样成立.不过,所有的y(t)和f(t,y(t))必须换成相应的近似向量函数. 第二,算法需要保存所有算过的的y的值. 这是很多分数阶微分方程的典型特征. 这将会导致一些问题,特别是当y的维数和分数阶有限元系统一样大时. 这种情况下,系统有临近的短期记忆,y(t)在过去一段时间长度的值可以忽略不计,以此来改善对存储空间的需求和计算效率.3. 数值结果为了说明这种方法的效率,我们分别考虑一个线性和一个非线性的算例. 讨论这些例子是因为他们解的逼近格式是可靠的,并且可以用其他数值方法求解. 这样我们就可以把用这种方法得到的结果和解析解以及其他数值方法的结果相比较.3.1例1线性方程在第一个例子中,我们考虑如下给出的线性方程:,0<α<2,(24) 且. (25) 初始条件仅当α>1是成立. (24),(25)的解是:,(26) 其中,. (27) 是Mittag–Leffler函数的阶.图1.α=0.75时y(t)的比较(O:解析解,X:本文数值方法)表1.α=0.75时h在不同值下函数y(t)的误差对不同的α和h可以得到很多结果,这里给出其中一些. 在各种情况下,我们另T=6.4s.考虑这个区间是因为它接近α=2的系统的时间. 这里图(1)比较了α=0.75时的解析解和二次数值方法. 在这种情况下,我们令h=0.1s. 注意到这两个结果几乎完全重合. 为了强调收敛性,表(1)列出了当α=0.75,h分别等于0.4,0.2,0.1,0.05和0.025时的结果. 注意到随着步长的缩小,误差也如期望一样的缩小了. 在大部分节点误差比R=E(2h)/E(h)都非常接近3.37,这表明误差阶为1.75(或E(h)=O()).图2. α不同时y(t)的比较(O:α=0.25,X:α=0.5,+:α=0.75,Δ:α=0.95,*:α=1.)图3.α=1.5时y(t)的比较(O:解析解,X:本文数值方法)表2.α=1.5时h在不同值下函数y(t)的误差图(2)展示了h=0.025,α分别等于0.25,0.5,0.75,0.95和1时y(t)的数值结果.因为解析解和数值结果完全一致,因此图中没用画出解析解. 当α=1时,精确解为y(t)=. 注意到随着α越来越接近1,数值解越收敛到解析解y(t)=,即在极限情况下,分数阶微分方程的解接近整数阶导数的解.更进一步地,我们给出了α>1的一系列结果.α<1和α>1的结果是分离的,因为y(t)的斜率在α=1处有一个跃迁.图3比较了y(t)在α=1.5,h=0.4时的解析解和数值解. 两个结果又一次几乎完全重合.为了突出收敛性,表2给出了α=1.5,h分别等于0.4,0.2,0.1,0.05和0.025时的数值解. 正如之前观察到的一样,在这种情况下,随着步长的减小,误差也随之减小. 在这种情况下,误差比接近5.5,这表明误差阶为2.5(或E(h)=O()). 这样一来,通过观察α<1和α>1的收敛结果,可以知道误差的收敛阶为1+α(或E(h)=O()),即误差的阶不仅依赖于h,还依赖于导数的阶α.图4. α不同时y(t)的比较(O:α=1.25,X:α=1.5,+:α=1.75,Δ:α=1.95,*:α=2.)图5.α=0.25,0.75,1.25,1.75时y(t)的比较(O:解析解,X:本文数值方法)表3.本文数值方法和文献(35)中y(t)的误差的比较.图4展示了h=0.025,α分别等于1.25,1.5,1.75,1.95和2时y(t)的数值结果.又一次,因为解析解和数值结果完全一致,因此图中没用画出解析解. 当α=2时,精确解为y(t)=cos(t). 注意到随着α越来越接近2,数值解逐渐收敛到整数阶导数的解. 图2和图4展示的收敛结果十分重要,因为他们说明了在极限情况下,分数阶微分方程和他们的解逼近整数阶微分方程以及他们的解析解. 表3比较了t=1.0文献(35)的解的误差和用本文数值方法在α=0.5和1.25,h=0.1,0.05,0.025的解的误差. 注意到本文的方法得到了更低阶的误差. 这是因为,这里的方法是一种高阶方法. 当α和h取其他值时这种趋势也能显现出来.3.2 例2.非线性方程在第二个例子中,我们考虑一个如下定义的非线性方程:(28) 且. (29) 和之前一样,第二个初始条件仅适用于α>1. (28)(29)的精确解在文献(35)中已给出,(30)注意到当α<1, 方程的解在t=0处的斜率趋近于无穷. 因此,他可能导致在t=0附近出现一个巨大的数值误差.表4. α=0.75和1.5,h取不同值下y(t) 的误差.表5. 文献(35)中y(t)的误差和用本文数值方法得到的y(t)的误差的比较上面给出了一些在不同α和h下的数值结果. 图5表示了h=0.05,α分别等于0.25,0.75,1.25,1.75时解析解和数值解的结果. 由它可以得到(1).解析解和数值解基本重合,当α取其他值是,可以得到同样的结果. (2)尽管在t=0处斜率非常大,但是方法给出了非常精确的结果. (3)正如预期的那样, 在t=1处,对所有的α,y(t)的值收敛到0.25. 表4列出了当α=0.75和1.5,h=0.1,0.05和0.025的数值解的误差. 注意到误差随着步长的减小而减小. 同样的趋势在α取其他值时也能观察到. 在尝试过的α的值中,误差比R=E(2h)/E(h)表明没有特定的收敛阶. 然而,当α=0.75和1.5时,收敛的误差的平均值接近12和15.表5比较了文献(35)中在t=1.0处解的误差和用这种方法在α=0.25和1.25,h=0.05时的误差. 这里我们用的是文献(35)中用Richardson外推法得到的值. 观察得到,本文的方法又一次给出了更小的误差. 当α=0.25时,这种方法给出了比Richardson外推法小得多的误差. 这可能是因为,当α等于0.25时y(t)在t=0附近的斜率改变非常迅速,并且线性方法不能精确地获得结果. 需要指出的是,这种数值方法对于α和h取其他值时同样给出了更加精确的结果.4.结论本文给出了一种经典的分数阶微分方程的数值逼近方法. 这里的分数阶微分方程是依据Caputo型分数阶导数给出的, 这种导数的性质可以把分数阶微分方程减弱为Volterra型积分方程. 时间区间被分成一组网格,通过3个连续节点的二次插值多项式逼近未知的和已知的函数y(t)和f(t,y(t)). 把这些多项式带入Volterra型积分方程可以得到一组代数方程,这种数值方法的提出就是用来解这些方程以及获取y(t)的解. 通过一个线性一个非线性的例子的解决,说明了这种数值方法的作用. 用这种方法得到结果和解析解以及其他数值方法的结果是一致的. 还得到一个结论就是结果随着步长的减小而收敛. 在极限情况小,当α逼近整数值,数值方法会得到一个整数阶系统的解. 结果还表明这种方法是数值稳定的.。

一类Caputo分数阶微分方程边值问题多解的存在性

一类Caputo分数阶微分方程边值问题多解的存在性郭彩霞;任玉岗;郭建敏【摘要】We investigate the existence and multiplicity of positive solutions for nonlinear Caputo fractional differential equation boundary value problem Dα0+u(t)+f(t,u(t))=0,t∈(0,1),u′(0)=u(1)=0{ , Where 1<α≤2,f:[0,+∞)×→[0,+∞)is continuous,and Dα0+is the standard Caputo differentiation.In the process of proof,we first transform it into integral equation,then differ-ential equation boundary value problem is further converted to discuss the problem of integral operator fixed point.Finally,by means of Leggett-Williams fixed point theorems on cone,ex-istence results of at least three positive solutions are obtained.The properties of the Green function and the conditions of the nonlinear term is very important.%研究一类Caputo分数阶微分方程边值问题:Dα0+u(t)+f(t,u(t))=0,t∈(0,1),u′(0)=u(1)=0{,【期刊名称】《广西科学》【年(卷),期】2016(023)004【总页数】4页(P374-377)【关键词】分数阶微分方程;边值问题;Leggett-Williams不动点定理【作者】郭彩霞;任玉岗;郭建敏【作者单位】山西大同大学数学与计算机科学学院,山西大同 037009;山西大同大学数学与计算机科学学院,山西大同 037009;山西大同大学数学与计算机科学学院,山西大同 037009【正文语种】中文【中图分类】O175.8分数阶微分方程在工程、化学、物理、生物等领域有着广泛应用,例如热传导领域和流体学领域[1-3],而且分数阶导数模型克服了经典整数阶微分模型理论与实验结果不吻合的缺点[4],因此研究分数阶微分方程边值问题有着重要的意义.近年来,大量文献报道微分方程[5-6]和分数阶微分方程[4,7-10]边值问题解的存在性.2005年,当1<α≤2时,Bai等[7]推导了分数阶微分方程边值问题目前研究分数阶微分方程边值问题的主要工具有锥拉伸与锥压缩不动点原理、Krasnoselskii不动点原理、Schauder不动点原理上下解等.本文利用Leggett-Williams不动点定理,参照文献[9]中的方法研究Caputo分数阶微分方程边值问题定义1.1[11] 一个连续函数u:(0,+)→的α阶Caputo导数定义为引理1.1[11] 令α>0,若u∈ACn[0,1]或u∈Cn[0,1],则引理1.2 令α∈(1,2],给定h∈C[0,1],则u′(0)=u(1)=0,证明由引理1.1可得,(1.1)式等价于方程,其中C1,C2∈.从而.由(1.2)式可知s.引理1.3 引理1.2中的G(t,s)有下列性质:(i)G(t,s)∈C([0,1]×[0,1],)且G(t,s)>0,t,s∈(0,1);证明 (i)~(iii)显然可得,只需证明(iv).令,则令γ,β,θ是锥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‖.(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.β(u1)<d,a<α(u2),d<β(u3),α(u3)<a.令E=C[0,1],其范数为‖u‖|.当时,定义E中锥P为P={u∈E:u在[0,1]上是非负的凹函数,‖u‖}.又定义锥P上的非负连续凹函数α,ψ和非负连续凸函数γ,β,θ为α(u)≤β(u),定理2.1 假设存在非负实数a,b和c使得,若f满足下列的条件:(H3)f(t,u(t))≤c αΓ(α),t∈[0,t3]∪[1-t3,1],u(t)∈[0,c].证明在锥P上定义算子A为首先,对任意u∈P,由(2.1)式和(2.2)式可知α(u)≤β(u),‖u‖≤γ(u). 若,则‖u‖≤c.又由(H3)得(1)若u∈Q(γ,β,a,c)且,则(2)若,由(H1)得(3)若u∈Q(γ,α,b,c),且,则(4)若,由(H2)得α(u1)>b,β(u2)<a,α(u3)<b,β(u3)>a.本文研究了一类Caputo分数阶微分方程边值问题多解的存在性.证明时,将微分方程边值问题转化为积分方程,进一步转化为讨论积分算子不动点的问题,然后通过运用Leggett-Williams不动点定理该分数阶微分方程边值问题至少有3个正解存在的结果,其中格林函数的性质和非线性项的条件至关重要.。

一类分数阶非线性系统解的性质

一类分数阶非线性系统解的性质田永强;钟守铭;包姣【摘要】给出了Caputo分数阶非线性微分方程的解的性质,并在此基础上给出一类含有边界条件的分数阶非线性系统的解.%Some properties on the concept of solution for fractional order differential equation and establish the solution for nonlinear fractional order system with boundary value problem are shown.【期刊名称】《科学技术与工程》【年(卷),期】2011(011)012【总页数】4页(P2633-2635,2639)【关键词】分数阶;Caputo分数阶导数;非线性系统【作者】田永强;钟守铭;包姣【作者单位】电子科技大学数学与科学学院,成都,611731;电子科技大学数学与科学学院,成都,611731;电子科技大学数学与科学学院,成都,611731【正文语种】中文【中图分类】O175.14分数阶微积分[1—3]自从17世纪在大数学家Leibniz,Euler,Lagrange,Able,Liouville 等人的推动下得到广泛的发展,并建立了分数阶微积分的理论基础。

近年来,在分数阶的发展中,fractional calcul已经代表着任意阶的积分和微分。

分数阶系统由于其包括整数阶系统,其实质是任意阶的系统,广泛地应用于图像处理、智能控制、机器人等各个分支,其地位的重要性也越来越引起人们的重视。

为了解决分数阶微分方程解的问题,最常见研究方法就是Laplace,Fourier变换及Taylor展式,数值方法等[1]。

关于分数阶非线性系统的研究,一是传统的非线性整数阶系统采用分数阶控制策略来研究;二是对非线性分数阶系统进行研究。

文献[4]采用不动点理论研究了分数阶带有线性边界条件的非线性系统解的性质文献[5]讨论了基于四点的一些非线性系统边界解的情况,给出了系统的解的表达式。

两点分数阶微分方程耦合系统边值问题的解

两点分数阶微分方程耦合系统边值问题的解曹竞文;胡卫敏【期刊名称】《江汉大学学报(自然科学版)》【年(卷),期】2014(42)3【摘要】讨论一类非线性分数阶微分方程耦合系统的两点边值问题,应用Green 函数将微分系统转化为等价的积分系统,应用不动点定理证明系统正解的存在性和唯一性,并给出系统无解的充分条件。

%Discusses a class of the two-point boundary value problems of nonlinear fractional differ⁃ential equation of coupling system,using the Green function,differential system can be converted to equivalent integral system,with the fixed point theorem ,the existence and uniqueness of positive solutions for system are abtained,sufficient conditions of no solutions are given.【总页数】4页(P23-26)【作者】曹竞文;胡卫敏【作者单位】伊犁师范学院数学与统计学院,新疆伊宁 835000;伊犁师范学院数学与统计学院,新疆伊宁 835000【正文语种】中文【中图分类】O175.8【相关文献】1.一类分数阶微分方程耦合系统边值问题解的存在性 [J], 薛益民;苏有慧;刘洁;苏莹2.非线性分数阶微分方程奇异两点边值问题的解 [J], 韩仁基;蒋威3.一类非线性分数阶微分方程耦合系统两点边值问题的可解性 [J], 曹竞文;胡卫敏4.分数阶微分方程耦合系统多点积分边值问题的解 [J], 张宁;张娣;史小艺5.一类非线性分数阶微分方程耦合系统边值问题的两个正解 [J], 彭钟琪;李媛;薛益民因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第55卷第4期 2017年7月 吉林大学学报(理学版) 

Journal of Jilin University(Science Edition) V01.55 No.4 

July 2017 

doi:10.13413/j.cnki.jdxblxb.2017.04.15 

一类非线性Caputo型分数阶 微分方程耦合系统的正解 

薛益民,苏莹,苏有慧 (徐州工程学院数学与物理科学学院,江苏徐州221018) 

摘要:利用Schauder不动点定理、Leray—Schauder抉择理论和Banach不动点定理,研究一 类含积分边值条件的非线性Caputo型分数阶微分方程耦合系统边值问题,得到了该耦合系 统正解存在性和唯一性的充分条件,并举例说明定理的适用性. 关键词:分数阶微分方程;耦合系统;正解;不动点定理 中图分类号:O175.8 文献标志码:A 文章编号:1671—5489【2017)04—0853—08 

Positive Solutions of a Coupled System of Nonlinear Caputo Type Fractional Differential Equations 

XUE Yimin,SU Ying,SU Youhui (School of Mathematics and Physical Science,Xuzhou Institute of Technology,Xuzhou 221018,Jiangsu Province,China) 

Abstract:We studied the boundary value problem of a coupled system of a class of nonlinear Caputo type fractional differential equations with integral boundary value conditions by applying Schauder fixed point theorem,Leray—Schauder choice theory and Banach fixed point theorem.We obtained sufficient conditions for the existence and uniqueness of positive solutions of the coupled system,and illustrated the applicability of the theorem by examples. Key words:fractional differential equation;coupled system;positive solution;fixed point theorem 

0 引 目 

’一2 < a <3

,’

。 

收稿日期:2016—08—29. 作者简介:薛益民(1977一),男,汉族,硕士,讲师,从事微分方程理论及应用的研究,E-mail:xueym@xzit.edu.cn. 基金项目:国家自然科学基金(批准号:11301454)、国家自然科学基金数学天元基金(批准号:11526177)、江苏省自然科学基金 (批准号:BK20151160)、江苏省高校自然科学基金(批准号:14KJB110025)、江苏省六大人才高峰项目(批准号:2013一JY-O03)和 安徽省教育厅重点项目(批准号:KJ2017A839). 854 吉林大学学报(理学版) 第55卷 其中: D。表示Caputo型分数阶导数;fe c([0,1]×[O,。。),[0,cx3)).受文献[6]启发,本文利用 Schauder不动点定理、Leray—Schauder抉择理论和Banach不动点定理研究如下非线性Caputo型分数 阶微分方程耦合系统正解的存在性和唯一性: D。U(£)+f(t, (£), (£))一o, o<t<1, D (£)+g(t, ( ), ( ))一0,0<t<1, “(0)一 (0):…一“( 2 (o)=== ( (o)一0, (1)一 J “( )d , (1) J 0 r] (0)一 (0)=…一 ‘ (0)一 ‘ (o)一0, (1):itI ( )ds, 

其中:n<d, n+1, ≥2,nE ,0<A<n;f,gE C([0,1]× × , ); D 表示Caputo型分数 阶导数. 

1 预备知识 定义1 Ll3 函数,:[0,Cx3)一 的a阶Caputo型导数定义为 D 厂( )一 南j。(£一 f ( , 其中n>O, :[a]+1,[a]为a的整数部分. 定义2E¨ 妇 函数 :[0,∞)一 的口阶Riemann—Liouville积分定义为 

( )一志j。(£一s) 厂(s) 其中:a>O;右端积分存在. 引理l 若a>。,则齐次分数阶微分方程 。。“(£)一。有唯一解 ( )一 . 引理2[¨ 若a>0,则Riemann—Liouville积分和Caputo型导数有如下性质: ㈤ 凳 

引理3E吲 V ( )∈c[o,1], <y< +1,y— ,p, ≥2, ∈ ,O<a<n,分数阶微分方程边值 问题 f D (£)+ ( )一0,0<t<1, l (0)一 (o)一…一“‘ (o)一U (0)一0, (1)= l (s)ds, 2 r1 ( ) 

有唯一解 ( )一f 1G (£,s) (s)ds,其中 

G (£,s)一 二 二 !Z二 ± 2二 丝二 2 二 2二 ( 一a)F(r+1) 

年 , o≤£≤ ≤1. — )F(y+1) ’ 、‘ ‘、~ 

,0≤S≤t≤1, (3) 

引理4(Schauder不动点定理)_16 设己,是Banach空间X的一个非空有界闭凸子集,P是U到其 自身的全连续映射,则至少存在一个aT.∈U,使得Px— . 引理5(Leray-Schauder抉择)[16 设V是Banach空间y的一个有界闭凸子集,E是 中相对开 球,且0∈E.若算子F: —V全连续,则下列叙述之一成立: 1)F在E中至少存在一个不动点; 2)存在U∈aE和∞∈(O,1),使得 一 F(M). 引理6(Banach压缩映射原ti)E" 设D是Banach空间E的闭子集,F:D—D是一个严格压缩 映射,即v z,yE D,l Fx--Fy I≤ l z— l成立,其中0<k<l,则F在E上有唯一不动点. 第4期 薛益民,等:一类非线性Caputo型分数阶微分方程耦合系统的正解855 2 主要结果 令X一{“( )l甜(£)E C(Eo,1], )),V uE X,定义范数l1 ll一 ]I (£)I,则(x,ll・I1)是一 个Banach空间.令Y一{ ( )I (£)E C(EO,1],哝)},V E Y,定义范数l1 l—s p—f ( )f. V( , )∈X×y,定义范数lJ( , )lI—ll Il+Il Il,则(x×Y,ll( , )l】)也是一个Banach空间. V( ,口)∈x×y,tE[O,1],定义算子T:X×y—x×y为 T(“, )( )一(T1( , )(£),T2(“, )( ))一 

(j 0 G0 )/( ( 口( j: ), ( ds). (4) 

则T的不动点即为耦合系统(1)的解. 为方便,记 K1一 丝 ±丛二 ( ~ ) r(d+1)’ Kz一 . ㈤ 

引理7设f,gEC(E0,1]× × , ),则式(4)定义的算子T:X×y—x×y是全连续的. 证明:由Gu(t,s), (t,s),f,g的连续性知,算子T是连续的,且T(X×y)c X×Y,即 T:X×y—X×Y.首先证明算子T一致有界.令ncX×Y有界,由f,g的连续性知,V( , )E 0, tE L0,1],存在L ,L2>O,使得 I f(t,“(£), (£))1≤L , l g(t,“(£),口( ))l≤L . 因此,V( , )Ef2,t∈EO,1],有 1 111(“, )( )I≤-f I (s, (s), (s))l d +Jl :: l厂(s, (s), (s))l ds+ 

l )lds+ l , )Ids≤ L ( + + ・ + ・F )一L。K . 类似可得 I Tz c… f≤Lz( + + ・ + 。 )一Lz K2. 

因此,算子T是一致有界的. 下面证明算子T:x×y—x×y等度连续.令0≤£ <£ ≤1,由式(3),(4)可得 

l T1(“, )(t2)一T1( , )(t1)l—l l(( (£z,5)一G (£1,s))f(s,“(s), (5))ds l≤ 

1( ( 一 r )+( 一 )).r(口+)\ (咒一 )口 … ‘ 。、。 /‘ 类似可得 l Tz( , )( z)一Tz( , )(£ )1≤ 鲁 ( (£ 一 )+( l一£ )). 

因此当tz一£ 时,有 Il丁l(“, )( 2)一Tl( , )( 1)II一0, Il ( , )( 2)一 ( , )( 1)ll一0. 从而当t2一 时,有 Il T(u, )( 2)一T(u, )( 1)Il—lI T1(“, )( 2)一T1( , )( 1)lI+ lJ T2( , )( )一 ( ,口)(z )ll一0. 所以,算子T:X×y—x×y是等度连续和一致有界的,由Arzela—Ascoli定理知算子T是全连续的. 定理1 设f,gE C([-0,1]× × , ),若存在常数a ,bi>0, 0,1,2,使得下列不等式成立: 1)J f(t,“, )j≤n。+日1 l I +n2 J l ,呀1,rl2 E(0,1);

相关文档
最新文档