一类二阶非线性微分方程奇异边值问题唯一整体解的精确渐近行为

合集下载

一类二阶非线性奇异边值问题的正解

一类二阶非线性奇异边值问题的正解

一类二阶非线性奇异边值问题的正解近十几年以来,研究非线性奇异边值问题的学术界越来越受到重视,非线性奇异边值问题具有重要的理论意义和现实意义,是一个值得深入研究的热点课题。

非线性奇异边值问题是指在一定条件下,当某一函数取极值时,至少有一个变量处于奇异边界(即满足一种奇异边界条件,如符号条件约束等)的优化问题,属于多变量函数的不等式约束优化问题。

非线性奇异边值问题存在着严格的数学限制,并且本质上属于非线性可微易变分析问题,其解的精确性和一致性极其重要。

在研究的过程中,学者们发现,有效研究非线性奇异边值问题的关键在于正确性,即给出的解满足所有约束条件,是真实有效的解。

不久前,学者们在数学模型分析中取得了重要进展,发展出了一类可以有效解决非线性奇异边值问题的新方法正解法。

正解法的核心思想是在优化的同时进行函数的精确解析,利用函数的数学特性引入变量的约束,从而有效解决非线性奇异边值问题,对于研究者来说是一个重要技术手段,也是一个值得挖掘学习的话题。

首先,非线性奇异边值问题的正解法是以建模和求解为主要内容,其目标是获得满足所有约束条件的正解,而正解可以理解为优化过程中函数满足偏导数等于零的可行解。

换句话说,正解就是寻找满足函数所有约束条件的可行解的过程,这是研究非线性奇异边值问题的基础。

其次,正解法的实质是运用数学方法解决非线性奇异边值问题,解决过程主要包括解析法和迭代法。

解析法是采用友元求解和极值分析技术,结合函数的数学特性,建立满足约束条件的数学模型,寻求最优解;迭代法是采用熵增准则,利用熵增的改进方向,进行非线性奇异边值的综合优化。

解决非线性奇异边值问题需要考虑到的因素不仅包括函数的平滑程度、优化后的可行性等,还需要考虑优化方案的可用性,因为正解法在解决非线性奇异边值问题时,它不仅考虑优化方案的求解,还要求其可用性。

可用性指的是优化方案是否能够满足实际的要求,例如满足所有约束条件,和各种物理特性的限制,这是非线性奇异边值问题的解决必备步骤。

一类二阶奇异微分方程边值问题正解的存在性

一类二阶奇异微分方程边值问题正解的存在性
M a . 2 11 r 0
文 章 编 号 : 0 8 0 7 ( 0 1 0 — 0 3 0 1 0 — 1 1 2 1 ) 20 2 — 8

类 二 阶 奇 异 微 分 方 程 边 值 问 题 正 解 的存 在 性 苏 源自恒 , 志 明 , 郭 白定 勇
( 州 大 学 数 学 与信 息科 学 学 院 , 东 广 州 5 0 0 ) 广 广 1 0 6
正 解 。 结 果 推 广 了现 有 文 献 中 的相 关 结 论 。 该
关 键 词 : 阶 奇 异 微 分 方 程 ; 解 ; 值 问题 ; 界 点 理 论 二 正 边 临 中 图 分 类 号 : 7 . O1 5 1 文 献 标 志码 : A
1 引 言 及 主 要 结 果
在微分 方 程研究 的诸 多方 向中 , 值 问题 正解 的存在 性具 有广泛 的应用背 景 和重要 的理论 意义 。 边 例
+ 。 。1 8 。 9 1年 , eetc i L 运用 打靶 法 在 边 值条 件 ( ) 0和 甜 M ) B rsyk 等 4 0= ( 一0下 , 研究 二 阶微 分 方 程 +( N一1 I=f u 正解 的存 在性 , 中 , ) t () 其 Ⅳ> 1 正整数 。2 0 为 0 5年 , o h u e等[改 进 了文 献 E ] B ner 5 4 的 方 法 , 与 文 献 E - 同的边 值 条 件 下 , 正 整 数 N> 1 推 广 到任 意 实 数 正 0 研 究 了二 阶微 分 方 程 在 4相 ] 将 , > ,
() 1
() 2
其 中 , E( , × ,() 一个 连续 函数 , 0 一0 故 方程 () £ M O +c ) c£是 3 P( ) , 1 在 一0处具 有奇 异性 。 在本 文 中 , 如下 的基本 假设 : 作

一类二阶奇异微分方程的迭代解

一类二阶奇异微分方程的迭代解

一类二阶奇异微分方程的迭代解一类二阶奇异微分方程的迭代解是指通过对方程逐步进行递归计
算得到的解。

这种方法通常适用于形如$y''=f(x)y$的方程,其中
$f(x)$在一些特殊点上为无穷或为零。

迭代解的得到需要先对方程进
行一些变形,通常需要借助于某些特殊函数或变量的引入。

通过对方
程进行迭代,每一次迭代都会得到一个更加精确的解,直到达到所需
的精度为止。

通常情况下,迭代解可以用级数形式表示,其中每一项
都可以通过递推关系得到。

迭代解的优点在于可以得到非常高的精度,并且在计算过程中可以控制误差的大小。

此外,迭代解方法还可以应
用于一些其他的特殊微分方程的求解中。

一类二阶非线性方程奇异边值问题

一类二阶非线性方程奇异边值问题

