§3.2 非齐次方程—纯强迫振动
一、定解问题
我们来考虑有界弦(或杆)的纯强迫振动
?
?
?
?
??
?
?
?
=
=
=
=
+=
=
=
=
=
0|
0|
0|
0|
),(
0
0
0
2
tt
t
lx
x
xxtt
u
u
u
u
txfuau <1>
<2>
<3>
二、求解思路
思路① :对于
??
??
?
=
=
+=
=
=
0|
0|
),(
0
0
2
tt
t
xxtt
u
u
txfuau
我们可以先由冲量定理求解:
??
??
?
=
=
=
=
=
),(|
0|
2
tt
t
xfv
v
vav
tt
t
xxtt
可由达朗贝尔公式求出);,( ttxv
∫= t dtxvtxu 0 );,(),( tt而
所以我们可以想到,对于:
?
?
?
?
??
?
?
?
=
=
=
=
+=
=
=
=
=
0|
0|
0|
0|
),(
0
0
0
2
tt
t
lx
x
xxtt
u
u
u
u
txfuau
也可以先用冲量原理求解
?
?
?
?
??
?
?
?
=
=
=
=
=
=
=
=
=
0|
0|
0|
0| 0
2
t
t
tt
t
lx
x
xxtt
v
v
v
v
vav
根据冲量原理,先求解:
式求出可由有界弦自由振动公其中 );,( ttxv
∫= t dtxvtxu 0 );,(),( tt这时
思路②:
我们考虑二阶非齐次的常微分方程的求解:
(*))()()()( xfyxQyxpxy =+′+′′对于
我们可采取常数变易法求解,即:
先求解对应的齐次方程:
0)()()( =+′+′′ yxQyxpxy
若它有通解:
)()()( 2211 xyCxyCxyg +=
则可令对应的齐次方程(*)有通解:
)()()()()( 2211 xyxCxyxCxyg +=
s
g
yxCxC
xy
从而求出、
求出并附加一个条件,即可代入将
)()(
(*))(
21
则对于定解问题<1>~<3>,我们也可先
考虑对应的齐次问题:
??
??
?
=
=
=
=
=
0|
0|
0
0
2
tt
t
xxtt
u
u
uau
再采用类似于常数变量法类似的方法求对
应的非齐次问题的特解。
三、求解:
1、对应的齐次条件
??
??
?
=
=
=
=
=
0|
0|
0
0
2
tt
t
xxtt
u
u
uau
[这一问题用分离变量法是易于求
得 本征值问题])(xX
??
???
=
=
=?′′
0)(
0)0(
0
lX
X
XX m
解得:
l
xnCxX
nln
nn
p
pm
sin)(
,...2,1,)( 2
=
=?=
)()(),( tTxXtxu =令
仿照常数变易法,令:
∑∞
=
=
1
sin)()(
n
n l
xntTtxu p, <4>
将之代入式<1>得:
),(sin)]()()([
1
2 txf
l
xntT
l
antT
n
nn =+
″∑∞
=
pp
由傅里叶级数的系数公式有:
2、求解 的方程)(tTn
)()()()( 2 tftTlantT nnn =+″ p
其中: apaa dlntfltf ln ∫= 0 sin),(2)(
再将式<4>代入初始条件<3>,得:
?
?
?
???
?
=′
=
∑
∑
∞
=
∞
=
1
1
0sin)0(
0sin)0(
n
n
n
n
l
xnT
l
xnT
p
p
<5>
<6>
比较两边的系数有:
??
???
=′
=
0)0(
0)0(
n
n
T
T
有常微分方程的常数变易法可求得关于方程
的解为:和的定解问题 ><>< 75)(tTn
ttptp dtl anfanltT t nn ∫ ?=
0
)(sin)()(
<7>
3、有界弦(杆)的纯强迫振动的解:
<8>
将<8>代入<4>,得定解问题<1>~<3>的解为:
xlndtl anfanltxu
n
t
n
pttpt
p sin])(sin)([),( 1 0∑ ∫
∞
=
?=
式给出。由其中 >< 6)(tfn
4、本征函数法
以上求解非齐次方程的方法,显然也适用于
求解带有其他齐次边界条件的各类方程。
其中主要步骤为:
(1) 用分离变量法求得对应的齐次问题(即对应的齐次方
程连同齐次边界条件)的本征函数。
(2)将未知函数 [或 等]按上面求得的本征函
数展开,其展开系数为另一变量的函数,代入非齐次方程
和初始条件(或另一变量的边界条件),得到关于时间因
子的常微分方程的初值条件(或另一单元函数的常微分方
程的边值问题),用常数变易法或拉氏变换法可求得其
解。
(3)将所求得的解代入未知函数的展开式中,即得到原定解
问题的解。这种分离变量的方法按其特点又叫本征函数
(或固有函数)法。
)( txu ,)( yxu ,
五、例题
求解定解问题:
??
??
?
=
==
=?
=
==
0|
0|,0|
sin
0
0
2
t
lxxxx
xxt
u
uu
tAuau w <1>
<2>
<3>
解:对应的齐次方程的本征值问题为
??
?
=′=′
=?′′
0)(,0)0(
0
lXX
XX m
求解得本征函数为:
,...2,1,0,cos)( == nl xnCxX nn p
令: ∑
∞
=
=
0
cos)(
n
n l
xntTu p
代入方程<1>和初始条件<3>得:
?
?
?
???
?
=′
=+′
∑
∑
∞
=
∞
=
0
0
2
0cos)0(
sincos)]()()([
n
n
n
nn
l
xnT
tAl xntTl xntT
p
wpp
<4>
比较等式两边的傅里叶系数有:
)0(
0)0(
sin)(
0
0 =
??
???
=
=′ n
T
AtT w
)321(
0)0(
0)()()( 2 …=
??
???
=
=+′ ,,n
T
tTl antT
n
nn
p
解之得:
,...)3,2,1(0)(;)cos1()(0 ==?= ntTtAtT nww
代入<4>式得: )cos1(),( tAtxu ww ?=
本节注意:
(1)以上方法也适于求解带有其他齐次边界条件的非齐次
方程的定解条件,其主要精神是:
<i> 先考虑对应的齐次问题,用分离变量法求得其齐次问
题的固有函数;
<ii> 将未知函数按本征函数展开,其展开系数为另一变
量的系数,代入原非齐次方程和初始条件(或另一变量的
边界条件),得另一变量的常微分方程定解问题;
<iii> 求常微分方程定解问题的解代入展开式得原定解问
题得解。这种分离变量的方法按其特点又叫本征函数。
(2) 对于有界弦的一般强迫振动:
??
??
?
==
==
=?
==
==
)(|,)(|
0|,0|
),(
00
0
2
xuxu
uu
txfuau
ttt
lxx
xxtt
yj
我们可以由迭加原理将之化为两个
以前讨论过的问题来求解,即:
?
?
?
?
??
?
?
?
=
=
=
=
=
=
=
=
=
)(|
)(|
0|
0|
0
0
0
2
xu
xu
u
u
uau
tt
I
t
I
lx
I
x
I
xxIttI
y
j
?
?
?
?
??
?
?
?
=
=
=
=
+=
=
=
=
=
0|
0|
0|
0|
),(
0
0
0
2
tt
II
t
II
lx
II
x
II
xxIIttII
u
u
u
u
txfuau
+
则:
III uuu +=
= 有界弦自由振动解
+
有界弦纯强迫振动解