微分方程第3章习题解
第3章 振动系统的运动微分方程题解
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)本题也可以用平面运动微分方程求解。
自动控制原理第三章习题解答
tp =
1−ξ 2
= e −π 0.6 /
1−0.62
= e −π 0.6 /
1−0.62
= 9 .5 %
π
1 − ξ ωn
2
=
π
1.6
= 1.96( s )
ts =
3-5
3.5
ξω n
=
3.5 = 2.92( s ) 1.2
设单位反馈系统的开环传递函数为
G ( s) =
0.4 s + 1 s ( s + 0.6)
s5 s4 s3 s2 s1 s0
1 12 35 3 20 25 16 80 3 3 5 25 10 25
有一对虚根,系统不稳定 3-13 已知单位反馈系统的开环传递函数
G ( s) =
K (0.5s + 1) s ( s + 1)(0.5s 2 + s + 1)
试确定系统稳定时的 K 值范围。 解:系统特征方程为
ε 0 ,试问 k1 应满足什么条件?
见习题 3-20 解答 3-2 设系统的微分方程式如下: (1)
&(t ) = 2r (t ) 0.2c
&&(t ) + 0.24c &(t ) + c(t ) = r (t ) (2) 0.04c
试求系统的单位脉冲响应 k(t)和单位阶跃响应 h(t)。已知全部初始条件为零。 解: (1) 因为 0.2 sC ( s ) = 2 R ( s ) 单位脉冲响应: C ( s ) = 10 / s 单位阶跃响应 h(t)
试求系统的超调量σ%、峰值时间tp 和调节时间ts。 解: h(t ) = 1 −
偏微分方程总复习和课后习题答案
一、基本概念
1. 偏微分方程的定义P1 2. 偏微分方程的阶数,线性、拟线性、完全非线性 偏微分方程的定义P10 3. 偏微分方程的适定性P23
二、方程的导出,分类与化简
三、公式的直接应用题
1. 2. 3. 4. 5. 达朗贝尔公式P36 公式P42 傅里叶(逆)变换P106 P110例 4.1.7结论 泊松公式P112
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at x a ( t ) 1 t d f ( , )d x a ( t ) 2a 0
1 2 u ( x t ) 3t xt 2
1 1 xa t C f1 ( x at ) ( x at ) ( )d 2 2a x0 2 1 1 xa t C f 2 ( x at ) ( x at ) ( )d 2 2 a x0 2
1 1 xat u [ ( x at ) ( x at )] ( )d 2 2a x a t
1 u ( x t ) x (1 a )t cos x sin at a
2 2 2
1 ( 7)
解:
2
1 22 1 x at x at x u ( x t ) 5 x t a t 2 (e e 2e ) 3 2a
1 ( 6)
解:
2 2u u 2 1 a f ( x , t ), x R ,t 0 2 2 t x u ( x, 0) ( x), u ( x, 0) ( x), x R1. t
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at x a ( t ) 1 t d f ( , )d x a ( t ) 2a 0
常微分方程习题3[1].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.求下列方程通解。
习题答案第3章
t s (2%) 4T
T ,为惯性环节的时间常数。 将已知数据 t s (2%)=15 秒代入上式,求得惯性环节的时间常数 T
则闭环传递函数为
3.75 秒。
WB ( s )
单位反馈系统的开环传递函数为
1 15s 1
WK ( s )
WB ( s) 1 1 WB ( s ) 15s
s0
劳斯表中出现 s 行为全 0 行,且无符号变化,则闭环系统临界稳定,有 2 对对称于原 点的特征根。可通过辅助方程得到。
3
令 F ( s ) 3s 18s 12 0
4 2
解得
3-9 系统如图 P3-4 所示,问 取何值系统方能稳定。
s1, 2 j 0.87 , s 3, 4 j 2.29 10 s ( s 1)
则 令 xc (t ) 0
xc (t ) L1 [ X c ( s )] 1 e t cos(3t )
可得 t m 0.94 s
阶跃响应的最大峰值 根据超调量的定义
x max (t m ) 1.37
%
调节时间 t s (5%)
x max (t m ) xc () 100% 37% x c ( )
0.2 sX c ( s ) 2 X r ( s )
又输入信号为 X r ( s )
1 ,则输出 s
X c ( s ) 10 s2
拉氏反变换后,得单位阶跃响应为
xc (t ) 10t
c (t ) 0.24 x c (t ) (2) 0.04 x
微分方程两侧同时取拉氏变换,得
5 4 3 2
(4) s 4 s 4 s 4 s 7 s 8s 10 0
微分几何(第三版)【梅向明_黄敬之_编】第三章课后题答案[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.证明挠曲线的主法线曲面与副法线曲面不是可展曲面。
理论力学(周衍柏 第二版)第3章习题解答
∑F ∑F
∑M
到最小时,
y
x
= N 2 − f1 = 0 ①
= f 2 + N1 − G1 − G2 = 0 ②
且梯子沿过 A 点平行于 z 轴的合力矩为零。即:
i
= G2 l cos θ + G1
l cos θ − f 2 l cos θ − N 2 l sin θ = 0 ③ 2
又因梯子是一个刚体。当一端滑动时,另一端也滑动,所以当梯与地面的倾角达
对于 C 球,它相对于过 D 点与 z 轴平行的轴的合力矩等于零。即:
= Tr sin(β − α ) − Gr sin β = 0 ②
tan β = 3 tan α
3.5 解 如题 3.5.1 图。
y A
o
f2
N2
N1
G2 G 1
f1 B x
题3.5.1图
梯子受到地面和墙的弹力分别为 N1 , N 2 ,受地面和墙的摩擦力分别为 f1 , f 2 。 梯子和人的重力分别为 G1 ,G2 且 G2 = 3G1 。设梯长为 l ,与地面夹角为 θ 。由于 梯子处于平衡,所以
2
=1
可求该切面的面积
⎛ y2 ⎞ ⎟ S ( y ) = πac⎜ − 1 ⎜ b2 ⎟ ⎝ ⎠
故积分
b b ⎛ y2 ⎞ 4 2 2 2 ⎜ ⎟ ρdy = πρab3c y dm y S dy y ac 1 = ⋅ ρ = π − 2 ⎟ ∫ ∫−b ( y ) ∫−b ⎜ b ⎠ 15 ⎝
同理可求
o
αα
y
T
B Tα
β β
β −α c
A r
题3.4.1图
Ox 轴竖直向下,相同的球 A 、 B 、 C 互切, B 、 C 切于 D 点。设球的重力大小
马知恩周义仓编常微分方程定性与稳定性方法部分习题参考解答
马知恩周义仓编常微分⽅程定性与稳定性⽅法部分习题参考解答第⼀章 基本定理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$$ ⽽零解是渐近稳定的.。
5(3)一阶线性微分方程
P( x)dx
, 求导 得
dy + P ( x ) y = Q( x ) dx
P( x)dx
′ = C′( x)e ∫ y
P( x)dx
+C(x)e ∫
[P( x)]
dy 将 y和 y ′代入原方程 , 得 + P ( x ) y = Q( x ) dx
C ′( x )e ∫
P ( x ) dx
y = e ∫2
P( x)dx
∫ P( x)dxdx + C] [∫ Q( x)e
y=e
∫ x2 1
2x
dx
x cos x ∫ x2 1dx 2 dx + C ∫ 2 e x 1
1 (C + sin x ) = 2 x 1
由初始条件 y x = 0 = 1
C = 1,
1 sin x 特解 y = 1 x2
y 原方程的解 tan = Ce x + (1 x ) 2
21
一阶微分方程
四,小结
(1)一阶线性微分方程 一阶线性微分方程
dy + P( x) y = Q( x) dx
y=e ∫
P( x)dx
∫ P( x)dxdx + C] [∫ Q( x)e
( n ≠ 0,1)
dy (2)伯努利微分方程 + P ( x ) y = Q( x ) y n 伯努利微分方程 dx
( x y ) = ∫0 f ( x )dx
3 2
x
(∫
x
0
3 ′ ydx ) = (x y)′
2
y
Q
y = 3 x y′
积分方程
y = x3
《流体力学》徐正坦主编课后答案第三章
第三章习题简答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ϕ。
微分方程数值解法(戴嘉尊)习题解答
+
R Lh
(eL( X
− x0 )
−1)
电子文档制作:成都信息工程学院 数学学院 杨韧 吴世良,2010 年 4 月
成都信息工程学院>>精品课程>>微分方程数值解
11、解:令 f(x,y)=-y+x+1
y y y x y x y x = + h(− + +1) = (1− h) + h( +1) = 0.9 + + 0.1
0.0988*1.0e-3
0.9
0.4973
0.4972
0.0640*1.0e-3
1
0.5002
0.5000
0.1773*1.0e-3
2.解:显然, y = e−x 是原初值问题的准确解。 由梯形公式得
整理可得: 于是:
yn+1
=
yn
+
h 2
[
f
(
xn
,
yn
)
+
f
(xn+1, yn+1)]
=
yn
+
h 2
成都信息工程学院>>精品课程>>微分方程数值解
微分方程数值解 习题解答
杨韧 吴世良(编)
成都信息工程学院 数学学院
二 O 一 O 年四月编写
电子文档制作:成都信息工程学院 数学学院 杨韧 吴世良,2010 年 4 月
目
成都信息工程学院>>精品课程>>微分方程数值解
录
第一章 常微分方程数值解 ......................................................................3 第二章 抛物型方程的差分方法 ..............................................................8 第三章 椭圆型方程的差分方法 ............................................................16 第四章 双曲型方程的差分方法 ............................................................25
微分方程-3
xdy − ydx ⎛ y⎞ = d⎜ ⎟ 2 x ⎝ x⎠
xdy + ydx = d (ln xy ) xy
xdx + ydy ⎞ ⎛1 = d ⎜ ln( x 2 + y 2 ) ⎟ x2 + y2 ⎠ ⎝2
xdy − ydx ⎛1 x + y⎞ = d ⎜ ln ⎟ 2 2 x − y ⎝2 x− y⎠
∂ ( µP ) ∂ ( µQ ) 1.公式法: Q = , ∂y ∂x ∂P ∂µ ∂Q ∂µ 两边同除 µ, µ +P =µ +Q ∂y ∂y ∂x ∂x ∂ ln µ ∂ ln µ ∂P ∂Q −P = − Q ∂x ∂y ∂y ∂x
特殊地: 求解不容易
∂µ ∂µ dµ a . 当µ只与x有关时; = 0, = , ∂y ∂x dx
u( x , y ) = ∫0 ( x − 3 xy )d x + ∫0 y 3 dy
x 3 2 y
x4 3 2 2 y4 = − x y + , 4 2 4
x 3 2 2 y − x y + = C. 原方程的通解为 4 2 4
4
4
江西理工大学理学院
2x y2 − 3x2 例2 求方程 3 dx + dy = 0的通解 . 4 y y
dy x2 + x3 + y 的通解. 例6 求微分方程 = − 1+ x dx
解1
1 dy + y = − x2 , 整理得 dx 1 + x
C . A 常数变易法: 对应齐次方程通解为 y = 1+ x x3 x4 C ( x) 设 y= . C ( x) = − − + C . 3 4 1+ x
常微分方程第三章习题解答
例题选讲
例1 证明方程(3.2)的线性无关解的个数不超 过n+1.它存在n+1个线性无关解吗?
证明 如果方程(3.2)有n+2个线性无关解:
Y1( x),K,Yn+2 ( x)
则Y1( x) − Yn+2 ( x),Y2 ( x) − Yn+2 ( x),K,Yn+1( x) − Yn+2 ( x)是方程(3.1) 的n + 1个解.由齐次线性方程组解的 结构知,这n + 1个解必然线性 相关, 于是, 存在不全为零的常数 C1 ,K, Cn+1使得 C1[Y1( x) − Yn+2 ( x)] + C2[Y2 ( x) − Yn+2 ( x)] + K + Cn+1[Yn+1( x) − Yn+2 ( x)] = 0
其中A是n× n实矩阵.
(b). 如果λ是矩阵 A的k重单特征根 , 则(3.3)有形式为
P( x)eλx的解 , 其中, P( x)是次数不超过 k − 1的多项式 向量函数 .
(c). 如果λ是矩阵A的k重特征根,则方程(3.3)存在
k个线性无关解
Y = ( R0 + R1 x + L + Rk−1 x k−1 )eλx 其中, R0 , R1 ,K, Rk−1由下列方程确定
x
+
1
⎜⎛ ⎜
−
1 1
1 −1
2 ⎟⎞2 − 2 ⎟ x2]
⎜⎝ 0 0 1⎟⎠ ⎜⎝ 0 0 0⎟⎠ 2 ⎜⎝ 0 0 0⎟⎠
⎜⎛1 0 0⎟⎞ ⎜⎛ x x 2x ⎟⎞ = e x[⎜ 0 0 1⎟ + ⎜ − x − x − 2x ⎟]
微分方程习题及答案
(2);
(3);
(4).
2、求连续函数,使得时有。
3、求以为通解得二阶微分方程、
4。某个三阶常系数微分方程有两个解与,求。
5、设有一个解为,对应齐次方程有一特解,试求:
(1)得表达式;
(2)该微分方程得通解.
6、已知可导函数满足关系式:
求。
7.已知曲线上原点处得切线垂直于直线,且满足微分方程,求此曲线方程.
5、长为6m得链条自桌上无摩察地向下滑动,设运动开始时,链条自桌上垂下部分长为1m,问需多少时间链条全部滑过桌面。
§7二阶常系数非齐次线性微分方程
1。求下列微分方程得通解
(1);
(2);
(3);
(4);
(5).
2。求下列微分方程得特解
(1);
(2)
3.设连续函数满足求。
4、一质量为得质点由静止开始沉入水中,下沉时水得反作用力与速度成正比(比例系数为),求此物体之运动规律、
(1);
(2).
5、 用适当得变换替换化简方程,并求解下列方程
(1);
(2)
(3)
(4)
6.求一曲线,使其任意一点得切线与过切点平行于轴得直线与轴所围城三角形面积等于常数、
7、设质量为得物体自由下落,所受空气阻力与速度成正比,并设开始下落时速度为0,求物体速度与时间得函数关系、
8。有一种医疗手段,就是把示踪染色注射到胰脏里去,以检查其功能。正常胰脏每分钟吸收掉染色,现内科医生给某人注射了0、3g染色,30分钟后剩下0。1g,试求注射染色后分钟时正常胰脏中染色量随时间变化得规律,此人胰脏就是否正常?
5。一链条悬挂在一钉子上,起动时一端离开钉子8m,另一端离开钉子12m,若不计摩擦力,求链条全部滑下所需时间。
常微分方程第三版课后答案
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也是方程的解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3-6
假设函数
f (x, y) 及 ∂f ∂y
都在区域 G 内连续,又 y
= ϕ (x, x0 , y0 ) 是方程
y′ = f (x, y)
满足初始条件
y0
=
ϕ(x0 , x0 ,
y0 ) 的解,试证
∂ϕ ∂y0
存在且连续,并写出其表达式。
证 1)因 f (x, y) 及 ∂f 都在区域 G 内连续,则 f 在 G 内满足局部利普希兹条件,故解 ∂y
y
=
⎪⎧0 ⎨
(x ≤ C)
3
。
⎪⎩(x − C) 2 (x > C) , C ≥ 0
评注:寻找解的存在唯一性定理中的条件所满足的区域,就是寻找 f (x, y) 连续和关于 y 满足利普希兹条件的区域,困难在于利普希兹条件的验证,除用定义外,还常用下面的结
论:
∂f 在 D 上存在且有界,则 f (x, y) 在 D 上关于 y 满足利普希兹条件。 ∂y
ϕ ( x1 ) > 0 ,由 ϕ ( x) 的连续性及ϕ ( x0 ) = 0 ,知必存在 x 0 , x0 ≤x 0 ≤ x1 ,使得
ϕ (x0 ) = 0 及ϕ (x) > 0 , x 0 ≤ x ≤ x1,
则有
x
∫ ϕ i ( x ) = ϕ i ( x 0 ) + x0 f ( x,ϕ i ( x))dx , i = 1,2 。
解 设 f ( x, y) = x 2 − y 2 ,显然,方程在 D 上满足解的存在唯一性定理,则
M = max f (x, y) = 4, a = 1, b = 1, ( x,y )∈D
所以
b
11
h = min(a, ) = min(1, ) = ,
M
44
方程过点 (−1,0) 的解的存在区间为: x + 1 ≤ 1 ,即 − 5 ≤ x ≤ − 3 。
式
ϕn (x)
− ϕ (x)
≤
MLn h n+1 ,在进行近似计算时,可以根据误差的要求,选取适当的逐 (n + 1)!
步逼近函数ϕ n ( x) 。
3-3
dy
讨论方程
=
3
1
y 3 在怎样的区域中满足解的存在唯一性定理的条件,并求通过点
dx 2
(0,0) 的一切解。
解
设
f
( x,
y)
=
3
1
y3
,则
x
∫ 而
ϕ(x) = ϕ1(x) − ϕ2 (x) =
[f
x0
(x,ϕ1 (x)) −
f
(x,ϕ 2 ( x))]dx
,其中 x 0 ≤
x≤
x1 。
由ϕ(x) = ϕ1 (x) − ϕ2 (x) > 0 及 f (x, y) 对 y 的不增性,知
x
∫ ϕ1(x) − ϕ2 (x) =
[f
x0
(x,ϕ1(x)) −
当 ε → 0+ 时,有 f (t) ≤ 0 。
因为 f (t) ≥ 0 ,即得 f (t) ≡ 0 ,从而
t
f (t) ≤ K exp( ∫α g (s)ds ), α ≤ t ≤ β
综上所述,不等式成立。 唯一性的证明。
设ϕ(t),ψ (t) 是初值问题 x′ = f (t, x), x(t0 ) = x0 的两个解,则有
t
t
∫ ∫ ϕ(t) = x0 + t0 f (ξ ,ϕ(ξ ))dξ ,ψ (t ) = x 0 +
f (ξ , ϕ (ξ )) dξ 。
t0
于是
t
t
∫ ∫ ϕ (t) −ψ (t) ≤ f (ξ ,ϕ (ξ )) − f (ξ ,ψ (ξ )) dξ ≤ L ϕ (ξ ) −ψ (ξ ) dξ ,
y′ = f (x, y) 的满足初始条件 y(x0 ) = y0 的解 y = y(x, x0 , y0 ) 对一切 x ≥ x0 有定义,试证
下列说法是等价的:
(1) 任 给 ε > 0 , 可 以 找 到 正 数 δ = δ (ε , x0 ) , 使 当 y0 ≤ δ 时 , 对 一 切 x ≥0 和 x0 ≥ x0 ,存在 δ 2 > 0 ,使 y(x0 , x0 , y0 ) ≤ δ 2 ,对一切 x ≥ x0 ,
有 y(x, x0 , y0 ) < ε ,因为方程的解 y = y( x, x0 , y0 ) 在 G 内连续依赖于初值 (x0 , y 0 ) ,
第三章 一阶微分方程的解的存在定理
3-1 求下列初值问题的近似解。
⎧dy ⎪
=
x
+
y2
1) 求初值问题 ⎨ dx
的第三次近似解;
⎪⎩ y(0) = 0
⎧dy ⎪
=
x
−
y2
2) 求初值问题 ⎨ dx
的第二次近似解。
⎪⎩ y(1) = 0
解 由解的存在唯一性定理知,1),2)中的初值问题的解分别在 (0,0) , (1,0) 的邻域
x
∫ y = y0 +
f (x, y)dx
x0
。
即
x
x
∫ ∫ φ = y0 +
x0 f (x, φ)dx 和 ψ = y0 + ∆y0 +
2
∂f
=
1
y
2 −
3
(y
≠
0)
,
∂y 2
故在 y ≠ 0 的任何区域上 ∂f 存在且连续,因而方程在这样的区域中满足解的存在唯一 ∂y
性定理的条件。
显然, y ≡ 0 是通过点 (0,0) 的一个解;
又由方程
dy
=
3
1
y3
得
dx 2
3
y = (x − C)2 。
所以通过点 (0,0) 的一切解为 y ≡ 0及
即有
t
f (t) ≤ w(t) ≤ K exp(∫α g(s)ds), α ≤ t ≤ β 。
t
∫ 2) K = 0 时,对任意ε > 0 ,由于 f (t ) ≤ f (s) g (s)ds , 所以 α t f (t) ≤ ε + ∫α f (s) g(s)ds,
由 1)有
t
f (t) ≤ ε exp( ∫α g (s)ds ),
4
4
4
设ϕ(x) 是初值问题
⎧dy ⎪
=
x2
−
y2
⎨ dx
⎪⎩ y(−1) = 0
的解, ϕ2 (x) 是第二次近似解,则
ϕ0 (x) = 0,
∫ ϕ 1 ( x) =
x x 2 dx
=
x3
1 +
,
−1
33
∫ ϕ 2 ( x) =
x [x2
x3 −(
+
1 ) 2 ]dx
=
−
x7
−
x4
+
x3
−
x
+
11
。
−1
33
63 18 3 9 42
在区间
x
+1
≤
1 4
上, ϕ 2
(x) 与ϕ(x)
的误差为
ϕ2 (x) − ϕ(x)
≤
ML2 h3 (2 + 1)!
=
4 L2 3!
1 43
,
取 ∂f = − 2 y ≤ 2 = L , ∂y
4⋅22 1 1
所以
ϕ2 (x) − ϕ(x) ≤
=。 3! 43 24
评注:需要掌握第 n 次近似解ϕ n (x) 和真正解 ϕ (x) 在区间 x − x0 ≤ h 内的误差估计公
t0
t0
其中 L 为利普希兹常数,由上面的不等式可知
0 ≤ ϕ(t) −ψ (t) ≤ 0,
因而有ϕ(t) ≡ ψ (t) 。
评注:格朗瓦耳不等式是微分方程中的重要不等式,表明积分不等式与其解的关系。用 格朗瓦耳不等式证明微分方程初值问题解的唯一性是一个很好的方法。
3-5 假定函数 f ( x, y) 于 (x0 , y0 ) 的邻域内是 y 的不增函数,试证初值问题
内存在且唯一。下面求它们的近似解。
1) φ0 (x) = 0 ,
x
x2
∫ φ1 ( x) =
xdx =
0
2
,
x
x4
x5 x2
∫ φ2 (x) =
[x +
0
]dx = +
4
20
2
,
∫ φ3 (x) =
x⎡
⎛ x5
x2
2
⎞
⎤
0
⎢x ⎢⎣
+
⎜⎜ ⎝
20
+
2
⎟⎟ ⎠
⎥dx = ⎥⎦
x2 x5 x8
x 11
++ +
y = ϕ ( x, x0 , y0 ) 在它的存在范围内对 x, x0 , y0 连续。 2)设由初值 (x0 , y0 ) 和 (x0 , y0 + ∆y0 ), ∆y0 足够小,所确定的解分别为
y = ϕ(x, x0 , y0 )
和
y = ψ (x, x0 , y0 + ∆y0 ) ,
则这两个解均满足积分方程
∂f 在 D 上存在且无界,则 f (x, y) 在 D 上关于 y 不满足利普希兹条件。 ∂y
其中 D 为某矩形区域。 3-4 证 明 格 朗 瓦 耳 (Gronwall) 不 等 式 : 设 K 为 非 负 常 数 , f (t) 和 g(t) 为 在 区 间