数学物理方程第二篇分离变量法

第二章 分离变量法

分离变量法是求解偏微分方程定解问题最常用的方法之一，它和积分变换

法一起统称为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 =， 或

, 1.i i i A T T i n λ=≤≤ （）

上式说明，正交矩阵T 的每一列都是实对称矩阵A 的特征向量，并且这n 个特征向量是相互正交的. 由于此结论在一定意义下具有普遍性，我们以定理的形式给出.

定理 设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 λ≤≤.

例 设n b R ∈，求解线性方程组 Ax b =.

解 由于向量组{1}i T i n ≤≤线性无关，故可做为n R 的一组基. 将,x b 按此组基分别展开为1

1 ,n

n

i i i i i i x x T b bT ====∑∑，则Ax b =等价于

1

1

n

n

i

i

i i

i i x AT bT ===∑∑，

或

1

1

n n

i i i

i i

i i x T bT λ===∑∑，

比较上式两边i T 的系数可得

1, 1i i i x b i n λ-=≤≤，

12( ... )n x x x x T =便是原问题的解.

例 设0n x R ∈,12()((),(),...,()), 0n n f t x t x t x t R t T =∈>. 求解非齐次常微分方程组

0(), (0)dx

Ax f t x x dt

=+=， （）

其中

'''

12((),(),...,()),0n dx x t x t x t t dt

T =>. 解 类似于上例，将0,,()x x f t 按基{1}i T i n ≤≤分别展开为

01

1

1

, , ()()n

n

n

i i i i

i i i i i x x T x x T f t f t T ======∑∑∑. 则（）等价于

0111()

() +(), (0), 1n

n n

i i i i i i i i i i i dx t T x t AT f t T x x i n dt =====≤≤∑∑∑， 或

011

()

(()()), (0),1n

n

i i i i i i i i i i dx t T x t f t T x x i n dt λ===+=≤≤∑∑， 比较上式两边i T 的系数可得

0()

()(), (0), 1i i i i i i dx t x t f t x x i n dt

λ=+=≤≤. （） （）是n 个一阶线性方程的初始值问题，很容易求出其解.请同学们给出解(),1i x t i n ≤≤的具体表达式.

2.1.2 一个二阶线性微分算子的特征值问题

在这一小节，我们讨论在本章常用的一些特征值问题. 代替上节的有限维线性空间n R 和n 阶实对称矩阵A ，在这儿要用到线性空间[0,]C l 的某个子空间

H 和该子空间上的二阶线性微分算子A . 一般地取

2{()[0,]()H X x C l X x =∈在0,x l =满足齐次边界条件}. （）

下面我们讨论二阶线性微分算子2

2d A dx

=-的特征值问题. 先取边界条件为

(0)0,()0X X l ==，设()X x H ∈是A 的特征函数，即()0X x ≠且满足

()()AX x X x λ=.

此问题等价于()X x 是下面问题的非零解

"()()0, 0

(0)()0 .X x X x x l X X l λ?+=<

==?

（） （）便是二阶线性微分算子2

2d A dx

=-的特征值问题，即要找出所有使得该

问题有非零解的λ. 下面求解特征值问题（）.

首先证明要使（）具有非零解，λ必须非负.

设)(x X 是相应于λ的一个非零解，用)(x X 乘（）中的方程，并在[]l ,0上积

分得

0)()()()("=+x X x X x X x X λ，

0)()()( 0 2 0

"

=+??

dx x X dx x X x X l

l

λ，

0)())(()()( 0

2 0

2'0'=+-??dx x X dx x X x X x X l

l

l λ.

由于0)()0(==l X X ，故有

2'2 0

()(())l

l

X x dx X x dx λ=??，

'2

2 0

(())()0l

l

X x dx

X x dx λ=≥??

. （）

当0λ=时，方程0)()("=+x X x X λ的通解为12()X x c c x =+. 利用边界条件

0)()0(==l X X 可得120c c ==，即()0X x =. 因此，0λ=不是特征值.

当0λ>时，方程0)()("=+x X x X λ的通解为

x C x C x X λλsin cos )(21+=. （） 利用边界条件0)()0(==l X X 确定常数21,C C 如下

10C =, l C l C λλsin cos 021+=,

或

0sin 2=l C λ.

由于要求（）中齐次微分方程的非零解，故2C 不能为零. 故有

0sin =l λ.

0>，从而有

πλn l = ， 1n ≥,

2

)(

l

n n πλ=，1n ≥ . 将n C C λ,,21代入到（）中，并略去任意非零常数2C 得

x l

n x X n π

sin

)(= ， 1n ≥ . 故特征值问题（）的解为

2)(

l n n πλ= ， x l

n x X n πsin )(= ，1n ≥ （） 注1 特征值问题是分离变量法的理论基础. 上面已求出特征值问题（）的

解为{ sin 1 }n x n l

π

≥. 在高等数学中知道，

在一定条件下区间[0 , ]l 的任一函数可按特征函数系{ sin 1 }n x n l

π

≥展开为Fourier 级数. 换言之，特征函数系

{ sin 1 }n x n l π≥是区间[0 , ]l 上满足一定条件的函数所成无穷维空间的一组基，

而且还是该空间上的一组正交基，即有0

sin

sin 0 , l

n m x xdx n m l l

ππ=≠?. 特征函数系{ sin

1 }n x n l

π

≥的这两个根本性质：正交性和完备性（基）

，和定理有限维空间n R 中相应结论很相似，只是现在的特征值和特征函数是无穷个. 另外，若改变（）中的边界条件，其相应的特征值和特征函数也会有所变化. 如将边界条件变为(0)0,'()0X X l ==，则特征值和特征函数分别为

2(21)(21)(

),()sin ,022n n n n X x x n l l

ππ

λ++==≥. 该特征函数系(21){ sin 1 }2n x n l π+≥也具有和特征函数系{ sin 1 }n x n l

π

≥类似

的性质，既正交性和完备性.此类问题的一般结果便是著名的Sturm —Liouville

定理，有兴趣的同学可参阅参考文献[1][4]-. 将以上的结果以定理的形式给出.

定理 [1],[4] 考虑二阶线性微分算子2

2d A dx

=-的特征值问题

"()()

()()0 , 0 ,

(0)0,()0 .

k m X x X x x l X X l λ?+=<

120......n λλλ≤<<<<→∞.

相应的特征函数系1{()}n n X x ≥在[0,]l 上是相互正交的. 且对于任一在区间[0,]l 上分段光滑的函数()f x ，可按特征函数系1{()}n n X x ≥展开为如下的Fourier 级数 1()()n n n f x f X x ∞

==∑，

其中Fourier 系数为

20

()(), 1()l

n

n l n f x X

x dx

f n X x dx

=

≥??.

为后面需要，下面再求解二阶线性微分算子2

2d A dx

=-带有周期边界条件的

特征值问题. 在偏微分方程教材中，习惯上用()θΦ表示周期函数，即考虑下面

二阶线性微分算子2

2d A dx

=-的周期边值问题

"()()0, () (2), .

θλθθθπθθ?Φ+Φ=-∞<<+∞

?

Φ=Φ+-∞<<+∞? （） 可证（）和以下问题等价

"''

()()0, 02(0) (2), (0) (2).

θλθθπ

ππ?Φ+Φ=≤≤??Φ=ΦΦ=Φ?? （） 和（）的证明相似易得（）中的特征值0≥λ. 当0λ=时，12()c c θθΦ=+, 由周期边界条件可得20c =. 所以0()1θΦ=为特征函数.

当0λ>时，方程通解为

θλθλθsin cos )(21c c +=Φ,

求导得

'()c c θΦ=-+.

由周期边界条件可得

112cos(2sin(2c c c c c c ?=+??

=-+?? 或

1212[1cos(2sin(20

sin(2[1cos(20.

c c c c π?--=??

+-=?? （）

由于要求非零解，故12,c c 不能同时为零. 因此，齐次方程组（）的系数矩阵行列式必为零，即

1cos(20π-=. 解之可得

2n n =λ，()cos sin .n n n c n d n θθθΦ=+

此时对每个正特征值2n n =λ，特征函数有二个，既θn cos ，θn sin . 总结所得结果为如下定理.

定理 考虑二阶线性微分算子2

2d A d θ

=-带有周期边界条件的特征值问题

"''

()()0, 02(0) (2), (0) (2).

θλθθπ

ππ?Φ+Φ=≤≤??Φ=ΦΦ=Φ?? 则该问题的特征值和特征函数分别为

00,λ=0()1;θΦ=2n n =λ，(){cos ,sin }, 1n n n n θθθΦ=≥.

§2?2 分离变量法

本节结合具体定解问题的求解来介绍分离变量法（method of separation of

variables ）. 所举例子仅限于一维弦振动方程，一维热传导方程混合问题以及平面上一些特殊区域上的位势方程边值问题. 对高维问题的处理放在其它章节中介绍.

以下多数例子均假定定解问题带有齐次边界条件. 否则,可利用边界条件齐次化方法转化之. 我们以弦振动方程的一个定解问题为例介绍分离变量法.

2． 弦振动方程定解问题

例求解两端固定弦振动方程的混合问题

2(,), 0, 0 (2.1)(0,)0, (,)0, 0 (2.2)(,0)(), (,0)(),0. tt xx t u a u f x t x l t u t u l t t u x x u x x x l ?ψ-=<<>==≥==≤≤ (2.3)

??

???

解 分四步求解.

第一步 导出并求解特征值问题. 即由齐次方程和齐次边界条件，利用变量分离法导出该定解问题的特征值问题并求解.

令)()(),(t T x X t x u =，并代入到齐次方程中得

0)()()()(''2''=-t T x X a x X t T ,

或

''''2()()