一类二阶非线性方程奇异边值问题
杨作东;毕卫萍
【期刊名称】《应用数学》
【年(卷),期】1997(10)1
【摘要】本文证明了方程div(|Du|p-2Du)+f(r.u(r),u'(r))=0(R1<r<R0)正对称解的存在性.这里f允许在u=0或。

u’=0处奇异.
【总页数】6页(P51-56)
【关键词】奇异边值问题;非线性方程;正对称解;边值问题
【作者】杨作东;毕卫萍
【作者单位】河南师范大学数学系
【正文语种】中文
【中图分类】O175.8
【相关文献】
1.一类奇异非线性二阶边值问题解的存在性 [J], 沈文国
2.一类二阶非线性奇异边值问题的正解 [J], 张艳红
3.一类奇异的二阶非线性常微分方程边值问题的正解 [J], 姜杰
4.一类二阶非线性奇异边值问题的多重正解 [J], 杨凤勉;张圣文
5.一类奇异二阶非线性方程两点边值问题正解存在性 [J], 杨作东
因版权原因,仅展示原文概要,查看原文内容请购买。

两类奇异的非线性微分方程边值问题解的存在性分析

两类奇异的非线性微分方程边值问题解的存在性分析

摘要分数阶微积分已有很长的历史. 早在1695年,在Leibniz 和L’Hospital 的往来书信中就已经提到了分数阶微分的概念. 在近三个世纪内,人们对分数阶微积分理论的研究主要集中在数学的纯理论领域. 然而在最近几十年内,许多学者纷纷指出分数阶微积分非常适合于刻画具有记忆和遗传性质的材料和过程,这些性质在经典模型中是常常被忽视的.如今,分数阶微分方程模型越来越多地被用于描述声学、热学系统、材料力学、信号处理、系统辨识、控制理论、机器人科学以及其它应用领域中的问题.本文的工作如下:第一部分是绪论,主要简要介绍了分数阶微积分和分数阶微分方程的研究历史和发展现状,以及分数阶微分方程正解存在性方面的研究工作.第二部分研究了一类奇异的非线性semipositone Sturm-Liouville 边值问题正解的存在性. 我们的主要方法是对非线性部分()f y 进行重新定义,使其转化成非奇异的p ositone 边值问题, 然后应用锥上的不动点定理以及泛函分析的知识证明该奇异非线性s emipositone Sturm-Liouville 边值问题的正解的存在性.第三部分讨论了一类奇异的非线性分数阶微分方程Dirichlet 边值问题正解的存在性. 我们的主要思想是重新研究非线性部分0(,(),())f t x t D x t β+,使其转化为非奇异的分数阶微分方程边值问题,然后再对每一个重新定义的非线性部分为0(,(),())n f t x t D x t β+(N n ∈)的边值问题,证明其存在正解n x ,最后通过紧集上函数列极限的性质给出原奇异非线性分数阶微分方程Dirichlet 边值问题的正解的存在性.关键词:微分方程,分数阶微分方程,边值问题,正解,奇异性,不动点.AbstractFractional calculus has a long history. As early as in 1695, the concept of fractional differential was already mentioned in the correspondence of Leibniz and L'Hospital. During the past three centuries, the research of fractional calculus theory was mainly concentrated in the pure theoretical field of mathematics. However, in the recent several decades many scholars in succession pointed out that fractional calculus is very suitable to characterize materials and processes with memory and hereditary properties, which were often neglected in the classical models .Nowadays, fractional differential equation models are increasingly used to describe the problems in acoustics, thermal systems, material mechanics ,signal processing, system identification, control theory, robotics and other applied fields .This thesis is divided as follows:The first part is an introduction, briefly presents the research history and development status of the fractional calculus and fractional differential equations, and some past research works about the existence of positive solutions of the fractional differential equations.The second part studies a singular nonlinear semipositone Sturm-Liouville boundary value problem. We redefine the nonlinear part ()f y , and make the singular boundary value problem transform into a nonsingular positone boundary value problem, and then prove the existence of a positive solution for the original singular nonlinear boundary value problem by using the cone fixed point theorem as well as knowledge of functional analysis.The third part discusses the positive solution existence for Dirichlet boundary value problem of a singular nonlinear fractional differential equation. We study itsnonlinear part 0(,(),())f t x t D x t β+, and have it transform into a nonsingular boundaryvalue problem, and then prove the existence of a positive solution n x for eachboundary value problem with redefined nonlinear part 0(,(),())n f t x t D x t β+(N n ∈), andfinally we give the existence of a positive solution for the original Dirichlet boundary value problem via the limit properties of a sequence of functions on compact sets. Keywords: Differential equation, fractional differential equation, boundary value problem, positive solution, singularity, fixed point.目录摘要 (1)Abstract........................................................................................................I I 第一章绪论 (1)1.1分数阶微积分的历史 (1)1.2分数阶微分方程的研究现状 (2)第二章带有奇异的非线性Semipositone Sturm-Liouville边值问题解的存在性52.1 引言 (5)2.2 预备知识 (6)2.3 主要结果 (7)第三章带有奇异的非线性分数阶微分方程Dirichlet边值问题正解的存在性 173.1 引言 (17)3.2 预备知识 (19)3.3 主要结果 (29)参考文献 (31)攻读硕士期间发表的论文 (34)后记 (35)第一章 绪论1.1分数阶微积分的历史牛顿和莱布尼茨发明的微积分是现代数学和古典数学的分水岭,数学的发展和应用自此发生了根本性的变化,分析、几何和代数一同成为数学的三个基本研究方向和工具.对大多数研究人员和工程师而言,分数阶微积分也许还是一个新奇的概念和数学工具,但它实际上早在300多年前就已被提出,和传统的整数阶微积分有着一样久远的历史.莱布尼茨最先引入/n n d y dx 来表示导数,也正是因为这个符号的出现,促使了L’Hospital 对分数阶导数问题的思考.1695年9月,L’Hospital 在写给莱布尼茨的信中问到:“一个函数()f x 的n 阶导数可以表示为()n n d f x dx ,如果当12n =时会有怎样的结果.” 莱布尼茨在回信中写道:“这显然是一个悖论,但总有一天会得出有用的结论.”由此,分数阶微积分诞生了,在之后300多年的学习研究过程中,莱布尼茨的这句话已经得到了验证,至少他说对了一半,尤其是在二十世纪,大量有关分数阶微积分的应用被人们所发现.尽管分数阶微积分有了这些应用以及一些数学背景,然而它的物理意义却很难把握,分数阶微积分的定义也不像整数阶微积分那样完善.历史上,莱布尼茨、欧拉、拉普拉斯、Lacroix 和傅里叶都曾对分数阶导数做出过重要贡献,其中,欧拉迈出了关键的第一步.他注意到,当n 时非整数时,幂函数a x 的导数na n d x dx在数学上有意义.1812年,拉普拉斯提出了积分形式的的函数()x T t t dt -⎰的分数阶导数的思想.1819年,Lacroix 第一次给出了1/21/21/2d x dx =1823年,Abel 在求解等时曲线的积分方程时,第一次使用分数阶算子并用分数阶微积分来表示该方程的解.1832年,Liouville 成功的应用了自己提出的分数阶导数的定义,解决了势理论问题,此后,Liouvile 发表的一系列文章使他成为分数阶微积分理论的创始人.19世纪末,物理学家Heaviside发表的一系列文章表明,分数阶算子可以应用于求解特定的整数阶微分方程,从数学角度看,他的方法并不严格,但被证明对电缆中电流的传输这类问题非常有效.后来Heaviside的结果被证明是正确的,但他的处理过程在数学上并不完善,直到1919年Bromwich才把这一工作完善,Heaviside的想法极大的促进了分数阶算子的发展,但当时分数阶微积分还没有被应用于科学和工程的物理和力学建模.20世纪40年代,力学家Scott和Gerasimov分别独立的提出了介于牛顿流体和胡克定律表征的分数阶导数模型.地理学家Caputo和Mainardi将分数阶微积分方法运用到复杂黏弹性和流变介质,发展了若干的力学模型,更为重要的是,Caputo发展了一个不同于传统的Riemann-Liouville分数阶导数的新定义(被人们称为Caputo定义),克服了前者的强奇异性,并且自然的将初始条件含在定义中,在解决实际问题时得到了非常广泛的应用.1965年,美国耶鲁大学的Mandelbrot教授提出分形的概念,并首次指出自然界和工程中存在大量分数维的事实,并且整体与部分之间存在自相似现象,他认为分数阶布朗运动与Riemann-Liouville提出的分数阶微积分定义有紧密的内在联系.从此,作为分形几何和分形动力学的基础,分数阶算子理论特别是分数阶微分方程的研究开始得到广泛关注,分数阶微积分的研究重点也逐渐从纯数学领域转移到其它学科领域.20世纪末至今,由于反常扩散、多孔介质力学、非牛顿流体力学、黏弹性力学等研究的需要,分数阶导数的研究和应用再度引起广泛重视,成为多个领域学者研究的热点.1.2分数阶微分方程的研究现状现实的世界本质上是分数阶的.过去用整数阶微积分描述自然界中的事物,但自然界中许多现象依靠传统的整数阶微积分是不能精确描述的,必须对传统的微积分进行拓展才能更好的描述与研究这样的现象.分数阶微分方程是扩展传统微积分学的一种直接方式,即允许微分方程中对函数的导数阶次选择分数,而不是现有的整数.目前,分数阶算子的定义主要有Riemann-Liouville型、Caputo型、Grunwald-Letnikov型,Weyl型、Erdelyi-Kober型、Riesz型以及Marchaud-Hadamard型分数阶微积分[]23,前面三种定义用的最为广泛,同时这三种定义还存在着一定关系:Riemann-Liouville型分数阶微分定义和Caputo型分数阶微分定义都是在Grunwald-Letnikov型分数阶微分定义的基础上进行改进而得到的,但是它们同时又有各自的侧重点,其中对于Riemanna-Liouville 型定义是从数学角度出发,在计算时初始条件是必要的,但这个定义在应用方面缺乏物理背景,是得它在很大程度上不能适用于具体模型;而对于Caputo 型定义,它正好弥补了Riemanna-Liouville 型定义可以很好的应用到数学模型中去,因为此时的初始条件恰好是整数阶的导数,这样连带的初始条件就可以很好的描述一些物理现象的特征,并能准确的对它进行计算,它还比较适用于拉普拉斯变换,有利于分数阶微分方程的进一步讨论随着分数阶微积分定义的出现,分数阶微分方程的求解就成为数学家至今仍在研究的主要课题,在分数阶微分方程的解析解研究方面:Agarwal []26,29,30研究了分数常微分方程解的存在性、唯一性;Miller 和Ross []21给出了一种分数阶微分方程的求解方法;Wyss []36等人给出了分数阶Black-Scholes 方程的一个完整解;Zhanbin Bai []40,41, Chuanzhi Bai []11等研究了分数微分方程正解的存在性;然而,由于分数阶微分方程的解析解以及基本解大多带有特殊函数(如多变量的Mittag-Lemer 函数),这些特殊函数的计算是相当困难的,而且并非所有的分数阶微分方程都能得到其解析解.因此,建立分数阶微分方程的数值方法是非常必要的,在分数阶微分方程的数值解研究方面:Diethelm []1314,等人对于Adams类型的分数阶微分方程,提出用预测校正法来得到微分方程的数值解并且讨论了分数阶非线性微分方程的求解问题,在特定初值和Riemann-Liouville 型分数阶微积分定义的条件下求解分数阶微分方程的数值解;Diethelm 和Ford []15在分数阶微积分的Caputo 定一下给出了给出了一种求解分数阶微积分的数值算法;Sayed []33等人对于线性分数阶微分方程给出了一种计算其近似的数值解的算法,该算法需要很大的计算量来得到计算权数;Agrawal []4等人在Caputo 型分数阶微积分的基础上,求解了分数阶漫射波方程,数值解在实际问题中得到了广泛的应用,数学家们给出了自己的解法,每种解法都随着计算机技术的快速发展得到了验证.在最近的十多年里,有关分数阶微分方程的论文和著作相继出现,在这些论文和著作当中,也有很大一部分文章是关注不同边值条件和不同阶数取值范围下的分数阶微分方程正解存在性和唯一性问题,在不同的边值条件和阶数条件范围下,可以采用不同方法来求解分数阶微分方程的解以及证明其正解的存在性.已知的求解方法中较多是采用各种推广的特殊函数来直接求解,其中Green 函数是研究的重点内容,不同的边值条件和阶数的取值范围会产生不同的Green函数以及相应的Green函数值的有所不同,进而导致在后续估计分数阶微分方程正解的存在的条件以及在证明正解存在性的方法上出现显著差别.本文主要利用非线性泛函分析中的不动点理论,把解的存在性转化为某个非线性算子不动点的存在性.研究了一类分数阶微分方程在边值条件下正解的存在性.第二章 带有奇异的非线性SemipositoneSturm-Liouville 边值问题解的存在性2.1 引言近年来,带有奇异的或非奇异的positone 问题(其中非线性项()f y 为非负值)的研究已引起了很多的学者的关注,详见文献[17,25,26,28].最近,文章[19,20]开始讨论了Semipositone 非奇异问题. 这里Semipositone 问题指的是非线性项()f y 可能在0y =处奇异并且f 可以取负值.本章主要研究了带有奇异的非线性Semipositone Sturm-Liouville 边值问题(2.1.1)解的存在性.0μ>这里是常数,1[0,1],q L ∈:(0,)f R ∞→连续并且在0y =处奇异, ,,,0:0.αβγδργβαγαδ≥=++>及在文献[27]中作者研究了带有奇异的Semipositone Dirichlet 边值问题 ()()(())0,01(0)(1)0;y t q t f y t t y y μ''+=<<⎧⎨==⎩解的存在性.受以上文献启发,本文讨论了带有奇异的Semipositone Sturm-Liouville 边值问题(2.1.1)解的存在性.本章主要利用锥上的不动点理论来建立边值问题解的存在性,本章第二部分首先介绍了一些基本定义和引理,给出我们后面用到的基本定理,第三部分是我们的主要定理并给出了(2.1.1)式当==1αγ,==0βδ特殊情形时的一个例子,边值问题(1)0,01(0)(1)0,0,0p q y y y t y y p q μ-''⎧++-=<<⎨==>>⎩ 当μ充分小时,有一个解()2[0,1](0,1),()0,0,1y C C y t t ∈⋂>∈且有()()(())0,01(0)(0)0,(1)(1)0;y t q t f y t t y y y y μαβγδ''+=<<⎧⎪'-=⎨⎪'+=⎩2.2 预备知识定理 2.2.1[]27(,)E E K E =⋅∈设是一个Banach 空间,是一个锥,,r R 都是常数且有r R <<0.{}(=:)R R A K K x E x R Ω⋂→Ω∈<假设:这里,A 是一个连续的紧映射并且假设下列条件成立:(1) (),[0,1)E r x A x x K λλ≠∈∈∂Ω⋂且, (2) ,E R Ax x x K >∈∂Ω⋂, 那么算子{}:A K x E r x R ⋂∈≤≤在集合上有一个不动点.引理 2.2.1[]27设{}[0,1]()0,[0,1]()[0,1]K y C y t t y t ∈≥∈=:并且是上的凸函数,如果y ,K ∈那么 01()(1),[0,1];=max ()t y t t t y t y y t ≤≤≥-∈这里. 引理2.2.2 1[0,1],()0,(0,1),q L q t t ∈>∈假设那么边值问题()()=0,01(0)(0)0,(1)(1)0;y t q t t y y y y αβγδ''+<<⎧⎪'-=⎨⎪'+=⎩(2.2.1) 的解0()(),[0,1]w t w t G t t C t ≤∈满足(,);G t t 其中(,)为边值问题 =0,(0)(0)0,(1)(1)0;y y y y y αβγδ''⎧⎪'-=⎨⎪'+=⎩其中1()(),01,1()(),01,t s s t G t s s t t s γδγβαργδγβαρ⎧+-+≤≤≤⎪⎪⎨⎪+-+≤≤≤⎪⎩(,)= 100[0,1]11max ()()()()()()t t t C s q s ds s q s ds t t ψϕψϕ∈⎧⎫=+⎨⎬⎩⎭⎰⎰ 记 ():,():,01t t t t t ϕγδγψβα=+-=+≤≤.证明:因为(2.2.1)式的解Green G t s 的函数(,)当t=s 的情形,1011()()()()()()()t t w t t s q s ds s t q s ds γδγβαγδγβαρρ=+-+++-+⎰⎰10()()1()()()()t t t t s q s ds s q s ds ϕψβαγδγρρρ=+++-⎰⎰ 10()()()()()()t t t t t t q s ds q s ds ϕψϕψρρ≤+⎰⎰ 所以有00()()()(,)t t w t C G t t C ϕψρ≤=. 引理2.2.3 :(0,)f R M ∞→>设的连续函数并且存在一个常数0,使得 ()0,f u M +≥(0,)u ∀∈∞,若边值问题*()()(()())0,01(0)(0)0,(1)(1)0;y t q t f y t t t y y y y μφαβγδ''⎧+-=<<⎪'-=⎨⎪'+=⎩(2.2.2) 211[0,1](0,1)()(),(0,1),y C C y t t t φ∈⋂>∈有一解且 ()=(),t Mw t φμ这里 *()(),0f v f v M v =+>.1()()()u t y t t φ=-那么 为(1.1.1).式的一个非负解证明:因为1()()()u t y t t φ''''''=-=*1()(()())()q t f y t t Mq t μφμ--+ []1()(()())()q t f y t t M Mq t μφμ=--++1()(()())()(())q t f y t t q t f u t μφμ=--=-所以有()u t ''=()(())q t f u t μ-,即1()()()u t y t t φ=-是(2.1.1).式的一个非负解2.3 主要结果假设下列条件成立:(H1):(0,)f R M ∞→>的连续函数并且存在一个常数0,使得()0,f u M +≥ (0,)u ∀∈∞.(H2)()()()f u M g u h u +=+,(0,)u ∀∈∞,其中0g ∞为(,)上正的连续递减函数 并且存在000()()(),0,0K g ab K g a g b a b >≤∀>>使得. h ∞为[0,)上的连续非负函数并且有0hg∞为(,)上的递增函数. (H3) 存在常数,(,)(1),L G t t Lt t ≤-使得 存在0,r MLC μ>使得00001,()11()r du K b gu MLC h r g r g r μμ>⎛⎫⎛⎫-+ ⎪⎪⎝⎭⎝⎭⎰ 1021002max 2(1)(),2(1)(),b t t q t dt t t q t dt ⎧⎫=--⎨⎬⎩⎭⎰⎰其中,100[0,1]11max ()()()()()()t t t C s q s ds s q s ds t t ψϕψϕ∈⎧⎫=+⎨⎬⎩⎭⎰⎰. (H4)11((1))(0,),()(,)2()((1))((1))a a Rg a a R a R r q s G s ds g R g a a R h a a R εμσεε--∈>≤-+-⎰存在有这里 00MLC Rμεε>≥是任意常数且满足1-,11[0,1]01()(,)sup ()(,)a aaat q s G s ds q s G t s ds ξξ--∈≤≤=⎰⎰满足.定理2.3.1 假设条件(H1)、(H2)、(H3)和 (H4)成立,那么边值问题(2.1.1) 式有一个解2[0,1](0,1)(0,1)()0.y C C t y t ∈⋂∈>且有当时证明:记0*0001:(1)m N m N a a R m ηηε⎧⎫=∈<<-⎨⎬⎩⎭且. 我们首先证明边值问题0,()()(()())0,011(0)(0),1(1)(1),m y t q t f y t t t y y m y y m N m μφαβγδ⎧''⎪+-=<<⎪⎪'-=⎨⎪⎪'+=∈⎪⎩(2.3.1) 对0m N ∀∈有一个解m y ,()0,()(),[0,1],,m m m y t y t t t r y R φ≥≥∈≤≤这里1()()(),()11()(),0.m f v M g v h v v mf vgh v v mm ⎧+=+≥⎪⎪=⎨⎪+≤≤⎪⎩欲证(2.3.1)式,我们接下来看下式*0,()()(()())0,011(0)(0),1(1)(1),m y t q t f y t t t y y m y y m N m μφαβγδ⎧''⎪+-=<<⎪⎪'-=⎨⎪⎪'+=∈⎪⎩(2.3.2) *1()()(),11()()(),01()(0),0.m f v M g v h v v m f v g h v v mm g h v m ⎧+=+≥⎪⎪⎪=+≤≤⎨⎪⎪+≤⎪⎩这里所以我们有*()0,(,).m f v v ≥∀∈-∞∞0,([0,1],){[0,1]()0,[0,1]()[0,1]}.m N E C y C y t t y t ∈=∈≥∈固定并且K=:且是上的凸函数()():[0,1]y t y t A K C →是边值问题(2.3.2)式的解当且仅当是算子1*01()()(()())(()())m Ay t G t s q s f y s s ds t t mμφϕψρ=-++⎰(,) (2.3.3)的不动点.:[0,1]A K C →由文献[27]知算子是连续的并且是全连续算子.:A K K →接下来验证*()0,(,)()0.m u K f v v Au t ∀∈≥∀∈-∞∞≥对,因为,所以有同时也容易看出()0A u t ''≤,(0,1),.t ∈ {}{}12=[0,1]:,[0,1]:.u C u r u C u R Ω∈<Ω=∈<设1,[0,1).y Ay y K λλ≠∈∈⋂∂Ω我们首先证明且=0.(0,1)=y Ay y Ay λλλλ≠∈当时,显然成立当时,假设成立,我们有*()()(()())0,01m y t q t f y t t t λμφ''+-=<< (2.3.4)00[0,1](0,1),(0,)()0,y t t t y t '∈∈≥由是凸函数可知,在区间上存在点使得当时有000(,1)()0,().t t y t t y t y r '∈≤==时有并且在处有0(,)()(1)=()(1)()(,),0,1(1)G t t y t t t y y t t t r w t G t t C L t t t ≥-≥-≤≤∈-因为以及,()0()()()=()1()1()MLC Mw t y t t y t y t y t r μμφ⎡⎤⎡⎤--≥-⎢⎥⎢⎥⎣⎦⎣⎦所以 0,r MLC μ>由于所以有0()()()10,0,1MLC y t t y t t r μφ⎡⎤-≥->∈⎢⎥⎣⎦() (2.3.5) *11()=()()()0m v f v g v h v v v m m≥+∞≤≤当时,,因为g 在(0,)上递减,所以当时,*1()()()()()m f v g h v g v h v m=+≤+,*(()())(()())(()()),(0,1)m f y t t g y t t h y t t t φφφ-≤-+-∈故有(0,1)(2.3.4)x ∈因此当时,由式我们有()()(()())(()())y x q x g y x x h y x x μφφ''≤-+--(()())()(()())1(()())h y x x q x g y x x g y x x φμφφ⎧⎫-=-+⎨⎬-⎩⎭(2.3.5)由式,我们有0(()())()()()11(()())MLC h y x x y x q x g y x r g y x x μφμφ⎧⎫-⎡⎤''≤-+⎨⎬⎢⎥-⎣⎦⎩⎭-() 00()1)1()())(MLC h r K gq x g y x rg r μμ⎧⎫≤-+⎨⎬⎩⎭(() (2.3.6) 不等式00()t t t t ≤两边从到积分得,00()()(())1)1()(t t MLC h r y t K g y t g q x dx rg r μμ⎧⎫'≤-+⎨⎬⎩⎭⎰() (2.3.7) =y Ay λ由知,1*01()()(,)()(()())(()())m y t Ay t G t s q s f y s s ds t t mμφϕψρ≤=-++⎰ 1*1(,)()(()())(()())m G t t q s f y s s ds t t mμφϕψρ≤-++⎰1*1(1)()(()())(()())m Lt t q s f y s s ds t t mμφϕψρ≤--++⎰ 因此++(0)(1)y y m mγδβαδβρρ++≤≤有,. 取++=max ,m mm γδβαδβηρρ⎧⎫++⎨⎬⎩⎭0(2.3.7)t ,对式两边从0到积分,以及由分部积分得00+00()1)1()()(rt mMLC duh r K g xq x dx g u r g r γδβρμμ+⎧⎫≤-+⎨⎬⎩⎭⎰⎰() 有 0000()11)11)()()(1mrt MLC duh r K g x x q x dx g u r g r t ημμ⎧⎫≤-+-⎨⎬-⎩⎭⎰⎰(()000(2.3.6()1t t t t t ≥类似的,如果我们对)式两边先从到积分,然后再对不等式两边从到积分得m1000()11)11)()()(rt MLC duh r K g x x q x dx g u r g r t ημμ⎧⎫≤-+-⎨⎬⎩⎭⎰⎰(()有 000()1)1()(mrMLC duh r K b g g u r g r ημμ⎧⎫≤-+⎨⎬⎩⎭⎰() (2.3.8) 其中0b 为条件(H3)中所定义,又因为由条件(H3)有00001()11()r du K b gu MLC h r g r g r μμ>⎛⎫⎛⎫-+ ⎪⎪⎝⎭⎝⎭⎰ 所以当η充分小时有0001()11()r du K b gu MLC h r g r g r ημμ>⎛⎫⎛⎫-+ ⎪⎪⎝⎭⎝⎭⎰ 与(2.3.8)式产生矛盾.接下来我们证明,当 2y K Ay y ∈⋂∂Ω>时,有,2K ⋂∂Ω因为当时, 有()(1),[0,1]y t t t R t ≥-∈.0(0,1)()()()1()(1)MLC t y y t t y t y t t t R r μφεε⎡⎤∈∈-≥-≥≥-⎢⎥⎣⎦当时, [],1t a a ∈-因此,当时,我们有 ()()(1),y t t a a R φε-≥-1()()(1)y t t a a R m φε-≥->由 有 []*(()())=(()())(()()),,1mf y t tg y t th y t t t a a φφφ--+-∈- 1*01()()(()())(()())m Ay G s q s f y s s ds mξμξφϕξψξρ=-++⎰(,) 1*()(()())a m a G s q s f y s s dsμξφ-≥-⎰(,)1(()())()(()())1(()())a ah y s s G s q s g y s s ds g y s s φμξφφ-⎧⎫-=-+⎨⎬-⎩⎭⎰(,)1((1))()1()((1))aa h a a R g R G q s ds g a a R εμξε-⎧⎫-≥+⎨⎬-⎩⎭⎰(,s)(),(2.9).Ay R y ξ≥=由条件(H4)知因此式成立211.2.1(\),(1),[0,1]m m m A y r y R y t t r t ∈ΩΩ≤≤≥-∈由定理知有不动点并且有.0(1)(1)()()m y t t r MLC t t Mw t t μμφ≥->-≥=因为m y 所以是边值问题(2.3.1)式的解.{}0[0,1].m m N y ∈下证是定义在区间上的有界,等度连续的函数族因为 *(()())(()())(()()),(0,1).m m m m f y t t g y t t h y t t t φφφ-≤-+-∈(()())()()(()())1(()())m m m m h y t t y x q x g y t t g y t t φμφφ⎧⎫-''≤-+⎨⎬-⎩⎭所以我们有-00()()1)1()())()m m MLC h R y x K g q x gy x rg R μμ⎧⎫''≤-+⎨⎬⎩⎭-(( (2.3.9) 0(),()()()1()1,(0,1)()m m m m m MLC Mw s r y R y s s y s y s s y s r μμφ⎡⎤⎡⎤≤≤-=-≥-∈⎢⎥⎢⎥⎣⎦⎣⎦又因为,(0,1),(0,)()0,(,1)()0,m m m m m t t y t t y t ''∈≥≤同时存在使得在区间上在区间上 (2.3.9)()m m t t t t ≤对式两边从到积分得00()()1)1()())(mt m t m y t MLC h R K g q x dx g y t r g R μμ'⎧⎫≤-+⎨⎬⎩⎭⎰(() (2.3.10)(2.1.9)()m m t t t t ≥另一方面式两边从到积分得00()()1)1()())(m tm t m y t MLC h R K g q x dx g y t r g R μμ'⎧⎫-≤-+⎨⎬⎩⎭⎰(()(2.3.11) 由(2.3.10)、(2.3.11)式可知'00()()1)1(),(0,1)())(m m y t MLC h R K g v t t g y t rg R μμ⎧⎫≤-+∈⎨⎬⎩⎭(()(2.3.12) 其中{}{}10max ,min ,()()t a t a v t q x dx =⎰,{}{}00010inf :sup :1m m a t m N t m N a <<∈≤∈<<.注:这里0,1()0m m m t y t '=是区间()上唯一的一点,满足,有{}0inf :0m t m N ∈>. 倘若不成立,那么存在0N 的子列,使得0m m t →∞→当时,有. 对(2.3.10)式两边从0m t 到积分可得()00000()1)1()()(()m m m m y t t MLC duh R du K g xq x dx g u r g R g u ημμ⎧⎫≤-++⎨⎬⎩⎭⎰⎰⎰() (2.3.13) 因为0,0m m m t η→∞→→当时,有,由(2.3.13)式可知m →∞当时有()0m m y t →,然而因为()m m y t 在区间[0,1]的最大值在m t 处取得,所以当0m m y →∞→当时,有这与()(1),[0,1]m y t t t r t ≥-∈矛盾 故有{}0:m 0m inf t N ∈>.类似的可以证明{}0sup :1m t m N ∈<. 定义映射0:[0,)[0,),()()y duI I y g u ∞→∞=⎰,显然{}0)m m N I y ∈(是有界的.(())(())m m I y t I y s -()()()()(())m m y t t my s sm y x dudx g u g y x '==⎰⎰00()1)1()(ts MLC h R K g v x dx rg R μμ⎧⎫≤-+⎨⎬⎩⎭⎰() 可知{}0)m m N I y ∈(是等度连续的. 由,[]1()I I R -在区间0,上的一致连续性以及11())(()=((()))((()))m m m m y t I y s I I y t I I y s ----知{}0[0,1]m m N y ∈是定义在区间上的有界,等度连续的函数族.由Arzela Ascoli -定理知,存在0[0,1],N N y C ∈的子列以及函数 m →∞当时m y y 有在区间[0,1]上一致收敛到同时有(0)(0)0y y αβ'-=(1)(1)0y y γδ'+=,r y R ≤≤,()(1),[0,1]y t t t r t ≥-∈,且有0()(1)(1)()()y t t t r MLC t t Mw t t μμφ≥->-≥=. 固定(0,1)t ∈,不失一般性,我们假设12t >,固定(0,1),x x t ∈>满足,对1[,]2s x ∀∈ 0()()()()1()1()MLC Mw s y s s y s y s y s r μμφ⎡⎤⎡⎤-=-≥-⎢⎥⎢⎥⎣⎦⎣⎦00(1)(1)112MLC MLC x s s r r r r μμ-⎡⎤⎡⎤≥--≥-⎢⎥⎢⎥⎣⎦⎣⎦选择0001(1)12MLC x n N r n r μ-⎡⎤∈<-⎢⎥⎣⎦使得,,设{}10:.N m N m n =∈≥ 当1,m y m N ∈时,由泰勒公式有12111()()()222x m m m y x y y x f s x s ds⎛⎫⎛⎫⎛⎫'''=+-+- ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭⎰=[]12111()(()())(()()()222x m m m m y y x q s g y s s h y s s s x ds μφφ⎛⎫⎛⎫⎛⎫'+-+-+-- ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭⎰因为(1)()m rs s y s R -≤≤,所以11,2my m N ⎧⎫⎛⎫'∈⎨⎬ ⎪⎝⎭⎩⎭为有界序列,[]0,1s ∈. 故112m m N y ∈⎧⎫⎛⎫'⎨⎬ ⎪⎝⎭⎩⎭有一个收敛的子列,不妨设子列的极限收敛到0r R ∈, 在1N 中当m →∞时,我们有,[]10211()()(()())(()()()22x y x y r x q s g y s s h y s s s x ds μφφ⎛⎫⎛⎫=+-+-+-- ⎪ ⎪⎝⎭⎝⎭⎰(2.3.14)对(2.3.14)式两边求二阶导有[]()()(()())(()())0y x q x g y x x h y x x μφφ''+-+-=所以有 [][]()()(()())(()())0,0,1y t q t g y t t h y t t t μφφ''+-+-=∈*()()(()())0,01y t q t f y t t t μφ''+-=<<.因此y 为(2.2.2)式的解并且有()()y t t φ>,(0,1)t ∈. 下面我们通过一个实例来给出了定理2.3.1的一个应用. 例:考虑边值问题(1)0,01(0)(1)0,p q y y y t y y μ-''⎧++-=<<⎨==⎩ (2.3.15) 这里0(0,)μμ∈且满足()100(1)12pp μμ++≤. (2.3.16)那么边值问题(1.3.15)式有一个解()0,(0,1)y y t t ≥∈且.我们将应用定理 2.3.1来证明,边值问题(2.3.15)式是(2.1.1)式当==1αγ,==0βδ的特殊情形. 设01,(),(),1p q M g y y h y y K -====,1L =,14a =,其中 100[0,1]11max ()()()()()()t t t C s q s ds s q s ds t t ψϕψϕ∈⎧⎫=+⎨⎬⎩⎭⎰⎰10[0,1]111max 112t t t sds sds t t ∈⎧⎫=+-=⎨⎬-⎩⎭⎰⎰. 11210021max 2(1),2(1)6b t t dt t t dt ⎧⎫=--=⎨⎬⎩⎭⎰⎰.01122r MLC r μμμ==<≤=取时有,,()1001112121121()11()ppp r p q du r gu r rp p MLC h r g r g r μμμ++⎛⎫- ⎪⎛⎫⎝⎭=-= ⎪+++⎛⎫⎛⎫⎝⎭-+ ⎪⎪⎝⎭⎝⎭⎰ (2.3.16)由式有()()00000011122222121()11()p pr du K b p p gu MLC h r g r g r μμμμμμ⎛⎫⎛⎫-- ⎪ ⎪⎝⎭⎝⎭≤<≤≤=++⎛⎫⎛⎫-+ ⎪⎪⎝⎭⎝⎭⎰ 取01111222MLC R R R μμε=>=-≥,当时有,1- . 最后当,1R q →∞>时有,13((1))320()((1))((1))333232pq p qp qR Rg a a R g R g a a R h a a R Rεεε-+-+⎛⎫⎪-⎝⎭=→-+-⎛⎫⎛⎫+ ⎪ ⎪⎝⎭⎝⎭. 因此定理2.3.1的条件(H1)、(H2)、(H3)、(H4)均满足,故边值问题(2.3.15)式有一个解()0,(0,1)y y t t ≥∈且.第三章 带有奇异的非线性分数阶微分方程D i r i c h 边值问题正解的存在性3.1 引言近年来,人们开始并越来越多的关注、研究分数阶微分方程,主要是因为分数阶微积分自身理论的发展以及在多种学科中的应用,例如物理学,化学,工程学等等,详见文献[8,23,30,31].分数阶微分方程的Dirichlet 边值问题是很多学者研究是的焦点,在文献[40]中作者研究了边值问题()(),()0(0)(1)0D y t f t y t y y α+===正解的存在性和多解性,这里[][]()12,0,1,0,f C α<≤∈∞为非负函数,Bai Zhanbin 通过Krasnosel’skii 不动点定理和Leggett-Williams 不动点定理得到了相关结论.在文献[32]中,Zhang 研究了边值问题()(2)()(),,0,01,1n D u t q t f u u u t n n αα-'+=<<-<≤, (3.1.1)(2)(2)(0)(0)(0)(1)0,n n u u u u --'===== (3.1.2)这里0D α+是标准的Riemann-Liouville 分数阶导数,q 可能在t=0处奇异,f 可能在(2)0,0,0n u u u -'===处奇异.在此基础上Goodrich 在文献[42]中研究了边值问题()0(),(),01,1,vD y t f t y t t n v n +-=<<-<≤ (3.1.3) ()0,02,i y i n =≤≤-(0) (3.1.4)00()0,12,t D y t n αα+=⎡⎤=≤≤-⎣⎦ (2.1.5)这里的3n >,可以看出边值条件(3.1.2)式是边值条件(3.1.4)和(3.1.5)式的特殊情形,文献[42]在Zhang 研究的基础上进一步阐述了Green 函数的有关性质Harnack-like 不等式,这是在锥上寻找正解存在性的一个重要性质.文献[28]中Agarwal 等研究了边值问题()(),(),()0,D u t f t u t D u t αμ+= (3.1.6) (0)(1)0,u u == (3.1.7)正解的存在性,这里12,0,1αμαμ<<>-≥,0D α+是标准的Riemann-Liouville 分数阶导数.f 是正的Caratheodory 函数并且(,,)0f t x y x =在处奇异,通过锥上的不动点定理以及函数列的相关性质证明了边值问题(3.1.6)、(3.1.7)式正解的存在性.本章主要在文献[28,42]的基础上研究下面的带有奇异的非线性分数阶微分方程Dirichlet 边值问题00()(,(),())0,01D x t f t x t D x t t αβ+++=<< (3.1.8) ()(0)0,02i x i n =≤≤- (3.1.9)01()0,02t D x t n μμ+=⎡⎤=≤≤-⎣⎦ (3.1.10)正解的存在性.这里1n n α-<≤,01βα<≤-,f 是正的Caratheodory 函数并且在[0,1],(0,)B B ⨯=∞⨯上满足Caratheodory 条件(([0,1]f C a r B ∈⨯,(,,)0f t x y x =在处奇异,0D α+是标准的Riemann-Liouville 分数阶导数.我们说函数f 在集合[0,1],(0,)B B R ⨯=∞⨯上满足Caratheodory 条件,如果函数f 满足下面三个条件:[]()(,,):0,1(,)a f x y x y B →∀∈是可测函数,成立, [](b)(,,):0,1f t B t →∈是连续的,a.e.成立, []1()0,1,c B L κκϕ∈对中的任一紧集,存在函数使得[](,,),0,1(,)f t x y t x y B κϕ≤∈∀∈a.e.,成立,函数[]0,1u C ∈称为边值问题(2.1.8)-(2.1.10)的一个正解,如果x 在区间(0,1)上有0x >,[]00,1D x C β+∈,[]100,1D x L μ+∈且满足边值条件(3.1.9)、(3.1.10)式和等式(3.1.8),对几乎所有的[]0,1t ∈成立.本文中假设函数f 满足下列条件:[]()1():0,1,(0,),H f Car C B B ∈⨯=∞⨯[]0lim (,,),..0,1,x f t x y a e t y +→=∞∈∀∈(3.1.11)并存在正整数m 满足(,,)(1)f t x y m t μ≥-,[]..0,1,(,)a e t x y ∈∈∀∈(3.1.12)()[]2():(,,)()()()(),..0,1,(,)H f t x y t q x p x y a e t x y B γω≤++∈∈∀∈这里[]()[)1()0,1,0,1,,0,1t L q C p C γω∈∈∈都是正的函数,其中q 单调递减,,p ω单调递增,且有[]10()((1))s q K s s ds αγ-<∞⎰,()1mK α=Γ+ (3.1.13)()()lim0x p x x xω→∞+= (3.1.14)因为(3.1.8)式是一个奇异方程,故我们定义1(,,)11(,,)0n f t x y x n f f t y x n n⎧≥⎪⎪=⎨⎪≤<⎪⎩有[]()[)**0,1,0,n f Car C B B ∈⨯=∞⨯,n ∈,由条件1()H 和2()H 可得[]*1(,,)()()()(),..0,1,(,)n f t x y t q p x y a e t x y B n γω⎛⎫≤++∈∈∀∈ ⎪⎝⎭ (3.1.15)[]1(,,)()()()(),..0,1,(,)n f t x y t q p x y a e t x y B n γω⎛⎫≤++∈∈∀∈ ⎪⎝⎭(3.1.16)接下来我们首先讨论一般的分数阶微分方程00()(,(),())0,01n D x t f t x t D x t t αβ+++=<< (3.1.17) 3.2预备知识定义2.1[]40空间[]0,1C 上的范数[]{}max ():0,1x x t t =∈,空间[]10,1L 上的 的范数1()Lxx t dt =⎰.定义 3.2[]40函数:(0,)(0)y R Riemann Liouville α∞→>-,阶数为的分数阶积分由以下公式给出:1001I ()()()()ty t t s y s ds ααα-+=-Γ⎰ (3.2.1) 上式右端在(0,)∞上有定义,其中10()s e s ds αα∞--Γ=⎰.定义3.3[]40函数:(0,)(0)y R Riemann Liouville α∞→>-,阶数为的分数阶微分由以下公式给出:0101()()()()n tn d y s D y t ds n dt t s ααα+-+⎛⎫= ⎪Γ--⎝⎭⎰ (3.2.2)上式右端在(0,)∞上有定义,其中[]1n α=+,[]α表示实数α的整数部分.引理3.1[]28关于分数阶微积分有如下性质:(1)00I ()()D y t y t αα++=,..(0,1]a e t ∈, 1(0,1)y L ∈, 0α>(2)如果0α>,0λ>,那么110()()D ttαλλαλλα---+Γ=Γ- (3) []10()(,1),0,1tt s s ds t B t βααβαβ----=-∈⎰,其中B 为Beta 函数1110(,)(1),0,0p q B p q x x dx p q --=->>⎰,()()(,)()p q B p q p q ΓΓ=Γ+(4) [][]1()(),0,1,()0,1.I I f t I f t t f t L αβαβ+=∈∈ 由性质(4)可知()111()()()(,)()()t s tt s s f d ds B t s f s ds αβαβττταβ+-----=-⎰⎰⎰引理3.2[]40设0α>,如果(0,1)(0,1)y C L ∈⋂,那么分数高阶微分方程0()0D y t α+=有唯一解1212()N N y t C t C t C t ααα---=+++,,1,2,,i C R i N ∈=,其中N 是大于或等于α的最小整数.引理 3.3[]40设0α>,如果(0,1)(0,1)y C L ∈⋂且关于α的分数阶导数0()(0,1)(0,1)D yt C L α+∈⋂,那么120012I ()()N N D u t u t C tC t C t ααααα---++=++++ (3.2.3)其中,,1,2,,i C R i N ∈=,N 是大于或等于α的最小整数.引理3.4[]42设[]0,1f L ∈,那么边值问题0()()0,01,1D u t f t t n n αα++=<<-<≤, ()(0)0,02i u i n =≤≤-01()0,02t D u t n ββ+=⎡⎤=≤≤-⎣⎦有唯一的解1()(,)()y t G t s y s ds =⎰,其中11111(1)(),01()(,)=(1),01()t s t s s t G t s t s t s ααμαααμαα-------⎧---≤≤≤⎪Γ⎪⎨-⎪≤≤≤⎪Γ⎩证明:由引理2.3知,边值问题的解为12120()I ()n n u t C t C t C t f t αααα---+=+++- (3.2.4)由边值条件(1.4)式知230n C C C ====,对上式两边求μ阶导数,由引理2.1以及边值条件(2.1.5)式知11010()1()()()()()tD u t C t t s y s ds μαμαμααμαμ----+Γ=--Γ-Γ-⎰ 当1t =时有,1110()10(1)()()()C s y s ds αμααμαμ--Γ=--Γ-Γ-⎰,故有 11101(1)()()C s y s ds αμα--=-Γ⎰ 1111001()(1)()()()()()t t y t s y s ds t s y s ds ααμααα----=---ΓΓ⎰⎰111111011((1)())()(1)()()()t tt s t s y s ds t s y s ds ααμαααμαα-------=---+-ΓΓ⎰⎰ 设[][]{}0,1:0,1X x C D x C β=∈∈,给空间X 赋以范数{}*max ,x x D x β=,由文献[14]知X 是Banach 空间.定义空间X 中的锥P ,[]{}:()0,0,1.P x X x t t =∈≥∈为了证明边值问题(3.1.9)、(3.1.10)、(3.1.17)有一个正解,我们首先通过公式定义锥上的算子n T ,10()(,)(,(),())n n T x G t s f s x s D x s ds β=⎰ (3.2.5)引理3.5如果条件1()H 和2H ()成立,那么:n T P P →是一个全连续算子. 证明:设x P ∈,因为[]*(0,1)n f Car B ∈⨯,所以[]10,1n f L ∈,故有10()(,)(,(),())n n T x G t s f s x s D x s ds β=⎰1110(1)(,(),())()n t s f s x s D x s ds ααμβα---=-Γ⎰ 101()(,(),())()t n t s f s x s D x s ds αβα---Γ⎰ []10()(,(),())0,1(,)0tn t s f s x s D x s ds C G t s αβ--∈≥⎰由以及,可知()()[]()()0,1,0n n T x t C T x t ∈≥ (3.2.6)接下来由引理3.1的性质(3)、(4)可知()()()()()()101nn tn n d D T x t t s T x s ds n dt βββ--⎛⎫=- ⎪Γ-⎝⎭⎰()()()111(,)(,(),())nn t nd t s G s f x D x d ds n dt ββτττττβ--⎛⎫=- ⎪Γ-⎝⎭⎰⎰=()()1111001s (1)(,(),())()nn tn d t s f x D x d ds n dt αβαμβτττττβα-----⎛⎫-- ⎪Γ-Γ⎝⎭⎰⎰()()110011()(,(),())()nn t s n d t s s f x D x d ds n dt βαβτττττβα---⎛⎫--- ⎪Γ-Γ⎝⎭⎰⎰()()111101s (1)(,(),())()nn t n d t s ds f x D x d n dt βααμβττττταβ-----⎛⎫=-- ⎪ΓΓ-⎝⎭⎰⎰ ()()()+101,(,(),())()nn t n d B n t f x D x d n dt αββαβττττταβ--⎛⎫--- ⎪ΓΓ-⎝⎭⎰()()1+1101,(1)(,(),())()nn n d t B n f x D x d n dt αβαμβαβττττταβ----⎛⎫=-- ⎪ΓΓ-⎝⎭⎰ ()()()+101,(,(),())()n n t n d B n t f x D x d n dt αββαβττττταβ--⎛⎫--- ⎪ΓΓ-⎝⎭⎰ ()()(11101,(1)(,(),())()i n n B n it f x D x d n αβαμβαβαβττττταβ----≤≤--∏-+=-ΓΓ-⎰()10(,(),())tn t f x D x d αββτττττ--⎫--⎪⎭⎰所以有()()()(11101(1)(,(),())i n nn iD T x t ts f s x s D x s ds n βαβαμβαβαβ----≤≤-∏-+=-Γ+-⎰()1(,(),())tn t s f s x s D x s ds αββ--⎫--⎪⎭⎰ (3.2.7)因此有[]0,1,:n n D T x C T P P β∈→.为了证明n T 是一个连续算子,设{}m x P ⊂是一个收敛序列而且有*lim 0m m x x →∞-=,可知*,m x M m ≤∀∈对,M 是一个正的常数,因为[]*(0,1)n f Car B ∈⨯,我们有[]lim (,(),())(,(),()),..0,1n m m n m f t x t D x t f t x t D x t a e t ββ→∞=∈由(2.1.15)、(2.1.16)式可知,10(,(),())()()()()n m m f t x t D x t t q p M M n βγω⎛⎫<≤++ ⎪⎝⎭(3.2.8)由Lebesgue 控制收敛定理有1lim (,(),())(,(),())0n m m n m f t x t D x t f t x t D x t dt ββ→∞-=⎰ (3.2.9)()10()()()()(,)(,(),())(,(),())n m n n m m n T x t T x t G t s f s x s D x s f s x s D x s ds ββ-=-⎰ ()1(1,)(,(),())(,(),())n m m n G s f s x s D x s f s x s D x s ds ββ≤-⎰10(,(),())(,(),())n m m n f s x s D x s f s x s D x s ds ββ≤-⎰()()()11101()()()()(1)(,(),())(,(),())i n n m n n m m n iD T x t D T x t ts f s x s D x s f s x s D x s dsn ββαβαμββαβαβ----≤≤-∏-+-≤--Γ+-⎰()()()101(,(),())(,(),())ti n nmm n it s f s xs D x s f s x s D x s dsn αβββαβαβ--≤≤-∏-++--Γ+-⎰()1012(,(),())(,(),())i n n m m n if s x s D x s f s x s D x s ds n ββαβαβ≤≤-∏-+≤-Γ+-⎰故有lim 0n m n m T x T x →∞-=,所以n T 是连续算子.最后,设P X Ω⊂是中的有界集,*,x x L ∀∈Ω≤有,L 是一个正的常数,由于[]*(0,1)n f Car B ∈⨯,所以存在[]10,1L ϕ∈使得()[]0(,(),())..0,1,n f t x t D x t t a e t x βϕ<≤∈∀∈Ω对[],0,1x t ∀∈Ω∈有,11()()(,)(,(),())(,(),())n n n L T x t G t s f s x s D x s ds f s x s D x s ds ββϕ=≤≤⎰⎰()11101()()(1)(,(),())i n n n iD T x t t s f s x s D x s ds n βαβαμβαβαβ----≤≤-∏-+=-Γ+-⎰()1(,(),())tn t s f s x s D x s ds αββ----⎰()012i n L i n αβϕαβ≤≤-∏-+≤Γ+-故()n T Ω是X 中的有界集,下证n T 是等度连续的,设1201t t ≤<≤,1111212101()()()()()(1)(,(),())()n n n T x t T x t t t s f s x s D x s ds αααμβα-----=--Γ⎰111120(()())(,(),())t n t s t s f s x s D x s dsααβ--+---⎰2112()(,(),())t n t t s f s x s D x s dsαβ---⎰1111111211201()(1)()(()())()()t t t s s ds t s t s s ds αααμααϕϕα------≤--+---Γ⎰⎰。

二阶椭圆型方程奇摄动边值问题的渐近解

二阶椭圆型方程奇摄动边值问题的渐近解

如 一
即可 ( f , Y ) 满足
a ∈ 2 塑 O y



一 >0 ;
可 ( , 0 ) 一0 ,
∈ 0 ;
可 ( O , Y ) = ( ) , 利用 虿 ( 0 ) :0 , ( 0 ) =0 , 得
O y 2一 e 一
三 。 ( ) 昙+ 6 ( ) G q 2 + c ( ) 器 舢 ) + e ( 训 ) 南 圳
其中系数 0 , b , c , d , e , ,都是有界区域 Q C R 上 的连续函数且 a ( x , Y ) >0 , / ( x , Y ) 0 . 若 存在两次连续可微函数 ( , Y ) 和q Y ( x , Y ) 使得 l ( ) I [ 一 皿 ] , V ( x , Y ) ∈Q,
本文在文f71的基础上进一步讨论一类二阶线性椭圆型方程的奇摄动边值问题利用边界层函数法89构造出问题的零次形式近似并应用椭圆型算子的最大值原理作出解的渐近估计
第 4 6 卷 第 1 期 2 0 1 3 年 3 月
数 学 研 究
J o u r n a l o f Ma t h e ma t i c a l S t u d y
式近似 , 并应用椭 圆型算子 的最大值 原理对 问题 的解作 出渐近估计.
关键词 奇摄动 ; 椭圆型方程; 边值 问题 ; 边 界层函数法 ; 最大值原理 O 1 7 5 . 2 5 文 献标识 号 A
中图分类 号
1 引 言
1 9 5 0年 , L e v ms o n [ 1 】 首先研究 了一类二阶椭 圆型方程的 D i r i c h l e t问题, 在构造 形式近 似式 的基础上, 应用椭 圆型方程 的极值原 理给 出精 确解 的渐近估计 . 随后涌现 出大批研 究 偏微分方程 奇摄动边 值 问题 的论文 [ 2 - 6 ] . 1 9 9 6年 , D e J a g e r 和江福汝 【 7 】 综合 阐述 了奇异 摄动理论和方法, 并对 二阶和 高阶偏微分方程 的奇摄动 问题作 了系统的研 究. 本文在文 【 7 ] 的基础上进一步讨论一类二 阶线 性椭 圆型方程 的奇摄动边值 问题, 利 用边 界层 函数法 【 8 , 9 j 构造 出问题 的零次形 式近似 , 并应 用椭 圆型算子 的最大值原理作 出解 的渐 近 估计 . 引理 1 【 7 J ( 椭圆型算子的最大值原理) 设

非线性分数阶微分方程奇异边值问题的唯一解

非线性分数阶微分方程奇异边值问题的唯一解
在 t=0或 者 t=1处有 奇性 。
( >0 R e n —iu ie 数 阶积分 定义 如下 : )i ma n Lo v l 分 l

y) 志 Jfs (d ( o —) ss ( y),
连续 函 数 Y ( 。 R 的 O 阶 : 0, 。) t
假设 右边 是 逐点定 义 在 ( ∞ )上 的 。 0, 定义 2 (/> ) ima nLo vl O 0 R e n .iu ie分数 阶微 分定 义如 下 : l
其性质 , 中2<o 3是 实数 , + 其 z <  ̄ 是标准 Re a nLov l i n —i ie型微分 , m ul 并利 用锥不 动点定理和 混合单调方 法证 明了奇异边值 问
题 解 的 唯 一 性 。 最 后 举 例 加 以说 明 。
关键 词
分数 阶微分 方程
奇异边 值 问题

2 1 S i eh E gg 0 c T c. n r 1 . .
非 线性分 数阶 微分 方程 奇异 边值 问题 的唯一 解
于 瑶
( 连 教 育 学 院 , 连 16 2 ) 大 大 10 1


研 究了非线性分 数阶微分 方程边值 问题 D + ()+ t )= 0<t ; ( )=Z1 u t , ) 0, ( <1 u o / )=U( )= , G en函数及 ( 0 0 的 r e
质有 着广 泛 的理论 和现 实 意义 ¨ 1 。 _ J
文 献 [ 的作 者 应 用 Kansl i 8] rsoes i 动 点 理 k s不 论 、 合单 调算 子 的 唯一 不 动 点 理 论 研 究 了非 线性 混
分 数 阶微 分方 程边 值 问题 D + ()+a ( , () =0, <t<1 M f t“ £ ) 0 ,

奇异微分方程边值问题的数值解法

奇异微分方程边值问题的数值解法

奇异微分方程边值问题的数值解法奇异微分方程(singular differential equation)是指微分方程中存在奇异点(singular point)的一类特殊微分方程。

这些奇异点通常是导致方程在一些点上不连续或无定义的地方。

差分法(finite difference method)是将微分方程转化为差分方程,并用差分方法进行逼近求解的一种方法。

它的基本思想是将区间离散化,将微分方程转化为一个线性方程组,然后通过求解线性方程组得到数值解。

差分法的步骤如下:1.将求解区间进行等距离离散化,将连续的问题转化为离散的问题。

2.将微分方程中的导数用中心差分或向前/向后差分表示,得到差分方程。

3.将边界条件转化为差分方程中的代数方程。

4.将离散化的差分方程和代数方程组成一个线性方程组,然后通过求解线性方程组得到数值解。

有限元法(finite element method)是一种将微分方程用虚位移法(variational principle)得到弱形式,然后通过离散化和近似求解的方法。

它的基本思想是将求解区域划分为有限个子区域,然后在每个子区域内选取适当的基函数,通过这些基函数的线性组合近似原方程。

有限元法的步骤如下:1.将求解区域划分为三角形或四边形的有限个子区域,每个子区域称为单元。

2.在每个单元内选取适当的基函数,通常为多项式函数。

3.将原方程化为弱形式,即将方程两边乘上一个测试函数,并在整个求解区域上进行积分。

4.在每个单元内进行积分近似,并通过对各个单元的积分进行求和,得到离散化的方程。

5.将边界条件转化为代数方程。

6.将离散化的方程和代数方程组成一个线性方程组,然后通过求解线性方程组得到数值解。

总结起来,奇异微分方程边值问题的数值解法包括差分法和有限元法。

这两种方法都需要将微分方程进行离散化,然后通过求解线性方程组得到数值解。

选择使用哪种方法主要取决于具体的问题和求解精度要求。

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

ς ≤ ς0 := min{ψ(T )(K(T ))2/(q−1),
m2θ0 c0
1/(q−1)
}.

u¯(t) = λF (t), t ≥ T, λ ∈ [λ0, ∞),
u(t) = ςF (t), t ≥ T, ς ∈ (0, ς0].
通过直接计算, 得到,
F
(t)
=
−q
2 −
1
k(t)(K
(t))−(q+1)/(q−1),
(b2)
lim
t→∞
d dt
K (t) k(t)
:= Ck ∈ [0, ∞), k(t) :=
t
b(t), K(t) = k(s)ds.
0
问题(1.1)来源于边界blow-up的非线性椭圆型问题解在边界附近的精确渐近行为的研究
中,见文献[1]- [17].
最近,S. Cano-Casanova, J. L´opez-G´omez [1] 得到了结果:
设b 满足
(B1) b ∈ C[0, ∞),b 在(0, ∞)上是正的单调递增函数;
g 满足(g1) 和(g2), 则问题(1.1) 存在唯一整体解ψ ∈ C2(0, ∞). 而且,如果g 还满足(g3)和
(g5) 存在p > 1 使得
C0
:=
lim
s→∞
g(s) sp

(0, ∞);
1国家自然科学基金资助课题10671169.
(1.2)
Γ(t) =
R t
ds A(s)
,
A(t) =
t
(p+1)/(p−1)
(b(τ ))1/(p+1)dτ
, t ∈ (0, R].
0
(1.3)
注意到,(g1) 和(g4)隐含(g2). 但文[1]未讨论ψ在无穷远处的精确渐近行为. 本文在权函数b满足(b1)和(b2)的条件下,应用摄动方法和上下解方法,构造新的上下解, 得到了ψ在无穷远处的精确渐近行为.
本文的主要结果是
定理1.1. 设b 满足(b1)- (b2), g 满足(g1) - (g4),则问题(1.1)的唯一整体解ψ ∈ C2(0, ∞)
满足
lim ψ(t)(K(t))2/(q−1) =
t→∞
2((q − 1)Ck + 2) c0(q − 1)2
1/(q−1)
.
(1.4)
特别,当g(u) = uq, b(t) = tσ,σ > 0 时,
[3] F. Cˆırstea, Y. Du, General uniqueness results and variation speed for blow-up solutions of elliptic equations, Proc. London Math. Soc. 91(2005) 459-482.
因此,
0
<
ς

hL
:=
lim inf
t→∞
h(t)

hM
:=
lim sup h(t)
t→∞

λ.
我们断定
lim h(t) 存在.
t→∞
反证. 假设hL < hM . 则存在两族数列{tn}, {sn}, n = 1, 2, · · ·, 使得
lim
n→∞
tn
=
lim
n→∞
sn
= ∞,
lim
n→∞
h(tn)
F (t) = b(t)(K(t))−2q/(q−1)Υ0(t), t ≥ T,
0 < u¯ = λb(t)g(u¯)Υ0(t)Υ1(t) ≤ b(t)g(u¯), t > T,
u¯(T ) ≥ ψ(T ), u¯(∞) = 0,
u = ςb(t)g(u)Υ0(t)Υ1(t) ≥ b(t)g(u), t > T, u(T ) ≤ ψ(T ), u(∞) = 0.
1
一类二阶非线性微分方程奇异边值问题唯一整体解的精确渐近行为
1
冯化冰, 张志军
烟台大学数学与信息科学学院,山东、烟台, 264005 E-mail: zhangzj@
摘要: 应用摄动方法并构造新的上下解,在权函数b满足适当的条件下,得到了一类半直线上 二阶非线性微分方程奇异边值问题唯一整体解在无穷远处的精确渐近行为.
当σ > 1时, Ck = 0.Leabharlann 22. 定理的证明
3
本节证明定理1.1. 先做一些预备. 考虑摄动问题
u = b(t)g(u), u(t) > 0, t > T, u(T ) = ψ(T ), u(∞) = 0. 这里,ψ 是问题(1.1)的唯一整体解,T > 0.
u ≥ b(t)g(u), u(t) > 0, t > T, u(T ) ≤ ψ(T ), u(∞) = 0.
(2.3)
引理2.1 在定理1.1 的假设条件下, 存在T > 0 和0 < ς0 < λ0, 使得对任意ς ∈ (0, ς0] 和λ ∈ [λ0, ∞), u¯(t) = λF (t) 和u(t) = ςF (t) 分别是问题(2.1) 的一个上解和下解, 这里,
[2] F. Cˆırstea, V. D. Rˇadulescu, Uniqueness of the blow-up boundary solution of logistic equations with absorbtion, C. R. Acad. Sci. Paris, S´er. I 335 (2002) 447-452.
即u¯ 和u 分别是问题(2.1) 的上下解, 并且(2.5)成立.
引理2.2 在引理2.1中,对任意给定ς ∈ (0, ς0] 和λ ∈ [λ0, ∞),
ςF (t) ≤ ψ(t) ≤ λF (t), ∀t ≥ T.
证明 我们只需证ψ(t) ≤ λF (t), ∀t ≥ T . 反证. 假设存在t0 > T 和λ ∈ [λ0, ∞), 使 得ψ(t0) > λF (t0). 由于ψ和F 在(0, ∞)是严格单调递减的且ψ(T ) ≤ λF (T ),则存在t1 ∈ [T, t0) 使得

h
(t)
F F
(t) (t)
+ 2h
(t)
F F
(t) (t)
+ h(t) =
b(t)g(h(t)F F (t)
(t))
.
再由(2.9), 对任意的n ≥ 1, 得到,
h(tn) ≥ h
(tn)
F F
(tn) (tn)
+ h(tn) =
b(tn
)g(h(tn)F F (tn)
(tn))
,
h(sn) ≤ h
1
2
b 还满足
(B2)
b0
:=
lim
t→0+
Γ(t)Γ (t) [Γ (t)]2

(0, ∞),
则ψ满足
lim
t→0+
ψ(t) Γ(t)
=
b−0 p/(p−1)
p+1 p−1
(p+1)/(p−1)
C0−1/(p−1),
这里,R > 0 是某一个常数,
mθ0 < Υ0(t) < M θ0, ∀t ≥ T

m c0ξq
< Υ1(t) <
M c0ξq
,
∀t ≥ T.
3
4
设λ 和ς 是满足下面条件的正常数
λ ≥ λ0 := max{ψ(T )(K(T ))2/(q−1),
M 2θ0 c0
1/(q−1)
},

ds
< ∞, a > 0, G(s) =
s
g(τ )dτ, s ≥ 0.
a G(s)
0
(g3) g∗(s) = g(s)/s 在(0, ∞)单调递增;
(g4) 存在q > 1 使得
c0
:=
lim
s→0+
g(s) sq

(0, ∞);
b 满足
(b1) b ∈ C1(0, ∞), b 在(0, ∞)上是正的单调递增函数;
矛盾. 从而结论获证.
定理1.1的证明 考虑辅助函数
h(t) :=
ψ(t) F (t)
,
t ∈ [T, ∞).
由引理2.2可知, 对任意给定ς ∈ (0, ς0] 和λ ∈ [λ0, ∞),
ς ≤ h(t) ≤ λ, ∀t ∈ [T, ∞).
(2.8)
4
5
(sn
)
F F
(sn) (sn)
+ h(sn) =
b(sn)g(h(sn)F F (sn)
(sn))
.
在不等式(2.10)和(2.11)中, 令n → ∞, 结合(2.6)和(2.7), 得到
hM

c0hqM θ0
和 hL ≤
c0hqL θ0
.
从而 这与假设hL < hM 矛盾. 因此,
hL = hM =
F (t) = (K(t))−2/(q−1).
(2.4)
而且,
u (t) < 0, u¯ (t) < 0, u (t) > 0, u¯ (t) > 0, t ≥ T.
证明 记
Υ0(t)
=
q
2 −1
q q
+ −
1 1

k
(t)K (t) k2(t)
, t>0

Υ1(t) =
(F (t))q g(ξF (t))
相关文档
最新文档