北邮数理方程课件第三章的分离变量法

合集下载

第三节分离变量法

第三节分离变量法

θ

1)
P3(cos θ)
=
1 (5 cos3 2
θ

3 cos θ)
·········
罗德利格(Rodrigues)公式

Pn(cos θ)
=
1 2nn!
dn d(cos θ)n
cos2 θ − 1 n

§ 3.5 例一
§ 3.5 例一
【问题】 电容率为ε的介质球置于均匀外电场E0中,求电势。
d dθ
(sin
θ
dΘ dθ
)
+
1 Φ sin2
θ
d2Φ dφ2
=
−n(n +
1)
(2)
式(2)两端乘以sin2 θ可得:
sin θ Θ(θ)
d dθ
(sin
θ
dΘ dθ
)
+ n(n
+
1) sin2 θ
=

1 Φ
d2Φ dφ2
=
m2
关于Φ(φ)的方程及其解:
1 Φ
d2Φ dφ2
=
−m2
Φ(φ) = Cm sin(mϕ) + Dm cos(mϕ)
ϕ(r,
θ,
φ)
=
(anmrn
n,m
+
bnm
1 rn+1
)Pnm(cos
θ)
cos(mφ)
+
(cnmrn
n,m
+
dnm
1 rn+1
)Pnm(cos
θ)
sin(mφ)
§ 3.3 拉普拉斯方程的通解
★ 拉氏方程在球坐标中的通解为

第三章 分离变量法

第三章 分离变量法
特征方程 r 2 pr q 0
p 4q r1, 2 , 2 (1) 有两个不相等的实根 ( p 2 4q 0) r 1 , r2
特征根
2
p
两个线性无关的特解
得齐次方程的通解为
y1 e ,
r1 x
y2 e ,
r2 x
y C1e
r1 x
C2e ;
第三章 分离变量法

分离变量法是求解偏微分方程最基本和 常用的方法。

理论依据:线性方程的叠加原理和 Sturm-Liouville 理论。 基本思想:将偏微分方程的求解化为相 应的常微分方程的求解.

主要内容
预备知识-常微分方程 带有两个第一(二)类齐次边界条件的齐次方程的问题 带有第三类齐次边界条件的齐次方程的问题 极坐标系下的变量分离法 非齐次问题的处理
r2 x
3.0 预备知识-常微分方程
(2) 有两个相等的实根 特解为
r1 x
( p 2 4q 0)
r 1x
r1 x
y1 e , y2 xe
齐次方程的通解为 y (C1 C 2 x )e
;
(3) 有一对共轭复根
特征根为 特解为
r1 i , r2 i , x x y1 e cos x, y2 e sin x,
( p 2 4q 0)
齐次方程的通解为 y e x (C1 cos x C 2 sin x ).
3.0 预备知识-常微分方程
y py qy 0
特征根的情况
r pr q 0
通解的表达式
2
r2 实根r1 r2 复根r1, 2 i

分离变量法3

分离变量法3

0 xa X X 0, X (0) X (a) 0
0 xa X X 0, X (0) X (a) 0
2 0
X 2 X 0
X (0) A B 0 A B0
X Ae x Be x X (a) Aea Be a 0 X 0
2 a n C n Dn ( x) cos xdx 0 a a
Cn e
nb a
Dn e

nb a
2 a n ( x) cos xdx a 0 a
Cn
nb 2 a n a ( x ) e ( x ) cos xdx 0 a a
n n , n 1,2,3, a n X n Bn sin x y a n
e
a y
2
Y Y 0
n 2 2 Yn 2 Yn 0 a
n y a

e Yn Cn
n y a
e Dn


n y a
e Dn
n y a
e Yn Cn
e Dn

n y a
n n n n y y y y n n a a a a un X nYn Cn Bn e Dn e cos x Cn e Dn e cos x a a
n n y y n a a sin u un x Cn e Dn e a n 1 n 1



n n y y n a a sin u x Cn e Dn e a n 1
u ( x,0) ( x) C n Dn sin

《分离变量法》课件

《分离变量法》课件
《分离变量法》 ppt课件
目 录
• 分离变量法简介 • 分离变量法的步骤 • 分离变量法的应用实例 • 分离变量法的优缺点 • 分离变量法的改进方向 • 分离变量法的未来展望
01
分离变量法简介
定义与特点
定义
分离变量法是一种求解偏微分方 程的方法,通过将多变量问题转 化为多个单变量问题,从而简化 求解过程。
交叉学科应用
探索分离变量法在交叉学 科中的应用,如生物医学 工程、环境科学等。
与其他方法的结合使用
与数值方法的结合
将分离变量法与其他数值方法(如有限元法、有限差分法等)结 合使用,形成更有效的数值计算方法。
与机器学习算法的结合
将分离变量法与机器学习算法相结合,用于数据分析和模式识别等 领域。
与优化算法的结合
工程与技术领域
分离变量法在解决工程与技术领域中的偏微 分方程问题方面具有优势,如结构分析、电 磁波传播、信号处理等。随着工程与技术的 不断发展,分离变量法有望在解决实际问题 中发挥更大的作用。
THANK YOU
感谢观看
立。
近似解
分离变量法得到的解是近似解,而 不是精确解。因为这种方法忽略了 各个变量之间的相互作用和影响。
数值稳定性
分离变量法在数值计算中可能存在 数值稳定性问题,例如数值误差的 累积和传播可能导致计算结果失真 或误差较大。
05
分离变量法的改进方向
算法优化
01
02
03
算法效率提升
通过改进算法结构,减少 计算复杂度,提高分离变 量法的计算速度。
精度控制
优化算法中的数值计算方 法,提高结果的精度和稳 定性,减少误差。
自适应调整
根据不同问题的特性,自 适应地调整算法参数,提 高分离变量法的适用性和 可靠性。

