常微分方程课件-中山版

常微分方程OrdinaryDifferentialEquation2014-2015学年第一学期刘汉泽修改[e-mail]•课程安排：计划上课18周（除去节假日、劳动周），从9月1日开始，单周4节；双周2节，上机。教材及参考资料•教材：常微分方程，(第三版）（2007年教育部精品教材），王高雄等（中山大学),高教出版社•参考书目：[1]常微分方程,东北师大数学系编，高教出版社[2]常微分方程讲义，王柔怀、伍卓群编，高教出版社[3]常微分方程教程，丁同仁等编，高教出版社[4]微分方程定性理论，张芷芬等编，科学出版社。教学安排•第1周——第18周，共54学时（含国庆等假期，实际课时更少）•考试安排：按学校、学院统一安排，•总成绩=平时（30%）+期末（70%），有小论文可以加分，一般每周四课代表收交作业，并统计作业情况。第一章绪论常微分方程是现代数学的一个重要分支，是人们解决各种实际问题的有效工具，它在几何、力学、物理、电子技术、航空航天、生命科学、经济领域等都有广泛的应用。随着计算技术和计算机的快速发展，常微分方程已经渗透到自然科学、社会科学、工程技术等学科的任何一个领域，正发挥着越来越大的作用。动力系统•Dynamicalsystemdescribestheevolutionofastateovertime••Curator:Dr.EugeneM.Izhikevich,Editor-in-ChiefofScholarpedia,thefreepeerreviewedencyclopedia第一章绪论主要内容•线性方程、二次方程、高次方程、指数方程、对数方程、三角方程和方程组•这些方程都是要把研究问题中的已知数和未知数之间的关系找出来，列出包含一个未知数或几个未知数的一个或者多个方程式，统称代数方程。•在实际工作中，常常出现一些特点和以上方程完全不同的问题•比如：某个物体在重力作用下自由下落，要寻求下落距离随时间变化的规律•火箭在发动机推动下在空间飞行，要寻求它飞行的轨道等•研究这些问题所建立的数学方程不仅与未知函数有关，而且与未知函数的导数有关，这就是我们要研究的微分方程•基本思想:把研究的问题中已知函数和未知函数之间的关系找出来，从列出的包含未知函数及其导数的一个或几个方程中去求得未知函数的表达式，即求解微分方程。微分方程的历史•微分方程差不多是和微积分同时产生•牛顿在建立微积分的同时，对简单的微分方程用级数来求解•瑞士数学家雅各布·贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论•法国数学家Poincare及前苏联数学家Lyapunov等对现代微分方程理论的建立做出了巨大的贡献与其他学科的关系•常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的•数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响•当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具1.1常微分方程模型•RLC电路•数学摆•人口模型•传染病模型•两生物种群生态模型•Lorenz方程RL电路基尔霍夫(Kirchhoff)第二定律在闭合回路中，所有支路上的电压的代数和等于零RLC电路数学摆人口模型•马尔萨斯(Malthus)假设：在人口自然增长的过程中，净相对增加率（单位时间内人口的净增长数与人口总数之比）是常数，记为r人口模型的改进•Verhulst：引入常数Nm（环境最大容纳量），假设：净相对增长率为))(1(mNtNrlogistic模型传染病模型•假设传染病传播期间其地区总人数不变，为常数n，开始时染病人数为x0，在时刻t的健康人数为y(t)，染病人数为x(t)•假设单位时间内一个病人能传染的人数与当时的健康人数成正比，比例系数为kSI模型易感染者：Susceptible已感染者：InfectiveSIS模型•对无免疫性的传染病，假设病人治愈后会再次被感染，设单位时间治愈率为muSIR模型（R：移出者（Removed））•对有很强免疫性的传染病，假设病人治愈后不会在被感染，设在时刻t的愈后免疫人数为r(t)，称为移出者，而治愈率l为常数两生物种群生态模型•意大利数学家沃特拉(Volterra)建立了一个关于捕食鱼与被食鱼生长情形的数学模型•假设在时刻t，被食鱼的总数为x(t)，而捕食鱼的总数为y(t)•假设单位时间内捕食鱼与被捕食鱼相遇的次数为bxy•捕食鱼的自然减少率同它们的存在数目y成正比Volterra被捕食-捕食模型两种群竞争模型Lorenz方程Lorenz吸引子，蝴蝶效应对初值的敏感性分形（fractal）吸引子总结•微分方程反映量与量之间的关系，与时间有关，是一个动态系统•从已知的自然规律出发，考虑主要因素，构造出由自变量、未知函数及其导数的关系史，即微分方程，从而建立数学模型•数学模型的建立有多种方式•研究微分方程的解和解结构的性质，检查是否与实际相吻合，不断改进模型•由微分方程发现或预测新的规律和性质1.2基本概念•1.2.1常微分方程基本概念定义（微分方程）联系自变量、未知函数及未知函数导数（或微分）的关系式称为微分方程，DE;2)1(xdxdy;0(2)ydxxdy;0)3(322xdtdxtxdtxd;sin35)4(2244txdtxddtxd;)5(zyzxz.0)6(2222uzyxyuxu例1：下列关系式都是微分方程微分方程如果在一个微分方程中，自变量的个数只有一个，则这样的微分方程称为常微分方程，ODE;2)1(xdxdy;0(2)ydxxdy;0)3(322xdtdxtxdtxd;sin35)4(2244txdtxddtxd都是常微分方程常微分方程如如果在一个微分方程中,自变量的个数为两个或两个以上,称为偏微分方程，PDE;)5(zyzxz.0)6(2222uzyxyuxu注:本课程主要研究常微分方程，同时把常微分方程简称为微分方程或方程偏微分方程如都是偏微分方程定义微分方程中出现的未知函数的最高阶导数或微分的阶数称为微分方程的阶数.2)1(xdxdy是一阶微分方程0(2)ydxxdy是二阶微分方程0)3(322xdtdxtxdtxd是四阶微分方程sin35)4(2244txdtxddtxd微分方程的阶如:)1(0),,dxdyy,F(x,nndxydn阶微分方程的一般形式为.,,,,,dxdyy,x,0),,dxdyy,F(x,是自变量是未知函数而且一定含有的已知函数是这里xydxyddxyddxydnnnnnn2)1(xdxdy是线性微分方程0(2)ydxxdysin35)4(2244txdtxddtxd线性和非线性0),,dxdyy,F(x,nndxyd如.,,,dxdyy阶线性方程则称其为的一次有理式及的左端为ndxydnn如果方程是非线性微分方程如0)3(322xdtdxtxdtxdn阶线性微分方程的一般形式111()()()(2)nnnnndydyaxaxyfxdxdx.)(),(),(1的已知函数是这里xxfxaxan不是线性方程的方程称为非线性方程微分方程的解定义:,),(满足条件如果函数Ixxy;)()1(阶的连续导数上有直到在nIxy,0))(),(),(,(:)2('xxxxFIxn有对.0),,dxdyy,F(x,(x)y上的一个解在为方程则称Idxydnn)(xy称为方程的显式解例.),(0ycosxysinx,y上的一个解在都是微分方程验证y证明:由于对sinx,yxsinycosx,y'(,),x故对有yyxsin0xsin.),(0ysinxy上的一个解在是微分方程故y.),(0yxcosy上的一个解在是微分方程同理y显式解与隐式解是方程的一个则称的解为方程所确定的隐函数如果关系式0),(,0),,dxdyy,F(x,Ix(x),y0),(yxdxydyxnn隐式解注:显式解与隐式解统称为微分方程的解，也叫微分方程的积分例如yxdxdy对一阶微分方程有显式解2211.yxyx和和隐式解:.122yx通解与特解定义如果微分方程的解中含有任意常数,且所含的相互独立的任意常数的个数与微分方程的阶数相同,则称这样的解为该方程的通解例如:为任常数2121,ccosx,sinxyccc.0y的通解是微分方程yn阶微分方程通解的一般形式为),,,(1nccxy.,,1为相互独立的任常数其中ncc注:使得行列式的某一邻域存在是指个独立常数含有称函数,),,,(,),,,(11nnccxnccxy0),,,(),,,()1(2)1(1)1('2'1'2121)1('nnnnnnnncccccccccccc.)(kkkdxd表示其中例.62y2y3cy''2321的通解是微分方程验证yyececexxxxxxecece2321'2cy证明:由于,4cy2321''xxxececexxxecece2321'''8cy故yy2y2y'')2(c2321xxxecece)8(c2321xxxecece)4(c22321xxxecece)32(c2321xxxecece6xe)c2cc2c(1111xecccc)22(-2222xecccc23333)228(86.62y2y3cy''2321的通解是微分方程故yyececexxx又3''3''1''3'2'1'321ccccccccc2222264xxxxxxxxxxeeeeeeeeee0.62y2y3cy''2321的解微分方程是故yyececexxx类似可定义方程的隐式通解如果微分方程的隐式解中含有任意常数,且所含的相互独立的任意常数的个数与微分方程的阶数相同,则称这样的解为该方程的隐式通解以后不区分显式通解和隐式通解,统称为方程的通解隐式通解也称为“通积分”在通解中给任意常数以确定的值而得到的解称为方程的特解例如.0ycosxysinx,y的特解都是方程y中分别取可在通解cosxsinxy21cc:,0,1c21得到c:,1,0c21得到csinx,ycosx.y定义定解条件为了从通解中得到合乎要求的特解,必须根据实际问题给微分方程附加一定的条件,称为定解条件求满足定解条件的求解问题称为定解问题常见的定解条件是初始条件,n阶微分方程的初始条件是指如下的n个条件:)1(01)1()1(000,,,,xxnnnydxydydxdyyy时当.1,,,,)1(0)1(000个常数是给定的这里nyyyxn当定解条件是初始条件时,相应的定解问题称为初值问题注1:n阶微分方程的初始条件有时也可写为)1(010)1()1(0000)(,,)(,)(nnnydxxydydxxdyyxy通常记为问题的解的初值问题也称满足条件阶微分方程求,)(,,)(,)(,0),,,,(:)1(010)1()1(0000CauchyydxxydydxxdyyxydxyddxdyyxFnnnnnn注2:0),,,,(nndxyddxdyyxF)