常微分方程初等解法的研究
- 格式:docx
- 大小:266.43 KB
- 文档页数:21
目录摘要 (I)关键词 (I)Abstract (I)Key words (I)1.前言 (1)2.常微分方程的求解方法 (1)2.1常微分方程变量可分离类型解法 (1)2.1.1直接可分离变量的微分方程 (3)2.1.2可化为变量分离方程 (3)2.2常数变易法 (7)2.2.1一阶线性非齐次微分方程的常数变易法 (7)2.2.2一阶非线性微分方程的常数变易法 (8)2.3积分因子法 (12)3.实例分析说明这几类方法间的联系及优劣 (14)3.1几个重要的变换技巧及实例 (14)3.1.1变为 (14)3.1.2分项组合法组合原则 (15)3.1.3积分因子选择 (15)参考文献 (16)致谢 (17)常微分方程初等解法及其求解技巧摘要常微分方程是微积分学的重要组成部分,广泛用于具体问题的研究中.求解常微分的问题,常常通过变量分离、两边积分,如果是高阶的则通过适当的变量代换,达到降阶的目的来解决问题.本文就是对不同类型的常微分方程的解法及其求解技巧的系统总结:先介绍求解常微分方程的几种初等解法,如变量分离法,常数变易法,积分因子法等,在学习过程中,通过对不同类型的方程求解,揭示常微分方程的求解规律.然后介绍几类方程求解中的变换技巧及规律,并通过实例来分析这几类方法之间的联系及优劣,从而能快速的找到最佳解法.关键词变量分离法常数变易法积分因子变换技巧Elementary Solution and Solving Skills of OrdinaryDifferential EquationAbstractOrdinary differential equations are important components of calculus and used extensively for the studies on specific issues. Ordinary differential equations are often resolved by the means of variable separation and both sides integral. If they are higher-order ones, we can reduce their order by proper variable substitution to solve this problem. This essay aims at concluding systematically the methods of different types of differential equations and its resoling skills. First of all, I’d would like to introduce several basic resolutions of differential equations, such as variable separation, constant threats, points factor, etc. In the process of learning, I’d like to reduce the law of resolving ordinary differential equations by resolving different types of equations. Then, we describe several equations resolutions and for transformation techniques and its laws,and we also analyze the advantages and disadvantages and connections by using the examples of these methods to be able to find the best solution quickly.Key wordsVariable separation; constant threats; points factor; transform techniques1.前言数学发展的历史告诉我们,300年来数学分析是数学的首要分支,而微分方程又是数学分析的心脏,它还是高等分析里大部分思想和理论的根源.人所共知,常微分方程从它产生的那天起, 就是研究自然界变化规律、研究人类社会结构、生态结构和工程技术问题的强有力工具.它的发展历史也是跟整个科学发展史大致同步的.现在,常微分方程在很多学科领域内有着重要的应用,自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性质的研究、化学反应稳定性的研究等.这些问题都可以转化为求常微分方程的解,或者化为研究解的性质的问题.常微分方程具有广泛的社会实践性,无论是在各类学科领域上,还是在实际生产生活中,都有举足轻重的作用.它所涉及范围之广,致使前人对它做了很深入的研究.应用常微分方程理论已经取得了很大的成就,但是,它现有的理论也还远远不能满足需要,还有待进一步的发展,使这门学科的理论更加完善.微分方程是表达自然规律的一种自然的数学语言.它从生产实践与科学技术中产生,而又成为现代科学技术中分析问题与解决问题的一个强有力的工具.人们在探求物质世界某些规律的过程中,一般很难完全依靠实验观测认识到该规律,反而是依照某种规律存在的联系常常容易被我们捕捉到,而这种规律用数学语言表达出来,其结果往往形成一个微分方程,而一旦求出方程的解,其规律则一目了然.所以我们必须能够求出它的解.常微分方程的初等解法,既是常微分方程理论中有自身特色的部分,也与实际问题密切相关;恰当对初等解法进行归类,能正确而又敏捷地判断一个给定的方程属于何种类型,从而能按照所介绍的方法进行分解.总之,常微分方程属于数学分析或基础数学的一个组成部分,在整个数学大厦中占据这重要位置,学好常微分方程基本理论与方法对进一步学习研究数学理论与实际应用均非常重要,因此本文对常微分方程的初等解法进行了简要归纳和分析,主要讨论变量分离方程,非恰当微分方程,线性微分方程,同时结合具体的实例,展示了初等解法在解题过程中的应用及其求解过程中的变换技巧和律.2.常微分方程的求解方法2.1常微分方程变量可分离类型解法定义1 如果一阶微分方程具有形式,则该方程称为可分离变量微分方程.若设,则可将方程化为.即将两个变量分离在等式两端.其特点是:方程的一端只含有的函数与,另一端只含有的函数与.对于该类程,我们通常采用分离变量的方法来处理。
常微分方程的初等解法1.常微分方程的基本概况1.1.定义:自变量﹑未知函数及函数的导数(或微分)组成的关系式,得到的便是微分方程,通过求解微分方程求出未知函数,自变量只有一个的微分方程称为常微分方程。
1.2.研究对象:常微分方程是研究自然科学和社会科学中的事物、物体和现象运动﹑演化和变化规律的最为基本的数学理论和方法。
物理﹑化学﹑生物﹑工程﹑航空﹑航天﹑医学﹑经济和金融领域中的许多原理和规律都可以描述成适当的常微分方程。
如牛顿运动规律、万有引力﹑能量守恒﹑人口发展规律﹑生态总群竞争﹑疾病传染﹑遗传基因变异﹑股票的涨伏趋势﹑利率的浮动﹑市场均衡价格的变化等。
对这些规律的描述﹑认识和分析就归结为对相应的常微分方程的理论和方法不仅广泛应用于自然科学,而且越来越多的应用于社会科学各个领域。
1.3.特点:常微分方程的概念、解法、和其它理论很多,比如,方程和方程组的种类及解法、解的存在性和唯一性、奇解、定性理论等等。
下面就方程解的有关几点简述一下,以了解常微分方程的特点。
求通解在历史上曾作为微分方程的主要目标,一旦求出通解的表达式,就容易从中得到问题所需要的特解。
也可以由通解的表达式,了解对某些参数的依赖情况,便于参数取值适宜,使它对应的解具有所需要的性能,还有助于进行关于解的其他研究。
1.4.应用:现在,常微分方程在很多学科领域内有着重要的应用,自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。
这些问题都可以化为求常微分方程的解,或者化为研究解的性质的问题。
应该说,应用常微分方程理论已经取得了很大的成就,但是,它的现有理论也还远远不能满足需要,还有待于进一步的发展,使这门学科的理论更加完善。
2.一阶的常微分方程的初等解法一阶常微分的初等解法包括变量分离方程与变量变换﹑可以化为变量分离方程的类型﹑线性微分方程与常数变易法﹑恰当微分方程与积分因子,下面我们就具体分析一阶常微分方程的初等解法。
二阶常系数常微分方程的初等解法求解技巧(word版可编辑修改)二阶常系数常微分方程的初等解法求解技巧(word版可编辑修改) 编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(二阶常系数常微分方程的初等解法求解技巧(word版可编辑修改))的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为二阶常系数常微分方程的初等解法求解技巧(word版可编辑修改)的全部内容。
二阶常系数常微分方程的初等解法求解技巧(word版可编辑修改)二阶常系数常微分方程的初等解法求解技巧郑燕,王俊霞太原师范学院数学系,山西晋中,030619摘要:本文总结介绍了三类二阶常系数常微分方程的初等解法求解技巧,分别是:特征根法;常数变易法;比较系数法.同时结合例题进行具体讲解.虽然当今社会关于二阶常微分方程初等解法求解技巧的研究已经获得了很大的成就,但它的已有理论仍然得不到求知者的满足,需要大家进一步发展,使之更加完善。
关键词:二阶常系数齐次线性微分方程;特征根法;常数变易法;比较系数法;二阶常系数非齐次线性微分方程.1。
预备知识(1.1)其中以及f(t)都是连续函数并且区间是a t b。
如果,则方程(1)就变成了(1.2)我们形如方程(1.2)的方程叫做二阶齐次线性微分方程,把方程(1。
1)叫做二阶非齐次线性微分方程.并且把方程(1.1)叫做方程(1.2)对应的齐次线性微分方程。
2.求解方法技巧2.1常数变易法常数变易法是将常数看作是的待定函数,然后求出非齐次线性方程的通解。
求解过程如下:设,是方程(1.2)的基本解组,则(2.1.1)是方程(1。
2)的通解。
将常数看作是t的待定函数,那么方程(2。
常微分方程讲解常微分方程第一章绪论在初等数学中,我们已经学过一些代数方程(如元个一次联立方程),并且用它们解决了一些有趣的应用问题,使我们初步体会到方程论(主要是设未知量、列方程和求解方程的方法)对于解决实际问题的重要性。
在解析几何与微积分中,我们又碰到一类不同的方程——方程的个数少于未知量的个数,也就是通常所说的函数方程。
例如,1) (设是自变量,则是未知函数);2),(设是自变量,则和是两个未知函数)。
这类函数方程与开头所说的代数方程相比,在概念上进了一步——确定自变量与因变量之间的函数关系。
利用这类方程可以解决一类新的问题,例如某些轨迹问题和极值问题等。
本课程所要讲述的方程与刚才说的那种函数方程又不一样,它们除了自变量和未知函数外,还包含了未知函数的导数(即微商)。
例如:1)(是自变量,是未知函数,是未知函数对的导数。
)2)(是自变量,是未知函数,是未知函数对的导数等等)。
这种联系着自变量、未知函数以及未知函数的导数(或微分)的关系式,数学上称之为微分方程。
其中未知函数的导数或微分是不可缺少的。
下面我们通过几个具体的例子,粗略地介绍常微分方程的一些物理背景和方程的建立问题,并讲述一些最基本的概念。
第一节微分方程:某些物理过程的数学模型在这一节中列举几个简单的实际例子,说明怎样从实际问题列成微分方程的问题。
例子虽然简单,但是从中能够简明地诱导出微分方程的一些基本概念,成为进一步探讨其他较复杂问题的借鉴。
掌握好这些例子,会有助于增进我们分析问题的能力。
例1 物体冷却过程的数学模型将某物体放置于空气中,在时刻时,测量得它的温度为,10分钟后测得温度为。
我们要求决定此物体的温度和时间的关系,并计算20分钟后物体的温度。
这里我们假定空气的温度保持为。
解为了解决上述问题,需要了解有关热力学的一些基本规律。
例如,热量总是从温度高的物体向温度低的物体传导的;在一定的温度范围内(其中包括了上述问题的温度在内),一个物体的温度变化速度与这一物体的温度和其所在介质温度的差值成比例。
一阶常微分方程初等解法研究1. 分离变量法:对于形如 dy/dx = f(x)g(y) 的方程,如果可以将变量 x 和 y 分离出来,即可写成 dx/f(x) = dy/g(y),然后对两边同时积分即可得到方程的解。
2. 齐次方程:对于形如 dy/dx = f(y/x) 的方程,我们可以令 v = y/x,然后通过代换和分离变量的方式将其转化为一阶线性方程,进而求解。
3. 线性齐次方程:对于形如 dy/dx + p(x)y = 0 的方程,我们可以通过乘以一个积分因子来将其转化为可分离变量的形式,进而求解。
4. 一阶线性方程:对于形如 dy/dx + p(x)y = q(x) 的方程,可以通过乘以一个积分因子来将其转化为一阶线性常微分方程组的形式,然后通过求解常微分方程组得到原方程的解。
5. 可分离变量的方程:对于形如 dy/dx = f(x)g(y) 的方程,如果可以将变量 x 和 y 分离出来,即可写成 dx/f(x) = dy/g(y),然后对两边同时积分即可得到方程的解。
以上是一阶常微分方程的初等解法研究。
这些方法广泛适用于各种类型的一阶常微分方程,能够通过简单的代数运算和积分求解方程,得到解析解。
但对于一些特殊类型的方程,可能需要借助其他方法求解,或者使用数值方法进行求解。
除了初等解法,还有一些其他的方法可以用于求解一阶常微分方程,如变量替换、常数变易法、特解叠加法等。
这些方法在特定情况下可以简化方程的求解过程,提高求解效率。
此外,对于更高阶的微分方程,可以利用一阶常微分方程的解法来进行逐步求解。
总结起来,一阶常微分方程初等解法的研究可以帮助我们理解微分方程的性质和求解方法,掌握这些解法对于解决实际问题和推导其他微分方程的解法都具有重要意义。
因此,研究一阶常微分方程的初等解法有着广泛的应用价值。
第二章 基本定理我们在第一章主要学习了初等积分法,掌握了几类常微分方程的解法.但是这些解法只适用于某些特殊的类型,很多其它的常微分方程不能用初等解法进行求解.1841年,法国数学家刘维尔(Liouville )证明了里卡蒂(Riccati )方程)0)(()()()(2≠++=x p x r y x q y x p dydx 除了某些特殊的类型外,一般不能用初等积分法求解.例如,很简单的里卡蒂方程22y x dxdy +=就不能用初等积分法求解.自然地,如果一个常微分方程不能用初等积分法求解,那么应该如何处理呢?是否存在解呢?如果存在解,它的解是否唯一呢?解的存在区间是什么呢?初值的微小误差对解有什么影响呢?这些问题在理论的研究和实际应用中,都有着重要的意义.本章将解决这些基本问题. 本章主要介绍解的存在唯一性定理、解的延展定理与比较定理、解对初值的连续依赖性定理以及解对初值的可微性定理,这些定理就回答了我们刚才的疑问,有效的处理解的存在性、唯一性、存在区间、初值对解的影响等问题,为我们使近似解法奠定理论基础,同时这些定理也是常微分方程理论的基础内容,对进一步的学习奠定基础.2.1 解的存在唯一性定理对于一般的常微分方程),(y x f dxdy = (2.1) 如果给出了初始条件00)(y x y =,我们就得到了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy (2.2) 这时,在什么样的条件下,柯西初值问题的解存在且唯一呢?解的存在区间是什么呢?我们有如下的解的存在唯一性定理.2.1.1 存在唯一性定理的叙述定理2.1(存在唯一性定理)如果方程(2.1)的右端函数),(y x f 在闭矩形区域b y y b y a x x a x R +≤≤-+≤≤-00002,:上满足如下条件:(1)在2R 上连续;(2)在2R 上关于变量y 满足李普希兹(Lipschitz )条件,即存在常数N ,使对于2R 上的任何一对点),(y x 和),(x 有不等式:y y N y x f y x f -≤-),(),(则初值问题(2.2)在区间],[0000h x h x +-上存在唯一解00)(),(y x x y ==ϕϕ 其中),(max ),,min(),(0y x f M Mb a h R y x ∈==. 在给出定理2.1的证明之前,我们先对定理2.1的条件和结论做些说明:1、在两个条件中,条件(2),即李普希兹条件比较难于验证,因为李普希兹常数N 难以确定.但是,我们可以将该条件加强,替换为:如果函数),(y x f 在闭矩形区域2R 关于y 的偏导数),(y x f y '存在且有界.这样,可以推出李普希兹条件成立.事实上,因为),(y x f y '有界,故设N y x f y ≤'),(,对2),(),,(R x y x ∈∀,由拉格朗日中值定理得:y y N y y x f y x f y x f y -≤-'=-),(),(),(ξ我们验证),(y x f y '在闭矩形区域2R 上有界也不容易,可以进一步将条件加强为:),(y x f y '在闭矩形区域2R 上连续.由闭区域上连续函数的性质知:),(y x f y '在闭矩形区域2R 上有界,所以李普希兹条件成立.因此,有如下的关系式:),(y x f y '在2R 上连续⇒),(y x f y '在2R 上存在且有界⇒李普希兹条件2、在定理2.1的结论中,解)(x y ϕ=的存在区间为],[0000h x h x +-,其中 ),(max ),,min(),(0y x f M Mb a h R y x ∈==.为什么解的存在区间不是],[00a x a x +-呢?这是因为我们研究问题的范围为闭矩形区域2R ,方程的解)(x y ϕ=不能超出2R 的范围,又因为),(max ),(y x f M Ry x ∈=,所以M y x f M ≤≤-),( 即 M dxdy M ≤≤- 由⎪⎩⎪⎨⎧=-=00)(y x y M dx dy 和⎪⎩⎪⎨⎧==00)(y x y M dx dy 得:001)()(y x x M x y +--=,002)()(y x x M x y +-= 因此)()()(21x y x y x y ≤=≤ϕ,即)(x y ϕ=夹在)(1x y 与)(2x y 之间.又,)(1x y 与)(2x y 在2R 上的存在区间为],[0000h x h x +-,故)(x y ϕ=的存在区间也是],[0000h x h x +-.2.1.2 存在性的证明首先,我们给出柯西初值问题(2.2)的等价转化,即求(2.2)的解)(x y ϕ=,等价于求解积分方程⎰+=xx d y f y y 0))(,(0ξξξ (2.3) 事实上,如果)(x y ϕ=是初值问题(2.2)的解,即有))(,()(x x f x ϕϕ='且00)(y x =ϕ从0x 到x 积分得:⎰+=xx d f y x 0))(,()(0ξξϕξϕ 即)(x y ϕ=是积分问题(2.3)的解.反过来,如果)(x y ϕ=是积分问题(2.3)的解,即有⎰+=xx d f y x 0))(,()(0ξξϕξϕ 则00)(y x =ϕ且))(,()(x x f x ϕϕ='即)(x y ϕ=是初值问题(2.2)的解.经过等价转化,我们将初值问题(2.2)的求解,转化为积分问题(2.3)的求解.下面用皮卡(Picard )逐次逼近来证明积分问题(2.3)的解的存在性,分为三个步骤:1、构造近似函数列{})(x n ϕ任取一个满足初值条件00)(y x y =的函数)(0x y ϕ=作为首项(初始项),并要求在2R 上的存在区间为:],[0000h x h x +-,简单起见,取00)(y x =ϕ,将它代入方程(2.3)的右端,所得到的函数用)(1x ϕ表示,并称为一次近似,即⎰+=xx d f y x 0))(,()(001ξξϕξϕ 再将)(1x ϕ代入方程(2.3)的右端就得到二次近似⎰+=xx d f y x 0))(,()(102ξξϕξϕ 序行此法,可以得到n 次近似⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ 为了保证上述的逐次逼近过程可以一直进行下去,必须有2))(,(R x x n ∈ϕ,即当],[0000h x h x x +-∈时,有,2,1)(0=≤-n b y x n ϕ 下面用数学归纳法证明b y x n ≤-0)(ϕ.显然,当],[0000h x h x x +-∈时,有b y y y x ≤=-=-0)(0000ϕ假设,当],[0000h x h x x +-∈时,有b y x n ≤--01)(ϕ,那么,对于)(x n ϕ有⎰-=-xx n n d f y x 0))(,()(10ξξϕξϕ 从而有b Mb M Mh x x M d f y x xx n n =≤≤-≤≤-⎰-00100))(,()(ξξϕξϕ 由数学归纳法知,当],[0000h x h x x +-∈时,有,2,1)(0=≤-n b y x n ϕ这样,我们就可以得到一个近似函数列{})(x n ϕ.2、证明近似函数列{})(x n ϕ在区间],[0000h x h x +-上一致收敛.由于无法得到{})(x n ϕ的通项公式,只知道首项和递推关系式,直接证明函数列{})(x n ϕ的收敛性比较困难,为此我们构造函数项级数+-++-+-)]()([)]()([)(1010x x x x x n n ϕϕϕϕϕ (2.4) 它的部分和是)()]()([)]()([)()(10101x x x x x x x S n n n n ϕϕϕϕϕϕ=-++-+=-+因此,证明{})(x n ϕ的收敛性转化为证明级数(2.4)的收敛性,下面我们证明级数(2.4)在区间],[0000h x h x +-上一致收敛.首先研究级数(2.4)的通项)(x n μ⎰=-xx d f x x 0))(,()()(001ξξϕξϕϕ 即⎰=-xx d y f y x 0),()(001ξξϕ 所以00010),()(x x M d y f y x x x -≤≤-⎰ξξϕ 因为⎰+=x x d f y x 0))(,()(001ξξϕξϕ,⎰+=x x d f y x 0))(,()(102ξξϕξϕ,所以 ⎰-≤-x x d f f x x 0))(,())(,()()(0112ξξϕξξϕξϕϕ由李普希兹条件,得 !2)()()()(200011200x x MN d x MN d N x x x x x x -=-≤-≤-⎰⎰ξξξξϕξϕϕϕ 下面用数学归纳法证明!)()(011n x x MN x x nn n n -≤---ϕϕ 显然,2,1=n 的时候,不等式成立(上面已经给出), 假设!)()(011n x x MN x x n n n n -≤---ϕϕ成立,那么对于1+n 的情形有 )!1(!)()())(,())(,()()(100111000+-=-≤-≤-≤-+--+⎰⎰⎰n x x MN d n x MN d N d f f x x n n x x n n xx n n x x n n n n ξξξξϕξϕξξϕξξϕξϕϕ由数学归纳法知,对一切自然数n ,均有!)()(011n x x MNx x nn n n -≤---ϕϕ 又00h x x ≤-,所以级数(2.4)的通项满足: !)(011n h MN v x n n n n -+=≤μ ( ,2,1=n ) 利用比式判别法,可知以n v 为通项的级数收敛,从而以)(x n μ为通项的级数(2.4)绝对收敛且一致收敛.又,每一个)(x n μ是连续的,所以级数(2.4)的和函数也是连续的,记为)(x ϕ,其存在区间也是],[0000h x h x +-.因此函数列{})(x n ϕ就收敛于)(x ϕ.3、证明)(lim )(x x n n ϕϕ∞→=是积分问题(2.3)的解,从而也是初值问题(2.2)的解.在⎰-+=x x n n d f y x 0))(,()(10ξξϕξϕ两端取极限,得到 ⎰-∞→∞→+=xx n n n n d f y x 0))(,(lim )(lim 10ξξϕξϕ 即⎰+=xx d f y x 0))(,()(0ξξϕξϕ 所以)(x ϕ是积分问题(2.3)的解,从而也是初值问题(2.2)的解.2.1.3 唯一性的证明下面我们证明解的唯一性.在证明唯一性之前,先介绍一个重要的不等式,即贝尔曼(Bellman )不等式.贝尔曼引理 设)(x y 为区间],[b a 上的非负连续函数,b x a ≤≤0.若存在,0≥δ 0≥k ,使得)(x y 满足不等式],[,)()(0b a x d y k x y xx ∈+≤⎰ττδ (2.5) 则有],[,)(0b a x e x y x x k ∈≤-δ证明 仅证明0x x ≥的情形,0x x ≤的情形类似.令)(x y 的原函数为⎰=xx d y x R 0)()(ττ,代入(2.5)得 δ≤-')()(x kR x R两边同时乘以积分因子)(0x x k e --,得)()(00)]()([x x k x x k e x kR x R e ----≤-'δ从0x 到x 积分得)()(00)(x x k x x k e e x kR -----≤δδ即)(0)(x x k e x kR -≤+δδ 由(2.5)知,)()(x kR x y +≤δ,所以],[,)(0b a x e x y x x k ∈≤-δ下面证明积分问题(2.3)的解的唯一性.假设积分问题(2.3)有两个解)(1x y 和)(2x y ,我们只需要证明:)(1x y )(2x y ≡,],[0000h x h x x +-∈事实上,因为⎰+=x x d y f y x y 0))(,()(101ξξξ,⎰+=xx d y f y x y 0))(,()(202ξξξ 所以有⎰-≤-xx d y f y f x y x y 0))(,())(,()()(2121ξξξξξ由李普希兹条件知⎰-≤-xx d y y N x y x y 0)()()()(2121ξξξ 令N k x y x y x y ==-=,0,)()()(21δ,由贝尔曼引理可知,0)(=x y ,即)(1x y )(2x y ≡. 这样,我们就完成了解的存在性与唯一性的证明.2.1.4 三点说明为了更好的理解和掌握解的存在唯一性定理,我们对该定理再做三点说明.1、在存在性的证明过程中,我们利用逐次逼近法构造了近似函数列{})(x n ϕ,其中首项为:00)(y x =ϕ,递推关系式为:⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ.该方法实际上给出了我们一种求初值问题(2.2)的近似解的方法,当用n 次近似解逼近精确解时,需要给出它的误差估计.事实上,有∑∑∞+=∞=+-≤-≤-101!)()()()(n k k k nk k k n k x x N N M x x x x ϕϕϕϕ 0)!1()(!)!1()(!10001010Nh n k k k n n k k k e n Nh N M k h N n Nh N M k h N N M +=+<≤+∞=+∞+=∑∑ 2、如果方程(2.1)是线性方程,即)()(x q y x p dxdy +-= 其中)(x p 和)(x q 在区间],[b a 上连续,这时,初值问题(2.2)在带型区域+∞<<-∞≤≤y b x a R ,:2满足定理2.1的条件.事实上,)()(),(x q y x p y x f +-=在2R 上连续,而且)(),(x p y x f y -='在2R 上也连续,所以),(y x f 关于变量y 满足李普希兹条件.这时,初值问题(2.2)的解存在且唯一,存在区间为],[b a .3、定理2.1中的李普希兹条件是保证解唯一的充分条件,那么这个条件是不是必要条件呢?回答是否定的,即李普希兹条件是解唯一的充分非必要条件.下面我们给出一个例子来说明李普希兹条件是解唯一的非必要条件,也就是说,即使李普希兹条件不成立,初值问题(2.2)的解也可能是唯一的.例1 试证方程0,ln ,0≠=⎩⎨⎧=y y y y dx dy 经过xOy 平面上任一点的解都是唯一的.证明 由00,ln ,0≠=⎩⎨⎧=y y y y dx dy 可得:0=y 或x Ce e y ±=. 任给xOy 平面上的一个点),(00y x ,只会对应0=y 或xCe e y ±=中的一个解,也就是说,过xOy 平面上任一点的解都是唯一的.但是,我们有0ln ln )0,(),(-==-y y y y x f y x f 因为+∞=→y y ln lim 0,所以找不到0>N ,使得 0)0,(),(-≤-y N x f y x f从而方程右端函数在0=y 的任何邻域上不满足李普希兹条件,但是初值问题(2.2)的解却是唯一的,这说明李普希兹条件是非必要条件.习 题 2.11.试判断方程y x dx dy tan =在区域 (1)π≤≤≤≤-y x R 0,11:1;(2)44,11:2ππ≤≤-≤≤-y x R上是否满足定理2.1的条件?2.讨论方程3123y dx dy =在怎样的区域中满足定理2.1的条件.并求通过)0,0(的一切解.3.试用逐次逼近法求方程2y x dxdy -=满足初值条件0)0(=y 的近似解: )(),(),(),(3210x x x x ϕϕϕϕ并在闭矩形区域11,11:2≤≤-≤≤-y x R 给出三次近似的误差估计.4.利用逐次逼近法求方程22x y dxdy -=适合初值条件1)0(=y 的近似解: )(),(),(210x x x ϕϕϕ并在闭矩形区域111,11:2≤-≤-≤≤-y x R 给出二次近似的误差估计.5.试证明定理2.1中的n 次近似解)(x n ϕ与精确解)(x ϕ有如下的误差估计式:10)!1()()(+-+≤-n n n x x n MN x x ϕϕ 6.在条形区域+∞<≤≤y b x a ,内,假设方程(2.1)的所有解都唯一,对其中任意两个解)(),(21x y x y ,如果有)()(0201x y x y <,则必有b x x x y x y ≤≤<021),()(.7.讨论方程323y dx dy = 解的唯一性.2.2 延展定理和比较定理由解的存在唯一性定理,我们知道,初值问题(2.2)的解在满足一定条件的情况下存在且唯一,但是解的存在区间不是],[00a x a x +-,而是],[0000h x h x +- 其中),(max ),,min(),(0y x f M Mb a h R y x ∈==.如果M 比较大的话,则解的存在区间就非常小,这对我们研究解的性质产生了很大的局限性,只能在很小的范围内有解,当x 超出这个范围时,解的情况就不清楚了.为了解决这个问题,我们有下面的延展定理.2.2.1 延展定理定理2.2(延展定理)如果方程(2.1)的右端函数在区域R R D ⨯⊂上连续,且关于变量y 满足局部的李普希兹条件,即对于D 内的任一闭矩形区域都满足李普希兹条件,则对任何一点D y x ∈),(00,初值问题(2.2)的解)(x y ϕ=可以向左右无限延展,直到))(,(x x ϕ任意接近区域D 的边界.在给出定理的证明之前,先对“))(,(x x ϕ任意接近区域D 的边界”进行说明.当区域D 有界时,积分曲线向左右延展可以任意接近;当区域D 无界时,积分曲线向左、右延展,或者任意接近区域D 的边界(边界存在的话),或者无限远离坐标原点.证明 首先证明区域D 有界的情形.设区域D 的边界为D D L -=(D 为D 的闭包).对于任意给定的正数ε,记L 的ε邻域为εU ,记L 的2ε邻域为2εU ,记L 的4ε邻域为4εU .则集合22εεU D D -=为闭集,且D D ⊂2ε,所以2εD 有界. 只要证明积分曲线可以到达2εD 的边界2εL ,由ε的任意性知,积分曲线就可以任意接近区域D 的边界L .事实上,以2εD 中的任意一点为中心,以4ε为半径的闭圆区域均包含在区域D 的内部.且在闭区域44εεU D D -=之内.从而,以2εD 中的任意一点为中心,以4221ε=a 为边长的正方形也在闭区域4εD 之内.记 ),(max 4),(1y x f M D y x ε∈= 则过2εD 的任意一点),(**y x 的积分曲线,必至少可在区间],[**h x h x +-上存在,其中)82,82min(),min(1111M M a a h εε==. 于是,过点),(00y x 的积分曲线)(x y ϕ=每向左或向右延展一次,其存在区间就伸长一个确定的正数h ,由于2εD 有界,)(x y ϕ=经过有限次延展后一定可以达到2εD的边界2εL .于是也就可以任意接近区域D 的边界L .其次考虑区域D 为无界的情形.这时,我们可以用闭圆区域,2,1},),{(222=≤+=n n y x y x S n与区域D 取交集,令n n S D D =,则 ∞==1n n D D .由于n D 为有界的区域,根据前面的证明,我们可知,过n D 内任一点的积分曲线能够任意接近n D 的边界.因此,过点),(00y x 的积分曲线)(x y ϕ=可以无限接近区域D 的边界.延展定理的证明,关键是第一步证明,也就是区域D 有界的时候,过点),(00y x 的积分曲线)(x y ϕ=向左向右延展的时候,一定要做等速延展,即延展步幅h 是不变的. 例1 试讨论方程2y dxdy=通过点)1,1(的解和通过点)1,3(-的解的存在区间. 解 该题目中研究问题的区域D 为整个坐标平面xOy .方程右端函数满足延展定理的条件.由2y dxdy=可以解得方程的通解为 xC y -=1代入1)1(=y 得:2=C .故通过点)1,1(的解为xy -=21 它可以向左无限延展,而当-→2x 时,+∞→y ,所以通过点)1,1(的解xy -=21的存在区间为)2,(-∞.代入1)3(-=y 得:2=C .故通过点)1,3(-的解为xy -=21它可以向右无限延展,而当+→2x 时,-∞→y ,所以通过点)1,3(-的解xy -=21的存在区间为),2(+∞.这个例子说明,尽管),(y x f 在整个坐标平面上满足延展定理的条件,解上的点))(,(x x ϕ也能无限接近区域D 的边界,但是延展的方向却不一定是无限向右和向左,可能是向上或向下,从而导致解的存在区间不是),(+∞-∞. 例2 试证明:对任意的0x 及满足条件100<<y 的0y ,方程221)1(y x y y dx dy ++-=的满足条件00)(y x y =的解)(x y y =在),(+∞-∞上存在.证明:令221)1(),(y x y y y x f ++-=,则222222)1(122),(y x x y y x y y x f y ++--++=' 显然),(),,(y x f y x f y '在xOy 平面上连续,满足解的存在唯一性条件及延展定理的条件,而1,0==y y 是),(y x f dxdy=的解, 因此,满足00)(y x y =,100<<y 的解存在,而且可以无限延展到xOy 平面的边界,且不能穿过1,0==y y ,故只能向左右无限延展,所以,)(x y y =在),(+∞-∞上存在.该例题说明,),(y x f 在整个坐标平面上满足延展定理的条件,当方程的解不能穿过1,0==y y 时,它就不能向上向下无限延展了,只能向左、向右延展,所以解的存在区间就是),(+∞-∞.在这里,1,0==y y 控制了解的延展方向,使它按照我们的要求进行延展,因此就有了下面的比较定理. 2.2.2 比较定理我们在使用延展定理的时候,通常会和比较定理配合使用,从而起到控制延展方向的作用.下面介绍一下比较定理.我们在考察方程(2.1)),(y x f dxdy=时,通常将右端函数),(y x f 进行放缩的处理,比如),(),(),(21y x F y x f y x F <<这时,我们可以同时考察),(1y x F dx dy =和),(2y x F dxdy = 我们有如下的比较定理:定理2.3 (第一比较定理)设定义在某个区域D 上的函数),(y x f ,),(1y x F 和),(2y x F 满足条件:(1)在D 满足解的存在唯一性定理及延展定理的条件,即在D 上连续,在D 上关于变量y 满足李普希兹条件;(2)在D 上有不等式),(),(),(21y x F y x f y x F <<设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ,⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=,)(1x y Φ=和)(2x y Φ=,则在它们的共同存在区间上有下列不等式:021),()()(x x x x x >Φ<<Φϕ 021),()()(x x x x x <Φ>>Φϕ证明 仅证当0x x >时,)()(2x x Φ<ϕ,其它的情形相类似. 由比较定理的条件(1),初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解在0x 的某一邻域内存在且唯一,分别记为)(x y ϕ=和)(2x y Φ=,它们满足0020)()(y x x =Φ=ϕ令)()()(2x x x h ϕ-Φ=,则0)()()(0020=-Φ=x x x h ϕ且0))(,())(,()()()(0002020020>-Φ='-Φ'='x x f x x F x x x h ϕϕ所以函数)(x h 在0x 的某一右邻域内是严格单调增加的.如果在0x x >时,0)(>x h 不是总成立,则至少存在一点01x x >,使得0)(1=x h ,且当10x x x <<时,0)(>x h ,因此在点1x 的左导数0)0(1≤-'x h ,这与0))(,())(,()()()(1112121121>-Φ='-Φ'='x x f x x F x x x h ϕϕ矛盾.因此当0x x >时,0)(>x h 总成立,即)()(2x x Φ<ϕ.比较定理的应用,关键是),(1y x F 和),(2y x F 的选取,因为初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解)(x y ϕ=的存在区间的延展,受到)(1x y Φ=和)(2x y Φ=的控制,即)(x y ϕ=夹在)(1x y Φ=和)(2x y Φ=之间.因此,我们必须能确定出)(1x y Φ=和)(2x y Φ=的存在区间,这就是我们选取),(1y x F 和),(2y x F 的标准,即⎪⎩⎪⎨⎧==001)(),(y x y y x F dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解)(1x y Φ=和)(2x y Φ=必须能够求得. 下面我们给出第二比较定理.定理2.4 (第二比较定理)设定义在某个区域D 上的函数),(y x f ,),(1y x F 和),(2y x F 满足条件:(1)在D 满足解的存在唯一性定理及延展定理的条件,即在D 上连续,在D 上关于变量y 满足李普希兹条件;(2)在D 上有不等式),(),(),(21y x F y x f y x F ≤≤设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ,⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=,)(1x y Φ=和)(2x y Φ=,则在它们的共同存在区间上有下列不等式:021),()()(x x x x x >Φ≤≤Φϕ 021),()()(x x x x x <Φ≥≥Φϕ习 题 2.21.设方程为),()(22y x f a y dxdy-= 假设),(y x f 及),(y x f y '在xOy 平面上连续,试证明:对于任意的0x 及a y <0,方程满足00)(y x y =的解都在),(+∞-∞上存在.2.指出方程2)1(2xy e y dxdy -=的每一个解的最大存在区间,以及当x 趋于这个区间的右端点时解的极限.3.讨论方程xx dx dy 1cos 12-= 解的存在区间.4.设),(y x f 在整个平面上连续有界,对y 有连续偏导数,试证明方程),(y x f dxdy=的任一解)(x y ϕ=在区间+∞<<∞-x 上有定义. 5.讨论方程212-=y dx dy 的通过点)0,0(的解,以及通过点)3,2(ln -的解的存在区间.6.在方程)(y f dxdy=中,如果)(y f 在),(+∞-∞上连续可微,且 )0(0)(≠<y y yf ,求证方程满足00)(y x y =的解)(x y 在区间),[0+∞x 上存在,且有0)(lim =+∞→x y x .2.3 解对初值的连续依赖性定理和解对初值的可微性定理通过前两节的存在唯一性定理和延展定理,加上比较定理,我们知道了初值问题(2.2)在什么样的条件下,解是存在的,是唯一的,而且存在区间比较小的时候,通过延展定理和比较定理可以将解的存在区间变大,从而在实际问题中可以达到我们的要求.但是,在实际问题中,还有一个问题需要解决,那就是误差问题.我们的初始条件00)(y x y =如果产生了微小的偏差,这个偏差对我们的初值问题(2.2)的解)(x y ϕ=会有什么影响呢?下面我们来解决这个问题. 我们在研究初值问题(2.2)的时候,习惯上把0x 和0y 当作常数来看待,这样初值问题(2.2)的解)(x y ϕ=被看作x 的函数.实际上,如果0x ,0y 变化,初值问题(2.2)的解)(x y ϕ=也会发生变化.例如方程xydx dy = 经过点),(00y x 的解为x x y y 0=,可以看作00,,y x x 的函数.对于一般的情形,初值问题(2.2)的解也可以看作00,,y x x 的函数,记为),,(00y x x y ϕ=,代入00)(y x y = 得:0000),,(y y x x =ϕ.如果我们的初始条件00)(y x y =发生了微小的误差,变为了**0)(y x y =,初值问题(2.2)的解也变化不大的话,称解连续依赖于初值.下面我们给出连续依赖性的严格定义.定义2.1 设初值问题⎪⎩⎪⎨⎧==**0)(),(y x y y x f dxdy的解),,(*0*0y x x y ϕ=在区间],[b a 上存在,如果对于任意给定的正数ε,存在正数δ (δ的选取与,ε**0,y x 有关),使得对于满足δδ<-<-*00*00,y y x x (2.2)的解),,(00y x x y ϕ=都在],[b a 上存在,且有],,[,),,(),,(*0*000b a x y x x y x x ∈<-εϕϕ则称初值问题(2.2)的解),,(00y x x y ϕ=在点),(*0*0y x 连续依赖于初值,0x 0y .定理2.4 (解对初值的连续依赖性定理)设),(y x f 在区域D 内连续,且关于变量y 满足李普希兹条件.如果D y x ∈),(*0*0,初值问题(2.2)有解),,(*0*0y x x y ϕ=,且当b x a ≤≤时,D y x x x ∈)),,(,(*0*0ϕ,则对任意的正数ε,存在0>δ,使对于满δδ<-<-*00*00,y y x x的任意),(00y x ,初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解),,(00y x x y ϕ=也在区间],[b a 上存在,且有εϕϕ<-),,(),,(*0*000y x x y x x证明 对于任意给定的正数ε,取εδ<<10,使得闭区域}),,(,),{(1*0*0δϕ≤-≤≤=y x x y b x a y x U整个含在区域D 内,这是可以做到的,因为区域D 是开区域,且当b x a ≤≤时,D y x x x ∈)),,(,(*0*0ϕ,所以,只要1δ的选取足够小,以曲线),,(*0*0y x x y ϕ=为中线,宽度为12δ的带形开区域U 就整个包含在区域D 内, 选取δ满足)(110a b N e M--+<<δδ其中N 为李普希兹常数,),(max ),(y x f M Uy x ∈=,同时还要求δ的选取,必须保证闭正方形δδ≤-≤-*0*02,:y y x x R含于带形开区域U 内.由存在唯一性定理知,对于任一200),(R y x ∈,初值问题(2.2)在0x 的某邻域上存在唯一解),,(00y x x y ϕ=,而且),,(00y x x y ϕ=在0x 的该邻域上可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(000000⎰+=而),,(*0*0y x x y ϕ=可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(*0*0*0*0*0*⎰+=对上述两式做差得:ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx x x )),,(,()),,(,(),,(),,(*0*000*00*0*000*⎰⎰-+-=-ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx xx )),,(,()),,(,(),,(),,(*0*000*00*0*000*0⎰⎰-+-≤-ττϕτττϕττϕτd y x f d y x f y x f y y x x xx |)),,(,(||)),,(,()),,(,(|0000*0*0*00**0⎰⎰+-+-≤δττϕττϕτδM d y x f y x f xx +-+≤⎰|)),,(,()),,(,(|00*0*0*0ττϕτϕδd y x y x N M xx |),,(),,(|)1(00*0*0*0-++≤⎰由贝尔曼引理,得εδδδϕϕ<<+≤+≤---1)(*0*000)1()1(),,(),,(*a b N x x N e M e M y x x y x x因此,只要在),,(00y x x y ϕ=有定义的区间上,就有εϕϕ<-),,(),,(*0*000y x x y x x .下面我们证明:),,(00y x x y ϕ=在区间],[b a 上有定义.事实上,因为εϕϕ<-),,(),,(*0*000y x x y x x即解),,(00y x x y ϕ=夹在εϕ+=),,(*0*0y x x y 和εϕ-=),,(*0*0y x x y 之间,而且,初值问题(2.2)满足延展定理的条件,所以,解),,(00y x x y ϕ=可以向左向右无限延展,直到无限接近区域D 的边界,于是,它在延展的时候,必须由直线a x =和直线b x =穿出区域U ,从而),,(00y x x y ϕ=在区间],[b a 上有定义.解对初值的连续依赖性说明,初值),(00y x 无法准确得到,但是我们能得到测量数据),(*0*0y x ,只要误差比较小,即δδ<-<-*00*00,y y x x .我们就可以用),(*0*0y x 代替),(00y x 去计算,得到初值问题的解),,(*0*0y x x y ϕ=,这个解可以非常接近真实解),,(00y x x y ϕ=,即εϕϕ<-),,(),,(*0*000y x x y x x .同理,如果方程的右端函数),(y x f 不能准确得到,只能得到),(y x f 的近似函数),(~y x f ,即)),((,),(),(~D y x y x f y x f ∈<-δ我们就可以用),(~y x f 代替),(y x f 去计算,得到初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 的解),,(00~y x x y ϕ=,那么),,(00~y x x y ϕ=能否代替),,(00y x x y ϕ=呢?我们有下面的解的连续依赖性定理.定理2.5 (解对被积函数的连续依赖性定理)在区域D 上,),(y x f 和),(~y x f 都连续,而且关于变量y 满足李普希兹条件, 若初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 在b x a ≤≤上有解),,(00~y x x y ϕ=,则对任意给定的正数ε,存在0>δ,只要),(y x f 满足)),((,),(),(~D y x y x f y x f ∈<-δ则初值问题(2.2)的解),,(00y x x y ϕ=在b x a ≤≤上存在,且有εϕϕ<-),,(),,(00~00y x x y x x .证明 由解的存在唯一性定理知,初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 的解),,(00~y x x y ϕ=存在,设其存在区间为],[b a ,且有⎰+=xx d y x f y y x x 0))],,(,([),,(00~~000~ξξϕξϕ而初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解),,(00y x x y ϕ=也存在,且可以表示为⎰+=xx d y x f y y x x 0))],,(,([),,(00000ξξϕξϕ则⎰⎰-=-xx xx d y x f d y x f y x x y x x 0))],,(,([))],,(,([),,(),,(0000~~0000~ξξϕξξξϕξϕϕ从而有⎰-≤-xx d y x f y x f y x x y x x 0|)),,(,()),,(,(|),,(),,(0000~~0000~ξξϕξξϕξϕϕ⎰-+-=xx d y x f y x f y x f y x f 0|)),,(,()),,(,()),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ ⎰-+-≤xx d y x f y x f y x f y x f 0|)),,(,()),,(,(||)),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ⎰+-≤xx d y x y x N 0)|),,(),,((|0000~ξδξϕξϕ ⎰-+-≤xx d y x y x N a b 0|),,(),,(|)(0000~ξξϕξϕδ由贝尔曼引理,得)(0000~)(),,(),,(a b N e a b y x x y x x --≤-δϕϕ取)(a b N e ab ---<εδ,则εϕϕ<-),,(),,(0000~y x x y x x .且解),,(00y x x y ϕ=在b x a ≤≤上存在. 例1 考虑方程,ln ,0≠=⎩⎨⎧-=y y y y dx dy 解的情况.解 显然1,1,0-===y y y 是方程的解,当1,1,0-≠≠≠y y y 时,有y y dxdyln -= 这时解得上半平面的通解为x Ce e y -=,下半平面的通解为xCe e y --=.可以看到,对于Ox 轴上的初值)0,(0x ,在任意有限闭区间上解对初值连续依赖,但是,在),0[+∞上,无论),(00y x ,00≠y 如何接近)0,(0x ,只要x 充分大,过),(00y x 的积分曲线就不能与过)0,(0x 的积分曲线(即0=y )任意接近了.这个例子说明,解在有限闭区间上对初值连续依赖,不能推广到无限区间,即,在无限区间上解对初值的连续依赖定理就不成立了.我们有时不仅要求解对初值连续依赖,而且还要知道解),,(00y x x y ϕ=对初值00,y x 的偏导数00,y x ∂∂∂∂ϕϕ是否存在.下面给出解对初值的可微性定理. 定理2.6 (解对初值的可微性定理)如果函数),(y x f 以及),(y x f y '在区域D 内连续,则初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解),,(00y x x y ϕ=在它有定义的区间上有连续偏导数00,y x ∂∂∂∂ϕϕ.并且有 ⎰-=∂∂'x x y d y x f e y x f x y x x 000)),,(,(00000),(),,(ττϕτϕ 及⎰=∂∂'xx y d y x f e y y x x 000)),,(,(000),,(ττϕτϕ 习 题 2.31.若函数),(y x f ,),(y x R 在区域D 内连续且满足李普希兹条件,设初值问题⎪⎩⎪⎨⎧=+=*0*0)(),(),(y x y y x R y x f dx dy 的解为),,(*0*0~y x x y ϕ=,存在区间为],[b a .对任意的正数ε,存在0>δ,使对于满足)),((,),(D y x y x R ∈<δ的),(y x R ,以及满足δδ<-<-*00*00,y y x x的任意),(00y x ,初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解),,(00y x x y ϕ=也在区间],[b a 上存在,且有εϕϕ<-),,(),,(*0*0~00y x x y x x 2.已知方程)sin(xy dxdy = 试求0000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x x y x x y 和0000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x y y x x y 3.设),,(00y x x ϕ是初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解,试证明0),(),,(),,(00000000=∂∂+∂∂y x f y y x x x y x x ϕϕ 2.4 欧拉折线法在第一章,我们介绍了方程的初等解法,即用微积分的知识求得常微分方程的函数解.但是绝大多数的方程不能用初等方法求解,在第二章的前三节中,我们给出了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 在什么样的条件下,解存在且唯一;在什么条件下,解的存在区间可以延展;在什么条件下连续依赖于初值;在什么条件下,解对初值是可微的.有了这些准备,我们就可以研究柯西初值问题的近似解.下面我们介绍求近似解的方法,欧拉折线法.假定函数),(y x f 在区域:+∞<<-∞≤≤y b x a ,上连续,且关于变量y 满足李普希兹条件,求柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 在区间],[0b x 上的近似解,我们采用的方法是:(1)等分区间],[0b x ,分点为n k kh x x k ,,1,0,0 =+=;小区间长度nx b h 0-=, (2)第一个小区间上用切线段逼近曲线:))(,(0000x x y x f y y -+=,(3)求出1x 所对应的纵坐标))(,(010001x x y x f y y -+=,(4)依次重复(2),(3)得到每个小区间上的线段,从而得到欧拉折线. 这样,我们就用欧拉折线作为柯西初值问题在区间],[0b x 近似解.欧拉折线法的前提是:柯西初值问题的解存在且唯一,而且解的存在区间是],[0b x .例1试用欧拉折线法,取步长1.0=h ,求初值问题⎪⎩⎪⎨⎧=+=1)1(22y y x dx dy 的解在4.1=x 时的近似值.解 令22),(y x y x f +=,2)1,1(=f ,这时12-=x y ,代入1.11=x 得:2.11=y ,65.2)2.1,1.1(=f ,这时2.1)1.1(65.2+-=x y , 代入2.12=x 得:465.12=y ,586225.3)465.1,2.1(=f ,这时465.1)2.1(586225.3+-=x y , 代入3.13=x 得:8236225.13=y ,0155990225.5)8236225.1,3.1(=f ,这时8236225.1)3.1(0155990225.5+-=x y ,代入4.14=x 得:53251824022.24=y 习 题 2.41. 试用欧拉折线法,取步长1.0=h ,求初值问题⎪⎩⎪⎨⎧=-=0)0(22y y x dx dy 的解在5.1=x 时的近似值.2.试用欧拉折线法,取步长1.0=h ,求初值问题 ⎪⎩⎪⎨⎧=+=2)1(22y y x dx dy 在区间]4.1,1[上的近似解.。
黎卡提方程的初等解法摘要:常微分方程有着深刻而生动的实际背景,它从生产实践与科学技术中产生,而又成为现代科学技术中分析问题与解决问题的一个强有力的工具,例如化学,生物学,电子技术等等都提出大量的微分方程问题,那么就需要探讨微分方程的求解问题,本文介绍了著名的黎卡提方程的,给出了黎卡提方程可用初等积分法求解的一些充分条件,使得黎卡提方程在满足一定条件下可以用初等解法求解,并给出一些特殊类型黎卡提方程的通解表示,最后举例对一些具体的黎卡提方程进行求解,及微分方程的应用举例。
关键词:黎卡提方程,变量方程,伯努利方程,线性方程0. 引言常微分方程是数学的一个重要分支,也是偏微分方程,变分法,控制论等数学分支的基础。
微分方程的理论和方法从17世纪末开始发展起来,很快成了研究自然现象的强有力的工具。
在17~18世纪,在力学,天文,物理和技术科学中,就已借助微分方程取得了巨大成就。
微分方程的首要问题是如何给定一个方程的通解或特解。
到目前为止,人们已经对许多微分方程得出了求解的一般方法。
例如一阶微分方程中的变量分离方程、线性方程等等。
求一个方程的解最自然的想法是用初等解法求解,即把微分方程的求解问题化为积分问题,但这是不容易做到的,能用初等解法求解的微分方程为数很少,绝大部分的微分方程都无法求出通解,黎卡提方程便是其中的一个。
意大利数学家黎卡提于 1724 年给出了它的特殊形式,后来引起许多学者的研究。
达朗贝尔在 1763年给出了它的一般形式,并首先称之为“黎卡提方程”;黎卡提方程'2y =p(x)y (x)y+r(x)q + 不同于线性微分方程'y =(x)y+r(x)q 之处是还多含一项2p(x)y ,但这就大大地改变了解的性质,即初等可积性丧失了,但在特殊情况下仍旧可以利用初等积分法进行求解。
文献[2]和[3]汇集了很多可积方程和可积性成果;60年代以来,《美国数学月刊》上又连续发表了多篇关于这方面的论文;近年来《数学通报》也发表了多篇关于这一内容的文章,如[2][4]及[5]。
数学毕业(学位)论文题目汇总一、数学理论1。
试论导函数、原函数的一些性质。
ﻫ2。
有界闭区域中连续函数的性质讨论及一些推广。
ﻫ3。
数学中一些有用的不等式及推广.4。
函数的概念及推广.ﻫ5。
构造函数证明问题的妙想。
6.对指数函数的认识。
ﻫ7。
泰勒公式及其在解题中的应用。
8。
导数的作用。
9。
Hilbert空间的一些性质。
ﻫ10。
Banach空间的一些性质。
ﻫ11。
线性空间上的距离的讨论及推广。
12。
凸集与不动点定理.ﻫ13。
Hilbert空间的同构.ﻫ14。
最佳逼近问题。
ﻫ15。
线性函数的概念及推广.ﻫ16.一类椭圆型方程的解.18.线性赋范空间上的模等价。
17。
泛函分析中的不变子空间。
ﻫ19.范数的概念及性质.20。
正交与正交基的概念。
22。
隐函数存在定理的再证明。
ﻫ23.线性空间的等距同构。
21。
压缩映像原理及其应用.ﻫ24。
列紧集的概念及相关推广。
25。
Lebesgue控制收敛定理及应用。
26。
Lebesgue积分与Riemann积分的关系。
27。
重积分与累次积分的关系.28。
可积函数与连续函数的关系。
29。
有界变差函数的概念及其相关概念。
ﻫ30。
绝对连续函数的性质。
31.Lebesgue测度的相关概念。
33。
可测函数的定义及其性质。
ﻫ34.分部积分公式的32。
可测函数与连续函数的关系。
ﻫ推广。
35。
Fatou引理的重要作用。
36.不定积分的微分的计算。
ﻫ37。
绝对连续函数与微积分基本定理的关系。
ﻫ38。
Schwartz 不等式及推广。
39。
阶梯函数的概念及其作用.40。
Fourier级数及推广。
ﻫ41.完全正交系的概念及其作用。
ﻫ42。
Banach空间与Hilbe rt空间的关系。
44。
数学分析中的构造法证题术,43。
函数的各种收敛性及它们之间的关系。
ﻫ45。
用微积分理论证明不等式的方法46.数学分析中的化归法47。
微积分与辩证法49。
在上有界闭域的D中连续函数的性质48. 积分学中一类公式的证明ﻫ51。
常微分方程初等积分法解法研究常微分方程及积分因子初等积分法解常微分方程的关键在于求解不定积分。
不定积分是解微分方程的主要手段,通过找到合适的积分因子,可以将一个一阶微分方程转化为一个可积的方程。
在本文中,将对常微分方程及积分因子进行研究。
dy/dx = f(x, y)其中,f(x,y)是已知函数。
解这个方程的方法之一就是通过积分来找到y。
我们需要将这个方程转化成一个可积的形式。
考虑一个形式为 dy/dx + P(x)y = Q(x) 的一阶常微分方程。
要将这个方程转化为可积的形式,需要找到一个因子M(x),使得通过乘以M(x)可以使得原方程的左侧变为一个可积的形式。
这个因子M(x)被称为积分因子。
要找到积分因子,通常通过求解方程 M(x) = 1/M dM/dx = P(x) 来确定。
最常见的积分因子是指数函数,即M(x) = e^(∫P(x)dx)。
通过乘以这个积分因子,原方程可以变为积分形式:d/dx (M(x)y) = M(x)Q(x)通过对上式两边进行不定积分,可以求解出y。
举个例子来说明。
考虑一阶常微分方程 dy/dx + xy = x^2、我们需要找到一个积分因子。
通过解方程 M(x) = 1/M dM/dx = x,可以得到M(x) = e^(1/2 x^2)。
d/dx (e^(1/2 x^2) y) = x e^(1/2 x^2)对上式两边不定积分,得到:e^(1/2 x^2) y = ∫x e^(1/2 x^2) dx通过不定积分求解上式,可以得到y。
通过求解积分因子,我们可以将一阶常微分方程转化为可积的形式。
这种方法适用于一阶线性常微分方程。
对于高阶常微分方程,可以通过转化为一组一阶微分方程来求解。
总结起来,常微分方程及积分因子的研究是通过寻找积分因子来将一阶常微分方程转化为可积的形式。
通过解不定积分,可以求解出未知函数。
初等积分法解常微分方程是一种常用的方法,对于一阶线性常微分方程特别适用。
2015届本科毕业论文(设计) 论文题目:常微分方程初等解法的研究学院:数学科学学院专业班级:数学与应用数学11-1班学生姓名:汤鹏指导老师:张新东副教授答辩日期:2015年5月5日新疆师范大学教务处目录引言 (1)1 常微分方程的定义及分类 (2)1.1 定义 (2)1.2 一阶线性微分方程 (2)1.3 一阶线性微分方程组 (2)2 一阶线性微分方程的解法 (4)2.1 分离变量法 (4)2.2 常数变易法 (5)2.3 全微分法 (6)2.4 参数法 (7)3 n阶常系数线性微分方程的解法 (9)3.1 单根的情形 (9)3.2 重根的情形 (10)4 常微分方程的应用 (11)4.1 人口动力学问题 (11)4.2 简谐运动 (11)4.3 电路理论 (12)4.4 MATLAB解常微分方程 (13)5 总结 (15)参考文献 (16)致谢 (17)常微分方程初等解法的研究摘要:本文主要对常微分方程的初等解法进行研究,使大家更深一步地了解常微分方程的分类、解法及其在其他领域的应用。
首先总结阐述常微分方程的定义和几种常见的类型,然后讲解了常微分方程的解法及方程组解的情况,最后讲述了常微分方程在以下四个方面的应用:动力学问题、简谐运动、电路理论及用MATLAB解常微分方程。
关键词:常微分方程;初等解法;方程组;动力学;MATLABResearch elementary solution of ordinary differentialequationsAbstract: This paper mainly elementary solution of ordinary differential equation is studied,make you a deeper understanding of classification,the ordinary differential equation solution and its application in other fields.Firstly summarizes the type describes the definition of ordinary differential equations and several common,then explain the ordinary differential equation solution and the solution of equations,and finally describes the application of ordinary differential equations in the following four aspects:dynamics,simple harmonic motion,boundary value problem and the solution of ordinary differential equation with MATLAB.Key words: Ordinary differential equations; The primary solution; Equations; Dynamics; MATLAB引言常微分方程是数学中的一个重要的方程之一。
常微分方程是人类在生活实践中得来的。
据史料记载它的的出现要比微积分还要早。
笛卡尔在光学问题上的研究由切线性质引出的镜面形状、伽利略研究自由落体运动等等[10]。
事实上,这些问题都要建立并求解微分方程。
本文首先给出了常微分方程的相关定义、分类及其解法,想让大家对常微分方程的相关知识进行整理和汇总,在此基础上应用实际应用例子,以体现常微分方程的重要作用。
1 常微分方程的定义及分类1.1 定义一般来说,微分方程就是联系自变量、未知函数以及未知函数的某些导数(或微分)之间的关系式[7]。
当未知函数中依赖于一个自变量时,相应的微分方程称为常微分方程[8]。
1.2一阶线性微分方程一阶线性微分方程的形式:)()(x q y x p dx dy+= (1.2.1) 其中()x p 和()x q 是区间b x a <<上的已知函数。
如果()0≡x q ,即y x p dx dy)(= (1.2.2) 则称其为一阶线性齐次方程。
如果0)(≠x q ,则称(1.2.1)式为一阶线性非齐次方程[8]。
1.3一阶线性微分方程组 一阶线性微分方程组:如果在一阶微分方程组中,函数i f ),,2,1)(,,,,(21n i y y y x n ΛΛ=关于n y y y ,,,21Λ 是线性的,一阶微分方程组可以写成:()()()()()()()()()()()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧++++=++++=++++=x f y x a y x a y x a dx dy x f y x a y x a y x a dx dy x f y x a y x a y x a dx dy n n nn n n n n n n n ΛM ΛΛ2211222221212112121111(1.3.1)则称(1.3.1)为一阶线性微分方程组。
为了方便记忆,可以把(1.3.1)写成向量的形式。
为此,记()()()()()()()()()⎝⎛⎪⎪⎪⎪⎪⎭⎫=x a x a a x a x a x a x a x a x a x A nn n n n n ΛM M M ΛΛ212222111211 及()()()()⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫=x f x f x f x F n M 21(1.3.1)可以写成向量形式()()x F Y x A dx dY+= (1.3.2) 如果(1.3.2)上()0≡x F ,方程组(1.3.2)可变成()Y x A dx dY= (1.3.3) 则称(1.3.3)为一阶线性齐次方程组。
2 一阶线性微分方程的解法2.1 分离变量法1)显示变量可分离方程的解法形如()()y x f dxdyϕ=(2.1.1) 的方程,称为显示变量可分离方程。
如果()0≠y ϕ,我们可将(2.1.1)写成()()dx x f y dy=ϕ这样变量就分离开了。
两边积分可得 ()()c dx x f y dy+=⎰⎰ϕ(2.1.2) 则称(2.1.2)是(2.1.1)的通解。
如果存在0y ,使()00=y ϕ,直接代入,可知0y y =也是(2.1.1)的解,可能不在通解中,必须予以补上[9]。
例1 求解方程323a yb a dx dy -=解 分离变量得dx a dy a y b 332=- 积分得()dx a a y b a y b b3323221=-- ()33232232a x a yb a y b b=-- 这曲线就是摆线[3]。
2)微分形式变量可分离方程的解法形如()()()()dy y N x M dx y N x M 2211=(2.1.3)是变量可分离方程的微分形式表达式。
这是,y x 和在方程中的地位是“平等”的。
1)当()001=y N ,则0y y =为方程(2.1.3)的解。
同理()002=x M ,则0x x =也是方程(2.1.3)的解。
2)当()()时021≠x M y N ,分离变量可得:()()()()dx x M x M dy y N y N 2112=对上式两端同时积分得(2.1.3)的通积分()()()()C dx x M x M dy y N y N +=⎰⎰2112例2 求解方程()()01122=-+-dy x y dx y x 解 首先,易见1,1±=±=x y 是方程的解,其次,当()()01122≠--y x 时,分离变量得01122=-+-y ydyx xdx 积分,得方程的通积分C y x ln 1ln 1ln 22=-+- (0≠C )或()()C y x=--1122(0≠C )2.2 常数变易法常数变易法主要针对是一阶线性非齐次方程。
即)()(x f y x p dxdy=+ 对于一阶齐次方程0)(=+y x p dx dy的通解是()⎰=-dxx p Ce y 使用常数变易法解决一阶非齐次方程的解,可将常数C 变易成函数()x C ,即令()()⎰=-dx x p e x C y (2.1.4) 为方程(1.2.1)的解,其中()x C 待定。
将(2.1.4)代入(1.2.1)中并积分可以得到()()()C e x f x C dxx p +⎰=⎰将其代入(2.1.4)得(1.2.1)的通积分公式为()()()()dx e x f e Ce y dx x p dx x p dx x p ⎰⎰⎰+⎰=-- 例3 求解方程2x x ydx dy +=解 先求对应齐次方程x y dx dy = 的通解是Cx y = 由常数变易法得,令()x x C y = 为原方程的解,并代入原方程有()()()2'x x C x C x x C +=+ 整理并积分得 ()C x x C +=221 代入原方程的通解为321x Cx y +=2.3全微分法如果微分形式的一阶方程()()0,,=+dy y x N dx y x M (2.3.1) 的左端恰好是一个二元函数()y x U ,的全微分,即()()()dy y x N dx y x M y x dU ,,,+= (2.3.2) 则称(2.3.1)是全微分方程[1]。
下面将以例题的形式介绍全微分法例4 求解方程 ()()046633222=+++dy y y x dx xy x解 因为xNxy y M ∂∂==∂∂12 所以原方程是全微分方程。
为了计算方便我们取0,000==y x 故方程的通积分为()⎰⎰=++yx C dy y dx xy x3022463即C y y x x =++422332.4参数法参数法主要针对的是一阶隐式微分方程()0,,'=y y x F(1)如果能解出'y ,就得到一个或者几个显示微分方程,能用初等积分很容易解出方程的通解。
(2)如果不能解出'y ,这就要用到参数法。
本文主要介绍一类可积分类型。
即,()0,'=y x F ()()0,'=y y F (1)首先讨论()0,'=y x F (2.4.1) 1)把方程(2.4.1)化成参数形式。
即,()()⎩⎨⎧==t y t x ψϕ't 为参数 (2.4.2) 2)对于(2.4.2)和沿着(2.4.1)的任何一条积分曲线恒满足基本关系式dx y dy '= 这样,把(2.4.2)代入上式并积分得()()C dt t t y +=⎰'ϕψ于是得到方程(2.4.1)的参数通解()()()⎪⎩⎪⎨⎧+==⎰C dt t t y t x 'ϕψϕ 同理,可以讨论()0,'=y y F (2.4.3)设其可以表示的参数形式()()⎩⎨⎧==t y t x ψϕ',由于dy y dx '1=,将其代入参数方程组中积分得 ()()C dt t t x +=⎰ψϕ' 从而(2.4.3)的参数形式通解为()()()⎪⎩⎪⎨⎧=+=⎰t y C dt t t x ϕψϕ' 例5 求解方程 ()'2'1y y x =+解 令()t y tan '=,有()t x sin =,原方程的参数形式为⎩⎨⎧==t y tx tan sin 由基本关系dx y dy '=有tdt dy sin =,积分可的C t y +-=cos 。