常微分方程 第二章.
常微分方程二阶线性微分方程习题课
二阶线性微分方程
例7 求微分方程 yy y2 y2 ln y 的通解.
解: y
0,
y
yy y2
y2
ln
y,
y
ln yx
y ,方程改写为: y
ln y
ln y,
令 z ln y z z 0,二阶常系数齐次线性方程
特征方程 2 1 0, 特征根 1.
通解 z C1e x C2e x ln y C1e x C2e x .
(3n 1)!
y( x)
x
x4
x7
x3n2
4! 7!
(3n 2)!
20
二阶线性微分方程
解 (2) 相应的齐次微分方程
y y y ex y(0) 1,y(0) 0
y y y 0, 特征方程 2 1 0
特征根
1,2
1 2
3 i, 2
非齐次方程的特解: y Ae x
将y, y, y 代入方程 A 1 ,
特征根的情况
通解的表达式
实根 1 2
实根 1 2
复根 1,2 i
y C1e1x C2e2 x y (C1 C2 x)e1x y e x (C1 cos x C2 sin x)
2
二阶线性微分方程
例2 求方程 y 4 y 4 y 0 的系数线性非齐次方程
f ( x) 2e x . 1.
(1) 求对应齐次方程的通解
特征方程 2 3 2 0,
特征根 1 1,2 2,
对应齐次方程通解 Y C1e x C2e2x
14
二线性微分方程
例10 设函数 y y( x)满足微分方程 y 3 y 2 y 2e x ,
例6 求 y(5) y(4) 2 y 2 y y y 0 的通解. 解 特征方程 5 4 2 3 2 2 1 0
《常微分方程》 (方道元 著) 课后习题答案 浙江大学出版社
= v0 + at.
dh dt |t=T
=0
2.一个湖泊的水量为V立方米,排入湖泊内含污染物A的污水量为V1 立方米/时,流入湖泊内不含污
0 不得超过 m 5 。试讨论湖泊中污染物A的浓度变化?
解:设污染物A的浓度为P(t),由题意可得 V P (t) + P (t)(V1 + V2 ) = P (0) = 5m
w
ω )e−s ds = y (x)。
4.考虑方程
w
.k
w
其中p(x)和q (x)都是以ω 为周期的连续函数,试证:
(1)若q (x) ≡ 0,则方程(2.4.23)的任一非零解以ω 为周期当且仅当函数p(x)的平均值 p ¯= 1 ω
ω
hd aw
答
dy + a(x)y ≤ 0, (x ≥ 0). dx
−
x 2y
= 0, y (0) = 1;
−2 ,令z = y 2 ,方程两边再乘以因子e−2x ,得到 (1)显然y ≡ 0是方程的解,当y = 0时,方程两边乘以 1 2y
方程的通解为 y = (Ce2x − x 1 2 − ) 4 8
hd aw
1 1
案 网
1.试求下列微分方程的通解或特解: √ dy − 4xy = x2 y ; (1) x dx
w
w
(3) y =
dy dx
1 1−x2 y = 1 + x, x ex + 0 y (t) dt; x4 +y 3 xy 2 ;
(4)
=
(5) 2xydy − (2y 2 − x)dx = 0;
(6) (y ln x − 2)ydx = xdy ;
常微分方程第二章
第二章 基本定理我们在第一章主要学习了初等积分法,掌握了几类常微分方程的解法.但是这些解法只适用于某些特殊的类型,很多其它的常微分方程不能用初等解法进行求解.1841年,法国数学家刘维尔(Liouville )证明了里卡蒂(Riccati )方程)0)(()()()(2≠++=x p x r y x q y x p dydx 除了某些特殊的类型外,一般不能用初等积分法求解.例如,很简单的里卡蒂方程22y x dxdy +=就不能用初等积分法求解.自然地,如果一个常微分方程不能用初等积分法求解,那么应该如何处理呢?是否存在解呢?如果存在解,它的解是否唯一呢?解的存在区间是什么呢?初值的微小误差对解有什么影响呢?这些问题在理论的研究和实际应用中,都有着重要的意义.本章将解决这些基本问题. 本章主要介绍解的存在唯一性定理、解的延展定理与比较定理、解对初值的连续依赖性定理以及解对初值的可微性定理,这些定理就回答了我们刚才的疑问,有效的处理解的存在性、唯一性、存在区间、初值对解的影响等问题,为我们使近似解法奠定理论基础,同时这些定理也是常微分方程理论的基础内容,对进一步的学习奠定基础.2.1 解的存在唯一性定理对于一般的常微分方程),(y x f dxdy = (2.1) 如果给出了初始条件00)(y x y =,我们就得到了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy (2.2) 这时,在什么样的条件下,柯西初值问题的解存在且唯一呢?解的存在区间是什么呢?我们有如下的解的存在唯一性定理.2.1.1 存在唯一性定理的叙述定理2.1(存在唯一性定理)如果方程(2.1)的右端函数),(y x f 在闭矩形区域b y y b y a x x a x R +≤≤-+≤≤-00002,:上满足如下条件:(1)在2R 上连续;(2)在2R 上关于变量y 满足李普希兹(Lipschitz )条件,即存在常数N ,使对于2R 上的任何一对点),(y x 和),(x 有不等式:y y N y x f y x f -≤-),(),(则初值问题(2.2)在区间],[0000h x h x +-上存在唯一解00)(),(y x x y ==ϕϕ 其中),(max ),,min(),(0y x f M Mb a h R y x ∈==. 在给出定理2.1的证明之前,我们先对定理2.1的条件和结论做些说明:1、在两个条件中,条件(2),即李普希兹条件比较难于验证,因为李普希兹常数N 难以确定.但是,我们可以将该条件加强,替换为:如果函数),(y x f 在闭矩形区域2R 关于y 的偏导数),(y x f y '存在且有界.这样,可以推出李普希兹条件成立.事实上,因为),(y x f y '有界,故设N y x f y ≤'),(,对2),(),,(R x y x ∈∀,由拉格朗日中值定理得:y y N y y x f y x f y x f y -≤-'=-),(),(),(ξ我们验证),(y x f y '在闭矩形区域2R 上有界也不容易,可以进一步将条件加强为:),(y x f y '在闭矩形区域2R 上连续.由闭区域上连续函数的性质知:),(y x f y '在闭矩形区域2R 上有界,所以李普希兹条件成立.因此,有如下的关系式:),(y x f y '在2R 上连续⇒),(y x f y '在2R 上存在且有界⇒李普希兹条件2、在定理2.1的结论中,解)(x y ϕ=的存在区间为],[0000h x h x +-,其中 ),(max ),,min(),(0y x f M Mb a h R y x ∈==.为什么解的存在区间不是],[00a x a x +-呢?这是因为我们研究问题的范围为闭矩形区域2R ,方程的解)(x y ϕ=不能超出2R 的范围,又因为),(max ),(y x f M Ry x ∈=,所以M y x f M ≤≤-),( 即 M dxdy M ≤≤- 由⎪⎩⎪⎨⎧=-=00)(y x y M dx dy 和⎪⎩⎪⎨⎧==00)(y x y M dx dy 得:001)()(y x x M x y +--=,002)()(y x x M x y +-= 因此)()()(21x y x y x y ≤=≤ϕ,即)(x y ϕ=夹在)(1x y 与)(2x y 之间.又,)(1x y 与)(2x y 在2R 上的存在区间为],[0000h x h x +-,故)(x y ϕ=的存在区间也是],[0000h x h x +-.2.1.2 存在性的证明首先,我们给出柯西初值问题(2.2)的等价转化,即求(2.2)的解)(x y ϕ=,等价于求解积分方程⎰+=xx d y f y y 0))(,(0ξξξ (2.3) 事实上,如果)(x y ϕ=是初值问题(2.2)的解,即有))(,()(x x f x ϕϕ='且00)(y x =ϕ从0x 到x 积分得:⎰+=xx d f y x 0))(,()(0ξξϕξϕ 即)(x y ϕ=是积分问题(2.3)的解.反过来,如果)(x y ϕ=是积分问题(2.3)的解,即有⎰+=xx d f y x 0))(,()(0ξξϕξϕ 则00)(y x =ϕ且))(,()(x x f x ϕϕ='即)(x y ϕ=是初值问题(2.2)的解.经过等价转化,我们将初值问题(2.2)的求解,转化为积分问题(2.3)的求解.下面用皮卡(Picard )逐次逼近来证明积分问题(2.3)的解的存在性,分为三个步骤:1、构造近似函数列{})(x n ϕ任取一个满足初值条件00)(y x y =的函数)(0x y ϕ=作为首项(初始项),并要求在2R 上的存在区间为:],[0000h x h x +-,简单起见,取00)(y x =ϕ,将它代入方程(2.3)的右端,所得到的函数用)(1x ϕ表示,并称为一次近似,即⎰+=xx d f y x 0))(,()(001ξξϕξϕ 再将)(1x ϕ代入方程(2.3)的右端就得到二次近似⎰+=xx d f y x 0))(,()(102ξξϕξϕ 序行此法,可以得到n 次近似⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ 为了保证上述的逐次逼近过程可以一直进行下去,必须有2))(,(R x x n ∈ϕ,即当],[0000h x h x x +-∈时,有,2,1)(0=≤-n b y x n ϕ 下面用数学归纳法证明b y x n ≤-0)(ϕ.显然,当],[0000h x h x x +-∈时,有b y y y x ≤=-=-0)(0000ϕ假设,当],[0000h x h x x +-∈时,有b y x n ≤--01)(ϕ,那么,对于)(x n ϕ有⎰-=-xx n n d f y x 0))(,()(10ξξϕξϕ 从而有b Mb M Mh x x M d f y x xx n n =≤≤-≤≤-⎰-00100))(,()(ξξϕξϕ 由数学归纳法知,当],[0000h x h x x +-∈时,有,2,1)(0=≤-n b y x n ϕ这样,我们就可以得到一个近似函数列{})(x n ϕ.2、证明近似函数列{})(x n ϕ在区间],[0000h x h x +-上一致收敛.由于无法得到{})(x n ϕ的通项公式,只知道首项和递推关系式,直接证明函数列{})(x n ϕ的收敛性比较困难,为此我们构造函数项级数+-++-+-)]()([)]()([)(1010x x x x x n n ϕϕϕϕϕ (2.4) 它的部分和是)()]()([)]()([)()(10101x x x x x x x S n n n n ϕϕϕϕϕϕ=-++-+=-+因此,证明{})(x n ϕ的收敛性转化为证明级数(2.4)的收敛性,下面我们证明级数(2.4)在区间],[0000h x h x +-上一致收敛.首先研究级数(2.4)的通项)(x n μ⎰=-xx d f x x 0))(,()()(001ξξϕξϕϕ 即⎰=-xx d y f y x 0),()(001ξξϕ 所以00010),()(x x M d y f y x x x -≤≤-⎰ξξϕ 因为⎰+=x x d f y x 0))(,()(001ξξϕξϕ,⎰+=x x d f y x 0))(,()(102ξξϕξϕ,所以 ⎰-≤-x x d f f x x 0))(,())(,()()(0112ξξϕξξϕξϕϕ由李普希兹条件,得 !2)()()()(200011200x x MN d x MN d N x x x x x x -=-≤-≤-⎰⎰ξξξξϕξϕϕϕ 下面用数学归纳法证明!)()(011n x x MN x x nn n n -≤---ϕϕ 显然,2,1=n 的时候,不等式成立(上面已经给出), 假设!)()(011n x x MN x x n n n n -≤---ϕϕ成立,那么对于1+n 的情形有 )!1(!)()())(,())(,()()(100111000+-=-≤-≤-≤-+--+⎰⎰⎰n x x MN d n x MN d N d f f x x n n x x n n xx n n x x n n n n ξξξξϕξϕξξϕξξϕξϕϕ由数学归纳法知,对一切自然数n ,均有!)()(011n x x MNx x nn n n -≤---ϕϕ 又00h x x ≤-,所以级数(2.4)的通项满足: !)(011n h MN v x n n n n -+=≤μ ( ,2,1=n ) 利用比式判别法,可知以n v 为通项的级数收敛,从而以)(x n μ为通项的级数(2.4)绝对收敛且一致收敛.又,每一个)(x n μ是连续的,所以级数(2.4)的和函数也是连续的,记为)(x ϕ,其存在区间也是],[0000h x h x +-.因此函数列{})(x n ϕ就收敛于)(x ϕ.3、证明)(lim )(x x n n ϕϕ∞→=是积分问题(2.3)的解,从而也是初值问题(2.2)的解.在⎰-+=x x n n d f y x 0))(,()(10ξξϕξϕ两端取极限,得到 ⎰-∞→∞→+=xx n n n n d f y x 0))(,(lim )(lim 10ξξϕξϕ 即⎰+=xx d f y x 0))(,()(0ξξϕξϕ 所以)(x ϕ是积分问题(2.3)的解,从而也是初值问题(2.2)的解.2.1.3 唯一性的证明下面我们证明解的唯一性.在证明唯一性之前,先介绍一个重要的不等式,即贝尔曼(Bellman )不等式.贝尔曼引理 设)(x y 为区间],[b a 上的非负连续函数,b x a ≤≤0.若存在,0≥δ 0≥k ,使得)(x y 满足不等式],[,)()(0b a x d y k x y xx ∈+≤⎰ττδ (2.5) 则有],[,)(0b a x e x y x x k ∈≤-δ证明 仅证明0x x ≥的情形,0x x ≤的情形类似.令)(x y 的原函数为⎰=xx d y x R 0)()(ττ,代入(2.5)得 δ≤-')()(x kR x R两边同时乘以积分因子)(0x x k e --,得)()(00)]()([x x k x x k e x kR x R e ----≤-'δ从0x 到x 积分得)()(00)(x x k x x k e e x kR -----≤δδ即)(0)(x x k e x kR -≤+δδ 由(2.5)知,)()(x kR x y +≤δ,所以],[,)(0b a x e x y x x k ∈≤-δ下面证明积分问题(2.3)的解的唯一性.假设积分问题(2.3)有两个解)(1x y 和)(2x y ,我们只需要证明:)(1x y )(2x y ≡,],[0000h x h x x +-∈事实上,因为⎰+=x x d y f y x y 0))(,()(101ξξξ,⎰+=xx d y f y x y 0))(,()(202ξξξ 所以有⎰-≤-xx d y f y f x y x y 0))(,())(,()()(2121ξξξξξ由李普希兹条件知⎰-≤-xx d y y N x y x y 0)()()()(2121ξξξ 令N k x y x y x y ==-=,0,)()()(21δ,由贝尔曼引理可知,0)(=x y ,即)(1x y )(2x y ≡. 这样,我们就完成了解的存在性与唯一性的证明.2.1.4 三点说明为了更好的理解和掌握解的存在唯一性定理,我们对该定理再做三点说明.1、在存在性的证明过程中,我们利用逐次逼近法构造了近似函数列{})(x n ϕ,其中首项为:00)(y x =ϕ,递推关系式为:⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ.该方法实际上给出了我们一种求初值问题(2.2)的近似解的方法,当用n 次近似解逼近精确解时,需要给出它的误差估计.事实上,有∑∑∞+=∞=+-≤-≤-101!)()()()(n k k k nk k k n k x x N N M x x x x ϕϕϕϕ 0)!1()(!)!1()(!10001010Nh n k k k n n k k k e n Nh N M k h N n Nh N M k h N N M +=+<≤+∞=+∞+=∑∑ 2、如果方程(2.1)是线性方程,即)()(x q y x p dxdy +-= 其中)(x p 和)(x q 在区间],[b a 上连续,这时,初值问题(2.2)在带型区域+∞<<-∞≤≤y b x a R ,:2满足定理2.1的条件.事实上,)()(),(x q y x p y x f +-=在2R 上连续,而且)(),(x p y x f y -='在2R 上也连续,所以),(y x f 关于变量y 满足李普希兹条件.这时,初值问题(2.2)的解存在且唯一,存在区间为],[b a .3、定理2.1中的李普希兹条件是保证解唯一的充分条件,那么这个条件是不是必要条件呢?回答是否定的,即李普希兹条件是解唯一的充分非必要条件.下面我们给出一个例子来说明李普希兹条件是解唯一的非必要条件,也就是说,即使李普希兹条件不成立,初值问题(2.2)的解也可能是唯一的.例1 试证方程0,ln ,0≠=⎩⎨⎧=y y y y dx dy 经过xOy 平面上任一点的解都是唯一的.证明 由00,ln ,0≠=⎩⎨⎧=y y y y dx dy 可得:0=y 或x Ce e y ±=. 任给xOy 平面上的一个点),(00y x ,只会对应0=y 或xCe e y ±=中的一个解,也就是说,过xOy 平面上任一点的解都是唯一的.但是,我们有0ln ln )0,(),(-==-y y y y x f y x f 因为+∞=→y y ln lim 0,所以找不到0>N ,使得 0)0,(),(-≤-y N x f y x f从而方程右端函数在0=y 的任何邻域上不满足李普希兹条件,但是初值问题(2.2)的解却是唯一的,这说明李普希兹条件是非必要条件.习 题 2.11.试判断方程y x dx dy tan =在区域 (1)π≤≤≤≤-y x R 0,11:1;(2)44,11:2ππ≤≤-≤≤-y x R上是否满足定理2.1的条件?2.讨论方程3123y dx dy =在怎样的区域中满足定理2.1的条件.并求通过)0,0(的一切解.3.试用逐次逼近法求方程2y x dxdy -=满足初值条件0)0(=y 的近似解: )(),(),(),(3210x x x x ϕϕϕϕ并在闭矩形区域11,11:2≤≤-≤≤-y x R 给出三次近似的误差估计.4.利用逐次逼近法求方程22x y dxdy -=适合初值条件1)0(=y 的近似解: )(),(),(210x x x ϕϕϕ并在闭矩形区域111,11:2≤-≤-≤≤-y x R 给出二次近似的误差估计.5.试证明定理2.1中的n 次近似解)(x n ϕ与精确解)(x ϕ有如下的误差估计式:10)!1()()(+-+≤-n n n x x n MN x x ϕϕ 6.在条形区域+∞<≤≤y b x a ,内,假设方程(2.1)的所有解都唯一,对其中任意两个解)(),(21x y x y ,如果有)()(0201x y x y <,则必有b x x x y x y ≤≤<021),()(.7.讨论方程323y dx dy = 解的唯一性.2.2 延展定理和比较定理由解的存在唯一性定理,我们知道,初值问题(2.2)的解在满足一定条件的情况下存在且唯一,但是解的存在区间不是],[00a x a x +-,而是],[0000h x h x +- 其中),(max ),,min(),(0y x f M Mb a h R y x ∈==.如果M 比较大的话,则解的存在区间就非常小,这对我们研究解的性质产生了很大的局限性,只能在很小的范围内有解,当x 超出这个范围时,解的情况就不清楚了.为了解决这个问题,我们有下面的延展定理.2.2.1 延展定理定理2.2(延展定理)如果方程(2.1)的右端函数在区域R R D ⨯⊂上连续,且关于变量y 满足局部的李普希兹条件,即对于D 内的任一闭矩形区域都满足李普希兹条件,则对任何一点D y x ∈),(00,初值问题(2.2)的解)(x y ϕ=可以向左右无限延展,直到))(,(x x ϕ任意接近区域D 的边界.在给出定理的证明之前,先对“))(,(x x ϕ任意接近区域D 的边界”进行说明.当区域D 有界时,积分曲线向左右延展可以任意接近;当区域D 无界时,积分曲线向左、右延展,或者任意接近区域D 的边界(边界存在的话),或者无限远离坐标原点.证明 首先证明区域D 有界的情形.设区域D 的边界为D D L -=(D 为D 的闭包).对于任意给定的正数ε,记L 的ε邻域为εU ,记L 的2ε邻域为2εU ,记L 的4ε邻域为4εU .则集合22εεU D D -=为闭集,且D D ⊂2ε,所以2εD 有界. 只要证明积分曲线可以到达2εD 的边界2εL ,由ε的任意性知,积分曲线就可以任意接近区域D 的边界L .事实上,以2εD 中的任意一点为中心,以4ε为半径的闭圆区域均包含在区域D 的内部.且在闭区域44εεU D D -=之内.从而,以2εD 中的任意一点为中心,以4221ε=a 为边长的正方形也在闭区域4εD 之内.记 ),(max 4),(1y x f M D y x ε∈= 则过2εD 的任意一点),(**y x 的积分曲线,必至少可在区间],[**h x h x +-上存在,其中)82,82min(),min(1111M M a a h εε==. 于是,过点),(00y x 的积分曲线)(x y ϕ=每向左或向右延展一次,其存在区间就伸长一个确定的正数h ,由于2εD 有界,)(x y ϕ=经过有限次延展后一定可以达到2εD的边界2εL .于是也就可以任意接近区域D 的边界L .其次考虑区域D 为无界的情形.这时,我们可以用闭圆区域,2,1},),{(222=≤+=n n y x y x S n与区域D 取交集,令n n S D D =,则 ∞==1n n D D .由于n D 为有界的区域,根据前面的证明,我们可知,过n D 内任一点的积分曲线能够任意接近n D 的边界.因此,过点),(00y x 的积分曲线)(x y ϕ=可以无限接近区域D 的边界.延展定理的证明,关键是第一步证明,也就是区域D 有界的时候,过点),(00y x 的积分曲线)(x y ϕ=向左向右延展的时候,一定要做等速延展,即延展步幅h 是不变的. 例1 试讨论方程2y dxdy=通过点)1,1(的解和通过点)1,3(-的解的存在区间. 解 该题目中研究问题的区域D 为整个坐标平面xOy .方程右端函数满足延展定理的条件.由2y dxdy=可以解得方程的通解为 xC y -=1代入1)1(=y 得:2=C .故通过点)1,1(的解为xy -=21 它可以向左无限延展,而当-→2x 时,+∞→y ,所以通过点)1,1(的解xy -=21的存在区间为)2,(-∞.代入1)3(-=y 得:2=C .故通过点)1,3(-的解为xy -=21它可以向右无限延展,而当+→2x 时,-∞→y ,所以通过点)1,3(-的解xy -=21的存在区间为),2(+∞.这个例子说明,尽管),(y x f 在整个坐标平面上满足延展定理的条件,解上的点))(,(x x ϕ也能无限接近区域D 的边界,但是延展的方向却不一定是无限向右和向左,可能是向上或向下,从而导致解的存在区间不是),(+∞-∞. 例2 试证明:对任意的0x 及满足条件100<<y 的0y ,方程221)1(y x y y dx dy ++-=的满足条件00)(y x y =的解)(x y y =在),(+∞-∞上存在.证明:令221)1(),(y x y y y x f ++-=,则222222)1(122),(y x x y y x y y x f y ++--++=' 显然),(),,(y x f y x f y '在xOy 平面上连续,满足解的存在唯一性条件及延展定理的条件,而1,0==y y 是),(y x f dxdy=的解, 因此,满足00)(y x y =,100<<y 的解存在,而且可以无限延展到xOy 平面的边界,且不能穿过1,0==y y ,故只能向左右无限延展,所以,)(x y y =在),(+∞-∞上存在.该例题说明,),(y x f 在整个坐标平面上满足延展定理的条件,当方程的解不能穿过1,0==y y 时,它就不能向上向下无限延展了,只能向左、向右延展,所以解的存在区间就是),(+∞-∞.在这里,1,0==y y 控制了解的延展方向,使它按照我们的要求进行延展,因此就有了下面的比较定理. 2.2.2 比较定理我们在使用延展定理的时候,通常会和比较定理配合使用,从而起到控制延展方向的作用.下面介绍一下比较定理.我们在考察方程(2.1)),(y x f dxdy=时,通常将右端函数),(y x f 进行放缩的处理,比如),(),(),(21y x F y x f y x F <<这时,我们可以同时考察),(1y x F dx dy =和),(2y x F dxdy = 我们有如下的比较定理:定理2.3 (第一比较定理)设定义在某个区域D 上的函数),(y x f ,),(1y x F 和),(2y x F 满足条件:(1)在D 满足解的存在唯一性定理及延展定理的条件,即在D 上连续,在D 上关于变量y 满足李普希兹条件;(2)在D 上有不等式),(),(),(21y x F y x f y x F <<设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ,⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=,)(1x y Φ=和)(2x y Φ=,则在它们的共同存在区间上有下列不等式:021),()()(x x x x x >Φ<<Φϕ 021),()()(x x x x x <Φ>>Φϕ证明 仅证当0x x >时,)()(2x x Φ<ϕ,其它的情形相类似. 由比较定理的条件(1),初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解在0x 的某一邻域内存在且唯一,分别记为)(x y ϕ=和)(2x y Φ=,它们满足0020)()(y x x =Φ=ϕ令)()()(2x x x h ϕ-Φ=,则0)()()(0020=-Φ=x x x h ϕ且0))(,())(,()()()(0002020020>-Φ='-Φ'='x x f x x F x x x h ϕϕ所以函数)(x h 在0x 的某一右邻域内是严格单调增加的.如果在0x x >时,0)(>x h 不是总成立,则至少存在一点01x x >,使得0)(1=x h ,且当10x x x <<时,0)(>x h ,因此在点1x 的左导数0)0(1≤-'x h ,这与0))(,())(,()()()(1112121121>-Φ='-Φ'='x x f x x F x x x h ϕϕ矛盾.因此当0x x >时,0)(>x h 总成立,即)()(2x x Φ<ϕ.比较定理的应用,关键是),(1y x F 和),(2y x F 的选取,因为初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解)(x y ϕ=的存在区间的延展,受到)(1x y Φ=和)(2x y Φ=的控制,即)(x y ϕ=夹在)(1x y Φ=和)(2x y Φ=之间.因此,我们必须能确定出)(1x y Φ=和)(2x y Φ=的存在区间,这就是我们选取),(1y x F 和),(2y x F 的标准,即⎪⎩⎪⎨⎧==001)(),(y x y y x F dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解)(1x y Φ=和)(2x y Φ=必须能够求得. 下面我们给出第二比较定理.定理2.4 (第二比较定理)设定义在某个区域D 上的函数),(y x f ,),(1y x F 和),(2y x F 满足条件:(1)在D 满足解的存在唯一性定理及延展定理的条件,即在D 上连续,在D 上关于变量y 满足李普希兹条件;(2)在D 上有不等式),(),(),(21y x F y x f y x F ≤≤设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ,⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=,)(1x y Φ=和)(2x y Φ=,则在它们的共同存在区间上有下列不等式:021),()()(x x x x x >Φ≤≤Φϕ 021),()()(x x x x x <Φ≥≥Φϕ习 题 2.21.设方程为),()(22y x f a y dxdy-= 假设),(y x f 及),(y x f y '在xOy 平面上连续,试证明:对于任意的0x 及a y <0,方程满足00)(y x y =的解都在),(+∞-∞上存在.2.指出方程2)1(2xy e y dxdy -=的每一个解的最大存在区间,以及当x 趋于这个区间的右端点时解的极限.3.讨论方程xx dx dy 1cos 12-= 解的存在区间.4.设),(y x f 在整个平面上连续有界,对y 有连续偏导数,试证明方程),(y x f dxdy=的任一解)(x y ϕ=在区间+∞<<∞-x 上有定义. 5.讨论方程212-=y dx dy 的通过点)0,0(的解,以及通过点)3,2(ln -的解的存在区间.6.在方程)(y f dxdy=中,如果)(y f 在),(+∞-∞上连续可微,且 )0(0)(≠<y y yf ,求证方程满足00)(y x y =的解)(x y 在区间),[0+∞x 上存在,且有0)(lim =+∞→x y x .2.3 解对初值的连续依赖性定理和解对初值的可微性定理通过前两节的存在唯一性定理和延展定理,加上比较定理,我们知道了初值问题(2.2)在什么样的条件下,解是存在的,是唯一的,而且存在区间比较小的时候,通过延展定理和比较定理可以将解的存在区间变大,从而在实际问题中可以达到我们的要求.但是,在实际问题中,还有一个问题需要解决,那就是误差问题.我们的初始条件00)(y x y =如果产生了微小的偏差,这个偏差对我们的初值问题(2.2)的解)(x y ϕ=会有什么影响呢?下面我们来解决这个问题. 我们在研究初值问题(2.2)的时候,习惯上把0x 和0y 当作常数来看待,这样初值问题(2.2)的解)(x y ϕ=被看作x 的函数.实际上,如果0x ,0y 变化,初值问题(2.2)的解)(x y ϕ=也会发生变化.例如方程xydx dy = 经过点),(00y x 的解为x x y y 0=,可以看作00,,y x x 的函数.对于一般的情形,初值问题(2.2)的解也可以看作00,,y x x 的函数,记为),,(00y x x y ϕ=,代入00)(y x y = 得:0000),,(y y x x =ϕ.如果我们的初始条件00)(y x y =发生了微小的误差,变为了**0)(y x y =,初值问题(2.2)的解也变化不大的话,称解连续依赖于初值.下面我们给出连续依赖性的严格定义.定义2.1 设初值问题⎪⎩⎪⎨⎧==**0)(),(y x y y x f dxdy的解),,(*0*0y x x y ϕ=在区间],[b a 上存在,如果对于任意给定的正数ε,存在正数δ (δ的选取与,ε**0,y x 有关),使得对于满足δδ<-<-*00*00,y y x x (2.2)的解),,(00y x x y ϕ=都在],[b a 上存在,且有],,[,),,(),,(*0*000b a x y x x y x x ∈<-εϕϕ则称初值问题(2.2)的解),,(00y x x y ϕ=在点),(*0*0y x 连续依赖于初值,0x 0y .定理2.4 (解对初值的连续依赖性定理)设),(y x f 在区域D 内连续,且关于变量y 满足李普希兹条件.如果D y x ∈),(*0*0,初值问题(2.2)有解),,(*0*0y x x y ϕ=,且当b x a ≤≤时,D y x x x ∈)),,(,(*0*0ϕ,则对任意的正数ε,存在0>δ,使对于满δδ<-<-*00*00,y y x x的任意),(00y x ,初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解),,(00y x x y ϕ=也在区间],[b a 上存在,且有εϕϕ<-),,(),,(*0*000y x x y x x证明 对于任意给定的正数ε,取εδ<<10,使得闭区域}),,(,),{(1*0*0δϕ≤-≤≤=y x x y b x a y x U整个含在区域D 内,这是可以做到的,因为区域D 是开区域,且当b x a ≤≤时,D y x x x ∈)),,(,(*0*0ϕ,所以,只要1δ的选取足够小,以曲线),,(*0*0y x x y ϕ=为中线,宽度为12δ的带形开区域U 就整个包含在区域D 内, 选取δ满足)(110a b N e M--+<<δδ其中N 为李普希兹常数,),(max ),(y x f M Uy x ∈=,同时还要求δ的选取,必须保证闭正方形δδ≤-≤-*0*02,:y y x x R含于带形开区域U 内.由存在唯一性定理知,对于任一200),(R y x ∈,初值问题(2.2)在0x 的某邻域上存在唯一解),,(00y x x y ϕ=,而且),,(00y x x y ϕ=在0x 的该邻域上可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(000000⎰+=而),,(*0*0y x x y ϕ=可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(*0*0*0*0*0*⎰+=对上述两式做差得:ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx x x )),,(,()),,(,(),,(),,(*0*000*00*0*000*⎰⎰-+-=-ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx xx )),,(,()),,(,(),,(),,(*0*000*00*0*000*0⎰⎰-+-≤-ττϕτττϕττϕτd y x f d y x f y x f y y x x xx |)),,(,(||)),,(,()),,(,(|0000*0*0*00**0⎰⎰+-+-≤δττϕττϕτδM d y x f y x f xx +-+≤⎰|)),,(,()),,(,(|00*0*0*0ττϕτϕδd y x y x N M xx |),,(),,(|)1(00*0*0*0-++≤⎰由贝尔曼引理,得εδδδϕϕ<<+≤+≤---1)(*0*000)1()1(),,(),,(*a b N x x N e M e M y x x y x x因此,只要在),,(00y x x y ϕ=有定义的区间上,就有εϕϕ<-),,(),,(*0*000y x x y x x .下面我们证明:),,(00y x x y ϕ=在区间],[b a 上有定义.事实上,因为εϕϕ<-),,(),,(*0*000y x x y x x即解),,(00y x x y ϕ=夹在εϕ+=),,(*0*0y x x y 和εϕ-=),,(*0*0y x x y 之间,而且,初值问题(2.2)满足延展定理的条件,所以,解),,(00y x x y ϕ=可以向左向右无限延展,直到无限接近区域D 的边界,于是,它在延展的时候,必须由直线a x =和直线b x =穿出区域U ,从而),,(00y x x y ϕ=在区间],[b a 上有定义.解对初值的连续依赖性说明,初值),(00y x 无法准确得到,但是我们能得到测量数据),(*0*0y x ,只要误差比较小,即δδ<-<-*00*00,y y x x .我们就可以用),(*0*0y x 代替),(00y x 去计算,得到初值问题的解),,(*0*0y x x y ϕ=,这个解可以非常接近真实解),,(00y x x y ϕ=,即εϕϕ<-),,(),,(*0*000y x x y x x .同理,如果方程的右端函数),(y x f 不能准确得到,只能得到),(y x f 的近似函数),(~y x f ,即)),((,),(),(~D y x y x f y x f ∈<-δ我们就可以用),(~y x f 代替),(y x f 去计算,得到初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 的解),,(00~y x x y ϕ=,那么),,(00~y x x y ϕ=能否代替),,(00y x x y ϕ=呢?我们有下面的解的连续依赖性定理.定理2.5 (解对被积函数的连续依赖性定理)在区域D 上,),(y x f 和),(~y x f 都连续,而且关于变量y 满足李普希兹条件, 若初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 在b x a ≤≤上有解),,(00~y x x y ϕ=,则对任意给定的正数ε,存在0>δ,只要),(y x f 满足)),((,),(),(~D y x y x f y x f ∈<-δ则初值问题(2.2)的解),,(00y x x y ϕ=在b x a ≤≤上存在,且有εϕϕ<-),,(),,(00~00y x x y x x .证明 由解的存在唯一性定理知,初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dxdy 的解),,(00~y x x y ϕ=存在,设其存在区间为],[b a ,且有⎰+=xx d y x f y y x x 0))],,(,([),,(00~~000~ξξϕξϕ而初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解),,(00y x x y ϕ=也存在,且可以表示为⎰+=xx d y x f y y x x 0))],,(,([),,(00000ξξϕξϕ则⎰⎰-=-xx xx d y x f d y x f y x x y x x 0))],,(,([))],,(,([),,(),,(0000~~0000~ξξϕξξξϕξϕϕ从而有⎰-≤-xx d y x f y x f y x x y x x 0|)),,(,()),,(,(|),,(),,(0000~~0000~ξξϕξξϕξϕϕ⎰-+-=xx d y x f y x f y x f y x f 0|)),,(,()),,(,()),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ ⎰-+-≤xx d y x f y x f y x f y x f 0|)),,(,()),,(,(||)),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ⎰+-≤xx d y x y x N 0)|),,(),,((|0000~ξδξϕξϕ ⎰-+-≤xx d y x y x N a b 0|),,(),,(|)(0000~ξξϕξϕδ由贝尔曼引理,得)(0000~)(),,(),,(a b N e a b y x x y x x --≤-δϕϕ取)(a b N e ab ---<εδ,则εϕϕ<-),,(),,(0000~y x x y x x .且解),,(00y x x y ϕ=在b x a ≤≤上存在. 例1 考虑方程,ln ,0≠=⎩⎨⎧-=y y y y dx dy 解的情况.解 显然1,1,0-===y y y 是方程的解,当1,1,0-≠≠≠y y y 时,有y y dxdyln -= 这时解得上半平面的通解为x Ce e y -=,下半平面的通解为xCe e y --=.可以看到,对于Ox 轴上的初值)0,(0x ,在任意有限闭区间上解对初值连续依赖,但是,在),0[+∞上,无论),(00y x ,00≠y 如何接近)0,(0x ,只要x 充分大,过),(00y x 的积分曲线就不能与过)0,(0x 的积分曲线(即0=y )任意接近了.这个例子说明,解在有限闭区间上对初值连续依赖,不能推广到无限区间,即,在无限区间上解对初值的连续依赖定理就不成立了.我们有时不仅要求解对初值连续依赖,而且还要知道解),,(00y x x y ϕ=对初值00,y x 的偏导数00,y x ∂∂∂∂ϕϕ是否存在.下面给出解对初值的可微性定理. 定理2.6 (解对初值的可微性定理)如果函数),(y x f 以及),(y x f y '在区域D 内连续,则初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解),,(00y x x y ϕ=在它有定义的区间上有连续偏导数00,y x ∂∂∂∂ϕϕ.并且有 ⎰-=∂∂'x x y d y x f e y x f x y x x 000)),,(,(00000),(),,(ττϕτϕ 及⎰=∂∂'xx y d y x f e y y x x 000)),,(,(000),,(ττϕτϕ 习 题 2.31.若函数),(y x f ,),(y x R 在区域D 内连续且满足李普希兹条件,设初值问题⎪⎩⎪⎨⎧=+=*0*0)(),(),(y x y y x R y x f dx dy 的解为),,(*0*0~y x x y ϕ=,存在区间为],[b a .对任意的正数ε,存在0>δ,使对于满足)),((,),(D y x y x R ∈<δ的),(y x R ,以及满足δδ<-<-*00*00,y y x x的任意),(00y x ,初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解),,(00y x x y ϕ=也在区间],[b a 上存在,且有εϕϕ<-),,(),,(*0*0~00y x x y x x 2.已知方程)sin(xy dxdy = 试求0000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x x y x x y 和0000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x y y x x y 3.设),,(00y x x ϕ是初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的解,试证明0),(),,(),,(00000000=∂∂+∂∂y x f y y x x x y x x ϕϕ 2.4 欧拉折线法在第一章,我们介绍了方程的初等解法,即用微积分的知识求得常微分方程的函数解.但是绝大多数的方程不能用初等方法求解,在第二章的前三节中,我们给出了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 在什么样的条件下,解存在且唯一;在什么条件下,解的存在区间可以延展;在什么条件下连续依赖于初值;在什么条件下,解对初值是可微的.有了这些准备,我们就可以研究柯西初值问题的近似解.下面我们介绍求近似解的方法,欧拉折线法.假定函数),(y x f 在区域:+∞<<-∞≤≤y b x a ,上连续,且关于变量y 满足李普希兹条件,求柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 在区间],[0b x 上的近似解,我们采用的方法是:(1)等分区间],[0b x ,分点为n k kh x x k ,,1,0,0 =+=;小区间长度nx b h 0-=, (2)第一个小区间上用切线段逼近曲线:))(,(0000x x y x f y y -+=,(3)求出1x 所对应的纵坐标))(,(010001x x y x f y y -+=,(4)依次重复(2),(3)得到每个小区间上的线段,从而得到欧拉折线. 这样,我们就用欧拉折线作为柯西初值问题在区间],[0b x 近似解.欧拉折线法的前提是:柯西初值问题的解存在且唯一,而且解的存在区间是],[0b x .例1试用欧拉折线法,取步长1.0=h ,求初值问题⎪⎩⎪⎨⎧=+=1)1(22y y x dx dy 的解在4.1=x 时的近似值.解 令22),(y x y x f +=,2)1,1(=f ,这时12-=x y ,代入1.11=x 得:2.11=y ,65.2)2.1,1.1(=f ,这时2.1)1.1(65.2+-=x y , 代入2.12=x 得:465.12=y ,586225.3)465.1,2.1(=f ,这时465.1)2.1(586225.3+-=x y , 代入3.13=x 得:8236225.13=y ,0155990225.5)8236225.1,3.1(=f ,这时8236225.1)3.1(0155990225.5+-=x y ,代入4.14=x 得:53251824022.24=y 习 题 2.41. 试用欧拉折线法,取步长1.0=h ,求初值问题⎪⎩⎪⎨⎧=-=0)0(22y y x dx dy 的解在5.1=x 时的近似值.2.试用欧拉折线法,取步长1.0=h ,求初值问题 ⎪⎩⎪⎨⎧=+=2)1(22y y x dx dy 在区间]4.1,1[上的近似解.。
常微分方程第二章第一讲
2.1.2 可化为变量分离方程的类型
引言 有的微分方程从表面上看,不是可分 离变量的微分方程,但是,通过适当的变量替 换,就可以很容易地化为“变量分离方程”, 在这里,介绍两类这样的方程。 1. 齐次方程
1)方程的类型
定义
dy y g ( ) (2.5) 的方程,称为齐次 dx x 微分方程,这里 g (u ) 是 u 的连续函数。 14
dy ( y) f ( x)dx C (2.2)
可以证明这就是方程(2.1)的通解.
2)如果存在 y0, ( y0 ) 0, 则方程( .1 使 2 )还有特解
y y0
(**)
微分方程(2.1)的所有解为:式(2.2)和(**).
注意:积分常数C 的相对任意性。
7
3.变量分离方程的解题步骤
即 1 , 2 1 ,
则 ON OM ,
PM 而 tan 2 , OP ON
_____ _____
则有 y'
y x x y
2 2
.
上述方程为齐次微分 方程,可用变量变换 法求解。
27
小结 1.变量分离方程的形状 dy f ( x) ( y )或M 1 ( x) N1 ( y ) dx M 2 ( x) N 2 ( y ) dy 0 dx 2.变量分离方程的求解:分离变量法 步骤:分离变量,两边积分,检查是否有遗漏的特解
2
(*)
23
分离变量,得 dX 1 u du 2 X 1 2u u 两边积分,得 ~ 2 2 ln X ln | u 2u 1 | C
即X (u 2u 1) C1 (C1 e ), 此外容易验证 u 2 2u 1 0 亦为方程(*)的解,因此方程(*)的通解为 X 2 (u 2 2u 1) C1, 其中C1为任意常数。
常微分第一二章
u2070C
例2.解析几何模型
求一曲 线线 上, 任 (x设 y,)一 处 在 点 的 曲 切线斜
该点横2倍 坐, 标并 的且该 (1曲 .,2)线经过点
解 (1)依据: 曲y线 y(x)导数的几 — — 何 d是 y意曲 义线
dx 在 (x,y)处的切线斜率
(2)数学模型
dy dx
2x .
4. 40学时的初步分配:
第一、二章:用6次课;(穿插课堂习题)
第 三 章:用4次课;
第 四 章:用6次课;
第 五 章:用4次课.
第二章 一阶微分方程的初等解法
一阶微分方程
六种解法(工具箱)
显式
dy f (x, y), d(x导数已解出
)
(M对( x称, y形)d式x) N ( x, y)dy 0
3. 重视应用,重视计算,重视解题格式; 4. 适度减弱教材中某些理论证明,补充工程中常用 的,工程师喜用的解题方法.
学习要求
1. 准确、熟练地掌握基本概念、基本解法,了解 相关的基本理论. 2. 初步学会由实际问题建立数学模型、求解、再 回到(解释、解决)实际问题的方法.
3. 认真听好课,及时预习和复习;上好习题课,按 质、按量及时完成作业.
y | x 1 2
(3)求解
通解: yx2 c
曲线族
特解: yx2 1 过(1,2)满足题意的一条曲线
例3.单摆运动数学模型
一根长为l 的细杆,一端联结一个质量为m 的球M, 另
一端悬挂在O点,若不计细杆的质量,在重力的作用下细
杆在某一铅直平面上摆动,求摆球的运动规律.
解 (1)物理依据:牛顿第二定律——Fma
画出方程 dy x2 y2 的线素场,并近似地描出积分曲
常微分方程(王高雄)第三版 2.1教学教材
(I)齐次方程
ddyxg(yx)
(II) 形如 ddyxfaa21xxbb12yycc12的方,程 其中 a1,b1,c1,a2,b2,c2为任意.常数
(I) 形如
dyg(y) dx x
(2.5)
方程称为齐次方程, 这里g(u)是u的连续函. 数
求解方法: 10 作变量代换(引入新变量)u y ,方程化为
x
du g(u)u, (这里d由 yx于 duu)
dx x
dx dx
20 解以上的变量分离方程
30 变量还原.
例4 求解方程 xdy 2xyy dx
(x0)
解: 方程变形为 dy2 yy dx x x
(x0)
这是齐次方程, 令u y 代入得 x
x du u 2 uu 即 x du 2 u
dx
为 (1)的情形,可化为变量分离方程求解.
解的步骤:
10解方 程 aa21xx 组 bb1 2yy cc1200,
得解yx
,
20 作变换 YXyx,方程化为
dY a1Xb1Y dX a2Xb2Y
g
(
Y X
)
30再经变 u换 Y,将以上方程化离 为方 变程 量分
X
40 求解
50 变量还原
dx
10 分离变量, 当 (y)0时 ,将 (2.1)写成
dy f (x)dx,
(y)
这样变量就“分离”开了.
20 两边积分得
dy
(y)f(x)d xc (2.2)
1 的某一原函数 f (x)的某一原函数 ( y)
由 (2.2)所确定 y的 (x,c)就 函 (2 为 .数 1)的.解
例:
分离变量:
丁同仁常微分方程教程第二章小结
−( Q
∂M ∂N − ) ∂y ∂x = −1 M
仅与 y 有关,故方程有积分因子 µ ( y ) = e ∫
− dy
= e− y
从而方程 −e − y + e− y ( x + 2 y )dy = 0 为全微分方程 其通解为: − ∫x e − y dx + ∫y ( x0 + 2 y )dy = c1
z =y−y1
Riccati方程不恒为零。观察其一 特解为y = ϕ ( x)
1
dy = p ( x) y 2 + q ( x ) y + r ( x), dx p( x)
习题课 目的:①对 1-5 节内容即初等积分法求解方法的归纳,并通过几个 典型例题的讲解,使大家认识到一阶方程解法的灵活性,多样性。 ②关于一阶方程的通解公式和通解结构定理, 也是一重要的
变量分离方程
∆=0
dy = f ( x)ϕ ( y ) dx
y = c(x)e ∫
− p( x)dx
线性方程
dy = p ( x) y + q ( x) dx
1 µ= ϕ1(y)
dy a x + b1 y + c1 = f( 1 ) dx a2 x + b2 y + c2
a1 b1 a2 b2
µ=
u= y x
从而 δ lim ce x →+∞
∫x0 a (τ ) dτ
≤ lim c e− k ( x − x0 ) = 0
δ x →+∞
再讨论前一项,首先,由积分性质
−
e
∫x0 a (τ ) dτ
x
∫
x
x0
常微分方程与动力系统第二章课后题参考答案
常微分方程与动力系统第二章习题参考答案 1.证明:因为()t Φ是线性齐次系统(LH )的一个基本解矩阵,由定理2.5知()t Φ在区间J 上满足矩阵微分系统()M LH ,即.()()()t A t t Φ=Φ,.1()()()A t t t -=ΦΦ所以由()A t 确定的线性齐次系统(LH )必唯一。
2.证明:因为()t ϕ,()t ψ分别是.()x A t x=和.()T x A t x =-的解,所以111()()()nk k k nnk k k a d t A t t dt a ϕϕϕϕ==⎛⎫ ⎪ ⎪== ⎪ ⎪⎪ ⎪⎝⎭∑∑ ,11211111122222*121()()()nn k k k n n kn kn n n nnk a a a a a a a d t A t t dta a a a ψψψψψψ==⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪ ⎪⎪ ⎪=-ψ=-=- ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪⎝⎭⎝⎭ ⎪⎝⎭∑∑ 因而1111112211(,)(,)(,),,nnk k k k k k nnkn k k nk k n n k a a d d d dt dt dt a a ψϕϕψψϕϕψϕψψϕψϕψϕ====⎡⎤⎛⎫⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥ ⎪-⎢⎥ ⎪⎪ ⎪⎢⎥ ⎪⎢⎥ ⎪ ⎪ ⎪⎢⎥=+= ⎪+⎢⎥ ⎪ ⎪ ⎪⎢⎥ ⎪-⎢⎥ ⎪⎪ ⎪⎢⎥⎪⎝⎭⎢⎥⎝⎭⎝⎭ ⎪⎣⎦⎢⎥⎝⎭⎣⎦∑∑∑∑ 11111111()0nnn n nnnnn n nnm m m m i ij j i ij j i mk k km k mk k km m m m m i j i j k k k k a a a a a a ϕψψϕϕψϕψϕψϕψ============-=+=-=-=∑∑∑∑∑∑∑∑∑∑∑∑所以(),()()()1nt t t t k kk ϕψϕψ≡≡∑=常数。
3.证明:设)t Φ(为系统.()x A t x=的一个基本解矩阵,则由定理2.11知[]1()Tt -Φ是系统.()Tx At x =-的基本解矩阵,由定理2.4知系统.()x A t x=满足初始条件00()x t x =的特解为100()))t t t x ϕ-=Φ(Φ(,[)0,0,t t ∈+∞由题可知)t Φ(与[]1()Tt -Φ在[)0,+∞上有界,从而由定理2.24知110()0k k t ∃=>和220()0k k t =>使得10120(),(),T t k t t t k t t -⎧Φ≤≤<+∞⎪⎨Φ≤≤<+∞⎪⎩,利用常数变易法公式(2.32),可知式.()()y A t y B t y=+的初始条件为00()y t y =的解满足1()()()()()()tt y t t t s B s y s ds ϕ-=+ΦΦ⎰因为1111()()(Ttttt---ΦΦ≤Φ所以12120()()(),tt y t k kx k k B s y s≤+≥⎰,利用格朗瓦尔不等式有12()120().tt k k B s dsy t k k x e⎰≤记12()12tt k k B s dsC k k e ⎰=设0()B t dt M +∞=<+∞⎰则()()tt B s ds B t dt M+∞≤=⎰⎰有1212k k MCk k e≤从而00(),y t C x t t ≤≥所以系统.()()y A t y B t y =+的一切解都在[)0,+∞上有界。
常微分方程线性方程全文
x(t)
4000000 20t
因此有
dx 100x 61.6, x(0) 0. dt 400000 2t
该方程有积分因子
(t) exp(
100 dt) (4000 0.02t)50
400000 2t
两边同乘以 (t) 后,整理得
d [x(4000 0.02t)50] 61.6(4000 0.02t)50 dt 积分得?
0
0
f (x)以 2 为周期 (令s 2 t )
x2 e a(2 s) f (s)ds x eat f (t)d t
0
2
1
c 1 e2a
0 eas f (s)ds
2
将c 的表达式代入通解,再一次利用f (x)的周期性得:
y
1 e 2 a
1
x 2 e a( x s) f (s)ds
将y和y代入原方程,
u( x)e p( x)dx p( x)u( x)e p( x)dx
p( x)u( x)e p( x)dx g( x)
u( x)e p( x)dx g( x)
8
u( x)dx g( x)e p( x)dx dx u( x) g( x)e p( x)dxdx C
解 原方程变形 dy 2sin y cos y 2cos2 y x 0
dx
22
2
1 2cos2
y
dy dx
tan
y 2
x0
2
d dx
tan
y 2
设tan y z 2
dz z x...... 一阶线性方程 dx
原方程通解: tan y C e x (1 x)
2
17
求 y 1 e y 的通解。
第二章 常微分方程
a0 a2
a1 a0 0 3 a3 a 2 a1 a5 a 4 a 3
a1 a n 3 a2 n 1, a0 a2
,
a1 a3 4 a5 a7
a0 0 0 a 2 a1 a0 a4 a3 a2 a 6 a5 a 4
a2 n 2 ,
0 0 an n 1 an
通解
n nt
yt c1 x e
或
1 1t
c2 x e
2 2t
cn x e
y=Yc
常数 c 由初始条件确定
第二章常微分方程——线性稳定性分析
四、线性稳定性分析方法
稳定性(stability)——系统的一种动态特性,指偏离定常 状态后能否自动返回该定常态的性质 ,系统抗干扰能力的
y e (c1 cos x c2 sin x)
x
重根
y e1x (c1 c2 x)
2。非其次方程特解:比较系数法
第二章常微分方程——二阶变系数方程
二、 二阶变系数方程的解法
1、级数解法
2 d y dy x 2 2 xF ( x ) G ( x ) y 0 dx dx
an 2 an ( n c k )( n c k )
(1) a0 (k 1) 1 a2n 2n 2 n!(n k )( n k 1)()(1 k ) 4 n! (n k 1)
n n
第一解
x ( 1) 2 y1 2k (k 1)a0 n 0 n! ( n k 1)
第二章常微分方程——线性稳定性分析
Routh指出,若采用如下的判定函数R i
第二章常微分方程
an (n c)(n c 1)xnc (F0 F1x F2 x2 ) an (n c)xnc
n0
n0
(G0 G1x G2 x2 ) an xnc 0
n0
第二章常微分方程——二阶变系数方程
首项xc的系数为0——指标方程
c2 (F0 1)c G0 0
第n项xn+c的系数为0 ——递推公式
rAs
)
dy dt
y
(rA
rAs )
[Qr (T )
Qr (Ts )]
第二章常微分方程——线性稳定性分析
将反应项与移热项线性展开
dx dt
1
rA cA
s
x
rA T
s
y
dy dt
rA cA
s
x
1
rA T
s
dQr dT
s
y
特征根方程
2 tr 0
detA I 0
从中可解出n个特征根和特征向量,构成基解矩阵
第二章常微分方程——一阶常系数方程组
通解 或
Y t e1t x 1 , e2t x 2 , ,ent x n
y t c1 x 1e1t c2 x 2e2t cn x nent
y=Yc 常数 c 由初始条件确定
y2
y c cc1
➢ 当c1-c2 为整数时,第二解为
y2
c
c
c2
y cc2
第二章常微分方程——二阶变系数方程
推导:设
y(x,c)
an不一定满足指标方程,将其代入
方程后有
x 2 d 2y dx 2
xF
(x
)
dy dx
G(x)y (c c1)(c c2)a0x c
《常微分方程指导与实验》第2章:一阶微分方程的解的存在定理
第二章 一阶微分方程的解的存在定理§2.1 一阶微分方程解的基本理论主要内容一 导数已解出方程初值问题解的存在唯一性定理 考虑导数已解出的一阶DE 的初值问题()()00,y f x y y x y '=⎧⎪⎨=⎪⎩(2.1)(2.2)这里()y x f ,是在闭矩形域R : a x x ≤-0,b y y ≤-0上的连续函数。
定义2.1 如果存在常数0>L ,使得对于所有的点()1,y x ,()2,y x R ∈,都有不等式()()2121,,y y L y x f y x f -≤-成立,则称函数()y x f ,在R 上关于y 满足李普希兹(Lipschitz )条件。
1定理2.1 (毕卡存在唯一性定理) 如果()y x f ,在R 上满足条件: 1)连续;2)关于y 满足李普希兹条件,则初值(2.1)和(2.2)在区间h x x ≤-0上存在唯一解()x y y =,其中()M b a h ,m in=,()y x f M R y x ,max ),(∈=。
注1 取数h 的意义。
注意到()y x f M R y x ,max ),(∈=,从而积分曲线()x y y =在任一点()()R x y x ∈,处的切线斜率()M x y ≤'。
于是从点()o y x p ,0引两条斜率分别为M 和M -的直线1l 和2l ,便知过点P 的积分曲线必限制在图2.1和图2.2的阴影区域内。
而直线1l 和2l 相交情形有如下两种可能。
(i )若相交成如图 2.1所示的情况,则a Mb>,积分曲线()x y y =在a x x ≤-0上不越出R ,从而应取a h =。
(ii )若相交成如图 2.2所示的情况,则a Mb >,积分曲线()x y y =在Mb x x ≤-0上不越出R ,从而应取Mb h =。
总之,取()M ba h ,min =,就是为了使初值问题(2.1)和(2.2)的解在h x x ≤-0上总存在。
常微分方程第二章 一阶微分方程的初等解法
du dx 1u2 x
两边积分得: ln u 1 u2 ln x ln c
整理后得 u 1 u2 cx
变量还原得 y 1 ( y )2 cx
x
x
du dx 1u2 x
最后由初始条件 y(1) 0,可定出c 1.
故初值问题的解为 y 1 (x2 1) 2
可2、化d为y 变a量1x 分b1 y离 方c1 法
由对数的定义有
y e p( x)dxc1
y e p( x)dxc1
即
y ec1e p(x)dx ce p(x)dx.
此外y 0也是方程的解,若在上式中充许c 0, 即知y 0也包括在上式中,
故方程的通解为
y ce p(x)dx , c为任常数.
例4
求初值问题
dy dx
y2
c os x的特解.
例:
y y sin x 0
并求满足条件的 y( ) 2 特解。
2
线性微分方程
例:
1、cos x dy y sin x cos2 x dx
二 伯努利(Bernoulli )方程
形如 dy p(x) y Q(x) yn dx
的方程,称为伯努利方程. 这里P(x), Q(x)为x的连续函数 。
故对应齐次方程通解为 y c(x 1)n
y
ce p(x)dx
ce
n dx x 1
c(x
1)n
其次应用常数变易法求非齐线性方程的通解,
令y c(x)( x 1)n为原方程的通解 , 代入得
dc(x) (x 1)n nc(x)(x 1)n1 nc(x)(x 1)n1 ex (x 1)n dx
解的步骤:
10
解方程组aa21xx
常微分方程教程+第二版+丁同仁+李承志+答案和练习第2章习题 第二章答案
习 题 2-1判断下列方程是否为恰当方程,并且对恰当方程求解: 1.0)12()13(2=++-dy x dx x解:13),(2-=x y x P , 12),(+=x y x Q ,则0=∂∂y P ,2=∂∂x Q , 所以 xQy P ∂∂≠∂∂ 即 原方程不是恰当方程.2.0)2()2(=+++dy y x dx y x解:,2),(y x y x P += ,2),(y x y x Q -=则,2=∂∂y P ,2=∂∂x Q 所以xQy P ∂∂=∂∂,即 原方程为恰当方程 则,0)22(=-++ydy xdy ydx xdx两边积分得:.22222C y xy x =-+ 3.0)()(=+++dy cy bx dx by ax (a,b 和c 为常数). 解:,),(by ax y x P += ,),(cy bx y x Q +=则,b y P =∂∂,b x Q =∂∂ 所以xQy P ∂∂=∂∂,即 原方程为恰当方程 则,0=+++cydy bxdy bydx axdx两边积分得:.2222C cy bxy ax =++ 4.)0(0)()(≠=-+-b dy cy bx dx by ax解:,),(by ax y x P -= ,),(cy bx y x Q -=则,b y P -=∂∂,b x Q =∂∂ 因为 0≠b , 所以xQ y P ∂∂≠∂∂,即 原方程不为恰当方程5.0sin 2cos )1(2=++udt t udu t解:,cos )1(),(2u t u t P += u t u t Q sin 2),(=则,cos 2u t t P =∂∂,cos 2u t x Q =∂∂ 所以xQy P ∂∂=∂∂,即 原方程为恰当方程则,0cos )sin 2cos (2=++udu udt t udu t两边积分得:.sin )1(2C u t =+ 6.0)2()2(2=++++dy xy e dx y e ye x x x解: xy e y x Q y e ye y x P x x x 2),(,2,(2+=++=,则,2y e y P x +=∂∂,2y e x Q x +=∂∂ 所以xQy P ∂∂=∂∂,即 原方程为恰当方程则,0])2()[(22=++++dy xy e dx y ye dx e x x x 两边积分得:.)2(2C xy e y x =++7.0)2(ln )(2=-++dy y x dx x xy解:,2ln ),(),(2y x y x Q x xy y x P -=+=则,1x y P =∂∂,1x x Q =∂∂ 所以xQy P ∂∂=∂∂,即 原方程为恰当方程则02)ln (2=-++ydy dx x xdy dx xy两边积分得:23ln 3y x y x -+.C = 8.),(0)(22为常数和c b a cxydy dx by ax =++解:,),(,),(22cxy y x Q by ax y x P =+=则,2by y P =∂∂,cy x Q =∂∂ 所以 当xQy P ∂∂=∂∂,即 c b =2时, 原方程为恰当方程则0)(22=++cxydy dx by dx ax两边积分得:233bxy ax +.C = 而当c b ≠2时原方程不是恰当方程.9.01222=-+-dt ts s ds t s 解:,),(,12),(22ts s s t Q t s s t P -=-= 则,212t s t P -=∂∂,212t s s Q -=∂∂ 所以xQ y P ∂∂=∂∂, 即原方程为恰当方程,两边积分得:C ts s =-2.10.,0)()(2222=+++dy y x yf dx y x xf 其中)(⋅f 是连续的可微函数.解:),(),(),(),(2222y x yf y x Q y x xf y x P +=+=则,2f xy y P '=∂∂,2f xy x Q '=∂∂ 所以xQy P ∂∂=∂∂, 即原方程为恰当方程,两边积分得:22()f xy dx C +=⎰,即原方程的解为C y x F =+)(22 (其中F 为f 的原积分).习 题 2-2. 1. 求解下列微分方程,并指出这些方程在平面上的有意义 的区域::(1)yx dx dy 2=解:原方程即为:dx x ydy 2= 两边积分得:0,2332≠=-y C x y .(2))1(32x y x dx dy += 解:原方程即为:dx xx ydy 321+=两边积分得:1,0,1ln 2332-≠≠=+-x y C x y .(3)0sin 2=+x y dxdy解: 当0≠y 时原方程为:0sin 2=+xdx y dy两边积分得:0)cos (1=++y x c .又y=0也是方程的解,包含在通解中,则方程的通解为0)cos (1=++y x c .(4)221xy y x dxdy+++=; 解:原方程即为:2(1)1dyx dx y=++ 两边积分得:c x x arctgy ++=22, 即 )2(2c x x tg y ++=. (5)2)2cos (cos y x dxdy= 解:①当02cos ≠y 时原方程即为:dx x y dy 22)(cos )2(cos = 两边积分得:2222sin 2tg y x x c --=. ②y 2cos =0,即42ππ+=k y 也是方程的解. (N k ∈) (6)21y dxdyx-= 解:①当1±≠y 时 原方程即为:xdx y dy =-21 两边积分得:c x y =-ln arcsin . ② 1±=y 也是方程的解.(7).yxe y e x dx dy +-=- 解.原方程即为:dx e x dy e y xy)()(--=+两边积分得:c e x e y x y ++=+-2222, 原方程的解为:c e e x y x y =-+--)(222.2. 解下列微分方程的初值问题.(1),03cos 2sin =+ydy xdx 3)2(ππ=y ;解:两边积分得:c yx =+-33sin 22cos , 即 c x y =-2cos 33sin 2因为 3)2(ππ=y , 所以 3=c .所以原方程满足初值问题的解为:32cos 33sin 2=-x y .(2).0=+-dy ye xdx x, 1)0(=y ;解:原方程即为:0=+ydy dx xe x,两边积分得:c dy y dx e x x=+-2)1(2, 因为1)0(=y , 所以21-=c , 所以原方程满足初值问题的解为:01)1(22=++-dy y dx e x x .(3).r d dr=θ, 2)0(=r ; 解:原方程即为:θd rdr=,两边积分得:c r =-θln , 因为2)0(=r , 所以2ln =c ,所以原方程满足初值问题的解为:2ln ln =-θr 即θe r 2=.(4).,1ln 2yx dx dy+= 0)1(=y ; 解:原方程即为:dx x dy y ln )1(2=+,两边积分得:3ln 3y y x x x c ++-=, 因为0)1(=y , 所以1=c ,所以原方程满足初值为:3ln 13y y x x x ++-=(5).321xy dxdyx=+, 1)0(=y ; 解:原方程即为:dx xx y dy 231+=, 两边积分得:c x y ++=--22121, 因为1)0(=y , 所以23-=c ,所以原方程满足初值问题的解为:311222=++yx .1. 解下列微分方程,并作出相应积分曲线的简图. (1).x dxdycos = 解:两边积分得:c x y +=sin . 积分曲线的简图如下:(2).ay dxdy=, (常数0≠a ); 解:①当0≠y 时,原方程即为:dx aydy= 积分得:c x y a +=ln 1,即 )0(>=c ce y ax②0=y 也是方程的解. 积分曲线的简图如下:y(3).21y dxdy-=; 解:①当1±≠y 时,原方程即为:dx y dy =-)1(2 积分得:c x yy+=-+211ln ,即 1122+-=x x ce ce y .②1±=y 也是方程的解.积分曲线的简图如下:(4).n y dx dy =, )2,1,31(=n ; 解:①0≠y 时,ⅰ)2,31=n 时,原方程即为 dx ydyn =, 积分得:c y n x n=-+-111.ⅱ)1=n 时,原方程即为dx ydy=积分得:c x y +=ln ,即)0(>=c ce y x.②0=y 也是方程的解.积分曲线的简图如下:4. 跟踪:设某A 从xoy 平面上的原点出发,沿x 轴正方向前进;同时某B 从点开始跟踪A ,即B 与A 永远保持等距b .试求B 的光滑运动轨迹.解:设B 的运动轨迹为)(x y y =,由题意及导数的几何意义,则有22yb ydx dy --=,所以求B 的运动轨迹即是求此微分方程满足b y =)0(的解.解之得:222222ln 21y b y b b y b b b x ----++=.5. 设微分方程)(y f dxdy=(2.27),其中f(y) 在a y =的某邻域(例如,区间ε<-a y )内连续,而且a y y f =⇔=0)(,则在直线a y =上的每一点,方程(2.27)的解局部唯一,当且仅当瑕积分∞=⎰±εa ay f dy)((发散). 证明:(⇒)首先经过域1R :,+∞<<∞-x a y a <≤-ε 和域2R :,+∞<<∞-x ε+≤<a y a内任一点(00,y x )恰有方程(2.13)的一条积分曲线, 它由下式确定00)(x x y f dyyy-=⎰. (*) 这些积分曲线彼此不相交. 其次,域1R (2R )内的所有 积分曲线c x y f dy +=⎰)(都可由其中一条,比如0)(c x y f dy+=⎰ 沿着 x 轴的方向平移而得到。
常微分方程-总复习
dy a1 x b1 y c1 dx a2 x b2 y c2
dx
x
2.3 恰当方程和积分因子 2.3.1 恰当方程 定义、判别方法、求解方法 2.3.2 积分因子 定义、特殊类型方程的积分因子的求法 2.4 一阶隐方程和参数表示
第三章 一阶微分方程解的存在定理
解的存在唯一性定理的内容及证明过程。
近似计算和误差估计;
解对初值的可微性
第四章 高阶微分方程
4.1 线性微分方程的一般理论 4.1.1 齐线性方程解的性质与结构 定理2-定理6 4.1.2 非齐线性方程与常数变易法 定理7 常数变易法 4.2 常系数线性方程的解法 4.2.2 复值函数与复值解 复值函数的运算性质、定理8、定理9
4.2.2 常系数齐线性方程和欧拉方程 欧拉待定指数函数法、根据特征根的性质确定 方程的基本解组、欧拉方程的求解 4.2.3 非齐线性方程-比较系数法
第五章 线性微分方程组
5.1 解的存在唯一性定理 5.1.1 记号和定义 将n阶线性微分方程的初值问题化为等价的微分 方程组的初值问题 5.1.2 存在唯一性定理 5.2 线性微分方程组的一般理论 5.2.1 齐线性微分方程组
定理2-定理6 定理1*定理2* 5.2.2 非齐线性微分方程组 定理7 定理8 常数变易公式
常微分方程
总复习
第一章 绪论
基本概念 常微分方程、偏微分方程、微分方程的阶 线性和非线性微分方程 解:隐式解、通解、特解 积分曲线
第二章 一阶微分方程
2.1 变量分离方程和变量变换 2.1.1 变量分离方程 2.1.2 可化为变量分离方程的类型 y 1) dy g
2) 2.2 线性方程与常数变易法 一阶齐线性微分方程、一阶非齐线性微分方程、 伯努利方程
常微分方程第2章
(u ) u
两边积分, 得
积分后再用
du
(u ) u
代替 u, 便得原方程的通解.
例1. 解微分方程 y
解: 令 u
y
y x
tan
y x
.
x u x u u tan u
, 则 y u x u , 代入原方程得
分离变量 两边积分 得
cos u sin u cos u
说明: 显然 x = 0 , y = 0 , y = x 也是原方程的解, 但在
例3. 在制造探照灯反射镜面时, 要求点光源的光线反 射出去有良好的方向性 , 试求反射镜面的形状. 解: 设光源在坐标原点, 取x 轴平行于光线反射方向, 则反射镜面由曲线 绕 x 轴旋转而成 . 过曲线上任意点 M (x, y) 作切线 M T,
dX X
积分得
arctan u
2 1 ln (1 u ) 2
ln C X
代回原变量, 得原方程的通解:
y 5 ln C ( x 1) ln 1 arctan x 1 2 x 1
y5
1
2
得 C = 1 , 故所求特解为
u
ln (1 e
x y
(1 e ) e 1 e
u
u
u
du
所求通解:
) y C ( C 为任意常数 )
解法2 分离变量
e
y
e C
y
x
即
(e C )e 1 0
x
(C<0 )
例4. 已知放射性元素铀的衰变速度与当时未衰变原
子的含量 M 成正比, 已知 t = 0 时铀的含量为 衰变过程中铀含量 M(t) 随时间 t 的变化规律.
常微分方程考研讲义第二章 一阶微分方程的初等解法
第二章、一阶微分方程的初等解法[教学目标]1. 理解变量分离方程以及可化为变量分离方程的类型(齐次方程),熟练掌握变量分离方程的解法。
2. 理解一阶线性微分方程的类型,熟练掌握常数变易法及伯努力方程的求解。
3. 理解恰当方程的类型,掌握恰当方程的解法及简单积分因子的求法。
4. 理解一阶隐式方程的可积类型,掌握隐式方程的参数解法。
[教学重难点] 重点是一阶微分方程的各类初等解法 ,难点是积分因子的求法以及隐式方程的解法。
[教学方法] 讲授,实践。
[教学时间] 14学时[教学内容] 变量分离方程,齐次方程以及可化为变量分离方程类型,一阶线性微分方程及其常数变易法,伯努利方程,恰当方程及其积分因子法,隐式方程。
[考核目标]1.一阶微分方程的初等解法:变量分离法、一阶线性微分方程的常数变易法、恰当方程与积分因子法、一阶隐方程的参数解法。
2.会建立一阶微分方程并能求解。
§1 变量分离方程与变量变换1、 变量分离方程1) 变量分离方程形如 ()()dy f x g y dx= (或1122()()()()0M x N y dx M x N y dy +=) (2.1) 的方程,称为变量分离方程,其中函数()f x 和()g y 分别是,x y 的连续函数.2) 求解方法如果()0g y ≠,方程(2.1)可化为,()()dy f x dx g y = 这样变量就分离开了,两边积分,得到 ()()dy f x dx c g y =+⎰⎰ (2.2) 把,()()dy f x dx g y ⎰⎰分别理解为1,()()f x y ϕ的某一个原函数. 容易验证由(2.2)所确定的隐函数(,)y x c ϕ=满足方程(2.1).因而(2.2)是(2.1)的通解.如果存在0y 使0()0g y =,可知0y y =也是(2.1)的解.可能它不包含在方程的通解(2.2)中,必须予以补上.3) 例题例1 求解方程dy x dx y=- 解 将变量分离,得到ydy xdx =-两边积分,即得 22222y x c =-+ 因而,通解为22x y c += 这里的c 是任意的正常数.或解出显式形式y =例2 解方程 2cos dy y x dx= 并求满足初始条件:当0x =时.1y =的特解.解 将变量分离,得到2cos dy xdx y = 两边积分,即得 1sin x c y -=+ 因而,通解为 1sin y x c=-+ 这里的c 是任意的常数.此外,方程还有解0y =.为确定所求的特解,以0x =.1y =代入通解中确定常数c ,得到 1c =-因而,所求的特解为 11sin y x =- 例3 求方程 ()dy P x y dx= (2.3) 的通解,其中()P x 是x 的连续函数.解 将变量分离,得到()dy P x dx y= 两边积分,即得 ln ()y P x dx c=+⎰ 这里的c是任意常数.由对数的定义,即有 ()P x dx c y e +⎰=即()P x dx c y e e ⎰=±令c e c ±= ,得到()P x dx y ce ⎰= (2.4)此外,0y =也是(2.3)的解.如果在(2.4)中允许0c =,则0y =也就包括在(2.4)中,因而,(2.3)的通解为(2.4),其中c 是任意常数.注: 1.常数c 的选取保证(2.2)式有意义.2.方程的通解不一定是方程的全部解,有些通解包含了方程的所有解,有些通解不能包含方程的所有解.此时,还应求出不含在通解中的其它解, 即将遗漏的解要弥补上.3.微分方程的通解表示的是一族曲线,而特解表示的是满足特定条件00()y x y =的一个解,表示的是一条过点00(,)x y 的曲线.2、可化为变量分离方程的类型1).形如 dy y g dx x ⎛⎫= ⎪⎝⎭(2.5) 的方程,称为齐次方程,这里的()g u 是u 的连续函数.另外,ⅰ)对于方程 (,)(,)dy M x y dx N x y = 其中函数(,)M x y 和(,)N x y 都是x 和y 的m 次齐次函数,即对0t >有(,)(,)m M tx ty t M x y ≡ (,)(,)m N tx ty t N x y ≡ 事实上,取1t x=,则方程可改写成形如(2.5)的方程. (1,)(1,)(1,)(1,)m m y y x M M dy x x y y dx x N N x x== ⅱ)对方程 (,)dy f x y dx= 其中右端函数(,)f x y 是x 和y 的零次齐次函数,即对0t >有(,)(,)f tx ty f x y =则方程也可改写成形如(2.5)的方程(1,)dy y f dx x= 对齐次方程(2.5)利用变量替换可化为变量分离方程再求解.令y u x =(2.6) 即y ux =,于是 dy du x u dx dx=+ (2.7) 将(2.6)、(2.7)代入(2.5),则原方程变为 ()du xu g u dx += 整理后,得到 ()du g u u dx x-= (2.8) 方程(2.8)是一个可分离变量方程,按照变量分离法求解,然后将所求的解代回原变量,所得的解便是原方程(2.5)的解.例4 求解方程dy y y tg dx x x=+ 解 这是齐次方程,以,y dy du u x u x dx dx ==+代入,则原方程变为 du xu u tgu dx +=+ 即du tgu dx x= (2.9) 分离变量,即有 dx ctgudu x=两边积分,得到 ln sin ln u x c=+这里的c是任意的常数,整理后,得到 sin u cx = (2.10)此外,方程(2.9)还有解0tgu =,即sin 0u =. 如果(2.10)中允许0c =,则sin 0u =就包含在(2.10)中,这就是说,方程(2.9)的通解为(2.10).代回原来的变量,得到原方程的通解为 sin ycx x =例5 求解方程(0).dyx y x dx +=<解 将方程改写为(0)dyyx dx x =< 这是齐次方程,以,ydyduu x u x dx dx ==+代入,则原方程变为dux dx =(2.11) 分离变量,得到dxx =两边积分,得到(2.11)的通解ln()x c =-+即2[ln()](ln()0)u x c x c =-+-+> (2.12) 这里的c 是任意常数.此外,(2.11)还有解0u =注意,此解不包括在通解(2.12)中.代回原来的变量,即得原方程的通解2[ln()](ln()0)y x x c x c =-+-+>及解0y =.原方程的通解还可表为2[ln()],ln()0,0,x x c x c y ⎧-+-+>=⎨⎩它定义于整个负半轴上. 注:1.对于齐次方程dy y g dx x ⎛⎫= ⎪⎝⎭的求解方法关键的一步是令y u x =后,解出y ux =,再对两边求关于x 的导数得dy du u x dx dx =+,再将其代入齐次方程使方程变为关于,u x 的可分离方程.2.齐次方程也可以通过变换x v y=而化为变量分离方程.这时x vy =,再对两边求关于y 的导数得dx dv v y dy dy =+,将其代入齐次方程dx x f dy y ⎛⎫= ⎪⎝⎭使方程变为,v y 的可分离方程 小结:这一讲我们主要讲解了一阶微分方程的可分离变量法和齐次方程的dy y g dx x ⎛⎫= ⎪⎝⎭形状的解法.而这一齐次方程通过变量替换任然可化为可分离方程,因而,一定要熟练掌握可分离方程的解法.2)形如 111222a xb yc dy dx a x b y c ++=++ (2.13) 的方程经变量变换化为变量分离方程,这里的121212,,,,,a a b b c c 均为常数.分三种情况来讨论(1)120c c ==情形.这时方程(2.13)属齐次方程,有 1122a x b y dy y g dx a x b y x +⎛⎫== ⎪+⎝⎭此时,令y u x=,即可化为变量可分离方程. (2)11220a b a b =,即1122a b a b =的情形.设1122a b k a b ==,则方程可写成22122222()()()k a x b y c dy f a x b y dx a x b y c ++==+++ 令22a x b y u +=,则方程化为 22()du a b f u dx=+ 这是一变量分离方程.(3)1112220,a b c c a b ≠及不全为零的情形.这时方程(2.13)右端的分子、分母都是,x y 的一次式,因此11122200a xb yc a x b y c ++=⎧⎨++=⎩ (2.14) 代表xy 平面上两条相交的直线,设交点为(,)αβ.显然,0α≠或0β≠,否则必有120c c ==,这正是情形(1)(只需进行坐标平移,将坐标原点(0,0)移至(,)αβ就行了,若令X x Y y αβ=-⎧⎨=-⎩(2.15) 则(2.14)化为112200a X b Y a X b y +=⎧⎨+=⎩ 从而(2.13)变为 1122a X bY dY Y g dX a X b Y X +⎛⎫== ⎪+⎝⎭(2.16) 因此,得到这种情形求解的一般步骤如下:(1)解联立代数方程(2.14),设其解为,x y αβ==;(2)作变换(2.15)将方程化为齐次方程(2.16);(3)再经变换Y u X=将(2.16)化为变量分离方程; (4)求解上述变量分离方程,最后代回原变量可得原方程(2.13)的解.上述解题的方法和步骤也适用于比方程(2.13)更一般的方程类型 111222a x b y c dy f dx a x b y c ⎛⎫+== ⎪++⎝⎭此外,诸如 ()dy f ax by c dx++ ()()0y xy dx xg xy dy += 2()dy x f xy dx =2dy y xf dx x ⎛⎫= ⎪⎝⎭ 以及(,)()(,)()0M x y xdx ydy N x y xdy ydx ++-=(其中,M N 为,x y 的齐次函数,次数可以不相同)等一些方程类型,均可通过适当的变量变换化为变量分离方程.例6 求解方程 13dy x y dx x y -+=+- (2.17) 解 解方程组 1030x y x y -+=⎧⎨+-=⎩得1, 2.x y == 令12x X y Y =+⎧⎨=+⎩ 代入方程(2.17),则有dY X Y dX X Y-=+ (2.18) 再令Y u X =即 Y uX = 则(2.18)化为2112dX u du X u u +=-- 两边积分,得 22ln ln 21X u u c=-+-+ 因此22(21)c X u u e +-=±记1,c e c ±= 并代回原变量,就得2212Y XY X c +-=221(2)2(1)(2)(1)y x y x c -+----=此外,易验证2210u u +-=即2220Y XY X +-=也就是(2.18)的解.因此方程(2.17)的通解为22262y xy x y x c +---=其中c 为任意的常数. 3、 应用举例例7 电容器的充电和放电如图(2.1)所示的R C -电路,开始时电容C 上没有电荷,电容两端的电压为零.把开关K 合上“1”后,电池E 就对电容C 充电,电容C 两端的电压C u 逐渐升高,经过相当时间后,电容充电完毕,再把开关K 合上“2”,这时电容就开始放电过程,现在要求找出充、放电过程中,电容C 两端的电压C u 随时间t 的变化规律.解 对于充电过程,由闭合回路的基尔霍夫第二定理,c u RI E += (2.19)对于电容C 充电时,电容上的电量Q 逐渐增多,根据C Q Cu =,得到 ()C C du dQ d I Cu C dt dt dt=== (2.20) 将(2.20)代入(2.19),得到c u 满足的微分方程 c c du RC u E dt+= (2.21) 这里R 、C 、E 都是常数.方程(2.21)属于变量分离方程.将(2.21)分离变量,得到C C du dt u E RC =-- 两边积分,得到 11ln C u E t c RC -=-+ 即 1112t t c RC RC C u E e ec e ---=±= 这里12c c e =±为任意常数.将初始条件:0t =时,0C u =代入,得到2c E =-.所以 1(1)t RC C u E e -=- (2.22)这就是R C -电路充电过程中电容C 两端的电压的变化规律.由(2.22)知道,电压C u 从零开始逐渐增大,且当t →+∞时,C u E →,在电工学中,通常称RC τ=为时间常数,当3t τ=时,0.95C u E =,就是说,经过3τ的时间后,电容C 上的电压已达到外加电压的95%.实用上,通常认为这时电容C 的充电过程已基本结束.易见充电结果C u E =.对于放电过程的讨论,可以类似地进行.例8 探照灯反射镜面的形状在制造探照灯的反射镜面时,总是要求将点光源射出的光线平行地射出去,以保证照灯有良好的方向性,试求反射镜面的几何形状.解 取光源所在处为坐标原点,而x 轴平行于光的反射方向,设所求曲面由曲线()0y f x z =⎧⎨=⎩ (2.23) 绕x 轴旋转而成,则求反射镜面的问题归结为求xy 平面上的曲线()y f x =的问题,仅考虑0y >的部分,过曲线()y f x =上任一点(,)M x y 作切线NT ,则由光的反射定律:入射角等于反射角,容易推知12αα=从而 OM ON =注意到 2dy MP tg dx NPα==及,,OP x MP y OM ===就得到函数()y f x =所应满足的微分方程式dy dx = (2.24) 这是齐次方程.由2.12知引入新变量x u y =可将它化为变量分离方程.再经直接积分即可求得方程的解.对于方齐次方程(2.24)也可以通过变换x v y =而化为变量分离方程也可由x yv =得dx dv v y dy dy=+代入(2.24)得到sgn dv v y v y dy+=+于是sgn dy y y =(2.25)积分(2.25)并代回原来变量,经化简整理,最后得2(2)y c c x =+ (2.26)其中c 为任意常数.(2.26)就是所求的平面曲线,它是抛物线,因此,反射镜面的形状为旋转抛物面22(2)y z c c x +=+ (2.27)小结: 本节我们主要讨论了一阶可分离微分方程和齐次微分方程的求解问题.将各种类型的求解步骤记清楚的同时要注意对解的讨论.§2 线性方程与常数变易法1、一阶线性微分方程 ()()()0dy a x b x y c x dx++= 在()0a x ≠的区间上可以写成 ()()dy P x y Q x dx=+ (2.28) 对于()a x 有零点的情形分别在()0a x ≠的相应区间上讨论.这里假设(),()P x Q x 在考虑的区间上是x 的连续函数.若()0Q x ≡,(2.28)变为 ()dy P x y dx= (2.3) 称为一阶齐线性方程.若()0Q x ≠,(2.28)称为一阶非齐线性方程.2、常数变易法(2.3)是变量分离方程,已在例3中求得它的通解为()P x dx y ce ⎰= (2.4)这里c 是任意的常数.下面讨论一阶非齐线性方程(2.28)的求解方法.方程(2.3)与方程(2.28)两者既有联系又有区别,设想它们的解也有一定的联系,在(2.4)中c 恒为常数时,它不可能是(2.28)的解,要使(2.28)具有形如(2.4)的解, c 不再是常数,将是x 的待定函数()c x ,为此令()()P x dx y c x e ⎰= (2.29)两边微分,得到 ()()()()()P x dx P x dx dy dc x e c x P x e dx dx ⎰⎰=+ (2.30) 将(2.29)、(2.30)代入(2.28),得到()()()()()()()()()P x dx P x dx P x dx dc x e c x P x e P x c x e Q x dx ⎰⎰⎰+=+ 即()()()P x dx dc x Q x e dx -⎰= 积分后得到()()()P x dx c x Q x e dx c -⎰=+⎰(2.31) 这里c是任意的常数..将(2.31)代入(2.29),得到()()()()()() =()P x dx P x dx P x dx P x dx P x dx y e Q x e dx c ce e Q x e dx --⎛⎫⎰⎰=+ ⎪⎝⎭⎰⎰⎰+⎰⎰(2.32) 这就是方程(2.28)的通解.这种将常数变易为待定函数的方法,通常称为常数变易法.实际上常数变易法也是一种变量变换的方法.通过变换(2.29)可将方程(2.28)化为变量分离方程.注: 非齐线性方程的通解是它对应的齐线性方程的通解与它的某个特解之和.例1 求方程1(1)(1)x n dy x ny e x dx++-=+的通解,这里的n 为常数. 解 将方程改写为(1)1x n dy n y e x dx x -=++ (2.33) 先求对应的齐次方程01dy n y dx x -=+ 的通解,得(1)n y c x =+令 ()(1)n y c x x =+ (2.34)微分之,得到 ()(1)(1)()n dy dc x x n x c x dx dx=+++ (2.35) 以(2.34)、(2.35)代入(2.33),再积分,得()x c x e c=+ 将其代入公式(2.34),即得原方程的通解(1)()n x y x e c=++ 这里c是任意的常数. 例2 求方程22dy y dx x y =-的通解. 解 原方程改写为 2dx x y dy y=- (2.36) 把x 看作未知函数,y 看作自变量,这样,对于x 及dx dy 来说,方程(2.36)就是一个线性方程了.先求齐线性方程2dx x dy y= 的通解为2x cy = (2.37)令2()x c y y =,于是2()2()dx dc y y c y y dy dy=+ 代入(2.36),得到 ()ln c y y c=-+ 从而,原方程的通解为 2(ln )x y cy =- 这里c是任意的常数,另外0y =也是方程的解. 特别的,初值问题00()()()dy P x y Q x dx y x y ⎧=+⎪⎨⎪=⎩ 的解为0000()()()=()x x s x x x P d P d P d x x y cee Q s e ds ττττττ-⎰⎰⎰+⎰例3 试证 (1)一阶非齐线性方程(2.28)的任两解之差必为相应的齐线性方程(2.3)之解;(2)若()y y x =是(2.3)的非零解,而 ()y y x =是(2.28)的解,则(2.28)的通解可表为 ()()y cy x y x =+,其中c 为任意常数.(3)方程(2.3)任一解的常数倍或两解之和(或差)仍是方程(2.3)的解.证 (1)设12,y y 是非齐线性方程的两个不同的解,则应满足方程使 1122()(1)()(2)dy py Q x dx dy py Q x dx=+=+ (1)—(2)有1212()()d y y p y y dx-=- 说明非齐线性方程任意两个解的差12y y -是对应的齐次线性方程的解.(2)因为 (()())()()(()()()()d cy x y x dy x d y x c p cy p y Q x p cy y Q x dx dx dx+=+=++=++ 故结论成立.(3)因为12121212()()()(),(),()d y y d y y d cy p cy p y y p y y dx dx dx+-==+=- 故结论成立. 3、Bernoulli 方程形如 ()()n dy P x y Q x y dx=+ ( 0,1n ≠) (2.38) 的方程,称为伯努利(Bernoulli )方程,这里(),()P x Q x 为x 连续函数.利用变量变换可将伯努利方程化为线性方程来求解.事实上,对于0y ≠,用n y -乘(2.38)两边,得到 1()()nn dy y y P x Q x dx --=+ (2.39) 引入变量变换1n z y -= (2.40)从而 (1)n dz dy n y dx dx-=- (2.41) 将(2.40)、2.41)代入(2.39),得到 (1)()(1)()dz n P x z n Q x dx=-+- (2.42) 这是线性方程,用上面介绍的方法求得它的通解,然后再代回原来的变量,便得到(2.38)的通解.此外,当0n >时,方程还有解0y =.例4 求方程26dy y xy dx x=-的通解 解 这是2n =时的伯努利方程,令 1z y -=,得2dz dy y dx dx -=- 代入原方程得到 6dz z x dx x=-+ 这是线性方程,求得它的通解为 268c x z x =+ 代回原来的变量y ,得到 2618c x y x =+ 或者 688x x c y -= 这是原方程的通解. 此外,方程还有解0y =.例5 求方程331dy dx xy x y =+的解 解 将方程改写为 33dx yx y x dy=+ 这是一个自变量为y ,因变量为x 的伯努利方程.解法同上.例6 求方程23y dy e x dx x+=的通解 这个方程只要做一个变换,令,y y du dy u e e dx dx==,原方程改写为 22231du x u u dx x x=+便是伯努利方程.小结;这次主要讨论了一阶线性微分方程的解法.其核心思想是常数变易法.即将非齐线性方程对应的齐线性方程解的常数变易为待定函数,使其变易后的解函数代入非齐次线性方程,求出待定函数()c x ,求出非齐次方程的解.我们还讨论了伯努利方程,求解过程为,先变换,将原方程化为非齐线性方程,再求解.§3 恰当方程与积分因子1、恰当方程的定义将一阶微分方程(,)dy f x y dx = 写成微分的形式(,)0f x y dx dy -=把,x y 平等看待,对称形式的一阶微分方程的一般式为(,)(,)0M x y dx N x y dy += (2.43)假设(,),(,)M x y N x y 在某区域G 内是,x y 的连续函数,而且具有连续的一阶偏导数. 如果存在可微函数(,)u x y ,使得(,)(,)du M x y dx N x y dy =+ (2.44)即 (,), (,)u u M x y N x y x y∂∂==∂∂ (2.45) 则称方程(2.43)为恰当方程,或称全微分方程.在上述情形,方程(2.43)可写成(,)0du x y ≡,于是(,)u x y C ≡就是方程(2.43)的隐式通解,这里C 是任意常数(应使函数有意义).2、 恰当方程的判定准则定理1设(,),(,)M x y N x y 在某区域G 内连续可微,则方程(2.43)是恰当方程的充要条件是, (,)M N x y G y x∂∂=∈∂∂ (2.46) 而且当(2.46)成立时,相应的原函数可取为000(,)(,)(,)x y x y u x y M s y ds N x t dt =+⎰⎰ (2.47) 或者也可取为000(,)(,)(,)yxy x u x y N x t dt M s y ds =+⎰⎰ (2.48) 其中00(,)x y G ∈是任意取定的一点.证明 先证必要性.因为(2.43)是恰当方程,则有可微函数(,)u x y 满足(2.45),又知(,),(,)M x y N x y 是连续可微的,从而有 22M u u N y y x x y x∂∂∂∂===∂∂∂∂∂∂ 下面证明定理的充分性,即由条件(2.46),寻找函数(,)u x y ,使其适合方程(2.45).从(2.47)可知 (,)u N x y y∂=∂ 0000(,)(,) =(,)(,) =(,)(,)(,)y y y x y y y y u M x y N x t dt x x M x y N x t dt M x y M x t dt M x y ∂∂=+∂∂++=⎰⎰⎰ 即(2.45)成立,同理也可从(2.48)推出(2.45).例1. 解方程 21()02x xydx dy y++= (2.49)解 这里21, =()2x M xy N y=+,则y x M x N ==,所以(2.49)是恰当方程.因为N 于0y =处无意义,所以应分别在0y >和0y <区域上应用定理2.3,可按任意一条途径去求相应的原函数(,)u x y .先选取00(,)(0,1)x y =,代入公式(2.47)有 22011()ln 22xyx x u xdx dy y y y =++=+⎰⎰再选取00(,)(0,1)x y =-,代入公式(2.47)有22011()()ln()22xyx x u x dx dy y y y -=-++=+-⎰⎰可见不论0y >和0y <,都有2ln ||2x u y y =+故方程的通解为2ln ||2x y y C +=.3、恰当方程的解法上述定理已给出恰当方程的解法,下面给出恰当方程的另两种常用解法. 解法1. 已经验证方程为恰当方程,从(,)x u M x y =出发,有2(,)(,)()()2x u x y M x y dx y y y φφ≡+=+⎰ (2.50)其中()y φ为待定函数,再利用(,)y u N x y =,有221()22x x y yφ'+=+ 从而1()y y φ'=于是有 ()ln ||y y φ=只需要求出一个(,)u x y ,因而省略了积分常数.把它代入(2.50)便得方程的通解为2ln ||2x u y y C =+= 解法2. 分项组合的方法 对(2.49)式重新组合变为21()02x xydx dy dy y++=于是 2()ln ||02x d y d y +=从而得到方程的通解为 2ln ||2x y y C +=4、积分因子的定义及判别对于微分形式的微分方程(,)(,)0M x y dx N x y dy += (2.43)如果方程(2.43)不是恰当方程,而存在连续可微的函数(,)0x y μμ=≠,使得 (,)(,)0M x y dx N x y dy μμ+= (2.51) 为一恰当方程,即存在函数(,)v x y ,使(,)(,)M x y dx N x y dy dv μμ+≡则称(,)x y μ是方程(2.43)的积分因子.此时(,)v x y C =是(2.51)的通解,因而也就是(2.43)的通解.如果函数(,),(,)M x y N x y 和(,)x y μ都是连续可微的,则由恰当方程的判别准则知道, (,)x y μ为(2.43)积分因子的充要条件是M Ny xμμ∂∂=∂∂ 即 ()M N N M x y y xμμμ∂∂∂∂-=-∂∂∂∂ (2.52) 5、积分因子的求法方程(2.52)的非零解总是存在的,但这是一个以μ为未知函数的一阶线性偏微分方程,求解很困难,我们只求某些特殊情形的积分因子. 定理2 设(,),(,)M M x y N N x y ==和(,)x y ϕϕ=在某区域内都是连续可微的,则方程(2.43) 有形如((,))x y μμϕ=的积分因子的充要条件是:函数(,)(,)(,)(,)(,)(,)y x x y M x y N x y N x y x y M x y x y ϕϕ-- (2.53)仅是(,)x y φ的函数,此外,如果(2.53)仅是(,)x y φ的函数((,))f f x y ϕ=,而()()G u f u du =⎰,则函数((,))G x y e ϕμ= (2.54) 就是方程(2.43)的积分因子.证明 因为如果方程(2.43)有积分因子()μμϕ=,则由(2.52)进一步知()()d M N N M d x y y xμϕϕμϕ∂∂∂∂-=-∂∂∂∂ 即y x x yM N d d N M μϕμϕϕ-=-由()μμϕ=可知左端是ϕ的函数,可见右端y x x yM N N M ϕϕ--也是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,于是,有()d f d μϕϕμ=, 从而 ()()f d G e e ϕϕϕμ⎰==反之,如果(2.53)仅是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,则函数(2.54)是方程(2.52)的解.事实上,因为()()()()G x y y x NM N M f e M N x yϕμμϕϕϕμ∂∂-=-=-∂∂ 因此函数(2.54)的确是方程(2.43)的积分因子.为了方便应用这个定理,我们就若干特殊情形列简表如下:例2. 解22(31)()0y xy dx xy x dy -++-= 解 这里2231,M y xy N xy x =-+=-,注意y x M N y x -=-所以方程不是恰当的,但是1y xM N Nx-=它仅是依赖与x ,因此有积分因子1dxx ex μ⎰≡=给方程两边乘以因子x μ=得到2223(3)()0xy x y x dx x y x dy -++-=从而可得到隐式通解22321122u x y x y x C ≡-+= 例3. 解方程2()(1)0xy y dx xy y dy ++++=解 这里2,1M xy y N xy y =+=++方程不是恰当的.但是1y xM N My-=-- 它有仅依赖于y 的积分因子 11dyy eyμ-⎰≡=方程两边乘以积分因子1y μ=得到 1()(1)0x y dx x dy y++++=从而可得到隐式通解21ln ||2u x xy y y C ≡+++= 另外,还有特解0y =.它是用积分因子乘方程时丢失的解. 例4. 解方程 223(2)()0y x y dx xy x dy +++=解 这里2232,M y x y N xy x =+=+,不是恰当方程.设想方程有积分因子()x y αβμμ=,其中α,β是待定实数.于是2112111()(2)y xM N y x x N y M x y y x x y x yαβαβαβαβαβαβ----⋅=⋅=--+- 只须取3,2αβ==.由上述简表知原方程有积分因子32x y μ=从而容易求得其通解为:446313u x y x y C ≡+=六、积分因子的其他求法以例4为例,方程的积分因子也可以这样来求:把原方程改写为如下两组和的形式:223()(2)0y dx xydx x ydx x dy +++=前一组有积分因子11yμ=,并且 21()()y dx xydy d xy y+= 后一组有积分因子21xμ=,并且 2321(2)()x ydx x dy d x y x+= 设想原方程有积分因子211()()xy x y y xαβμ==其中α,β是待定实数.容易看出只须3,2αβ==,上述函数确实是积分因子,其实就是上面找到一个.例5. 解方程 1212()()()()0M x M y dx N x N y dy += 其中1M ,2M ,1N ,2N 均为连续函数.解 这里12()()M M x M y =,12()()N N x N y =.写成微商形式就形式上方程是变量可分离方程,若有0y 使得20()0M y =,则0y y =是此方程的解;若有0x 使得10()0N x =,则0x x =是此方程的解;若21()()0M y N x ≠,则有积分因子211()()M y N x μ=并且通解为1212()()()()M x N y u dx dy N x M y ≡+⎰⎰ 例6、试用积分因子法解线性方程(2.28).解 将(2.28)改写为微分方程[()()]0P x y Q x dx dy +-= (2.55)这里()(),1M P x y Q x N =+=-,而()M Ny xP x N∂∂-∂∂=- 则线性方程只有与x 有关的积分因子()P x dxe μ-⎰=方程(2.55)两边乘以()P x dxe μ-⎰=,得()()()()()0P x dx P x dx P x dxxP x e ydx e dy Q x e dx ---⎰⎰⎰-+= (2.56) (2.56)为恰当方程,又分项分组法()()()()0P x dx P x dxd ye Q x e dx --⎰⎰-=因此方程的通解为()()()P x dx P x dxye Q x e dx c --⎰⎰-=⎰即()()[()]P x dxP x dxy e Q x e dx c -⎰⎰=+⎰与前面所求得的结果一样.注:积分因子一般不容易求得可以先从求特殊形状的积分因子开始,或者通过观察法进行“分项分组”法求得积分因子.§4 一阶隐方程与参数表示1、一阶隐方程一阶隐式微分方程的一般形式可表示为:(,,)0F x y y '=如果能解出(,)y f x y '=,则可化为显式形式,根据前面的知识求解. 例如方程2()()0y x y y xy ''-++=,可化为y x '=或y y '=但难以从方程中解出y ',或即使解出y ',而其形式比较复杂,则宜采用引进参数的方法求解.一般隐式方程分为以下四种类型:1) (,)y f x y '= 2) (,)x f y y '= 3) (,)0F x y '= 4)(,)0F y y '=2、求解方法Ⅰ)可以解出y (或)x 的方程1) 讨论形如(,)y f x y '= (2.57)的方程的解法,假设函数(,)f x y '有连续的偏导数,引进参数y p '=,则方程(2.57)变为 (,)y f x p = (2.58)将(2.58) 的两边对x 求导数,得到 f f dpp x y dx∂∂=+∂∂ (2.59) 方程(2.59)是关于,x p 的一阶微分方程,而且属于显式形式.若求得(2.59)的通解形式为(,)p x c ϕ=,将其代入(2.58),于是得到(2.57)通解为 (,(,))y f x x c ϕ=若求得(2.59)的通解形式为(,)x p c ψ=,于是得到(2.57)的参数形式的通解为(,)((,),)x p c y f p c p ψψ=⎧⎨=⎩其中p 为参数, c 是任意常数.若求得(2.59)的通解形式为(,,)0x p c Φ=,于是得到(2.57)的参数形式的通解为(,,)0(,)x p c y f x p Φ=⎧⎨=⎩其中p 为参数, c 是任意常数. 例1 求方程3()20dy dyx y dx dx+-= 的解 解 令dyp dx=,于是有32y p xp =+ (2.60) 两边对x 求导数,得到 2322dp dp p p x p dx dx=++ 即 2320p dp xdp pdx ++= 当0p ≠时,上式有积分因子p μ=,从而32320p dp xpdp p dx ++= 由此可知4234p xp c += 得到42223344c pc x p p p -==-将其代入(2.60),即得43342()c p y p p-=+故参数形式的通解为22334(0) 212c x p p p c y p p ⎧=-⎪⎪≠⎨⎪=-⎪⎩当0p =时,由(2.60)可知0y =也是方程的解.例2 求方程22()2dy dy x y xdx dx =-+的解. 解 令dy p dx =,得到222x y p xp =-+ (2.61) 两边对x 求导数,得到2dp dp p px p x dx dx =--+ 或 (2)(1)0dpp x dx--= 由10dpdx -=,解得p x c =+,于是得到方程的通解为222x y cx c =++ (2.62)由20p x -=,解得2xp =,于是得到方程的一个解为24x y = (2.63)特解(2.63)与通解(2.62)中的每一条积分曲线均相切,因此称为方程的奇解.2) 讨论形如(,)dyx f y dx= (2.64) 的方程的求解方法,方程(2.64)与方程(2.57)的求解方法完全类似,假定函数(,)f y y ' 有连续偏导数.引进参数dyp dx=,则(2.64) 变为 (,)x f y p = (2.65) 将(2.65) 的两边对y 求导数,得到1f f dpp y x dy∂∂=+∂∂ (2.66) 方程(2.66))是关于,y p 的一阶微分方程,而且属于显式形式.设其通解为 (,,)0y p c Φ= 则(2.64)的通解为(,,)0(,)y p c x f y p Φ=⎧⎨=⎩Ⅱ)不显含y (或)x 的方程 3) 讨论形如(,)0F x y '= (2.67) 的方程的解法.记dyp y dx'==,此时(,)0F x p =表示的是xp 平面上的一条曲线,设曲线用参数形式表示为()x t ϕ=,()p t ψ= (2.68) 由于dy pdx =,进而()()dy t t dt ψϕ'= 两边积分,得到()()y t t dt c ψϕ'=+⎰ 于是得到方程(2.67)参数形式的解为()()()x t y t t dt c ϕψϕ=⎧⎪⎨'=+⎪⎩⎰c 是任意常数.例3 求解方程3330x y xy ''+-=解 令y p tx '==,则由方程得 331t x t =+, 2331t p t =+ 于是 23339(12)(1)t t dy dt t -=+ 积分得到 23333329(12)314(1)2(1)t t t y dt c c t t -+=+=+++⎰ 故原方程参数形式的通解为: 3332313142(1)t x t t y c t ⎧=⎪+⎪⎨+⎪=+⎪+⎩4) 讨论形如(,)0F y y '= (2.69) 的方程,其解法与方程(2.67)的求解方法类似.记dy p y dx '==,此时(,)0F y p =表示的是yp 平面上的一条曲线,设曲线用参数形式表示为()y t ϕ=,()p t ψ= 由关系式dy pdx =可知 ()()t dt t dx ϕψ'=,于是0p ≠时,有 ()()t dx dt t ϕψ'=, ()()t x dt c t ϕψ'=+⎰ 故方程(2.69)的参数形式的通解 ()()()t x dt c t y t ϕψϕ'⎧=+⎪⎨⎪=⎩⎰c 是任意常数.此外,不难验证,若(,0)0F y =有实根y k =,则y k =也是方程的解. 例4 求解方程 22(1)(2)y y y ''-=-. 解 令2y yt '-=,则有222(1')y y y t -=由此可以得2'1y t =-,1y t t=+ 代入1dx dy p=,得到 222111(1)1dx dt dt t t t=-+=-- 积分,得到1x c t=+ 故原方程参数形式的通解为 11x c t y t t ⎧=+⎪⎪⎨⎪=+⎪⎩其中c 是任意常数.此外, 当0y '=时原方程变为24y =,于是2y =±也是方程的解. 例5求解方程y '=解 令y p '=,则有p =,取,(,)22p tgt t ππ=∈-,则sin sec tgt x t t === 由dy pdx =得到cos sin dy tgt tdt tdt == 所以cos y t c =-+故原方程参数形式的通解为 sin cos x t y t c =⎧⎨=-+⎩其中c是任意常数.。
常微分方程教案(王高雄)第二章
第二章目录内容提要及其它 (1)第二章一阶微分方程的初等解法(初等积分) (2)第一节变量分离方程与变量变换 (2)一、变量分离方程 (2)二、可化为变量分离方程的类型 (6)1、齐次方程 (6)2、可化为变量分离方程 (7)三、应用例题选讲 (10)第二节线性方程与常数变易法 (11)第三节恰当方程与积分因子 (15)一、恰当方程 (15)二、积分因子 (20)第四节一阶隐含方程与参数表示 (23)一、可以解出y(或x)的方程 (24)二、不显含y(或x)的方程 (25)本章小结及其它 (27)内容提要及其它授课题目(章、节)第二章:一阶微分方程的初等解法教材及主要参考书(注明页数)教材:常微分方程(第三版),王高雄等,高等教育出版社,2006年,p30-74主要参考书:[1]常微分方程,东北师范大学微分方程教研室编,高等教育出版社,2005,p1-70[2]常微分方程教程,丁同仁等编,高等教育出版社,1991,p1-20[3]偏微分方程数值解法(第2版),陆金甫关治,清华大学出版社,2004,p1-12[4]常微分方程习题解,庄万主编,山东科学技术出版社,2003,p28-169[5]微分方程模型与混沌,王树禾编著,中国科学技术大学出版社,1999,p15-158[6]差分方程和常微分方程,阮炯编著,复旦大学出版社,2002,p38-124目的与要求:掌握变量分离方程、齐次方程、线性方程、伯努利方程和恰当方程的解法.理解变量变换思想方法和积分因子方法,并能应用于求解一些特殊的常微分方程.掌握四类典型的一阶隐方程的解法.能熟练求解变量分离方程、齐次方程、线性方程、伯努利方程、恰当方程和四类典型的一阶隐方程.领会变量变换思想方法和积分因子方法,并能应用于求解一些特殊的常微分方程.教学内容与时间安排、教学方法、教学手段:教学内容:第1节变量分离方程与变量变换;第2节线性方程与常数变易法;第3节恰当方程与积分因子;第4节一阶隐方程与参数表示:可以解出(或y x)的方程、不显含(或y x)的方程.时间安排:8学时教学方法:讲解方法教学手段:传统教学方法与多媒体教学相结合。
常微分方程21变量分离方程与变量变换
1
dx
1
x
两边积分得: 2y 2 lnxc1
整理后得通解为:
y
(ln
x
4 c1)2
(ln
4 cx )2
,
其c中 ec1,由于 y32函 x1在 x数 0无意 , 义
故此解 x只 0或 x在 0之一中.有意义
2此 020/1/2外 3 还y有 0解 ,这个 常微分解 方程 未包含,在 应通 补.解 上中
( y)
由 (2.2)所确定 y的 (x,c)就 函 (2 为 .数 1)的.解
2020/1/23
常微分方程
例:
dyx2 y2 1
dx
分离变量:
dy y2 1
x2dx
两边积分:
dy
y2
1
x2d
xC
arc tayn1x3 C 3
2020/1/23
常微分方程
注: 若存 y0,使 在 (y0)0,则 yy0也(是 2.1)的,解 可能
解: 设在t雪 时球 刻的v(体 t)表 , 积 面 为 s(积 t)则 , 为
2020/1/23
dv(t) ks(t) dt
根据球体的体积和表面积的关系得
122
s(t)(4)333v3(t) 常微分方程
12
引入新常 (数 4)333k,再利用题中条件得
dvk(4)13332v32
解:
解方程组
x y 1 0 x y 3 0
得 x1,y2,
令Xx1,Yy2代入方程得
dY X Y
1
Y X
dX X Y 1 Y
令u Y ,得 X du 1u2
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
两边同时关于 x 求不定积分:
dy( x ) ( y( x )) ( x )dx
写出通解:(结果含有一个任意常数)
( x) ( x) C
® 需要注意 ( y ) 的零点 y y0有可能
是方程的解,不可遗漏,见以下分析。
例1 解方程 y (sinx )(siny ) .
0
(e x y)dx ( x 2 sin y)dy 0
是恰当方程,并求它的通解.
解: 这里M ( x, y) e y, N ( x, y) x 2sin y. N ( x, y ) M ( x, y) , 所以 1 x y 故所给方程是恰当方程.
x
由于所求函数 u( x, y)满足 u u x e y, x 2 sin y, x y
du( x, y) M ( x, y)dx N ( x, y)dy,
从而(1)为恰当方程。
这时, 取( x0 , y0 ) R, 则
u( x, y)
( x, y )
( x0 , y0 )
M ( x, y)dx N ( x, y)dy
y y0
M ( x, y0 )dx N ( x, y)dy,
二、恰当方程的求解
1 不定积分法
10 判断M ( x, y)dx N ( x, y)dy 0是否为恰当方程 , 若是进入下一步 .
20 求u ( x, y ) M ( x, y )dx ( y ),
u 3 由 N ( x, y )求 ( y ). y 例1 验证方程
N N M [ M ( x, y )dx] 0. x y x x y
于是, (7)右端的确只含有 y, 积分之得
( y) [ N M ( x, y)dx]dy, y d ( y )
dy N M ( x, y )dx y
sin x x y y 4.如下方法: M ( x, y ) N ( x, y) 由于 , y x 由数学分析曲线积分与路径无关的定理知:
M ( x, y)dx N ( x, y)dy为某函数u( x, y)的全微分,
即有函数u( x, y),使
故
(7 )
u ( x, y ) M ( x, y )dx [ N M ( x, y)dx]dy, (8) y 即u( x, y)存在, 从而(1)为恰当方程 。
注:若(1)为恰当方程,则其通解为
M ( x, y)dx [ N y M ( x, y)dx]dy c, c为任常数
N ( x, y ) M ( x, y ) y , cos x 2 xe x y
故所给方程是恰当方程.
由于M ( x, y), N ( x, y)在 全平面上连续 ,
故取( x0 , y0 ) (0,0),则
M ( x, y ) y cos x 2 xe
y
N ( x, y ) sin x x 2e y 2,
即应满足
(4)
u ( x, y ) M ( x, y )dx ( y ).
u u ( x, y ) M ( x, y )dx ( y ). N ( x, y), (6) y 这里 ( y)是y的任意可微函数 ,
d ( y ) N M ( x, y )dx 因此 dy y
ydx xdy 1 x y d (ln ). 2 2 x y 2 x y
例2 求方程 (3x 2 6xy2 )dx (6x 2 y 4 y3 )dy 0 的通解. 解:
这里M ( x, y) 3x2 6xy 2 , N ( x, y) 6x2 y 4 y3 ,
d ( y ) x x 2 sin y dy d ( y ) 2 sin y 即 dy
( y) 2 cos y, x 故 u( x, y) e yx 2 cos y.
积分后得: 从而方程的通解为
e yx 2 cos y c.
x
2 分组凑微法 采用“分项组合”的方法,把本身已构成全微 分的项分出来,再把余的项凑成全微分.
由偏导数的定义 , 只要将y看作常数 , 将e x y对x积分得
u ( x, y ) (e y )dx ( y ) e yx ( y).
x
x
u( x, y) e yx ( y).
x
对u( x, y)关于y求偏导数, 得 ( y)应满足的方程为
( x, y )
u( x, y)
( 0, 0 )
x
M ( x, y)dx N ( x, y)dy
y
y
2 xdx 0 (sin x x 2e y 2)dy
M ( x,0)dx N ( x, y)dy
0 x 0
0 2
x y sin x x 2 (e y 1) 2 y. y sin x x 2e y 2 y. . 故通解为: y sin x x 2e y 2 y c, c为任常数
2 方程为恰当方程的充要条件 定理1 设函数M ( x, y)和N ( x, y)在一个矩形区
域R中连续且有连续的一阶 偏导数, 则方程
M ( x, y)dx N ( x, y)dy 0,
为恰当方程的充要条件是
(1)
(2).
M ( x, y ) N ( x, y ) , y x
证明 “必要性”设(1)是恰当方程, 则有函数u( x, y),使得
下面选择 ( y),使u同时满足(6),即 u d ( y ) M ( x, y)dx N dy y y
(7 )
下面证明 (7)的右端与x无关, 即对x的偏导数常等于零 [ N M ( x, y )dx] 事实上 x y N [ M ( x, y )dx] x x y
N ( x, y ) M ( x, y ) , 所以 12 xy x y
故所给方程是恰当方程. 把方程重新“分项组合”得
3x dx 4 y dy (6xy dx 6x ydy) 0
2 3 2 2
即
dx3 dy4 (3 y 2dx2 3x 2dy2 ) 0 d ( x y 3x y ) 0
1 2 1 2 2 1 2 d sin x d x y d y 0 2 2 2 或写成 d (sin2 x x 2 y 2 y 2 ) 0,
故通解为:
sin 2 x x 2 y 2 y 2 c,
由初始条件 y(0) 2, 得 c 4,
故所求的初值问题的解为:
dy 解 分离变量 sin xdx, sin y
csc y dy sin xdx
ln | csc y cot y | cos x C1 ,
y 1 cos y cos x C1 cos x C1 tan e 即 e 2 sin y y cos x tan Ce 是原方程的通解。 2
N ( x, y ) M ( x, y ) , 2 xy x y
故所给方程是恰当方程.
2
把方程重新“分项组合”得
2
cos x sin xdx ( xy dx x ydy) ydy 0,
即
1 2 d 1 x2 y2 d 1 y2 0, d sin x 2 2 2
故
M ( x, y ) N ( x, y) . y x
M ( x, y ) N ( x, y ) , “充分性” 若 y x
则需构造函数 u( x, y),满足
du( x, y) M ( x, y)dx N ( x, y)dy,
u M ( x, y ), (5) x u N ( x, y), (6) y 从(5)出发, 把y看作参数, 解这个方程得
dy 1 y2 例2 求微分方程 的通解. 2 dx (1 x ) xy
解 分离变量 y dx 2 dy ( 1 y 0) 2 2 1 y x(1 x ) y dx dy 两边积分得 2 2 1 y x(1 x )
y 1 x x 1 y 2 dy x(1 x 2 ) dx
u u du( x, y ) dx dy M ( x, y)dx N ( x, y)dy x y u u N ( x, y ) M ( x, y ), 故有 y x 2 2 N u M u . 从而 , x xy y yx
2u 2u 2u 2u , 由于 和 都是连续的 , 从而有 yx xy yx xy
§2.2
变量分离方程
一、 可分离变量的微分方程 一阶微分方程的形式: dy f ( x , y ) 或 y f ( x , y ) dx 可用分离变量法解的一阶方程的一般形式 dy (1) ( x ) ( y ) dx
具体解法
分离变量: 设 ( y ) 0
dy ( x )dx ( y)
3 4 2 2
或写成
故通解为:
x y 3x y c, c为任常数 。
3 4 2 2
例3 验证方程
(cosx sin x xy2 )dx y(1 x2 )dy 0,
是恰当方程,并求它满足初始条件y(0)=2的解.
解: 这里M ( x, y) cos x sin x xy2 , N ( x, y) y(1 x 2 ),
第二章 初等积分法
§2.1 恰当方程
1 恰当方程的定义 u( x, y ), 使得 定义1 若有函数
du( x, y) M ( x, y)dx N ( x, y)dy