§2.2 初边值问题的分离变量法
数学物理方法分离变量法

1 0.5
-0.5
2.5
5
7.5
10 12.5
15
-1
波腹 波节
15
6、 分离变量法概要:
(1)将偏微分方程化简为常微分方程(U=XT) (2)确定固有值和固有函数(利用边界条件) (3)确定形式解(级数形式解) (4)确定级数解中的待定常数(利用初始条件)
16
例:求解
utt a2uxx 0
R(0)
欧拉方程
30
第二步:求解周期本征值问题和欧拉方程
0 ( ) ( 2 )
nn
n2
( )
an
cos
n
bn
sபைடு நூலகம்n
n
n 0,1, 2,
2R R R 0
X(0) 0
X(l) 0
C1 0 C2 sin l 0
非零解 C2 0
sin l 0
n2 2
l2
n 1,2,3
则X(x)的一族非零解为
n x
X ( x) C2 sin l
C2是积分常数
上解称为满足边界条件的固有解(特征解),λ称为
固有值(特征值),sin函数称为固有函数(特征函数)。
(3)当r1、2 i时,y(x) ex (c1 cos x c 2 sin x)
7
本方程特征方程r2+λ=0,由上面结论知,方程的解与 λ的不同取值有关,分情况讨论:
(1) 0
X ( x) C1e x C2e x
X(0) 0
C1 C2 0
x, t 是相互独立的变量
T '' a2T 0
分离变量法

