微分方程与解
第1讲 微分方程与解
微分方程
什么是微分方程?它是怎样产生的?这是首先要回答的问题.
300多年前,由牛顿(Newton,1642-1727)和莱布尼兹(Leibniz,1646-1716)所创立的微积分学,是人类科学史上划时代的重大发现,而微积分的产生和发展,又与求解微分方程问题密切相关. 这是因为,微积分产生的一个重要动因来自于人们探求物质世界运动规律的需求. 一般地,运动规律很难全靠实验观测认识清楚,因为人们不太可能观察到运动的全过程. 然而,运动物体(变量) 与它的瞬时变化率(导数) 之间,通常在运动过程中按照某种己知定律存在着联系,我们容易捕捉到这种联系,而这种联系,用数学语言表达出来,其结果往往形成一个微分方程. 一旦求出这个方程的解,其运动规律将一目了然. 下面的例子,将会使你看到微分方程是表达自然规律的一种最为自然的数学语言. 例1 物体下落问题
设质量为m 的物体,在时间t =0时,在距地面高度为H 处以初始速度v (0) = v 0垂直地面下落,求ss 此物体下落时距离与时间的关系.
解 如图1-1建立坐标系,设为t
. 于是物体下落的速度为
加速度为
质量为m 的物体,在下落的任一时刻所受到的外力有重力mg 和空气阻力,当速度不太大时,空气阻力可取为与速度成正比. 于是根据牛顿第二定律 F = ma (力=质量×加速度) 可以列出方程
其中k > 0为阻尼系数,g 是重力加速度.
(·= ) (1.1)
(1.1)式就是一个微分方程,这里t 是自变量,x 是未知函数,是未知函数对t 导
数. 现在,我们还不会求解方程(1.1),但是,如果考虑k =0的情形,即自由落体运动,此
时方程(1.1)可化为
将上式对t 积分两次得
(1.2)
其中
和
(1.3)
是两个独立的任意常数,它是方程(1.2)的解.
一般说来,微分方程就是联系自变量、未知函数以及未知函数的某些导数之间的关
系式. 如果其中的未知函数只与一个自变量有关,则称为常微分方程;如果未知函数是两个或两个以上自变量的函数,并且在方程中出现偏导数,则称为偏微分方程. 本书所介绍的都是常微分方程,有时就简称微分方程或方程. 例如下面的方程都是常微分方程
(1.4)
(1.5)
(·= ) (1.6)
(′= ) (1.7)
在一个常微分方程中,未知函数最高阶导数的阶数,称为方程的阶. 这样,一阶常
微分方程的一般形式可表为
如果在(1.8)中能将y ′解出,则得到方程
或
(1.9)
(1.10) (1.8)
(1.8)称为一阶隐式方程,(1.9)称为一阶显式方程,(1.10)称为微分形式的一阶方程. n 阶隐式方程的一般形式为
n 阶显式方程的一般形式为
(1.12) (1.11)
在方程(1.11)中,如果左端函数F 对未知函数y 和它的各阶导数y ′,y ″,…,y (n ) 的全体而言是一次的,则称为线性常微分方程,否则称它为非线性常微分方程. 这样,一个以y 为未知函数,以x 为自变量的n 阶线性微分方程具有如下形式:
(1.13)
显然,方程(1.4)是一阶线性方程;方程(1.5)是一阶非线性方程;方程(1.6)是二阶线性方程;方程(1.7)是二阶非线性方程. 通解与特解
微分方程的解就是满足方程的函数,可定义如下. 定义1. 1 设函数
在区间I 上连续,且有直到n 阶的导数. 如果把
代入方程(1.11),得到在区间I 上关于x 的恒等式,
则称
为方程(1.11)在区间I 上的一个解.
这样,从定义1.1可以直接验证:
1. 函数y = x2+C是方程(1.4)在区间(-∞,+∞)上的解,其中C 是任意的常数. 2. 函数
是方程(1.5)在区间(-1,+1)上的解,其中C 是任意常
数. 又方程(1.5)有两个明显的常数解y =±1,这两个解不包含在上述解中. 3. 函数立的任意常数. 4. 函数
是方程(1. 7) 在区间(-∞,+∞)上的解,其中
和
是独立的
是方程(1.6)在区间(-∞,+∞)上的解,其中
和
是独
任意常数.
这里,我们仅验证3,其余留给读者完成. 事实上,在(-∞,+∞)上有
所以在(-∞,+∞)上有
从而该函数是方程(1.6)的解.
从上面的讨论中,可以看到一个重要事实,那就是微分方程的解中可以包含任意常数,其中任意常数的个数可以多到与方程的阶数相等,也可以不含任意常数. 我们把n 阶常微分方程(1.11)的含有n 个独立的任意常数C 1,C 2,…,C n 的
解
,称为该方程的通解,如果方程(1.11)
的解
不包含任
意常数,则称它为特解. 由隐式表出的通解称为通积分,而由隐式表出的特解称为特积分.
由上面的定义,不难看出,函
数
分别是方程(1.4),(1.5)和(1.6)的通解,函数
和是方程(1.7)
的通积分,而函数y =±1是方程(1.7)的特解. 通常方程的特解可对通解中的任意常数以定值确定,这种确定过程,需要下面介绍的初始值条件,或简称初值条件.
初值问题 例 1中的函数(1.3)显然是方程(1.2)的通解,由于
和
是两个任意常数,这表明
方程(1.2)有无数个解,解的图像见下面的图a 和图b 所示.
图a (C 1>固定,C 2>0) 图b (C 1=0,C2>0) 而实际经验表明,一个自由落体运动仅能有一条运动轨迹. 产生这种多解性的原因是因为方程(1.2)所表达的是任何一个自由落体,在任意瞬时t 所满足的关系式,并未考虑运动的初始状态,因此,通过积分求得的其通解(1.3)所描述的是任何一个自由落体的运动规律. 显然,在同一初始时刻,从不同的高度或以不同初速度自由下落的物体,应有不同的运动轨迹. 为了求解满足初值条件的解,我们可以把例1中给出的两个初始值条件,即
初始位置 x (0)= H 初始速度
代入到通解中,推得
于是,得到满足上述初值条件的特解为
(1.14)
它描述了初始高度为H ,初始速度为v 0的自由落体运动规律. 求微分方程满足初值条件的解的问题称为初值问题. 于是我们称(1.14)是初值问题
的解.
对于一个n 阶方程,初值条件的一般提法是
其中
是自变量的某个取定值,而
(1.15)
是相应的未知函数及导数的给
定值. 方程(1.12)
的初值问题常记为
初值问题也常称为柯西(Cauchy )问题. 对于一阶方程,若已求出通解
(1.16)
,只要把初值条件
代入通解中,得到方程
从中解出C ,设为
.
,代入通解,即得满足初值条件的解
对于n 阶方程,若已求出通解得到n 个方程式
后,代入初值条件(1.15),
如果能从(1.17)
式中确定出
.
例2 求方程
(1.17)
,代回通解,即得所求初值问题的
的满足初值条件 解 方程通解为
求导数后得
将初值条件代入,得到方程组
的解.
解出
和
得
故所求特解为
积分曲线
为了便于研究方程解的性质,我们常常考虑解的图象. 一阶方程(1.9)
的一个特解
的图象是xoy 平面上的一条曲线,称为方程(1.9)的积分曲线
,而通解的图象是平面上的一族曲线,称为积分曲线族. 例如,方程(1.4)的通解
+C是xoy 平面上的一族抛物曲线. 而
是过点(0,0) 的一条积分曲线. 以后,为了叙
述简便,我们对解和积分曲线这两个名词一般不加以区别. 对于二阶和二阶以上的方程,
也有积分曲线和积分曲线族的概念,只不过此时积分曲线所在的空间维数不同,我们将在第4章详细讨论.
最后,我们要指出,本书中按习惯用
而
分别代表,
本
分别代表
本节要点:
1.常微分程的定义,方程的阶,隐式方程,显式方程,线性方程,非线性方程. 2.常微分方程解的定义,通解,特解,通积分,特积分. 3.初值问题及初值问题解的求法. 4.解的几何意义,积分曲线.
作业:
练习1.1 1, 2.
1.指出下列方程的阶数,是否是线性方程:
(1) (2) (3)
(4) (5) (6)
2.验证给出函数是否为相应方程的解
(1) ,
,(C 为任意常数)
(2) , ,(C 为任意常数)
(3) (4)
,
,
答案:
1.(1)一阶,非线性 (2)一阶,非线性 (3)四阶,线性 (4)三阶,非线性 (5)二阶,非线性
2.(1)是 (2)是
(6)一阶,非线性 3)不是 (4)是
(