三. 试用拉格朗日方程建立弹簧振子的运动微分方程,并求出其振动周期(已知:弹簧的倔强系数为K ,物块的质量为m )。
四. 长l 2,质量为m 的均匀棒,其上端A 靠在光滑的墙上,下端则固联一不能伸长的线BC ,线的上端固结于墙上C 点,C 点与A 点在同一垂直线上,棒与墙所成的角度为α,线与墙所成的角度为β,如果ABC 平面为与墙垂直的铅垂面。求平衡时αβ与之间的关系。(用刚体平衡方程求解)。
三. 解:系统自由度1=S ,取q=x,系统的动能2
'2
1x m T = 系统的势能221kx V = =-=V T L 2'2
1x m -221kx 代入拉氏方程:0)(=??-??x l x l dt d ,得: 0''=+kx x m 0''=+
∴x m k x 令m
k w =2,则w 为弹簧振子简谐振动的圆频率。 k
m W T ππ22==∴
四. 解:αβ
cos 0(1)0
sin 0(2)0sin 2cos 0(3)()0yi xi
B i T mg F N T F mgl N l m F ββαα?-==?-==??-==?
∑∑∑ )
1()2((3)N tg N mg tg mg ββ==得:。代入式得
sin 2cos 0mgl mg tg l αβα-?= 即:202tg tg tg tg αβαβ-=∴=
习题2.5 2.ydy x xdy ydx 2=- 。 解: 2x ,得: ydy x xdy ydx =-2 c y x y d +-=221 即c y x y =+2 2 1 4. xy x y dx dy -= 解:两边同除以x ,得 x y x y dx dy - =1 令u x y = 则dx du x u dx dy += 即 dx du x u dx dy +=u u -=1 得到 ()2ln 2 1 1y c u -=, 即2 ln 21?? ? ??-=y c y x 另外0=y 也是方程的解。 6.()01=-+xdy ydx xy 解:0=+-xydx xdy ydx x d x y x d y y d x -=-2 得到c x y x d +-=??? ? ??2 21
即 c x y x =+2 2 1 另外0=y 也是方程的解。 8. 32 x y x y dx dy += 解:令 u x y = 则: 21u x u dx du x u dx dy +=+= 即2 1u x dx du x = 得到22x dx u du = 故c x u +-=-11 即 21 1x x c y += 另外0=y 也是方程的解。 10. 2 1?? ? ??+=dx dy dx dy x 解:令 p dx dy = 即p p x 2 1+= 而 p dx dy =故两边积分得到 c p p y +-=ln 2 12 因此原方程的解为p p x 21+=,c p p y +-=ln 212 。 12.x y xe dx dy e =?? ? ??+-1 解: y x xe dx dy +=+1
常微分方程与动力系统第二章习题参考答案 1.证明:因为()t Φ是线性齐次系统(LH )的一个基本解矩阵,由定理2.5知()t Φ在区间J 上满足矩阵微分系统()M LH ,即. ()()()t A t t Φ=Φ, . 1 ()()() A t t t -=ΦΦ所以由()A t 确定的线性齐次系统(LH )必唯一。 2.证明:因为()t ?,()t ψ分别是. ()x A t x = 和. ()T x A t x =-的解,所以 11 1 () ()()n k k k n nk k k a d t A t t dt a ????==?? ? ?== ? ? ? ??? ∑∑ , 11211111122222* 121 ()()()n n k k k n n kn k n n n nn k a a a a a a a d t A t t dt a a a a ψψψψψψ==?????? ? ? ? ? ? ?=-ψ=-=- ? ? ? ? ? ? ????? ??? ∑∑ 因而 1111 112 2 1 1 (,)(,)(,),,n n k k k k k k n n kn k k nk k n n k a a d d d dt dt dt a a ψ??ψψ ??ψ?ψ ψ?ψ?ψ?====?? ?? ?????????? ?-?? ? ? ??? ??? ? ? ???=+= ?+?? ? ? ??? ?-?? ? ? ??? ????? ???? ??????? ?? ∑∑∑∑ 11 111 1 1 1()0 n n n n n n n n n n n n m 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 ?ψψ??ψ?ψ?ψ?ψ== === = == == = = -= += -=-=∑∑∑∑∑∑∑ ∑∑∑∑∑所以 (),() ()()1 n t t t t k k k ?ψ?ψ≡≡ ∑=常数。 3.证明:设)t Φ(为系统. ()x A t x = 的一个基本解矩阵,则由定理2.11知 [ ]1 () T t -Φ是系统. ()T x A t x =-的基本解矩阵,由定理 2.4知系统. ()x A t x = 满足初始条件00()x t x =的特解为1 00()))t t t x ?-=Φ(Φ(,[) 0,0,t t ∈+∞由题可 知)t Φ(与[ ]1 () T t -Φ在[)0,+∞上有界,从而由定理2.24知110()0 k k t ?=>
习题2.2 求下列方程的解。 1.dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 2 1e x -(x x cos sin +)+c] =c e x -21 (x x cos sin +)是原方程的解。 2.dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ?-dt 3 (?e t 2 e -? -dt 3c dt +) =e t 3- (5 1e t 5+c) =c e t 3-+5 1e t 2 是原方程的解。 3.dt ds =-s t cos +21t 2sin 解:s=e ?-tdt cos (t 2sin 2 1?e dt dt ?3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为:dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy +1212--y x x =0 解:原方程可化为:dx dy =-1212+-y x x ?=-dx x x e y 1 2(c dx e dx x x +?-221) )21(ln 2+=x e )(1 ln 2?+--c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 234xy x x += 解:dx dy 234xy x x += =23y x +x y 令 x y u = 则 ux y = dx dy =u dx du x + 因此:dx du x u +=2u x 21u dx du = dx du u =2 c x u +=33 1 c x x u +=-33 (*) 将x y u =带入 (*)中 得:3433cx x y =-是原方程的解.
常微分方程习题2.1 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为: x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 22 2 22 2 2 3 22 3 2 )1(1)1)(1(),0(ln 1ln 2 1ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+ ? + =+)故原方程的解为(即两边积分得故分离变量得显然
.0;0;ln ,ln ,ln ln 0 110000 )1()1(4===-==-+=-++=-=+≠===-++x y c y x xy c y x xy c y y x x dy y y dx x x xy x y xdy y ydx x 故原方程的解为即两边积分时,变量分离是方程的解,当或解:由:
1 拉格朗日方程的应用及举例 拉格朗日方程有以下几个特点:(1)拉格朗日方程适用于完整系统,可以获得数目最少的运动微分方程,即可以建立与自由度数目相同的n 个方程,是一个包含n 个二阶常微分方程组,方程组的阶数为2n 。求解这个方程组可得到以广义坐标描述的系统运动方程。(2)拉格朗日方程的形式具有不变性。对于任意坐标具有统一的形式,即不随坐标的选取而变化。特别是解题时有径直的程序可循,应用方便。(3)所有的理想约束的约束反力均不出现在运动微分方程中。系统的约束条件愈多,这个特点带来的便利越突出。(4)拉格朗日方程是以能量的观点建立起来的方程,只含有表征系统运动的动能和表征主动力作用的广义力,避开了力、速度、加速度等矢量的复杂运算。(5)拉格朗日方程不但可以建立相对惯性系的运动,还可以直接建立相对非惯性系的动力学方程,只要写出的动能是绝对运动的动能即可,至于方程所描述的运动是对什么参考系的运动,则取决于所选的广义坐标。 纵观拉格朗日方程,看出分析力学在牛顿力学的基础上,提出严密的分析方法,从描述系统的位形到建立微分方程都带有新的飞跃。我们还应看到,虽然拉格朗日方法在理论上和应用上都有重要的价值,但是,牛顿力学的价值并未降低,特别是它的几何直观性和规格化的方法使人乐于应用,由于计算机的广泛使用,牛顿一欧拉方法又有所发展。我们将会看到,用拉格朗日方程求解,在获得数量最少的运动微分方程时,其求导过程有时过于繁琐,并有较多的耦合项。 应用拉格朗日方程建立动力学方程时,应首先建立以广义坐标q 和广义速度q 表示的动能函数和广义力Q 。为此,首先讨论动能的计算和广义力的计算,在此基础上,再讨论拉格朗日方程的应用。 一、动能的计算 对于系统的动能,可以写出关于广义速度q 的齐次函数的表达式。在实际计算中,应用理论力学的有关知识就可以建立以广义坐标和广义速度所表达的动能函数。 例1-1 已知质量为m ,半径为r 的均质圆盘D ,沿OAB 直角曲杆的AB 段只滚不滑。圆盘的盘面和曲杆均放置在水平面上。已知曲杆以匀角速度ω1绕通过O 点的铅直轴转动,试求圆盘的动能。 解:取广义坐标x 和?,x 为圆盘与曲杆接触点到曲杆A 点的距离,?为曲杆OAB 的转角,? = ω1t 。 应用柯尼希定理求圆盘的动能。为此,先求圆盘质心C 的速度和相对于质心平动坐标
习题2.2 求下列方程的解 1. dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 21 e x -(x x cos sin +)+c] =c e x -2 1 (x x cos sin +)是原方程的解。 2. dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ? -dt 3 (?e t 2 e -?-dt 3c dt +) =e t 3- (5 1 e t 5+c) =c e t 3-+5 1 e t 2 是原方程的解。 3. dt ds =-s t cos + 21t 2sin 解:s=e ? -tdt cos (t 2sin 2 1 ?e dt dt ? 3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为: dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy + 1212 --y x x =0 解:原方程可化为: dx dy =-1212 +-y x x ? =-dx x x e y 2 1 2(c dx e dx x x +? -2 21) ) 2 1(ln 2 + =x e )(1ln 2 ?+- -c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 2 3 4xy x x += 解: dx dy 2 3 4 xy x x += =2 3y x + x y 令 x y u = 则 ux y = dx dy =u dx du x + 因此:dx du x u += 2 u x 2 1u dx du = dx du u =2 c x u +=3 31 c x x u +=-33 (*) 将 x y u =带入 (*)中 得:3 4 3 3cx x y =-是原方程的解.
第二章:动力学系统的微分方程模型 利用计算机进行仿真时,一般情况下要给出系统的数学模型,因此有必要掌握一定的建立数学模型的方法。在动力学系统中,大多数情况下可以使用微分方程来表示系统的动态特性,也可以通过微分方程可以将原来的系统简化为状态方程或者差分方程模型等。在这一章中,重点介绍建系统动态问题的微分方程的基本理论和方法。 在实际工程中,一般把系统分为两种类型,一是连续系统;其数学模型一般是高阶微分方程;另一种是离散系统,它的数学模型是差分方程。 §2.1 动力学系统统基本元件 任何机械系统都是由机械元件组成的,在机械系统中有3种类型的基本机械元件:惯性元件、弹性元件和阻尼元件。 1 惯性元件:惯性元件是指具有质量或转动惯量的元件,惯量可以定义为使加速度(或角加速度)产生单位变化所需要的力(或力矩)。 惯量(质量)= ) 加速度(力(2 /) s m N 惯量(转动惯量)= ) 角加速度(力矩(2/) s rad m N ? 2 弹性元件:它在外力或外力偶作用下可以产生变形的元件,这种元件可以通过外力做功来储存能量。按变形性质可以分为线性元件和非线性元件,通常等效成一弹簧来表示。 对于线性弹簧元件,弹簧中所受到的力与位移成正比,比例常数为弹簧刚度k 。 x k F ?= 这里k 称为弹簧刚度,x ?是弹簧相对于原长的变形量,弹性力的方向总是指向弹簧的原长位移,出了弹簧和受力之间是线性关系以外,还有所谓硬弹簧和软弹簧,它们的受力和弹簧变形之间的关系是一非线性关系。 3 阻尼元件:这种元件是以吸收能量以其它形式消耗能量,而不储存能量,可以形象的表示为一个活塞在一个充满流体介质的油缸中运动。阻尼力通常表示为: α x c R = 阻尼力的方向总是速度方向相反。当1=α,为线性阻尼模型。否则为非线性阻 尼模型。应注意当α等于偶数情况时,要将阻尼力表示为: ||1--=αx x c R 这里的“-”表示与速度方向相反
《常微分方程》第三版答案 习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时c=1 特解为y= e 2 x . 2. y 2 dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。解:y 2dx=-(x+1)dy 2 y dy dy=-1 1+x dx 两边积分: - y 1=-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解x=0,y=1时c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+31 x x + y y 21+dy=31 x x +dx 两边积分:x(1+x 2 )(1+y 2 )=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +-
令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即ln(y 2+x 2)=c-2arctg 2 x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为:tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32+=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c.
习 题 3-1 复摆重P ,对质心的回转半径为C ρ,质心距转动轴的距离为a ,复摆由水平位置无初速地释放,列写复摆的运动微分方程。 解:系统具有一个自由度,选复摆转角?为广义坐标,原点及正方向如如题4-1图所示。 复摆在任意位置下,根据刚体绕定轴转动微分方程 O O M J =? 其中 )(22 a g P J C O += ρ 得到复摆运动微分方程为 ?? ρcos )(22 Pa a g P C =+ 或 0cos )(22 =-+?? ρga a C 3-2均质半圆柱体,质心为C ,与圆心O 1的距离为e ,柱体半径为R ,质量为m ,对质心的回转半径为C ρ,在固定平面上作无滑动滚动,如题3-2图所示,列写该系统的运动微分方程。 解:系统具有一个自由度,选θ为广义坐标。 半圆柱体在任意位置的动能为: 222 1 21ωC C J mv T += 用瞬心法求C v : 2222*2)cos 2()(θθθ Re R e CC v C -+== θω = 2 C C m J ρ= 故 222222 1)cos 2(21θρθθ C m Re R e m T +-+= 系统具有理想约束,重力的元功为 题3-1图 题3-2图
θθδd mge W sin -= 应用动能定理的微分形式 W dT δ= θθθρθθd mge m Re R e m d C sin 21)cos 2(2122222-=?? ????+-+ θθθθθθθθθθ ρd m g e d m R e d m R e d R e m C s i n s i n c o s 2)(2222-=+-++ 等式两边同除dt , θθθθθθθθθθ ρ s i n s i n c o s 2)(2222m g e m R e m R e R e m C -=+-++ 0≠θ ,等式两边同除θ 故微分方程为 0s i n s i n )c o s 2(2222=+++-+θθθθρθ m g e m R e Re R e m C ① 若为小摆动θθ≈sin ,1cos ≈θ,并略去二阶以上微量,上述非线性微分方程可线性化,系统微摆动的微分方程为 0])[(22=++-θθρge r R C 要点及讨论 (1)本题也可以用平面运动微分方程求解。系统的受力图与运动分析图如图(b )所示。列写微分方程 ??? ??--=-=-=④③② θ θθρsin )cos (2Ne e R F m mg N y m F x m C C C 上述方程包含C x ,C y ,θ ,F ,N 五个未知量,必须补充运动学关系才能求解。建立质心坐标与广义坐标θ之间的关系 ?? ?-=-=θθ θcos sin e R y e R x C C , ???=-=θθθθθ sin cos e y e R x C C 所以 ?????+=+-=⑥ ⑤22cos sin sin cos θθθθθθθθθ e e y e e R x C C 运动学方程式⑤⑥与方程②③④联立,消去未知约束力N ,F ,就可以得到与式①相同的系统运动微分方程。 因为在理想约束的情况下,未知约束力在动能定理的表达式中并不出现,所以用动能定理解决已知力求运动的问题更简便、直接。 (2)本题也可用机械能守恒定律求解。 系统的动能 222222 1)c o s 2(21θρθθ C m Re R e m T +-+=
上海交通大学 致远学院 2016年秋季学期 《常微分方程与动力系统》课程教学说明 一.课程基本信息 1.开课学院(系):致远学院 2.课程名称:《常微分方程与动力系统》 (An Introducation to Differential Equations and Dynamical Systems) 3.学时/学分:48学时/ 3学分 4.先修课程:数学分析、高等代数、空间解析几何;或线性代数、高等数学。 5.上课时间:星期五 6-8节(12:55-15:40) 6.上课地点:东下院 101 7.期末考试时间:2017-01-(02-13)考试周 8.任课教师:肖冬梅, xiaodm@https://www.wendangku.net/doc/757925788.html, 9.办公室及电话:数学楼2305,54743151转2305 10.助教:何鸿锦,hehongjin000@https://www.wendangku.net/doc/757925788.html, 11.答疑(office hour):星期三晚18:30 – 20:30,数学楼2305室二.课程主要内容(如何可以,请提供中英文) 除期中考试2学时+习题课2学时外,其余全是课堂教学 第一章基本概念(3学时) 主要内容: 1.1什么是微分方程?什么是常微分方程?常微分方程的分类 1.2什么是常微分方程解?什么是特解?什么是通解? 1.3常微分方程建模:初始值问题和边界值问题 1.4关于常微分方程和解的几何看法:向量场、积分曲线 重点与难点:常微分方程和解的几何观点,方向场和积分曲线的作图 第二章一阶常微分方程的初等解法(6学时) 主要内容: 2.1 变量分离法 2.2 一阶线性常微分方程 2.3 全微分方程(或恰当方程)和积分因子 2.4 替代法和某些可解的常微分方程 重点与难点:全微分方程和积分因子,变换的技巧 第三章基本理论(8学时) 主要内容:
常微分方程 2.1 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为: x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 2 2 2 22 2 2 3 22 3 2 )1(1)1)(1(),0(ln 1ln 21ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+ ? + =+) 故原方程的解为(即两边积分得故分离变量得显然 .0;0;ln ,ln ,ln ln 0 110000 )1()1(4===-==-+=-++=-=+≠===-++x y c y x xy c y x xy c y y x x dy y y dx x x xy x y xdy y ydx x 故原方程的解为即两边积分时,变量分离是方程的解,当或解:由:
10ln 1ln ln 1ln 1,0 ln 0 )ln (ln :931:8. cos ln sin ln 0 7ln sgn arcsin ln sgn arcsin 1 sgn 11,)1(,,,6ln )1ln(2 11 11,11,,,0 )()(:5332 2 22 2 22 2 22 2 c dx dy dx dy x y cy u d u u dx x x y u dx x y dy x y ydx dy y x x c dy y y y y dx dy c x y tgxdx ctgydy ctgxdy tgydx c x x x y c x x u dx x x du x dx du dx du x u dx dy ux y u x y y dx dy x c x arctgu dx x du u u u dx du x u dx du x u dx dy ux y u x y x y x y dx dy dx x y dy x y e e e e e e e e x y u u x y x u u x y x y y x x x +===+=+-===-?-=--+-=-=+-===-=+?=+?=?=--=+===-+=+-=++ =++-++=++===+-==-++-+-- 两边积分解:变量分离:。 代回原变量得:则有:令解:方程可变为:解:变量分离,得 两边积分得:解:变量分离,得::也是方程的解。 另外,代回原来变量,得两边积分得:分离变量得:则原方程化为: 解:令:。两边积分得:变量分离,得:则令解:
常微分方程(第三版) 答案
常微分方程习题答案 2.1 1.?Skip Record If...?,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 ?Skip Record If...??Skip Record If...?并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: ?Skip Record If...?3 ?Skip Record If...? 解:原式可化为: ?Skip Record If...??Skip Record If...??Skip Record If...? ?Skip Record If...? 12.?Skip Record If...? 解?Skip Record If...??Skip Record If...? ?Skip Record If...? 15.?Skip Record If...? ?Skip Record If...?16.?Skip Record If...? 解:?Skip Record If...? ?Skip Record If...?,这是齐次方程,令?Skip Record If...? 17. ?Skip Record If...? 解:原方程化为?Skip Record If...? 令?Skip Record If...? 方程组?Skip Record If...??Skip Record If...? 则有?Skip Record If...? 令?Skip Record If...? 当?Skip Record If...?当?Skip Record If...? 另外 ?Skip Record If...? ?Skip Record If...?
习题3.1 1 求方程dx dy =x+y 2通过点(0,0)的第三次近似解; 解: 取0)(0=x ? 20020012 1)()(x xdx dx y x y x x x ==++=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x +=+=++=???? dx x x x y x x ])20 121([)(252003+++=?? = 118524400 1160120121x x x x +++ 2 求方程dx dy =x-y 2通过点(1,0)的第三次近似解; 解: 令0)(0=x ? 则 20020012 1)()(x xdx dx y x y x x x ==-+=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x -=-=-+=???? dx x x x y x x ])20 121([)(252003--+=?? =118524400 1160120121x x x x -+- 3 题 求初值问题: ?????=--=0 )1(22y y x dx dy R :1+x ≤1,y ≤1 的解的存在区间,并求解第二次近似解,给出在解的存在空间的误差估计; 解: 因为 M=max{22y x -}=4 则h=min(a,M b )=4 1 则解的存在区间为0x x -=)1(--x =1+x ≤4 1 令 )(0X ψ=0 ; )(1x ψ=y 0+?-x x x 0)0(2dx=31x 3+31;
)(2x ψ =y 0+])3131([2132?-+-x x x dx=31x 3-9x -184x -637x +4211 又 y y x f ??),(2≤=L 则:误差估计为:)()(2x x ψ-ψ≤32 2 )12(*h L M +=2411 4 题 讨论方程:31 23y dx dy =在怎样的区域中满足解的存在唯一性定理的条件, 并求通过点(0,0)的一切解; 解:因为y y x f ??),(=3221-y 在y 0≠上存在且连续; 而312 3y 在y 0 σ≥上连续 由 3123y dx dy =有:y =(x+c )23 又 因为y(0)=0 所以:y =x 2 3 另外 y=0也是方程的解; 故 方程的解为:y =?????≥00023 x x x 或 y=0; 6题 证明格朗瓦耳不等式: 设K 为非负整数,f(t)和g(t)为区间βα≤≤t 上的连续非负函数,
习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+3 1 x x + y y 21+dy=3 1 x x +dx 两边积分:x(1+x 2 )(1+y 2 )=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0
解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c. 9.x(lnx-lny)dy-ydx=0 解:原方程为: dx dy =x y ln x y
拉格朗日方程的应用及举例 拉格朗日方程有以下几个特点:(1)拉格朗日方程适用于完整系统,可以获得数目最少的运动微分方程,即可以建立与自由度数目相同的n个方程,是一个包含n个二阶常微分方程组,方程组的阶数为2n。求解这个方程组可得到以广义坐标描述的系统运动方程。(2)拉格朗日方程的形式具有不变性。对于任意坐标具有统一的形式,即不随坐标的选取而变化。特别是解题时有径直的程序可循,应用方便。(3)所有的理想约束的约束反力均不出现在运动微分方程中。系统的约束条件愈多,这个特点带来的便利越突出。(4)拉格朗日方程是以能量的观点建立起来的方程,只含有表征系统运动的动能和表征主动力作用的广义力,避开了力、速度、加速度等矢量的复杂运算。(5)拉格朗日方程不但可以建立相对惯性系的运动,还可以直接建立相对非惯性系的动力学方程,只要写出的动能是绝对运动的动能即可,至于方程所描述的运动是对什么参考系的运动,则取决于所选的广义坐标。 纵观拉格朗日方程,看出分析力学在牛顿力学的基础上,提出严密的分析方法,从描述系统的位形到建立微分方程都带有新的飞跃。我们还应看到,虽然拉格朗日方法在理论上和应用上都有重要的价值,但是,牛顿力学的价值并未降低,特别是它的几何直观性和规格化的方法使人乐于应用,由于计算机的广泛使用,牛顿一欧拉方法又有所发展。我们将会看到,用拉格朗日方程求解,在获得数量最少的运动微分方程时,其求导过程有时过于繁琐,并有较多的耦合项。 应用拉格朗日方程建立动力学方程时,应首先建立以广义坐标q和广义速度q 表示的动能函数和广义力Q。为此,首先讨论动能的计算和广义力的计算,在此基础上,再讨论拉格朗日方程的应用。 一、动能的计算 对于系统的动能,可以写出关于广义速度q 的齐次函数的表达式。在实际计算中,应用理论力学的有关知识就可以建立以广义坐标和广义速度所表达的动能函数。 例1-1已知质量为m,半径为r的均质圆盘D, 沿OAB直角曲杆的AB段只滚不滑。圆盘的盘面和曲 杆均放置在水平面上。已知曲杆以匀角速度 1绕通过 O点的铅直轴转动,试求圆盘的动能。 解:取广义坐标x和 ,x为圆盘与曲杆接触点到 曲杆A点的距离, 为曲杆OAB的转角, = 1t。 应用柯尼希定理求圆盘的动能。为此,先求圆盘质心C 的速度和相对于质心平动坐标标准
习题2.2 求下列方程的解 1.dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 2 1e x -(x x cos sin +)+c] =c e x -21 (x x cos sin +)是原方程的解。 2.dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ?-dt 3 (?e t 2 e -? -dt 3c dt +) =e t 3- (5 1e t 5+c) =c e t 3-+5 1e t 2 是原方程的解。 3.dt ds =-s t cos +21t 2sin 解:s=e ?-tdt cos (t 2sin 2 1?e dt dt ?3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为:dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy +1212--y x x =0 解:原方程可化为:dx dy =-1212+-y x x ?=-dx x x e y 21 2(c dx e dx x x +?-221) )21(ln 2+=x e )(1 ln 2?+--c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 234xy x x += 解:dx dy 234xy x x += =23y x +x y 令 x y u = 则 ux y = d x d y =u dx du x + 因此:dx du x u +=2u x 21u dx du = dx du u =2 c x u +=33 1 c x x u +=-33 (*) 将x y u =带入 (*)中 得:3433cx x y =-是原方程的解.
动力系统的概念 这一章是对于事实的调查,而且来源于应用于全书的动力系统理论。我们的主要目的是为后面的章节确定固定使用的常用符号和专业术语,并且回想一些常常在课本的前言中不被讨论的理论的一些方面。为了更容易的阅读,我们保持讨论时采用非专业术语,并尽可能地避免技术上的符号和观点。然而许多遗漏的细节可以从研究生使用的动力系统的课本的前言中找到,一些更加先进的课题仅仅在研究性的文章中涉及到。在某些情况下,我们将提供一些在更深的章节中关于这个主题的参考。另外,我们鼓励读者使用附录A 和B 作为基于不同的几何和函数分析的参考。 流量,映射,动力系统 对于任意的集合P ,一个变换群:P P t F →中的任意的一个参数t 属于实数,如果 ()x x F =0对于所有的x 属于集合P ,并且s t s t F F F ο=+对于任意的 ,t s , 属于实数都成立, 则被称为一个流。这两个属性表明t F 和它的逆t F -是不可以转化的。这一组合t (,)p F 叫做基于空间P 的一个连续的动力系统。换句话说,一个连续的动力系统包括一个可能状态集合和唯一决定将来状态)(x F t 的当前的状态函数x 的变化规则。通过x 这一点的变化轨迹是集 合)()(x F U x t R t ∈=γ。一个固定点的流是一个点x 且x x F t =)(对于任意的R t ∈都成立。 这个流的一个周期的轨迹就是通过这一点x 对于那些存在的正数T,并且满足x x F T =)(的这 样的轨迹。 如果用以上所说的映射族t F 定义只需0≥t ,且对于所有的t ,s 满足()x x F =0和s t s t F F F ο=+,则t F 叫做半流形。注:半流形通常是不可逆的,动力系统的一个典型的特征是在无穷大的空间中是确定的。 当有单独向映射P P f →:且存在()f P ,时,离散动力系统是确定的。这样的系统还有一些性质即通过f 的迭代次数可以得出唯一的当前状态决定所有的将来状态()(),...,2x f x f 。这时x 的取值范围是确定的在集合()()Y Z n n x f x ∈=γ中,其中
常微分方程第三版答 案
习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+3 1 x x + y y 21+dy=3 1 x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。
5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2 x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c.
Xxxxxx U N I V E R S I T Y 《微分方程定性理论》实践报告 所属学院:理学院 专业班级:应用数学 姓名: 学号:xxxxxxxxxxx 实践课题:动力系统综述 实践成绩: 任课教师:
动力系统综述 随着数学知识的不断扩充及科学技术的不断发展,动力系统被广泛应用于工程、力学、生态等各大领域,推动着社会的发展。动力系统是随时间而演变的系统。 随着数学知识的不断扩充及科学技术的不断发展,动力系统被广泛应用于工程、力学、生态等各大领域,推动着社会的发展。动力系统是随时间而演变的系统。对于含参数的系统,当参数变化并经过某些临界值时,系统的定性性态,如平衡点或周期运动的数目和稳定性等会发生突然变化,这种变化称为分叉[2]。 分叉理论主要研究当参数在分叉值附近变化时,系统轨线的拓扑结构或定性性态将如何变化。近几十年来,动力系统的分叉理论被系统而深入的研究,并得到了迅猛的发展,且广泛应用于物理、化学、生物、工程等研究领域中,分叉问题的研究己成为非线性动力系统研究的重点和难点之一。 1动力系统简介 动力系统的研究起源于牛顿的经典力学理论.假设空间R n 的一个质点M 在时刻t 的坐标为),,,(21n x x x x =并且己知质点M 此时的运动速度为))(,),(),(()(21x v x v x v x v n =,并且只与坐标x 有关.那么质点M 的运动方程为: )(x v dt dx = (1) 这个方程是一个自治的微分方程.更进一步如果方程(1)满足微分方程解的存在和唯一性定理的条件,那么对任何的初值条件00)(x t x =,则方程存在唯一解),,()(00x t t t =?。 我们称x 取值的空间n ?为相空间,而称((t , x )的取值空间“n ???”为增广相空间.按照微分方程的几何意义,方程(1)定义了增广相空间中的一个向量场.解的几何意义为增广相空间中经过点),(00x t 的唯一的积分曲线[1]. 2 动力系统在力学中的应用 稳定性是系统的一个重要特性。对系统运动稳定性分析是系统与控制论的一个重要组成部分,一个实际的系统必须是稳定的,不稳定的系统是不能付诸于工程实施的。 设系统的向量状态方程为: 0,)(),,(00≥==t x t x t x f x (2.1) 式中:x 为n 维状态向量;),(??f 为n 维向量函数。