常微分方程第一章

常微分方程讲解常微分方程第一章绪论在初等数学中，我们已经学过一些代数方程（如元个一次联立方程），并且用它们解决了一些有趣的应用问题，使我们初步体会到方程论（主要是设未知量、列方程和求解方程的方法）对于解决实际问题的重要性。

在解析几何与微积分中，我们又碰到一类不同的方程——方程的个数少于未知量的个数，也就是通常所说的函数方程。

例如，1） (设是自变量，则是未知函数)；2），（设是自变量，则和是两个未知函数）。

这类函数方程与开头所说的代数方程相比，在概念上进了一步——确定自变量与因变量之间的函数关系。

利用这类方程可以解决一类新的问题，例如某些轨迹问题和极值问题等。

本课程所要讲述的方程与刚才说的那种函数方程又不一样，它们除了自变量和未知函数外，还包含了未知函数的导数(即微商)。

例如:1）（是自变量，是未知函数，是未知函数对的导数。

）2）（是自变量，是未知函数，是未知函数对的导数等等）。

这种联系着自变量、未知函数以及未知函数的导数（或微分）的关系式,数学上称之为微分方程。

其中未知函数的导数或微分是不可缺少的。

下面我们通过几个具体的例子，粗略地介绍常微分方程的一些物理背景和方程的建立问题，并讲述一些最基本的概念。

第一节微分方程：某些物理过程的数学模型在这一节中列举几个简单的实际例子，说明怎样从实际问题列成微分方程的问题。

例子虽然简单，但是从中能够简明地诱导出微分方程的一些基本概念，成为进一步探讨其他较复杂问题的借鉴。

掌握好这些例子，会有助于增进我们分析问题的能力。

例1 物体冷却过程的数学模型将某物体放置于空气中，在时刻时，测量得它的温度为，10分钟后测得温度为。

我们要求决定此物体的温度和时间的关系，并计算20分钟后物体的温度。

这里我们假定空气的温度保持为。

解为了解决上述问题，需要了解有关热力学的一些基本规律。

例如，热量总是从温度高的物体向温度低的物体传导的；在一定的温度范围内（其中包括了上述问题的温度在内），一个物体的温度变化速度与这一物体的温度和其所在介质温度的差值成比例。

常微分方程课程总结第一章 绪论§1.2微分方程的基本概念（1）常微分方程偏微分方程微分方程：凡含有未知函数的导数或微分的方程叫微分方程。

常微分方程：未知函数为一元函数的微分方程。

()(),dyaxy a dxdy p x y Q x dx=+=为常数 偏微分方程：未知函数为多元函数，从而出现偏导数的微分方程。

()22,22242u uf x y x y u u y x ∂∂+=∂∂∂∂=∂∂（2）线性与非线性一般n 阶线性微分方程具有形式：（等式左面全是一次有理整式）()(1)11()()()().n n n n y a x y a x y a x y f x --'++++=（3）解和隐式解微分方程的解：代入微分方程能使方程成为恒等式的函数. 隐式解：Φ（x,y ）=0 （4）通解和特解通解：微分方程的解中含有任意常数,且任意常数的个数与微分方程的阶数同.） 特解： 确定了通解中任意常数以后的解. 初始条件：用来确定任意常数的条件.初值问题: 求微分方程满足初始条件的解的问题.（5）积分曲线：微分方程任一特解的图形都是一条曲线,称为微分方程的积曲线。