《分离变量法》课件

《分离变量法》课件
法的计算效率。
06
总结与展望
总结
内容回顾
详细梳理了分离变量法的基本概 念、应用场景、实施步骤和注意 事项,帮助学习者全面理解这一
方法。
案例分析
通过具体的案例分析,展示了分离 变量法在解决实际问题中的应用, 加深学习者对方法的理解和掌握。
互动问答
鼓励学习者在课程结束前提出疑问 ,并对常见问题进行了解答,有助 于巩固学习效果。
展望
新应用领域
实践应用建议
探讨分离变量法在未来可能的应用领 域,如人工智能、大数据分析等,为 学习者提供新的思路和方向。
为学习者提供将分离变量法应用于实 际问题的建议和指导,帮助他们更好 地实现学以致用。
方法改进
介绍分离变量法的最新研究进展和可 能的改进方向,激发学习者进一步探 索和研究。
谢谢您的聆听
02
分离变量法的原理
原理概述
分离变量法是一种求解偏微分方程的 方法,通过将多个变量分离,将复杂 的偏微分方程简化为一系列简单的常 微分方程,从而求解。
该方法适用于具有多个变量的偏微分 方程,特别是当各变量之间相互独立 时。
数学模型建立
首先,需要建立偏微分方程,并确定变量 的个数。
然后,通过适当的变换,将偏微分方程转 化为全微分方程。
求解过程
通过分离变量法,可以将 $u(x, t) = X(x) T(t)$,从而将波动方程 转化为 $X''(x) = -lambda X(x)$ 和 $T''(t) = -omega^2 T(t)$, 其中 $lambda$ 和 $omega$ 是常数。
应用实例二:化学反应动力学模型
总结词
描述化学反应速率
THANKS

数学物理方程分离变量法研究生高校本科生PPT课件

数学物理方程分离变量法研究生高校本科生PPT课件

l n , (n 1, 2,.....)
从而得到了固有值为
n
n2 2
l2
,
(n 1, 2,.....)
相应的固有函数为
(11) 固有值问题
Xn (x)
B sin
n
l
x
,
(n 1, 2,.....)
X ''(x) X (x) 0 (6)
(12)
X
(0)
0,
X
(l)
0
(10)
第12页/共97页
ut |t0
x)
Cn sin
n1
n
l
x
n a
(x) Dn
n1
l
sin
n
l
x
第14页/共97页
§2.2.1 齐次方程定解问题的解法
u
ut
|t0 ( |t0
x)
Cn
n1
sin
n
l
x
n a
(x) Dn
n1
l
sin
n
l
x
这两式正好是(x)和 (x)关于 sin
n
l
x的正弦展开。
根据Fourier级数展开法则(见下页附录),便可得到
§2.2.1 齐次方程定解问题的解法
iii)求方程满足边界条件的特解。
设u(x,t) X (x)T (t) (4)
为了求出T (t),把(11)式代入(7)式,得 T ''(t) a2T(t) 0 (7)
T
''(t)
n2 2a2
l2
T
(t)
0
其通解为一对共轭复根,即
n

数理方程分离变量法

数理方程分离变量法

