微分方程第3章习题解
- 格式:pdf
- 大小:195.81 KB
- 文档页数:12
45 / 2045习 题3-1 复摆重P ,对质心的回转半径为C ρ,质心距转动轴的距离为a ,复摆由水平位置无初速地释放,列写复摆的运动微分方程。
解:系统具有一个自由度,选复摆转角ϕ为广义坐标,原点及正方向如如题4-1图所示。
复摆在任意位置下,根据刚体绕定轴转动微分方程O O M J =ϕ其中)(22a gP J C O +=ρ 得到复摆运动微分方程为ϕϕρcos )(22Pa a gP C =+ 或0cos )(22=-+ϕϕρga a C3-2均质半圆柱体,质心为C ,与圆心O 1的距离为e ,柱体半径为固定R ,质量为m ,对质心的回转半径为C ρ,在平面上作无滑动滚动,如题3-2图所示,列写该系统的运动微分方程。
解:系统具有一个自由度,选θ为广义坐标。
半圆柱体在任意位置的动能为:222121ωC C J mv T +=题3-1图题3-2图46 / 2046用瞬心法求C v :2222*2)cos 2()(θθθ Re R e CC v C -+== θω = 2CC m J ρ= 故2222221)cos 2(21θρθθ Cm Re R e m T +-+=系统具有理想约束,重力的元功为 θθδd mge W sin -=应用动能定理的微分形式W dT δ=θθθρθθd mge m Re R e m d C sin 21)cos 2(2122222-=⎥⎦⎤⎢⎣⎡+-+ θθθθθθθθθθρd mge d mRe d mRe d R e m C sin sin cos 2)(2222-=+-++ 等式两边同除dt ,θθθθθθθθθθρ sin sin cos 2)(2222mge mRe mRe R e m C -=+-++ 0≠θ ,等式两边同除θ故微分方程为0sin sin )cos 2(2222=+++-+θθθθρθmge mRe Re R e m C ① 若为小摆动θθ≈sin ,1cos ≈θ,并略去二阶以上微量,上述非线性微分方程可线性化,系统微摆动的微分方程为0])[(22=++-θθρge r R C要点及讨论(1)本题也可以用平面运动微分方程求解。
习题3.11. 试用变量分离法求下列一阶微分方程的解. (1);dy x dx y=- 解: 分离变量得ydy xdx =-,两边积分得原方程的通解为222.x y C += (2)2cos .dyy x dx= 解: 分离变量得21cos dy xdx y =,两边积分得原方程的通解为1sin x C y-=+.0y =也是原方程的解. (3)2;dyxy dx= 解: 分离变量得12dy xdx y=,两边积分得原方程的通解为2ln y x C =+或2.x y Ce = (4) 22(1)(1);xy x dy y dx +=+ 解: 分离变量得2211(1)y dy dx y x x =++,即221()11y xdy dx y x x =-++. 两边积分得2211ln(1)ln ln(1)ln 22y x x C +=-++,通解为222(1)(1).x y Cx ++= (5)24;y x dye dx-= 解:分离变量得24yx e dy e dx --=,积分得241124y x e e C ---=-+,通解为422x y e e C ---=.(6)22;1dy dx x =- 解: 分离变量得221dy dx x =-,积分得微分方程的通解为1ln .1x y C x -=++ (7) 21;2dy y dx -= 解: 分离变量得221dy dx y =-,积分得原方程的通解为1ln 1y C y -=+.另外,1y =±也是解.(8)2cos .x dye y dx= 解: 分离变量得2sec x ydy e dx =,积分得原方程的通解为tan .x y e C =+另外,2y k ππ=+也是解.2. 作适当的变量变换求解下列方程. (1)2();dyx y dx=+ 解:令u x y =+,原方程变形为21du u dx -=,分离变量得211du dx u=+,积分得 arctan u x C =+,原方程的通解为tan().x y x C +=+(2)21;()dy dx x y =+ 解: 令u x y =+,原方程变形为211du dx u -=,分离变量得221u du dx u =+, 积分得arctan u u x C -=+,原方程的通解为arctan()y x y C =++.(3)21;21dy x y dx x y -+=-+ 解: 令210210x y x y -+=⎧⎨-+=⎩得11,.33x y =-=作代换11,33x X y Y =-=+,原方程变为齐次方程22dY X Y dX X Y -=-,再令Y u X =,该齐次方程变为212du uu X dX u-+=-, 分离变量得2121121u dX u u X -=-+,两端积分得21ln(1)ln ln 2u u X C -+=+,原方程的通解为2222(31)(31)(31)(31)(31).y y x x Cx x ---+++=+ (4)5;2dy x y dx x y -+=-- 解:令u x y =-,原方程变形为512du u dx u +-=-,分离变量得(2)7u du dx -=-,原方程的通解为2(2)14x y x C --=-+.(5)22(1)(41)81;dyx y xy dx=+++++ 解: 原方程即2(41)2dy x y dx =+++,作代换,令41u x y =++,方程变为21(1)24duu dx-=+,分离变量得2149du dx u =+,原方程的通解为2tan(6)(41).3x C x y +=++ (6) 6252222dy y x dx xy x y-=+; 解: 原方程即3322321()()232d y y x dx xy x -=+,令3y u =,方程变为齐次方程2221232du u x dx xu x -=+ 再令u v x =,后一方程又变为22116v dv dx v v x+=--,积分得 22ln(6)[ln(3)ln(2)]ln ln 5v v v v x c --+--+=+整理并代换变量得原方程的解散为:373315(3)(2)y x y x Cx -+=.(7) 322323.32dy x xy x dx x y y y++=+- 解:原方程即2222(231)(321)dy x x y dx y x y ++=+-,亦即222222()231()321d y x y d x x y ++=+- (1) 令22,u x v y ==,(1)式可变为231321dv u v du u v ++=+- (2) 作代换11,55u v ξη=-=-,(2)式变为2332d d ηξηξξη+=+ (3)作代换z ηξ=,(3)式变为2332dz z z d zξξ++=+,分离变量得22321z dz d z ξξ+=-- (4) (4)式两端积分得231ln(1)ln 2ln ln 21z z C z ξ--+=-++,整理并代回变量得原方程的通解为 22522(2)().y x C x y -+=+3. 已知0()()1(0)xf x f t dt x =≠⎰,试求函数()f x 的一般表达式.解:原方程变形为1()()xf t dt f x =⎰,两端求导得2()()()f x f x f x '=-,并由已知式子可知lim ()x f x →=∞。
3-1 设系统的微分方程式如下:(1) )(2)(2.0t r t c= (2) )()()(24.0)(04.0t r t c t c t c=++ 试求系统闭环传递函数Φ(s),以及系统的单位脉冲响应g(t)和单位阶跃响应c(t)。
已知全部初始条件为零。
解:(1) 因为)(2)(2.0s R s sC =闭环传递函数ss R s C s 10)()()(==Φ 单位脉冲响应:s s C /10)(= 010)(≥=t t g单位阶跃响应c(t) 2/10)(s s C = 010)(≥=t t t c(2))()()124.004.0(2s R s C s s =++ 124.004.0)()(2++=s s s R s C `闭环传递函数124.004.01)()()(2++==s s s R s C s φ单位脉冲响应:124.004.01)(2++=s s s C t e t g t 4sin 325)(3-= 单位阶跃响应h(t) 16)3(61]16)3[(25)(22+++-=++=s s s s s s C t e t e t c t t 4sin 434cos 1)(33----=3-2 温度计的传递函数为11+Ts ,用其测量容器内的水温,1min 才能显示出该温度的98%的数值。
若加热容器使水温按10ºC/min 的速度匀速上升,问温度计的稳态指示误差有多大解法一 依题意,温度计闭环传递函数11)(+=ΦTs s 由一阶系统阶跃响应特性可知:o o T c 98)4(=,因此有 min 14=T ,得出 min 25.0=T 。
视温度计为单位反馈系统,则开环传递函数为Tss s s G 1)(1)()(=Φ-Φ=⎩⎨⎧==11v TK !用静态误差系数法,当t t r ⋅=10)( 时,C T Ke ss ︒===5.21010。
解法二 依题意,系统误差定义为 )()()(t c t r t e -=,应有 1111)()(1)()()(+=+-=-==ΦTs TsTs s R s C s R s E s e C T sTs Ts ss R s s e s e s ss ︒==⋅+=Φ=→→5.210101lim )()(lim 203-3 已知二阶系统的单位阶跃响应为)1.536.1sin(5.1210)(2.1o tt et c +-=-试求系统的超调量σ%、峰值时间tp 和调节时间ts 。
微分方程习题§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.求下列方程通解。
§4.直纹面和可展曲面1. 证明曲面r =}32,2,31{2432v u u uv u v u +++是可展曲面.证法一: 已知曲面方程可改写为r =},2,{432u u u +v }32,,31{2u u ,令()a u =},2,{432u u u ,()b u =}32,,31{2u u ,则r =()a u + v ()b u ,且()b u ≠0,这是直纹面的方程 ,它满足(',,')a b b =23226412334013u u u u u u =0 ,所以所给曲面为可展曲面。
证法二:证明曲面的高斯曲率为零。
(略)2。
证明曲面r={cosv-(u+v)sinv, sinv+(u+v)cosv,u+2v}是可展曲面。
证法一: 曲面的方程可改写为 r=()a v + u ()b v ,其中()a v ={cosv-vsinv,sinv+vcosv, 2v},()b v ={-sinv, cosv,1} ,易见()b v ≠0,所以曲面为直纹面,又因为(',,')a b b =2sin cos 2cos sin 2sin cos 1cos sin 0v v v v v v v v vv ------=0,所以所给曲面为可展曲面。
证法二:证明曲面的高斯曲率为零。
(略)3.证明正螺面r={vcosu,vsinu,au+b}(a ≠0)不是可展曲面。
证法一:原曲面的方程可改写为r=()a u + v ()b u ,其中()a u ={0,0,au+b},()b u ={cosu,sinu,0}.易见()b u ≠0, 所以曲面为直纹面,又因为(',,')a b b =00cos sin 0sin cos 0au u u u -=a ≠0.故正螺面不是可展曲面。
证法二:证明曲面的高斯曲率为零。
(略)4.证明挠曲线的主法线曲面与副法线曲面不是可展曲面。
马知恩周义仓编常微分⽅程定性与稳定性⽅法部分习题参考解答第⼀章 基本定理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$$ ⽽零解是渐近稳定的.。
第三章习题简答3-1 已知流体流动的速度分布为22y x u x -= ,xy u y 2-=,求通过1,1==y x 的一条流线。
解:由流线微分方程yx u dyu dx =得dy u dx u x y =则有 dy y x xydx )(222-=-两边积分可得C y y x yx +-=-3322即0623=+-C y x y将x=1,y=1代入上式,可得C=5,则 流线方程为05623=+-y x y3-3 已知流体的速度分布为⎭⎬⎫==-=-=tx x u ty y u y x 00εωεω(ω>0,0ε>0)试求流线方程,并画流线图。
解:由流线微分方程yx u dyu dx =得dy u dx u x y =则有 tydy txdx 00εε-=两边积分可得C y x +-=22流线方程为C y x =+223-5 以平均速度s m v /5.1=流入直径为D=2cm 的排孔管中的液体,全部经8个直径d=1mm 的排孔流出,假定每孔出流速度依次降低2%,试求第一孔与第八孔的出流速度各为多少?题3-5图解:由题意得:v 2=v 1(1-2%),v 3=v 1(1-2%)2,…,v 8=v 1(1-2%)7 根据质量守恒定律可得282322212832144444dv d v d v d v D v Q Q Q Q Q πππππ⋅+⋅⋅⋅+⋅+⋅+⋅=⋅+⋅⋅⋅+++=sm d vD v v d v v v v d D v /4.80)98.01(001.002.002.05.1)98.01()98.01(98.01)98.01(4)(448228221812832122=-⨯⨯⨯=--⋅=∴--⋅=+⋅⋅⋅+++⋅=⋅πππ则 v 8=v 1(1-2%)7=80.4×(1-2%)7=69.8m/s3-6 油从铅直圆管向下流出。
管直径cm d 101=,管口处的速度为s m v /4.11=,试求管口处下方H=1.5m 处的速度和油柱直径。
第三章 一阶微分方程的解的存在定理3-1 求下列初值问题的近似解。
1) 求初值问题⎪⎩⎪⎨⎧=+=0)0(2y yx dx dy 的第三次近似解;2) 求初值问题⎪⎩⎪⎨⎧=-=0)1(2y yx dx dy 的第二次近似解。
解 由解的存在唯一性定理知,1),2)中的初值问题的解分别在)0,0(,)0,1(的邻域内存在且唯一。
下面求它们的近似解。
1) 0)(0=x φ,2)(201x xdx x φx==⎰,220]4[)(25042x x dx x x x φx+=+=⎰,⎰⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++=xdx x x x x φ02253220)(440016*********x x x x +++=。
2) 0)(0=x φ,212)(211-==⎰x xdx x φx,3011412620212)(2351222--++-=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛--=⎰x x x x dx x x x φx。
评注:逐次逼近函数序列00)(y x y =,⎰-+=xx n n dx x y x f y x y 0))(,()(10,在实际中有广泛的应用。
利用此序列求近似解时,须验证初值问题的解存在唯一,否则求出的结果可能并不是我们想要的近似解。
3-2 设1,11:),(≤≤+∈y x D y x ,求初值问题⎪⎩⎪⎨⎧=--=0)1(22y yx dxdy 的解的存在区间,并求第二次近似解,给出在解的存在区间的误差估计。
解 设22),(y x y x f -=,显然,方程在D 上满足解的存在唯一性定理,则1,1,4),(max ),(====∈b a y x f M Dy x ,所以41)41,1min(),min(===M b a h , 方程过点)0,1(-的解的存在区间为:411≤+x ,即4345-≤≤-x 。
设)(x ϕ是初值问题 0)1(22⎪⎩⎪⎨⎧=--=y yx dx dy 的解,)(2x ϕ是第二次近似解,则,0)(0=x ϕ313)(3121+==⎰-x dx x x xϕ,4211931863])313([)(34712322+-+--=+-=⎰-x x x x dx x x x xϕ。
3t15t=e ( e +c)5=c e 3t +15e 2t 是原方程的解ds 13. =-s cost + sin2tdt 2cos tdt 13dt解:s=e ( sin2t e dt c )=esint( sin t coste sin t dt c) sin tsint sint= e( sin tee c )常微分方程 习题 2.2求下列方程的解1. dy = y sin x dx解: y=e ( sinxe dx c)x1 x=e x [- e x (sinx cos x )+c]= ce sint sint 1 是原方程的解。
4.dy xy e x x n,n 为常数. dx n解:原方程可化为:dy xy e x x n dx n方程的解。
=c e(sinx cos x )是原yendxx x ( e x x e n n dx n xdx c)2.dx+3x=e 2tnxx (ec)dt 解:原方程可化为:dx=-3x+edt是原方程的解 .所以:3dtx=ee2te 3dt5.dy +1 22x y 1=0 dx x 2dt c)ds23P(x) ,Q(x) (x 1)3 x1P(x)dxee=(x+1) 2((x 21) c)即: 2y=c(x+21+)(x+14) 为方程的通解。
8.d dy x =x y y 33dx x+y 1 2 解: xy 2dy y yP(y)dy P(y)dy( e Q(y)dy c) =y( 1*y 2dy c)y3= y cy23即 x=y +cy 是方程的通解 ,且 y=0也是方程的解。
2解:原方程可化为:dy dx1x 22xy 1x7.dy 2y (x 1)3dx x 1 解:dy 2y(x 1)3 dx x 1 (x 1)2(ln x 2e方程的通解为:ln x 2 1( e x dx c)1= x 2(1 ce x )P(x)dx P(x)dxy=e ( e Q(x)dx c) =(x+1)(=(x+1)((x 11)2 *(x+1)3dx+c) (x+1)dx+c) 是原方程的解.x=edx c )2则P(y)=y 1,Q(y) y 2方程的通解9. dy ay x 1,a 为常数 dx x x解:(P x) a ,Q(x) x 1xP(x )dxeedx方程的通解为:y=(x)dx P (x)dx(e Q(x)dx=xa(1 x+1dx+c)x a时,x 方程的通解为11.dy xy x 3y 3 dx 解:dy xyx 3y 3dx 两边除以3y c)d 3y xy 2 x 3 ydxdy2( xy 2 x 3)y=x+ln/x/+c当 y=cx+xln/x/-1当 a 1时, 方程 的通解为a 0,1时,方程的通解为y=cxa x 1 +-1- a adx 令y 2 z dz 2( xz x 3) dx P(x) 2x,Q(x) 2x 3 epx dx e2xdxe x 2 方程的通解为:z= e dx( e dxQ(x)dx c)10.x d d x y y x 3解:d dy x 1x y x 3P(x) 1,Q(x) x =e =xx(e x (2x 3)dx c) 22ce x1故方程的通解为y :2(x 2 ce x 1) 1,且y 0也是方程的解。