下载文档原格式
1 / 16习题２-１判断下列方程是否为恰当方程,并且对恰当方程求解： １．0)12()13(2=++-dy x dx x解：13),(2-=x y x P , 12),(+=x y x Q ,则0=∂∂y P ,2=∂∂xQ, 所以 x Q y P ∂∂≠∂∂ 即 原方程不是恰当方程． ２．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=∂∂xQ所以x Q y 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 xQ =∂∂ 所以x Q y 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 xQ=∂∂ 因为 0≠b , 所以x Q y P ∂∂≠∂∂,即 原方程不为恰当方程５．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 xQ=∂∂ 所以x Q y P ∂∂=∂∂,即 原方程为恰当方程则,0cos )sin 2cos (2=++udu udt t udu t2 / 16两边积分得：.sin )1(2C u t =+ ６．0)2()2(2=++++dy xy e dx y e ye xxx解： xy e y x Q y e ye y x P xxx2),(,2,(2+=++=,则,2y e y P x +=∂∂,2y e xQx +=∂∂ 所以x Q y P ∂∂=∂∂,即 原方程为恰当方程则,0])2()[(22=++++dy xy e dx y ye dx e xxx两边积分得：.)2(2C xy e y x=++７．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 x y两边积分得：23ln 3y x y x -+.C = ８．),(0)(22为常数和c b a cxydy dx by ax =++解：,),(,),(22cxy y x Q by ax y x P =+=则,2by y P =∂∂,cy xQ =∂∂ 所以 当x Q y P ∂∂=∂∂,即 c b =2时, 原方程为恰当方程则0)(22=++cxydy dx by dx ax两边积分得：233bxy ax +.C = 而当c b ≠2时原方程不是恰当方程．９．01222=-+-dt ts s ds t s 解：,),(,12),(22ts s s t Q t s s t P -=-= 则,212t s t P -=∂∂,212tss Q -=∂∂ 所以x Q y P ∂∂=∂∂, 即原方程为恰当方程,两边积分得：C ts s =-2．3 / 1610．,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 xQ '=∂∂ 所以x Q y P ∂∂=∂∂, 即原方程为恰当方程,两边积分得：22()f xy dx C +=⎰,即原方程的解为C y x F =+)(22<其中F 为f 的原积分>．习题２-2. 1. 求解下列微分方程,并指出这些方程在平面上的有意义 的区域:：〔１〕yx dx dy 2= 解：原方程即为：dx x ydy 2= 两边积分得：0,2332≠=-y C x y ．〔２〕)1(32x y x dx dy += 解：原方程即为：dx xx ydy 321+=4 / 16两边积分得：1,0,1ln 2332-≠≠=+-x y C x y ．〔３〕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 ．〔４〕221xy y x dx dy +++=；解：原方程即为：2(1)1dyx dx y =++ 两边积分得：c x x arctgy ++=22, 即 )2(2c x x tg y ++=． 〔５〕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 ∈〕 〔６〕21y dxdyx-= 解：①当1±≠y 时 原方程即为：xdx y dy =-21 两边积分得：c x y =-ln arcsin ． ②1±=y 也是方程的解.〔７〕．yxe y e x dx dy +-=- 解．原方程即为：dx ex dy e y xy)()(--=+5 / 16两边积分得：c e x e y x y++=+-2222, 原方程的解为：c ee x y xy=-+--)(222.2. 解下列微分方程的初值问题． 〔１〕,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 ． 〔２〕．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．〔３〕．r d dr=θ, 2)0(=r ； 解：原方程即为：θd rdr=,两边积分得：c r =-θln ,因为2)0(=r , 所以2ln =c ,所以原方程满足初值问题的解为：2ln ln =-θr 即θe r 2=．〔４〕．,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 ++-=6 / 16〔５〕．321xy dxdyx=+, 1)0(=y ； 解：原方程即为：dx xx y dy 231+=, 两边积分得：c x y ++=--22121, 因为1)0(=y , 所以23-=c ,所以原方程满足初值问题的解为：311222=++yx ．１． 解下列微分方程,并作出相应积分曲线的简图． 〔１〕．x dxdycos = 解：两边积分得：c x y +=sin ． 积分曲线的简图如下：〔２〕．ay dxdy=, 〔常数0≠a 〕； 解：①当0≠y 时,原方程即为：dx ay dy = 积分得：c x y a +=ln 1, 即 )0(>=c cey ax②0=y 也是方程的解． 积分曲线的简图如下：7 / 16〔３〕．21y dxdy-=； 解：①当1±≠y 时,原方程即为：dx y dy =-)1(2 积分得：c x yy+=-+211ln ,即 1122+-=x x ce ce y ．②1±=y 也是方程的解．积分曲线的简图如下：〔４〕．n y dx dy =, )2,1,31(=n ； 解：①当0≠y 时, ⅰ〕2,31=n 时,原方程即为 dx y dy n =,积分得：c y n x n=-+-111．8 / 16ⅱ〕1=n 时,原方程即为dx ydy= 积分得：c x y +=ln ,即 )0(>=c ce y x．②0=y 也是方程的解．积分曲线的简图如下：4. 跟踪：设某A 从xoy 平面上的原点出发,沿x 轴正方向前进；同时某9 / 16B 从点开始跟踪A,即B 与A 永远保持等距b ．试求B 的光滑运动轨迹．解：设B 的运动轨迹为)(x y y =,由题意与导数的几何意义,则有22yb ydx dy --=,所以求B 的运动轨迹即是求此微分方程满足b y =)0(的解．解之得：222222ln21y 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 轴的方向平移而得到。
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 x N ∂∂+N xu ∂∂⇔ 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)) 由已知条件上式恒成立，故原命题得证。
习题2.3解：卫j ，赳=1 ..y；x则 W =ex所以此方程是恰当方程。
凑微分，x 2dx -2ydy (ydx xdy)二 0得:1 x 3 xy _ y 2 二 C3(y 「3x 2)dx 「(4y 「x)dy 二 0.:MN ._ 1 , — 1._y ；x则如」cy ex所以此方程为恰当方程。
凑微分，ydx xdy - 3x 2dx - 4ydy = 0得 x 3 -xy 2y 2 二 C型 二 2x(x-y)2-2x 2(x-y) _ 2xy:x(x_y)4 (x_y)31、 验证下列方程是恰当方程，并求出方程的解。
1. 2(x y)dx (x - 2y)dy = 02.解:3.y 21 1 x 2EYdx [「E ]dy，解:.:M 2y(x-y)2-2y 2(x-y)(-1) 2xy3-y(x-y)4「 (x - y)(x-y)4则也=卫ex cy1 x 2y (x-y)2因此此方程是恰当方程。
2.:u y 1 2 _.x (x-y) x2.:u 1 x■7(x-y)2(2)(1) 做x 的积分，1u2dx dx (y)-In x 「(y)(3) 做y 的积分，孔-(-1)y 2 (x-y)2y d (y):y(x-y)2dy=-2xy y 2 d(y)(x-y)2dydy (y)二y2 - 2xy 1 X2 -2xy y2 y (x - y)2(x - y)2(x-y)24、i(1)dy=l ny_y yy u 二- 一- In x In y - y =In x—y x 故此方程的通解为In 1二C2(3xy2 2x3)dx 3(2x2y y2)dy = 0 解:.:M :N12xy ,12xy .:M :N■:y :x则此方程为恰当方程。
凑微分，6xy 2dx 4x 3dx 6x 2ydy 3y 2dy = 03d(x 2y 2) d(x 4) d(x 3) = 0得：x 43x 2y 2y 3= Cd(-cos x)+d (sin- )+dx+d(- 1)=0yx y所以，d(sin y -cos- +x -1 )=0 xy y 故所求的解为sin$ -cos- +x -- =Cx y y求下列方程的解：2 26. 2x(y e x -1)dx+ e x dy=0解：迴二 2x e x2, 岂=2x e"cyex所以，迴二型，故原方程为恰当方程2 2又 2xy e x dx-2xdx+ e x dy=0解:sin x yyM=y 12cos 丄cos* -冷 xy sin y x 1 、+ r )dy=0ycos y+1xcos yxcM 鋼 _ 1 2y sin xy cN1x2sin；:x yy 所以=cN :x 因为 1 . sin xdx-cos yxcos y +-y3sin$x x x,故原方程为恰当方程yy1 y x2cos 丄 dx+dx+ — cos- dy-sin xdy+丄 dy=0y yyyxx 13 cos------------- 2yy xy xcy ex所以，d(ye"-x2)=0故所求的解为ye"-x2=C7.(e x+3y2)dx+2xydy=0解：e x dx+3y2dx+2xydy=0e x x2dx+3x2y2dx+2x3ydy=0所以，d e x( x2-2x+2)+d( x 3y2)=0即 d [e x( X2-2X+2)+ x3y2]=0故方程的解为e x( x2-2X+2)+ x3y2=C8.2xydx+( x2+1)dy=0解：2xydx+ x2dy+dy=Od( x2y)+dy=0即d(x2y+y)=0故方程的解为x2y+y=C9、ydx -xdy = x2 y2 dx解：两边同除以x2 y2得yd x ~ x d y -dxx + yf 、即，d arctg — = dx< y丿故方程的通解为argtg - =x + c ly丿10、ydx - X y3 dy =0解：方程可化为：ydx了叽ydyy/ \即， d — = ydyly丿故方程的通解为：-=-^2c 即：2x=yy 2cy 2同时，y=0也是方程的解。
常微分方程与动力系统第二章习题参考答案 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,+∞上有界。
第二章 一阶微分方程的初等解法§2.1 变量分离方程与变量变换习题2.1求下列方程的解 1．xy dxdy2=，并求满足初始条件：1,0==y x 的特解． 解 分离变量，得到xdx ydy 2=，两边积分，即得C x y ~ln 2+=，因而，通解为 2x Ce y =，这里C 是任意常数．此外，方程还有解0=y ．由10==x y 得1=C ，特解2x e y =．2．0)1(2=++dy x dx y ，并求满足初始条件：1,0==y x 的特解． 解 分离变量，得到12+-=x dx y dy ，两边积分，即得C x y ~1ln 1++-=-，因而，通解为Cx y ++=1ln 1，这里C 是任意常数．此外，1-=x 和0=y 是两条积分曲线．由10==x y 得1=C ，特解11ln 1++=x y ．3．yx xy y dx dy 321++=． 解 分离变量，得到)1(122x x dx y ydy +=+，两边积分，即得C xx y ~1ln )1ln(222++=+，所以得通解222)1)(1(Cx y x =++，这里0>C 是任意正常数．4．0)1()1(=-++xdy y ydx x ．解 分离变量，得到dx xx dy y y +=-11，两边积分，即得C x x y y ~ln ln ++=-，因此得通解C xy y x =+-ln ，这里C 是任意常数．另有特解0=x 和0=y ．5．0)()(=-++dx y x dy x y ．解 变形得x y x y dx dy +-=，这是齐次方程，设x y u =，得dxdu x u dx dy +=，代入原方程得 11+-=+u u dx du xu ，分离变量得 x dx du u u -=++211，两边积分，即得 C x u u +-=++ln )1ln(21arctan 2，即C y x x y =++)ln(21arctan 22，这里C 是任意常数．6．22y x y dxdy x -+=．解 变形得 21sgn ⎪⎭⎫ ⎝⎛-+=x y x x y dx dy ，这是齐次方程，设x y u =，得dx du x u dx dy +=，代入原方程得 21sgn u x dxdux-=，分离变量积分，即得 C x x u +=ln sgn arcsin ，即C x x xy+=ln sgn arcsin． C x u u +-=++ln )1ln(21arctan 2，即C y x x y =++)ln(21arctan 22，这里C 是任意常数． 7．0cot tan =-xdy ydx ．解 分离变量，得到xdx ydy tan cot =，两边积分，即得C x y ~cos ln sin ln +-=，所以通解为C y x =sin cos ，这里0≠C 的任意常数．另有特解πk y =，Z k ∈及2ππ+=k x ，Z k ∈，这只须在通解表达式中允许0=C 即可，故通解为C y x =sin cos ，这里C 是任意常数．8．032=++ye dx dy x y ．解 分离变量，得到dx e e ydy x y32-=，两边积分，即得C e e x y ~3232+-=--，得到通解C eey x=--2323，这里C 是任意常数．9．0)ln (ln =--ydx dy y x x ．解 变形得y y x x dy dx )ln (ln -=，令u y x =，则dydu y u dy dx +=，代入方程并分离变量得，ydy u u du =-)1(ln ，两边积分，即得C y u ~ln 1ln ln +=-，或1ln +=Cy u ，回代原变量有，1ln+=Cy yx ，或1+=Cy ye x ，这里0≠C 的任意常数．另有特解满足01ln =-u ，即ey x =，这只须在通解表达式中允许0=C 即可，故通解为1+=Cy ye x ，这里C 是任意常数．10．y x e dxdy-=． 解 分离变量，得到dx e dy e xy=，积分得C e e xy +=，这里C 是任意常数．作适当的变量变换求解下列方程（11—17） 11．2)(y x dxdy+=． 解 设y x u +=，则dx dy dx du +=1，原方程化为C x u u dxdu+=⇒+=arctan 12，即通解为 C x y x +=+)a r c t a n(，这里C 是任意常数． 12．2)(1y x dx dy +=． 解2)(y x dydx+=，由上题，注意到这里的x 和y 相当于上题的y 和x ，得到方程的通解为 C y y x +=+)a r c t a n(，这里C 是任意常数． 13．1212+-+-=y x y x dx dy ． 解 由⎩⎨⎧=+-=+-012,012y x y x 得⎪⎪⎩⎪⎪⎨⎧=-=31,31y x ，令⎪⎪⎩⎪⎪⎨⎧-=+=,31,31y Y x X 就有Y X Y X dX dY 22--=，这是齐次方程，令X Y u =，有u dXduX dX dY +=，代入方程后分离变量，X dX du u u u =+--)1(2212，得到C X u u ~ln 2)1ln(2+-=+-，回代变量得C y x y xy x =-++-22即为原方程的通解，这里C 是任意常数．14．25--+-=y x y x dx dy ． 解 令u y x =-，则dx du dx dy =-1，代入方程得 27--=u dx du ，分离变量并积分得，C x u u =+-1442，即C y x y xy x =+++-410222为方程的通解，这里C 是任意常数．15．18)14()1(22+++++=xy y x dxdy． 解 变形为2)14(2+++=y x dx dy ，令u y x =++14，则dxdydx du 41+=，代入原方程得942+=u dx du ，分离变量解之得，C x u +=632arctan ，回代原变量并变形化简，得到通解 )14(2)6t a n (3++=+y x C x ，这里C 是任意常数．16．2252622yx xy x y dx dy +-=． 解 变形为232632)2(3)(xxy x y dx y d +-=，令x y u 3=，则原方程化为 1262+--=u u u dx du x ，解之得537)2()3(Cx u u =+-，即153373)2()3(Cx x y x y =+-为方程的通解，这里C 是任意常数．17．yy y x xxy x dx dy -+++=32232332． 解 变形为123132)()(222222-+++=y x y x x d y d ，令⎪⎩⎪⎨⎧==22,yY x X ，原方程变为123132-+++=Y X Y X dX dY ，由⎩⎨⎧=-+=++0123,0132Y X Y X ，得到⎩⎨⎧-==1,1Y X ．设⎩⎨⎧+=-=1,1Y v X u ，则有v u v u dv du 2332++=，再令s v u =，得到dvds v s dv du +=，于是23)1(32+-=s s dv ds v ，解得C v s s =-+65)1)(1(，逐步回代变量，得原方程的通解为C y x y x =--+52222)2)((，这里C 是任意常数．18．证明方程)(xy f dxdyy x =经变换u xy =可化为变量分离方程，并由此求解下列方程： （1）xdy dx y x y =+)1(22；（2）222222y x y x dx dy y x -+=． 证明 令u xy =，则得dx dy x y dx du +=，代入原方程得]1)([+=u f u dxdu x 是变量分离方程．（1）中221)(y x xy f +=，所以)2(2+=u u dxdux，分离变量求解得 C u x u ++=)]2(ln[2arctan224，即得原告方程的通解C y x x xy ++=)]2(ln[2arctan2224．（2）中2222)(u u u f -+=，所以224udx du x -=，分离变量求解得 C x u u +=-ln 43123，即得原告方程的通解 C x y x xy +=-ln 431233． 19．已知0,1)()(0≠=⎰x dt t f x f x，试求函数)(x f 的一般表达式．解 变形后等式两边对x 求导，有 ])(1[])(['='⎰x f dt t f x，即 )()()(2x f x f x f '-=，解得)(21)(C x x f +±=，由1)()1(1=⎰dt t f f ，得0=C ，所以xx f 21)(±=．20．求具有性质)()(1)()()(s x t x s x t x s t x -+=+的函数)(t x ，已知)0(x '存在．解 因为)0(x '存在，故)(t x 在0=t 连续，即)0()(lim 0x t x x =→．由)()(1)()()(s x t x s x t x s t x -+=+，令0=s 就有)0()(1)0()()(x t x x t x t x -+=，得到0)0(=x ．ss x s x t x t x s t x s x t x s x t x s t x s t x )()()(1)(1)()()(1)()()()(2⋅-+=--+=-+，令0→s 取极限，由于右边的极限为)0()](1[2x t x '+，故左边的极限存在，从而得到函数)(t x 满足的方程， )0()](1[)(2x t x t x '+='，解之得 C t x t x +'=)0()(arctan ，或])0(t a n [)(C t x t x +'=．由0)0(=x ，推出Z k k C ∈=,π，所以])0(tan[)(πk t x t x +'=，Z k ∈．21．求一曲线，使它的切线介于两坐标轴之间的部分被切点分成相等的部分． 解 由习题 1.2—9（4），知曲线)(x f y =应满足的方程0=+'y y x ，即xy dx dy -=，分离变量解之得，C x y ~ln ln +-=，或C xy =为所求的曲线．22．在图（2.1）所示的C R -电路中，设10=E 伏，100=R 欧，01.0=C 法，而开始时电容C 上没有电荷，问：（1）当开关K 合上“1”后，经过多长时间电容C 上的电压5=C u 伏？（2）当开关K 合上“1”后，经过相当长的时间（如1分钟后）开关K 从“1”突然转至“2”，试求C u 的变化规律，并问经过多长时间5=C u 伏？解 （1）由例7，)1(1t RCC eE u --=，将10=E ，100=R ，01.0=C 代入，有)1(10t C e u --=，由)1(105te --=，反解出)(6931.02ln s t ≈=，即经过约6931.0秒，电容C 上的电压5=C u 伏．（2）同样由例7，t RCC Eeu 1-=，代入具体数值有t C e u -=10，由te -=105，同样得到)(6931.02ln s t ≈=，即经过约6931.0秒，电容C 上的电压5=C u 伏．23．求出习题1.2第9题（1）所确定的曲线，其中4πα=．解 由习题1.2—9（1），ααt a n t a n y x x y y -+='，代入4πα=得y x x y y -+='，这是齐次方程，令u x y =，则dx du x u dx dy +=，代入得2211u u dx du x -+=，解出C y x xy++=)ln(arctan 222即为所求曲线．24．证明满足习题1.2第9题（7）所给条件的曲线是抛物线族． 证明 由习题1.2—9（7），0(>='k kx y 常数)，解之得C kx y +=221，这是抛物线族，顶点在),0(C ，对称轴为y 轴．§2.2 线性方程与常数变易法习题2.2求下列方程的解： 1．x y dxdysin +=． 解 首先，求齐次线性方程y dxdy=的通解，从dx y dy =得到齐次方程通解x ce y =，令xe x c y )(=为方程的解，代入得x ex c xsin )(-='，即cx x e x c x ~)co s (si n 21)(++-=-，故原方程的通解为x e c x x y ~)cos (sin 21++-=，其中c ~为任意常数． 2．t e x dt dx 23=+．解 由03=+x dtdx ，解出t ce x 3-=，设t e t c x 3)(-=是原方程的解，代入原方程得，t e t c 5)(='，故c e t c t ~51)(5+=，所以原方程的通解为t t e e c x 2351~+=-，其中c ~为任意常数．3．t t s dt ds 2sin 21cos +-=． 解 由t s dt ds cos -=，解得tce s sin -=，设t e t c s s i n )(-=是原方程的解，代入原方程得，t e t c t 2sin 21)(sin ='，得c t e t c t ~)1(sin )(sin +-=，所以通解为1sin ~sin -+=-t e c x t ，其中c ~为任意常数．4．n x e y x n dx dy n x ,=-为常数．解 由0=-y xndx dy ，解得n cx y =，设n x x c y )(=是原方程的解，代入原方程得，x e x c =')(，即ce x c x ~)(+=，所以通解为n x x c e y )~(+=，这里c ~为任意常数． 5．01212=--+y x xdx dy ． 解 由0212=-+y xx dx dy ，解得x e cx y 12=，设xe x x c y 12)(=是原方程的解，代入原方程得，ce x c e xx c x x ~)(1)(112+=⇒='--，所以通解x e x c x y 122~+=，这里c ~为任意常数．6．234xyy x dx dy +=． 解 原方程即231y x y x dx dy +=，这是2-=n 的Bernoulli 方程，令3y z =，就有，233x z xdx dz +=，解这个一阶线性方程得通解为)ln 3(3c x x z +=，即)ln 3(33c x x y +=，这里c 为任意常数．7．3)1(12+++=x x y dx dy ．解 由12+=x y dx dy ，得2)1(+=x c y ，令2)1)((+=x x c y 为原方程的解，代入原方程得，1)(+='x x c ，即c x x c ~)1(21)(2++=，所以原方程通解为24)1(~)1(21+++=x c x y ，其中c ~为任意常数．8．3y x ydx dy +=． 解 变形为21y x ydy dx +=，把x 看作未知函数，y 看作自变量，对于y 及dy dx 来说，这是一个线性方程．先解对应的齐线性方程x ydy dx 1=，得cy x =，其次把c 看作)(y c ，即设y y c x )(=为变形后方程的解，代入变形后的方程得y dy y dc =)(，得到c y y c ~21)(2+=，从而原方程的通解为y c y x ~213+=，其中c~为任意常数． 9．xx x ay dx dy 1++=． 解 先解xay dx dy =，得a cx y =，设ax x c y )(=为原方程的解，代入原方程得，11)(++='a x x x c ，即 ⎪⎪⎪⎩⎪⎪⎪⎨⎧=++=+-≠≠+--=-0,~ln ,1,~1ln ,0,1,~)11()(a c x x a c x x a a c x a a xx c a ， 所以原方程通解为⎪⎪⎪⎩⎪⎪⎪⎨⎧=++=+-≠≠+--=0,~ln ,1,~1ln ,0,1,~11a c x x a x c x x a a x c a a x y a ，其中c ~为任意常数． 10．3x y dxdyx =+． 解 先解y dx dy x +，得xc y =，设x x c y )(=为原方程的解，代入原方程得，3)(x x c =' ，即c x x c ~41)(4+=，所以原方程通解为x c x y ~413+=，这里c~为任意常数． 11．33y x xy dxdy =+．解 这是3=n 的Bernoulli 方程，令2-=y z 代入有322x xz dxdz -=，解这个一阶线性方程得通解为122++=x cez x ，即1)1(222=++x ce x y 为原方程的通解，这里c 为任意常数．另有特解0=y ．12．xdy ydx x y =-)2ln (．解 变形为y xy x x dx dy 2ln 2-=，这是2=n 的Bernoulli 方程，令1-=y z 代入有 x xz x dx dz ln 2+-=， 解这个一阶线性方程得通解为241ln 21cx x z ++=，即1)41ln 21(2=++cx x y 为原方程的通解，这里c 为任意常数．另有特解0=y ．13．dx x y xydy )2(22-=． 解 变形为yy x dx dy 211-=，这是1-=n 的Bernoulli 方程，令2y z =代入有 12-=z xdx dz ，解这个一阶线性方程得通解为2cx x z +=，即22cx x y +=，这里c 为任意常数．14．23x x e dx dy y +=． 解 设u e y=，则dx dy u dx dy e dx du y ==，代入原方程得2213u xu x dx du +=，这是2=n 的Bernoulli 方程，令1-=u z 代入有 213xz x dx dz --=，解这个关于z 的一阶线性方程得通解为3~21xcx z +-=，回代原变量得原方程的通解322)(x e x c y =-，其中c 为任意常数．15．331yx xy dx dy +=． 解 变形为33x y yx dydx +=，把x 看作未知函数，y 看作自变量，对于y 及dy dx 来说，这是一个3=n 的Bernoulli 方程．令2-=x u ，有322y yu dydu--=，解这个一阶线性方程得通解为122+-=-y ce u y ，即得原方程的通解1222+-=--y ce x y ，这里c 为任意常数．16．⎰+=xx dt t y e y 0)(．解 两边求导得一阶线性方程x e y dxdy+=，解之得通解x e c x y )(+=，从原方程知道有初始条件10==x y ，代入通解表达式中得1=c ，故原积分方程的解为xe x y )1(+=．17．设函数)(t ϕ于+∞<<∞-t 上连续，)0(ϕ'存在且满足关系式)()()(s t s t ϕϕϕ=+，试求此函数．解 由于ss t st s t st s t 1)()()()()()()(-⋅=-=-+ϕϕϕϕϕϕϕ，且)0(ϕ'存在，故在该式中令0→s 取极限就有，)0()()(ϕϕϕ'='t t ，解得t ce t )0()(ϕϕ'=．若0)(≡t ϕ，则是解；若)(t ϕ不恒为零，则由)0()()0()(ϕϕϕϕt t t =+=得1)0(=ϕ，由此得1=c ，所以t e t )0()(ϕϕ'=．18．如图所示的L R -电路，试求：（1）当开关1K 合上10秒后，电感L 上的电流；（2）1K 合上10秒后再将2K 合上，求2K 合上20秒后，电感L 上的电流． 解 （1）由Kirchhoff 第二定律得，E dtdIL I R =+1，把101=R ，2=L ，50=E 代入得到微分方程255=+I dtdI，初始条件0=t 时，0=I ．解之得t e I 555--=，当10=t 时，5055--=eI 约为5安培．（2）由Kirchhoff 第二定律得，E dt dILRI =+，其中320201020102121=+⨯=+=R R R R R ，2=L ，50=E 代入得25310=+I dt dI ，初始条件0=t 时，5=I ．解之得)3(25310t e I --=，当20=t 时，)3(253200--=eI 约为7.5安培．19．试求图示的L R -电路电感上电流)(t I 的变化规律，并解释其物理意义，设0=t 时，0=I ．解 由Kirchhoff 第二定律得，E dtdILRI =+，即t L U I L R dt dI m ωsin =+，初始条件为00==t I，求出其通解为)sin(222ϕωω-++=-t L R U ce I m t LR，其中RL ωϕ=tan ，20πϕ<<．由初始条件得，ϕωsin 222L R U c m +=，所以，)]sin([sin 222ϕωϕω-++=-t e L R U I t LR m ．其物理意义是：当t 增大时，第一项逐渐衰减而趋于零（称为暂时电流），事实上很快就消失而不起作用．而第二项就起着重要作用（称为稳定电流）．稳定电流是一个周期函数，其周期与电动势的周期相同，而相角相差ϕ-．20．试证：（1）一阶非齐线性方程（2.28）的任两解之差必为相应的齐线性方程（2.3）之解；（2）若)(x y y =是（2.3）的非零解，而)(~x y y =是（2.28）的解，则方程（2.28）的通解可表为)(~)(x y x cy y +=，其中c 为任意常数； （3）方程（2.3）任一解的常数倍或任两解之和（或差）仍是（2.3）的解． 证明 （1）设)(1x y y =，)(2x y y =是方程（2.28）)()(x Q y x P dxdy+= 的任意两个解，即)()(11x Q y x P dx dy +=，)()(22x Q y x P dxdy+=，由此得到 ))(()]()([)]()([)(21212121y y x P x Q y x P x Q y x P dxdy dx dy dx y y d -=+-+=-=-，所以)()(21x y x y y -=是齐线性方程（2.3）：y x P dx dy)(=之解． （2）由于)(x y y =是（2.3）的非零解，故)()()(x y x P dxx dy =，而)(~x y y =是（2.28）的解，即)()(~)()(~x Q x y x P dx x y d +=，所以)]()(~)([)()()(~)())(~)((x Q x y x P x y x cP dxx y d dx x dy c dx x y x cy d ++=+=+)()](~)()[(x Q x y x cy x P ++=，所以)(~)(x y x cy y +=是（2.28）的解，其中含有一个任意常数c ，故是方程（2.28）的通解，其中c 为任意常数．（3）设)(1x y y =，)(2x y y =都是方程（2.3）的解，即)()()(11x y x P dx x dy =， )()()(22x y x P dxx dy =， 因此有)]()[()()()())((1111x ky x P x y x kP dxx dy k dx x ky d ===，)]()()[()()()()()()()]()([21212121x y x y x P x y x P x y x P dxx dy dx x dy dx x y x y d ±=±=±=±，所以，方程（2.3）任一解的常数倍或任两解之和（或差）仍是（2.3）的解． 21．求解习题1.2第9题（5）和（6）．解 （5）的方程为 2x y x y ='-，或变形为x y xy -='1，这是一阶线性方程． 先解对应的齐次方程y xy 1='，得到cx y =，设原方程的解为)(x xc y =，代入原方程得1)(-='x c ，即c x x c +-=)(，故所求曲线方程为cx x y +-=2，其中c 为任意常数．（6）的方程为 x y x y ='-2，或变形为2121-='y x y ，这是一阶线性方程．同样先解对应的齐次方程y xy 21='，得到x c y =，设原方程的解为x x c y )(=，代入原方程得xx c 21)(-='，即c x x c +-=)(，故所求曲线方程为x c x y +-=，其中c 为任意常数．22．求解下列方程：（1）01)1(2=+-'-xy y x ；（2）0)12()1(322=+--'-x y x y x x ； （3）0sin cos sin 3=--'x y x x y ．解 （1）先解xy y x ='-)1(2，得12-=x cy ，设方程的解为1)(2-=x x c y ，代入方程得)1sgn(1)(2232---='-x x x c ，推出1)(2-+=x c x x y 为原方程的通解（需分1>x ，1-<x 及1<x 三种情形分别求解后再统一），这里c 为任意常数．（2）先解y x y x x )12()1(22-='-，得到12-=x cxy ，设原方程的解为1)(2-=x x x c y ，代入原方程得 1)1()(22---='x x x x c ，即c x x c +-=11)(2，所以原方程的通解为12-+=x cxx y ，这里c 为任意常数．（3）先解0cos sin =-'y x x y ，得到x c y tan =，设原方程的通解为x x c y tan )(=，代入原方程得x x c sin )(='，即c x x c +=co s )(，所以通解x c x y tan sin +-=，这里c 为任意常数．§2.3 恰当方程与积分因子习题2.3验证下列方程是恰当方程，并求出方程的解： 1．0)2()(2=-++dy y x dx y x ． 证明 y x N y x M 2,2-=+=，所以xNy M ∂∂==∂∂1，即所给方程是恰当方程． 改写方程为02)(2=-++ydy xdy ydx dx x ，即0)31(23=-+y xy x d ，得原方程的解为c y xy x =-+2331，其中c 为任意常数． 2．0)4()3(2=---dy x y dx x y ． 证明 )4(,32x y N x y M --=-=，所以xNy M ∂∂==∂∂1，即所给方程是恰当方程． 改写方程为043)(2=--+ydy dx x xdy ydx ，即0)2(23=--y x xy d ，得原方程的解为c y x xy =--232，其中c 为任意常数．3．0])(1[]1)([2222=--+--dy y x x y dx x y x y ．证明 2222)(1,1)(y x x y N x y x y M --=--=，所以x N y x xy y M ∂∂=-=∂∂3)(2，即所给方程是恰当方程．改写方程为 0)(222=+---y dyx dx y x dy x dx y ， 凑为0)()()()(2=-+----y dyx dx y x y x xyd xy d y x ， 即0)ln ln (=-+-y x y x xy d ，得原方程的通解为c y x yx xy=-+-ln ln ，其中c 为任意常数．4．0)2(3)23(22232=+++dy y y x dx x xy ． 证明 )2(3,)23(22232y y x N x xy M +=+=，所以xNxy y M ∂∂==∂∂12，即所给方程是恰当方程．改写方程为034)(623=+++dy y dx x xdy ydx xy ，即0)3(3422=++y x y x d ，得原方程的解为c y x y x =++34223，其中c 为任意常数．5．0)1sin cos 1()1cos sin1(222=+-++-dy yy x y x x y x dx x y x y y x y ． 证明 由于2221sin cos 1,1cos sin 1yy x y x x y x N xyx y y x y M +-=+-=， 所以，x Ny x y x y x y x y x y x y xy M ∂∂=--+-=∂∂cos sin 1sin cos 13232，即所给方程是恰当方程． 改写方程为01cos sin 222=++-+-dy ydx x y x ydx xdy y x y xdy ydx ， 即0)1cos (sin=-+-y x y x x y d ，得原方程的解为c yx y x x y =-+-1cos sin ，其中c 为任意常数．求下列方程的解：6．0)1(222=+-dy e dx ye x x x ．解 改写方程为02)2(22=-+⋅xdx dy e y dx xe x x ，即0)(22=-x ye d x ，所以得到原方程的通解c x ye x =-22，这里c 为任意常数．7．02)3(2=++xydy dx y e x． 解 由于xy N y e M x2,32=+=，故y xNy y M 2,6=∂∂=∂∂． 因为xN xNy M 2=∂∂-∂∂只与x 有关，所以方程有只与x 有关的积分因子 2ln 22x eexdxx ==⎰=μ，以2x =μ乘方程两边得，0233222=++ydy x dx y x dx e x x，即0)()(223=+xe d x y x d ，故得原方程的通解为c e x x y x x=+-+)22(223，这里c 为任意常数．8．0)1(22=++dy x xydx ．解 改写为0)2(2=++dy dy x xydx ，凑微分得0))((22=++dy dy x x yd ，得原方程的通解c y y x =+2，其中c 为任意常数．9．dx y x xdy ydx )(22+=-． 解 以22y x +除方程两边，有dx yx xdy ydx =+-22，即dx y xd =)(arctan ，得到原方程的通解为c x yx+=arctan，这里c 为任意常数． 10．0)(3=+-dy y x ydx ． 解 改写为dy y xdy ydx 3=-，得ydy yxdy ydx =-2，即)21()(2y d y x d =，所以得到原方程的通解c y y x +=221，或cy y x +=321，其中c 为任意常数． 11．0)1(=+--xdy dx xy y ．解 由x N xy y M =--=,1，得1,1=∂∂-=∂∂x Nx y M ，由于1-=∂∂-∂∂Nx Ny M 与y无关，故方程有只与x 有关的积分因子x dxe e --=⎰=)1(μ，以x e -乘方程两边有，0)1(=+----xdy e dx xy y e x x ，分组得，0])([=--+---dx e dx xye xdy ydx ex x x，凑微分得0])1[(=+-x e xy d ，即得方程的通解为xce xy =+1，这里c 为任意常数．12．0)(2=--xdy dx x y ．解 由x N x y M -=-=,2，得1,1-=∂∂=∂∂xN y M ，由于x N x Ny M 2-=∂∂-∂∂只与x 有关，故方程有积分因子2)2(1x edxx =⎰=-μ，以21x乘方程两边并组合变形有， dx xxdyydx =-2， 即dx x y d =-)(，得到方程的通解为c x xy-=-，或)(x c x y -=，这里c 为任意常数．13．0)2(=++xdy dx y x ．解 改写为02=++xdy ydx xdx ，显然有积分因子x ，故以x =μ乘方程两边有，0])([222=++dy x x yd dx x ，即0)()31(23=+y x d x d ，得到通解c y x x =+2331，其中c为任意常数．14．0)cos()]sin()cos([=+++++dy y x x dx y x y x x ． 解 改写为 0)sin())(cos(=++++dx y x dy dx y x x ， 即 0)sin())(sin(=+++dx y x y x xd ，或0)]sin([=+y x x d ，所以原方程的通解为c y x x =+)sin(，其中c 为任意常数． 15．0)cos sin ()sin cos (=++-dy x x x y dx x x x y ．解 x x x y N x x x y M c o s s i n ,s i n c o s+=-=，则x x x x y xNsi n cos cos -+=∂∂，x yMcos =∂∂，由于1=-∂∂-∂∂M x Ny M 与x 无关，故方程有积分因子y dy e e =⎰=μ，以y e 乘方程两边并分项组合有，0sin cos )1()cos sin cos (=+-+⋅+-xdy ye xdx y e dy e x x xdx xe xdx e y y y y y ，或写为0)(sin )(sin )1()](cos ))(cos (cos [=+-+++y y y y e xd y x d y e e xd x x xd xdx e ，即0sin )](sin )(sin )[1()cos (=++-+xdy e e xd x d e y xe x d yyyy， 也即0)1(sin )sin ()1()cos (=-+-+y xd e x e d y xe x d y y y ，故0)sin )1(cos (=-+x e y xe x d yy，得到方程的通解为c e x y x x y=-+]sin )1(cos [，这里c 为任意常数．16．0)53()24(3=+++xdy ydx y xdy ydx x ．解 改写方程为 05324342=+++dy xy dx y dy x xydx ，可看出y x 2=μ是一个积分因子，用它乘方程两边有053244352423=+++dy y x dx y x ydy x dx y x ，分项组合就有，0)()(5324=+y x d y x d ，故方程的通解为c y x y x =+5324，其中c 为任意常数．17．试导出方程0),(),(=+dy y x N dx y x M 分别具有形为)(y x +μ和)(xy μ的积分因子的充要条件．解 设)(y x +μ是0),(),(=+dy y x N dx y x M 的积分因子⇔0),()(),()(=+++dy y x N y x dx y x M y x μμ是恰当方程 ⇔xN y M ∂∂=∂∂)()(μμ ⇔NM y x xNy M y x d y x d -+∂∂-∂∂-=++)()()()(μμ ⇔NM y x xNy M -+∂∂-∂∂-)()(μ应为y x +的函数)(y x f +． 又设)(xy μ是0),(),(=+dy y x N dx y x M 的积分因子⇔0),()(),()(=+dy y x N xy dx y x M xy μμ是恰当方程⇔xN y M ∂∂=∂∂)()(μμ ⇔xMyN xy xNy M xy d xy d -∂∂-∂∂=)()()()(μμ ⇔)()()(xy xMyN xy xNy M ϕμ=-∂∂-∂∂． 18．设),(y x f 及yf∂∂连续，试证方程0),(=-dx y x f dy 为线性方程的充要条件是它有仅依赖于x 的积分分子．证明 “⇒”即例4．“⇐” 若方程0),(=-dx y x f dy 有反依赖于x 的积分因子，则)()(x p yfx N y M -=∂∂-=∂∂-∂∂ 仅与x 有关，所以)()()(),(x Q y x p dy x p dy y fy x f +==∂∂=⎰⎰，其中)(x Q 是x 的任意连续函数．从而方程为0)]()([=+-dx x Q y x p dy ，即)()(x Q y x p dxdy+=是线性方程． 19．试证齐次方程0),(),(=+dy y x N dx y x M 当0≠+yN xM 时有积分因子yNxM +=1μ．证明 将方程两端同乘以N 1，得0=+dy dx N M ，即0)(=+dy dx xy g ． 设xyu =，则ux y =，从而0)(=++xdu udx dx u g ，或0])([=++xdu dx u u g ，这是可分离变量方程，取积分因子])([11u u g x +=μ，则有0])(1[ln =++⎰du uu g x d ，得到通解为c du u u g x =++⎰)(1ln ，其中c 为任意常数．积分因子为yN xM xy N M x Nu u g x N N +=+⋅=+⋅==1][11])([1111μμ，通解可写为c xyd yN xM N x =++⎰)(ln ，c 为任意常数．20．设函数)(u f ，)(u g 连续、可微且)()(u g u f ≠，试证方程0)()(=+dy xy xg dx xy yf有积分因子1)])()([(--=xy g xy f xy μ．证明 以1)])()([(--=xy g xy f xy μ乘以方程两边，有0)]()([)()(=-+xy g xy f xy dyxy xg dx xy yf ，或0)]()([))((=--+y dy xy g xy f xy xdy ydx xy f ，即0)ln )]()([)((=--⎰y du u g u f u u f d )(xy u =，因而方程0)()(=+dy xy xg dx xy yf 有积分因子1)])()([(--=xy g xy f xy μ．21．假设方程（2.43）中的函数),(,),(y x N y x M 满足关系)()(y Mf x Nf xNy M -=∂∂-∂∂ 其中)(,)(y g x f 分别为x 和y 的连续函数，试证方程（2.43）有积分因子))()(exp(⎰⎰+=dy y g dx x f μ．证明 因为yMy g M y M y M M y ∂∂+=∂∂+∂∂=∂∂μμμμμ)()(，xNx f N x N x N N x ∂∂+=∂∂+∂∂=∂∂μμμμμ)()(， 所以0][)]()([)()(=∂∂-∂∂+-=∂∂-∂∂xNy M x Nf y Mg N x M y μμμμ， 即)()(N xM y μμ∂∂=∂∂，从而方程0=+Ndy Mdx μμ是恰当方程，故方程（2.43）有积分因子))()(exp(⎰⎰+=dy y g dx x f μ．22．求出Bernoulli 方程的积分因子． 解 Bernoulli 方程为n y x Q y x P dx dy )()(+= )1,0(≠n ，以ny n --)1(乘方程两边，并令u yu=-1，化为关于dx du u ,的一阶线性方程)()1()()1(x Q n u x P n dxdu-+-=，后者有积分因子⎰--dx x P n e )()1(，从而Bernoulli 方程的积分因子⎰-=--dxx P n ne yn )()1(1μ． 23．设),(y x μ是方程（2.43）的积分因子，从而求得可微函数),(y x U ，使得)(Ndy Mdx dU +=μ．试证),(~y x μ也是方程（2.43）的积分因子的充要条件是)(),(~U y x μϕμ=，其中)(t ϕ是t 的可微函数． 证明 “⇐”若),(y x μ是方程（2.43）的积分因子，且)(Ndy Mdx dU +=μ，则))(()()()()(~0⎰==+=+=dU U d dU U U Ndy Mdx Ndy Mdx ϕϕϕμμ， 所以c dU U =⎰)(ϕ为（2.43）的通解，故)(),(~U y x μϕμ=亦是方程（2.43）的积分因子，其中)(t ϕ是t 的可微函数．“⇒”设dV Ndy Mdx =+22μμ，则M xV2μ=∂∂，N y V 2μ=∂∂． 由Ndy Mdx dU μμ+=，则M xUμ=∂∂，N y U μ=∂∂，得到yU x UN M y V x V ∂∂∂∂==∂∂∂∂，所以，0=∂∂∂∂∂∂∂∂yU x U yV x V，因而存在函数)(t Φ，使得)(U V Φ=，由此得 ))(())(()(Ndy Mdx U Ndy Mdx U dU U dV +Φ'=+Φ'=Φ'=μμμ )(2N d y M d x+=μ， 得到)()(2U U μϕμμ=Φ'=．24．设),(,),(21y x y x μμ是方程（2.43）的两个积分因子，且21μμ不恒为常数，求证c =21μμ（任意常数）是方程（2.43）的通解．证明 由于21,μμ是方程（2.43）的两个积分因子，由上题结论)(21U ϕμμ=，其中)(2Ndy Mdx dU +=μ，这里)(t ϕ是t 的可微函数．由于)(21U ϕμμ=不恒为常数，故有)(U ϕ'不恒为零，由此在c U =)(ϕ两边微分得0)()(2=+'Ndy Mdx U μϕ，因此得到，0)(2=+Ndy Mdx μ，所以c U =)(ϕ是方程（2.43）的解，又c U =)(ϕ中含有一个任意常数，故c U =)(ϕ即c =21μμ（任意常数）是方程（2.43）的通解．25．假设第19题中微分方程还是恰当的，试证它的通解可表为c y x yN y x xM =+),(),(（c 为任意常数）． 证明 由于方程是恰当的，故11=μ即是一个积分因子，而由第19题yNxM +=12μ也是积分因子，且yN xM +=21μμ不恒为常数，所以由第24题所证结论，就知道它的通解为c y x yN y x xM =+),(),(，c 为任意常数．§2.4 一阶隐方程与参数表示习题2.4求解下列方程： 1．y y x '+='13． 解 解出31y y x ''+=，设p y ='，方程为31ppx +=，两边对y 求导，有 dydpp p p )23(134+-=， 即dp p p dy )23(23+-=，所以c p pp ++=2232，因此得原方程的通解为 ⎪⎪⎩⎪⎪⎨⎧++=+=c p py pp x 223,11223（p 为参数），c 为任意常数．2．0)1(33='--'y x y ．解 设31t y ='-，则31t y -='，21t t x -=．由dt t tt dx y dy )21)(1(23---='=，得c t t t y ++-=52521，从而通解为 ⎪⎪⎩⎪⎪⎨⎧++-=-=c t t t y t tx 522521,1，（t 为参数），c 为任意常数． 3．y e y y ''=2．解 设p y ='，则原方程为pe p y 2=，两边对x 求导有，dxdppep p p)2(+=，即dp e p dx p )2(+=，解得c e p x p ++=)1(，所以通解为⎪⎩⎪⎨⎧=++=ppep y c e p x 2,)1( （p 为参数），c 为任意常数． 4．a y y 2)1(2='+（a 为常数）． 解 解出212y a y '+=，p y ='，则原方程为212pay +=，两边对x 求导有， dx dpp p a p ⋅⋅+-=2)1(222，或dp p a dx 22)1(4+-=，解得c p a papx +-+-=arctan 2122，所以通解为 ⎪⎪⎩⎪⎪⎨⎧+=+-+-=2212,arctan 212p a y c p a p ap x （p 为参数），c 为任意常数． 5．122='+y x ．解 设t x cos =，t y sin ='，则由tdt dx y dy 2sin -='=，得c t t y ++-=2s i n 4121，所以通解为⎪⎩⎪⎨⎧++-==c t t y t x 2sin 4121,cos （t 为参数），c 为任意常数． 6．22)2()1(y y y '-=-'．解 令yt y ='-2，则有t t y -=1，所以21t y +='，dt ty dy dx 21-='=，由此解出c tx +=1，于是求得通解为⎪⎪⎩⎪⎪⎨⎧-=+=tt y c t x 1,1 （t 为参数），或消去参数t 得c x c x y +--=1，c 为任意常数．习题2.5求下列方程的解：1．1cos sin =+x dx dyx y ． 解 原方程为x y x dx dy sec tan +⋅-=，由y x dxdy⋅-=tan ，得到x c y cos =，设原方程的解是x x c y cos )(=，代入原方程得出x x c 2sec )(='，即c x x c +=tan )(，因此原方程的通解为x c x y cos sin +=，c 为任意常数．2．ydy x xdy ydx 2=-． 解 方程两边同乘以21x ，有ydy xxdy ydx =-2，凑微分得0)21(2=+x yy d ，故得通解c xyy =+221，这里c 为任意常数． 3．1sin 4-=-x e dxdy y ．解 改写为0sin 4=-+xdx dx e dy e yy，两边乘以xe 并凑微分得0)sin 4(=-⎰xdx e e e d x y x ，所以c xdx e e e xy x +=⎰sin 4，即c x x e e e xy x +-=)cos (sin 2，其中c 为任意常数．4．xyx ydx dy -=． 解 这是齐次方程．设x y u =，则ux y =，dxdu x u dx dy +=，代入原方程化为 uuu dx du x-=1． 分离变量求解得，c uux 22ln =+，即2)ln 21(y c y x -=，这里c 为任意常数．5．0)(22=-+dy e x dx y xye yx yx ．解 变形为yxe x y x y dx dy -⎪⎭⎫ ⎝⎛+=2，这是齐次方程．设x y u =，则ux y =，dx du x u dx dy +=，代入化原方程为u e u dxdux 12-=，分离变量求解得x c e u ln 1=+-，即c e x y x=+ln ，其中c 为任意常数．6．0)1(=-+xdy ydx xy ． 解 变形为02=-=yxdy ydx xdx ，凑微分得0)21(2=+y xy d ，所以原方程的通解为c yxy =+221，其中c 为任意常数． 7．0)2()122(=-++-+dy y x dx y x ．解 变形为2)(1)(2-+-+-=y x y x dx dy ，设u y x =+，则dxdudx dy =+1，代入原方程后得，21-+-=u u dx du ，解之得c x u u +-=+-1ln 3，即c y x y x +++=+1ln 2，这里c 为任意常数．8．32x y x y dx dy +=． 解 这是2=n 的Bernoulli 方程．令1-=y z ，有311xz x dx dz --=，解这个一阶线性方程，得x cxz +=21，即x c x y +=211，这里c 为任意常数． 9．23-+=x y dxdy． 解 先解y dxdy3=，得到x ce y 3=，设原方程的解是x e x c y 3)(=，代入原方程后得，x e x x c 3)2()(--='，所以c e x x c x +--=-3)53(91)(，得到x ce x y 3)53(91+--=是原方程的通解，这里c 为任意常数．10．21⎪⎭⎫ ⎝⎛+=dx dy dx dy x ．解 设p dxdy=，则p p x +=1，两边对y 求导得dy dp p p )11(12+-=，从中就可解出c p p y +-=ln 212，所以通解为⎪⎪⎩⎪⎪⎨⎧+-=+=c p p y p px ln 21,12（p 为参数），c 为任意常数．11．312+++-=y x y x dx dy ． 解 改写并分项组合，有0)1()3()(2=+-+=+dx x dy y ydx xdy ，凑微分，得0)21331(22=--++x x y y xy d ， 所以，c x x y y xy =--++2221331是方程的通解，这里c 为任意常数． 12．x y xe dxdy e =+-)1(． 解 原方程即y x xe dx dy +=+1，设u y x =+代入方程得u xe dxdu=，这是分离变量方程．解出c e x u =+-221，即得原方程的通解为c e x y x =++-)(221，这里c 为任意常数．13．02)(22=-+xydy dx y x ．解 设xy N y x M 2,22-=+=，则y x N y y M 2,2-=∂∂=∂∂，x N x Ny M 2-=∂∂-∂∂仅与x 有关，故方程有积分因子2)2(1xedxx =⎰=-μ，用它乘方程两边并分项组合有， 0222=-+x xydy dx y dx ，即0)(2=-x y x d ，所以c xy x =-2，或cx y x =-22是原方程的通解，其中c 为任意常数． 14．1++=y x dx dy． 解 由y dxdy =，解得x ce y =，设原方程的解为x e x c y )(=，代入有xe x x c -+=')1()(，即c ex x c x++-=-)2()(，所以通解为)2(+-=x ce y x ，其中c 为任意常数．15．xy e dx dy x y+=．解 设u x y =，则ux y =，dx du x u dx dy +=，代入化原方程为u e dxdu x =，分离变量解之得c e x u=+-ln ，即c ex xy=+-ln ，其中c 为任意常数．16．y e dxdyx -=++21)1(． 解 分离变量得112+=--x dx e dy y ，两边积分得通解12+=-x c e y，c 为任意常数． 17．0)1()(2=++-dy x y dx y x ．解 改写为1111-+-+=y xxy x dx dy ，这是1-=n 的Bernoulli 方程．设2y z =，则原方程化为一阶线性方程xx z x dx dz +-+=1212，解之得2)1(12+++=x c x z ，因此得原方程的通解为22)1(12+++=x c x y ，这里c 为任意常数．18．0)1(24322=-+dy y x dx y x ．解 )1(2,4322-==y x N y x M ，则y x xN y x y M 226,8=∂∂=∂∂，y M x Ny M 21-=-∂∂-∂∂只与y 有关，故有积分因子yedyy 1)21(=⎰=-μ，用它乘以方程并分项组合有02)24(21213232=-+-dy y dy y x dx y x ，凑微分得，0)434(21233=-y y x d ，所以通解为c y y x =-212333，或c y y x =-)3(3，其中c 为任意常数．19．0422=+-⎪⎭⎫⎝⎛x dx dy y dx dy x ．解 解出y y x y '+'=2)4(2，设p y ='，则原方程为pp x y 2)4(2+=，两边对x 求导有，dxdpp x p p p )221(2422-++=， 或dp pxdx =，解得cp x =，所以通解为⎪⎩⎪⎨⎧+==)4(2,2p c y cp x （p 为参数）， 或消去参数p ，得2242c x cy +=，c 为任意常数．另外还有042=-p ，或x y 2±=也是解．20．1]1[22=⎪⎭⎫⎝⎛-dx dy y ．解 令t dxdysin =，代入方程有1)sin 1(22=-t y ，即t y sec ±=． 由于dt tdt t t t y dy dx 2cos 1sin tan sec ±=±='=，所以1tan c t x +±=，得到原方程的通解为⎩⎨⎧±=+±=ty c t x sec ,tan 1 （t 为参数）， 消去参数t 得1)(22++=c x y ，其中c 为任意常数．21．0)1()1(=-++dy yxe dx e yx yx ． 解 设u y x =，则yu x =，dy duy u dy dx +=，代入原方程化简得uu e u e dy du y ++-=1，分离变量求解得c y e u u=+)(，即c yex yx=+是原方程的通解，其中c 为任意常数．22．0324223=-+dy yx y dx y x ． 解 设42233,2y x y N y x M -==，则x N y x y M ∂∂=-=∂∂46，故为恰当方程． 由于)()(2),(323y yx y dx y x y x u ϕϕ+=+=⎰，其中)(y ϕ是y 的待定可微函数，再由422423)(3y x y y y x y u -='+-=∂∂ϕ，得到21)(y y ='ϕ，即有y y 1)(-=ϕ，因此得到方程的通解为c y yx y x u =-=1),(32，即322cy y x =-，这里c 为任意常数．23．0)1(2=++-dy y x ydx ．解 变形为dy y xdy ydx )1(2+=-，看出有积分因子21y =μ，用21y =μ乘以方程两边并凑微分得)1()(y y d y xd -=，即得方程的通解是c yy y x +-=1或cy y x =+-12，这里c 为任意常数．24．0)]([22=-+-xdy dx y x x y ．解 变形为dx y x x xdy ydx )(22+=-，看出有积分因子221yx +=μ，用它乘以方程两边并凑微分得)2()(arctan 2x d y x d =，得方程的通解是c x y x +=2arctan 2，或其等价形式)2tan(2c x y x +=，其中c 为任意常数． 25．0=-+x e dxdydx dy． 解 设p dxdy=，则p e p x +=，两边对y 求导有dy dp e p p )1(1+=，即dp e p dy p )1(+=，由此得到c e p p y p +-+=)1(212，所以方程的通解为 ⎪⎩⎪⎨⎧+-+=+=ce p p y e p x pp )1(21,2 （p 为参数），c 为任意常数． 26．0)()32(2232=++++dy y x dx y y x xy ． 解 设2232,32y x N y y x xy M +=++=，则x x Ny x x y M 2,222=∂∂++=∂∂，1=∂∂-∂∂N xNy M 与y 无关，所以方程有积分因子x dx e e =⎰=1μ，以之乘方程的两边，分项组合得到0)3()2(2322=++++dy e y dx e y dy e x dx ye x dx xye x xxxx，即。
第二章 基本定理我们在第一章主要学习了初等积分法，掌握了几类常微分方程的解法.但是这些解法只适用于某些特殊的类型，很多其它的常微分方程不能用初等解法进行求解.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[上的近似解.。
第二章、一阶微分方程的初等解法[教学目标]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 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=- 解 将变量分离，得到两边积分，即得因而，通解为22x y c += 这里的c 是任意的正常数.或解出显式形式例2 解方程并求满足初始条件：当0x =时.1y =的特解.解 将变量分离，得到两边积分，即得因而，通解为这里的c 是任意的常数.此外，方程还有解0y =.为确定所求的特解，以0x =.1y =代入通解中确定常数c ，得到 1c =-因而，所求的特解为例3 求方程()dy P x y dx= （2.3） 的通解，其中()P x 是x 的连续函数.解 将变量分离，得到两边积分，即得这里的c 是任意常数.由对数的定义，即有即令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 >有 事实上，取1t x=，则方程可改写成形如(2.5)的方程. ⅱ)对方程(,)dy f x y dx =其中右端函数(,)f x y 是x 和y 的零次齐次函数，即对0t >有则方程也可改写成形如(2.5)的方程对齐次方程（2.5）利用变量替换可化为变量分离方程再求解.令y u x =（2.6） 即y ux =，于是 dy du x u dx dx=+ （2.7） 将（2.6）、（2.7）代入（2.5），则原方程变为整理后，得到 ()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 tgu dx x= （2.9） 分离变量，即有两边积分，得到 这里的c 是任意的常数，整理后，得到sin u cx = （2.10）此外，方程（2.9）还有解0tgu =,即sin 0u =. 如果（2.10）中允许0c =，则sin 0u =就包含在（2.10）中，这就是说，方程（2.9）的通解为（2.10）.代回原来的变量，得到原方程的通解为例5 求解方程(0).dy x y x dx +=<解 将方程改写为 这是齐次方程，以,y dy du u x u x dx dx ==+代入，则原方程变为du xdx =（2.11） 分离变量，得到两边积分，得到（2.11）的通解即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 =.原方程的通解还可表为它定义于整个负半轴上.注：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）属齐次方程，有 此时，令y u x=，即可化为变量可分离方程. （2）11220a b a b =，即1122a b a b =的情形. 设1122a b k a b ==，则方程可写成 令22a x b y u +=，则方程化为这是一变量分离方程.（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）化为从而（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）更一般的方程类型此外，诸如以及（其中,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）化为两边积分，得因此记1,c e c ±=并代回原变量，就得此外，易验证即也就是（2.18）的解.因此方程（2.17）的通解为其中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）分离变量，得到两边积分，得到即这里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 ，则由光的反射定律：入射角等于反射角，容易推知从而注意到及,,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 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、一阶线性微分方程在()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 dxy 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 dxy c x e ⎰= （2.29）两边微分，得到()()()()()P x dxP x dx dy dc x e c x P x e dx dx⎰⎰=+ （2.30） 将（2.29）、（2.30）代入（2.28），得到 即 积分后得到()()()P x dxc x Q x e dx c -⎰=+⎰（2.31） 这里c 是任意的常数..将（2.31）代入（2.29），得到()()()()()() =()P x dxP x dx P x dx P x dx P x dxy 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 dyx ny e x dx++-=+的通解，这里的n 为常数. 解 将方程改写为(1)1x n dy n y e x dx x -=++ （2.33） 先求对应的齐次方程 的通解，得令 ()(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），再积分，得 将其代入公式（2.34），即得原方程的通解 这里c 是任意的常数. 例2 求方程22dy y dx x y=-的通解. 解 原方程改写为2dx x y dy y=- （2.36） 把x 看作未知函数，y 看作自变量，这样，对于x 及dxdy来说，方程（2.36）就是一个线性方程了.先求齐线性方程 的通解为2x cy = （2.37）令2()x c y y =,于是 代入（2.36），得到 从而，原方程的通解为这里c 是任意的常数,另外0y =也是方程的解. 特别的，初值问题 的解为例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 是非齐线性方程的两个不同的解，则应满足方程使 （1）—（2）有说明非齐线性方程任意两个解的差12y y -是对应的齐次线性方程的解.（2）因为 故结论成立.（3）因为12121212()()()(),(),()d y y d y y d cy p cy p y y p y y dx dx dx+-==+=- 故结论成立.3、Bernoulli 方程形如()()n dyP x y Q x y dx=+ （ 0,1n ≠） （2.38） 的方程，称为伯努利（Bernoulli ）方程，这里(),()P x Q x 为x 连续函数.利用变量变换可将伯努利方程化为线性方程来求解.事实上，对于0y ≠，用n y -乘（2.38）两边，得到1()()n n dyy y P x Q x dx--=+ （2.39） 引入变量变换1n z y -= （2.40） 从而(1)ndz dyn y dx dx-=- （2.41） 将（2.40）、2.41）代入（2.39），得到(1)()(1)()dzn P x z n Q x dx=-+- （2.42） 这是线性方程，用上面介绍的方法求得它的通解，然后再代回原来的变量，便得到（2.38）的通解.此外，当0n >时，方程还有解0y =. 例4 求方程26dy yxy dx x=-的通解 解 这是2n =时的伯努利方程，令 1z y -=,得 代入原方程得到这是线性方程，求得它的通解为 代回原来的变量y ，得到 或者这是原方程的通解. 此外，方程还有解0y =. 例5 求方程331dy dx xy x y =+的解 解 将方程改写为这是一个自变量为y ，因变量为x 的伯努利方程.解法同上.例6 求方程23y dy e x dx x+=的通解 这个方程只要做一个变换，令,y ydu dyu e e dx dx==，原方程改写为 便是伯努利方程.小结;这次主要讨论了一阶线性微分方程的解法.其核心思想是常数变易法.即将非齐线性方程对应的齐线性方程解的常数变易为待定函数，使其变易后的解函数代入非齐次线性方程，求出待定函数()c x ，求出非齐次方程的解.我们还讨论了伯努利方程，求解过程为，先变换，将原方程化为非齐线性方程，再求解.§3 恰当方程与积分因子1、恰当方程的定义 将一阶微分方程 写成微分的形式把,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 ≡,于是就是方程(2.43)的隐式通解,这里C 是任意常数(应使函数有意义). 2、 恰当方程的判定准则定理1设(,),(,)M x y N x y 在某区域G 内连续可微,则方程(2.43)是恰当方程的充要条件是, (,)M Nx y G y x∂∂=∈∂∂ (2.46) 而且当(2.46)成立时,相应的原函数可取为 00(,)(,)(,)xyx y u x y M s y ds N x t dt =+⎰⎰ (2.47)或者也可取为0(,)(,)(,)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 是连续可微的,从而有下面证明定理的充分性,即由条件(2.46),寻找函数(,)u x y ,使其适合方程(2.45).从(2.47)可知即(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)有 再选取00(,)(0,1)x y =-,代入公式(2.47)有 可见不论0y >和0y <,都有故方程的通解为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)便得方程的通解为 解法2. 分项组合的方法 对(2.49)式重新组合变为于是 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 ，使则称(,)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 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）进一步知 即由()μμϕ=可知左端是ϕ的函数，可见右端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）的解.事实上，因为因此函数（2.54）的确是方程（2.43）的积分因子.为了方便应用这个定理，我们就若干特殊情形列简表如下： 例2.解22(31)()0y xy dx xy x dy -++-=解 这里2231,M y xy N xy x =-+=-,注意 所以方程不是恰当的,但是 它仅是依赖与x ，因此有积分因子 给方程两边乘以因子x μ=得到 从而可得到隐式通解例3. 解方程2()(1)0xy y dx xy y dy ++++=解 这里2,1M xy y N xy y =+=++方程不是恰当的.但是 它有仅依赖于y 的积分因子 11dyy eyμ-⎰≡=方程两边乘以积分因子1y μ=得到 1()(1)0x y dx x dy y++++= 从而可得到隐式通解另外，还有特解0y =.它是用积分因子乘方程时丢失的解. 例4. 解方程 223(2)()0y x y dx xy x dy +++=解 这里2232,M y x y N xy x =+=+，不是恰当方程.设想方程有积分因子()x y αβμμ=，其中α，β是待定实数.于是只须取3,2αβ==.由上述简表知原方程有积分因子 从而容易求得其通解为： 六、积分因子的其他求法以例4为例，方程的积分因子也可以这样来求：把原方程改写为如下两组和的形式：前一组有积分因子11y μ=，并且 后一组有积分因子21xμ=，并且 设想原方程有积分因子其中α，β是待定实数.容易看出只须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 ≠，则有积分因子 并且通解为例6、试用积分因子法解线性方程（2.28）.解 将（2.28）改写为微分方程[()()]0P x y Q x dx dy +-= （2.55）这里()(),1M P x y Q x N =+=-，而 则线性方程只有与x 有关的积分因子方程（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）为恰当方程，又分项分组法 因此方程的通解为 即与前面所求得的结果一样.注：积分因子一般不容易求得可以先从求特殊形状的积分因子开始，或者通过观察法进行“分项分组”法求得积分因子.§4 一阶隐方程与参数表示1、一阶隐方程一阶隐式微分方程的一般形式可表示为:如果能解出(,)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 dp p x y dx∂∂=+∂∂ (2.59) 方程(2.59)是关于,x p 的一阶微分方程,而且属于显式形式.若求得(2.59)的通解形式为(,)p x c ϕ=,将其代入(2.58),于是得到(2.57)通解为 若求得(2.59)的通解形式为(,)x p c ψ=,于是得到(2.57)的参数形式的通解为 其中p 为参数, c 是任意常数.若求得(2.59)的通解形式为(,,)0x p c Φ=,于是得到(2.57)的参数形式的通解为 其中p 为参数, c 是任意常数.例1 求方程3()20dy dy x y dx dx +-= 的解 解 令dy p dx=,于是有32y p xp =+ (2.60) 两边对x 求导数,得到即 2320p dp xdp pdx ++=当0p ≠时,上式有积分因子p μ=,从而由此可知得到将其代入(2.60),即得故参数形式的通解为当0p =时,由(2.60)可知0y =也是方程的解.例2 求方程22()2dy dy x y x dx dx =-+的解. 解 令dy p dx=,得到222x y p xp =-+ (2.61) 两边对x 求导数,得到 2dp dp p p x p x dx dx =--+ 或 (2)(1)0dp p x dx--= 由10dp dx -=,解得p x c =+,于是得到方程的通解为222x y cx c =++ (2.62)由20p x -=,解得2x p =,于是得到方程的一个解为24x y = (2.63) 特解(2.63)与通解(2.62)中的每一条积分曲线均相切,因此称为方程的奇解.2) 讨论形如 (,)dy x f y dx= (2.64) 的方程的求解方法,方程(2.64)与方程(2.57)的求解方法完全类似,假定函数(,)f y y ' 有连续偏导数.引进参数dy p dx=,则(2.64) 变为 (,)x f y p = (2.65)将(2.65) 的两边对y 求导数,得到1f f dp p y x dy∂∂=+∂∂ (2.66) 方程(2.66))是关于,y p 的一阶微分方程,而且属于显式形式.设其通解为 则(2.64)的通解为Ⅱ）不显含y （或）x 的方程3) 讨论形如(,)0F x y '= (2.67) 的方程的解法. 记dy p y dx'==,此时(,)0F x p =表示的是xp 平面上的一条曲线,设曲线用参数形式表示为()x t ϕ=,()p t ψ= (2.68) 由于dy pdx =,进而两边积分,得到于是得到方程(2.67)参数形式的解为 c 是任意常数.例3 求解方程3330x y xy ''+-=解 令y p tx '==,则由方程得于是 23339(12)(1)t t dy dt t -=+ 积分得到故原方程参数形式的通解为：4) 讨论形如(,)0F y y '= (2.69) 的方程,其解法与方程(2.67)的求解方法类似.记dy p y dx'==,此时(,)0F y p =表示的是yp 平面上的一条曲线,设曲线用参数形式表示为 由关系式dy pdx =可知 ()()t dt t dx ϕψ'=,于是0p ≠时,有 故方程(2.69)的参数形式的通解c 是任意常数.此外,不难验证,若(,0)0F y =有实根y k =,则y k =也是方程的解. 例4 求解方程 22(1)(2)y y y ''-=-.解 令2y yt '-=,则有由此可以得 代入1dx dy p=，得到 积分,得到1x c t =+故原方程参数形式的通解为其中c 是任意常数.此外, 当0y '=时原方程变为24y =,于是2y =±也是方程的解.例5 求解方程y '=解 令y p '=,则有p =,取,(,)22p tgt t ππ=∈-,则sin sec tgt x t t === 由dy pdx =得到所以cos y t c =-+故原方程参数形式的通解为其中c 是任意常数.。
习题２-4１．求解下列微分方程：（１）yx xy y --='22；解：令ux y =，则原方程化为uu u dx du x --=+212，即x dxdu u u =--122，积分得：c x u u u +=--+-ln 1ln 2111ln2 还原变量并化简得：3)()(y x c x y +=-（２）4252--+-='y x x y y ；解：由⎩⎨⎧=--=+-042052y x x y 得 ⎩⎨⎧-==21y x令2,1+=-=y v x u ， 则有vu u v du dv --=22，由第一题的结果知此方程解为3)()(v u c u v +=-， 还原变量并化简得：.)1(33++=+-y x c x y（３）14212-+++='y x y x y ；解：令y x v 2+=， 则1212121-++=+=v v dx dy dx dv ， 即1214-+=v v dx dv ，此方程为变量分离方程， 分离变量并积分得：c x v v +=+-14ln 8321，还原变量并化简得：c y x x y =++--184ln 348． （４）xy y x y -='33．解：①当0≠y 时，方程两边同时乘以32--y ,则233222--+-='-xy x y y ， 令2-=y z ， 则322x xz dxdz-=, 此方程为一阶线性方程，由公式得：122++=x ce z x还原变量得：122)1(2-++=x ce y x . ②0=y 也是方程的解．2. 利用适当的变换，求解下列方程： （１）)cos(y x y -='；解：令y x u -=，则u dx dy dx du cos 11-=-=， ①当1cos ≠u 时，有dx udu =-cos 1， 即 dx u du=2sin 22，两边积分得：c x uctg +=221还原变量化简得：2sin 2sin 22cos yx c y x x y x -+-=-． ②当1cos =u 时，即πk x y 2+=)(Z k ∈也是方程的解． （２）0)()3(22=+++dv uv u du v uv ； 解：方程两边同时乘以u 则原方程化为：0)()3(2322=+++dv v u u du uv v u ，即 0)()3(2232=+++vdv u du uv dv u vdu u 此方程为全微分方程，则原方程的解为：c v u v u =+22321． （３）)2(2)3(222yx y x dx dy y x -=++；解：原方程即为324222222++-=y x x y xdx ydy ，令u y v x ==22,，则324++-=v u vu dv du ,由⎩⎨⎧=++=-03024v u v u 得⎩⎨⎧-=-=21v u ， 令⎩⎨⎧+=+=21v n u m ，则有n m n m dn dm +-=24令z n m=，则zn m =， 124+-=+=z z z n dn dz dn dm ， 则有1)2)(1(+--=z z z n dn dz ，此方程为变量分离方程， 分离变量并积分得：n c zz ln 2)1(ln32+=--，还原变量并化简得：322222)32()1(-+-=+-y x c y x .（４）yy y x xxy x dx dy 8237323223-+-+=． 解：原方程即为823732222222-+-+=y x y x xdx ydy ，令22,x v y u ==,则823732-+-+=u v u v dv du ，由⎩⎨⎧=++=-+08230732u v u v ⎩⎨⎧==⇒21v u ， 令⎩⎨⎧-=-=21v n u m ， 则m n m n dn dm 2332++=，令z n m=，可将方程化为变量分离形方程， n dn dz zz =-+)2223(2，两边积分得：c n z z z +=---+ln 1ln 2111ln 432， 还原变量并化简得：)3()1(22522-+=--y x c y x ．3. 求解下列微分方程： （１）．2241xy y --='； 解：令xy z =， 则原方程可化为：)41(12-+-=z z x dx dz ， ①当21≠z 时，即21≠xy 时方程为x dxdz z =--2)21(1 ，此方程为变量分离方程， 两边积分得：c x z +=-ln 211还原变量并化简得：cxx x x y ++=ln 121； ②当21=z 时，xy 21=是方程的特解． （２）．1222++='xy y x y x ； 解：原方程即为：221x x y y y ++='， 令xy z =，则2)1(1+=z xdx dz ，此方程为变量分离方程， 分离变量积分得：c x z +=+-ln 11， 还原变量并化简得：cxx x x y +--=ln 11． 4. 试把二阶微分方程0)()(=+'+''y x q y x p y 化为一个黎卡提方程． 解：令⎰=udxe y ， 则⎰='udxue y ，+⎰=''udxe u y 2⎰'udxe u ，代入原方程可得：=+'+''y x q y x p y )()(+⎰udxe u 2⎰'udxe u ＋)()(x q ue x p udx+⎰⎰udxe ＝０，即有：0)()(2=++'+x q u x p u u ，此方程为一个黎卡提方程．5. 求一曲线，使得过这一曲线上任一点的切线与该点向径的夹角等于45．解：设此曲线为)(x y y =，由题意得：1451==+-tg xy dx dy x y dx dy ，化简得：y x y x dx dy -+=， 此方程为齐次方程，解之得：c y x x y arctg =+-)ln(2122．6. 探照灯的反光镜（旋转面）应具有何种形状，才能使点光源发射的光束反射成平行线束？解：取点光源所在处为坐标原点，而x 轴平行于光的反射方向，建立三维坐标系．设所求曲面由曲线⎩⎨⎧==0)(z x f y 绕x 轴旋转而成,则求反射镜面问题归结为求 xy 平面上的曲线y=f(x)的问题．由题意及光的反射定律，可得到函数)(x f y =所应满足的微分方程式：22yx x ydx dy ++=，此方程为齐次方程， 解之得：)2(2x c c y +=，（其中c 为任意正常数）．)2(2x c c y +=就是所求的平面曲线，它是抛物线，因此反射镜面的形状为旋转抛物面)2(22x c c z y +=+．习题２-5１．求解下列微分方程：（１）．0)()23(2232=++++dy y x dx y xy y x ；解：方程两边同乘xe33， 则)33()369(233323323=++++dy y e dx y e dy x e xydx e ydx x e x x x x x ，此方程为全微分方程，即 c y e y x e x x =+33233． （２）．0)2(2=-+-dy e xy ydx y ；解：方程两边同乘y e y 21， 则 0)12(22=-+dy yxe dx e y y即01)2(22=-+dy ydy xe dx e yy 此方程为全微分方程，即有 c y xe y =-ln 2 ．（３）．0)3()63(2=+++dy xyy x dx y x ；解：方程两边同乘 xy ， 则0)3()63(232=+++dy y x dx x y x即 0)36()3(232=+++dy y xdx dy x ydx x 此方程为全微分方程，即有c x y y x =++2333 ．（４）．22()0ydx x y x dy -++=; 解：方程两边同乘221y x +， 则 022=-+-dy yx xdyydx ， 此方程为全微分方程，即 c y yxarctg=- （5）．0)1(2223=-+dy y x dx xy ；解：方程两边同乘21y ， 则0)1(222=-+dy y x xydx ， 此方程为全微分方程，即c y x y=+21. （6）．0)1(=-+xd y dx xy y ；解：方程两边同乘21y ， 则0)1(2=-+dy y xdx y xdx ， 此方程为全微分方程，即c x y x =+221. （7）0)(2223=-+dy xy x dx y ；解：方程两边同乘y x 21， 则 02)2(22=+-dy y dy x y dx x y ， 此方程为全微分方程，即 c y xy =+-ln 22（8）．0)c o s2(=++dy y y ctgy e dx e xx解：方程两边同乘y sin ， 则02sin )cos sin (=++ydy yc ydy e ydx e x x ，此方程为全微分方程，即 11cos cos 2sin 224xe y y y y c -+=. 2. 证明方程（5.1）有形如)),((y x φμμ=的积分因子的充要条件是)),((y x f yP P x Q Q xQy P φ=∂∂-∂∂∂∂-∂∂，并写出这个积分因子。
第二章、一阶微分方程的初等解法[教学目标]1. 理解变量分离方程以及可化为变量分离方程的类型（齐次方程），熟练掌握变量分离方程的解法。
2. 理解一阶线性微分方程的类型，熟练掌握常数变易法及伯努力方程的求解。
3. 理解恰当方程的类型，掌握恰当方程的解法及简单积分因子的求法。
4. 理解一阶隐式方程的可积类型，掌握隐式方程的参数解法。
[教学重难点] 重点是一阶微分方程的各类初等解法，难点是积分因子的求法以及隐式方程的解法。
[教学方法] 讲授，实践。
[教学时间] 14学时[教学内容] 变量分离方程，齐次方程以及可化为变量分离方程类型，一阶线性微分方程及其常数变易法，伯努利方程，恰当方程及其积分因子法，隐式方程。
[考核目标]1.一阶微分方程的初等解法：变量分离法、一阶线性微分方程的常数变易法、恰当方程与积分因子法、一阶隐方程的参数解法。
2.会建立一阶微分方程并能求解。
§1 变量分离方程与变量变换1、变量分离方程1) 变量分离方程形如（2.1）的方程，称为变量分离方程.2) 求解方法(2.1)可化为，这样变量就分离开了,两边积分,得到（2.2）.容易验证由（2.2 2.1）.因而（2.2）是（2.1）的通解.2.1）的解.可能它不包含在方程的通解（2.2）中，必须予以补上.3) 例题例1解将变量分离，得到两边积分，即得因而，通解为.或解出显式形式例2 解方程.解将变量分离，得到两边积分，即得因而，通解为.因而，所求的特解为例3 求方程（2.3）.两边积分，即得.由对数的定义，即有即（2.4）2.3）的解.如果在（2.4（2.4 2.3）的通解为（2.4.注: 1.(2.2)式有意义.2.方程的通解不一定是方程的全部解,有些通解包含了方程的所有解，有些通解不能包含方程的所有解.此时，还应求出不含在通解中的其它解,即将遗漏的解要弥补上.3..2、可化为变量分离方程的类型1）.形如（2.5）.另外,ⅰ)对于方程(2.5)的方程.ⅱ)对方程则方程也可改写成形如(2.5)的方程对齐次方程（2.5）利用变量替换可化为变量分离方程再求解.（2.6）（2.7）将（2.6）、（2.7）代入（2.5），则原方程变为整理后，得到（2.8）方程（2.8）是一个可分离变量方程，按照变量分离法求解，然后将所求的解代回原变量，所得的解便是原方程（2.5）的解.例4解即（2.9）分离变量，即有两边积分，得到（2.10）此外，方程（2.9如果（2.102.10 2.9）的通解为（2.10）.代回原来的变量，得到原方程的通解为例5解将方程改写为（2.11）分离变量，得到两边积分，得到（2.11）的通解即.此外，（2.11注意，此解不包括在通解（2.12）中.代回原来的变量，即得原方程的通解原方程的通解还可表为它定义于整个负半轴上.注：1.2..可分离方程.而这一齐次方程通过变量替换任然可化为可分离方程，因而，一定要熟练掌握可分离方程的解法.2）形如（2.13）.分三种情况来讨论（1.这时方程（2.13）属齐次方程，有.（2.这是一变量分离方程.（3.这时方程（2.13（2.14）1）（只需进行坐标（2.15）则（2.14）化为从而（2.13）变为（2.16）因此，得到这种情形求解的一般步骤如下：(1)解联立代数方程（2.14(2)作变换（2.15）将方程化为齐次方程（2.16）；(3)2.16）化为变量分离方程；(4)求解上述变量分离方程，最后代回原变量可得原方程（2.13）的解.上述解题的方法和步骤也适用于比方程（2.13此外，诸如以及.例6求解方程（2.17）解解方程组代入方程（2.17），则有（2.18）再令即则（2.18）化为两边积分，得因此此外，易验证即也就是（2.18）的解.因此方程（2.17）的通解为.3、 应用举例例7 电容器的充电和放电如图（2.1电荷，电容两端的电压为零.1.解 对于充电过程，由闭合回路的基尔霍夫第二定理，（2.19）（2.20） 将（2.20）代入（2.19（2.21）.方程（2.21）属于变量分离方程.将（2.21）分离变量，得到两边积分，得到即.所以（2.22）.由（2.22）知道，电压已达到外加电压的95%..易见对于放电过程的讨论，可以类似地进行.例8探照灯反射镜面的形状在制造探照灯的反射镜面时，总是要求将点光源射出的光线平行地射出去，以保证照灯有良好的方向性，试求反射镜面的几何形状.解（2.23）,仅,从而注意到（2.24）这是齐次方程.由2.12.再经直接积分即可求得方程的解.对于方齐次方程（2.242.24）得到于是（2.25）积分（2.25）并代回原来变量，经化简整理，最后得（2.26）.（2.26）就是所求的平面曲线，它是抛物线，因此，反射镜面的形状为旋转抛物面（2.27）小结: 本节我们主要讨论了一阶可分离微分方程和齐次微分方程的求解问题.将各种类型的求解步骤记清楚的同时要注意对解的讨论.§2 线性方程与常数变易法1、一阶线性微分方程（2.28）..2.28）变为（2.3）称为一阶齐线性方程.2.28）称为一阶非齐线性方程.2、常数变易法（2.3）是变量分离方程，已在例3中求得它的通解为（2.4）.下面讨论一阶非齐线性方程（2.28）的求解方法.方程(2.3)与方程(2.28)两者既有联系又有区别,设想它们的解也有一定的联系,在(2.4),它不可能是(2.28)的解,要使(2.28)具有形如(2.4)的解是常数,为此令（2.29）两边微分，得到（2.30）将（2.29）、（2.30）代入（2.28），得到即积分后得到（2.31） ..将（2.31）代入（2.29），得到（2.32）这就是方程（2.28）的通解.这种将常数变易为待定函数的方法，通常称为常数变易法.实际上常数变易法也是一种变量变换的方法.通过变换（2.29）可将方程（2.28）化为变量分离方程.注: 非齐线性方程的通解是它对应的齐线性方程的通解与它的某个特解之和. 例1. 解 将方程改写为（2.33）先求对应的齐次方程的通解，得令（2.34） 微分之，得到（2.35） 以（2.34）、（2.35）代入（2.33），再积分，得将其代入公式（2.34），即得原方程的通解.例2 .解原方程改写为（2.36）2.36）就是一个线性方程了.先求齐线性方程的通解为（2.37）于是代入（2.36），得到从而，原方程的通解为ln)y,.特别的，初值问题的解为例3 试证 （1）一阶非齐线性方程（2.28）的任两解之差必为相应的齐线性方程（2.3）之解；（22.28）的解，则（2.28）.（3）方程（2.3）任一解的常数倍或两解之和（或差）仍是方程（2.3）的解. 证 （1（1）—（2）有.（2）因为故结论成立.（3故结论成立.3、Bernoulli方程形如（（2.38）.利用变量变换可将伯努利方程化为线性方程来求解. 2.38）两边，得到（2.39）引入变量变换（2.40）从而（2.41）将（2.40）、2.41）代入（2.39），得到（2.42）（2.38）的通解.例4解得代入原方程得到这是线性方程，求得它的通解为或者这是原方程的通解.例5解将方程改写为.解法同上.例6便是伯努利方程.小结;这次主要讨论了一阶线性微分方程的解法.其核心思想是常数变易法.即将非齐线性方程对应的齐线性方程解的常数变易为待定函数，使其变易后的解函数代入非.我们还讨论了伯努利方程，.§3 恰当方程与积分因子1、恰当方程的定义将一阶微分方程写成微分的形式（2.43）.使得即则称方程(2.43)为恰当方程,或称全微分方程.在上述情形,方程(2.43)于是就是方程(2.43)的隐式通解,(应使函数有意义).2、恰当方程的判定准则定理1,则方程(2.43)是恰当方程的充要条件是而且当(2.46)成立时,相应的原函数可取为或者也可取为.证明 先证必要性.因为(2.43)是恰当方程,(2.45),,从而有下面证明定理的充分性,即由条件(2.46),使其适合方程(2.45).从(2.47)可知即(2.45)成立,同理也可从(2.48)推出(2.45). 例1. 解方程(2.49)解所以(2.49)是恰当方程.2.3,可按任意一条途代入公式(2.47)有代入公式(2.47)有都有3、恰当方程的解法上述定理已给出恰当方程的解法,下面给出恰当方程的另两种常用解法.解法1.已经验证方程为恰当方程,,有(2.50),有于是有因而省略了积分常数.把它代入(2.50)便得方程的通解为解法2.分项组合的方法对(2.49)式重新组合变为于是从而得到方程的通解为4、积分因子的定义及判别对于微分形式的微分方程（2.43）如果方程（2.432.43）的积分因子.(2.51)的通解，因而也就是（2.43）的通解.,则由恰当方程的判别准则知道即5、积分因子的求法方程(2.52)定理2的，则方程（（2.53）2.53（2.54）就是方程（2.43）的积分因子.证明 因为如果方程（2.43 2.52）进一步知即于是，有从而反之，如果（2.532.54）是方程（2.52）的解.事实上，因为因此函数（2.54）的确是方程（2.43）的积分因子.为了方便应用这个定理，我们就若干特殊情形列简表如下：例2.解注意所以方程不是恰当的,但是从而可得到隐式通解例3.解.但是从而可得到隐式通解它是用积分因子乘方程时丢失的解.例4.解方程.设想方程有积分因子.于是由上述简表知原方程有积分因子从而容易求得其通解为：六、积分因子的其他求法以例4为例，方程的积分因子也可以这样来求：把原方程改写为如下两组和的形式：设想原方程有积分因子..例5. 解方程.并且通解为例6、试用积分因子法解线性方程（2.28）.解将（2.28）改写为微分方程（2.55）方程（2.55（2.56）（2.56）为恰当方程，又分项分组法因此方程的通解为即与前面所求得的结果一样.注：积分因子一般不容易求得可以先从求特殊形状的积分因子开始，或者通过观察法进行“分项分组”法求得积分因子.§4 一阶隐方程与参数表示1、一阶隐方程一阶隐式微分方程的一般形式可表示为:则可化为显式形式,根据前面的知识求解.例如方程而其形式比较复杂,则宜采用引进参数的方法求解.一般隐式方程分为以下四种类型:2、求解方法1)讨论形如（2.57）,则方程(2.57)变为将(2.58) ,得到方程(2.59),而且属于显式形式.若求得(2.59)将其代入(2.58),于是得到(2.57)通解为若求得(2.59)于是得到(2.57)的参数形式的通解为.若求得(2.59)于是得到(2.57)的参数形式的通解为.例1 的解解,得到即,从而由此可知得到将其代入(2.60),即得 故参数形式的通解为,由(2.60).例2. 解,得到或(2.63)特解(2.63)与通解(2.62)中的每一条积分曲线均相切,因此称为方程的奇解.2)讨论形如的方程的求解方法,方程(2.64)与方程(2.57)的求解方法完全类似,有连续偏导数.则(2.64) 变为将(2.65) ,得到方程(2.66)),而且属于显式形式.设其通解为则(2.64)的通解为3) 讨论形如的方程的解法.,设曲线用参数形进而两边积分,得到于是得到方程(2.67)参数形式的解为.例3解则由方程得于是积分得到故原方程参数形式的通解为：4) 讨论形如的方程,其解法与方程(2.67)的求解方法类似.,设曲线用参数,有故方程(2.69)的参数形式的通解. 此外,不难验证,.例4 求解方程解则有由此可以得积分,故原方程参数形式的通解为.此外,. 例5解则故原方程参数形式的通解为.友情提示：方案范本是经验性极强的领域，本范文无法思考和涵盖全面，供参考!最好找专业人士起草或审核后使用。
《常微分方程》作业参考答案一．求解下列方程1．x c y cos =2．通解为：x x c y sin cos +=3．dx x x dy 122-= ⎰⎰--=122)1(xx d dy 2ln 1y x c =-+ 1)0(==c y 2ln |1|1y x ∴=-+4．'(1)ln(1)y yyy x x x -=++ 令 xuy x yu =⇔= (1)ln(1)dyduu x u u u dx dx ∴=+=+++故 (1)ln(1)dux u u dx =++(1)ln(1)du dx u u x =++ ln(1)ln(1)d u dxu x +=+ln ln(1)ln ln u x c ∴+=+ ln(1)u cx +=cx e u =+1 cx e x y=+∴1 )1(-=cx e x y5. 可分离变量方程，通解为.)1)(1(222cx y x =++6．齐次方程，通解为 c x x yx y =++ln 422sin .7．全微分方程，通解为 .64224c y y x x =+-8．.0222=++y dx dyx dx y d9. 解为 .)3(3x x y -=10. 通解为 .2sin 222c y x y x =++11．方程为 .011222=+-y x dx dyx dx y d12．通解为 ).tan(21c x c y +=二．1．通解为：c e e x y +=2212. 通解为： t t e c c e c z y x 2321123101210⎪⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛-3.0)0(0==y y 2121x y =52220121x x y += 4. x uN y uM ∂∂=∂∂ xu N x N u y u M y M u ∂∂+∂∂=∂∂+∂∂ 令 u y x =+22 y u d u d y u 2⋅=∂∂∴ x ud u d x u 2⋅=∂∂ u d u d x x N u u d u d y y M u 22+∂∂=+∂∂ ⎪⎪⎭⎫ ⎝⎛∂∂-∂∂=-∴y M x N u u d u d x y )(2故满定充要条件的表达式为：)(22y x xy y M xN +=--∂∂∂∂ϕ 5.)(2122y x v +=)(*dtdv)(22s x +-≤∠0 022≠+s x ∴（0.0）渐近稳定 6.一次近似方程为：⎪⎩⎪⎨⎧+=--=y x dtdy y x dt dx 32 特征方程为：012=++λλ 3-=∴∆＜0 P ＝1＞0 ∴0)Re(0)Re(21<<λλ, 则（0.0）局部渐过稳定. 7.01032=--λλ 5,221=-=λλx B x B x A x A y o 2sin )(2cos )(101*1+++=为x x y y y 2cos 10'3"=-- 之特解,±2λ不是特征根5=a 是特征方程的单根 x o e c x c x c x y 52122)(++=∴*故其通解为： 215221y y e c ec y x x +++=-8.特征根为：2.1.1321==-=λλλ 11-=λ所属的特征向量为：⎪⎪⎪⎭⎫ ⎝⎛-=532α12=λ所属的特征向量为：⎪⎪⎪⎭⎫ ⎝⎛=111β13=λ所属的特征向量为：γ⎪⎪⎪⎭⎫ ⎝⎛=101通解为：t t t e c e c e c z y x 2321101111531⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛-9.0:)0(=o y y 2121x y =52220121x x y -= 10.特征方程为：01072=++λλ07>=p 010>=g 0>∆故 (0.0)为稳定结点11.1.一次近似方程为：⎪⎩⎪⎨⎧-=--=yx y x t d y d dt x d 0222=++∴λλ0)Re(1<λ 0)Re(2<λ ∴（0.0）为局部渐近稳定 2.)(2122y x v +=. )1)((2222)(-++=*y x y x l dt dv 故122<+y x 0<∴dtdv 故（0.0）局部渐近稳定. 12. 1.,00=y ,31),(3020001x dx x dx y x f y y x x==+=⎰⎰ .63131)91(),(730620102x x dx x x dx y x f y y x x+=+=+=⎰⎰ 2. ,),(22y x y x f += ∴ ,5),(max ),(==∈y x f M Dy x ,42max max ),(),(L y y f D y x D y x ===∂∂∈∈ .5252,1min ,min =⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=m b a h则 .7564)52(32145)()(322=⋅⋅⋅≤-x y x y 13. 系数阵为 ,110111110⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡- 特征方程为 .0)1()det(2=--=-λλλE A E A λ-的初等因子为 2)1(,-λλ，通解为.101010101112321t t e t c e c c z y x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛14.证：设 [).),0()(..,0+∞∈∀≤>∃x M x f t s M .则[)+∞∈∀,0x ，有 .)1()(0)(0000M y e M y ds e Me y x y x x xx s x+≤-+=+≤--⎰[]),,0()(0x C x y ∈ ∴ [].,0,)(..,00x x M x y t s M ∈≤>∃令 {},,max 0M y M K += ∴ [).,0,)(+∞∈∀≤x K x y15.通解为 .)21(221xx e x x x c e c y -++=16．,2=α 特解为 ,1x y = 通解为 ).ln 21(221x x x c x c y +-+=。
习题2.31、验证下列方程是恰当方程，并求出方程的解。
1. 0)2()(2=-++dy y x dx y x 解： 1=∂∂yM，x N ∂∂=1 . 则xNy M ∂∂=∂∂ 所以此方程是恰当方程。
凑微分，0)(22=++-xdy ydx ydy dx x 得 ：C y xy x =-+23312． 0)4()3(2=---dy x y dx x y解： 1=∂∂yM，1=∂∂x N . 则xNy M ∂∂=∂∂ . 所以此方程为恰当方程。
凑微分，0432=--+ydy dx x xdy ydx 得 C y xy x =+-2323． 0])(1[]1)([2222=--+--dy y x x y dx xy x y解： 3422)(2)()1)((2)(2y x xyy x y x y y x y y M -=-----=∂∂ 3422)(2)()(2)(2y x xyy x y x x y x x x N -=-----=∂∂ 则yNx M ∂∂=∂∂ .因此此方程是恰当方程。
x y x y x u 1)(22--=∂∂ （1） 22)(1y x x y y u --=∂∂ （2） 对（1）做x 的积分，则)(1)(22y dx x dx y x y u ϕ+--=⎰⎰ =---yx y 2)(ln y x ϕ+ （3） 对（3）做y 的积分，则dy y d y x y y x y y u )()(2)()1(22ϕ+--+---=∂∂ =dy y d y x y xy )()(222ϕ+-+- =22)(1y x x y -- 则11)(21)(2)(1)(2222222-=-+--=-----=y y x y xy x y y x xy y y x x y dy y d ϕ y y dy yy -=-=⎰ln )11()(ϕyx xyx y y x y xy y x y y y x y x y u --=--+-=-+---=ln ln ln ln 222 故此方程的通解为C yx xyx y =-+ln 4、 0)2(3)23(22232=+++dy y y x dx x xy解：xy yM12=∂∂，xy x N 12=∂∂ . xNy M ∂∂=∂∂ . 则此方程为恰当方程。