nx nat nat u( x, t ) ( An cos Bn sin ) sin l l l n 1
现在要求出叠加系数 An和 Bn 满足初始条件
u
ut
t 0
( x)
t 0
( x)
0 x l
n a n ut ( x, 0) Bn sin x ( x) l l n 1
2 2
得固有值:
n n 2 a
(n 1,2,.....)
n x 固有函数: X ( x) sin n a
而 于是有
(n 1,2,.....)
Yn ( y ) An e
un ( x, y) ( An e
n y a
n y a
Bne
n y a
n y a
Bn e
叠加得
n x )sin a
n x )sin a
u( x, y) (An e
n 1

n y a Bnen Nhomakorabea y a
为确定叠加系数,将 u ( x,

y)
代入非齐次边界条件
n x ( An Bn ) sin a u0 n 1 n n b b ( A e a B e a ) sin n x U n n a n 1
T X " 2 aT X
2 T a T 0
X " X 0
n ( (2n2l1) )2 , n 1,2,
X n sin (2n2l1) x,
固有 值 (特 征值) 问题
Tn (t ) Cne

(2 n 1)2 2 a 2 4l 2

分离变量法

分离变量法

0.5
1 2 1.5 t1
0.5
00
1 0.8
0.6 0.4 x 0.2
欲使关于V(x,t)的定解问题可分离变量,W(x)要满足:
4 W 8 0
9
9
W (0) 0,W ( ) 2
求解得: W (x) x2
原问题变为:
Vtt
4 9
Vxx
(0
x
,t
0)
V x0 0,V x 0
V t0 sin 3x,Vt t0 0
0.5
00
1 0.8
0.6 0.4 x 0.2
(2). 当 0 时
X Ax B
A B 0
(3).当 0 时
X (x) Acos x B sin x
A 0, B sin L 0
9
1
0.5 n 0
0.5
1 2 1.5 t1
0.5
00
1 0.8
0.6 0.4 x 0.2
sin L 0
n
n2 2
L2(n 1, 2,3 )X Nhomakorabean
(
x)
Bn
sin
n
L
x
(10)
注:对于参数λ的某些值,问题(8),(9)的非平 凡解存在,称这种λ值为固有值(本征值);同时 称相应的非平凡解X(x)为固有函数(本征函数); 求解固有值和固有函数的问题称为固有值问题 (本征值问题)。
分离变量的核心问题是固有值问题(本征值问题)!
(2n
1) at
2L
4、一般解为:
u
x,
t
n0
An
cos
(2n
1) 2L
at
Bn

数理方程分离变量法

数理方程分离变量法

2 解的性质
un (x,t)
(Cn
cos n a t
l
Dn
sin
n a t)sin
l
n
l
x
An
cos(nt
n )sin
n
l
x
其中:
An Cn2 Dn2
n
n a
l
n
arctan
Dn Cn
x=x0时:
un (x0,t)
An sin
n
l
x0 cos(nt
n )
t=t0时:
un (x,t0 )
An
cos10(2n11 1)t
第12页/共82页
sin
(2n 1)
10
x
例2求下列定解问题
2u
t
2
a2
2u x2
,
u(0, t )
0,
u(l , t ) x
0,
u(x, 0)
x2
2lx,
u(x, 0) t
0,
0 x l,t 0 t0 0 xl
解: u(x,t) X (x)T (t)
X (0) 0,
X (10) 0
n n2 2 /100 , n 1,2,3,
X n (x)
Bn
sin
n
10
x
T 104T 0 Tn 100 n2 2Tn 0 Tn Cn cos10nt Dn sin10nt
于是得到一系列分离变量形式的特解
un
X nTn
Bn
sin
n
10
x(Cn
(Cn cos10nt
X (0) 0 X (10) 0
X X

分离变量-PPT精选文档

分离变量-PPT精选文档
第1部分 直角坐标系中的 分离变量法
适用条件
• 奇次问题 • 只有一个边界条件是非奇次的问题 • 超过一个非奇次边界条件的问题
主要内容
• • • • • 分离变量法 直角坐标系中热传导方程的分离 有限大物体的一维奇次问题 半无限大物体的一维奇次问题 非奇次问题分解成简单问题
2.1 分离变量法
数学模型 物理模型
d Γ 2 a Γ 0 d
2 d X x 2 X x 0 2 x
时间变量函数T(τ)满足微分方程
d Γ 2 a Γ 0 d
求解
Γ e
a 2
物理意义:
Γ 0
L x
TF (x )
0
求解
T x , X x Γ x
空间变量函数X(β,x) 特征值问题 时间变量函数(τ)的解:
2 d X x 2 X x 0 2 x
0 x L
x 0
Γ e
a 2
k 1
X h 0 1X x
L 0
2
1 ' T x , e X , x X , x x 'dx ' m m F N m 1 m
2 a m

2.4 半无限大物体的一维奇次问题
物理模型
数学模型
h1 , 0
x TF
T x , T 1 T x , T x a
空间变量函数X(x)满足微分方程
2 d X x 2 X x 0 2 x
0 x L
x 0
X 0 x
k
X hX0 x

北邮数理方程课件 第三章 分离变量法.

北邮数理方程课件 第三章 分离变量法.

第三章 分离变量法3。

2 基础训练3.2.1 例题分析例1 解下列定解问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂-==∂∂=><<∂∂=∂∂====0,20,00,0020022222t t lx x t u lx x u x uu 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)中边界条件,得00A 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 ln πβ从而得到一系列固有值与固有函数2224)12(ln n πλ+= ),2,1,0(2)12(sin)( =+=n x ln B x X n n π与这些固有值相对应的方程(3)的通解为),2,1,0(2)12(sin 2)12(cos )( =+'++'=n tlan D t l a n C t T n nn ππ于是,所求定解问题的解可表示为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 D33202)12(322)12(sin )2(2ππ+-=+-=⎰n l xdxln 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),(0332πππ++⨯+-=∑∞=例2 演奏琵琶是把弦的某一点向旁边拨开一小段距离,然后放手任其自由振动。

数学物理方程课件:8.1分离变量法介绍

数学物理方程课件:8.1分离变量法介绍

X (0) 0 和 X '(l) 0.
C1 0 和
C2 cos l 0
(k
1)2 2 2 l2
(2k
1)2 4l 2
2
k 0,1,2,3
X
(x)
C2
sin
(2k
1)x 2l
C.
(2k 1)2 2a2
T ''
4l 2
T 0;
T (t) Acos (2k 1)at B sin (2k 1)at ,
nat l
Bn
sin
nat l
) cos
nx l
.
n0 n 1,2,3
D.
u(x,t) A0 B0t
n1
(
An
cos
na l
t
Bn
sin
na l
t
)
cos
nx l
.
由初始条件:
u t0 (x)
A0
n1
An
sin
nx l
(x),
A0
1 l
l 0
(
)d ,
An
2 l
l 0
( ) cos n l
X
(
x)
C2
sin
nx l
:本征值
:本征函数
C2是积分常数。
X ''X 0;
:本征值方程
C.
T ''
n2 2a2 l2
T
0;
T (t) Acos nat B sin nat ,
l
l
A、B 是积分常数。
un (x,t)
( An
cos
nat l

3 分离变量法

3 分离变量法
X ( x) A cos x B sin x ,
根据(3.1.6)得
X (0) A 0 ,
X (l ) B sin l 0 ,
因为要求非零解,所以 sin l 0 ,因而得 l n , 则
n n , (n 1, 2, ) l
(3.1.1a) (3.1.1b) (3.1.1c)
其中 ( x), ( x) 为已知函数. 为求解混合问题的解, 考虑物理中的简谐波, 简谐波 前行遇到固定边界反射回来为同频率的反向波,与原来的波叠加形成驻波,各点 的频率相同, 只是振幅不同. 驻波可表示为 u( x, t ) X ( x)T (t ) , 只含有自变量 x 的 函数和 t 的函数,即具有可分离变量的形式. 下面试求方程(3.1.1a)的非平凡解, 使它满足齐次边界条件(3.1.1b),设解为
X 0 ( x) A0 x B0 ,
3 分离变量法
分离变量法是求解常微分方程的一个重要解法 . 对于某些典型区域上的常 系数线性偏微分方程,分离变量法同样是一个非常重要的解法. 通过使用分离变 量法把偏微分方程转化为常微分方程,解出这些常微分方程的解,则得到所要求 的偏微分方程的解 . 本章重点介绍分离变量法及特征函数法的基本思想和具体 步骤. 这里我们所讨论的问题主要是有界区域上的弦振动问题、热传导问题以及 圆域上的拉普拉斯方程问题 . 本章所需要的数学工具是常系数线性常微分方程 理论和傅里叶(Fourier)级数理论相关的知识.
u( x, t ) X ( x)T (t )
将(3.1.2)代入方程(3.1.1a),两边同除以 a 2 XT ,得到,
T (t ) X ( x) 2 a T (t ) X ( x)

分离变量法

分离变量法

第二章 分离变量法一 齐次偏微分方程的分离变量法1 有界弦的自由振动(1) 考虑两端固定的弦振动方程的混合问题⎪⎪⎪⎩⎪⎪⎪⎨⎧====><<∂∂=∂∂==)(|),(|0),(),0(0,0,01022222x u x u t l u t u t l x x u a t u t t φϕ ① 这个定解的特点是:偏微分方程是齐次的,边界条件是齐次的。

求解这样的方程可用叠加原理。

类似于常微分方程通解的求法先求出其所有线性无关的特解,通过叠加求定解问题的解。

所谓),(t x u 具有分离变量的形式,即)()(),(t T x X t x u =把)()(),(t T x X t x u =带入方程①中,可得到常微分方程定解为:),(t x u =∑∞=1),(n n t x u =l x n l t an D l t an C n n n πππ∑∞=+1sin )sin cos (其中:⎰=l n dx l x n x l C 0sin )(2πϕ,⎰=l n dx lx n x an D 0sin )(2πφπ 2离变量法的解题步骤可以分成三步:(一) 首先将偏微分方程的定解问题通过分离变量转化为常微分方程的定解问题。

(二) 确定特征值与特征函数。

(三) 求出特征值和特征函数后,再解其它的常微分方程,将所得的解与同一特征值报骊应的特征函数相乘得到所有分离变量的特解。

3 有限长杆上的热传导设有一均匀细杆,长为l ,比热为c ,热传导系数为k ,杆的侧面是绝缘的,在杆的一端温度保持为0度,另一端杆的热量自由散发到周围温度是0的介质中,杆与介质的热交换系数为0k ,已知杆上的初温分布为)(x ϕ,求杆上温度的变化规律,也就是要考虑下列问题:0,0,22222><<∂∂=∂∂t l x xu a t u (2.18) 0),(,0),0(=+∂∂=t l hu xt l u t u ),( (2.19) )()0,(x x u ϕ= (2.20) 其中ρc k a =2,00>=k k h注意到此定解问题中方程和边界条件均是齐次的,因此仍用分离变量法来求解。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第三章 分离变量法3。

2 基础训练3.2.1 例题分析例1 解下列定解问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂-==∂∂=><<∂∂=∂∂====0,20,00,0020022222t t lx x t u lx x u x uu 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)中边界条件,得00A 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 ln πβ从而得到一系列固有值与固有函数2224)12(l n n πλ+=),2,1,0(2)12(sin)(Λ=+=n x ln B x X n n π与这些固有值相对应的方程(3)的通解为),2,1,0(2)12(sin 2)12(cos )(Λ=+'++'=n tlan D t l a n C t T n nn ππ于是,所求定解问题的解可表示为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 D33202)12(322)12(sin )2(2ππ+-=+-=⎰n l xdxln lx x l C l n 故所求的解为x ln t l a n n l t x u n 2)12(sin 2)12(cos)12(132),(0332πππ++⨯+-=∑∞=例2 演奏琵琶是把弦的某一点向旁边拨开一小段距离,然后放手任其自由振动。

设弦 长为l ,被拨开的点在弦长的1n (0n 为正整数)处,拨开距离为h ,试求解弦的振动,即求解定解问题⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧=⎪⎪⎩⎪⎪⎨⎧≤--≤≤===<<=-====0|)1()()10(|0||)0(000000002t t t l x x xx tt u l n n l l x l h n x l hxn u u u l x u a u解:将)()(),(t T x X t x u =代入原方程及边界条件得''20T a T μ+= (1)''0(0)()0X X X X l μ⎧+=⎨==⎩(2)解(2)第一式可得x D x C x X μμsin cos )(+=由(2)的第二式得222ln n πμ=,Λ,3,2,1,sin)(==n lxn D x X n n π将μ代入(1)并解得latn B l at n A t T n n n ππsincos)(+= 1(,)(cossin )sin n n n n at n at n x u x t A B l l lπππ∞==+∑ 由初始条件得⎪⎪⎩⎪⎪⎨⎧≤--≤≤=∑∞=)1()()10(cos 00001l n n l l x l h n x l hxn l at n A n n π 0sin cos 1=∑∞=lxn l at n B l a n n nπππ 所以0=n B00222010100sin)1(2)(sin [200n n n n h n dxn l l x l h dx l xn l hx n l A l n n n πππ-=--+=⎰⎰ 从而∑∞=-=1020220sin cos sin 1)1(2),(n l xn l at n n n n n h n t x u ππππ 例3 求解细杆的导热问题,杆长l ,两端保持零度,初始温度分布20/)(|t x l bx u t -==.解:该问题的定解问题为⎪⎪⎩⎪⎪⎨⎧==-=====0||)(|0202l x x t xxt u u l x l bx u u a u (1)令)()(),(t T x X t x u =, 代入(1)第一式可得,0)()(''=+x X x X λ (2)0)()(2'=+t T a t T λ (3)由(2)得()sin X x A B =+ (4) 由(1)第三式可得0)()0(=t T X ,0)()(=t T l X(0)()0X X l λ≠==因为所以由0)0(=X 得0=A ,由()0X l B ==,0≠B 得),3,2,1(222Λ==n l n n πλ, 于是有sin n n n x X B lπ=,2222()(1,2,3,)n n a t a tl n n n T t C e C e n πλ--===L ,因此222(,)sinn a t ln n n xu x t C elππ-=, ∑∞=-=1sin),(222n lta n n lxn eC t x u ππ 将2)(l x l bx -作Fourier 展开得 ∑∞==-12sin )(n n l xn B l x l bx π 其中),3,2,1(]cos 1[4sin )(23302Λ=-=-=⎰n n n bdxl xn t x l bx l B l n πππ于是),3,2,1(]cos 1[433Λ=-==n n n b B C n n ππ因此x lk e k blxn e n n b t x u k l ta k n l ta n πππππππ∑∑∞=-∞=--=-=1)12(3313312sin)12(8sin ]cos 1[4),(22222222例4 在矩形域 b y a x ≤≤≤≤0,0内求Laplace 方程022222=∂∂+∂∂=∇y ux u u (1) 的解,使其满足边界条件⎪⎩⎪⎨⎧========0000b y y y y a x x u u Ay u u (2)(3)解:令 )()(),(y Y x X y x u =,代入式(1),有0)()(=-''x X x X λ (4)0)()(=+''y Y y Y λ (5)又由边界条件(3)得0)()0(='='b Y Y (6)当0<λ时,式(5)的通解为yyeC e C y Y λλ---+=21)(由式(6)有 021=+-C C021=+----bbeC e C λλ由此得 021==C C ,即式(5)、(6)无非零解.当0=λ时,式(5)的通解为01)(A y A y Y +=1)(A y Y ='由 00)()0(1=⇒='='A b Y Y ,得 常数)()(00A y Y =.当0>λ时,式(5)的通解为y B y A y Y λλsin cos )(+= 0cos sin )(=+-='y B y A y Y λλλλ由 0)0(='Y 得0=B ,由 0)(='b Y 得0sin =b A λλ,得πλλn b ==或0,即 ),2,1(222Λ==n b n πλ.由此可见,本征值为 ),2,1,0(222Λ==n b n πλ本征函数为 ),2,1,0(cos)(Λ==n y bn A y Y n n π将λ的值代入式(4),解得x D C X 000+=),2,1()(Λ=+=-n eD eC x X bx n n bx n n n ππ故问题的一般解为y b n e D e C x D C y Y x X y Y x X y x u n b xn n b x n n n n n πππcos )()()()(),(100100∑∑∞=-∞=⎪⎪⎭⎫ ⎝⎛+++=+= (7)由边界条件 00==x u 得到()0cos10=++∑∞=byn D C C n n n π 一个无穷级数等于零,说明各项系数均为零,故),2,1(000Λ==+=n D C C n n (8)又由Ay u a x ==得Ay y b n e D e C a D n b an n b a n n =⎪⎪⎭⎫ ⎝⎛++∑∞=-πππcos 10 将Ay 展开成Fourier 余弦级数,并比较系数有221221200Abb b A Aydy b a D b ==⎪⎭⎫ ⎝⎛=⎰故 aAbD 20=)1(cos 2cos 2220-==+⎰-πππππn n Ab dy b y n Ay b eD eC b ba n n ba n n (9) 从式(8)和(9)中解得),2,1()1(cos ,)1(cos 2222Λ=--=-=n ban shn n Ab D b a n sh n n Ab C n n ππππππ代入式(7)并整理得by n b x n sh ba n sh n n Ab x a Ab y x u n πππππcos 1cos 22),(122∑∞=-+= (10) 例5 带电云与大地之间的静电场近似匀强静电场,其电场强度0E 是垂直的.水平架设的输电线处在这个静电场中.输电线是导体圆柱.柱面由于静电感应出现感应电荷,圆柱附近的静电场也就不再是匀强的了.不过,离圆柱“无限远”处的静电场仍保持匀强,现研究导体圆柱怎样改变了匀强静电场(即讨论导线附近的电场分布).解:化成定解问题,取柱轴为z 轴,设导线“无限长”,那么场强和电势都与z 无关,只需在x ,y 平面上讨论.如图3-2所示,圆柱在x ,y 平面的截面是圆周为半径)a a y x (222=+作为静电场的边界,所以我们采用极坐标.柱外空间无电荷,电势满足二维Laplace 方程02222=∂∂+∂∂y ux u ,化成极坐标为 )(01122222a u u u >=∂∂+∂∂+∂∂ρϕρρρρ (1)边界条件:导体中的电荷不再移动,说明导体中电势相同,又因为电势具有相对意义,可以把导体的电势当作零,故0==a u ρ (2)“无穷远”处也为一个边界(圆内则考虑圆心点),“无穷远”处静电场仍为匀强静电场0E ,由于选取了x 轴平行0E ,故有00,y x E E E ==即 000,cos uE u E x E xρϕ∂-==-=-∂ 因此有0cos u E ρρϕ→∞→- (3)图3-2 输电线对带电云和大地之间电场的影响分离变量,令)()(),(ϕρϕρΦ==R u u代入方程(1),得0=Φ+Φ''λ (4)0222=-+R d dR d R d λρρρρ (5)因为极角具有周期性,)2,(),(πϕρϕρ+和应表示一个点,同一处的u 应该相同,故有),()2,(ϕρπϕρu u =+即 )()()2()(ϕρπϕρΦ=+ΦR R所以有 )()2(ϕπϕΦ=+Φ (6) 方程(6)称为自然周期条件.方程(4),(6)构成本征值问题,解之ϕϕϕm B m A sin cos )(+=Φ即),2,1,0(2Λ==m m λ方程(5)可以写成02222=-+R m d dR d R d ρρρρ (7) 为欧拉方程.作变换 te =ρ化成常系数线性微分方程,其通解为)0(1)0(ln 0>+==+=m DC R mD C R mmρρρ于是得到极坐标系中Laplace 方程的本征解ρϕρln ),(000D C u +=⎪⎪⎭⎫ ⎝⎛++=m m mmm m m D C m B m A u ρρϕϕϕρ1)sin cos (),( 一般解应叠加∑∞=⎪⎪⎭⎫ ⎝⎛++++=1001)sin cos (ln ),(m m m mmm m D C m B m A D C u ρρϕϕρϕρ (8) 由边界条件(2),有01)sin cos (ln 100=⎪⎭⎫ ⎝⎛++++∑∞=m m m m m m m a D a C m B m A a D C ϕϕ一个Fourier 级数为零,各系数为零,即01,0ln 00=+=+mmm m a D a C a D C 由此m m m a C D a D C 200,ln -=-=于是将解化简为∑∞=⎪⎪⎭⎫ ⎝⎛-++=1201)sin cos (ln),(m m m m m m a m B m A a D u ρρϕϕρϕρ (9) 再由边界条件(3),对于∞→ρ略去aρln及mρ1项,即01lim (cos sin )cos m m m m A m B m E ρρϕϕρϕ∞→∞=+=-∑比较系数100(1)0m m A E A m B =-=≠=代入方程(9),导体周围的电势分布2000(,)lncos cos a u D E E aρρϕρϕϕρ=-+ (10)例6 长为l 的理想传输线,一端接于电动势为 t v ωsin 0的交流电源,另一端开路,求解线上的稳恒电振荡.解:经历交流电的许多周期后,初始条件所引起的自由振荡衰减到可以认为已经消失,这时的电振荡完全是由交流电源引起的,所以叫稳恒振荡.因此是没有初始条件的问题:⎪⎩⎪⎨⎧==⎪⎭⎫ ⎝⎛==-==0,10002l x t j x xx tt u e v u LC a u a u ω 为了计算方便,将电动势 t v ωsin 0写成 tj e v ω0,最后将得到的解取虚部.由于振荡完全由交流电源引起,当然可以认为振荡的周期与交流电源相同,即令t j e x X t x u ω)(),(=代入方程得0)()()(22=''-t j t j e x X a x X e j ωωω即0)()()(2=⎪⎭⎫⎝⎛+''x X a x X x X ω其通解为x aj x aj BeAex X ωω-+=)(故有t j x aj x aj e BeAet x u ωωω][),(-+=由 t j x e v u ω00== 得0v B A =+ (1) 及 0==l x u 得 0=+-l aj l aj BeAe ωω(2)从式(1),(2)中解出⎪⎭⎫ ⎝⎛-=-=⎪⎭⎫ ⎝⎛=-=--l a e jv ev B l a e jv ev A laj laj l aj l aj ωωωωωωsin 21sin 21020020带入解的表达式,得t j t j x l aj x l aj tj l x aj l x aj e a l x l a v l a j e e ev l a j e e ev t x u ωωωωωωωωωωω⎪⎭⎫ ⎝⎛-=⎪⎭⎫ ⎝⎛-=⎪⎭⎫⎝⎛--=------sin )(sinsin 2][sin 2][),(0)()(0)()(0取虚部,并以 LCa 1=代入,得传输线内稳恒的电振荡 t LCl x l LC v t x u ωωωsin sin )(sin ),(0-=例7 试解出具有放射衰变的热传导方程2220x u u a Ae x tα-∂∂-+=∂∂ 已知边界条件为0|0,|0x x l u u ====初始条件为00|t u T ==(常数) 解 令2221,Ab B a a ==,定解问题可以化为 222000||0|x x x l t u u b Be t x u u u T α-===⎧∂∂=+⎪∂∂⎪⎪==⎨⎪=⎪⎪⎩由于对应的齐次问题具有第一类边界条件,故令11(,)()sin ,()sin xn n n n n n u x t T t x Be f t x l l αππ∞∞-====∑∑ 代入上述方程和初始条件得21101()()sin ()sin (0)sin n nn n n nn n b n n T t T t x f t x l l l n T x T l ππππ∞∞==∞=⎧⎡⎤⎛⎫'+=⎪⎢⎥ ⎪⎪⎝⎭⎢⎥⎣⎦⎨⎪=⎪⎩∑∑∑ 即2000()()()22(0)sin [1(1)]n n n l n n n b T t T t f t l T n T T d l l n ππξξπ⎧⎛⎫'+=⎪ ⎪⎪⎝⎭⎨⎪==--⎪⎩⎰ (1)(2) 其中22021(1)22()sin n ll n e n n B f t Be d l l l n l ααξππξξπα--⎡⎤--⎣⎦==⎛⎫+ ⎪⎝⎭⎰(3) 求解(1)式得到22222()()n b t l n n n l T t C e f t n bππ⎛⎫- ⎪⎝⎭=+ (4)将(2)(3)、式代入(4)式得02221(1)221(1)n lnn e T B C n b n n l απππα-⎡⎤--⎣⎦⎡⎤=---⎣⎦⎛⎫+ ⎪⎝⎭2202221(1)22()1(1)1n l n b n b t t nl l n e T B T t e e n b n n l απππππα-⎛⎫⎛⎫-- ⎪ ⎪⎝⎭⎝⎭⎡⎤⎡⎤--⎣⎦⎢⎥⎡⎤=----⎣⎦⎢⎥⎛⎫⎣⎦+ ⎪⎝⎭故得原定解问题的解为22022121(1)22(,)1(1)1sin n l n b n b t t n l l n e T B n u x t e e x n b n l n l αππππππα-⎛⎫⎛⎫∞-- ⎪ ⎪⎝⎭⎝⎭=⎧⎫⎪⎪⎡⎤⎡⎤--⎪⎪⎣⎦⎢⎥⎡⎤=----⎨⎬⎣⎦⎢⎥⎛⎫⎪⎪⎣⎦+ ⎪⎪⎪⎝⎭⎩⎭∑ 即2202121(1)22(,)1(1)1sin n l n n t t n al al n e T A n u x t e e x n n l n l αππππππα-⎛⎫⎛⎫∞-- ⎪ ⎪⎝⎭⎝⎭=⎧⎫⎪⎪⎡⎤⎡⎤--⎪⎪⎣⎦⎢⎥⎡⎤=----⎨⎬⎣⎦⎢⎥⎛⎫⎪⎪⎣⎦+ ⎪⎪⎪⎝⎭⎩⎭∑ 例8 在环形域)0(22b a b y x a <<≤+≤内求解下列定解问题.⎪⎪⎩⎪⎪⎨⎧=∂∂=<+<-=∂∂+∂∂=+=+0,0),(12222222222222b y x a y x n u u b y x a y x y u x u 解:由于求解区域式环形区域,所以我们选用平面极坐标系,利用直角坐标系与极坐标系之间的关系⎩⎨⎧==θρθρsin cos y x 可将上述定解问题用极坐标θρ,表示:⎪⎪⎩⎪⎪⎨⎧=∂∂=<<=∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂==0,0,2cos 12112222b a u u b a uu ρρρρθρθρρρρρ (1)(2) 这是一个非齐次方程附有齐次边界条件的定解问题.采用固有函数法,并注意到圆域内Laplace 方程所对应的固有函数,可令问题(1)--(2)的解的形式为]sin )(cos )([),(0∑∞=+=n n n n B n A u θρθρθρ代入式(1)并整理得到θρθρρρρρθρρρρρ2cos 12sin )()(1)(cos )()(1)(202222=⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡-'+''+⎥⎦⎤⎢⎣⎡-'+''∑∞=n n n n n n n n B n B B n A n A A比较两端关于θθn n sin ,cos 的系数,可得2222212)(4)(1)(ρρρρρρ=-'+''A A A (3))2(0)()(1)(22≠=-'+''n A n A A n n nρρρρρ (4) 0)()(1)(22=-'+''ρρρρρn nnB n B B (5)再由条件(2)得0)()(0)()(='=='=b B a B b A a A nn nn (6)(7)方程(4)与(5)都是齐次的欧拉方程,它们的通解分别为nn nnn n n n n n d c B d c A --'+'=+=ρρρρρρ)()(其中n nn n d c d c '',,,都是任意常数.由条件(6)与(7)可得 0)()2(0)(≡≠≡ρρn n B n A下面的任务就是要确定2()A ρ.方程(3)是一个非齐次的欧拉方程,利用待定系数法可求得它的一个特解42)(ρρ=*A所以,它的通解为422212)(ρρρρ++=-C C A由条件(6)确定21,C C ,得442244244661)2(2b a b a b a C b a b a C +--=++-= 因此42442244244662)2(2)(ρρρρ++--++-=-ba b a b a b a b a A 原定解问题的解为θρρρθρ2cos ])()2()2[(1),(4442224426644b a b a b a b a ba u +--+++-=- 例9 求解一端固定,一端作周期运动t ωsin 的弦的振动问题.⎪⎪⎩⎪⎪⎨⎧≤≤====><<=====lx u u t u u t l x u a u t t t l x x xx tt 00,0sin ,00,00002ω解法一:令 ),(),(),(t x W t x v t x u += 取 t lxt x lt t t x W ωμμμsin )()()(),(112=+-= (1) 则定解问题转化为⎪⎪⎩⎪⎪⎨⎧-=====-xl x v x v t l v t v t x l v a v t xx tt ωωω)0,(,0)0,(0),(),0(sin 22(2)又令 ),(),(),(t x v t x v t x v I II+=其中 ),(),,(t x v t x v I II分别满足⎪⎪⎩⎪⎪⎨⎧-=====I I I I IIxl x v x v t l v t v v a v t xx tt ω)0,(,0)0,(0),(),0(2 (3)⎪⎪⎪⎩⎪⎪⎪⎨⎧====+=I I I I I I I I I II I 0)0,()0,(0),(),0(sin 22x v x v t l v t v t x l v a v t xx tt ωω (4)由分离变量法求解式(3),得 x ln l at n an l t x v n nπππωsin sin 2)1(),(122∑∞=I-=(5)用固有函数法求解问题(4),即设lx n t f t x ll x n t v t x v n n n n πωωπsin)(sin sin)(),(121∑∑∞=∞=I I==其中t n dx l x n t x l l t f n l n ωπωπωωsin 2)1(sin sin 2)(2102+-==⎰ 代入方程(4),得⎪⎩⎪⎨⎧=='==+''),2,1(0)0(,0)0()()()(222Λn v v t f t v l n a t v n nn n nπ 由参数变易法,得ττπωττπωτπωπττπωτπωπd l t a n l t a n n a n l d lt a n n a n lt v t n tn n ⎰⎰⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡⎪⎭⎫⎝⎛---⎪⎭⎫ ⎝⎛-+--=--=++021021)(cos )(cos 212)1()(sin sin 2)1()( 因此x l n t t t t a n l t x v n n n n n n πωωωωωωωωπωsin sin sin sin sin )()1(),(1221⎥⎦⎤⎢⎣⎡---++-=∑∞=+I I(6)其中 lan n πω= 原定解问题的解为t lxt x v t x v t x u ωsin ),(),(),(++=I I I 解法二:取 t al xa t x W ωωωsin sin sin),(=则原问题化为⎪⎪⎩⎪⎪⎨⎧-====''=''====ωωωa l x a v v v v v a v t t t l x x xx tt sin sin ,00,00002 (7)注意到方程和边界条件同时齐次化了. 用分离变量法解方程(7),得l xn l at n l l n al a t x v n n n nn ππββααωπωsin sin sin 1sin 11sin ),(1⎪⎪⎭⎫ ⎝⎛-=∑∞= 其中 ⎪⎭⎫⎝⎛-=⎪⎭⎫ ⎝⎛+=l n al n an n πωβπωα, 原定解问题的解为t al xa t x v t x u ωωωsin sin sin),(),(+= 应当指出,同两种方法得到的定解问题的解在形式上不一样,但可以证明它们是等价的,这是由定解问题解的唯一性决定的.例10 求下列定解问题⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂===><<+∂∂=∂∂====0,00,000022222t t l x x t u u B u u t l x A x u a tu (1)(2)(3)的解,其中A ,B 均为常数.解:这个定解问题的特点是:方程及边界条件都是非齐次的.根据上述原则,首先应将边界条件化成齐次的.由于问题中方程(1)的自由项及边界条件(2)都与t 无关,所以我们有可能通过一次代换将方程及边界条件都变成齐次的.具体做法如下:令 )(),(),(x W t x V t x u += 代入方程(1),得A x W x V a t V +⎥⎦⎤⎢⎣⎡''+∂∂=∂∂)(22222为了使这个方程及边界条件同时化成齐次的,选)(x W 满足⎪⎩⎪⎨⎧===+''==BW W A x W a l x x ,00)(02 (4)问题 (4)是一个二阶常系数线性非齐次常微分方程的边值问题,它的解可以通过两次积分求得:x l B a Al x a A x W ⎪⎭⎫⎝⎛++-=22222)( 求出函数)(x W 之后,再由问题(1)-(3)可知函数),(t x V 为下列定解问题⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂-===><<∂∂=∂∂====0),(00,000022222t t l x x t V x W V V V t l x x V a tV (5)(6)(7)的解.采用分离变量法,可得式(5)满足齐次边界条件(6)的解为 x l n t l a n D t l a n C t x V n nn πππsin sin cos ),(1∑∞=⎪⎭⎫ ⎝⎛+=(7) 利用式(7)中第二个条件可得0=n D . 于是定解问题(5)-(7)的解可表示为∑∞==1sin cos),(n n x ln t l a n C t x V ππ 代入式(7)中第一个条件,得∑∞==-1sin)(n n x ln C x W π 即∑∞==⎪⎭⎫ ⎝⎛+-1222sin 22n n x l n C x l B a Al x a A π 由傅氏级数的系数公式可得πππππππn B n a Al n n a Al xdx l n x l B aA xdx l n x l a A xdxl n x l B a Al x a A l C ll l n cos 22sin 2sin sin 222222233220220220222⎪⎪⎭⎫ ⎝⎛++-=⎪⎭⎫ ⎝⎛+-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+-=⎰⎰⎰ (8) 因此,原定解问题的解为∑∞=+⎪⎭⎫ ⎝⎛++-=1222sin cos 22),(n n x l n t l a n C x l B a Al x a A t x u ππ其中n C 由式(8)确定.例11在扇形区域内求下列定解问题200||0|()au u u u f θθαρθ===⎧∇=⎪==⎨⎪=⎩的解。

相关文档
最新文档