常微分方程在数学建模中的应用 这里介绍几个典型的用微分方程建立数学模型的例子. 一、人口预测模型 由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去，则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 例1（ 马尔萨斯 （Malthus ） 模型） 英国人口统计学家马尔萨斯（1766—1834）在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是：在人口自然增长过程中,净相对增长（出生率与死亡率之差）是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型. 解 设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理（因人口总数很大,可近似地这样处理,此乃离散变量连续化处理）,据马尔萨斯的假设,在t 到t t ?+时间段内,人口的增长量为 t t rN t N t t N ?=-?+)()()(, 并设0t t =时刻的人口为0N ,于是 ?????==． ， 00)(d d N t N rN t N 这就是马尔萨斯人口模型,用分离变量法易求出其解为 )(00e )(t t r N t N -=, 此式表明人口以指数规律随时间无限增长. 模型检验：据估计1961年地球上的人口总数为9 1006.3?,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3?=N ,02.0=r ,于是 ) 1961(02.09 e 1006.3)(-?=t t N . 这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人 口大约每35年翻一番,而上式断定34.6年增加一倍（请读者证明这一点）． 但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地（事实上,地球表面还有80%被水覆盖）,我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改. 例2（逻辑Logistic 模型） 马尔萨斯模型为什么不能预测未来的人口呢？这主要是地
常微分方程的实际应用 于萍 摘要：常微分方程在当代数学中是极为重要的一个分支，它的实用价值很高，应用也很广泛，本文主要介绍常微分方程在几何、机械运动、电磁振荡方面的应用，并举例说明，体会常微分方程对解决实际问题的作用，在解决实际问题过程中通常是建立起实际问题的数学模型，也就是建立反映这个实际问题的微分方程，求解这个微分方程，用所得的数学结果解释实际问题，从而预测到某些物理过程的特定性质，以便达到能动地改造世界，解决实际问题的目的。 关键字：常微分方程，几何，机械运动，电磁振荡，应用
Abstract: Nomal differential equation is an important part of math at it has a high practical value. This thesis shows the use in geometry, mechaics and electrothermal and makes some examples. Also, it summarizes the normal move of dealing with practical problems by the normal differential equation. Normal, we set up the maths matic model of the problem, solute the normal differentical equation make the use of the result to explain practical problems and make a forecast of some special character of physical process. Key: Normal differetial equation geometry mechanics electrothermal use
第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )() (=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、 xy dx dy = 解：当0≠y 时，有 xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(112 12 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(12 12 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有 dy y N y Q dx x P x M ) () ()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(2 2 =-+-dy x y dx y x 解：当0)1)(1(2 2 ≠--y x 时，有 dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(2 2 =--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(2 2 为常数C C y x =--。 ⑵可化为变量可分离方程的方程： ①、形如 )(x y g dx dy = 解法：令x y u =，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到
微分方程在高中物理中的应用 高中阶段，我们经常会遇到一些需要定性分析的物理问题，其实如果我们应用高等数学 的知识，可以把其中一些问题进行定量的分析。 例如，质量为m 的物体从高度H 自由下落，所受阻力f 与速度v 成正比，g 为重力加速 度这是我们平时常见的一类问题。但我们只知道速度V 最终会趋近于某一数值v0。下面我 进行一下定量分析。 根据题目所给信息，可列出动力学方程 mg-kv=ma ① a=dv/dt ② 结合①式可得mg-kv=mdv/dt 这里移项可得dt=mdv/(mg-kv)③ 两边同时积分便可的到 V=mg(ce*(-kt/m)+1)/k 又∵自由下落，可得t=0时v=.0 ∴v=mg(1-e*(-kt/m))/k ④ 由④式知，当t 趋近于正无穷时，e*(-kt/m)=0， 此时v=mg/k ⑤ 若按照正常思路，当物体受力平衡时，mg=kv,此时也能得到⑤式的结论。 而在高考中，更为常见的是在电磁场中的同类问题，我们不妨看一下下面这一道例题 （2012·山东理综）如图所示，相距为L 的两条足够长的光滑平行金属导轨与水平面的夹 角为θ，上端接有定值电阻，匀强磁场垂直于导轨平面，磁感应强度为B 。将质量为m 的导 体棒由静止释放，当速度达到v 时开始匀速运动，此时对导体棒施加一平行于导轨向下的 拉力，并保持拉力的功率为P ，导体棒最终以2v 的速度匀速运动。导体棒始终与导轨垂直 且接触良好，不计导轨和导体棒的电阻，重力加速度为g ，下列选项正 确的是 A ．P =2mg sin θ B ．P =3mg sin θ C ．当导体棒速度达到v /2时加速度为12 g sin θ D ．在速度达到2v 以后匀速运动的过程中，R 上产生的焦耳热等于拉力 所做的功 我们根据题目也可以列出动力学方程 Mgsin θ-B*2L*2V/R=ma ① a=dv/dt ② 同样可以解得v=(mgR sin θ/B*2L*2)(1-e*(-B*2L*2t/mR))③ 从③式可以看出当t 趋近于正无穷时，v=mgR sin θ/B*2L*2即B*2L*2v/R=mg sin θ转化而来。 所以题目中所说当速度到达V 时开始匀速运动存在明显错误。应改为近似于做匀速直线运 动。
微分方程应用 1 引言 常微分方程的形成与发展和很多学科有着密切的联系，例如力学、天文学、物理学等.数学的其他分支的快速发展，产生出很多新兴学科，这些新兴学科的产生都对常微分方程的发展有着深刻的影响，而且当前计算机的快速发展更是为常微分方程的应用及理论研究提供了非常有力的工具. 数学解决实际问题就必须建立模型，而数学建模就是把数学语言描述实际现象的过程.利用数学去解决各类实际问题时，建立数学模型是十分重要的一步，但是也是最困难的一步.建立数学模型的过程，是把错综复杂的实际问题简化、抽象为合理的数学结构的过程.要通过大量调查、收集相关数据资料，观察和研究实际对象的固有特征和内在规律，抓住问题的主要矛盾，建立起反映实际问题的数量关系，然后利用数学的理论和方法去分析和解决问题. 因此本文先简要介绍了如何建立微分方程模型，并通过具体的实例来简单地介绍了微分方程在数学建模中的应用. 2 数学模型简介 通常我们把现实问题的一个模拟称为模型．如交通图、地质图、航空模型和建筑模型等．利用字母、数学及其它数学符号建立起来的等式或不等式以及图表、图象、框图等来模拟现实的模型称为数学模型．数学模型在实际生活中经常碰到，如求不规则图形的面积，可建立定积分的数学模型，求变化率的问题可建立导数模型，统计学中抽样调查，买彩票中奖的概率问题等等．学会建立数学模型对解决实际生活问题会有很大的帮助． 建立数学模型是沟通摆在面前的实际问题与数学工具之间联系的一座必不可少的桥梁．随着科学技术的进步，特别是电子计算机技术的迅速发展，数学已经渗透到从自然科学技术到工农业生产建设，从经济生活到社会生活的各个领域．一般地说，当实际问题需要我们对所研究的现实对象提供分析、预报、决策、控制等方面的定量结果时，往往都离不开数学的应用，而建立数学模型则是这个过程的关键环节． 3 常微分方程模型 3．1 常微分方程的简介
常微分方程解题方法总结 来源：文都教育 复习过半， 课本上的知识点相信大部分考生已经学习过一遍 . 接下来， 如何将零散的知 识点有机地结合起来， 而不容易遗忘是大多数考生面临的问题 . 为了加强记忆， 使知识自成 体系，建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴， 他强调读 书要“由薄到厚、由厚到薄”，对同学们的复习尤为重要 . 以常微分方程为例， 本部分内容涉及可分离变量、 一阶齐次、 一阶非齐次、 全微分方程、 高阶线性微分方程等内容， 在看完这部分内容会发现要掌握的解题方法太多， 遇到具体的题 目不知该如何下手， 这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形 式将常微分方程中的解题方法加以总结，一目了然，便于记忆和查询 . 常微分方程 通解公式或解法 ( 名称、形式 ) 当 g( y) 0 时，得到 dy f (x)dx ， g( y) 可分离变量的方程 dy f ( x) g( y) 两边积分即可得到结果； dx 当 g( 0 ) 0 时，则 y( x) 0 也是方程的 解 . 解法：令 u y xdu udx ，代入 ，则 dy 齐次微分方程 dy g( y ) x dx x u g (u) 化为可分离变量方程 得到 x du dx 一 阶 线 性 微 分 方 程 P ( x)dx P ( x) dx dy Q(x) y ( e Q( x)dx C )e P( x) y dx
伯努利方程 解法：令 u y1 n，有 du (1 n) y n dy ， dy P( x) y Q( x) y n（n≠0,1）代入得到du (1 n) P(x)u (1 n)Q(x) dx dx 求解特征方程：2 pq 三种情况： 二阶常系数齐次线性微分方程 y p x y q x y0 二阶常系数非齐次线性微分方程 y p x y q x y f ( x) (1)两个不等实根：1, 2 通解： y c1 e 1x c2 e 2x (2) 两个相等实根：1 2 通解： y c1 c2 x e x (3) 一对共轭复根：i , 通解： y e x c1 cos x c2 sin x 通解为 y p x y q x y 0 的通解与 y p x y q x y f ( x) 的特解之和. 常见的 f (x) 有两种情况： x （ 1）f ( x)e P m ( x) 若不是特征方程的根，令特解 y Q m ( x)e x；若是特征方程的单根，令特 解 y xQ m ( x)e x；若是特征方程的重根， 令特解 y*x2Q m (x)e x； （2）f (x) e x[ P m ( x) cos x p n ( x)sin x]
淮北师范大学 2013届学士学位论文 常微分方程数值解法的误差分析 学院、专业数学科学学院数学与应用数学 研究方向计算数学 学生姓名李娜 学号 20091101070 指导教师姓名陈昊 指导教师职称讲师 年月日
常微分方程数值解法的误差分析 李娜 （淮北师范大学数学科学学院，淮北，235000） 摘要 自然界与工程技术中的很多现象，往往归结为常微分方程定解问题。许多偏微分方程问题也可以化为常微分方程问题来近似求解。因此，研究常微分方程的数值解法是有实际应用意义的。数值解法是一种离散化的数学方法，可以求出函数的精确解在自变量一系列离散点处的近似值。随着计算机计算能力的增强以及数值计算方法的发展，常微分方程的数值求解方法越来越多，比较成熟的有Euler 法、后退Euler法、梯形方法、Runge—Kutta方法、投影法和多步法，等等．本文将对这些解的误差进行分析,以求能够得到求解常微分数值解的精度更好的方法。 关键词：常微分方程, 数值解法, 单步法, 线性多步法, 局部截断误差
Error Analysis of Numerical Method for Solving the Ordinary Differential Equation Li Na (School of Mathematical Science, Huaibei Normal University, Huaibei, 235000) Abstract In nature and engineering have many phenomena , definite solution of the problem often boils down to ordinary differential equations. So study the numerical solution of ordinary differential equations is practical significance. The numerical method is a discrete mathematical methods, and exact solution of the function can be obtained in the approximation of a series of discrete points of the argument.With the enhanced computing power and the development of numerical methods,ordinary differential equations have more and more numerical solution,there are some mature methods. Such as Euler method, backward Euler method, trapezoidal method, Runge-Kutta method, projection method and multi-step method and so on.Therefore, numerical solution of differential equation is of great practical significance. Through this paper, error of these solutions will be analyzed in order to get a the accuracy better way to solve the numerical solution of ordinary differential. Keywords:Ordinary differential equations, numerical solution methods, s ingle ste p methods, l inear multi-step methods, local truncation error
常微分方程及其应用
第5章常微分方程及其应用 习题5.2 1．求下列各微分方程的通解： （1）?Skip Record If...?；（2）?Skip Record If...?； （3）?Skip Record If...?；（4）?Skip Record If...?； （5）?Skip Record If...?；（6）?Skip Record If...?． 2．求下列各微分方程满足所给初始条件的特解： （1）?Skip Record If...?，?Skip Record If...?；（2）?Skip Record If...?，?Skip Record If...?； （3）?Skip Record If...?，?Skip Record If...?；（4）?Skip Record If...?，?Skip Record If...?； （5）?Skip Record If...?，?Skip Record If...?；（6）?Skip Record If...?，?Skip Record If...?． 5.3 可降阶微分方程及二阶常系数线性微分方程 案例引入求微分方程?Skip Record If...?的通解． 解两边积分，得?Skip Record If...? 两边再积分，得?Skip Record If...? 所以，原方程的通解为?Skip Record If...?，其中?Skip Record If...?为任意常数． 5.3.1 可降阶微分方程 仅供学习与交流，如有侵权请联系网站删除谢谢20
1. 形如?Skip Record If...?的微分方程 特点：方程右端为已知函数?Skip Record If...?． 解法：对?Skip Record If...?连续积分?Skip Record If...?次，即可得含有 ?Skip Record If...?个任意常数的通解． 2. 形如?Skip Record If...?的微分方程 特点：方程右端不显含未知函数?Skip Record If...?． 解法：令?Skip Record If...?，则?Skip Record If...?．于是，原方程可化为?Skip Record If...?．这是关于?Skip Record If...?的一阶微分方程．设其通解为?Skip Record If...?，即?Skip Record If...?．两边积分，即可得原方程通解?Skip Record If...?，其中?Skip Record If...?为任意常数． 3. 形如?Skip Record If...?的微分方程 特点：方程右端不显含自变量?Skip Record If...?． 解法：令?Skip Record If...?，则?Skip Record If...?．于是，原方程可化为?Skip Record If...?．这是关于?Skip Record If...?的一阶微分方程．设其通解为?Skip Record If...?，即 ?Skip Record If...?．分离变量，得?Skip Record If...?．然后两边积分，即可得原方程通解 ?Skip Record If...?，其中?Skip Record If...?为任意常数．例5-7求微分方程?Skip Record If...?的通解． 解两边积分，得?Skip Record If...? 仅供学习与交流，如有侵权请联系网站删除谢谢20
目录 1 引言 (1) 2 二阶常系数常微分方程的几种解法 (1) 2.1 特征方程法 (1) 2.1.1 特征根是两个实根的情形 (2) 2.1.2 特征根有重根的情形 (2) 2.2 常数变异法 (4) 2.3 拉普拉斯变化法 (5) 3 常微分方程的简单应用 (6) 3.1 特征方程法 (7) 3.2 常数变异法 (9) 3.3 拉普拉斯变化法 (10) 4 总结及意义 (11) 参考文献 (12)
二阶常微分方程的解法及其应用 摘要：本文通过对特征方程法、常数变易法、拉普拉斯变换法这三种二阶常系数常微分方程解法进行介绍，特别是其中的特征方程法分为特征根是两个实根的情形和特征根有重根的情形这两种情况，分别使用特征值法、常数变异法以及拉普拉斯变换法来求动力学方程，现今对于二阶常微分方程解法的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题方面卓有成效。应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。 关键词：二阶常微分方程;特征分析法；常数变异法；拉普拉斯变换
METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper introduces the solution of the characteristic equation method, the method of variation of parameters, the Laplasse transform method the three kind of two order ordinary differential equations with constant coefficients, especially the characteristic equation method which is characteristic of the root is the two of two real roots and characteristics of root root, branch and don't use eigenvalue method, method of variation of constants and Laplasse transform method to obtain the dynamic equation, the current studies on solution of ordinary differential equations of order two has made many achievements, especially in the aspect of solving the problem of two order linear differential equation with constant coefficients very fruitful. Application of the theory of ordinary differential equations has made great achievements, however, the existing theory it is still far from meeting the need, needs further development, to make the discipline theory more perfect. Keywords:second ord er ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言 数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方程
i.常微分方程初值问题数值解法 常微分方程初值问题的真解可以看成是从给定初始点出发的一条连续曲线。差分法是常微分方程初值问题的主要数值解法，其目的是得到若干个离散点来逼近这条解曲线。有两个基本途径。一个是用离散点上的差商近似替代微商。另一个是先对微分方程积分得到积分方程，再利用离散点作数值积分。 i.1 常微分方程差分法 考虑常微分方程初值问题：求函数()u t 满足 (,), 0du f t u t T dt =<≤ （i.1a ） 0(0)u u = (i.1b) 其中(,)f t u 是定义在区域G : 0t T ≤≤, u <∞上的连续函数，0u 和T 是给定的常数。我们假设(,)f t u 对u 满足Lipschitz 条件，即存在常数L 使得 121212(,)(,), [0,]; ,(,)f t u f t u L u u t T u u -≤-?∈∈-∞∞ （i.2） 这一条件保证了（i.1）的解是适定的，即存在，唯一，而且连续依赖于初值0u 。 通常情况下，（i.1）的精确解不可能用简单的解析表达式给出，只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler 法。为此，首先将求解区域[0,]T 离散化为若干个离散点： 0110N N t t t t T -=<< <<= （i.3） 其中n t hn =，0h >称为步长。 在微积分课程中我们熟知，微商（即导数）是差商的极限。反过来，差商就是微商的近似。在0t t =处，在（i.1a ）中用向前差商 10()()u t u t h -代替微商du dt ，便得 10000()()(,())u t u t hf t u t ε=++ 如果忽略误差项0ε，再换个记号，用i u 代替()i u t 便得到 1000(,)u u hf t u -= 一般地，我们有 1Euler (,), 0,1, ,1n n n n u u hf t u n N +=+=-方法： （i.4） 从(i.1b) 给出的初始值0u 出发，由上式可以依次算出1,,N t t 上的差分解1,,N u u 。
《常微分方程》课程大纲 一、课程简介 课程名称：常微分方程学时/学分：3/54 先修课程：数学分析，高等代数,空间解析几何，或线性代数（行列式，矩阵与线性方程组，线性空间F n，欧氏空间R n，特征值与矩阵的对角化）, 高等数学（多元微积分，无穷级数）。 面向对象：本科二年级或以上学生 教学目标：围绕基本概念与基本理论、具体求解和实际应用三条主线开展教学活动，通过该课程的教学，希望学生正确理解常微分方程的基本概念，掌握基本理论和主要方法，具有一定的解题能力和处理相关应用问题的思维方式，如定性分析解的性态和定量近似求解等思想，并希望学生初步了解常微分方程的近代发展，为学习动力系统学科的近代内容和后续课程打下基础。 二、教学内容和要求 常微分方程的教学内容分为七部分，对不同的内容提出不同的教学要求。(数字表示供参考的相应的学时数，第一个数为课堂教学时数，第二个数为习题课时数) 第一章基本概念（2，0） （一）本章教学目的与要求： 要求学生正确掌握微分方程，通解，线性与非线性，积分曲线，线素场（方
向场），定解问题等基本概念。本章教学重点解释常微分方程解的几何意义。 （二）教学内容： 1．由实际问题：质点运动即距离与时间关系（牛顿第二运动定律），放射性元素衰变过程，人口总数发展趋势估计等，通过建立数学模型，导出微分方程。 2．基本概念（常微分方程，偏微分方程，阶，线性，非线性，解，定解问题，特解，通解等）。 3．一阶微分方程组的几何定义，线素场（方向场），积分曲线。 4．常微分方程所讨论的基本问题。 第二章初等积分法（4，2） （一）本章教学目的与要求： 要求学生熟练掌握分离变量法，常数变易法，初等变换法，积分因子法等初等解法。 本章教学重点对经典的几类方程介绍基本解法，勾通初等积分法与微积分学基本定理的关系。并通过习题课进行初步解题训练，提高解题技巧。 （二）教学内容： １. 恰当方程（积分因子法）; ２. 分离变量法 ３. 一阶线性微分方程（常数变易法） ４. 初等变换法（齐次方程，伯努利方程，黎卡提方程）
应 用 题（每题10分） 1、设()f x 在(,)-∞∞上有定义且不恒为零，又()f x '存在并对任意,x y 恒有()()()f x y f x f y +=，求()f x 。 2、设()()()F x f x g x =，其中函数(),()f x g x 在(,)-∞∞内满足以下条件 ()(),()(),(0)0,()()2x f x g x g x f x f f x g x e ''===+= （1）求()F x 所满足的一阶微分方程； （2）求出()F x 的表达式。 3、已知连续函数()f x 满足条件320 ()3x x t f x f dt e ??=+ ??? ?，求()f x 。 ; 4、已知函数()f x 在(0,)+∞内可导，()0,lim ()1x f x f x →+∞ >=，且满足 1 1 0()lim ()h x h f x hx e f x →? ?+ ?= ? ?? ? ，求()f x 。 5、设函数()f x 在(0,)+∞内连续，5 (1)2 f = ，且对所有,(0,)x t ∈+∞，满足条件 1 1 1 ()()()xt x t f u du t f u du x f u du =+? ??，求()f x 。 6、求连续函数()f x ，使它满足10 ()()sin f tx dt f x x x =+?? 。 7、已知可微函数()f t 满足 31() ()1()x f t dt f x t f t t =-+?，试求()f x 。 8、设有微分方程 '2()y y x ?-=， 其中21 ()01 x x x ??。试求在(,)-∞∞内的连续函 数()y y x =使之在(,1)-∞和()1,+∞内部满足所给方程，且满足条件(0)0y =。 9、设位于第一象限的曲线()y f x = 过点122?? ? ? ?? ，其上任一点(,)P x y 处的法线与y 轴的交点为Q ，且线段PQ 被x 轴平分。 （1）求曲线()y f x =的方程； （2）已知曲线sin y x =在[0,]π上的弧长为l ，试用l 表示曲线()y f x =的弧长s 。 ' 10、求微分方程(2)0xdy x y dx +-=的一个解()y y x =，使得由曲线()y y x =与直线 1,2x x ==以及x 轴所围成的平面图形绕x 轴旋转一周的旋转体体积最小。 11、设曲线L 位于xOy 平面的第一象限内，L 上任一点M 处的切线与y 轴总相交，交点记
各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u］=dx/x 两端积分,得∫du/［φ(u)-u］=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce－ ∫P(x)dx,再令y=u e－∫P(x)dx代入原方程解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e－∫P(x)dx［∫Q(x)e∫P(x)dx dx+C］ 即y=Ce－∫P(x)dx +e－ ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1 y(n-2)= ∫［∫f(x)dx+C1］dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程
令y ’=p 则y ”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C 1) 即dy/dx=φ(y,C 1),即dy/φ(y,C 1)=dx,所以∫dy/φ(y,C 1)=x+C 2 5.二阶常系数齐次线性微分方程解法 一般形式:y ”+py ’+qy=0,特征方程r 2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y ”+py ’+qy=f(x) 先求y ”+py ’+qy=0的通解y 0(x),再求y ”+py ’+qy=f(x)的一个特解y*(x) 则y(x)=y 0(x)+y*(x)即为微分方程y ”+py ’+qy=f(x)的通解 求y ”+py ’+qy=f(x)特解的方法: ① f(x)=P m (x)e λx 型 令y*=x k Q m (x)e λx ［k 按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m (x)的m+1个系数 ② f(x)=e λx ［P ｌ(x)cos ωx+P n (x)sin ωx ］型 令y*=x k e λx ［Q m (x)cos ωx+R m (x)sin ωx ］［m=max ﹛ｌ,n ﹜,k 按λ+i ω不是特征方程的根或是特征方程的单根依次取0或1］再代入原方程,分别确定Q m (x)和R m (x)的m+1个系数
第5章 常微分方程及其应用 习题5.2 1．求下列各微分方程的通解： （1）02 =+ydy dx x ； （2）0ln =-'y y y x ； （3）0)()(2 2 =-++dy y x y dx x xy ； （4）03=-'xy y ； （5）x e y y =-'2； （6）x x y y cos tan +='． 2．求下列各微分方程满足所给初始条件的特解： （1）y x e y -='2，0)0(=y ； （2） 011=+-+dy x y dx y x ，1)0(=y ； （3）x y y cos =-'，0)0(=y ； （4）x x y y sec tan =-'，0)0(=y ； （5）x x x y y sin = + '，1)(=πy ； （6）()0122 =+-+dx x xy dy x ，0)1(=y ． 5.3 可降阶微分方程及二阶常系数线性微分方程 案例引入 求微分方程x y 6=''的通解． 解 两边积分，得12 36C x xdx y +=='? 两边再积分，得 () 21312 3C x C x dx C x y ++=+= ? 所以，原方程的通解为213 C x C x y ++=，其中21C C 、为任意常数． 5.3.1 可降阶微分方程 1. 形如)() (x f y n =的微分方程 特点：方程右端为已知函数)(x f ． 解法：对)() (x f y n =连续积分n 次，即可得含有n 个任意常数的通解． 2. 形如),(y x f y '=''的微分方程 特点：方程右端不显含未知函数y ．
解法： 令)(x p y ='，则)(x p y '=''．于是，原方程可化为),(p x f p ='．这是关于 p p ',的一阶微分方程．设其通解为),()(1C x x p ?=，即),(1C x y ?='．两边积分，即可 得原方程通解21),(C dx C x y +=? ?，其中21C C 、为任意常数． 3. 形如),(y y f y '=''的微分方程 特点：方程右端不显含自变量x ． 解法：令)(y p y ='，则dy dp p dy dp y dx dy dy dp y ='=?= ''．于是，原方程可化为 ),(p y f p p ='．这是关于p p ',的一阶微分方程．设其通解为),()(1C y y p ψ=，即 ),(1C y dx dy ψ=．分离变量，得 dx C y dy =),(1 ψ．然后两边积分，即可得原方程通解 21) ,(C x C y dy +=? ψ，其中21C C 、为任意常数． 例5-7 求微分方程x x y cos sin -='''的通解． 解 两边积分，得12sin cos )cos (sin C x x dx x x y +--=-=''? 两边再积分，得()2 1 1 2cos sin 2sin cos C x C x x dx C x x y +++-=+--=? 第三次积分，得()322 121sin cos 2cos sin C x C x C x x dx C x C x x y ++++=+++-= ? 所以，原方程的通解为322 1sin cos C x C x C x x y ++++=，其中321C C C 、、为常数． 例5-8 求微分方程0='-''y y x 的通解． 解 令)(x p y ='，则)(x p y '=''．原方程可化为0=-'p p x ，即01 =-'p x p ．这是关于p p ',的一阶线性齐次微分方程．其通解为： x C e C e C x p x dx x 1ln 11 1222)(==?=，即x C y 12='．两边积分，即得原方程通解
第八章 常微分方程的数值解法 一．内容要点 考虑一阶常微分方程初值问题:?????==0 0)() ,(y x y y x f dx dy 微分方程的数值解：设微分方程的解y (x )的存在区间是[a,b ]，在[a,b ]内取一系列节 点a= x 0< x 1<…< x n =b ,其中h k =x k+1-x k ；(一般采用等距节点，h=(b-a)/n 称为步长)。在每个节点x k 求解函数y(x)的近似值:y k ≈y(x k )，这样y 0 , y 1 ,...,y n 称为微分方程的数值解。 用数值方法，求得f(x k )的近似值y k ，再用插值或拟合方法就求得y(x)的近似函数。 (一)常微分方程处置问题解得存在唯一性定理 对于常微分方程初值问题:?????==0 0)() ,(y x y y x f dx dy 如果： (1) 在B y y A x x 00≤-≤≤,的矩形内),(y x f 是一个二元连续函数。 (2) ),(y x f 对于y 满足利普希茨条件，即 2121y y L y x f y x f -≤-),(),(则在C x x 0≤≤上方程?????==0 0)() ,(y x y y x f dx dy 的解存在且唯一，这里C=min((A-x 0),x 0+B/L),L 是利普希茨常数。 定义：任何一个一步方法可以写为),,(h y x h y y k k k 1k Φ+=+，其中),,(h y x k k Φ称为算法的增量函数。 收敛性定理：若一步方法满足： (1)是p 解的. (2) 增量函数),,(h y x k k Φ对于y 满足利普希茨条件. (3) 初始值y 0是精确的。则),()()(p h O x y kh y =-kh =x -x 0,也就是有 0x y y lim k x x kh 0h 0 =--=→)( (一)、主要算法 1．局部截断误差 局部截断误差：当y(x k )是精确解时，由y(x k )按照数值方法计算出来的1~ +k y 的误差y (x k+1)- 1~ +k y 称为局部截断误差。 注意：y k+1和1~ +k y 的区别。因而局部截断误差与误差e k +1=y (x k +1) -y k +1不同。 如果局部截断误差是O (h p+1),我们就说该数值方法具有p 阶精度。
常微分方程数值解法 【作用】微分方程建模是数学建模的重要方法，因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题，大体上可以按以下几步： 1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。 2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。 3. 运用这些规律列出方程和定解条件。基本模型 1. 发射卫星为什么用三级火箭 2. 人口模型 3. 战争模型 4. 放射性废料的处理通常需要求出方程的解来说明实际现象，并加以检验。如果能得到解析形式的解固然是便于分析和应用的，但是我们知道，只有线性常系数微分方程，并且自由项是某些特殊类型的函数时，才可以得到这样的解，而绝大多数变系数方程、非线性方程都是所谓“解不出来” 的于是对于用微分方程解决实际问题来说，数值解法就是一个十分重要的手段。 1. 改进Euler 法： 2. 龙格—库塔( Runge—Kutta )方法： 【源程序】 1. 改进Euler 法： function [x,y]=eulerpro(fun,x0,x1,y0,n);%fun 为函数，(xO, x1)为x 区间，yO 为初始值，n 为子 区间个数 if nargin<5,n=5O;end h=(x1-xO)/n; x(1)=xO;y(1)=yO; for i=1:n x(i+1)=x(i)+h; y1=y(i)+h*feval(fun,x(i),y(i)); y2=y(i)+h*feval(fun,x(i+1),y1); y(i+1)=(y1+y2)/2; end 调用command 窗口 f=i nlin e('-2*y+2*x A2+2*x') [x,y]=eulerpro(f,O,,1,1O) 2 x +2x , (0 < x < , y(0) = 1 求解函数y'=-2y+2 2. 龙格—库塔( Runge—Kutta )方法： [t,y]=solver('F',tspan ，y0) 这里solver为ode45, ode23, ode113，输入参数F是用M文件定义的微分方程y'= f (x, y)右端的函数。tspan=[t0,tfinal]是求解区间，y0是初值。 注：ode45和ode23变步长的，采用Runge-Kutta算法。 ode45表示采用四阶-五阶Runge-Kutta算法，它用4阶方法提供候选解，5阶方法控制误差，是一种自适应步长(变步长)的常微分方程数值解法，其整体截断误差为(△ 口人5解 决的是Nonstiff(非刚性)常微分方程。