()()

X x T t X x a T t = . 上式左端是x 的函数而右端是t 的函数，要二者相等，只能等于同一常数.令此常数为-λ，则有

λ-=)()("x X x X ， "2()

()

T t a T t λ=- ，

上面的第一个方程为

0)()("=+x X x X λ .

利用齐次边界条件（），并结合0)(≠t T 得

0)()0(==l X X .

由此便得该定解问题的特征值问题为

"()()0, 0

(0)()0.

X x X x x l X X l λ?+=<

其解为

特征值：2(

) , 1 ;n n n l πλ=≥特征函数： ()sin , 1 .n n X x x n l

π

=≥

第二步 正交分解过程. 即将初值和自由项按特征函数系{}1()n n X x ≥展成Fourier 级数，并将),(t x u 也用特征函数{}1()n n X x ≥表出.

11()()sin

n n n n n n x X x x l

π

???∞∞

====∑∑ , （） 1

1

()()sin

n n n n n n x X x x l

π

ψψψ∞

∞

====∑∑, （） 11(,)()()()sin

n n n n n n f x t f t X x f t x l

π

∞

∞

====∑∑, （） 1

1

(,)()()()sin

n n n n n n u x t T t X x T t x l

π

∞

∞

====∑∑ （） 这里n ?，n ψ和)(t f n 分别为)(x ?，)(x ψ和),(t x f 的Fourier 系数，具体表示如下

02()sin l n n d l l π

??ααα=

?, 02()sin l n n d l l π

ψψααα=?,

02()(,)sin l n n f t f t d l l

π

ααα=?,

而)(t T n 为待定函数.

第三步 待定系数法. 即先将),(t x f 和),(t x u 的Fourier 级数代入到（）中，导出关于)(t T n 满足的常微分方程. 再利用初值条件（）得出)(t T n 满足的初始条

件.

假设中的级数可逐项求导，并将（）和（）代入到（）中得

"2

"111

()()()()()()n n

n n

n n n n n T t X

x a

T t X

x f t X x ∞

∞

∞

===-=∑∑∑,

"2

1

1

1

()()()(())()()n n

n

n

n

n n n n n T t X

x a

T t X

x f t X x λ∞

∞

∞

===--=∑∑∑,

"2

1

1

(()())()()()n

n n n n n n n T t a T t X x f t X x λ∞

∞

==+=∑∑. （）

由于Fourier 展式是唯一的，比较（）两端)(x X n 系数得

"2

()()(), 1.n n n n T t a T t f t n λ+=≥ （）

在（）中令0=t 并结合（）得

()(0)()()n n n n n n x T X x X x ??∞

∞

====∑∑ （）

比较（）两端)(x X n 系数得

(0), 1.n n T n ?=≥ （） 类似地可得

'(0), 1.n n T n ψ=≥ （） 结合（），（）和（）便得出关于)(t T n (1)n ≥满足的二阶常系数非齐次方程初始值问题

"2

'

()()(), 0

(0), (0).

n n n n n n n n T t a T t f t t T T λ?ψ?+=>??==?? （） 第四步 求解关于)(t T n 的定解问题（），并将其结果代入到（）中即可. 为简单起见，我们设()0,1n f t n =≥. 将n λ代入到（）中可得方程的通解为

t l

a

n d t l a n c t T n n n ππsin cos

)(+= ， 利用初始条件确定常数,n n c d 如下

'(0), (0)n n n n n

n a

T c T d l

π?ψ====. 故有

()cos

sin n n n l n a n a T t t t l n a l

ψππ?π=+ . 最后将上式代入到（）中便得定解问题（）—（）的解为

1

2(,)()sin cos sin l n n n a n u x t d t x l

l l l

πππ

?ααα∞

==∑?

1

2

()sin

sin sin l

n n n a n d t x n a l l l

πππ

ψαααπ∞

=+∑

?

注1 利用分离变量法求解（）—（），需要假设在（）中可通过无穷求和号

∑逐项求导. 而通过∑号求导要对无穷级数加某些条件，

在这里就不做专门讨论了. 今后遇到此类问题，我们均假设一切运算是可行的，即对求解过程只作形式上的推导而不考虑对问题应加什么条件. 通常称这样得出的解为形式解. 验证形式解是否为真解的问题，属于偏微分方程正则性理论的范围. 一般地讲，偏微分方程定解问题的解大多数是以无穷级数或含参变量积分形式给出的. 对这两类函数可微性的研究需要较深的数学知识，也有一定的难度，有兴趣的同学可查阅参考文献[1]和[2]. 我们约定：本书只求定解问题的形式解.

注2 当(,)0f x t =时，由（）可以看出：两端固定弦振动的解是许多简单振动(,)()sin

n n n u x t T t x l π=的叠加，当(11)k kl

x x k n n

==≤≤-时，对任意的时刻t ，(,)0n k u x t =，即(,)n u x t 在振动的过程中有(1)n +个点永远保持不动，所以称这样

的振动为驻波,而k x 称为该驻波的节点.显然当21

(11)2k x l k n n

+=

≤≤-时sin 1x =，在这些点上振幅最大，称这些点为驻波的腹点. 因此，求特征函数实际上就是求由偏微分方程及边界条件所构定的系统所固有的一切驻波. 利用由系统本身所确定的简单振动来表示一些复杂的振动，便是分类变量法求解波动问题的物理解释.

注3 例的求解方法也叫特征函数法（eigenfunction method ），现已成为固定模式，也具有普适性. 初学者似乎会感到有些繁琐，但随着进一步的学习，同学们就会熟练掌握这一方法. 特征函数法的关键之处是求解偏微分方程定解问题相应的特征值问题，而基本思想就是笛卡尔（Descartes ）坐标系的思想.如在三维空间3R 中，每个向量可由基{,,}i j k 的线性组合表出，两个向量

111222 , a i b j c k a i b j c k αβ=++=++

相等当且仅当在基{,,}i j k 下两个向量的坐标相等. 既121212 , , a a b b c c ===. 与此相类似，在例求解中也是比较方程或初始条件两边()n X x 的系数而得到（）. 与三维空间3R 相比较，例中特征函数系{ sin

1 }n x n l

π

≥相当于3R

中的基{,,}i j k ，而{ T () 1 }n t n ≥也就相当于上面的111{,,}a b c ，即定解问题的解

关于基函数{ sin

1 }n x n l

π

≥的坐标. 因此，在具有可数基的无穷维空间中，特征函数法也称为待定系数法. 例 设有一均匀细弦，其线密度为ρ. 若0x =端为自由端，x l =端固定.初始速度和初始位移分别为零，并受到垂直于弦线的外力作用，其单位长度所受外力为sin t ω. 求此弦的振动. 解 所求定解问题为

21 sin , 0, 0 (0,)0, (,)0, 0

(,0)0, (,0)0, 0.

tt xx x t u a u t x l t u t u l t t u x u x x l ρω-?-=<<>?

==≥??==≤≤? （）

利用特征函数法求解该问题.

情形1 非共振问题，即22, 0n a n ωλ≠≥. 该定解问题的特征值问题为

"

'

()()0, 0

(0)0, ()0.

X x X x x l X X l λ?+=<

2(21)(

)2n n l πλ+= ， (21)()cos 2n n X x x l

π

+= ，0n ≥ 将1sin t ρω-按特征函数{}0)(≥n n x X 展开成Fourier 级数得

11

sin ()()n n n t f t X x ωρ

∞

==∑ , （

021214()sin sin sin sin 2(21)l n n n f t t d t f t l l n ωπααωωρπρ

+=

==+?. 令

0(,)()()n n n u x t T t X x ∞

==∑ （）

完全类似例的求解过程可得，对于任意0, ()n n T t ≥满足下面问题

"2

'

()()sin , 0

(0)0, (0)0.

n n n n n n T t a T t f t t T T λω?+=>??==?? （） 初值问题中齐次方程的通解为

12()cos sin n T t c c =+，

而非齐次方程的一个特解为

22

()sin n

n n f T t t a ωλω

=

-. 因此，（）的通解为

1222

()cos sin sin n

n n f T t c c t a ωλω=++

-. （）

由初始条件可确定出

120, c c ==

. 最后将所得到的()n T t 代入到（）中便得（）的解.

情形2 共振问题，即存在某个0,n ≥ 使得22n a ωλ=.

不妨假设220a ωλ=. 此时，在情形1中求解所得到的{ T () 1 }n t n ≥不变. 当

0n =时，要求解以下问题

"2

000'

00()()sin , 0

(0)0, (0)0.

T t T t f t t T T ωω?+=>??==?? （） （）中齐次方程通解为

012()cos sin T t c t c t ωω=+.

为求得非齐次方程的一个特解，要将（）中方程的自由項换为0i t f e ω，而求以下问题的一个特解

"2000()(). i t T t T t f e ωω+=

令()i t T t Ate ω=并代入到上面非齐次方程中可得 02f i

A ω

=-

，故有 00()sin cos 22f t f t

T t t i t ωωωω

=

-， 取其虚部便得（）中方程的一个特解为

00()Im(())cos 2f t

T t T t t ωω

==-.

结合以上所得结果便可得到（）中方程的通解为

0012()cos sin cos 2f t

T t c t c t t ωωωω

=+-，

由初始条件确定出 01220, 2f

c c ω==，由此可得

0002()sin cos 22f f t

T t t t ωωωω

=-.

将()n T t 代入到（）中便得在共振条件下（）的解为

001

021

12(,)()()

()()()()

(sin cos )cos ()()

222 (,)(,) .

n n n n n n n n n u x t T t X x T t X x T t X x f f t t t x T t X x l u x t u x t π

ωωωω∞

=∞

=∞

===+=-+=+∑∑∑

可以证明: 2(,)u x t 是有界的. 而在1(,)u x t 的表达式中取 2k k t π

ω

=

，则

1(,)u x t 中的基本波函数cos

2x l

π

的振幅0()k T t 当k 逐渐变大时将趋于无穷大，

最终要导致弦线在某一时刻断裂，这种现象在物理上称为共振. 注意到在上面求解

过程中我们取周期外力的频率ω

等于系统的第一固有频率函数分量上发生共振. 一般地讲，当周期外力的频率ω很接近或等于系统的某

个固有频率时，系统都会有共振现象发生，即弦线上一些点的振幅将随着时间的增大而不断变大，导致弦线在某一时刻断裂.

2． 热传导方程定解问题 例 求解下面热方程定解问题

20, 0, 0 (0,), (,)sin , 0(,0)0, 0.t xx x u a u x l t u t u u l t t t u x x l ω?=<<>?

==≥??=≤≤?

（）

解 利用特征函数法求解.

首先将边界条件齐次化,取0(,)sin w x t u x t ω=+，并令w u v -=，则（）转化为

20cos , 0, 0 (0,)0, (,)0, 0(,0), 0.

t xx x v a v x t x l t v t v l t t v x u x l ωω?-=-<<>?

==≥??=-≤≤?

（）

利用分离变量法可得（）的特征值问题为

"()()0, 0

(0)0, '()0.

X x X x x l X X l λ?+=<

==? 特征值和特征函数分别为

2

(21)(

)2n n l

πλ+= 0≥n ,

(21)()sin

2n n X x x l

π

+= 0≥n . 将(,)cos f x t x t ωω=-，

0)(u x -=?按特征函数{}0)(≥n n x X 展成Fourier 级数得 0cos ()()n n n x t f t X x ωω∞

=-=∑ ，

02(21)()(1)cos sin cos 2l n n n f t t d f t l l

πωαωααω+=

-=?, 其中

1

22

8(1)(12)n n l f n ωπ+-=+.

00

n n n u X ?∞

=-=∑，

其中

00

042(21)()sin 2(12)l n u n u d l l n π

?ααπ

-+=

-=+? . 令

0(,)()(), n n n v x t T x X x ∞

==∑

并将（）代入到（）中的方程得

'2

"00

()()()()cos ()n n

n n

n n n n n T t X

x a

T t X

x f tX x ω∞

∞

∞

===-=∑∑∑ ,

'20

(()())()cos ()n

n n

n

n n n n T t a T t X

x f tX x λω∞

∞

==+=∑∑.

在（）中令0=t 并结合（）得

()(0)()()n n n n n n x T X x X x ??∞

∞

====∑∑.

比较上面两式中特征函数()n X x 的系数便得

'2()()cos , 0(0).

n n n n n n T t a T t f t t T λω??+=>??=?? （） （）是一阶常系数常微分方程初值问题.齐次方程通解为

t a n n Ce t T λ2

)(-=.

令()cos sin n T t A t B t ωω=+，并利用待定系数法求特解可得

2242242

()cos sin n n n

n n n a f f T t t t a a λωωωωλωλ=+++,

故有

22242242()cos sin n a t

n n n

n n n

a f f T t Ce t t a a λλωωωωλωλ-=++++ （） 在上式中代0t =得

2242n n

n n a f C a λ?ωλ=++,

22

42

n n

n n

a f C a λ?ωλ=-+ . 最后将代入到中便得的解为

0(21)(,)()sin

2n n n v x t T t x l

π

∞

=+=∑ . 故（）的解为

),(),(),(t x w t x v t x u +=

0 (,)sin v x t u x t ω=++

其中)(t T n 由（）给出.

2． 平面上位势方程边值问题

考虑矩形域上Poisson 方程边值问题

1212(,), , (,)(), (,)(), (,)(), (,)(), .

xx yy u u f x y a x b c y d u a y g y u b y g y c y d u x c f x u x d f x a x b +=<<<

==≤≤??==≤≤? （）

我们假设0)()(21==x f x f 或0)()(21==y g y g . 否则，利用边界条件齐次化方法化非齐次边界条件为齐次边界条件. 当然，也可以利用叠加原理将（）分解为二个问题，其中一个关于x 具有齐次边界条件，而另一个关于y 具有齐次边界条件. 例 求解Dirichlet 问题

0, 02, 0 1 (0,)0, (2,)0, 01(,0)1, (,1)(1), 0 2.xx yy u u x y u y u y y u x u x x x x +=<<<

==≤≤??==-≤≤?

（）

解 令)()(),(y Y x X y x u =并将其代入到（）中齐次方程得

0)()()()(""=+y Y x X y Y x X ,

λ-=-=)

