微分方程第2章习题解
常微分方程第三版习题答案
常微分方程第三版习题答案常微分方程是数学中的一个重要分支,它研究的是描述自然界中变化规律的方程。
在学习常微分方程的过程中,习题是非常重要的一部分,通过解习题可以加深对理论知识的理解和应用能力的培养。
本文将为大家提供《常微分方程第三版》习题的部分答案,希望能对大家的学习有所帮助。
1. 习题一1.1 解:首先,我们根据题意列出方程:$\frac{dy}{dt} = 2y + t^2$这是一个一阶线性常微分方程,我们可以使用常数变易法来求解。
令$y = u(t)e^{2t}$,则$\frac{dy}{dt} = \frac{du}{dt}e^{2t} + 2ue^{2t}$将上述结果代入原方程,得到:$\frac{du}{dt}e^{2t} + 2ue^{2t} = 2(u(t)e^{2t}) + t^2$化简得到:$\frac{du}{dt}e^{2t} = t^2$两边同时除以$e^{2t}$,得到:$\frac{du}{dt} = t^2e^{-2t}$对上式两边同时积分,得到:$u = -\frac{1}{4}t^2e^{-2t} + C$将$u$代入$y = u(t)e^{2t}$,得到最终的解:$y = (-\frac{1}{4}t^2e^{-2t} + C)e^{2t}$1.2 解:首先,我们根据题意列出方程:$\frac{dy}{dt} = \frac{t}{y}$这是一个一阶可分离变量的常微分方程,我们可以通过分离变量来求解。
将方程变形,得到:$ydy = tdt$对上式两边同时积分,得到:$\frac{1}{2}y^2 = \frac{1}{2}t^2 + C$解得:$y^2 = t^2 + C$由于题目中给出了初始条件$y(0) = 1$,将初始条件代入上式,得到:$1 = 0 + C$解得:$C = 1$将$C$代入$y^2 = t^2 + C$,得到最终的解:$y^2 = t^2 + 1$2. 习题二2.1 解:首先,我们根据题意列出方程:$\frac{dy}{dt} = 2ty + t^2$这是一个一阶线性常微分方程,我们可以使用常数变易法来求解。
偏微分方程总复习和课后习题答案
一、基本概念
1. 偏微分方程的定义P1 2. 偏微分方程的阶数,线性、拟线性、完全非线性 偏微分方程的定义P10 3. 偏微分方程的适定性P23
二、方程的导出,分类与化简
三、公式的直接应用题
1. 2. 3. 4. 5. 达朗贝尔公式P36 公式P42 傅里叶(逆)变换P106 P110例 4.1.7结论 泊松公式P112
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at x a ( t ) 1 t d f ( , )d x a ( t ) 2a 0
1 2 u ( x t ) 3t xt 2
1 1 xa t C f1 ( x at ) ( x at ) ( )d 2 2a x0 2 1 1 xa t C f 2 ( x at ) ( x at ) ( )d 2 2 a x0 2
1 1 xat u [ ( x at ) ( x at )] ( )d 2 2a x a t
1 u ( x t ) x (1 a )t cos x sin at a
2 2 2
1 ( 7)
解:
2
1 22 1 x at x at x u ( x t ) 5 x t a t 2 (e e 2e ) 3 2a
1 ( 6)
解:
2 2u u 2 1 a f ( x , t ), x R ,t 0 2 2 t x u ( x, 0) ( x), u ( x, 0) ( x), x R1. t
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at x a ( t ) 1 t d f ( , )d x a ( t ) 2a 0
第二章习题解答
第二章2-3 设系统传递函数为342)(2++=s s s G 初始条件0/)0(,1)0(=-=dt dc c 。
求单位阶跃输入r (t)=1(t)时,系统的输出响应c (t)。
【解】系统传递函数与微分方程是一一对应的,故通过传递函数先求出微分方程,然后通过拉氏变换的方法求解微分方程。
系统对应的微分方程为 4()3()2()c c t c t r t ++= 在给定的非零初始条件下,进行拉氏变换22(43)()(0)(0)4(0)s s C s sc c c s++---=整理后2221()(43)(43)s C s s s s s s +=-++++部分分式展开后,拉氏反变换111223242/35/25/6()[()][][](43)(43)13255326t t s c t L C s L L s s s s s s s s e e -----+==-=-+++++++=-+2-4 在图2-48中,已知G (s) 和H (s)两方框对应的微分方程分别为()2()5()4()3()6()c t c t e t b t b t c t +=+=图2-48 习题2-4系统结构框图且初始条件为零,试求传递函数C (s)/R (s)。
【解】求出每个方框的传递函数,利用反馈等效的方法求C(s)/R(s)。
根据定义可得 5()2G s s =+,6()43H s s =+ 255()5()25(43)10075(2)56()1()()(2)(43)30411361(2)(43)C s G s s s s R s G s H s s s s s s s +++====+++++++++2-5 图2-49是由电阻、电容和运算法放大器组成的无源网络和有源网络,试列写以V in (t)为输入量,V out (t)为输出量的传递函数。
(a) (b )(c) (d)图2-49 习题2-5电路图【解】(a) 1211211,1RZ R Z C s RC s C s===+ 22112121211()1()11Z C s RC s G s R Z Z R C C s RC s C s +===+++++(b ) 21122211R Z R Z R Cs R Cs ===+ 2222111211()1R Z R Cs R G s Z R R R Cs +=-==-+ (c) 32321123232321()(1)1()1()1R R R R Cs Cs Z R Z R R Cs R R Cs R R Cs++==+==++++ 323232211132(1)()11()()1R R Cs R R Cs R Z R Cs G s Z R R R R Cs ++++=-=-=-++ (d)本题和(b)、(c)做法图通,因为反馈通路有接地的部分。
第二章习题解题过程和参考答案
第二章习题解题过程和参考答案第二章习题解题过程和参考答案2-1 试建立题2-1图所示各系统的微分方程 [其中外力)(t f ,位移)(t x 和电压)(t u r为输入量;位移)(t y 和电压)(t u c为输出量;k (弹性系数),μ(阻尼系数),R (电阻),C (电容)和m (质量)均为常数]。
解:2-1(a) 取质量m 为受力对象,如图,取向下为力和位移的正方向。
作用在质量块m 上的力有外力f(t),重力mg ,这两个力向下,为正。
有弹簧恢复力[]0)(y t y k +和阻尼力()dy t dtμ,这两个力向上,为负。
其中,0y 为0)(=t f 、物体处于静平衡位置时弹簧的预伸长量。
根据牛顿第二定理F ma ∑=,有[]22()()()()dy t d y t f t mg k y t y m dt dtμ+-+-= 其中:0ky mg =代入上式得22)()()()(dt t y d mdt t dy t ky t f =--μ整理成标准式:22()()()()d y t dy t m ky t f t dt dtμ++=μμ()f t[()k y t +()dy t dt或也可写成:22()()1()()d y t dy t k y t f t dt m dt m mμ++=它是一个二阶线性定常微分方程。
2-1(b) 如图,取A 点为辅助质点,设该点位移为()Ax t ,方向如图。
再取B 点也为辅助质点,则该点位移即为输出量()y t ,方向如图A 点力平衡方程:1()()[()()][]AAdx t dy t k x t x t dt dtμ-=- ① B 点力平衡方程:2()()()[]Adx t dy t k y t dt dtμ=- ②由①和②:12[()()]()A k x t x t k y t -= 得:21()()()Akx t x t y t k=-二边微分:21()()()Adx t k dx t dy t dt dt k dt=-③将③代入②:221()()()()[]k dx t dy t dy t k y t dt k dt dtμ=--整理成标准式:1221()()()k k k dy t dx t y t k dt dtμ++=或也可写成:()A t AB1211212()()()()k k k dy t dx t y t dt k k k k dtμ+=++它是一个一阶线性定常微分方程。
常微分方程标准答案-一二章
习题1.24. 给定一阶微分方程2dyx dx=, (1). 求出它的通解; (2). 求通过点()1,4的特解; (3). 求出与直线23y x =+相切的解; (4). 求出满足条件102ydx =⎰的解;(5). 绘出(2),(3),(4)中的解得图形。
解:(1). 通解显然为2,y x c c =+∈;(2). 把1,4x y ==代入2y x c =+得3c =,故通过点()1,4的特解为23y x =+;(3). 因为所求直线与直线23y x =+相切,所以223y x cy x ⎧=+⎨=+⎩只有唯一解,即223x c x +=+只有唯一实根,从而4c =,故与直线23y x =+相切的解是24y x =+;(4). 把2y x c =+代入12ydx =⎰即得5c =,故满足条件12ydx =⎰的解是253y x =+;(5). 图形如下:-1.5-1-0.500.51 1.512345675. 求下列两个微分方程的公共解:242422,2y y x x y x x x y y ''=+-=++--解:由2424222y x x x x x y y +-=++--可得()()222210y x xy -++=所以2y x =或212y x =--,2y x =代入原微分方程满足,而212y x =--代入原微分方程不满足,故所求公共解是代入原微分方程不满足。
6. 求微分方程20y xy y ''+-=的直线积分曲线。
解:设所求直线积分曲线是y kx b =+,则将其代入原微分方程可得2200010k b k xk kx b k b k b k k -=⎧+--=⇒⇒====⎨-=⎩或所以所求直线积分曲线是0y =或1y x =+。
8. 试建立分别具有下列性质的曲线所满足的微分方程:(2). 曲线上任一点的切线介于两坐标轴之间的部分等于定长l ; (5). 曲线上任一点的切线的纵截距等于切点横坐标的平方。
郑大微分方程一、二、三章习题解答
习题解答 习题1.11.一质量为m 的物体,从高度0s 处以初速度0v 铅直向上抛出.设空气的阻力与速度成正比,试求物体的运动规律所满足的微分方程,并写出初始条件.解 如图建立坐标系,设时刻t 时物体的高度为()x t .因物体所受的合力为f mg k dx dt =+,方向向下,由Newton 第二定律可得22d x dx m mg k dt dt=--(0k >为常数), 化简后可得微分方程220d x k dxg dt m dt++=. 若令dx dt v =,则得速度v 满足的一阶微分方程0dv kv g dt m++=, 相应的初始条件为0(0)v v =.2.一高温物体在C20的恒温介质中冷却.设在冷却过程中降温速度与物体和其所在介质的温度差成正比.已知物体的初始温度为0u ,试求物体的温度)(t u 所满足的微分方程,并写出初始条件.解 设时刻t 时物体的温度为)(t u ,由Newton 冷却定律可得微分方程(20)duk u dt=--(0k >为常数), 相应的初始条件为0(0)u u =.3.已知曲线上任一点的切线的纵截距是切点的横坐标和纵坐标的等差中项,试求这曲线所满足的微分方程.解 如图所示建立坐标系,设所求曲线为()y y x =,曲线上的坐标为(,)x y ,过点(,)x y 处的切线上的坐标为(,)X Y ,则切线方程为'()Y y y X x -=-,易见该切线在纵轴上的截距为'b y xy =-.由条件可知'2x yy xy +-=, 整理可得微分方程11'02y y x -+=. 习题1.24.求下列两个微分方程的公共解:24'2y y x x =+-,242'2y x x x y y =++--.解 公共解当然满足关系式2424222y x x x x x y y +-=++--,化简,得22()[2()1]0y x y x -+-=.所以2y x =和212y x =-可能是两个方程的公共解.进一步验证可知前者是公共解,而后者不是.5.求微分方程2''0y xy y +-=的直线解. 解 设直线解为y ax b =+,则2()0a a x ax b +-+=.比较同次幂系数得a b =,2a a =,故0a b ==,或1a b ==.亦即所求的直线解为0y =或1y x =+.6.试求下列曲线族所满足的微分方程:(1)2y cx x =+, (2)12x x y c e c xe =+; (3)平面上的一切圆.解 (1)从2y cx x =+,'2y c x =+消去c 可得微分方程2'0xy x y --=.(2)从12x x y c e c xe =+,12'(1)x x y c e c x e =++,12"(2)x x y c e c x e =++消去12,c c 可得微分方程"2'0y y y -+=.(3)从222()()x a y b c -+-=,()()'0x a y b y -+-=, 21'()"0y y b y ++-=,3'"()"'0y y y b y +-=消去,,a b c 可得微分方程22[1(')]"'3'(")0y y y y +-=.7.给定微分方程22234'x y y xy ,证明其解曲线关于坐标原点(0,0)O 成中心对称的曲线,也是此微分方程的解曲线.证明 设00(,)x y 是方程22234'x y y xy 的积分曲线上任意一点,根据题意,我们只需证明00(,)x y --也是方程22234'x y y xy 的解即可.事实上,设()y y x =为任意积分曲线,00(,)x y 为其上任一点,则2223000004'()()()x y x y x x y x .又设1()y y x =为与积分曲线()y y x =关于坐标原点成中心对称的曲线,则010()()y x y x =--,010'()'()y x y x =-.代入上式,得2223010100104'()()()x y x y x x y x ,即2223010100104()'()()()()x y x y x x y x ,即00(,)x y --也是方程22234'x y y xy 的解.习题1.31. 试用图像法作出如下微分方程的方向场和积分曲线的略图: (1)||'xy xy y; (2)2)1(' y y ; (3)xy y 1'; (4)22'y x y =-. 解 (1)当0xy >时,(,)1f x y º,即在第一、第三象限内任何点的方向斜率均为1;当0xy <时,(,)1f x y ?,即在第二、第四象限内任何点的方向斜率均为1-.由此不难画出方向场及积分曲线的略图.(2)易见2(,)(1)f x y y =-满足解的存在唯一性条件.考察等斜线2(1)y k -=(0k ³),即 1y k =?.当0k =时,1y =(容易验证它是一条积分曲线)。
微分方程习题及答案
微分方程习题§1 基本概念1. 验证下列各题所给出的隐函数是微分方程的解.(1)y x y y x C y xy x -='-=+-2)2(,22(2)⎰'=''=+y 0 222t -)(,1e y y y x dt2..已知曲线族,求它相应的微分方程(其中21C , ,C C 均为常数)(一般方法:对曲线簇方程求导,然后消去常数,方程中常数个数决定求导次数.)(1)1)(22=++y C x ;(2)x C x C y 2cos 2sin 21+=.3.写出下列条件确定的曲线所满足的微分方程。
(1)曲线在()y x , 处切线的斜率等于该点横坐标的平方。
(2)曲线在点P ()y x ,处的法线x 轴的交点为Q,,PQ 为y 轴平分。
(3)曲线上的点P ()y x ,处的切线与y 轴交点为Q , PQ 长度为2,且曲线过点(2,0)。
§2可分离变量与齐次方程1.求下列微分方程的通解(1)2211y y x -='-;(2)0tan sec tan sec 22=⋅+⋅xdy y ydx x ;(3)23xy xy dxdy =-;(4)0)22()22(=++-++dy dx y y x x y x .2.求下列微分方程的特解(1)0 ,02=='=-x y x y e y ;(2)21 ,12==+'=x y y y y x 3. 求下列微分方程的通解(1))1(ln +='xy y y x ; (2)03)(233=-+dy xy dx y x .4. 求下列微分方程的特解(1)1 ,022=-==x y yx xy dx dy ; (2)1 ,02)3(022==+-=x y xydx dy x y .5. 用适当的变换替换化简方程,并求解下列方程(1)2)(y x y +=';(2))ln (ln y x y y y x +=+'(3)11+-='yx y (4)0)1()1(22=++++dy y x xy x dx xy y6. 求一曲线,使其任意一点的切线与过切点平行于y 轴的直线和x 轴所围城三角形面积等于常数2a .7. 设质量为m 的物体自由下落,所受空气阻力与速度成正比,并设开始下落时)0(=t 速度为0,求物体速度v 与时间t 的函数关系.8. 有一种医疗手段,是把示踪染色注射到胰脏里去,以检查其功能.正常胰脏每分钟吸收掉%40染色,现内科医生给某人注射了0.3g 染色,30分钟后剩下0.1g ,试求注射染色后t 分钟时正常胰脏中染色量)(t P 随时间t 变化的规律,此人胰脏是否正常?9.有一容器内有100L 的盐水,其中含盐10kg ,现以每分钟3L 的速度注入清水,同时又以每分钟2L 的速度将冲淡的盐水排出,问一小时后,容器内尚有多少盐?§3 一阶线性方程与贝努利方程1.求下列微分方程的通解(1)2x xy y =-'; (2)0cos 2)1(2=-+'-x xy y x ;(3)0)ln (ln =-+dy y x ydx y ;(4))(ln 2x y y y -='; (5)1sin 4-=-x e dxdy y 2.求下列微分方程的特解(1)0 ,sec tan 0==-'=x yx x y y ;(2)1|,sin 0==+'=x y xx x y y 3.一 曲线过原点,在) ,(y x 处切线斜率为y x +2,求该曲线方程.4.设可导函数)(x ϕ满足方程⎰+=+ x0 1sin )(2cos )(x tdt t x x ϕϕ,求)(x ϕ. 5.设有一个由电阻Ω=10R ,电感H L 2=,电流电压tV E 5sin 20=串联组成之电路,合上开关,求电路中电流i 和时间t 之关系.6.求下列贝努利方程的通解(1) 62y x xy y =+' (2)x y x y y tan cos 4+='(3)0ln 2=-+y x x dydx y (4)2121xy x xyy +-='§4 可降阶的高阶方程1.求下列方程通解。
微分方程课后习题答案
微分方程课后习题答案微分方程是数学中的重要分支,它研究的是描述自然现象中变化规律的方程。
在学习微分方程的过程中,课后习题是巩固知识、提高技能的重要途径。
本文将为大家提供一些微分方程课后习题的答案,希望能够帮助大家更好地理解和掌握微分方程的知识。
1. 一阶线性微分方程题目:求解微分方程 dy/dx + y = 2x解答:这是一个一阶线性微分方程,我们可以使用常数变易法来求解。
首先,将方程改写为 dy/dx = 2x - y设 y = u(x) * v(x),其中 u(x) 是未知函数,v(x) 是待定函数。
将 y = u(x) * v(x) 带入方程,得到 u(x) * v'(x) + u'(x) * v(x) = 2x - u(x) * v(x)整理得 u(x) * v'(x) + u'(x) * v(x) - u(x) * v(x) = 2x根据乘积法则,有 (u(x) * v(x))' = 2x对上式两边同时积分,得到 u(x) * v(x) = x^2 + C,其中 C 是常数。
然后,我们需要求解 u(x) 和 v(x)。
由于 v(x) 是待定函数,我们可以选择 v(x) = e^(-x),这样 v'(x) = -e^(-x)。
将 v(x) = e^(-x) 带入 u(x) * v'(x) + u'(x) * v(x) - u(x) * v(x) = 2x,得到 u'(x) * e^(-x) = 2x对上式两边同时积分,得到 u(x) * e^(-x) = x^2 + C将 u(x) * e^(-x) = x^2 + C 代入 y = u(x) * v(x),得到 y = (x^2 + C) * e^x所以,原微分方程的通解为 y = (x^2 + C) * e^x,其中 C 是常数。
2. 二阶线性常系数齐次微分方程题目:求解微分方程 d^2y/dx^2 + 2dy/dx + 2y = 0解答:这是一个二阶线性常系数齐次微分方程,我们可以使用特征方程法来求解。
常微分方程第四版课后练习题含答案
常微分方程第四版课后练习题含答案第一章:常微分方程基本概念和初值问题1.2 课后练习题1.2.1(1)y′=2y+3,y(0)=1,求解y(t);(2)y′+ty=1,y(0)=0,求解y(t)。
解答:(1)该微分方程为一阶线性常微分方程,其通解为$$y(t)=Ce^{2t}-\\frac{3}{2}$$代入初始条件y(0)=1,可得$$C=\\frac{5}{2}$$所以$$y(t)=\\frac{5}{2}e^{2t}-\\frac{3}{2}$$(2)首先设$u(t)=e^{\\frac{t^2}{2}}y(t)$,则$u'(t)=e^{\\frac{t^2}{2}}(y'+ty)$。
代入原方程可得$$u'(t)=e^{\\frac{t^2}{2}}$$对其积分得$$u(t)=\\int e^{\\frac{t^2}{2}} dt +C=\\frac{\\sqrt{2\\pi}}{2}erf\\frac{t}{\\sqrt{2}}+C$$其中$erf(x)=\\frac{2}{\\sqrt{\\pi}}\\int_0^x e^{-t^2} dt$称为误差函数。
进一步解得$$y(t)=e^{-\\frac{t^2}{2}}u(t)-ue^{-\\frac{t^2}{2}}=-\\frac{\\sqrt{2\\pi}}{2}erf\\frac{t}{\\sqrt{2}}e^{-\\frac{t^2}{2}}$$ 代入初始条件y(0)=0即可得到最终解答。
第二章:一阶线性微分方程2.2 课后练习题2.2.1求下列方程的通解:(1)(2x+1)y′+y=1;(2)(x−1)y′−y=2x;(3)$(2+\\cos x)y'-y=2-x\\cos x$。
解答:(1)该微分方程为一阶线性常微分方程,设方程的通解为$y=Ce^{-\\int \\frac{1}{2x+1} dx}+\\frac{1}{2x+1}$。
常微分方程与动力系统第二章课后题参考答案
常微分方程与动力系统第二章习题参考答案 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曲面的概念1.求正螺面r ={ u v cos ,u v sin , bv }的坐标曲线.解 u-曲线为r ={u 0cos v ,u 0sin v ,bv 0 }={0,0,bv 0}+u {0cos v ,0sin v ,0},为曲线的直母线;v-曲线为r ={0u v cos ,0u v sin ,bv }为圆柱螺线.2.证明双曲抛物面r ={a (u+v ), b (u-v ),2uv }的坐标曲线就是它的直母线。
证 u-曲线为r ={ a (u+0v ), b (u-0v ),2u 0v }={ a 0v , b 0v ,0}+ u{a,b,20v }表示过点{ a 0v , b 0v ,0}以{a,b,20v }为方向向量的直线;v-曲线为r ={a (0u +v ), b (0u -v ),20u v }={a 0u , b 0u ,0}+v{a,-b,20u }表示过点(a 0u , b 0u ,0)以{a,-b,20u }为方向向量的直线。
3.求球面r =}sin ,sin cos ,sin cos {ϑϕϑϕϑa a a 上任意点的切平面和法线方程。
解 ϑr=}cos ,sin sin ,cos sin {ϑϕϑϕϑa a a -- ,ϕr =}0,cos cos ,sin cos {ϕϑϕϑa a -任意点的切平面方程为00cos cos sin cos cos sin sin cos sin sin sin cos cos cos =------ϕϑϕϑϑϕϑϕϑϑϕϑϕϑa a a a a a z a y a x即 xcos ϑcos ϕ + ycos ϑsin ϕ + zsin ϑ - a = 0 ; 法线方程为ϑϑϕϑϕϑϕϑϕϑsin sin sin cos sin cos cos cos cos cos a z a y a x -=-=- 。
常微分方程教程+第二版+丁同仁+李承志+答案和练习第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 轴的方向平移而得到。
马知恩周义仓编常微分方程定性与稳定性方法部分习题参考解答
马知恩周义仓编常微分⽅程定性与稳定性⽅法部分习题参考解答第⼀章 基本定理1设有 $$\bex \frac{\rd \bbx}{\rd t}=\bbf(t,\bbx),\quad \bbx(t_0)=\bbx^0,\quad (t_0,\bbx^0)\in \bbR\times \bbR^n. \eex$$试证: 若 $\bbf\in C^1(G)$, 则在 $(t_0,\bbx^0)$ 的领域内, 此 Cauchy 问题的解存在惟⼀.证明: 由 $f\in C^1(G)$ 蕴含 $f\in C(G)$ 且在 $G$ 内适合 Lipschitz 条件知有结论.2试讨论下列⽅程解的存在区间:(1) $\dps{\frac{\rd y}{\rd x}=\frac{1}{x^2+y^2}}$;(2) $\dps{\frac{\rd y}{\rd x}=y(y-1)}$.解答:(1) 由 $\dps{\frac{\rd x}{\rd y}=x^2+y^2}$ 的解的存在区间有限知 $y$ 有界, ⽽由解的延拓定理, 原⽅程解的存在区间为 $\bbR$.(2) 直接求解有 $\dps{y=\frac{1}{1-\frac{y_0-1}{y_0}e^x}}$, ⽽a.当 $0\leq y_0\leq 1$ 时, 原⽅程解的存在区间为 $\bbR$;b.当 $y_0<0$ 时, 原⽅程解的存在区间为 $\dps{\sex{\ln\frac{y_0}{y_0-1},\infty}}$;c.当 $y_0>1$ 时, 原⽅程解的存在区间为 $\dps{\sex{-\infty,\ln\frac{y_0}{y_0-1}}}$.3 设有⼀阶微分⽅程式 $$\bex \frac{\rd x}{\rd t}=(t-x)e^{tx^2}. \eex$$ 试证: 过任⼀点 $(t_0,x_0)\in\bbR^2$ 的右⾏解的存在区间均为 $[t_0,+\infty)$.证明: 由 $$\bex \frac{\rd x}{\rd t}=(t-x)e^{tx^2}=\left\{\ba{ll} <0,&x>t,\\ >0,&x<t \ea\right. \eex$$ 知解在 $\sed{x>t}$ 内递减,在 $\sed{x<t}$ 内递增. 当 $x_0>t_0$ 时, 在 $$\bex \sed{(t,x);t\in\bbR, t_0<x<x_0} \eex$$ 内应⽤解的延伸定理知解定与$\sed{x=t}$ 相交, 之后解递增, 在 $$\bex \sed{(t,x);t\in\bbR,x<t} \eex$$ 内应⽤延伸定理及⽐较定理即知结论.4设有⼀阶⽅程 $\dps{\frac{\rd x}{\rd t}=f(x)}$, 若 $f\in C(-\infty,+\infty)$, 且当 $x\neq 0$ 时有 $xf(x)<0$. 求证过 $\forall\(t_0,x_0)\in\bbR^2$, Cauchy 问题的右⾏解均在 $[t_0,+\infty)$ 上存在, 且 $\dps{\lim_{t\to+\infty}x(t)=0}$.证明: 由题意, $$\bex f(x)\left\{\ba{ll} >0,&x<0,\\ <0,&x>0. \ea\right. \eex$$ ⽽由 $f$ 的连续性, $f(0)=0$. 于是当 $x_0=0$ 时,由解的唯⼀性知 $x=0$. 当 $x_0>0$ 时, 在 $$\bex \sed{(t,x);t\in\bbR,0<x<x_0} \eex$$ 内应⽤延伸定理及惟⼀性定理知 $x(t)$ 递减趋于 $0$. 当 $x_0<0$ 时, 在 $$\bex \sed{(t,x);t\in\bbR,x_0<x<0} \eex$$ 内应⽤延伸定理及惟⼀性定理知 $x(t)$ 递增趋于 $0$.5若 $\bbf(t,\bbx)$ 在全空间 $\bbR\times\bbR^n$ 上连续且对 $\bbx$ 满⾜局部 Lipschitz 条件且 $$\bex \sen{\bbf(t,\bbx)}\leq L(r),\quad r=\sqrt{\sum_{i=1}^n x_i^2},\quad \bbx=(x_1,\cdots,x_n)^T, \eex$$ 其中 $L(r)>0, r>0$, 且 $$\bee\label{1.5:1}\int_a^{+\infty}\frac{\rd r}{L(r)}=+\infty,\quad a>0. \eee$$ 试证: 对 $\forall\ (t_0,\bbx^0)\in\bbR\times\bbR^n$, Cauchy 问题的解均可对 $t$ ⽆限延拓.证明: 由解的延伸定理, 仅须证明在任何有限区间 $-\infty<\alpha<t<\beta<+\infty$ 上, $\bbx(t)$ 有界. 为此, 令 $y(t)=\sen{\bbx(t)}$,则 $$\beex \bea \frac{\rd y(t)}{\rd t}&=2\bbx(t)\cdot\frac{\rd \bbx(t)}{\rd t} =2\bbx(t)\cdot \bbf(t,\bbx(t)),\\\sev{\frac{\rd y(t)}{\rd t}} &\leq 2\sqrt{y(t)}\cdot L\sex{\sqrt{y(t)}},\\ \frac{\rd \sqrt{y(t)}}{L\sex{\sqrt{y(t)}}}&\leq \rd t,\\ \int_\alpha^\beta \frac{\rd \sqrt{y(t)}}{L\sex{\sqrt{y(t)}}} &\leq \int_\alpha^\beta \rd t=\beta-\alpha. \eea \eeex$$ 这与\eqref{1.5:1} ⽭盾 (事实上, 当 $\alpha,\beta\gg 1$, $|\alpha-\beta|\ll 1$ 时, 不等式右端可任意⼩, ⽽不等式左端有积分发散知可⼤于某⼀正常数).6设有微分⽅程 $$\bex \frac{\rd \bbx}{\rd t}=\bbf(t,\bbx), \eex$$ $\bbf\in C(G\subset \bbR\times\bbR^n)$, 试证: 若对$\forall\ (t_0,\bbx^0)\in G$, Cauchy 问题的解都存在唯⼀, 则解必对初值连续依赖.证明: 参考[家⾥蹲⼤学数学杂志第134期, 常微分⽅程习题集, 第1600页].7 试在定理 1.1 的假设下, 利⽤ Gronwall 引理直接证明解对初始时刻 $t_0$ 的连续依赖性.证明: 参考定理 1.7 的证明.8 设有⼀阶 Cauchy 问题 $$\bex \frac{\rd y}{\rd x}=x^2+(y+1)^2,\quad y(0)=0. \eex$$ 试利⽤⽐较定理证明, 若设解的右⾏饱和区间为 $[0,\beta)$, 则 $\dps{\frac{\pi}{4}\leq \beta\leq 1}$.证明: 仅须注意到当 $0\leq x\leq 1$ 时, $$\bex (y+1)^2\leq x^2+(y+1)^2\leq 1+(y+1)^2. \eex$$ 再利⽤⽐较定理即知结论.第⼆章 动⼒系统的基本知识1试证明: $\Omega_P=\vno$ 的充要条件是 $L_P^+$ 趋于⽆穷.证明: $\ra$ ⽤反证法. 若 $L_P^+$ 不趋于⽆穷, 则 $$\bex \exists\ M>0, t_n\nearrow +\infty,\st \sen{\mbox{ $\varphi$}(P,t_n)}\leq M. \eex$$ 由 Weierstrass 定理, $$\bex \exists\ \sed{t_n'}\subset \sed{t_n},\st \mbox{ $\varphi$}(P,t_n)\to Q,\eex$$ ⽽ $Q\in \Omega_P$, 这是⼀个⽭盾. $\la$ 亦⽤反证法. 若 $\Omega_P\neq \vno$, ⽽设 $Q\in \Omega_P$, 则 $$\bex\exists\ t_n\nearrow+\infty,\st \mbox{ $\varphi$}(P,t_n)\to Q. \eex$$ 这与 $L_P^+$ 趋于⽆穷⽭盾.2试证明: 若 $\Omega_P$ 仅含惟⼀奇点 $P^*$, 则当 $t\to+\infty$ 时必有 $L_P^+$ 趋向于 $P^*$.证明: ⽤反证法. 设 $$\bee\label{2.2:1} \exists\ \ve_0>0,\ t_n\nearrow+\infty, \st \sen{\mbox{ $\varphi$}(P,t_n)-P^*}\geq\ve_0. \eee$$ 则(1)若 $\sed{t_n}$ 有有界的⼦列, 则适当抽取⼦列 $\sed{t_n'}$ 后有 $$\bex \mbox{ $\varphi$}(P,t_n')\to Q. \eex$$ 于是 $Q\in\Omega_P=\sed{P^*}$. 这与 \eqref{2.2:1} ⽭盾.(2)若 $\sed{t_n}$ ⽆有界的⼦列, 则 $\dps{\lim_{n\to\infty}\mbox{ $\varphi$}(P,t_n)=\infty}$, ⽽ $\infty\in\Omega_P=\sed{P^*}$, ⼜是⼀个⽭盾.3试证明: 若 $\Omega_P$ 有界且 $\Omega_P$ ⾮闭轨, 则 $\forall\ R\in \Omega_P$, $\Omega_R$ 与 $A_R$ 必均为奇点.证明: ⽤反证法证明 $\Omega_R$ 为奇点集, $A_R$ 为奇点集类似可证. 设 $\Omega_R$ 含有常点. 由 $R\in \Omega_P$ 及$\Omega_P$ 为不变集知 $L_R\subset \Omega_Q$. 于是按引理 2.3, $L_R$ 为闭轨线, $L_R=\Omega_R\subset \Omega_P$. 这与 $\Omega_P$ ⾮闭轨⽭盾.4试证明: ⼀系统的圈闭奇点的集合是⼀闭集.证明: 全体奇点的集合为 $$\bex \sed{\bbx^*\in G; \bbf(\bbx^*)=\mbox{ $0$}}. \eex$$ 由 $\bbf$ 的连续性即知结论.5 若 $L_P^+$ 有界且 $\Omega_P$ 仅由奇点构成, 能否断定 $\Omega_P$ 仅含⼀个奇点?解答: 不能断定. 仅能说 $\Omega_P$ 为由奇点构成的连通闭集或闭轨线.6 设 $O(0,0)$ 是⼀平⾯⾃治系统的惟⼀奇点, 且是稳定的, 全平⾯没有闭轨线. 试证: (1) 此系统的任⼀轨线必负向⽆界; (2) 任⼀有界的正半轨闭进⼊奇点 $O$.证明:(1) ⽤反证法. 若有⼀轨线负向有界, 则在定理 2.8 中, 由全平⾯没有闭轨线知 (3),(4) 不成⽴; 由 $O$ 为惟⼀奇点知 (1),(2),(5) 不成⽴. 这是⼀个⽭盾.(2) 对有界正半轨⽽⾔, 定理 2.8 中仅有 (1),(2),(5) 可能成⽴. 若 (1),(2) 成⽴, 则结论已证; ⽽由全平⾯没有闭轨线知 (5) 不成⽴.第三章 稳定性理论1 讨论⽅程 $$\bee\label{3.1:1} \sedd{\ba{ll}\frac{\rd x_1}{\rd t}=x_2,\\ \frac{\rd x_2}{\rd t}=-a^2\sin x_1\ea} \eee$$ 零解的稳定性.解答: 选取 $$\bex V(\bbx)=\frac{x_2^2}{2}+a^2(1-\cos x_1), \eex$$ 则 $V$ 在原点的⼀邻域内是正定的, 且沿 \eqref{3.1:1} 的轨线有 $$\bex \dot V(\bbx)=V_{x_1}x_1'+V_{x_2}x_2'=0. \eex$$ 由此, 零解是稳定的, 但不是渐近稳定的.2 证明⽅程 $\dps{\frac{\rd x}{\rd t}=-x+x^2}$ 的零解是指数渐近稳定的, 但不是全局渐近稳定的.证明: 解该微分⽅程有: $$\bex \ba{ccc} -\frac{1}{x^2}\frac{\rd x}{\rd t}=\frac{1}{x}-1,&\frac{\rd y}{\rd t}=y-1\\sex{y=\frac{1}{x}},&\frac{\rd z}{\rd t}=-e^{-t}\ \sex{z=e^{-t}y},\\ z=e^{-t}+C,&y=Ce^t+1,&x=\frac{1}{1+Ce^t}. \ea \eex$$由此, 原微分⽅程的解为 $$\bex x=0,\mbox{ 或 }x(t)=\frac{1}{1+Ce^t}. \eex$$ 取初值 $(t_0,x_0),\ x_0\neq 0$, 有 $$\bexx(t,t_0,x_0)=\frac{x_0}{1+e^{t-t_0}(1-x_0)}. \eex$$ 故当 $|x_0|<1$ 时, $$\bex |x(t,t_0,x_0)|\leq \sev{\frac{1}{x_0}-1}e^{-(t-t_0)}. \eex$$ 这说明零解是指数渐近稳定的. 但由于从 $(t_0,1)$ 出发的解 $x(t,t_0,1)=1$ 不趋于零解, ⽽零解不是全局渐近稳定的.3 在相空间 $\bbR^n$ 中给出 $\dps{\frac{\rd \bbx}{\rd t}=\bbf(t,\bbx),\ \bbf(t,0)=0}$ 的零解稳定、渐近稳定、不稳定的⼏何解释.解答: 零解是稳定的 $\lra\ \forall\ \ve>0,\ \exists\ \delta>0,\ \forall\ P\in B_\delta,\ L_P^+\subset B_\ve$; 零解是渐进稳定的$\lra\ \exists\ U\ni O,\ \forall\ P\in U,\ L_P^+\to 0$; 零解是不稳定的 $\lra\ \exists\ \ve_0>0,\ \exists\ P_n\to0, \stL_{P_n}^+\bs B_\ve\neq \vno$.4判断下列系统零解的稳定性:(1) $\dps{\sedd{\ba{ll} \frac{\rd x_1}{\rd t}=mx_2+\alpha x_1(x_1^2+x_2^2),\\ \frac{\rd x_2}{\rd t}=-mx_1+\alphax_2(x_1^2+x_2^2); \ea}}$;(2) $\dps{\frac{\rd^2x}{\rd t^2}+\sex{\frac{\rd x}{\rd t}}^3+f(x)=0,}$ 其中 $xf(x)>0\ (x\neq 0), f(0)=0$;(3) $\dps{\frac{\rd^2x}{\rd t^2}-\sex{\frac{\rd x}{\rd t}}^2sgn\sex{\frac{\rd x}{\rd t}}+x=0}$.解答:(1) 取 $$\bex V=x_1^2+x_2^2, \eex$$ 则 $V$ 正定, 且沿微分⽅程的轨线有 $$\bex \dot V=2\alpha(x_1^2+x_2^2)\sedd{\ba{lll} \mbox{正定},&\alpha>0,\\ 0,&\alpha=0,\\ \mbox{负定},&\alpha<0. \ea} \eex$$ 于是当 $\alpha>0$ 时, 由定理 3.3, 零解是不稳定的; 当 $\alpha=0$ 时, 由定理 3.1, 定理是稳定的; 当 $\alpha<0$ 时, 由定理 3.1, 零解是渐近稳定的.(2) 令 $\dps{x_1=x,x_2=\frac{\rd x}{\rd t}}$, 则 $$\bex \frac{\rd x_1}{\rd t}=x_2,\quad \frac{\rd x_2}{\rd t}=-x_2^3-f(x_1). \eex$$ 取 $$\bex V=\frac{x_2^2}{2}+\int_0^{x_1}f(t)\rd t, \eex$$ 则 $V$ 正定, 且沿微分⽅程的轨线有 $\dot V=-x_2^4\leq 0.$再 $$\bex \sed{\bbx;\dot V(\bbx)=0}=\sed{0}, \eex$$ 我们据定理 3.2 知零解是渐近稳定的.(3) 令 $\dps{x_1=x,x_2=\frac{\rd x}{\rd t}}$, 则 $$\bex \frac{\rd x_1}{\rd t}=x_2,\quad \frac{\rd x_2}{\rd t}=x_2^2sgn(x_2)-x_1. \eex$$ 取 $$\bex V=\frac{x_1^2+x_2^2}{2}, \eex$$ 则 $V$ 正定, 且沿微分⽅程的轨线有 $\dot V=x_2^2|x_2|$是正定的. 我们据定理 3.3 知零解是不稳定的.5 若存在有⽆穷⼩上界的正定函数 $V(t,\bbx)$, 它沿着 $$\bex (3.3.1)\quad \frac{\rd\bbx}{\rd t}=\bbf(t,\bbx),\quad \bbf(t,0)=0 \eex$$ 解曲线的全导数 $\dot V(t,\bbx)$ 负定, 证明 (3.3.1) 的零解是渐近稳定的.证明: 仅须注意到存在正定函数 $W(x)$, $W_1(x)$ 使得 $$\bex W(\bbx)\leq V(t,\bbx)\leq W_1(\bbx). \eex$$ ⽽可仿照定理 3.1 的证明.6 讨论 $\dps{\frac{\rd x}{\rd t}=\frac{g'(t)}{g(t)}x}$ 零解的稳定性, 其中 $\dps{g(t)=\sum_{n=1}^\infty \frac{1}{1+n^4(t-n)^2}}$. 能否得到零解渐近稳定的结果? 为什么?解答: 直接求解有 $$\bex x(t)=\frac{x_0}{g(t_0)}{g(t)}, \eex$$ ⽽由 $$\bex |x(t)|\leq\frac{|x_0|}{g(t_0)}\sez{2+\sum_{n\neq [t],[t]+1}\frac{1}{1+n^4(t-n)^2}} \leq \frac{|x_0|}{g(t_0)}\sez{2+\sum_{n=1}^\infty\frac{1}{n^4}} \eex$$ 知零解是稳定的; 由$$\bex |x(k)|=\frac{|x_0|}{g(t_0)}\sez{1+\sum_{n\neq k}\frac{1}{n^4(k-n)^2}}\geq \frac{|x_0|}{g(t_0)} \eex$$ 知零解不是渐近稳定的.7证明 $\dps{\frac{\rd x}{\rd t}=-\frac{x}{t+1}}$ 的零解是渐近稳定的, 但不存在有⽆穷⼩上界的正定函数 $V(t,x)$, 使得 $\dotV(t,x)$ 负定 (该习题表明习题 5 中渐近稳定性定理中的条件不是必要的).证明: 直接求解有 $$\bex x(t)=\frac{x_0}{1+t}. \eex$$ ⽽零解是渐近稳定的.。
常微分方程2.2习题参考解答
习题2.2求下列方程的解1.dxdy =x y sin +解:y=e ⎰dx (⎰x sin e ⎰-dx c dx +)=e x [-21e x -(x x cos sin +)+c]=c e x -21(x x cos sin +)是原方程的解。
2.dt dx +3x=e t2解:原方程可化为:dt dx =-3x+e t 2所以:x=e ⎰-dt 3(⎰e t 2e -⎰-dt 3c dt +)=e t 3-(51e t 5+c)=c e t 3-+51e t 2是原方程的解。
3.dt ds =-s t cos +21t 2sin 解:s=e ⎰-tdt cos (t 2sin 21⎰e dt dt ⎰3c +)=e t sin -(⎰+c dt te t t sin cos sin )=e t sin -(c e te t t +-sin sin sin )=1sin sin -+-t ce t 是原方程的解。
4.dx dy nx x e y n x =-,n 为常数.解:原方程可化为:dx dy nx x e y nx +=)(c dx e x e e y dx x n n x x n +⎰⎰=⎰-)(c e x x n +=是原方程的解.5.dx dy +1212--y x x =0解:原方程可化为:dx dy =-1212+-y x x ⎰=-dx x x e y 212(c dx e x x+⎰-221))21(ln 2+=x e )(1ln 2⎰+--c dx ex x =)1(12x ce x +是原方程的解.6.dx dy 234xy x x +=解:dx dy 234xy x x +==23y x +x y 令x y u =则uxy =dx dy =u dx du x +因此:dx du x u +=2u x 21udx du =dxdu u =2c x u +=331c x x u +=-33(*)将x y u =带入(*)中得:3433cx x y =-是原方程的解.3332()21()227.(1)12(1)12(),()(1)1(1)(())1(1)dx P x dx x P x dx dy y x dx x dy y x dx x P x Q x x x e e x e Q x dx c x +--=++=+++==++⎰⎰==+⎰⎰++⎰⎰P(x)dx 232解:方程的通解为: y=e =(x+1)(*(x+1)dx+c) =(x+1)((x+23221(1)()211,()(())dy y x c dy y dx x y dx x y dy y yQ y y ye yQ y dy c -+++==+=⎰⎰==⎰⎰+⎰⎰2243P(y)dy P(y)dy P(y)dy 1)dx+c)=(x+1) 即:2y=c(x+1)+(x+1)为方程的通解。
常微分方程答案 一二章
习题1.24. 给定一阶微分方程2dyx dx=, (1). 求出它的通解; (2). 求通过点()1,4的特解;(3). 求出与直线23y x =+相切的解; (4). 求出满足条件102ydx =⎰的解;(5). 绘出(2),(3),(4)中的解得图形。
解:(1). 通解显然为2,y x c c =+∈ ;(2). 把1,4x y ==代入2y x c =+得3c =,故通过点()1,4的特解为23y x =+;(3). 因为所求直线与直线23y x =+相切,所以223y x cy x ⎧=+⎨=+⎩只有唯一解,即223x c x +=+只有唯一实根,从而4c =,故与直线23y x =+相切的解是24y x =+;(4). 把2y x c =+代入12y d x =⎰即得53c =,故满足条件12ydx =⎰的解是25y x =+; (5). 图形如下:-1.5-1-0.500.51 1.512345675. 求下列两个微分方程的公共解:242422,2y y x x y x x x y y ''=+-=++--解:由2424222y x x x x x y y +-=++--可得()()222210y x xy -++=所以2y x =或212y x =--,2y x =代入原微分方程满足,而212y x =--代入原微分方程不满足,故所求公共解是代入原微分方程不满足。
6. 求微分方程20y xy y ''+-=的直线积分曲线。
解:设所求直线积分曲线是y kx b =+,则将其代入原微分方程可得2200010k b k xk kx b k b k b k k -=⎧+--=⇒⇒====⎨-=⎩或所以所求直线积分曲线是0y =或1y x =+。
8. 试建立分别具有下列性质的曲线所满足的微分方程:(2). 曲线上任一点的切线介于两坐标轴之间的部分等于定长l ; (5). 曲线上任一点的切线的纵截距等于切点横坐标的平方。
数理方程第一章、第二章习题全解
u( 0 , t) = u( l, t) = 0 现考虑初始条件,当冲量 k 作用于 x = c处时, 就相当于在这点 给出了一个初速度 , 我们考虑以 c点为中心 , 长为 2δ的一小段弦 ( c δ, c + δ) , 设弦是均匀的 , 其线密度为 ρ, 则这 一小段 弦的质量 为 2δρ, 受冲击时速度为 ut ( x, 0) , 由动量定理得
h c
x
l
h -
c(
l
-
x)
(0 ≤ x ≤ c) ( c < x ≤ l)
ut ( x, 0) = ψ( x ) = 0
则 u( x, t) 是下列定解问题的解 :
utt - a2 uxx = 0
( 0 < x < l, t > 0)
u( x, 0) = φ( x ) , ut ( x, 0 ) = ψ( x )
2 .4 习题全解
1. 设弦的两端固定于 x = 0 及 x = l, 弦的初始位称如图 2 2 所 示,初速度为零, 又设有外力作用, 求弦作横向振动时的位移函数 u( x, t) 。
解 如图 2 2 所示, 弦作横向振动时初始条件为
62
数学物理方程与特殊函数导教·导学·导考
图2 2
u( x, 0) = φ( x ) =
5. 若 F( z) , G( z) 是任意两个二次连续可微函数 , 验证
u = F( x + at ) + G( x - at )
满足方程
2u t2
=
a2
2x2u。
解 作自变量代换ξ= x + at,η= x - at, 由复合函数求导法则
有
所以 于是
u t
微分方程数值解法(李荣华3版)第二章习题答案(大)
第二章习题课(2007.4.28)习题1.求两点边值问题22sin , 0142(0)0, (1)0xLu u u x u u ππ⎧''=-+=<<⎪⎨⎪'==⎩(1.1)的线性有限元解函数(区间等距剖分成2段或3段),要求在计算总刚度矩阵和总荷载向量时,所涉及的定积分用两种方法: 1. 精确求解;2. 用中矩形公式近似计算。
解:第一步:写出原问题(1.1)的等价变分形式(基于虚功原理)试探函数空间和检验函数空间均为:11(){ |(), ()0 }E H I u u H I u a =∈=.在(1.1)的第一个式子两边同时乘以检验函数空间1()E H I 中的任意元素v ,再在区间(0,1)I =上积分,可得21112sin42xu vdx uvdx vdx ππ''-+=⎰⎰⎰ (1.2)其中111011[(1)(1)(0)(0)]u vdxu v dx vu u v dx v u v u u v dx'''''-=-''''=--''=⎰⎰⎰⎰分部积分(1.3)将(1.3)代入(1.2),可得211()2sin42xu v uv dx vdx ππ''+=⎰⎰记21010(,)()4()2sin 2a u v u v uv dx x f v vdxππ⎧''=+⎪⎪⎨⎪=⎪⎩⎰⎰ 则可以得到原问题(1.1)的等价变分问题:求1()E u H I ∈,使得1(,)(), ()Ea u v f v v H I =∀∈. (1.4)第二步:线性有限元空间的构造1.网格剖分(这里以等距剖分3段为例)2.一次Lagrange 有限元空间的定义1{ ():|(),1,2,3, (0)0 }E i h h h e i h V u C I u P e i u =∈∈==.3. Lagrange 节点基函数的构造113, [0,]312()23, [,]330,x x x x x φ⎧∈⎪⎪⎪=-∈⎨⎪⎪⎪⎩在别处 ; 21231, [,]332()33, [,1]30,x x x x x φ⎧-∈⎪⎪⎪=-∈⎨⎪⎪⎪⎩在别处; 3232, [,1]()30,x x x φ⎧-∈⎪=⎨⎪⎩ 在别处.4.空间E hV 中元素的(整体)表示记 (), 1,2,3i h i u u x i ==,则对E hh u V ∀∈,有31()()h j j j u x u x φ==∑ (1.5)第三步:写出线性有限元方程将原变分问题(1.4)中1()EHI 的试探函数子空间和检验函数子空间均取为E h V ,则可以得到原问题(1.1)的近似变分问题:求 E hhu V ∈,使得 (,)(), E h h h h h a u v f v v V =∀∈. (1.6)利用(1.5)并将 h v 取为(), 1,2,3i x i φ=则上述近似变分问题等价于求123,,u u u R ∈,使得31(,)(), 1,2,3j j i i j a u f i φφφ===∑⇔ 31(,)(), 1,2,3j i j i j a u f i φφφ===∑⇔ 31(,)(), 1,2,3i j j i j a u f i φφφ===∑ 写成矩阵形式AU b =其中111213212223313233(,)(,)(,)(,)(,)(,)(,)(,)(,)a a a A a a a a a a φφφφφφφφφφφφφφφφφφ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,123u U u u ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦, 123()()()f b f f φφφ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦其中(a ) 精确求解以11(,)a φφ和1()f φ的计算为例:212211110122222223311111031222222233103(,)[()]4[()][()]44[3(3)][(3)(23)]44a dxdx dxx dx x dx πφφφφππφφφφππ'=+''=+++=++-+-=⎰⎰⎰⎰⎰1221(,)(,)a a φφφφ==,1331(,)(,)a a φφφφ==,22(,)a φφ=2332(,)(,)a a φφφφ==,33(,)a φφ=11101233103()2sin2 2sin (3)2sin (23)22xf dxx x x dx x dx πφφππ==+-=⎰⎰⎰(b )中矩形公式近似求解中矩形公式:()()()2baa bg x dx b a g +≈-⎰.以11(,)a φφ和1()f φ的计算为例:222222112221111(,)[3(3)][(3)(23)]34634211 (9)(9)3163162 (9)316a ππφφπππ≈++-+-=+++=+ 111111162()2sin (3)2sin (23)32632222 sin sin32438f ππφππ≈+-=+习题2.导出下面边值问题1122(), ()(), ()()d du Lu p qu f a x bdx dx u a u a u b u b αβαβ⎧=-+=<<⎪⎨⎪''+=+=⎩ (2.1)的线性有限元方程。
第2章 平稳过程习题答案
第二章 平稳过程1.指出下面所给的习题中,哪些是平稳过程,哪些不是平稳过程?(1)设随机过程Xt e t X -=)(,t >0,其中X 具有在区间),0(T 中的均匀分布 解:∵ 该随机过程的数学期望为∴ 该随机过程不是平稳过程。
(2)设随机过程}),({+∞<<-∞t t X 在每一时刻的状态只取0或1数值,而在不同时刻的状态是相互独立的,且对任意固定的t 有1}0)({ }1)({p t X P p t X P -====其中10<<p解:∵ 该随机过程的数学期望为p t X P t X P t EX t m X ==⋅+=⋅==}0)({0}1)({1)()((常数)该随机过程的自相关函数为:}0)()({0}1)()({1)]()([),(=+⋅+=+⋅=+=+ττττt X t X P t X t X P t X t X E t t R X2}1)({}1)({p t X P t X P ==+==τ 结果与t 无关∴ 该随机过程是平稳随机过程。
(3)设}1,{≥n X n 是独立同分布的随机序列,其中j X 的分布列为定义 ∑==nj jn XY 1,试对随机序列}1,{≥n Y n ,讨论其平稳性。
解:∵ 0211211}1{)1(}1{1=⋅-⋅=-=-+=⋅=j j j X P X P EX ∴ ∑∑=====nj jnj j n EXXE EY 110)((常数)又因为随机序列n Y 的自相关函数。
⎥⎥⎦⎤⎢⎢⎣⎡=+=+∑∑=+=n j m n k k j Y X X E m n Y n EY m n n R 11)()(),( m 为自然数⎰≠--=-===---TTtT xt xtx const e Tte Tt dx T e t EX t m 00]1[111)()(⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛+=∑∑∑==+=n j n k m n k k j j X X X E 111 ⎥⎥⎦⎤⎢⎢⎣⎡+⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛=∑∑∑∑∑∑=+===+==n j m n k k j n j j nj mn k k j nj j X X E X E X X X E 11211121 ∑∑=+=⋅+=nj mn k k jn EX EXEY 112n n n n DY EY DY EY =+==22)(∵ []∑∑∑∑======-==⎪⎪⎭⎫ ⎝⎛=nj n j n j j j j j n j j n np EX EX EX DX X D DY 1112221)( 即)(),(m R np m n n R Y Y ≠=+ ∴ 该随机过程不是平稳过程。
常微分方程课后习题部分答案 (2)
18. 设),(y x f 及连续,试证方程0),(=-dx y x f dy 为线性方程的充要条件是它有仅依赖于x 的积分因子.证:必要性 若该方程为线性方程,则有)()(x Q y x P dx dy += , 此方程有积分因子⎰=-dx x P e x )()(μ,)(x μ只与x 有关 .充分性 若该方程有只与x 有关的积分因子)(x μ .则0),()()(=-dx y x f x dy x μμ为恰当方程 , 从而dxx d y y x f x )()),()((μμ=∂-∂ ,)()(x x y f μμ'-=∂∂ , )()()()()()()()(x Q y x P x Q y x x x Q dy x x f +=+'-=+'-=⎰μμμμ . 其中)()()(x x x P μμ'-= .于是方程可化为0))()((=+-dx x Q y x P dy 即方程为一阶线性方程.20.设函数f(u),g(u)连续、可微且f(u)≠g(u),\,试证方程yf(xy)dx+xg(xy)dy=0 有积分因子u=(xy[f(xy)-g(xy)])1-证:在方程yf(xy)dx+xg(xy)dy=0两边同乘以u 得:uyf(xy)dx+uxg(xy)dy=0 则y uyf ∂∂=uf+uy y f ∂∂+yf y u ∂∂=)(g f xy f -+)(g f xy y f y-∂∂-yf 222)()(g f y x y g xy y f xy g f x -∂∂+∂∂+- =2)(g f xy y f gy y g yf-∂∂-∂∂=2)(g f x y xy xy f g y xy xy g f -∂∂∂∂-∂∂∂∂ =2)(g f xy f g xy g f-∂∂-∂∂ 而x uxg ∂∂=ug+ux x g ∂∂+xg x u ∂∂=)(g f xy g -+)(g f xy x g x -∂∂- xg 222)()(g f y x x g xy x f xy g f y -∂∂-∂∂+-=2)(g f xy x xy xy f xg x xy xy g xf-∂∂∂∂-∂∂∂∂=2)(g f xy f g xy g f -∂∂-∂∂ 故y uyf ∂∂=xuxg ∂∂,所以u 是方程得一个积分因子 21.假设方程(2.43)中得函数M (x,y )N(x,y)满足关系xN y M ∂∂-∂∂= Nf(x)-Mg(y),其中f(x),g(y)分别为x 和y 得连续函数,试证方程(2.43)有积分因子u=exp(⎰dx x f )(+⎰dy y g )()证明:M(x,y)dx+N(x,y)dy=0 即证x uN y uM ∂∂=∂∂)()(⇔u y M ∂∂+M y u ∂∂=u xN ∂∂+N x u ∂∂⇔ u(y M ∂∂-x N ∂∂)=N xu ∂∂- M y u ∂∂⇔u(y M ∂∂-x N ∂∂)=Ne ⎰⎰+dy y g dx x f )()(f(x) -M e ⎰⎰+dy y g dx x f )()(g(y)⇔u(y M ∂∂-x N ∂∂)=e ⎰⎰+dy y g dx x f )()((Nf(x)-Mg(y)) 由已知条件上式恒成立,故原命题得证。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
∂( μ(xy)M ) = ∂( μ(xy)N )
∂y
∂x
即
μ(xy)(∂M − ∂N ) = N ∂μ(xy) − M ∂μ(xy)
∂y ∂x
∂x
∂y
µ(xy)(∂M − ∂N ) = ( yN − xM ) dµ(xy) ,
∂y ∂x
d (xy)
∂M ∂N −
∂y
∂x
dµ ( xy)
=
⋅
1
= g(xy) ,
µ(x, y) =
1
。
xM (x, y) + yN (x, y)
方法 3 用定义求积分因子。
由积分因子的定义,只需证明二元函数 µ(x, y) =
1
满足
xM (x, y) + yN (x, y)
∂(µM ) ∂(µN )
=
即可。为此,我们计算
∂y
∂x
∂( M )
∂(µM ) xM + yN
=
∂y
∂y
仅依赖于 x 的积分因子。 证 必要性。若方程 dy − f (x, y)dx = 0 为线性方程,则方程可写为
dy − (P(x) y + Q(x))dx = 0,令
M = −(P(x) y + Q(x)) , N = 1 ,
∂M ∂N
−
∂M
∂y
由题有 连续,
∂x = −P(x) ,
∂y
N
由定理 2-2 的结论 1 方程有积分因子 e∫ −P( x) dx ,仅依赖于 x 。
x m{[M (1,u) + N (1, u)u]dx + xN (1,u)du} = 0 ,
可以看出上方程为可分离变量的方程,只要给上式乘以积分因子
µ(x, y)
=
1 x m+1[M (1, u) + uN (1,u)]
=
1 xM (x, y) +
yN(x, y)
,
方程就可变量分离,即化为恰当方程,因此,齐次方程的积分因子是
yN − xM d(xy) µ(xy)
因此方程具有形如 µ( xy) 的积分因子的充要条件是
∂M ∂N −
∂y ∂x = g (xy) 。 yN − xM
评注:利用对称形式的微分方程的系数容易判断方程是否具有特殊形式的积分因子,从 而给出求积分因子的思路。
2-6 设 f (x, y) 及 ∂f 连续,试证方程 dy − f (x, y)dx = 0 为线性方程的充要条件是它有 ∂y
分离变量,可求得
f (x) = ± 1 , 2(x + C)
1
代入原方程可得 C = 0 ,从而 f (x) 的一般表达式为 f (x) =
。
2x
评注:本题中常数的确定不能直接通过所给积分方程得到,而是需将通解代回原方程来 确定。
2-2 求具有性质 x(t + s) = x(t) + x(s) 的函数 x(t) ,已知 x′(0) 存在。 1 − x(t)x(s)
有积分因子 µ = (xy[ f (xy) − g(xy)])−1 。
证 方法 1 用积分因子定义证明。
令 M = yf (xy) , N = xg(xy)
∂(M ⋅ μ) − ∂(N ⋅ μ)
∂y
∂x
f ′( f − g) − ( f ′ − g ′) f g ′( f − g ) − ( f ′ − g ′)g
解 由题得
dx = xy + x3 y3 , dy 这是以 x 为未知函数和以 y 为自变量的迫努利方程,则有 x −3 dx = x −2 y + y 3 ,
dy
令 z = x −2 ,
dz = −2 yz − 2 y3 , dy
而 dz = −2 yz 的解为 z = C~e− y 2 。 dy
充分性。设方程 dy − f (x, y)dx = 0 有仅依赖于 x 的积分因子 µ(x) ,即
µ(x)dy − µ(x) f (x, y)dx = 0
为恰当方程,有
∂(−µ(x) f (x, y)) = dµ(x) ,
∂y
dx
− µ(x) ∂f (x, y)) = dµ(x) ,
∂y
dx
∂f (x, y) 1 dµ(x)
解 由导数的定义可得
x′(t) = lim x(t + s) − x(t)
s →0
s
x(s) + x 2 (t)x(s) = lim
s→0 [1 − x(t)x(s)]s
1 + x 2 (t) x(s)
= lim
⋅,
s→0 1 − x(t)x(s) s
显然可得 x(0) = 0 ,故
x′(t) = [1 + x 2 (t)]⋅ lim x(s) − x(0) = x′(0) ⋅[1 + x 2 (t)]
g
x
(
y x
)
=
−
y x2
g′
=
1 N2
(M x N
− NxM),
y gy ( x)
=
1 x
g′ =
1 N2
(M y N
− MN y ) ,
∂(µM )
− ∂(µN )
=
−N2
xy x2
g′+ N2
y g′ x
=
N2(y x
−
y)g′ x
= 0,
∂y
∂x
(xM + gN )2
(xM + gN )2
因而 µ 是齐次方程的积分因子。
=
( xM
1 +
yN ) 2
[ ∂M ∂y
( xM
+
yN )
−
∂(xM + ∂y
yN )
M
]
=
( xM
1 + yN)2
[ yN
∂M ∂y
−
yM
∂N ∂y
− NM ],
∂( N )
∂(µN) = xM + yN
∂x
∂x
=
( xM
1 +
yN ) 2
[ ∂N ∂x
(xM
+
yN ) −
∂(xM + ∂x
yN )
采用常数变易法,令 z = C~( y)e−y2 代入 dz = −2zy − 2 y3 中得 dy
C~( y) = − y 2 e y 2 + e y 2 + C , 故 z = − y 2 + 1+ Ce−y2 ,
从而原方程的解为
x2 (− y 2 +1 + Ce−y2 ) = 1 。 评注:在微分方程中,变量 x 与 y 具有同等的地位,对同一个方程,既可以就 y 求解, 也可以就 x 进行求解,如果方程 dy = f (x, y) 就 y 求解比较困难,可以尝试将原方程变化
dx − dy = d ln x ,
xy
y
所以原方程变为
1
x
{(M (x, y)x + N (x, y) y)d ln(xy) + (M (x, y)x − N (x, y) y)d ln } = 0 。
2
y
用 µ(x, y) =
1
乘上式两边,得
M (x, y)x + N (x, y) y
1 d ln(xy) + 1 M (x, y)x − N (x, y) y d ln x = 0 ,
子的充要条件是
∂[μ(x + y)M (x, y)] ∂[ μ(x + y)N (x, y)]
=
。
∂y
∂x
则有
μ(x +
y)( ∂M
−
∂N ) =
N
∂μ(x +
y) − M
∂μ(x +
y)
,
∂y ∂x
∂x
∂y
即
μ(x +
y)( ∂M
−
∂N ) =
N
dμ(x +
y)
−M
dμ(x +
y)
,
∂y ∂x
d(x + y)
M (x, y)dx + N (x, y)dy = 0
∂M
中的函数满足关系
− ∂N
=
Nf (x) − Mg ( y) ,其中
f (x), g( y) 分别为
x和
y 的连续函
∂y ∂x
∫ ∫ 数,试证方程 M (x, y)dx + N (x, y)dy = 0 有积分因子 µ = exp( f (x)dx + g ( y)dy) 。
=−
,
∂y
µ(x) dx
上式右端仅为 x 的函数,令其为 P (x) ,积分上式,得
f (x, y) = P(x) y + Q(x) ,
故该方程为线性方程。 评注:一阶线性方程一般用常数变易法求解,此例给出了线性方程的又一种求解方法, 即积分因子法。
2-7 设函数 f (u), g(u) 连续、可微且 f (u) ≠ g(u) ,试证方程 yf (xy)dx + xg(xy)dy = 0
N]
=
( xM
1 + yN)2
[ xM
∂N ∂x
−
xN
∂M ∂x
−
NM ] ,
∂(µM ) − ∂(µN )
∂y
∂x
= x(NM x − MN x ) + y(NM y − MN y ) , (xM + yN )2