第三章 分离变量法
用分离变量法求解各种有界问题。
目的:
2、分离变量法的解题步骤;
3、本征值问题:
特殊函数微分方程、 分离变量?? →?
??
?
=+?
=?
0
04
uu
u
l
1、有界弦自由振动的解;
中心内容:
§3.1 有界弦的自由振动
?
?
?
?
?
??
?
?
?
?
??
???
=
=
??
?
=
=
=
=
=
=
=
)(
)(
0|
0|
0
0
2
xu
xu
u
u
uau
lxt
t
lx
x
xxtt
y
j
<1>
<2>
<3>
考虑长为两端固定的弦的自由振动:l
解题思路:
1、当弦一端固定时,其自由振动可看为反射波;
2、当弦两端固定时,其自由振动会形成驻波;
)(2cos1 lgp xtAu ?=正波:
反波: )(2cos2 lgp
xtAu +=
驻波:
l
ppg xtAuuu 2cos2cos2
21 =+= )()( tTxX=
我们设定解条件<1>---<3>的特解为:
)()(),( tTxXtxu =
一、分离变量
则由式<1>,得:
m=′′=′′ XXTaT2
即:
?
?
?
=?′′
=?′′
0
0
2TaT
XX
m
m
则由式<2>,得:
??
?
=
=
0)()(
0)()0(
tTlX
tTX
即:
??
?
=
=
0)(
0)0(
lX
X
同理由式<3>可得:
??
?
=′′
=
)()0()(
)()0()(
xTxX
xTxX
y
j
不成立此二式,任意函数是和由于 )()( xx yj
二、本征值问题
1、先考虑定解问题
??
??
?
??
?
=
=
=?′′
0)(
0)0(
0
lX
X
XX m <4>
<5>
(1) 定义:
件的是常数,只能在边界条式中,m>< 4
值,称些特定的限制下取特定的值,这 m
称为本征值问题。
和相应本征值的问题征函数。求本征值
的非零解称为本方程相应于不同的
下的本征值;在边界条件为方程
m
m ><
><><
4
54
2、求解
得:
???
???
?
=
==
xlnxX
nlnk
n
p
p
sin)(
...321
本征函数:
,,,,本征值:
三、关于T(t)方程的通解
将本征值: 2
22
l
n pm ?=
代入方程: 02 =?′′ TaT m
得: 0)()( 2
222
=?′′ tTlnatT p
此方程通解为:
tl anBtl anAtT nnn pp sincos)( ′+′=
由特解:
)()(),( tTxXtxu =
得式<1>满足式<2>的特解为:
)()(),( tTxXtxu nnn =
l
xnt
l
anBt
l
anA
nn
ppp sin)sincos( += <6>
nnnnnn CBBCAA
′=′= ,式中
四、有界弦的自由振动的解
由迭加原理,将式<6>的解迭加,得:
∑∞
=
+=
1
sin)sincos(),(
n
nn l
xnt
l
anBt
l
anAtxu ppp
将上式代如式<3>,得:
?
?
?
???
?
=
=
∑
∑
∞
=
∞
=
1
1
sinsin)(
sin)(
n
n
n
n
l
xn
l
anBx
xlnAx
ppy
pj
则:
???
???
?
=
=
∫
∫
l
n
l
n
dl ananB
dl anlA
0
0
)(2
)(2
apayp
apaj
此即是定解条件<1>---<3>的解
五、解的物理意义
<6>式可从新改写为:
l
ntNtxu
nnnn
ptdw sin)cos(),( ?=
其中
,222 nnn BAN +=
,lnn ptw =
n
n
n A
Btg 1?=d
代表一个驻波可见 ),( txun
布代表弦上各点的振幅分l xnNn psin
是初位相是位相因子, nnnt ddw )cos( ?
是弦振动的固有频率nw
这些点振幅为零,保持不动,称为节点
0sinsin)...10( ==== pp mlxnnmnmlx mm ,、、当
02 1sinsin
)...21(2 12
=?=
=?=
pp klxn
nklnkx
k
k 时,、当
这些点振幅达到最大值,称为腹点
∑∞
=
=
1
),(
n
nutxu弦的振动 ,表示一系列振幅不
同,频率不同,位相不同的驻波的迭加。
六、小结
1、分离变量法
我们称 的驻波为基波; 的驻波为二次驻波1=n 2=n
分离变量法的主要精神是:把未知函数按自
变量(包括多个自变量)的单元函数分开,
[如:令 ]从而将偏微分
方程的问题化为解常微分方程的问题。
)()()( tTxXtxu =,
2、解题步骤
(1) 对齐次方程和齐次边界条件分离变量;
(2) 解关于空间因子的常微方程的本征值
问题;
(3)求其它常微分方程的解,与本征函数
相乘,得到特解。
(4) 迭加,由初始条件或非齐次边界条
件确定迭加系数,而最后得到所求定解
问题的解。
七、例题分析
解:其定解问题可表示为:
?
?
?
?
??
?
?
?
=
=
??
?
=
=
<<<<=+
=
=
=
=
)(|
0|
0|
0|
0,0,0
0
0
xfu
u
u
u
byaxuu
y
by
axx
xx
yyxx <7>
<8>
<9>
<10>
考察在矩形薄板内稳定状态的温度分布,
板的两对边绝热,而其余的两边一边温度保持
为零,另一边的温度由 规定。
)( xf
)()(),( yYxXyxu =令
分别代入<7>、<8>、<9>式,得:
0)(
0)(
0)0(
0
0
=
??
?
=′
=′
=+′′
=?′′
bY
aX
X
YY
XX
m
m
解本征值问题:
??
?
=′=′
=?′′
0)(,0)0(
0
aXX
XX m
<11>
得:
xanAxX
nan
nn
p
pm
cos)(
...210,2
22
=
=?= 、、、
的方程得:代入将 )(2
22
yYan pm ?=
02
22
=?′′ YanY p
其通解为:
??
???
≠+=+
=+
= )0(,)(
)0(,
)(
00
nFyanEshyanDshyanCch
nDyC
yY ppp
其中 DCthnaFCDE 122 , ?=+= p
??
???
≠=+
==+
)0(0)(
)0(000
nFbanEsh
nDbC
p
故有:
)0(,00 ≠?=?= EbFbDC
由边界条件 有:0)( =by
因此:
???
???
?
≠?
=?
=
)0()(
)0(
)(
0
0
nbyanshE
nb ybD
yYn p
由式<11>得:
???
???
?
≠?
=?
=
=
)0()(cos
)0(2
)()(),(
0
nbyanxshana
nb yba
yYxXyxu
n
nnn
pp
其中 nnn EAaDAa == ,2 000
于是,我们有:
∑
∑
∞
=
∞
=
?+?=
=
1
0
0
)(cos2
),(),(
n
n
n
n
byanxshanab yba
yxuyxu
pp
将非齐次边界条件式<10>代入上边的解,得:
∑∞
=
?+=
1
0 )(cos
2)( n n a
bnxsh
a
naaxf pp
故
...3,2,1,cos)(2 =?= ∫ ndxaxnxf
a
bnasha
a
o
n
p
p
dxxfaa
a
o
∫= )(20
故此问题的解为:
∑∞
=
?
+?=
1
*0 cos
)(
2),( n n xa
n
a
bnsh
byansh
ab ybayxu pp
p
其中: dxaxnxfaa
a
o
n
p∫= cos)(2*