微分方程第3章习题解

合集下载

第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)本题也可以用平面运动微分方程求解。

自动控制原理第三章习题解答

自动控制原理第三章习题解答
σ % = e −πξ /
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[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章

习题答案第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]

微分几何(第三版)【梅向明_黄敬之_编】第三章课后题答案[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. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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) 为 在 区 间
相关文档
最新文档