∞
kπ a kπ a ⎞ kπ ⎛ t + Bk sin t ⎟ sin x = ∑ ⎜ Ak cos l l ⎠ l k =1 ⎝
∞
使得
9
k π a Bk kπ ψ ( x ) = ut ( x , 0) = ∑ sin x l l k =1
∞
kπ ϕ ( x ) = u ( x , 0) = ∑ Ak sin x l k =1
X ( x) = C .
情形(C)
λ >0
其通解为
X ( x) = C1 cos λ x + C2 sin λ x,
由边界条件推出 C2 = 0,
22
再由 X ′(l ) = C1 λ sin λ l = 0 知道为了使 必须
C1 ≠ 0,
sin λ l = 0.
λ l = kπ ,
kπ λ = λk = 2 , l
sin λ l = 0.
λ l = kπ ,
kπ λ = λk = 2 , l
2 2
于是有
(k = 1,2,3," ).
本征值
(k = 1,2 ,3," ).
X k ( x) = Ck sin
kπ x, l
(k = 1,2," )
本征 函数
7
k 2π 2 λ = λk = 2 , l
(k = 1,2 ,3," ) 代入另一个方程可得
u ( x, t ) = X ( x)T (t )
把变量形式的解代入方程可得
XT ′′ = a 2 X ′′ T
即
T ′′(t ) X ′′( x) = 2 a T (t ) X ( x)
以及
第二章分离变量法

第⼆章分离变量法第⼆章分离变量法§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 uu 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 的函数,只有它们均为常数时才能相等。
分离变量法的解题步骤总结

n
如果电势不依赖于方位角,则
B n () r [ A r ( ] P ( c o s ) n n 1 ) r n 0 S
n nS
几种常见的边界条件
导体为等位体; ˆ 均匀场: E i E , E zE r c o s z 0 0 0 S
B m n m (r) [A r (n ] P c o s )c o sm n ( 1 ) r n 0m 0 S
n n mS n
D m n m [ C r (n ]P c o s )s inm n ( 1 ) r n 0m 0 S
n n mS
r S
柱的轴心或球心处若没有线电荷或点电荷,则 电位为有限值; 若电荷分布于有限区域,则无穷远处电位趋近 于零; 周期边界条件。
分离变量法的解题步骤总 结
解题步骤
确定求解区域,写出电势所满足的方程(一般 为Laplace方程)和边界条件(包括物理边界 条件和自然边界条件); 根据边界的形状选取坐标系; 写出通解(其中包含有待定系数); 把边界条件代入通解中,确定待定系数,从而 得到问题的解; 对问题进行讨论。
二维情况下直角坐标系通解形式
xy , A x BC D y A s i n k x B c o s k x C s i n h ky D c o s h ky 1 x 1 x 1 y 1 y A s i n h k x B c o s h k x C s i n ky D c o s ky 2 x 2 x 2 y 2 y
分离变量法

1.2
分离变量法的物理意义
令 Nn =
2 A2 n + Bn ,
αn = arctan 混合问题 (1) 的解的每一项可化为 un (x, t) = Nn sin
Bn , An
nπ anπ x sin t + αn . l l
un (x, t) 是振动元素。对于弦上任意一点 x , un (x, t) 描述了这一 anπ nπ , 频 率 ωn (x) = ,初 点 的 简 谐 振 动 , 其 振 幅 an (x) = Nn sin l l l n−1 相位为 αn 。于是 ,当 x = 0, , . . . , l, l 时,振幅 an (x) = 0 ;当 n n l 3l 2n − 1 x= , ,..., l 时,振幅 an (x) = ±Nn 达到最大。因此弦的振动 2n 2n 2n 可以看作一系列具有特定频率的驻波的叠加。 特别地,考虑定解问题 utt − a2 uxx = A (x) sin ωt, (x, t) ∈ (0, l) × (0, +∞) , u (x, 0) = ut (x, 0) = 0, x ∈ [0, l] , u (0, t) = u (l, t) = 0, t ∈ [0, +∞) . x 我们可得
4
将 u (x, t) , f (x, t) , ϕ (x) , ψ (x) 均按特征函数系展开:
∞
u (x, t) =
n=1 ∞
Tn (t) sin fn (t) sin
n=1 ∞
nπ x, l nπ x, l
f (x, t) = ϕ (x) =
n=1 ∞
ϕn sin ψn sin
n=1
nπ x, l nπ x, l
分离变量法

12.2.3直角坐标系分离变量例题分析
上面我们已经研究的例题12.2.1讨论的是两个边界点均
为第一类齐次边界条件的定解问题.下面讨论的例
题12.2.2是既有第一类,也有第二类齐次边界条件的定解 问题;而例题12.2.3讨论的是均为第二类齐次边界条件的 定解问题,注意到本征值和本征函数的区别.
(2) 若
,则(12.2.22)的解为
由(12.2.23)得 ,应该排出
只能得到无意义的解
(3)若
,则方程的解是 ,
由(12.2.23)则
注意到 .在
且要得到非零解,只有
条件下, 可以是任意常数.条件 ,即
故得到本征值为
相应的本征函数是
系数B可以在求通解时考虑进去,故此将系数认为是
归一化的 .
y y
C1er1x C1erx
C2er2x C2 xerx
y
e
x
(C1
cos
x
C2
sin
x)
求解(12.2.5),将 三种可能逐一加以分析
(1)
(12.2.5)的解为
和
由(12.2.7)确定,即有
由此解出
(2)、
被排除 方程(12.2.5)的解是
和
由(12.2.7)确定,即
第十二章 分离变量法
本章中心内容
用分离变量法求解各种有界问题;
本章基本要求
掌握有界弦的自由振动解及其物理意义
着重掌握分离变量法的解题思路、 解题步骤及其核心问题---本征值问题
第十二章 分离变量法
问题的引入
uutt
分离变量法在解初边值问题中的应用

分离变量法在解初边值问题中的应用
分离变量法是一种常用的解决初边值问题的方法,它可以将复杂的微分方程组分解为一组简单的微分方程,从而解决初边值问题。
分离变量法的基本思想是将复杂的微分方程组分解为一组简单的微分方程,每个微分方程
只包含一个未知函数的一个变量,而其他变量都是常数。
这样,每个微分方程都可以单独
解决,而不需要考虑其他变量的影响。
分离变量法在解决初边值问题时,首先要将复杂的微分方程组分解为一组简单的微分方程,然后将每个微分方程的解表示为一个未知函数的一个变量的函数,最后将这些函数组合起来,得到最终的解。
分离变量法在解决初边值问题时,可以有效地减少计算量,提高计算效率,并且可以得到
更加精确的解。
因此,分离变量法在解决初边值问题中有着重要的应用。
分离变量法-2

T (t) a2T (t) 0
X ( x) X ( x) 0
三、在右列方程组中,解出非零的 X ( x) 。
以下的任务:
X ( x) X ( x) 0
X (0) X (l) 0
确定λ 取何值时 ,方程 X ( x) X ( x) 0 有满足条件
l
t sin n
l
x
(n 1,2,3)
由上式所确定的 u( x, t) , 确实是
2u
t
2
a2
2u x2
,
0 xl,t0
u |x0 0 , u |xl 0 , t 0
u
|t
0
(
x
)
;
u (x) , t
t0
其中, T(t) 0 。盖由于 T(t) 0 ,则 u( x, t ) 0 , 所涉及的解,显然
不是我们所需要的(零解!!!)。
由此可见,只有 X (0) X (l) 0 。将此结果与所得到的常微分方
程
中的第二个方程(关于X )联立
X ( x) X ( x) 0 X (0) X (l) 0
(
t
)
C
n
cos
a
n
l
t
Dn
sin a n
l
t
,
(n 1,2,3)
将 Xn( x) 和 Tn (t ) 一并代入 u( x, t) X ( x)T(t) , 经整理后得
un
(
x,
t
)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
上页 下页 返回
由X ( l ) hX ( l ) 0得:
B( cos l h sin l ) 0
B 0 cos l h sin l 0
(2.9) (2.10) (2.11)
即,应满足: 令 v l
tan l h
则
v tan v lh
由边界条件(2.3)、(2.4),得
X (0) 0, X ( l ) hX ( l ) 0
下面讨论固有值问题(2.6)、(2.7):
(ⅰ)当 0 时,只有平凡解 X 0 ;
(ⅱ)当 0时,
X A cos x B sin x
(2.8)
由X (0) 0 得: A 0
X m X n )dx ( n m ) X m X n dx 0 0 ( X m X n 0
X m X n )dx (Xm Xn
0 l
l
l
Xm X n )dx ( X m Xn Xm Xn )dx 0 ( X m X n 0
T V V T a V T 2 V aT
2 2 T a T 0,
V V 0
x2 y2 l 2
Helmholts方程
代入边界条件, 得 V
11
0
上页 下页 返回
V V 0 得偏微分方程固有值问题: (I) V x2 y2 l 2 0
J n() 有无穷多个正根,从小到大分别记为:
n) (m ( m 1,2,)
15
上页 下页 返回
为使 R l
0, 只需取
( n) m
l
( n) m
( m 1,2,) , ,
2
( n) ( n) 相应的有, Rm (r ) J n ( m r ) l
2 2 R ( ) R ( ) ( n ) R( ) 0 2
可用幂级数解法, 得到原点处有界的解为
R( ) J n ( )
2
n
( 1) k 0 k !( n k )! 2
k
2k
(2.20) 0 ( m n) 0 sin m x sin n xdx M n 0 ( m n) (2.21)
l
7
上页 下页 返回
1 Ak Mk
0 ( x ) sin
0 () sin
l
l
k xdx
(2.22)
代入(2.18)得:
1 u( x , t ) k 1 M k k d e
¶ Ô È Î Ò à p 0¾ õ Ê Ç Ò » Ö Â Ê Õ Á ²
µ Ä ¡ £µ ± ( x ) Î ª Ó Ð ½ ç º ¯ Ê ý Ê ±£ ¬± í ´ ï Ê ½ (2.23)Ö Ð µ Ä ¼ ¶ Ê ý ¹ × Ó Ú
13
上页 下页 返回
3 当 0 时,方程的通解为 ○
A cos B sin ;
是周期为 2 只有当 n , 即 n2 (n 1, 2, )时,
的函数
n An cos n Bn sin n ;
综上所述, n n2 ( n 0,1, 2, )是固有值问题 (II)的 固有值,相应的固有函数是
其中h为正常数。 分离变量法求解: 令 u( x, t ) X ( x )T (t ) 代入(2.1),得: 即:
2
XT a 2 X T
T X 2 aT X
上页 下页 返回
T a 2T 0
(2.5) (2.6) (2.7)
X X 0
0
l
i .e ., sin m x sin n xdx 0 ( m n)
0
l
(2.20)
当m n k时 记M k X k dx
2 0
l
1 cos 2 k x M k sin k xdx dx 0 0 2 tan k l l sin 2 k l l l 2 4 k 2 2 k 1 tan 2 k l k tan k l h l h Mk 0 (2.21) 2 2 2( h k )
§2
初边值问题的 分离变量法
2015-3-6
1.一个空间变量的情形
ut a 2 u xx 0 ( t 0, 0 x l ), t 0 : u ( x ), x 0 : u 0, x l : u x hu 0, ( 2.1) ( 2.2) ( 2.3) ( 2.4)
(2.13)
可用图解法求解(2.13) :
4
上页 下页 返回
f
f tan v
求解( 2.13)得: v v k ( k 1,2) 且v k 满足
v2
O
v1
2
f
v lh
v
1 ( k ) v k k 2
从而(2.11)的解,即, (2.6)(2.7)的固有值为:
l
l
l
Xm Xn ) 0 (Xm Xn
X m ( l )( hX n ( l )) ( hX m ( l )) X n ( l ) 0
( n m ) X m X n dx 0
0
9
l
上页 下页 返回
当m n时 ( n m ) 0 所以 X m X n dx 0
0
2 r R rR ( r ) R 0 2
由自然周期条件,得
( ) ( 2 )
0 下面求解固有值问题 (II) ( ) ( 2 )
1 当 0 时,方程没有非零周期解; ○ 2 当 0 时,方程的周期解为 A ; ○ 0 0
18
上页 下页 返回
3.混合问题古典解的存在性
Ó É Ó Ú ± í ´ ï Ê ½ (2.23)Ö Ð º ¬ Ó Ð Ò ò × Ó e µ ±t Ê ±£ ¬¼ ¶ Ê ý
a 2 k t
£ ¬Ò ò ´ Ë ¶ Ô È Î Ò à 0£ ¬
k 1
p k
e
a 2 k t
k 1 a 2 k t
sin k x
(2.18)
是满足泛定方程(2.1)及边界条件(2.3)、(2.4)的解。 为使(2.18)满足初始条件(2.2),须成立:
( x ) Ak sin k x
k 1
(2.19)
由于固有函数系X k sin k x在[0, l ]上正交,即
n) 从而 (Ⅰ)的固有值为 (m ,相应的固有函数为
( n) U nm Rm ( r ) n () n) n) (m (m a nm J n ( r ) cos n bnm J n ( r ) sin n l l
将 代入T 的方程,得T
(n) m
( n) m
作极坐标变换 x Байду номын сангаас r cos , y r sin
1 1 Vrr r Vr r 2 V V 0, 0 r l ,0 2 V ( l , ) 0 (齐次边界条件) V (0, ) (自然边界条件) (自然周期条件) V ( r , ) V ( r , 2 )
a 2 k t
sin k x
(2.23)
固有函数系X k sin k x在[0, l ]上的正交性证明:
证明:
n X n 0 Xn m X m 0 Xm
(1) ( 2)
上页 下页 返回
8
用 X m , X n 分别乘以(1)、(2)式,相减后在[0, l ]上积分:
代入初始条件,得
n) n) (m (m a nm J n ( l r ) cos n bnm J n ( l r ) sin n ( r , ) n 0 m 1
(m j ) 0, l 2 r Jn( r) Jn( r ) dr l ( n) 0 l l J ( n 1 m ) 0, ( m j ) 2
l 2 l
10
上页 下页 返回
2.圆形区域上的热传导问题
u t a 2 ( uxx u yy ), x 2 y 2 l 2 , t 0 u( x , y ,0) ( x , y ) u 2 2 2 0 x y l (2.25)
解: 令 u( x , y, t ) T ( t )V ( x , y ), 代入方程, 得
图2.1
vk k l
2
( k 1,2,)
5
上页 下页 返回
相应的固有函数为:
vk X k Bk sin k x Bk sin x ( k 1,2,) l k时,T满足:T a 2 k T 0
(2.15)
解得:
n An cos n Bn sin n (n 0,1, 2, )
下面讨论R(r): 由自然边界条件和齐次边界条件,得
R(0) , R( l ) 0
14
上页 下页 返回
2 2 r R rR ( r )R 0 讨论固有值问题 (Ⅲ) R(0) , R( l ) 0 令 r , 并将 n2 代入方程得n阶Bessel方程
( n) m ( n) k