第二章 一阶微分方程的初等解法§2.1 变量分离方程与变量变换2.1.1、变量分离方程)()(y x f dxdyϕ= ⎰⎰+=c dx x f y dy )()(ϕ 2.1.2、可化为变量分离方程的类型1.形如)(x y g dx dy =，称为齐次微分方程，令u ＝xy ,即y ＝ux ，于是dx dy ＝x dx du ＋u ，代入原方程，变形为x dx du ＋u ＝g （u ），整理得dx du ＝xuu g -)(2.形如222111c x b x a c x b x a dx dy ++++= 的方程也可经变量变换化为变量分离方程（1）常数）(212121k c c b b a a ===，方程化为dxdy ＝k ，有通解c kx y += （2）≠==k b b a a 212121c c 情形，令u ＝y b x a 21+，这时有dx du ＝dx dy b a 22+＝2122c u c ku b a +++是分离变量方程 （3）2121b b a a ≠情形，若21c c 、不全为零，方程右端分子、分母都是x 、y 的一次多项式，因此111c x b x a ++＝0，222c y b x a ++＝0，交点（),βα，令X ＝x －α，Y ＝y －β，化为011=+Y b X a ， 022=+Y b X a 。

常微分方程第二版答案第一章【篇一：常微分方程第一章】程1.1学习目标：1. 理解微分方程有关的基本概念, 如微分方程、方程阶数、解、通解、初始条件、初值问题等的定义和提法. 掌握处理微分方程的三种主要方法: 解析方法, 定性方法和数值方法.2. 掌握变量分离法，用变量替换将某些方程转化为变量分离方程, 掌握一阶线性方程的猜测检验法, 常数变易法和积分因子法, 灵活运用这些方法求解相应方程, 理解和掌握一阶线性方程的通解结构和性质.3. 能够大致描述给定一阶微分方程的斜率场, 通过给定的斜率场描述方程解的定性性质; 理解和掌握欧拉方法, 能够利用欧拉方法做简单的近似计算.4. 理解和掌握一阶微分方程初值问题解的存在唯一性定理, 能够利用存在唯一性定理判别方程解的存在性与唯一性并解决与之相关的问题, 了解解对初值的连续相依性和解对初值的连续性定理, 理解适定性的概念.5. 理解自治方程平衡点, 平衡解, 相线的概念, 能够画出给定自治方程的相线, 判断平衡点类型进而定性分析满足不同初始条件解的渐近行为.6. 理解和掌握一阶单参数微分方程族的分歧概念, 掌握发生分歧的条件, 理解和掌握各种分歧类型和相应的分歧图解, 能够画出给定单参数微分方程族的分歧图解, 利用分歧图解分析解的渐近行为随参数变化的状况.7. 掌握在给定的假设条件下, 建立与实际问题相应的常微分方程模型, 并能够灵活运用本章知识进行模型的各种分析.1.2基本知识： (一) 基本概念1. 什么是微分方程:联系着自变量、未知函数及它们的导数（或微分）间的关系式（一般是指等式），称之为微分方程. 2. 常微分方程和偏微分方程:(1) 如果在微分方程中，自变量的个数只有一个，则称这种微分方程为常微分方程,dy2dyd2ydy()?t?y?0. ?b?cy?f(t)例如 , dtdtdtdt2(2) 如果在微分方程中，自变量的个数为两个或两个以上，则称这种微分方程为偏?2t?t?2t?2t?2t?4微分方程. 例如 , . ???02222?t?x?x?y?z本书在不特别指明的情况下, 所说的方程或微分方程均指常微分方程.3. 微分方程的阶数: 微分方程中出现的未知函数最高阶导数的阶数.例如，d2ydy?b?cy?f(t) 是二阶常微分方程； 2dtdt?2t?t?2t?2t?2t?4与是二阶偏微分方程. ???02222?t?x?x?y?z4. n阶常微分方程的一般形式：dydnyf(t,y,,...,n)?0,dtdtdydnydydnydnyn)是t,y,,...,n的已知函数，而且一定含有n的这里f(t,y,dtdtdtdtdt 项；y是未知函数，t是自变量. 5. 线性与非线性：dydnydydny,...,n)?0的左端是y及,...,n的一次有理式，（1）如果方程f(t,y,dtdtdtdtdydny,...,n)?0为n阶线性微分方程. 则称f(t,y,dtdt（2）一般n阶线性微分方程具有形式：dnydn?1ydy?a(t)?...?a(t)?an(t)y?f(t)1n?1nn?1dtdtdt这里a1(t),…, an(t),f(t)是t的已知函数.（3）不是线性方程的方程称为非线性方程. （4）举例：d2ydy?cy?f(t)是二阶线性微分方程；方程2?bdtdtd2?g方程2?sin??0是二阶非线性微分方程；ldt方程(dy2dy)?t?y?0是一阶非线性微分方程. dtdt6. 解和隐式解：dydny,...,n)?0后，能使它变为恒等式，则如果将函数y??(t)代入方程f(t,y,dtdt）?0决定的隐函数y??(t)是称函数y??(t)为方程的解. 如果关系式?（t,y方程的解，则称?（t,y）?0为方程的隐式解. 7. 通解与特解：把含有n个独立的任意常数c1,c2,...,cn的解 y??(t,c1,c2,...,cn)称为n阶方程dydnyf(t,y,,...,n)?0的通解. 其中解对常数的独立性是指，对?及其 n?1阶导数dtdtd?dn?1?,...,n?1关于n个常数 c1,c2,...,cn的雅可比行列式不为0, 即 dtdt ???c1????c1???(n?1)?c1???c2????c2???(n?1)?c2??????cn????cn??0.??(n?1)??cn为了确定微分方程一个特定的解，通常给出这个解所必须满足的条件，称为定解条件.dydny,...,n)?0的初始条件是常见的定解条件是初始条件, n阶微分方程f(t,y,dtdtdydn?1y(1)(n?1)?y0,...,n?1?y0指如下的n个条件：t?t0，y?y0,，这里dtdt(1)(n?1)是给定的n+1个常数. 求微分方程满足定解条件的解，就是所谓t0,y0,y0,...,y0定解问题. 当定解条件为初始条件时，相应的定解问题称为初值问题. 把满足初始条件的解称为微分方程的特解. 初始条件不同，对应的特解也不同.(二) 解析方法1．变量分离方程形如dy?f(t)?(y)的方程为变量分离方程，其中f(t),?(y)分别为t,y的连续函数.dt方程解法如下：若?(y)?0，则dy?f(t)dt?(y)dy??(y)??f(t)dt?c上式确定方程的隐式通解. 如果存在y0，使得??y0??0，则y?y0也是方程的解. 2. 可化为变量分离方程的方程(1) 齐次方程dyy?g()的方程为齐次方程，g?u?为u的连续函数. dttydydu?t?u，从而原方程变为解法如下：做变量替换u?，即y?ut，有tdtdtdudug(u)?ut?u?g(u)，整理有?，此为变量分离方程，可求解. dtdtt形如 (2) 形如dya1t?b1y?c1的方程, 其中a1??a2,?b1,?b2,?c1,?c2为常数. ?dta2t?b2y?c2?a1b1c1???k的情形. a2b2c2此时方程化为dy?k,可解得y?kt?c. dt?a1a2b1b2?0,即a1b1??k的情形: a2b2ku?c1dudy?a2?b2?a2?b2dtdtu?c2令 u?a2t?b2y, 则有此为变量分离方程. ?a1b1a2b2?0的情形y. t对c1?c2?0的情况, 直接做变量替换u?当c1,c2不全为零, 求 ? ?a1t?b1y?c1?0的解为?a2t?b2y?c2?0?t??. ??y???t?t??令 ? , 则方程组化为y?y???原方程化为3．一阶线性微分方程?a1t?by1?0. ?at?by?0?22dya1t?byy??g()的齐次方程可求解. dta2t?byt(1) 一般形式：a(t)dy?b(t)y?c(t)?0，若a(t)?0，则可写成 dtdy?p(t)y?qt(的形式). dtp(t)dtdy,?c为任意常数. ?p(t)y，通解为ce?(2) 一阶齐次线性微分方程：dtdy?p(t)y?q(t)，q(t)?0. (3) 一阶非齐次线性微分方程：dt性质1 必有零解 y?0;性质2 通解等于任意常数c与一个特解的乘积; 性质3 任意两个解的线性组合也是该微分方程的解. (5) 非齐次线性微分方程的性质性质1 没有零解;性质2 非齐次方程的解加上对应齐次方程的解仍为非齐次方程的解; 性质3 任意两个非齐次方程的解的差是相应齐次方程的解. (6) 一阶非齐次线性微分方程的解法：(i) 猜测-检验法对于常系数的情形，即 p(t) 为常数, 此时方程为(4) 齐次线性微分方程的性质dy?ay?q(t), a为常数. dt对应齐次方程的通解为ce, 只需再求一个特解, 这时根据q(t)为特定的函数,bt可猜测不同的形式特解. 事实上, 当q(t)?ae, a,b为给定常数, 且b?a 时at可设待定特解为ce, 而当b?a时, 可设特解形式为cte, 后代入方程可确定待定常数c. 当q(t)为cosat,??sinat或它们的线性组合时, 其中a为给定常数. 这时可设待定特解为bcosat?csinat代入方程后确定b,?c的值. 当btbtq(t)具有多项式形式a0tn?a1tn?1???an?1t?an, 其中a0,?a1,??an 为给定常数且a0?0, 这时可设待定特解为b0t?bt1nn?1???bn?1t?bn代入方程可求得bi,?i?0,1?,??,n的值. 对于q(t)有上述几种线性组合的形式, 则可设待定特解是上述形式特解的线性组合. (ii) 常数变易法: 令y?c(t)e?p(t)dt，代入方程，求出c(t)后可求得通解为【篇二：常微分课后答案2.1】>1.dy?2xy,并求满足初始条件：x=0,y=1的特解. dx解：对原式进行变量分离得1dy?2xdx,两边同时积分得：lny?yc?1,故它的特解为y?ex。

常微分⽅程第⼀章第⼀章⼀阶微分⽅程1.1学习⽬标：1. 理解微分⽅程有关的基本概念, 如微分⽅程、⽅程阶数、解、通解、初始条件、初值问题等的定义和提法. 掌握处理微分⽅程的三种主要⽅法: 解析⽅法, 定性⽅法和数值⽅法.2. 掌握变量分离法，⽤变量替换将某些⽅程转化为变量分离⽅程, 掌握⼀阶线性⽅程的猜测检验法, 常数变易法和积分因⼦法,灵活运⽤这些⽅法求解相应⽅程, 理解和掌握⼀阶线性⽅程的通解结构和性质.3. 能够⼤致描述给定⼀阶微分⽅程的斜率场, 通过给定的斜率场描述⽅程解的定性性质; 理解和掌握欧拉⽅法, 能够利⽤欧拉⽅法做简单的近似计算.4. 理解和掌握⼀阶微分⽅程初值问题解的存在唯⼀性定理, 能够利⽤存在唯⼀性定理判别⽅程解的存在性与唯⼀性并解决与之相关的问题, 了解解对初值的连续相依性和解对初值的连续性定理, 理解适定性的概念.5. 理解⾃治⽅程平衡点, 平衡解, 相线的概念, 能够画出给定⾃治⽅程的相线, 判断平衡点类型进⽽定性分析满⾜不同初始条件解的渐近⾏为.6. 理解和掌握⼀阶单参数微分⽅程族的分歧概念, 掌握发⽣分歧的条件, 理解和掌握各种分歧类型和相应的分歧图解, 能够画出给定单参数微分⽅程族的分歧图解, 利⽤分歧图解分析解的渐近⾏为随参数变化的状况.7. 掌握在给定的假设条件下, 建⽴与实际问题相应的常微分⽅程模型, 并能够灵活运⽤本章知识进⾏模型的各种分析.1.2基本知识： (⼀) 基本概念1. 什么是微分⽅程:联系着⾃变量、未知函数及它们的导数（或微分）间的关系式（⼀般是指等式），称之为微分⽅程. 2. 常微分⽅程和偏微分⽅程:(1) 如果在微分⽅程中，⾃变量的个数只有⼀个，则称这种微分⽅程为常微分⽅程,例如 )(22t f cy dt dy b dt y d =++, 0)(2=++y dtdyt dt dy .(2) 如果在微分⽅程中，⾃变量的个数为两个或两个以上，则称这种微分⽅程为偏微分⽅程. 例如 0222222=??+??+??zTy T x T , t T x T ??=??422. 本书在不特别指明的情况下, 所说的⽅程或微分⽅程均指常微分⽅程. 3. 微分⽅程的阶数: 微分⽅程中出现的未知函数最⾼阶导数的阶数. 例如，)(22t f cy dt dyb dty d =++ 是⼆阶常微分⽅程； 0222222=??+??+??zTy T x T 与t T x T ??=??422是⼆阶偏微分⽅程. 4. n 阶常微分⽅程的⼀般形式：(,,,...,)0n n dy d yF t y dt dt=,这⾥(,,,...,)n n dy d y F t y dt dt 是,,,...,n n dy d y t y dt dt 的已知函数，⽽且⼀定含有n n d y dt的项；y 是未知函数，t 是⾃变量. 5. 线性与⾮线性：（1）如果⽅程(,,,...,)0n n dy d y F t y dt dt =的左端是y 及,...,n n dy d ydt dt的⼀次有理式，则称(,,,...,)0n n dy d yF t y dt dt=为n 阶线性微分⽅程. （2）⼀般n 阶线性微分⽅程具有形式：1111()...()()()n n n n n n d y d y dy a t a t a t y f t dt dt dt---++++= 这⾥1()a t ,…, ()n a t ,()f t 是t 的已知函数.（3）不是线性⽅程的⽅程称为⾮线性⽅程. （4）举例：⽅程)(22t f cy dt dyb dt y d =++是⼆阶线性微分⽅程；⽅程0sin 22=+φφl gdtd 是⼆阶⾮线性微分⽅程；⽅程0)(2=++y dtdy t dt dy 是⼀阶⾮线性微分⽅程. 6. 解和隐式解：如果将函数()y t ?=代⼊⽅程(,,,...,)0n n dy d yF t y dt dt=后，能使它变为恒等式，则称函数()y t ?=为⽅程的解. 如果关系式,0t yΦ=（）决定的隐函数()y t ?=是⽅程的解，则称,0t yΦ=（）为⽅程的隐式解. 7. 通解与特解：把含有n 个独⽴的任意常数n c c c ,...,,21的解 12(,,,...,)n y t c c c ?=称为n 阶⽅程(,,,...,)0n n dy d yF t y dt dt =的通解. 其中解对常数的独⽴性是指，对?及其 1n -阶导数11,...,n n d d dt dt--关于n 个常数 n c c c ,...,,21的雅可⽐⾏列式不为0, 即 1212(1)(1)(1)120n n n n n nc c c c c c c c c ?---??????'''??????≠??????L L M M L M L.为了确定微分⽅程⼀个特定的解，通常给出这个解所必须满⾜的条件，称为定解条件.常见的定解条件是初始条件, n 阶微分⽅程(,,,...,)0n n dy d yF t y dt dt =的初始条件是指如下的n 个条件： 1(1)(1)00001,,...,n n n dy d y t t y y y y dt dt---====，，这⾥(1)(1)0000,,,...,n t y y y -是给定的n+1个常数. 求微分⽅程满⾜定解条件的解，就是所谓定解问题. 当定解条件为初始条件时，相应的定解问题称为初值问题. 把满⾜初始条件的解称为微分⽅程的特解. 初始条件不同，对应的特解也不同.(⼆) 解析⽅法1．变量分离⽅程形如()()dyf t y dt=的⽅程为变量分离⽅程，其中(),()f t y ?分别为,t y 的连续函数.⽅程解法如下：若()0y ?≠，则()()()()dyf t dt y dyf t dt cy ??==+??上式确定⽅程的隐式通解. 如果存在0y ，使得()00y ?=，则0y y =也是⽅程的解. 2. 可化为变量分离⽅程的⽅程(1) 齐次⽅程形如 ()dy yg dt t=的⽅程为齐次⽅程，()g u 为u 的连续函数. 解法如下：做变量替换y u t =，即y ut =，有dy dut u dt dt=+，从⽽原⽅程变为()du t u g u dt +=，整理有()du g u udt t-=，此为变量分离⽅程，可求解. (2) 形如111222a tb yc dy dt a t b y c ++=++的⽅程, 其中121212,,,,a a b b c c , 为常数. ●111222a b c k a b c ===的情形. 此时⽅程化为,dyk dt=可解得y kt c =+. ●11220,a b a b =即1122a bk a b ==的情形: 令 22,u a t b y =+ 则有 122222ku c du dya b a b dt dt u c +=+=++ 此为变量分离⽅程. ●11220a b a b ≠的情形对120c c ==的情况, 直接做变量替换y u t =. 当12,c c 不全为零, 求 11122200a t b y c a t b y c ++=??++=?的解为t y αβ=??=?. 令 T t Y y αβ=-??=-?, 则⽅程组化为112200a T bY a T b Y +=??+=?. 原⽅程化为12()a T bY dY Yg dT a T bY T+==+的齐次⽅程可求解. 3．⼀阶线性微分⽅程(1) ⼀般形式：()()()0dya tb t yc t dt++=，若()0a t ≠，则可写成()()dyP t y Q t dt=+的形式. (2) ⼀阶齐次线性微分⽅程：()dyP t y dt =，通解为(),P t dt ce c ? 为任意常数.(3) ⼀阶⾮齐次线性微分⽅程：()()dyP t y Q t dt=+，()0Q t ≠.(4) 齐次线性微分⽅程的性质性质1 必有零解 0y =;性质2 通解等于任意常数c 与⼀个特解的乘积; 性质3 任意两个解的线性组合也是该微分⽅程的解. (5) ⾮齐次线性微分⽅程的性质性质1 没有零解;性质2 ⾮齐次⽅程的解加上对应齐次⽅程的解仍为⾮齐次⽅程的解; 性质3 任意两个⾮齐次⽅程的解的差是相应齐次⽅程的解.(6) ⼀阶⾮齐次线性微分⽅程的解法：(i) 猜测-检验法对于常系数的情形，即 ()P t 为常数, 此时⽅程为()dyay Q t dt=+, a 为常数. 对应齐次⽅程的通解为atce , 只需再求⼀个特解, 这时根据()Q t 为特定的函数,可猜测不同的形式特解. 事实上, 当()BtQ t Ae =, ,A B 为给定常数, 且B a ≠时可设待定特解为BtCe , ⽽当B a =时, 可设特解形式为BtCte , 后代⼊⽅程可确定待定常数C . 当()Q t 为cos ,sin At At 或它们的线性组合时, 其中A 为给定常数. 这时可设待定特解为cos sin B At C At +代⼊⽅程后确定,B C 的值. 当()Q t 具有多项式形式1011n n n n a t a t a t a --++++L , 其中01,,n a a a L 为给定常数且00a ≠, 这时可设待定特解为1011n n n n b t b t b t b --++++L 代⼊⽅程可求得,0,1,,i b i n = L 的值. 对于()Q t 有上述⼏种线性组合的形式, 则可设待定特解是上述形式特解的线性组合. (ii) 常数变易法: 令()()P t dty c t e ?=，代⼊⽅程，求出()c t 后可求得通解为()()(())P t dtP t dty e Q t e dt c -?=+.(iii) 积分因⼦法: ⽅程改写为()()dyP t y Q t dt-=, 将()P t dt e µ-?=, 乘⽅程两端得 ()()()()()P t dt P t dtP t dt dy e e P t y Q t e dt---?-= 即 ()()()()P t dtP t dt d ye Q t e dt--?=, 从⽽通解为 ()()()P t dt P t dt ye Q t e dt c --?? =+?，即 ()()(())P t dt P t dt y e Q t e dt c -??= +?.注意, ⾮齐次线性微分⽅程通解的结构是: ⾮齐次线性微分⽅程的通解等于其对应的齐次线性微分⽅程的通解加上⾮齐次线性微分⽅程的⼀个特解.4. 伯努利(Bernoulli)⽅程. 形如()()n dyP t y Q t y dt=+的⽅程, 其中 n 是常数且0,1,(),()n P t Q t ≠ 是连续函数, 称为伯努利⽅程. 伯努利⽅程可通过变量替换 1nz y-=化为(1)()(1)()dyn P t z n Q t dt=-+-, 这是关于未知函数z 的线性⽅程, 可求其通解.(三) 定性⽅法与数值⽅法：1. 斜率场：⼀阶微分⽅程(,)dyf t y dt =的解()y t ?=代表ty 平⾯上的⼀条曲线，称之为微分⽅程的积分曲线. 微分⽅程(,)dyf t y dt=的通解()y t ?=，c 对应于ty 平⾯上的⼀族曲线，称之为微分⽅程的积分曲线族. 满⾜初始条件00()y t y =的特解就是通过点00(,)t y 的⼀条积分曲线. ⽅程(,)dy f t y dt =的积分曲线上的每⼀点(,)t y 处的切线斜率dydt刚好等于函数(,)f t y 在这点的值. 也就是，积分曲线的每⼀点(,)t y 以及这点上的切线斜率dydt恒满⾜⽅程；反之，如果在⼀条曲线每点上其切线斜率刚好等于函数(,)f t y 在这点的值，则这⼀条曲线就是⽅程的积分曲线. 这样，可以⽤(,)f t y 在ty 平⾯的某个区域D 内定义过各点的⼩线段，其斜率为(,)f t y ，⼀般称这样的⼩线段为斜率标记. ⽽对ty 平⾯上D 内任⼀点(,)t y , 有这样⼀个⼩线段与之对应, 这样在D 内形成⼀个⽅向场, 称为斜率场. 斜率场是⼏何直观上描述解的常⽤⽅法2. 欧拉⽅法：求微分⽅程初值问题00(,)()dyf t y dty t y== 的解，可以从初始条件00()y t y =出发，按照⼀定的步长t ? 依照某种⽅法逐步计算微分⽅程的近似解()n n y y t =, 这⾥0n t t n t =+?这样求出的解称为数值解. 利⽤欧拉公式10(,),n n n n n y y f t y t t t n t +=+? =+?,可求初值问题的近似解，这种⽅法称为欧拉⽅法.欧拉⽅法具有⼀阶误差精度 .如果我们先⽤欧拉公式求出近似解,再利⽤梯形公式进⾏校正, 得到的近似解将具有2阶误差精度,具体为预测: 1(,)n n n n y y f t y t +=+?,校正: 11,11[(,)()]2n n n n n n y y f t y f t y t ++ +=++?, 这种⽅法称为改进的欧拉⽅法.(四) 解的存在性、唯⼀性及解对初值的连续相依性1. 利普希茨（lipschitz ）条件: 函数(,)f t y 称为在区域2D ?R 内关于y 满⾜利普希茨条件，是指如果存在常数0L >，使得不等式1212(,)(,)f t y f t y L y y -≤-对于所有的12(,),(,)t y t y D ∈都成⽴, 其中L 称为利普希茨常数. 2. 基本定理(1) 解的存在性定理: 设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << < 和函数()y t , 定义于区间00(,)t t εε-+内，是初值问题00(,)()dyf t y dt y t y== 的解. (2) 解的唯⼀性定理: 设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <且关于y 满⾜利普希茨条件. 如果00(,)t y D ∈并且12(),()y t y t 是初值问题00(,)()dyf t y dt y t y==在区间00(,)t t εε-+内的两个解，那么对任意的00(,)t t t εε∈-+,12()()y t y t =,即解是唯⼀的.注记1: 存在性定理和唯⼀性定理结合在⼀起称为初值问题解的存在唯⼀性定理，叙述如下：设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << < 和函数()y t , 定义于区间00(,)t t εε-+内，是初值问题00(,)()dyf t y dt y t y== 的唯⼀解. 因⽽当我们判断初值问题解的存在唯⼀性时，要检查(,)f t y 需要满⾜的条件.注记2: 由于利普希茨条件较难检验，常⽤(,)f t y 在2{(,):,}D t y a t b c y d =∈ ≤≤ ≤≤R上对y 有连续偏导数来代替. 事实上，如果在D 上y f ??存在且连续，则yf在D 上有界. 设在D 上L yf≤??, 这时 2121212(,())(,)(,)f t y y y f t y f t y y y yθ?+--=-?21y y L -≤,其中 12(,),(,),01t y t y D θ∈ <<. 但反过来满⾜利普希茨条件的函数(,)f t y 不⼀定有偏导数存在. 例如(,)||f t y y = 在任何区域内都满⾜利普希茨条件，但它在0y =处没有导数.(3) 解对初值的连续相依性定理设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <条件. 如果00(,)t y D ∈，00(,,)y t t y ?=是初值问题00(,)()dyf t y dt y t y==在区间00(,)t h t h -+内的解，其中 0h >,那么，对任意给定的0>ε，必能找到正数(,)0h δδε=>，使得当2220000t t y y δ-+-<（）（）时，初值问题00(,)()dyf t y dt y t y==的解00(,,)y t t y ?=在区间00(,)t h t h -+内也有定义，并且0000|(,,),,|,t t y x t y ??ε-<（） 00(,)t t h t h ∈-+. (4) 解对初值的连续性定理设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <条件. 如果00(,)t y D ∈，00(,,)y t t y ?=是初值问题00(,)()dyf t y dt y t y==的解, 那么00(,,)t t y ?作为00,,t t y 的三元函数在它存在的范围内是连续的.3. 初值问题的适定性当⼀个微分⽅程初值问题的解存在, 唯⼀并且解连续的依赖于初始条件时, 我们称该问题是适定的. 那么, 对于常微分⽅程初值问题00(,)()dyf t y dt y t y==, 只要在00(,)t y 所在的区域内,(,)f t y 连续并且关于y 满⾜利普希茨条件, 则该初值问题是适定的.(五) ⾃治⽅程的平衡点与相线1. ⾃治⽅程当⼀阶微分⽅程(,)dy f t y dt =的右端项只是y 的函数⽽与⾃变量t ⽆关, 即()dy f y dt=时, 称为⾃治⽅程.2. 平衡解与平衡点对⾃治⽅程()dyf y dt=⽽⾔, 若()0f y =有解0y y =, 则称 0()y t y ≡是⽅程的平衡解, ⽽点0y 称为⽅程的⼀个平衡点. 3. 相线相线是仅仅对⾃治⽅程()dyf y dt=⽽⾔的⼀种简化的斜率场. ⾃治⽅程的斜率场在⽔平直线上的斜率标记是⼀样的, 这样只要知道⼀条竖直直线上的斜率标记,我们就可以知道整个斜率场. 因⽽, 在⼀个竖直的直线上, 我们⽤向上的箭头表⽰正的导数, ⽤向下的箭头表⽰负的导数. 对于导数为零的点, ⽤实⼼圆点来标记它, 则形成该⾃治⽅程的相线. 4. 画相线的基本步骤 (1) 画出y -线(竖直线),(2) 找到并在y -线上标记平衡点,不连续点或定义域外的点 (3) 找到()0f y >的区间, 在这些区间上画上向上的箭头, (4) 找到()0f y <的区间, 在这些区间上画上向下的箭头.5. 初值问题0(),(0)dyf y y y dt= =解的渐近⾏为 (1) 趋向于平衡点, 如01()(1),2f y y y y =- =;(2) 在⽆限时间内趋于⽆穷, 如0(),1f y y y = =;(3) 在有限时间内趋于⽆穷(爆破), 如20(),1f y y y = =;(4) 在有限时间内停⽌(导数趋于⽆穷), 如 01(),1f y y y=- =. 6. 平衡点的分类对于⾃治⽅程()dyf y dt=, 如果()f y 在(,)-∞+∞ 内连续, 那么它的解当t 增加时要么(在有限或⽆限时间⾥)趋于+∞或-∞, 要么渐近趋于平衡点. 因⽽,平衡点在⾃治⽅程的研究中起着重要的作⽤. (1) 汇对于初值接近0y 的解, 当t 增加时, 都渐近趋于0y . 对于这样的平衡点0y , 我们称之为汇, 它是稳定的. (2) 源对于初值接近0y 的解, 当t 增加时, 都远离0y . 对于这样的平衡点0y , 我们称之为源,它是不稳定的. (3) 结点既不是源也不是汇的平衡点, 我们称之为结点,它也是不稳定的. 7. 判断平衡点类型的线性化⽅法 1. 如果0y 是⾃治⽅程()dyf y dt=的⼀个平衡点, 即0()0f y =, 那么 (1) 0y 是源当且仅当()f y 在0y 附近严格单调增加; (2) 0y 是汇当且仅当()f y 在0y 附近严格单调递减. 2. (线性化定理) 如果0y 是⾃治⽅程()dyf y dt=的⼀个平衡点, 即0()0f y =, 并且()f y 是连续可微的, 那么 (1) 若0()0f y '> 则0y 是源; (2) 若0()0f y '<, 则0y 是汇;(3) 若0()0f y '=, 则需要进⼀步的信息决定其类型.(六) 分歧⼀阶微分⽅程解的渐近⾏为随参数变化发⽣了类型的变化, 我们称之为分歧现象(或分⽀, 分叉).1. 分歧发⽣的条件对于单参数微分⽅程族()(,)dyf y f y dtµµ==, 0µµ=是⼀个分歧值的必要条件是: 存在平衡点0y , 使得 0000(,)(,)0ff y y yµµ?==?. 这样我们要找分歧点可以通过求解⽅程组 (,)0(,)0f y fy y µµ=??=?, 得到解 00(,)y µ,0µ为可能的分歧值, ⽽0y 是可能发⽣分歧的平衡点. 2. 分歧图解与分歧类型分歧图解是y µ 平⾯上⽅程在分歧值附近的所有相线的图, ⽤以强调当参数经过分歧值时相线所经历的变化. (1) 鞍结点分歧在分歧图解(图1-1)中, 当µ从左到右经过分歧值0µ时, ⽅程的平衡点从两个变为⼀个再变为不存在, 这种分歧⼀般称之为鞍结点分歧. 这类分歧图解在分歧值附近是抛物线的形状(2) 在分歧图解(图1-2)中,当µ从右到左经过分歧值0µ=时, ⽅程的平衡点由三个变为⼀个, 这种分歧⼀般称之为⾳叉分歧.图 1-1 鞍结点分歧图 1-2 ⾳叉分歧图 1-3 跨越分歧图 1-4 复合分歧(3) 在分歧图解(图1-3)中, 当0µ= 时, ⽅程有⼀个平衡点; 当0µ≠ 时, ⽅程有两个平衡点. 0µ=是⼀个分歧值. 虽然在分歧值的两侧⽅程都有两个平衡点，但平衡点的稳定性会改变. 当0µ> 时, 0y =是⼀个汇,它是稳定的; 当0µ<时, 0y =是⼀个源,它是不稳定的.这类分歧⼀般称为跨越分歧.(4) 在分歧图解(图1-4)中, 当 µ从左到右变化时,相应的⽅程平衡点依次由⼀个变为两个,三个,两个再变回⼀个, 这种分歧⼀般称之为复合分歧.(七) ⼀阶微分⽅程的应⽤1. 增长和衰减问题设 ()S t 为正在增长或衰减的某研究对象的总量. 如果假设它随时间的变化率dS dt与当前数⽬成正⽐, 其⽐例系数为 k , 则有dS kS dt =, 或 0dS kS dt-=. 设()S t 可微, 因⽽是连续函数. Malthus ⼈⼝模型满⾜上述微分⽅程, 虽然对⼈⼝问题,()S t 是离散的, 只能取整数值, 但该模型系统在⼀定情况下提供了很好的近似对某⼀⽣物种群进⾏研究时, 该⽣物种群的增长往往受资源和环境的限制, 引进参量N , 称为最⼤承载量, ⽤以表⽰⾃然资源和环境条件所能容纳的最⼤数量, 并且假定（1）当基数很⼩时，增长率与当前数成正⽐；（2）当基数很⼤，达到资源和环境不能承受的时候，数量开始减少，即增长率为负的.此时⽅程可改写为(1)dS Sk S dt N=-, 称为具有增长率k 和最⼤承载量N 的Logistic 模型，该模型最早由荷兰⽣物学家 Verhulst在1838年提出.2. 温度问题⽜顿冷却定律(亦适应于加热的情况)说明物体的温度随时间的变化率与物体所处的周围环境的温差成正⽐, 设 T 是物体的温度, T 是所处环境的温度, 那么物体温度随时间的变化率为dTdt, ⽜顿冷却定律可表⽰为 ()dTk T T dt=--, 其中k 是正的⽐例系数, ⽽负号表⽰在冷却过程中, 物体温度 T ⼤于周围环境温度T , 变化率0dT dt <. 在加热过程中0dT dt>, 此时T T <. 3. 稀释问题⼀容器最初容纳0V 升盐⽔溶液, 其中含盐 a 克. 每升含盐 b 克的盐⽔溶液以e 升/分的速度注⼊,同时, 搅拌均匀的溶液以f 升/分的速度流出, 问在任何时刻 t , 容器中的含盐量.设Q 为任何时刻容器中的含盐量. Q 的变化率dQdt等于盐的注⼊率减去流出率. 盐的注⼊率是 be 克/分. 要决定流出率, ⾸先计算在时刻t , 容器中的溶液的体积, 它等于最初的体积0V 加上注⼊的体积 et 后减去流出的体积ft . 因此, 在任⼀时刻t , 盐⽔的体积是 0V et ft +-. 在任何时刻的浓度是0Q V et ft +-, 由此得流出率为 0QfV et ft+-/分.于是得到微分⽅程0dQ Qf be dt V et ft =-+-, 即 0dQ fQ be dt V et ft+=+-, 这是⼀个⼀阶线性⽅程.4. 电路⼀个简单的 RC 回路是包含有电阻R (欧姆), 电容C (法拉)和电源V (伏特),如图1-5.图1-5 RC 电路图1-6 RL 电路由电路学知识，C 的电压()v t 与电阻R 的电压之和应为电源的电压()V t . 电路中的电流I (安培)为 ()dQ dCv t dv I Cdt dt dt ===, 其中 Q 为电量从⽽R 处的电压为 dvRI RC dt=, 由此我们可以建⽴RC 电路的模型如下：()dv RC v V t dt +=, 即 ()dv V t vdt RC-=. 对于⼀个包含有电阻R (欧姆), 电感L (亨利)和电源V (伏特)的RL 回路,如图1-6. 电路中的电流应满⾜的基本⽅程为 dI R VI dt L L+=.(⼋) 种群⽣态学中的模型设()y t 表⽰⼀个⽣物种群的数量, t 为时间, 最简单的种群模型是 Malthus 模型dyky dt=. Malthus 模型的解()(0)kty t y e =预测了种群数量的指数增长.由于种群数量⼤的时候，对资源的竞争加剧，因此单位增长率会随种群数⽬增⼤⽽减⼩，因此更为合理的假设是()dyyf y dt = (*) 这⾥()f y 是单位增长率，因为dy dt 为增长率，y 是种群数量, ⽽()/dyf y y dt=. 当考虑种群数量的变化时.对()f y ⽽⾔, 其代数形式并不重要, ⽽关键是其单调性, 凸凹性, 这样我们可以对其进⾏⼤致分类:(1) 若()f y 在[0,)+∞上是递减的，称(*)为 Logistic 型； (2) 若()f y 在[0,)+∞上是先增后减的，称(*)为 Allee 效应型；(3) 若()f y 在[0,)+∞上是递减再递增最后递减的，称(*)为 Hysteresis 型.1.3典型例题：例1考虑微分⽅程3220dyy y y dt=--, 问 (1) y 为何值时, ()y t 将保持不变?(2) y 为何值时, ()y t 将增加?(3) y 为何值时, ()y t 将减少?解: 因为当0dy dt =时, ()y t 将保持不变; 当0dy dt >时, ()y t 将增加; 当0dy dt<时, ()y t 将减少. 由3220dyy y y dt=--知,(1) 当32200y y y --=, 即0,4,5y y y = =-=时, ()y t 将保持不变. (2) 当32200y y y -->, 即40y -<< 或5y > 时, ()y t 将增加. (3) 当32200y y y --<, 即4y <- 或05y << 时, ()y t 将减少. 例2假定在鄱阳湖中⼀种鱼类的数量()S t 随时间的变化按Logistic 模型增长, 增长率为k , 最⼤承载量为N , 即有(1)dS Sk S dt N=-. 如果每年要从湖中捕获⼀定量的鱼, 试按下述不同情形对模型做适当修改,(1) 每年捕获10吨?(2) 每年捕获总量的三分之⼀? (3) 捕获量与总量的平⽅根成正⽐?解: (1)(1)10dS Sk S dt N =--. (2) 1(1)3dS S k S S dt N =--.(3) (1)dS Sk S l S dt N=--, 其中 l 是捕获量与总量平⽅根的⽐例系数.例3求解⽅程dy tdt y=- 解：变量分离得 ydt tdy =-.两边积分 22222y t c =-+. 通解为 22t y c +=, c 为任意正常数.例4求解⽅程231dy y dx xy x y+=+ 解：变量分离得221(1)ydy dxy x x =++, 两边积分2221()1(1)1ydy dx xdx y x x x x ==-+++.即22111ln(1)ln ||ln(1)22y x x c +=-++, 1c 为任意常数, 整理得222(1)(1)y x cx ++=, 12c c e =为任意正的常数.例5求解⽅程tan dy yxy dx x -=. 解: 将⽅程改写为 tan dy y ydx x x=+, 这是齐次⽅程,做变量替换yu x =，即y ux =，有dy du x u dx dx=+，从⽽原⽅程变为tan dux u u u dx+=+ 即tan du udx x=利⽤分离变量法求得 sin u cx =, 代回原变量得通解为sin ycx x =, c 为任意常数例6 求解⽅程22dyx y x y dx=+-.解: ⽅程改写为2sgn 1()dy y y x dx x x=+?- 令u=y x ，则y ux =，从⽽2sgn 1du x u u x u dx+=+?- 当210u -≠时，2sgn 1du xdx xu =-, arcsin sgn ln u x x c =?+, 即 arcsinsgn ln yx x c x=?+, c 为任意常数．此外，还有解210u -=，即22y x =．例7求解⽅程13dy x y dx x y -+=+- 解: 解⽅程组 1030x y x y -+=??+-=?的解为12x y =??=?. 令 12X x Y y =-??=-? , 则原⽅程化为 dY X Y dX X Y -=+.令 Y u X =,则可化为变量分离⽅程 21,12dX u du X u u +=-- 解得 222Y XY X c --=, 代回原变量有22262y xy x y x c +---=, c 为任意常数. 例8求解⽅程2()dyy b t dt-=, 其中 (1) 2()1b t t t =++, (2) 4()tb t e = (3) 2()3tb t e = (4) ()cos3b t t =(5) 422()3cos31t t b t e e t t t =+++++解: 对应齐次⽅程的通解为 2ty ce =, 下⾯⽤猜测-检验法求特解(1) 设 21y At Bt C =++ 代⼊221dyy t t dt-=++, 有 2222()1At B At Bt C t t +-++=++解得 1,1,12A B C =- =- =-, 从⽽21112y t t =---, 原⽅程的通解为 22112ty ce t t =---, c 为任意常数. (2) 设 42ty Ae = 代⼊ 42t dy y e dt-=, 有44442t t t Ae Ae e -=解得 12A =, 从⽽4212ty e =, 原⽅程的通解为 2412tt y ce e =+, c 为任意常数.(3) 不能设2tAe 形式的特解, 因为它是相应齐次⽅程的解,不可能是⾮齐次⽅程的解,设 23ty Ate = 代⼊22t dyy e dt-=, 有 2222223t t t t Ate Ae Ate e +-=解得 3A =, 从⽽233ty te =, 原⽅程的通解为2223(3)t t ty ce te c t e =+=+, c 为任意常数.(4) 设 4cos3sin 3y A t B t =+ 代⼊2cos3dyy t dt-=, 有 3sin33cos32(cos3sin3)cos3A t B t A t B t t -+-+=有 2310320A B A B -+-=??--=?, 解得 23,1313A B =- =,从⽽423cos3sin 31313y t t =-+, 原⽅程的通解为 223cos3sin 31313ty ce t t =-+, c 为任意常数.(5) 根据叠加原理, 由前⾯4个⼩题知⽅程有特解422512313cos3sin 31213132t t y e te t t t t =+-+---原⽅程的通解为242212313cos3sin 31213132t t t y ce e te t t t t =++-+---,c 为任意常数.例9求⽅程22dy y dx x y =-的通解. 解: 将⽅程改写为222dx x y x y dy y y-==-. 求齐次线性微分⽅程2dx x dy y=, 得通解为2x cy =. (常数变易法) 令 2()x c y y =代⼊原⽅程得()1,()ln ||dc y c y y c dy y=- =-+, 从⽽可得原⽅程的通解为2(ln ||)x y y c =-+, c 为任意常数.例10 求⽅程26dy yty dt t=-的通解. 解: 此为 2n =的伯努利⽅程. 令 1z y -=可得6dz z t dt t=-+，此为线性⽅程可求通解为 268c t z t =-+, 代回原变量得2618c t y t =-+, 即688t t c y -=, c 为任意常数. 此外, 原⽅程还有解0y =. 例11 ⽤积分因⼦法求解⽅程32(1)1dy y t dt t =+++. 解: ⽅程改写为 32(1)1dy y t dt t -=++, 积分因⼦为 221()(1)dtt t e t µ- -+?==+,乘⽅程两端得 23(1)2(1)1dyt t y t dt--+-+=+, 即2(1)1d t y t dt -+=+, 有 421(1)(1)2y t c t =+++, c 为任意常数.例12 若()f t 连续且0()()10tf t f s ds t = , ≠?, 试求函数()f t 的⼀般表达式.解: 设0()()tF t f s ds =?, 则()F t 可导且()()F t f t '=, 这样有1,dFFFdF dt dt= =, 得 2()2,()2F t t c F t t c =+ =±+, ⼜(0)0F =, 得0c =. 从⽽ ()2F t t =±,进⽽ 1()()2f t F t t'==±. 例13 求具有性质 ()()()1()()y t y s y t s y t y s ++=- 的函数 ()y t , 已知(0)y '存在.解: ⾸先令 0s =, 由已知可得 ()(0)()1()(0)y t y y t y t y +=-,化简有 2(0)(1())0y y t +=, 知 (0)0y =. 由函数的导数定义。

微分方程的基本概念第一章常微分方程微分方程是数学理论（尤其是微积分）与实践相结合的重要途径之一。

它是研究许多自然科学、工程技术、生物技术、农业、经济和许多其他问题的有力工具。

因此，微分方程具有重要的应用价值。

本章主要介绍常微分方程的一些基本概念和几种常见微分方程的一些基本解。

下面我们通过两个具体例题来说明微分方程的基本概念。

示例一曲线通过点（1,2），曲线任意点m（x，y）的切线斜率为2x。

求出曲线方程。

解设所求的曲线方程为y=y(x),则根据导数的几何意义可知，未知函数y=y(x)应满足下面的关系：阿迪？2X，（1）DX，当x=1，y=2，也就是说，y（1）=2（2）到（1）dy?2x两端积分，得dx2y=2xdx?x?c(3)其中c是任意常数。

如果y（1）=2代，C=1代（3），即得所求曲线方程y?x?1（4）例2质量为M的粒子只有在重力的作用下才能从静止状态自由下落。

试着找到它的运动方程解在中学阶段就已经知道，从高度为h处下落的自由落体，离地面高度s的变化规律为s=h－程．二百一十二gt，其中g为重力加速度．这个规律是怎么得到的呢？下面我们给出推导过2m?取质点下落的铅垂线为s轴，它与地面的交点为原点，并规定正?h向朝上．设质点在时刻t的位置在s(t)（如图1-1）．因为质点只受方向向下的重s(t)力的作用(空气阻力忽略不计)，由牛顿第二定律f=ma，得D2S（T）M=-mg。

参见图1-11D2S（T），即=G（5）2dt根据质点由静止状态自由下降的假设，初始速度为0，所以s=s(t)还应满足下列条件s|t=0=h，Ds | t=0=0，（6）DT积分方程（6）的两边，和ds(t)=－g?dt=－gt+c1，(7)dt两边再积分，得s（t）=？（？gt？c1）dt=－12gt+c1t+c2，(8)2其中c1，c2均为任意常数．将条件（7）代入方程（8）和（9）中，得到C1=0和C2=h。

然后运动方程为s(t)=－十二gt+h．(9)2上述两个例子中的关系式（1）和（5）中，都含有未知函数的导数，自变量也都只有一个，且方程都附加有一定的条件。

第一章绪论1.1 引言常微分方程是现代数学的一个重要分支，是人们解决各种实际问题的有效工具。

微分方程的理论和方法从17世纪末开始发展起来，很快成了研究自然现象的强有力工具，在17到18世纪，在力学、天文、科学技术、物理中，就已借助微分方程取得了巨大的成就。

1864年Leverrer根据这个方程预见了海王星的存在，并确定出海王星在天空中的位置。

现在，常微分方程在许多方面获得了日新月异的应用。

这些应用也为常微分方程的进一步发展提供了新的问题，促使人们对微分方程进行更深入的研究，以便适应科学技术飞速发展的需要。

研究常微分方程常用数值解是数学工作者的一项基本的且重要的工作。

在国内外众多数学家的不懈努力，使此学科基本上形成了一套完美的体系。

微分方程的首要问题是如何求一个给定方程的通解或特解。

到目前为止，人们已经对许多微分方程得出了求解的一般方法。

由于在生产实际和科学研究中所遇到的微分方程问题比较复杂，使这些问题的解即使能求出解析表达式，也往往因计算量太大而难于求出，而对于一些典型的微分方程则可以运用基本方法求出其解析解，并可以根据初值问题的条件把其中的任意常数确定下来。

由于求通解存在许多困难，人们就开始研究带某种定解条件的特解。

首先是Cauchy对微分方程初始解的存在惟一性进行了研究。

目前解的存在惟一性、延拓性、大范围的存在性以及解对初始解和参数的延续性和可微性等理论问题都已发展成熟。

与此同时，人们开始采取各种近似方法来求微分方程的特解，例如求微分方程数值解的Euler折线法、Runge-Kutta法等，可以求得若干个点上微分方程的近似解。

最后，由于当代高科技的发展为数学的广泛应用和深入研究提供了更好的手段。

用计算机结合Matlab软件求方程的精确解、近似解，对解的性态进行图示和定性、稳定性研究都十分方便有效。

本章先介绍常微分的一般概念、导出微分方程的一些典型例子及求解微分方程的思路分析。

从而得到常微分方程的常用数值解法。