常微分方程王高雄1

常微分方程OrdinaryDifferentialEquation,教材及参考资料,教材：常微分方程(第三版），王高雄等,高教出版社参考书目：1常微分方程,东北师大数学系编，高教出版社2常微分方程讲义，王柔怀、伍卓群编，高教出版社3常微分方程及其应用，周义仓等编，科学出版社4微分方程定性理论，张芷芬等编，科学出版社,教学安排,第1周第18周，共72学时第5周周一(4.6)清明节，第8周周五(5.1)劳动节，第16周周一(6.22)端午节放假，实际授课时66学时）授课内容：第1章第5章考试安排：在结课后一周考试总成绩=平时（20%）+期末（80%），有小论文可以加分，每周五课后按学号尾数单双号交替交作业答疑地点:25教15楼尚月强办公室,第一章绪论,常微分方程是现代数学的一个重要分支，是人们解决各种实际问题的有效工具，它在几何、力学、物理、电子技术、航空航天、生命科学、经济领域等都有广泛的应用。随着计算技术和计算机的快速发展，常微分方程已经渗透到自然科学、社会科学、工程技术等学科的任何一个领域，正发挥着越来越大的作用。,动力系统,Dynamicalsystemdescribestheevolutionofastateovertime/article/History_of_dynamical_systems,PhilipHolmes(2007),Scholarpedia,2(5):1843Curator:Dr.EugeneM.Izhikevich,Editor-in-ChiefofScholarpedia,thefreepeerreviewedencyclopedia,第一章绪论,线性方程、二次方程、高次方程、指数方程、对数方程、三角方程和方程组。这些方程都是要把研究问题中的已知数和未知数之间的关系找出来，列出包含一个未知数或几个未知数的一个或者多个方程式。,在实际工作中，常常出现一些特点和以上方程完全不同的问题比如：某个物体在重力作用下自由下落，要寻求下落距离随时间变化的规律火箭在发动机推动下在空间飞行，要寻求它飞行的轨道等研究这些问题所建立的数学方程不仅与未知函数有关，而且与未知函数的导数有关，这就是我们要研究的微分方程,基本思想:把研究的问题中已知函数和未知函数之间的关系找出来，从列出的包含未知函数及其导数的一个或几个方程中去求得未知函数的表达式，即求解微分方程,微分方程差不多是和微积分同时先后产生的牛顿在建立微积分的同时，对简单的微分方程用级数来求解瑞士数学家雅各布贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论法国数学家Poincare及前苏联数学家Lyapunov等对现代微分方程理论的建立做出了巨大的贡献,常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具,1.1常微分方程模型,RLC电路数学摆人口模型传染病模型两生物种群生态模型Lorenz方程,RL电路,基尔霍夫(Kirchhoff)第二定律,在闭合回路中，所有支路上的电压的代数和等于零,RLC电路,数学摆,数学摆是系于一根长度为的线上而质量为的质点M.在重力作用下,它在垂直于地面的平面上沿圆周运动.如图所示.试确定摆的运动方程.,解:Newton第二定律:,取反时针运动方向为计量摆与铅垂线所成的角的正方向.则由Newton第二定律,得到摆的运动方程为,附注1:如果研究摆的微小振动,即当比较小时,可以取的近似值代入上式,这样就得到微小振动时摆的运动方程:,附注2:假设摆是在一个有粘性的介质中作摆动,如果阻力系数为则摆的运动方程为:,附注3:假设摆还沿着摆的运动方向受到一个外力F(t)的作用,则摆的运动方程为:,人口模型,马尔萨斯(Malthus)假设：在人口自然增长的过程中，净相对增加率（单位时间内人口的净增长数与人口总数之比）是常数，记为r,人口模型的改进,Verhulst：引入常数Nm（环境最大容纳量），假设：净相对增长率为,logistic模型,传染病模型,假设传染病传播期间其地区总人数不变，为常数n，开始时染病人数为x0，在时刻t的健康人数为y(t)，染病人数为x(t)假设单位时间内一个病人能传染的人数与当时的健康人数成正比，比例系数为k,SI模型易感染者：Susceptibley(t)已感染者：Infectivex(t),SIS模型,对无免疫性的传染病，假设病人治愈后会再次被感染，设单位时间治愈率为mu,SIR模型（R：移出者（Removed）,对有很强免疫性的传染病，假设病人治愈后不会在被感染，设在时刻t的愈后免疫人数为r(t)，称为移出者，而治愈率l为常数,两生物种群生态模型,意大利数学家沃特拉(Volterra)建立了一个关于捕食鱼与被食鱼生长情形的数学模型假设在时刻t，被食鱼的总数为x(t)，而捕食鱼的总数为y(t)假设单位时间内捕食鱼与被捕食鱼相遇的次数为bxy捕食鱼的自然减少率同它们的存在数目y成正比捕食鱼的自然增长率同它们它们的存在数目y及被被捕食鱼x成正比,Volterra被捕食-捕食模型,两种群竞争模型,相互竞争同一资源：,Lorenz方程,Lorenz吸引子,系统总体稳定性和局部不稳定性:在吸引子外的一切运动都趋向（吸引）到吸引子，一切到达吸引子内的运动都互相排斥,蝴蝶效应,一只南美洲亚马孙河流域热带雨林中的蝴蝶，偶尔扇动几下翅膀，可以在两周以后引起美国德克萨斯州的一场龙卷风,对初值的敏感性,电影蝴蝶效应,当一个人小时候受到微小的心理刺激，长大后这个刺激会被放大,分形（fractal）,一个粗糙或零碎的几何形状，可以分成数个部分，且每一部分都（至少近似地）是整体缩小后的形状。,总结,微分方程反映量与量之间的关系，与时间有关，是一个动态系统从已知的自然规律出发，考虑主要因素，构造出由自变量、未知函数及其导数的关系史，即微分方程，从而建立数学模型数学模型的建立有多种方式研究微分方程的解和解结构的性质，检查是否与实际相吻合，不断改进模型由微分方程发现或预测新的规律和性质,1.2基本概念与常微分方程的发展史,1.2.1常微分方程基本概念,定义（微分方程）联系自变量、未知函数及未知函数导数（或微分）的关系式称为微分方程,例1：下列关系式都是微分方程,微分方程,如果在一个微分方程中，自变量的个数只有一个，则这样的微分方程称为常微分方程,都是常微分方程,常微分方程,如,如果在一个微分方程中,自变量的个数为两个或两个以上,称为偏微分方程,注:本课程主要研究常微分方程，同时把常微分方程简称为微分方程或方程,偏微分方程,如,都是偏微分方程,定义微分方程中出现的未知函数的最高阶导数或微分的阶数称为微分方程的阶数.,是一阶微分方程,是二阶微分方程,是四阶微分方程,微分方程的阶,如:,n阶微分方程的一般形式为,是线性微分方程,线性和非线性,如,如果方程,是非线性微分方程,如,n阶线性微分方程的一般形式,不是线性方程的方程称为非线性方程,微分方程的解,定义,称为方程的显示解,例,证明:,显式解与隐式解,隐式解,注:显式解与隐式解统称为微分方程的解,例如,有显式解,和隐式解:,通解与特解,定义如果微分方程的解中含有任意常数,且所含的相互独立的任意常数的个数与微分方程的阶数相同,则称这样的解为该方程的通解,例如:,n阶微分方程通解的一般形式为,注:,例,证明:,由于,故,又,类似可定义方程的隐式通解,如果微分方程的隐式解中含有任意常数,且所含的相互独立的任意常数的个数与微分方程的阶数相同,则称这样的解为该方程的隐式通解,以后不区分显式通解和隐式通解,统称为方程的通解,隐式通解也称为“通积分”,在通解中给任意常数以确定的值而得到的解称为方程的特解,例如,定义,定解条件,为了从通解中得到合乎要求的特解,必须根据实际问题给微分方程附加一定的条件,称为定解条件,求满足定解条件的求解问题称为定解问题,常见的定解条件是初始条件,n阶微分方程的初始条件是指如下的n个条件:,当定解条件是初始条件时,相应的定解问题称为初值问题,注1:n阶微分方程的初始条件有时也可写为,注2:,例（P19）,解,由于,且,解以上方程组得,积分曲线和方向场,积分曲线,一阶微分方程,称为微分方程的积分曲线,方向场,在方向场中,方向相同的点的几何轨迹称为等斜线.,所规定的方向场,常微分方程求解的几何意义：在方向场中寻求一条曲线，使这条曲线上每一点切线的方向等于方向场中该点的方向。,例画出方程,的方向场。,等倾线方程,即,也就是说，方向场中每点的方向与该点等倾线垂直。,x,y,o,微分方程组,驻定与非驻定，动力系统,驻定（自治）,非驻定（非自治）,相空间、奇点和轨线,不含自变量、仅由未知函数组成的空间称为相空间积分曲线在相空间中的投影称为轨线称为平衡解（驻定解、常数解），奇点、平衡点,垂直等倾线、水平等倾线,1.2.2常微分方程的发展史,300多年前，Newton与Leibniz奠定微积分基本思想的同时,就正式提出了微分方程的概念.17世纪末到18世纪，常微分方程研究的中心问题是如何求出通解的表达式