()

()()(""y Y y Y x X x X , "()()0, 0 2

(0)0, (2)0. X x X x x X X λ?+=<

0)()("=-y Y y Y λ （） （）便是（）的特征值问题，其解为

2)2(πλn n =， x n x X n 2sin )(π

=， 1≥n .

将n λ代入到（）中得

0)()("=-y Y y Y n λ ， （） 该方程有两个线性无关解y n e 2

π，y n e

2

π-

. 由于2n sh

y π，2

n ch y π也是（）的解且线性无关，故（）通解为

y n ch d y n sh

c y Y n n n 2

2)(π

π+=. 令

1

1

(,)()()()sin 222

n n n n n n n n n u x y X x Y y c sh

y d ch y x πππ

∞

∞

====+∑∑ 则),(y x u 满足（）中方程和关于x 的齐次边界条件. 利用关于y 的边界条件可如下确定n c ，n d ，

∑∞

==12

sin

1n n x n d π

, ))1(1(2

2sin

12220n n n d n d --=?=?π

ααπ. （） x n n ch d n sh

c x x n n n ∑∞

=+=-12

sin )22()1(π

ππ,

22233

0216(1)164(1)(1)sin ,

2222n n

n n n n n n c sh d ch d n ππππααααπ----+=-=?

2

2)

)1(1(22)1(416)1(16332

2

ππ

ππππn sh n ch

n n sh n n c n n

n

n -------= . （） 故（）解为

1

(,)()sin ,222

n n n n n n u x y c sh

y d ch y x πππ∞

==+∑ （） 其中n c ，n d 由（）和（）确定.

对于圆域，扇形域和圆环域上的Poisson 方程边值问题，求解方法和矩形域

上的定解问题无本质区别，只是在此时要利用极坐标.

同学们自己可验证：令θρcos =x ，θρsin =y 作自变量变换，则有

θθρρρρ

ρu u u u u yy xx 2

11+

+

=+.

令)()(),(θρθρΦ=R u ，将其代入到极坐标下的Laplace 方程中得

"'"2

1

1

()()()()()()0R R R ρθρθρθρρ

Φ+

Φ+

Φ=,

"'"2

1

1(()())()()()0R R R ρρθρθρ

ρ+

Φ+

Φ=,

"'"

2

1

()()

()

1

()

()

R R R ρρθρ

λθρρ+

Φ=-=-Φ,

故有

0)()("=Φ+Φθλθ, （） 0)()()('"2=-+ρλρρρρR R R . （） 方程（）结合一定的边界条件便得相应定解问题的特征值问题,而 （）是欧拉（Euler ）方程. 对（）作自变量变换s e =ρ可得

s e =ρ ， ρln =s ,

'1

s dR dR ds R d ds d ρρρ

==,

2222'''2222211()ss s d R d R ds dR d s R R d ds d ds d ρρρρρ

=+=-. 将以上各式代入到（）得

''0ss R R λ-= . （） 例 求下面扇形域上Dirichlet 问题

22220, 0, 0, 4

(,0)0, 0 2

(0,)0, 0 2 (,), 4. xx yy u u x y x y u x x u y y u x y xy x y ?+=>>+

=≤≤??=≤≤?

?=+=? （）

的有界解.

解 令θρcos =x ，θρsin =y 作自变量变换，（）转化为

2110, 0, 0 2 2

(,0)0, (,)0, 022(2,)2sin 2, 0.2u u u u u u ρρρθθ

πθρρρπρρρπθθθ?

++=<<<

?

==≤≤??

?

=≤≤??

（）

令)()(),(θρθρΦ=R u 代入到（）中的方程，并结合边界条件可得

"()()0, 0<

(0)0, (/2)0.θλθθππ?Φ+Φ=?Φ=Φ=? （）

0)()()('"2=-+ρλρρρρR R R . （） （）便是（）的特征值问题.

求解特征值问题（）可得

22

4)2/(n n n ==ππλ ， θθn n 2sin )(=Φ ， 1≥n .

将n λ代入到（）中，并令s e =ρ作自变量变换可得

"240ss R n R -=,

2222()ns ns n n n n n n n R c e d e c d ρρρ--=+=+.

由于是求（）的有界解，故有∞<)0(R ，即0=n d . 从而有

n n n c R 2)(ρρ=.

上面求出的(,)()()n n n u R ρθρθ=Φ对每个1n ≥都满足（）中的方程和齐次边界条件，由叠加原理得

∑∑∞

=∞

==Φ=1

21

2sin )()(),(n n n n n n n c R u θρθρθρ, （）

也满足（）中的方程和齐次边界条件.为使（）中的非齐次边界条件

(2,)2sin u θθ=得以满足，在（）中令2ρ=得

21

2sin 22sin 2n n n c n θθ∞

==∑， （）

比较上式两边特征函数θθn n 2sin )(=Φ的系数得 11

2

c =

， 1)( 0≠=n c n . 将1c ，1)(≠n c n 代入到（）中便得（）的解为 θρθρ2sin 2

1),(2

=u . 例 求解圆域上Dirichlet 问题

2110, 0, 02 (,)(), 02.u u u a u a ρρρθθ

ρθπρρθ?θθπ?++=<<≤

)()(),(θρθρΦ=R u 并代入到（）中的方程可得

"()()0

() (2).θλθθπθ?Φ+Φ=?Φ=Φ+? （）

0)()()('"2=-+ρλρρρρR R R . （） （）是定解问题（）的特征值问题. 由定理知（）的解为

2, ()cos sin , 0n n n n n c n d n n λθθθ=Φ=+≥.

将n λ代入到（）中可得(要利用自然边界条件(0,)u θ<∞)

00)(c R =ρ,n n n c R ρρ=)(，1≥n

利用叠加原理可得（）的如下形式解

∑∞

=++=10)sin cos (),(n n n n n d n c c u θθρθρ . （）

根据边界条件)(),(θ?θ=a u 得

01

()(cos sin )n n n n c a c n d n ?θθθ∞

==++∑,

其中

2001()2c d π

?ττπ

=

?, ?=π

τττ?π20cos )(1d n a c n n ,

?=π

τττ?π20

sin )(1d n a d n n .

将以上各式代入到（）中便得（）的解为

2 2 0

011

1(,)()()(()cos cos 2n n u d n d n a

π

π

ρρθ?ττ?τττθπ

π∞

==

+∑?

?

)sin sin )(1

2 0

?

+

π

θτττ?π

n d n . （）

注4 利用等式)Re()(cos 1

)(1

∑∑∞

=-∞

==-n in n n n e c n c τθτθ可将（）化为如下形式

222220

1

()()

(,),22cos()

a u d a a πρ?τρθτπ

ρρθτ-=

+--?

（）

式（）称为圆域上调和函数的Poisson 公式. 在后面学习中还将用其它方法导出它.

注5 在例和例中，如果方程中自由项),(θρf 不为零，若),(θρf 特殊，可用函数代换将自由项化为零而转化齐次方程. 对于一般的),(θρf ，要利用特征函数方法求解.

注6 上面例—例几个定解问题的求解思想和主要过程，是伟大的数学家和物理学家Fourier 给出的，详细内容见参考文献[5]. 在这部著名论著中，Fourier 首次利用偏微分方程来研究热问题，并系统地介绍了分离变量法的基本思想和主要步骤.

结合本节所举例子，请同学们小结一下在本章所学过的特征值问题，二阶常系数非齐次常微分方程和欧拉方程的求解方法. 习 题 二 1. 设有如下定解问题