常微分方程 第二章.
常微分方程二阶线性微分方程习题课

二阶线性微分方程
例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,+∞上有界。
- 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