二阶常系数线性微分方程的应用举例
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第七章常微分方程7.13 二阶常系数线性微分
方程的应用举例
数学与统计学院
赵小艳
解 受力分析 例1 (弹簧的机械振动)
如图,弹簧下挂一物体.设在垂直方向有一随时间变化的外力
作用在物体上,物体将受外力驱使而上下振动,求物体的振动规律.
pt H t f sin )(1=x x
o )(1t f ;sin )()1(1pt H t f =外力;)2(kx f -=弹性力v f μ-=0)3(介质阻力,ma F =由t x kx pt H t x m d d d d μ--=sin 22可得.t x d d μ-= 设振动开始时刻为0,t 时刻物体离开平衡位置
的位移为x (t ).
.0,000====t t x x d 还应满足初始条件:
.0,000====t t x x d 还应满足初始条件:2m t x kx pt H t x m d d d d μ--=sin 22
可得m m pt h x t x t x sin 2222=++ωδd d d d 强迫振动的微分方程
2m m m pt h x t x t x sin 2222=++ωδd d d d 强迫振动的微分方程 对应齐次方程: 02222=++x t x t
x ωδd d d d 自由振动的微分方程 其特征方程: 0222=++ωδλλ.
,222221ωδδλωδδλ---=-+-=.
0)1(22>-ωδ.)(2)(12222t
t e C e C x ωδδωδδ-+----+=齐次方程的通解为
.0)(→∞→t x t 时,当此时物体运动按指数函数规律衰减.
t x O
对应齐次方程: 02222=++x t x t
x ωδd d d d 齐次方程的特征方程: 0
222=++ωδλλ.
0)2(22=-ωδ.
)(21t e t C C x δ-+=齐次方程的通解为物体运动仍按指数规律衰减.
.0)3(22<-ωδ.
222,1δωδλ-±-=i .)sin cos (222221t C t C e x t δωδωδ-+-=-齐次方程的通解为
t x O
02222=++x t x t x ωδd d d d .
0)2(22=-ωδ.
)(21t e t C C x δ-+=齐次方程的通解为物体运动仍按指数规律衰减. .0)3(22<-ωδ.
222,1δωδλ-±-=i .)sin cos (2222
21t C t C e x t δωδωδ-+-=-齐次方程的通解为.)sin(22ϕδωδ+-=-t Ae t .arctan ,2
12221是两个任意常数其中C C C C A =+=ϕ一方面,物体运动的振幅随时间的增大而减少,
t x
O t x O
pt h x t x t
x sin 2222=++ωδd d d d 非齐次方程的特解
,0≠μ设2221ωδδλ-±-=,,不是特征值ip ,
sin cos 21pt A pt A x +=*可令特解为⎩⎨⎧=--=+-h pA A p pA A p 122221222)(02)(δωδω得⎪⎪
⎩⎪⎪⎨⎧+--=+--=⇒222222222222214)()(4)(2p p p h A p p ph A δωωδωδpt p p ph pt p
p p h x cos 4)(2sin 4)()
(222222222222δωδδωω+--+--=*)sin(ψ-=pt B .2arctan ,⎪⎫ ⎛==p h B δψ其中
pt h x t x t
x sin 2222=++ωδd d d d pt ph pt p h x cos 2sin )
(22δω--=*)sin(ψ-=pt B 非齐次方程的通解为.),sin()sin(2
2ωδψϕδωδ<-++-=-其中pt B t Ae x t 暂态项
稳态项 ).sin(,1ψ-=≈>>*pt B x x t 时当取决于外界强迫力的作用
,很小时当阻尼μ222224)(p
p h B δω+-=振幅.22p h -≈ω.)sin(2
2ϕδωδ+-=-t Ae x t
pt h x t x t
x sin 2222=++ωδd d d d 2
221ωδδλ-±-=,非齐次方程的通解为.),sin()sin(2
2ωδψϕδωδ<-++-=-其中pt B t Ae x t 暂态项 稳态项
).sin(,1ψ-=≈>>*pt B x x t 时当取决于外界强迫力的作用
,很小时当阻尼μ222224)(p
p h B δω+-=振幅.
22p h -≈ω,相差不大时与弹簧的固有频率当外界强迫力的频率m
k p =ω