习题3.1
1.考察长为l 的均匀细杆的导热问题,若
(1)杆的两端温度保零度; (2)杆的两端均绝热;
(3)杆的一端为恒温零度,另一端绝热,而初始温度分布均为)(x ?; 试用分离变量法求解在这三种情况下的杆的导热问题的解。 解:(1)该问题的数学模型为
??
???=>==><<=)()0,(0,0),(),0(0
,0,2x x u t t l u t u t l x u a u xx t ? 其D a =2
Step1:分离变量:令)()(),(t T x X t x u =,代入齐次方程及齐次边界条件有:
)()()()0()()()()(2
==''='t T l X t T X t T x x a t T x X
由于0)(,0)(≠≠t T x X 所以有:
)()0()
()()
()(2
==='=''?
l X X t T a t T x X x X λ
整理得
)()(0
)()0(0
)()(2
=-'===-''t T a t T l X X x X x X λλ
Step2:求解特征值问题
?
?
?===-''0)()0(0
)()(l X X x X x X λ 讨论:若0=λ,则0)(=''x X
此时Dx C x X +=)(
将0)()0(==l X X 代入得0==D C ,于是0)(=x X
∴ λ=0不合适,舍去。
若0>λ时方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ λ±=2,1r
∴ x
x
De
Ce
x X λλ-+=)(
将0)()0(==l X X 代入得
?????=+=+-0
0l
l De Ce
D C λλ 得C=-D=0
∴ 0)(=x X ∴ 0>λ也不合适,舍去。 若0<λ时方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ i r λ-±=2,1
∴ x D x C x X λλ-+-=s i n c o s )( 将0)0(=X 代入有C=0 将0)(=l X 代入有0sin
=-l D λ
∵ 0≠D ∴ 0s i n =-l λ ∴
πλn l =- ∴ ,2,1,)(
2
=-=n l
n n πλ
此时 ,2,1 sin )(==n x l
n D x X n n π
Step3:将2
)(
l
n n πλ-=代入关于)(t T 的常微分方程有
0)()(
)(2
=+'t T l a n t T n n π
∴ ,2,1 )(2
)(==-n e
C t T t
l a n n n π
∴ t
l a n n n n n n e
C l
x n D t T x X t x u 2
)(s i n
)()(),(ππ-==
x l
n e
a t
l
a n n ππsin
2
)(-= ,2,1=n
其中n a 待求。 Step4:叠加
∑
∑
∞
=∞
=-=
=
1
1
)(
sin
),(),(2
n n t
l
a n n n x l
n e
a t x u t x u ππ
由于)()0,(x x u ?=
所以:∑
∞
==
1
sin
)(n n x l
n a x π?
∴ ,2,1,s i n
)(20
==
?
n x d x l
n x l
a l
n π?
综上知该热传导问题的解为:
∑
∞
=-=
1
)(
sin
),(2
n t
l
a n n x l
n e
a t x u ππ
其中 ?
==
l
n n x d x l
n x l
a 0
,2,1,s i n
)(2 π?
(2)该问题的数字模型为
??
?
??=>==><<=)()0,(0,0),(),0(0
,0,2x x u t t l u t u t l x u a u x x xx t ? 其D a =2
Step1:分离变量:令)()(),(t T x X t x u =,并代入齐次方程及齐次边界条件中有:
)()()()0()()()()(2
='='''='t T l X t T X t T x x a t T x X
由于0)()(≠t ,T x X 故上面方程可化为 所以有:
)()0()
()()
()(2
='='='=''?
l X X t T a t T x X x X λ
整理得
)()0(0)()(0
)()(2='='=-'=-''l X X t T a t T x X x X λλ
Step2:求解下面的特征值问题
??
?='='=-''0
)()0(0
)()(l X X x X x X λ 讨论:若0=λ,则方程变为0)(=''x X
这样 x D C x X 00)(+=
∵ 0)(D x X =' 而0)()0(='='l X X ∴ 00=D 此时有 0)(C x X = (取≠0C 0即可)
若0>λ时 则方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ λ±=2,1r
∴ x
x
De
Ce
x X λλ-+=)(
x
x
e
D e C x X λλλλ-
-=')(
∵ 0)0(='X ∴ 0=-D C 0)(='l X ∴ 0=--l
l
e
D e C λλλλ
∴ C=-D=0 此时0)(≡x X ∴ 0>λ不合适
若0<λ时,方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ i r λ-±=2,1
∴ x D x C x X λλ-+-=s i n c o s )(
x D
x C
x X λλλλ--+---='cos
sin
)(
∵ 0)0(=X ∴ 0=D
∵ 0)(='l X ∴ 0s i n =-l C λ ∵ 0≠C ∴ 0s i n =-l λ ∴
,2,1,)(
2
=-==-n l
n n l n πλπλ
∴ ,2,1,c o s
)(==n l
x n C x X n n π
综上所述。
该问题的特征值为 2
)(
l n n πλ-=
特征函数为 2,1,0,c o s )(==n x l
n C x X n n π
Step3:将2
)(
l
n n πλ-=代入关于)(t T 的微分方程求)(t T 有:
0)()(
)(2
=+'t T l a n t T n n π
∴ ,2,1,0 )(2
)(==-n e
D t T t
l a n n n π
∴ t
l
a n n n n n n e
xD l
n C t T x X t x u 2
)(
cos
)()(),(ππ-==
,2,1,0,cos
2
)(==-n xe
l
n a t
l
a n n ππ
Step4:叠加,原方程的解为:
∑
∑
∞
=∞
=-=
=
1
1
)(
cos
),(),(2
n n t
l
a n n n x l
n e
a t x u t x u ππ
∵ )()0,(x x u ?=
∴ ∑∑∞
=∞
=+
==1
1
0c o s 22c o s )(n n n
n x l n a a x l n a x π
π
? 其中 d x x l
a l
)(1
0?=
? ,2,1,cos
)(20
==
?
n xdx l
n x l
a l
n π?
(3)该问题的数字模型为:
??
?
??=>==><<=)()0,(0,0),(,0),0(0
,0,2x x u t t l u t u t l x u a u x xx t ? 其中D a =2
Step1:分离变量:令)()(),(t T x X t x u =,代入齐次方程及齐次边界条件中有:
)()()()0()()()()(2
='=''='t T l X t T X t T x x a t T x X
由于0)(,0)(≠≠t T x X 整理上面方程有 所以有:
)()0()
()()
()(2
='=='=''?
l X X t T a t T x X x X λ
整理得
)()(0
)()0(0
)()(2
=-'='==-''t T a t T l X X x X x X λλ
Step2:求解下面的特征值问题
??
?='==-''0
)()0(0
)()(l X X x X x X λ
讨论:若0=λ,则0)(=''x X ∴ Dx C x X +=)(
这样D x X =')(
∵ 00)0(=?=C X 00)(=?='D l X ∴ 0)(=x X
∴ 0=λ不合适,舍去。
若0>λ时 则方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ λ±=2,1r
∴ x
x
De
Ce
x X λλ-+=)(
x
x
e
D e C x X λλλλ-
-=')(
∵ 0)0(=X ∴ 0=+D C 0)(='l X ∴ 0=--l
l
e
D e C λλλλ
∴ C=-D=0
∴ 这说明0>λ也不合适
若0<λ时,方程0)()(=-''x X x X λ的特征方程为
02
=-λr ∴ i r λ-±=2,1
∴ x D x C x X λλ-+-=s i n c o s )(
x D
x C
x X λλλλ--+---='cos
sin
)(
∵ 0)0(=X ∴ 0=C
∵ 0)(='l X ∴ 0c o s =-l C λ ∵ 0≠C ∴ 0c o s =-l λ
∴
,2,1,0,)21
(
2=+
-=n l
n n πλ
∴ x l
n D x X n n π21s i n
)(+
=
Step3:将2)21
(
πλl
n n +
-=代入关于
)(t T 的微分方程有:
0)()21
(
)(2=+
+'t T a l n t T n
n π
得 ,2,1,0 )(2)21(
==+-n e
C t T t
a l
n n n π
∴ t
a l
n n
n n n n e xC l n D t T x X t x u 2)21(
21sin
)()(),(ππ+-+
==
x l
n e
a t
a l
n n ππ21
sin
2)21(
+=+-
Step4:叠加,原方程的解为:
∑∑∞
=∞
=+
-+=
=
1
1
)21(
21sin
),(),(2n n t
a l
n n
n
x l
n e
a
t x u
t x u ππ
∵ )()0,(x x u ?=
∴ ∑∞
=+
=
1
21s i n
)(n n
x l
n a
x π? 其中 2,1,0,,21
s i n
)(2
0=+
=
?n dx x l
n x l a l
n π?
2.今有一弦,其两端被钉子钉紧作自由振动,其初始位移为:
??
?≤≤-≤≤=2
1 )21 )(x x h(x o hx x ?
初速度为零,试求弦振动方程的解(其中h 为常数) 解:该问题的数字模型为
??
?
??====><<=0
)0,(),()0,(0
),(),0(,2 2x u x x u t l u t u o t x o u a u t xx tt ? 由本节推导知:
∑∞
=???
??
?+=
1sin sin cos ),(n n n x l n at l n b t l a n a t x u πππ 其中2=l
)()0,(x x u ?= x l
n a x h n ∑
∞
==
∴1
s i n
)(π?
?
=∴l
n xdx l
n x l
a 0
sin
)(2π?xdx n x ?
=
2
2
sin
)(π?
xdx n x xdx n x ?
?
+
=2
1
1
02sin
)(2
sin
)(π?π?
xdx n x h xdx n hx ?
?
-+
=
2
1
1
2sin
)2(2
sin
π
π
??
-??
? ??-+-
=2
11
2cos
)2(22cos 2x n d x n h x n xd n h o
πππ
π
??
?
?
?
?+--??????--
=?
?
2
1
1
2cos
12
2cos )2(22cos
1
2cos 2xdx n x n x n h
xdx n x n x n h π
πππ
ππ ???
??
?+
-??????-
-=122sin 202012sin 22x n n n h
x n n n h πππ
πππ ??? ??
--
=
2sin 0)(42
sin
)
(422
ππππn n h
n n h 2
sin )(82
ππn n h = ,...2,1=n 又 ∵ x l n t l a n l a n b t l a n l a n a t x u n n
n t π
ππππs i n c o s s i n ),(1???
???
????
??+??? ??-=
∑∞
=
∵ 0)0,(=x u t ,∴ ∑∞
=??
?
??=1
s i n 0h h
x l n l a
n b
ππ,∴ 0=n b ∴ 原方程的解为
x l
n t l
a n a t x u h n ππsin
cos
),(1
∑
∞
==
∑∞
==
1
2
2
sin
2
cos
2
sin
)
(8h x n t a n n n h
ππππ
3.求解下面的定解问题
(1)??
?
??====><<=0),(),0(0)0,(,sin 3)0,(0
,0,2t u t u x u x x u t x u a u t xx tt ππ
解:根据本节推可知(其中π=l )
∑
∞
=+=
1
sin
)sin
cos
(),(n n n x l
n t l
a n
b t l
a n a t x u πππ
nx nat b nat a
n n n
sin sin cos 1∑∞
=+=
∑∞
=+-=
1
sin ]cos sin [),(n n n
t nx nat na b nat naa
t x u
x x u sin 3)0,( = ,∑∞
==
∴1
sin sin 3 h n
nx a
x
2,3,n , 0a ,3n 1===∴a 0)0,(=x u t ,∑∞
===
∴1
0sinh 0 h n
x na b
, 1,2,n 0 ==∴n b
x at x at a t x u sin cos 3sin cos ),(1==∴
(2)???
??====><<=0
)0,(,)0,(0),(),0(0
,0 ,3
2x u x x u t u t u t x a a u t x xx tt ππ 解:取π=l
Step1:分离变量:令),()(),(t T x X t x u =并代入齐次方程和齐次边界条件有:
)()()()(2
t T x X a t T x X ''=''
0)()()()0(==t T X t T X π
由于 0)(,0)(≠≠t T x X ∴ 上面方程变形为 λ=?
''='')
()()
()(2
t T a
t T x X x X 0)()0(='=πX X
整理有??
?='==-''0
)()0(0)()(πλX X x X x X
0)()(2
=-''t T a t T λ
Step2:求解特征值问题
?
?
?='==-''0)()0(0)()(l X X x X x X λ
由第1题知 22
)21(21+-=??
??
?
?
??+-=n l n n πλ )21sin(21
sin
)(x n C x l n C x X n n n +=+
=π, ,2,1,0=n Step3:将2)2
1
(+-=n n λ代入关于T (t )的微分方程
0)()2
1()(2
2=++''t T a n t T n n
其特征方程为 0)2
1(2
22
=+
+a n γ
,i n a )2
1( 2,1+
±=∴γ t n a B t n a A t T n n n )2
1(sin )2
1(cos )( +
++=∴,n=0,1,2,…
)()(),(t T x X t x u n n n =??
????
++++
=t n a B t n a A x n C n n n )21(sin )21(cos )21
sin( x n t n a b t n a a n n )21
s i n ()21
(s i n )21
(c o s +???
???+++=
Step4:叠加
∑∞
==
),(),(n n t x u t x u x n t n a n a a n n ??? ??+?
????
???? ??+??? ??+=
∑∞
=21sin 21sin 21cos 0 ∵ ∑∞
=??
???? ??+??? ??+-=
021sin 21),(n n t t n a n a a t x u x n t n a n a b n
??? ?
?
+?????? ?
?
+??? ??
++21sin
21cos 21
∵ 3
)0,(x x u =,∴ ∑∞
=??? ??
+=0
3
21s i n n n x n a x
所以
xdx n x xdx n x a n 2
12sin
221sin 2
3
3
+=
??? ?
?+=
?
?
π
π
ππ
x n d x n 2
1
2cos
1222 0 3+??? ??+-=
?ππ ??
?
?
??+-++-
=?
π
π
π 0
2
32
12cos
302
)12(cos )12(4xdx n x x n x n
x n d x n xdx n x n 2
12sin
)12(242
12cos
)12(12
2
2
2
++=
++=
?
?
π
π
π
π
??
?
???+-++=?
π
π
π 0
22
2
12s i n
202
12s i n )12(24
xdx n x x n x n
??
????++?+??? ??
++=?ππππ 0 22212c o s 122221s i n )12(24
dx n xd n n n ??
?????
?????????+-+++??? ??++=?
π
ππππ 0
22
2
1
2c o s 0212c o s 12421s i n )12(24
xdx
n x n x n n n ??????++?+-??? ?
?
++=
0212s i n 12212421s i n )12(24
22
ππππx n n n n n ?
???????? ??
++-??? ??++=
ππππ21s i n )12(821s i n )12(2422
n n n n πππ??? ?
?+?
????
?
+-
+=21s i n )12(8
)12(24
22
2
n n n , n =0, 1, 2,… 又∵ 0)0,(=x u t ,x n n a b n n
??? ??
+??? ??+=
∑∞
=21sin 2100
,∴ 0=n b n =0, 1, 2,… ∴ 原方程的解为
x n t n a a
t x u n n
??? ?
?
+??? ??+=
∑∞
=21s i n 21c o s ),(0
t n a n n n n ??? ??+??? ??+?
????
?
+-
+=
∑
∞
=21c o s 21s i n )12(8
)12(24
22
2
πππx n ??
? ??+21s i n (3)??
?
??==≤≤=><<=0),1(),0(10 ,0)0,(0
,10 ,4N
t u t u x x u t x u u xx t
解:令 0),(),(N t x v t x u += 则:??
?
??==-=><<=0),1(),0()0,(0
,10 ,40t v t v N x v t x v v xx t
由于 a=2,l =1,根据第1题结果有:
x l
n e
a
t x v t l a n n n
ππsin
),(2
1
??
? ??-∞
=∑=
()x n e
a
t
n n n
ππsin 221
-∞
=∑=
其中 x d x l
n x l
a l
n π?s i n
)(2 0
?
=
??-=-=1
01
0s i n 21
s i n
2x d x n N x d x n N ππ
[]1cos 201)cos(1200-=
??? ??
--=ππππn n N x n n N = ,2,1,0 201
2)12(40
=??
???=+=+-k k
n k n k N π ∴ ()x n e
a t x v t
n n n ππsin ),(221
-∞
=∑
=
()x k e
k N t
k k ππ
π)12sin()12(42
2
1240
++-=
+-∞
=∑
∴ ),(),(0t x v N t x u +=∑∞
=+-++-
=0
)12(40
0)12s i n ()
12(1
42
2k t
k x k e
k N N ππ
π
(4)??
?
??
??≤≤==≤≤=≤≤-=<<<<=+10 0),1(),0(10 0)1,(10 )1()0,(1
0,10 0y y u y u x x u x x x x u y x u u yy xx
解:step1:分离变量,令)()(),(y Y x X y x u =,并代入齐次方程和齐次边界 条件中有:
0)()()()(=''+''y Y x X y Y x X 0)()1()()0(=+y Y X y Y X
∵ 0)(,0)(≠≠y Y x X 上面方程进行整理有:
λ?''-
='')
()()
()(y Y y Y x X x X
0)1()0(==X X
即 ?
??===-''0)1()0(0)()(X X x X x X λ
由第1题知特征值为()22
π
πλn l n n -=??
? ??-=
特征函数为x n C l
x n C x X n n n ππsin sin )(==,n =1,2,……
step3:将()2π
λn n -=代入至于Y (y )的常微分方程有
0)()()(2
=-''y Y n y Y π
其特征方程为0)(2
2=-πn r ∴ πn r ±=2,1
∴ y n n y n n n e B e A y Y ππ-+=)( n =1,2,……
∴ []y n n y n n n n n n e B e A x n C y Y x X y x u πππ-+==s i n )()(),([]x n e b e a y n n y n n πππs i n -+= step4叠加:原方程的解为:
[]∑∑∞
=-∞
=+=
=
1
1
sin ),(),(n y
n n y
n n
n n
x n e
b e
a
y x u
y x u πππ
∵ )1()0,(-=x x x u ,∴ []∑∞
+=
-1
s i n )1(n n n
x n b a
x x π
∴ ()
xdx n x x
l
b a l
n n πsin 2 0
2
?-=
+()
xdx n x x
l
πsin 2 0
2
?
-=()
x n d x x
n l
ππ
cos 2 0
2
?--
=
()??
????--
--=?x d x n x x n x x n l πππc o s 1201c o s )(2 0 2
()x n d x n l
ππs i n 12)(2
2
?-=
??
??
??--=
?x d x n x n x n l
πππ 0
2s i n 01
s i n )12()(4
[]
1)1()(401c o s )(433--==n
n x n n πππ 即 []
1)
1()
(43
--=
+n
n n n b a π (1)
又 ∵ 0)1,(=x u 即[]∑∞
=-+=
1
sin 0n n n n n
x n e
b e a
ππ
π
∴ 0=+-π
πn n n n e
b e a (2) 联立(1)、(2)有 3
3
)
()1(1)(2)
(1)1()(2πππππ
π
n n sh e
n n sh e
a n
n n
n n --=
---
=--
3
3
)
()
1(1)(2)
(1
)1()(2πππππ
π
n n sh e
n n sh e
b n
n n
n n --=
---
=
∴ 原方程的解为
[][]
∑∞
=--??
?
????---?--=13
3sin )()()1(12)()()1(12),(n y n n n y n n n x n e n n sh e e n n sh e y x u πππππππππ [
][]∑∞
=---?--=
1
3
sin )
()
()
1(12n y
n n y
n n n
x n e
e
e
e
n sh n πππππ
ππ
[
]
∑
∞
=?---=1
3
sin )1()
()()
1(14n n
x n y shn n sh n ππππ
4.求阻尼波动问题的解
??
?
??≥==≤≤==><<=+0 0),(),0(0 )()0,(,0)0,(0,0 ,2t t l u t u l x x g x u x u t l x u c au u t xx t tt 解:Step1:分离变量,令)()(),(t T x X t x u =,并代入齐次方程和齐次边界条件中有
)()()()()()(2
t T x X c t T x aX t T x X ''='+''
0)()()()0(==t T l x t T X
由于0)(≠x X ,0)(≠t T ,于是上面方程变为:
λ )
()
()()
()(2
?'+''=
''t T c t T a t T x X x X
0)()0(==l X X
整理得下面的常微分方程有:
??
?===-'+''=-''
0)()0(0
)()()( 0)()(2l X X t T c t T a t T x X x X λλ Step2:求解下面的特征值问题
?
?
?===-'' 0)()0(0
)()(l X X x X x X λ 由第1题的结果有特征值为2
)(l
n n πλ-=
特征函数为 ,2,1 sin )(==n x l
n c x X n n π
Step3:将2
)(
l
n n πλ-=代入关于)(t T 的常数分方程有:
0)()(
)()(2
=+'+''t T l
c n t T a t T n n π
上式是关于)(t T 的二阶常数系线性齐次常微分方程,其特征方程为
0)(
2
2
=++l
c n ar r
π,∴2
)
2(2
2
2,1l
c n a
a r π-±
-=
当2
2)2(
l
c n a π>时 t
r n t
r n n e
B e A t T 21)(+=
其中2
)
2(2
2
1l
c n a
a r π-±-=
,2
)
2(2
2
2l
c n a a r π-±
-=
当2
2
)2(
l
c n a π=时 at n n n e
t B A t T 2
1
)()(-+=
当2
2)2(
l
c n a π<时 ]2
)
2(s i n
2)
2(c o s
[)(2
2
2
2
2
1t l
c m a
B t l
c n a
A e
t T n n at
n ππ-+-=-
Step4:叠加 当2
2)2(
l
c n a π>时
[
]∑
∞
=+=
1
21e
e
sin
),(n t
r t
r Bn An x l
n Cn t x u π[]x l
n bn a
n t
r t
r n πsin
e
e
1
21∑∞
=+=
∵0)0,(=x u ∴∑∞
==+0
0sin
)(n n n l
x n b a π
∴0=+n n b a (1) 又∵)()0,(x g x u =
而∑
∞
=+=
1
21sin
][),(21n t
r n t
r n t x l
n e
r b e
r a t x u π,∴∑
∞
=+=
1
21sin
][)(n n n x l
n r b r a x g π
∴?
=
+l
n n x l
n x g l
r b r a 0
21sin
)(2π (2)
(1)两边同乘2r 有022=+n n b r a r (3) (2)-(3)得:?
?
?
-=
?-=
l
l
n xdx l
n x g l
l
c n a
x l
n x g l
r r a 0
2
2
2
1sin
)(2)
2(1sin
)(21πππ
(1)两边同乘以1r 有011=+n n b r a r (4) (2)-(4)得:?
?
?--
=?
-=
l
l
n xdx l
n x g l
l
c n a
xdx l
n x g l
r r b 0
2
2
1
2sin
)(2)
2(1sin
)(21πππ
∴∑
?
?
∞
=?--
?-=
1 0
2
2
2
2
sin
])(2
)
2(
1sin
)(2
)
2(
1[
),(21n l
l
t
r t
r x
l
n e
xdx l
n sim
x g l l c n a xe
l
n x g l l
c n a t x u πππππ
∑
?
∞
=-??-=
1
2
2
sin ][sin
)(2
)
2(
121n t
r t
r l
x
l n e
e
xdx l
n x g l l c n a πππ
∑
?
∞
=--
---??-=
1
)2(
2
1)2(
2
1
2
1 0
2
2
sin
][sin
)(2)
2(
1
2
2
2
2
n t
l
c n a t
l
c n a at l
l
x n e e
e
xdx l
n x g l
l
c n a πππππ
∑?
∞
=---?
?=
1
2
2
2
2
2
1 0
sin
)2(
2
1)
2(
2sin
)(2
n at l
x
l
n t l
c n a sh l
c n a e
xdx l
n x g l ππππ
令?
=
l
n x l
n x g l
D 0
sin
)(2π , 2
)2(
l
c n πγ=
则当γ>2a 时,原方程的解为
∑
∞
==
1
sin
)(),(n n x l
n t T D t x u π
其中?
=l
n xdx l
n x g l
D 0
sin
)(2π,t a
sh a
e t T at
2
2)(2
2
2
1γ
γ
--=
-
同理可得
γπ==2
2
)2(
l
c n a
时,原方程有解:∑∞
==
1
sin
)(),(n n x l
n t T
D
t x u π
其中 ?
=
l
n xdx l
n x g l
D 0
sin
)(2π,at te
t T 2
1)(-=
当γπ=<2
2
)2(
l
c n a 时,∑
∞
==
1
)sin(
)(),(n n x l
n t T D t x u π
其中 ?
=
l
n xdx l
n x g l
D 0
sin
)(2π, t a
a
e
t T at
2
sin
2)(2
2
2
1--=
-γγ
5.均匀细杆长为l ,在0=x 固定,而另一端受着一个沿杆长方向的力Q ,如果在开始一瞬间,突然停止这个力的作用,求杆的纵振动。
解:假设细杆的左端点固定,于是该问题的数字桂型为:
?
????
??===??===><<=??
0)0,()0,(
0),(),0( 0,0, 0
0 2x u E Qx dx E Q dx x u x u t l u t u t l x u a u t x x x xx tt σσ
其中E 为细杆的杨氏模量,σ为细杆的横截面积
Step1:分离变量:令)()(),(t T x X t x u =,并代入齐次方程和齐次边界条件中有:
)()()()(2
t T x X a t T x X ''='',0)()()()0(='=t T l X t T X
因为0)(≠x X ,0)(≠t T ,故上面式子可变为
λ=''='')
()()
()(2
t T a t T x X t X (待常数)
0)()0(='=l X X
得常微分方程有:
0)()(=-''x X x X λ
0)()0(='=l X X 0)()(2
=-''t T a t T λ
Step2:求解特征值问题
?
?
?='==-''0)()0(0
)()(l X X x X x X λ 由第1题的结果可知特征值为2)21(πλl
n n +
-=
特征函数
x l n C x X n n π21
sin )(+
= n =0,1,2,……
Step3:将2)21
(
πλl n n +
-=代入关于
)(t T 的常微分方程求)(t T 。 0)()21
(
)(2=+
+''t T a l
n t T n n π
其特征根为 i a l
n r π212,1+
±
=
∴at l n B at l
n A t T n
n n ππ21
sin
21
cos
)(+++
= n =0,1,2,……
∴x l n at l n B at l n a t T x X t x u n n n n n πππ21sin
]21sin
21
cos
[)()(),(++
++
==
n =0,1,2,……
Step4:叠加:原方程的解为:
∑∑∞
=∞
=++++
=
=
21sin
]21sin
21
cos
[),(),(n n n
n
n
x
l
n at l
n b at l
n a
t x u
t x u πππ
其中n n b a ,得求
∵σ
E Qx x u =
)0,( ∴
∑∞
=+
=
21sin
n n
x
l
n a
E Qx πσ
∴?
+=
l
n xdx l n E Qx l
a 0
21
sin
2πσ
?
+
-=+=
l
n
n l E Q
xdx l n x lE Q 0
2
2)2
1
()1(221
sin
2π
σ
πσ
n =0,1,2,……
∵∑∞
=+?+++++
?
-=
21sin
]21
cos
2121
sin
21)1([),(n n
n
t x l
n at l
n a l
n b
at l
n a l
n a
t x u πππππ
∵0),(==t x u t
即 x l
n a l
n b
n n
ππ21
sin
2100
++
=
∑∞
= ∴0=n
b n =0,1,2,……
∴原方程的解为
x
l
n at l
n n E lQ
t x u n n ππσ
π21sin
21cos
)
12()
1(8),(0
2
2
+++-=
∑
∞
=
6.长为2l 的均匀细杆,被作用在两端的压力压缩式)1(2ε-l ,在t =0时,把这个载荷移去,试证,若x =0是杆的中点,则在t 时刻,坐标为x 的杆的截面位移),(t x u 由下式确定:
at l
n x l
n n l
t x u n n 2)12(cos
2)12(sin
)
12()
1(8),(0
2
1
2
π
ππ
ε+++-=
∑
∞
=+
证明:参见第1章,该问题的数学模型为:
?
????
??=-=-=≥==><<=
0)0,( )0,(0
0),(),0(0,0 ,2x u x x l
l x u t t l u t u t l x u a u t x xx tt εε
Step1:分离变量:令)()(),(t T x X t x u =并代入齐次方程和齐次边界条件有,
)()()()(2
t T x X a t T x X ''=''
0)()()()0(='=t T l X t T X
由于0)(,0)(≠≠t T x X ,故上面的方程可变为
λ=''=
'')
()()
()(2
t T a t T x X x X
0)()0(='=l X X
整理即得两个常微分方程
??
?=-''='==-''0)()( 0
)()0(0)()(2
t T a t T l X X x X x X λλ Step2:求解特征值问题
??
?='==-''
0)()0(0
)()(l X X x X x X λ 由第1题的结果有特征值为2)21(
πλl
n n +
-=
特征函数为 x l
x
n C x X n n ππ21s i n
)(+
= n =0,1,2,……
Step3:将2)21
(
πλl
n n +
-=代入关于)(t T 的常微分方程求)(t T .0)()21
(
)(2//
=+
+t T a l
n t T n
n π
at l n B at l n A t T n
n n ππ21sin
21
cos
)(+++
=∴
x l n at l n b at l
n a t T x x t x U n n n n n πππ21sin
]21sin
21
cos
[)()(),(++++
==∴
step4:叠加:原方程的解为
∑∑∞
=∞
=++++
=
=
21sin
]21sin
21
cos
[),(),(n n n
n
n
x l
n at l
n b at l
n a
t x u
t x u πππ
0)0,(=x u t 得 0=n b
x x u ε-=)0,( 即∑∞
=+
=
-0
21sin
n n
x
l
n a
x πε
?
+
-=
∴l
n xdx l
n x l
a 0
21
sin
2πε?
+-=l
xdx l
n x l 0
212sin
2πε
?
++?
=
l
x l
n xd n l l
212cos
)12(22ππ
ε]212cos
212cos
[)12(40
0?+-++=
l
l
x l
n x l
n x n πππ
ε
]212sin
)12(2)2
1cos([1
240l
x l
n n l n l n ππ
πε++-
+
?+=
1
2
2
2
2
)
1()12(8])2
1[sin()12(8+-+=
+
+-
=n n l n n l π
εππ
ε
∴原方程的解为
x l
n at l
n n l
t x u n πππ
ε21sin
21cos
)
12(18),(0
2
2
++
+=
∑
∞
=
x l
n at l
n n l
n πππ
ε2)12(sin
2)12(cos
)
12(180
2
2
+++=
∑
∞
= [证毕]
7.求下列高维高波动的解
?????????
?
?========><<=0
),1,,(),0,,(0),,1,(),,0,(0),,,1(),,,0(0
)0,,,(sin sin sin )0,,,(0,1,,0,2t y x u t y x u t z x u t z x u t z y u t z y u z y x u z y x z y x u t z y x u a u t xx tt πππ 解:step1:分离变量:令)()()()(),,,(t T z Z y Y x X t z y x u = 代入齐次微分方程有:
][)(2
T Z XY ZT Y X YZT X a t T XYZ ''+''+''=''
两边同除以)()()()(2
t T z Z y Y x X a 有:
Z
Z Y Y X
X T a T ''+
''+
''=
''2
要使上式成立,上式右边每一项必须是常数
令
,,
,
γβα=''=''=''Z
Z Y Y X
X 且μγβα=++
于是有4个常微分方程
0)()(=-''x X x X α 0)()(=-''y Y y Y β 0)()(=-''z Z z Z γ 0)()(2
=-''t T a t T μ
<<电磁场与电磁波>>读书报告 姓 名: 学 院: 学 号: 专 业: 题 目:分离变量法在求静态场的解的应用 成 绩: 二〇一四年四月 Xxx 工程学院 电子工程类
一.引言 分离变量法是在数学物理方法中应用最广泛的一种方法。在求解电磁场与电磁波的分布型问题和边值型问题有很重要的应用。分布型问题是指已知场源(电荷分布、电流分布)直接计算空间各点和位函数。而边值型问题是指已知空间某给定区域的场源分布和该区域边界面上的位函数(或其法向导数),求场内位函数的分布。求解这两类问题可以归结为在给定边界条件下求解拉普拉斯方程或泊松方程,即求解边值问题。这类问题的解法,例如镜像法,分离变量法,复变函数法,格林函数法和有限差分法,都是很常用的解法。这里仅对在直角坐标系情况下的分离变量法作简单介绍。 二.内容 1.分离变量法的特点: 分离变量法是指把一个多变量的函数表示成几个单变量函数乘积,从而将偏微分方程分离为几个带分离常数的常微分方程的方法,属于解析法的一种。它要求要求所给边界与一个适当的坐标系的坐标面重合.在此坐标系中,待求偏微分方程的解可表示成三个函数的乘积,每一函数仅是一个坐标的函数。我们仅讨论直角坐标系中的分离变量法. 2.推导过程: 直角坐标系中的拉普拉斯方程: 222 222 0 x y z ??? ??? ++=??? 我们假设是三个函数的乘积,即
(,,)()()()x y z X x Y y Z z ?= 其中X 只是x 的函数,同时Y 是y 的函数Z 是z 的函数,将上式带入拉普拉斯方程,得 然后上式同时除以XYZ ,得 0X Y Z X Y Z '''''' ++= 上式成立的唯一条件是三项中每一项都是常数,故可分解为下列三个方程: 即 α,β,γ为分离常数,都是待定常数,与边值有关但不能全为实数或全为虚数 。 由上式得2220αβγ++=,下面以X ”/X =α2式为例,说明X 的形式与α的关系 当α2=0时,则 当α2 <0时,令α=jk x (k x 为正实数),则 或 当α2 >0时,令α=k x ,则 或 a ,b ,c ,d 为积分常数,由边界条件决定Y(y)Z(z)的解和X(x)类似。 3解题步骤 1,2λα =±00 ()X x a x b =+12()x x jk x jk x X x b e b e -=+12()sin cos x x X x a k x a k x =+12()x x k x k x X x d e d e -=+12() s x x X x c hk x c chk x =+
第三章行波法与积分变换法 」 分离变量法,它是求解有限区域内定解问题常用的一种方法。 J 行波法,是一种针对无界域的一维波动方程的求解方法。 」 积分变换法,一个无界域上不受方程类型限制的方法。 作如下代换; X at, X at 利用复合函数求导法则可得 同理可得 2 a 2(£ 代入(1)可得 =0o u(x,t) F( ) G( ) F(X at) G(X at) 这里F,G 为二阶连续可微的函数。再由初始条件可知 F(X ) G(X ) (X ), aF (X ) aG (X ) (X ). X 2 u -2 )(」 2 2」 2 u ~2 先对求积分,再对 求积分,可得u(X,t)d 的一般形式 § 3.1 一维波动方程的达朗贝尔 (D 'alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: 2 u 下 u 2 2 u a 2 , X (X), u 0, (1) (X ),- (2) 2 ■4), (3)
由(3)第二式积分可得 1 X F(x) G(x) - 0 (t)dt C , a 0 利用(3)第一式可得 所以,我们有 1 1 x at u(x,t) [ (x at) (x at)] (t)dt 2 2a x at 此式称为无限弦长自由振动的达朗贝尔公式。 二、 特征方程、特征线及其应用 考虑一般的二阶偏微分方程 AU xx 2BU xy CU yy DU x EU y Fu 0 称下常微分方程为其特征方程 A(dy)2 2Bdxdy C(dx)2 0。 由前面讨论知道,直线x at 常数为波动方程对应特征方程的积分曲线, 称为特征线。已知,左行波F(x at)在特征线x at G 上取值为常数值F(CJ , 右行波G(x at)在特征线x at C 2上取值为常数值G(C 2),且这两个值随着特 征线的移动而变化,实际上,波是沿着特征线方向传播的。称变换( 2)为特征 变换,因此行波法又称特征线法。 注:此方法可以推广的其他类型的问题。 三、 公式的物理意义 由 U(x,t) F (x at) G(x at) 其中F(x at)表示一个沿x 轴负方向传播的行波,G(x at)表示一个沿x 轴正方 向传播的行波。达朗贝尔公式表明:弦上的任意扰动总是以行波形式分别向两个 方向传播出去,其传播速度为a 。因此此法称为行波法。 四、 依赖区间、决定区域、影响区域 F(x) 1 2(X ) 2a (t)dt G(x) (x) 1 x 2a o (t)dt (4)
第三章 分离变量法 3。2 基础训练 3.2.1 例题分析 例1 解下列定解问题: ???? ?????=??-==??=><?=??====0,20,00,002 002 2222t t l x x t u lx x u x u u t l x x u a t u (1) 解:分离变量,即令 (,)()()u x t X x T t = (2) 代入方程((1)中第一式),得 0)()(2=+''t T a t T λ (3) 0)()(=+''x X x X λ (4) 其中λ为分离常数。(2)式代入边界条件((1)中第二式),得 0)()0(='=l X X (5) 相应的本证值问题为求 ?? ?='==+''0 )()0(0 )()(l X X x X x X λ (6) 的非零解.下面针对λ的取值情况进行讨论: (1)当0λ<时,(6)式中方程的通解是 ()X x Ae =+ (7)
其中A ,B 为积分常数,(7)代入(6)中边界条件,得 00 A B Ae +=???-+=?? (8) 由(8)得A=B=0,得X (x )=0,为平凡解,故不可能有0λ<。 (2) 当0λ=时,(6)式中方程的通解是 ()X x Ax B =+ 由边界条件得A=B=0,得X (x )=0,为平凡解,故也不可能有0λ=。 (3)当 02 >=βλ时,上述固有值问题有非零解.此时式(6)的通解为 x B x A x X ββsin cos )(+= 代入条件(6)中边界条件,得 0cos ,0==l B A β 由于 0≠B ,故 0cos =l β,即 ),2,1,0(212Λ=+= n l n πβ 从而得到一系列固有值与固有函数 2 2 24)12(l n n πλ+= ),2,1,0(2)12(sin )(Λ=+=n x l n B x X n n π 与这些固有值相对应的方程(3)的通解为 ),2,1,0(2)12(sin 2)12(cos )(Λ=+'++'=n t l a n D t l a n C t T n n n ππ 于是,所求定解问题的解可表示为 x l n t l a n D t l a n C t x u n n n 2)12(sin 2)12(sin 2)12(cos ),(0πππ+??? ? ? +++=∑∞ = 利用初始条件确定其中的任意常数n n D C ,,得 0=n D
第二章 分离变量法 偏微分方程定解问题常用解法,分离变量法。 解常微分方程定解问题时,通常总是先求出微分方程的特解,由线性无关的特解叠加出通解,而后用定解条件定出叠加系数 一阶线性偏微分方程的求解问题,基本方法也是转化为一阶线性常微分方程组的求解问题 对于二阶以及更高阶的偏微分方程定解问题,情况有些不同:即使可以先求出通解,由于通解中含有待定函数,一般来说,很难直接根据定解条件定出,因此,通常的办法就是把它转化为常微分方程问题 §2.1 有界弦的自由振动 什么是分离变量法?使用分离变量法应具备那些条件? 下面通过两端固定的弦的自由振动问题来说明。 定解问题:考虑长为l ,两端固定的弦的自由振动,其数理方程及定解条件为 .0 ),(u ),(u 0, ,0u ,0u 0, l,0 ,0 t 0022 222l x x x t t x x u a t u t t l x x ≤≤==>==><?=??====ψ? 分析: 1. 方程和边界条件都是齐次的,求这样的问题可用叠加原理。 2. 我们知道,在解常微分方程定解问题时,通常总是先求出微分方程的特解, 由线性无关的特解叠加出通解,而后用定解条件定出叠加系数。 启发:能否运用类似求常微分方程定解问题的方法求偏微分方程?也既是能否先找出满足齐次方程及齐次边界条件的足够多的特解,再用其作线性组合使其满足初始条件。 由分析,我们现在试求方程的变量分离形式: )()(),(t T x X t x u = 的非零解。 将),(t x u 代入方程,可得 ) () ()()()()()()(2''''' '2 ' 'x T a x T x X x X t T x X a t T x X = ?= 此式中,左端是关于x 的函数,右端是关于t 的函数。因此,左端和右端相等,就必须等于一个与t x ,无关的常数。设为λ-,则有
2019年数学物理方程-第二章分离变量法.doc
第二章 分离变量法 分离变量法是求解偏微分方程定解问题最常用的方法之一,它和积分变换 法一起统称为Fourier 方法. 分离变量法的本质是把偏微分方程定解问题通过变量分离,转化为一个所谓的特征值问题和一个常微分方程的定解问题,并把原定解问题的解表示成按特征函数展开的级数形式. 本章介绍两个自变量的分离变量法,更多变量的情形放在其他章节中专门讨论. §2?1 特征值问题 2.1.1 矩阵特征值问题 在线性代数中,我们已学过线性变换的特征值问题. 设A 为一n 阶实矩阵,A 可视为n R 到自身的线性变换。该变换的特征值问题(eigenvalue problem )即是求方程: ,n Ax x x R λ=∈, (1.1) 的非零解,其中C λ∈为待定常数. 如果对某个λ,问题(1.1)有非零解n x R λ∈,则λ就称为矩阵A 的特征值(eigenvalue),相应的n x R λ∈称为矩阵A 的特征向量(eigenvector). 一般来讲,特征值问题(1.1)有不多于n 个相异的特征值和线性无关的特征向量. 但可证明: 任一n 阶矩阵都有n 个线性无关的广义特征向量,以此n 个线性无关的广义特征向量作为n R 的一组新基,矩阵就能够化为Jordan 标准型. 若A 为一n 阶实对称矩阵,在线性代数中有一个重要结果,即存在一个正交矩阵T 使得 1T AT D -=, (1.2) 其中D =diag 12(,,...,)n λλλ为实对角阵. 设12[ ... ]n T T T T =,i T 为矩阵T 的第i 列向量(1)i n ≤≤,则式(1.2)可写为如下形式 1212 [ ... ][ ... ]n n A T T T T T T D =, 或 , 1.i i i A T T i n λ=≤≤ (1.3) 上式说明,正交矩阵T 的每一列都是实对称矩阵A 的特征向量,并且这n 个特征向量是相互正交的. 由于此结论在一定意义下具有普遍性,我们以定理的形式给出. 定理1.1 设A 为一n 阶实对称矩阵,考虑以下特征值问题 ,n Ax x x R λ=∈, 则A 的所有特征值为实数,且存在n 个特征向量,1i T i n ≤≤,它们是相互正交的(正交性orthogonality ),可做为n R 的一组基(完备性completeness ). 特征值问题在线性问题求解中具有重要的意义,下面举例说明之. 为简单起见,在下面两个例子中取A 为n 阶非奇异实矩阵,故A 的所有特征值非零,并且假设A 有n 个线性无关的特征向量,i T 相应的特征值为, 1i i n λ≤≤. 例1.1 设n b R ∈,求解线性方程组 Ax b =. 解 由于向量组{1}i T i n ≤≤线性无关,故可做为n R 的一组基. 将,x b 按此
§2.2 有限杆上的热传导 定解问题:一均匀细杆,长为l ,两端坐标为l x x == ,0。杆的侧面绝热,且在端点0=x 处温度为零,而在l x = 处杆的热量自由发散到周围温度为0的介质中。初始温度为)(x ?,求杆上的温度变化情况,即考虑下定解问题: .0 ),(u 0, ,0hu ,0u 0, l,0 ,0002 2 2l x x t x u t x x u a t u t l x x ≤≤=>=+??=><<=??-??===? 仍用分离变量法求解。此定解问题的边界条件为第三类边界条件。类似§2.1中步骤,设)()(),(t T x X t x u =,代入上面的方程可得 ?????=+=+?-==. 0)()(,0)()() ()()()( 2 ' '22'2 2'''x X x X t T a t T x T a x T x X x X βββ 从而可得通解 x B x A x X ββsin cos )(+= 由边界条件知 .0)()(,0)0('=+=l hX l X X 从而 ?? ???-=?=+=.tan 0sin cos , 0h l l h l A βββββ 令 αγ γαβγ=?- ==tan 1 ,hl l 上方程的解可以看作曲线γtan 1=y ,αγ=2y 交点的横坐标,显然他们有无穷多个,于是方程有无穷多个根。用下符号表示其无穷多个正根 ,,21n γγγ 于是得到特征值问题的无穷个特征值
1,2,3...) (n ,2 2 2== l n n γβ 及相应的特征函数 x B x X n n n βsin )(= 再由方程0)()(22'=+t T a t T β, 可得 t a n n n e A t T 2 2)(β-=, 从而我们得到满足边界条件的一组特解 x e C t x u n t a n n n ββsin ),(2 2-= 由于方程和边界条件是齐次的,所以 ∑∞ =-=1 sin ),(2 2n n t a n x e C t x u n ββ 仍满足此方程和边界条件。 下面研究一下其是否满足初始条件。 )(sin 1 x x C n n n ?β=∑∞ = 可以证明}{sin x n β在区域[0,l]上具有正交性,即 ?≠=l m n xdx x 0 n m ,0sin sin ββ 证明: ) )((sin cos cos sin ))((2)sin()()sin()( ) (2)sin()(2)sin( ))cos()(cos(2 1sin sin 00=+--- =+-+---+=++- --=--+- =??m n m n m n n m n m m n m n m n m n m n m n m n m n m n m n l m n m n l m n l l l l l l l l dx x x xdx x ββββββββββββββββββββββββββββββββββββ 完成。 令 ?=l n n n xdx x L 0 ,sin sin ββ 于是, ?= l n n n xdx x L C 0 sin )(1β ?
分离变量法 分离变量法是近年来发展较快的思想方法之一.高考数学试题中,求参数的范围常常与分类讨论、方程的根与零点等基本思想方法相联系.其中与二次函数相关的充分体现数形结合及分类思想方法的题目最为常见.与二次函数有关的求解参数的题目, 相当一部分题目都可以避开二次函数,使用分离变量,使得做题的正确率大大提高.随着分离变量的广泛使用,越来越多的压轴题都需要使用该思想方法. 分离变量法:是通过将两个变量构成的不等式(方程)变形到不等号(等号)两端,使两端变量各自相同,解决有关不等式恒成立、不等式存在(有)解和方程有解中参数取值范围的一种方法.两个变量,其中一个范围已知,另一个范围未知. 解决问题的关键: 分离变量之后将问题转化为求函数的最值或值域的问题.分离变量后,对于不同问题我们有不同的理论依据可以遵循.以下定理均为已知x 的范围,求a 的范围: 定理1 不等式()()f x g a ≥恒成立?[]min ()()f x g a ≥(求解()f x 的最小值);不等式 ()()f x g a ≤恒成立?[]max ()()f x g a ≤(求解()f x 的最大值). 定理2 不等式()()f x g a ≥存在解?[]max ()()f x g a ≥(求解()f x 的最大值);不等式 ()()f x g a ≤存在解?[]min ()()f x g a ≤(即求解()f x 的最小值). 定理3 方程()()f x g a =有解?()g a 的范围=()f x 的值域(求解()f x 的值域). 解决问题时需要注意:(1)确定问题是恒成立、存在、方程有解中的哪一个;(2)确定是求最大值、最小值还是值域. 再现性题组: 1、已知当x ∈R 时,不等式a+cos2x<5-4sinx 恒成立,求实数a 的取值范围。 2、若f(x)=2 33x x --在[1,4]x ∈-上有()21f x x a ≥+-恒成立,求a 的取值范围。 3、若f(x)=233x x --在[1,4]x ∈-上有2 ()251f x x a a ≥+--恒成立,求a 的取值范围。 4、若方程42210x x a -+=有解,请求a 的取值范围 5、已知32 11132 y x ax x = -++是(0,)+∞上的单调递增函数,则a 的取值范围是( ) .0A a <.22B a -≤≤.2C a <.2D a ≤ 6、求使不等式],0[,cos sin π∈->x x x a 恒成立的实数a 的范围。 再现性题组答案: 1、解:原不等式4sin cos 25x x a ?+<-+当x ∈R 时,不等式a+cos2x<5-4sinx 恒成立max a+5>(4sinx+cos2x)?-,设f(x)=4sinx+cos2x 则 22f(x)= 4sinx+cos2x=2sin x+4sinx+1=2(sinx 1)+3 --- ∴a+5>3a<2-∴
第三章 行波法与积分变换法 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ? ? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 22 2 2 2 22 2))((,ηηξξηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2 2222222ηηξξ ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对η求积分,再对ξ求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知
). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3) 由(3)第二式积分可得 C dt t a x G x F x += -?0)(1)()(ψ, 利用(3)第一式可得 .2 )(21)(21)(,2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 二、特征方程、特征线及其应用 考虑一般的二阶偏微分方程 02=+++++Fu Eu Du Cu Bu Au y x yy xy xx 称下常微分方程为其特征方程 0)(2)(22=+-dx C Bdxdy dy A 。 由前面讨论知道,直线常数=±at x 为波动方程对应特征方程的积分曲线,称为特征线。已知,左行波)(at x F +在特征线1C at x =+上取值为常数值)(1C F ,右行波)(at x G -在特征线2C at x =-上取值为常数值)(2C G ,且这两个值随着特征线的移动而变化,实际上,波是沿着特征线方向传播的。称变换(2)为特征变换,因此行波法又称特征线法。 注:此方法可以推广的其他类型的问题。 三、公式的物理意义 由 )()(),(at x G at x F t x u -++= 其中)(at x F +表示一个沿x 轴负方向传播的行波, )(at x G -表示一个沿x 轴正方向传播的行波。达朗贝尔公式表明:弦上的任意扰动总是以行波形式分别向两个 方向传播出去,其传播速度为a 。因此此法称为行波法。
七、分离变量法 分离变量法是近年来发展较快的思想方法之一.高考数学试题中,求参数的范围常常与分类讨论、方程的根与零点等基本思想方法相联系.其中与二次函数相关的充分体现数形结合及分类思想方法的题目最为常见.与二次函数有关的求解参数的题目, 相当一部分题目都可以避开二次函数,使用分离变量,使得做题的正确率大大提高.随着分离变量的广泛使用,越来越多的压轴题都需要使用该思想方法. 分离变量法:是通过将两个变量构成的不等式(方程)变形到不等号(等号)两端,使两端变量各自相同,解决有关不等式恒成立、不等式存在(有)解和方程有解中参数取值范围的一种方法.两个变量,其中一个范围已知,另一个范围未知. 解决问题的关键: 分离变量之后将问题转化为求函数的最值或值域的问题.分离变量后,对于不同问题我们有不同的理论依据可以遵循.以下定理均为已知x 的范围,求a 的范围: 定理1 不等式()()f x g a ≥恒成立?[]min ()()f x g a ≥(求解()f x 的最小值);不等式()()f x g a ≤恒成立?[]max ()()f x g a ≤(求解()f x 的最大值). 定理2 不等式()()f x g a ≥存在解?[]max ()()f x g a ≥(求解()f x 的最大值);不等式()()f x g a ≤存在解?[]min ()()f x g a ≤(即求解()f x 的最小值). 定理3 方程()()f x g a =有解?()g a 的范围=()f x 的值域(求解()f x 的值域). 解决问题时需要注意:(1)确定问题是恒成立、存在、方程有解中的哪一个;(2)确定是求最大值、最小值还是值域. 再现性题组: 1、已知当x ∈R 时,不等式a+cos2x<5-4sinx 恒成立,求实数a 的取值范围。 2、若f(x)=2 33x x --在[1,4]x ∈-上有()21f x x a ≥+-恒成立,求a 的取值范围。 3、若f(x)=233x x --在[1,4]x ∈-上有2 ()251f x x a a ≥+--恒成立,求a 的取值范围。 4、若方程42210x x a -+= 有解,请求a 的取值范围 5、已知32 11132 y x ax x = -++是(0,)+∞上的单调递增函数,则a 的取值范围是( ) .0A a < .22B a -≤≤ .2C a < .2D a ≤ 6、求使不等式],0[,cos sin π∈->x x x a 恒成立的实数a 的范围。 再现性题组答案: 1、解:原不等式4sin cos 25x x a ?+<-+当x ∈R 时,不等式a+cos2x<5-4sinx 恒成立max a+5>(4sinx+cos2x)?-,设f(x)=4sinx+cos2x 则 22f(x)= 4sinx+cos2x=2sin x+4sinx+1=2(sinx 1)+3 --- ∴a+5>3a<2-∴
第五讲补充常微分方程求解相关知识。
第二章 分离变量法 偏微分方程定解问题常用解法,分离变量法。 解常微分方程定解问题时,通常总是先求出微分方程的特解,由线性无关的特解叠加出通解,而后用定解条件定出叠加系数 一阶线性偏微分方程的求解问题,基本方法也是转化为一阶线性常微分方程组的求解问题 对于二阶以及更高阶的偏微分方程定解问题,情况有些不同:即使可以先求出通解,由于通解中含有待定函数,一般来说,很难直接根据定解条件定出,因此,通常的办法就是把它转化为常微分方程问题 (第六讲) §2.1 有界弦的自由振动 什么是分离变量法?使用分离变量法应具备那些条件? 下面通过两端固定的弦的自由振动问题来说明。 定解问题:考虑长为l ,两端固定的弦的自由振动,其数理方程及定解条件为 .0 ),(u ),(u 0, ,0u ,0u 0, l,0 ,0 t 0022 222l x x x t t x x u a t u t t l x x ≤≤==>==><?=??====ψ? 分析: 1. 方程和边界条件都是齐次的,求这样的问题可用叠加原理。 2. 我们知道,在解常微分方程定解问题时,通常总是先求出微分方程的特解, 由线性无关的特解叠加出通解,而后用定解条件定出叠加系数。 启发:能否运用类似求常微分方程定解问题的方法求偏微分方程?也既是能否先找出满足齐次方程及齐次边界条件的足够多的特解,再用其作线性组合使其满足初始条件。 由分析,我们现在试求方程的变量分离形式: )()(),(t T x X t x u = 的非零解。 将),(t x u 代入方程,可得 ) () ()()()()()()(2''''' '2 ' 'x T a x T x X x X t T x X a t T x X = ?= 此式中,左端是关于x 的函数,右端是关于t 的函数。因此,左端和右端相等,就必须等于一个与t x ,无关的常数。设为λ-,则有
第三章 分离变量法 3。2 基础训练 3.2.1 例题分析 例1 解下列定解问题: ???? ?????=??-==??=><?=??====0,20,00,002 002 2222t t l x x t u lx x u x u u t l x x u a t u (1) 解:分离变量,即令 (,)()()u x t X x T t = (2) 代入方程((1)中第一式),得 0)()(2=+''t T a t T λ (3) 0)()(=+''x X x X λ (4) 其中λ为分离常数。(2)式代入边界条件((1)中第二式),得 0)()0(='=l X X (5) 相应的本证值问题为求 ? ? ?='==+''0)()0(0 )()(l X X x X x X λ (6) 的非零解.下面针对λ的取值情况进行讨论: (1)当0λ<时,(6)式中方程的通解是 ()X x Ae =+ (7) 其中A ,B 为积分常数,(7)代入(6)中边界条件,得 00 A B Ae +=???-+=?? (8) 由(8)得A=B=0,得X (x )=0,为平凡解,故不可能有0λ<。 (2) 当0λ=时,(6)式中方程的通解是 ()X x Ax B =+ 由边界条件得A=B=0,得X (x )=0,为平凡解,故也不可能有0λ=。 (3)当 02 >=βλ时,上述固有值问题有非零解.此时式(6)的通解为
x B x A x X ββsin cos )(+= 代入条件(6)中边界条件,得 0cos ,0==l B A β 由于 0≠B ,故 0cos =l β,即 ),2,1,0(21 2Λ=+= n l n πβ 从而得到一系列固有值与固有函数 2 2 24)12(l n n πλ+= ),2,1,0(2)12(sin )(Λ=+=n x l n B x X n n π 与这些固有值相对应的方程(3)的通解为 ),2,1,0(2)12(sin 2)12(cos )(Λ=+'++'=n t l a n D t l a n C t T n n n ππ 于是,所求定解问题的解可表示为 x l n t l a n D t l a n C t x u n n n 2)12(sin 2)12(sin 2)12(cos ),(0πππ+??? ? ? +++=∑∞ = 利用初始条件确定其中的任意常数n n D C ,,得 0=n D 3 32 02)12(322)12(sin )2(2ππ+- =+-=?n l xdx l n lx x l C l n 故所求的解为 x l n t l a n n l t x u n 2)12(sin 2)12(cos )12(132),(0 3 3 2 π ππ++?+- =∑∞ = 例2 演奏琵琶是把弦的某一点向旁边拨开一小段距离,然后放手任其自由振动。设弦 长为l ,被拨开的点在弦长的0 1 n (0n 为正整数)处,拨开距离为h ,试求解弦的振动,即求解定解问题
第三章 行波法与积分变换法 (第十三讲 ) 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ?? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 222222 22))((,ηηξξ ηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2222222 2ηηξξ ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对求积分,再对求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知 ). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3) 由(3)第二式积分可得 C dt t a x G x F x +=-?0 )(1)()(ψ, 利用(3)第一式可得
.2 )(21)(21)(,2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 例 求解柯西问题: ?????+∞≤≤-∞==+∞≤≤-∞>=-+==.,0,3,,0,03202 x u x u x y u u u y y y yy xy xx 解:其特征方程为 0)(32)(22=--dx dxdy dy 由此可得特征线方程为 d y x c y x =+=-3 因此作变换 ?? ?+=-=y x y x μξ, 3 从而可得 η ξ???u 2=0 从而有 )()3(),(y x G y x F y x u ++-= 由初始条件可得 )()3(3)()3(' ' 2=+-=+x G x F x x G x F 所以有 C x G x F =-)(3)3(, 从而可得 C x x G C x x F +=-=4 3)(4 9)3(2 2
第二章 分离变量法 分离变量法是求解偏微分方程定解问题最常用的方法之一,它和积分变换 法一起统称为Fourier 方法. 分离变量法的本质是把偏微分方程定解问题通过变量分离,转化为一个所谓的特征值问题和一个常微分方程的定解问题,并把原定解问题的解表示成按特征函数展开的级数形式. 本章介绍两个自变量的分离变量法,更多变量的情形放在其他章节中专门讨论. §2?1 特征值问题 2. 矩阵特征值问题 在线性代数中,我们已学过线性变换的特征值问题. 设A 为一n 阶实矩阵, A 可视为n R 到自身的线性变换。该变换的特征值问题(eigenvalue problem )即是求方程: ,n Ax x x R λ=∈, () 的非零解,其中C λ∈为待定常数. 如果对某个λ,问题()有非零解n x R λ∈,则λ就称为矩阵A 的特征值(eigenvalue),相应的n x R λ∈称为矩阵A 的特征向量(eigenvector). 一般来讲,特征值问题()有不多于n 个相异的特征值和线性无关的特征向量. 但可证明: 任一n 阶矩阵都有n 个线性无关的广义特征向量,以此n 个线性无关的广义特征向量作为n R 的一组新基,矩阵就能够化为 Jordan 标准型. 若A 为一n 阶实对称矩阵,在线性代数中有一个重要结果,即存在一个正交矩阵T 使得 1T AT D -=, () 其中D =diag 12(,,...,)n λλλ为实对角阵. 设12[ ... ]n T T T T =,i T 为矩阵T 的第i 列向量(1)i n ≤≤,则式()可写为如下形式 1212 [ ... ][ ... ]n n A T T T T T T D =, 或
第二章 分离变量法 §2.1 有界弦的自由振动 为了了解什么是分离变量法以及使用分离变量法应该具备什么条件,我们选取两端固定的弦的自由振动问题为例,通过具体地求解逐步回答这些问题。 讨论两端固定的弦的自由振动,归结求解下列定解问题: 22222000,0,0 (2.1)0,0,0 (2.2)(),(),0 (2.3)x x l t t u u a x l t t x u u t u u x x x l t ?ψ====???=<<>?????==>????==≤≤??? 这个定解问题的特点是:偏微分方程是线性齐次的,边界条件也是齐次的。求解这样的问题,可以运用叠加原理。我们知道,在求解常系数线性齐次常微分方程的初值问题时,是先求出足够多个特解(它们能构成通解),再利用叠加原理作这些特解的线性组合,使满足初始条件。这就启发我们,要解问题(2.1~2.3),先寻求齐次方程(2.1)的满足齐次边界条件(2.2)的足够多个具有简单形式(变量被分离的形式)的特解,再利用它们作线性组合使满足初始条件(2.3)。 这种思想方法,还可以从物理模型得到启示。从物理学知道乐器发出的声音可以分解成各种不同频率的单音,每种单音,振动时形成正弦曲线,其振幅依赖于时间t ,即每个单音可以表示成
(,)()sin u x t A t x ω= 的形式,这种形式的特点是:u (x ,t )中的变量x 与t 被分离出来。 根据上面的分析,现在我们就试求方程(2.1)的分离变量形式 (,)()()u x t X x T t = 的非零解,并要求它满足齐次边界条件(2.2),式中X (x ),T (t )分别表示仅与x 有关及仅与t 有关的待定函数。 由(,)()()u x t X x T t =得 2222()(),()()u u X x T t X x T t x t ??''''==?? 代入方程(2.1)得 2()()()()X x T t a X x T t ''''= 或 2()()()() X x T t X x a T t ''''= 这个式子左端仅是x 的函数,右端仅是t 的函数,只有它们均为常数时才能相等。令此常数为-λ,则有 2()()()() X x T t X x a T t λ''''==- 这样我们得到两个常微分方程: 2()()0T t a T t λ''+= (2.4) ()()0X x X x λ''+= (2.5) 再利用边界条件(2.2),由于u (x ,t )=X (x ) T (t ),故有 (0)()0,()()0X T t X l T t == 但T (t )不恒等于零,因为如果T (t )≡0,则u (x ,t )=0,这种解称
第三章 贝塞尔函数 对两个自变量的情形,在第二章中比较系统地介绍了分离变量法的基本思想 以及求解偏微分方程定解问题的主要步骤. 本章讨论多于两个自变量的情形,其求解过程和两个自变量情形基本相同,区别仅在于特征值问题的求解要用到一类特殊函数—贝塞尔(Bessel )函数. 本章前两节围绕一类特征值问题的求解,比较系统地介绍二阶常微分方程的幂级数解法,以及Bessel 函数的一些基本性质. 第三节介绍多于两个自变量情形的分离变量法. §3?1 二阶线性常微分方程的幂级数解法 3.1.1 常系数线性方程的基解组 在高等数学中,同学们已学过常微分方程的一些求解方法. 对于常系数线性常微分方程,只要求出特征方程的根,就很容易写出齐次方程的基解组,由此可得齐次方程通解表达式. 例1.1 求解下列齐次微分方程 (1) '''320y y y -+=. (2) '''4130y y y ++=. (3) '''440y y y ++=. 解 (1) 特征方程为 2320λλ-+=, 特征根为121,2,λλ== 故基解组为 2{, }x x e e . (2)特征方程为 24130λλ++=, 特征根为1223, 23i i λλ=-+=--,是一对共轭复数,基解组为(23)(23){, }i x i x e e -+--, 这两个解为复值函数. 为得到实值函数的基解组,利用齐次微分方程解的线性性质得 2(23)(23)1 cos3 (+ )2x i x i x e x e e --+--=, 2(23)(23)1 sin 3 ( )2x i x i x e x e e i --+--=-, 这两个实值函数22cos3, sin3x x e x e x --也是方程(2)的解,由此得方程(2)的基解组为 22{cos3, sin3}x x e x e x --. (3)特征方程为 2440λλ++=,
分离变量法例题 例:两块半无限大接地平行于xz 平面的导体板,一块位于y = 0,另一块位于y = d ;平行板的有限端x = 0处被与之绝缘并保持常电势φ0的导体板封闭,如图所示。求导体板间的电势。 解:对于本问题,求解区域是x > 0的两平行板之间,区域内无电荷分布,因此电势满足拉普拉斯方程。区域的边界在y = 0、y = d 、x = 0、及x → ∞处。本问题实际是一个二维问题,即静电势与z 无关。因此,本定解问题: 20??= ( x > 0,0 < y < d ) (1) 0x ?→∞= (2) 00y ? == (x > 0) (3) 0y d ?== (x > 0) (4) 00x ??== (0 < y < d ) (5) (2)的条件是我们通常的选择。实际上(2)、(3)、(4)、(5)为边界条件。 因本问题为二维问题,(),x y ??=。(1)在直角坐标系中可写成: 2222 0x y ????+=?? 分离变量法的核心是将多维函数分解为多个一维函数的乘积。令 (,)()()x y X x Y y ?= 将其带入上式得: 2222d d 0d d X Y Y X x y += 将x 变量项和y 变量项整理为: 22221d 1d d d X Y X Y x y =- 上式坐标仅是x 的函数,而右边仅是y 的函数。这样,我们就将变量分离了。在上面的方程中,两半无限大接地平行导体板 y = d
对任意x 和y 成立,方程两边必等于常数。即: 222221d 1d d d X Y k X Y x y ==- (3-3-3) 式中k 为实数常数,称为分离常数。为什么我们将常数写为k 2而不是-k 2,后面我们将清楚这一点。上式可分为两个微分方程: 2221d d X k X x = 2221d d Y k Y y =- 我们知道上面的微分方程k 为非零时的解为: kx kx X Ae Be -=+ sin()cos()Y C ky D ky =+ 若k = 0,根据边界条件只能得出零解,因此,k 为非零值。式中A 、B 、C 、D 为积分常数,由边界条件确定。这样,我们得到: [][sin()cos()]kx kx Ae Be C ky D ky ?-=++ 由边界条件(2),我们得到A = 0及k > 0。这就是为什么我们将常数写为k 2而不是-k 2的原因。它可使电势φ在x 方向单调地增加或单调地减少而不是振荡。由边界条件(3),我们得到D = 0。而边界条件(4)给出: sin()0kd = 由此式及条件k > 0,我们得到: ,1,2,3,n k n d π==? 我们不取n = 0的原因是因为它给出的是零解。因此,我们得到对应n 值的电势解: (,)sin(),1,2,3,n x d n n n x y B e y n d ππ?-==? 其中C 已并入B n 。因拉普拉斯方程是线性方程,任何解的线性叠加也是方程的解。因此,我们将所有n 值的解叠加起来得到了更为一般的解: 1(,)sin()n x d n n n x y B e y d ππ?∞-==∑ 式中B n 为常数。此解满足边界条件(2)、(3)、(4)。由边界条件(5), 我们有 01sin( )n n n B y d π?∞==∑ 上式是一在[0,d ]区间展开的正弦傅里叶级数,其系数B n 为: