狄利克雷问题(10.30)

  • 格式:docx
  • 大小:120.23 KB
  • 文档页数:5

狄利克雷问题
狄利克雷问题(Dirichlet problem )就是在给定边界条件的区域D 内求解拉普拉斯方程的问题,即
20, ()()()u u x f x x D ∇==∈∂
当区域D 是一个长为L ,宽为l 的矩形时,即
(){},:0, 0D x y x l y L =≤≤≤≤
此时狄利克雷问题就变成了求解带有若干边界条件的二阶偏微分方程,即
0xx yy u u +=
1212(,0)(), (,)(), (0,)(), (,)()u x f x u x L f x u y g y u l y g y ====
根据叠加原理,分别求解当120g g ==和120f f ==时方程的解,其解的加和即为原方程的解。

同时120g g ==与120f f ==这两种情形是等价的,因此只需要求解其中一个问题,并改变相应的变量就可以得到另一个解。

因此我们这里求解120g g ==的情形,即
0xx yy u u +=
12(,0)(), (,)(), (0,)(,)0u x f x u x L f x u y u l y ====
根据变量分离法,首先假设函数(,)u x y 可以写作
(,)()()u x y X x Y y =
代入微分方程可得
()''()"()()X x Y t X x Y y =-
移项整理得
2''()''()0()()
X x Y y v X x Y y =-≡-< (可以证明,当上式中的比例常数为非负数时,方程只有0解,与题设不符) 于是有
2"()()X x v X x =-
2"()()Y y v Y y =
对于方程2"()()0X x v X x +=,其对应的特征根方程有两个共轭的复根vi ±,因此其通解为
012()(cos sin )
X x e C vx C vx =+
将边界条件(0,)(,)0u y u l y ==代入微分方程可得 (0)()0X X l ==
所以
10, sin 0C vl ==
进而有
, n vl n v l
ππ==
所以 2()sin n X x C x l
π= 当21C =时,就得到了特征值22(/)v n l π=对应的特征函数,即
()sin n X x x l
π= 对于方程2"()()0Y y v Y y -=,其对应的特征根方程有两个不同的实根v ±,因此其通解为
34()vy vy Y y C e C e -=+
由于双曲余弦函数
cosh 2
vy vy
e e vy -+= 双曲正弦函数
sinh 2
vy vy
e e vy --= 所以函数()Y y 可以表示为
()cosh sinh n n Y y vy vy αβ=+ 将n v l
π=代入上式得 ()cosh sinh n n n n Y y y y l l
ππαβ=+ 所以1(,0)()u x f x =
(,)sin cosh sinh n n n n n n u x y x y y l l l π
ππαβ⎛⎫=+ ⎪⎝
⎭ 由上式可知,对于每个不同的n ,都有一个(,)n u x y 与之对应。

而1(,0)()u x f x =,2(,)()u x L f x =,由于1()f x ,2()f x 都是给定的函数,所以此处得到的还不是方程的解。

构造函数(,)U x y 为一切(,)n u x y 的加和,即
11(,)(,)sin
cosh sinh n n n n n n n n U x y u x y x y y l l l πππαβ∞∞==⎛⎫==+ ⎪⎝
⎭∑∑ 由初始条件1(,0)()u x f x =,2(,)()u x L f x =可得
()111()sin cosh 0sinh 0sin n n n n n n n f x x a x l l ππαβ∞
∞===+=∑∑
其中
n n a α=
211()sin cosh sinh sin n n n n n n n L n L n f x x b x l l l l ππππαβ∞∞==⎛⎫=+= ⎪⎝
⎭∑∑ 其中
cosh
sinh n n n n L n L b l l
ππαβ=+ 可以解得 n n a α=
cosh csch coth sinh n n n n n n L b a n L n L l b a n L l l l πππβπ-=
=- 所以
1(,)sin
cosh csch coth sinh n n n n n n n L n L n U x y x a y b a y l l l l l πππππ∞=⎡⎤⎛⎫=+- ⎪⎢⎥⎝⎭⎣⎦
∑ 其中
102()sin d l n n a f x x x l l
π=⎰ 202()sin d l n n b f x x x l l
π=
⎰ 因为120f f ==时的情形为: 0xx yy u u +=
12(,0)(,)0, (0,)(), (,)()u x u x L u y g y u l y g y ====
由于此情形与120g g ==时的情形完全等价,所以此情形下的解为
1(,)sin cosh csch coth sinh n n n n n n n l n l n V x y y c x d c x L
L L L L πππππ∞=⎡⎤⎛⎫=+- ⎪⎢⎥⎝⎭⎣⎦∑
其中
102()sin d L n n c g y y y L L π=⎰ 202()sin d L n n d g y y y L L π=⎰
所以,方程最终的解为 (,)(,)(,)u x y U x y V x y =+。