X"Y
XY
X"X0
流
ny
ny
n (na)2,n1,2, 固有
程
Yn(y)Anea Bnea
Xn(x) sinna x,
值 （特
图
unXn(x)Yn(y)
征值） 问题
u u(x, y)
u (x ,y ) X n (x )Y n (y )
设形式解为：
u(x,y)X (x)Y(y)
代入上述泛定方程,得到
特点：边界条件 均非齐次
u x 0 p 0 u x a P(0 y b )
uy 0u 0
u U(0xa ) y b
让 ( x , y ) 和 w ( x , y ) 分别满足拉普拉斯方程,并各有
一组齐次边界条件，即
xx x
0
y
y
p
0
0
xa
P
y 0 0
0 yb
w w
xx
x0
【解】 设 u(x,t)X(x)T(t)并代入方程得
XTa2XT0
X (0)T (t) 0
X (l)T (t) 0
X(0) 0
X(l) 0
u |t0 ( x ) ut t0 (x)
u a2u
tt
xx
ux |x0ux |xl=0
分 离 变 量
uT(t)X(x)
X(0)X(l)0
T a2T
注：该解称为古典解，在求解中我们假设无穷级数是收敛的。
如上的方法称为分离变量法，是齐次发展方程求解的一个 有效方法。下面对该方法的步骤进行总结。
u |t0 ( x ) ut t0 (x)
u a2u
tt
xx
u|x0u|xl=0
分 离 变 量
uT(t)X(x)
X(0)X(l)0
T a2T
X" X
【解】杆上温度满足下列泛定方程和定解条件
ut a2uxx 0, 0 x l, t 0 u x00, ux xl 0
u t0(x)
分析：方程与边界条件均为齐次，用分离变量法，根据分离变量法流程，分析如下
试探解 u(x,t)X(x)T(t)
代入方程和边界条件得 固有值问题
XX 0
X(0)0 X(l)0
和常微分方程
方程左边是傅里叶正弦级数,这就提示我们把右边 的展开为傅里叶正弦级数，然后比较傅里叶系数，得
An2 l0l()sinnld
Bnn 2a0l()sin nl d
按上述公式计算出系数 A n 和 B n ，则可得原问题的解：
u(x,t) n 1(Anconsla tBnsin nla)tsin
nx
l
Ta2T0
u|t0(x)
u a2u
t
xx
u|x0ux |xl=0
分 离 变 量
uT(t)X(x)
T a2T
X" X
Ta2T0
X(0)X(l)0
X"X0
流 程
Tn(t)Cne(2n14)l222a2t
n ((2n2l1))2,n1,2, 固有
Xn sin(2n2l1)x,
值 （特
图
unTn(t)Xn(x)
Ta2T0
X"X0
流程 TnAncosnlatBnsinnlat
n
(
n l
)2
固有
X
n
sin
n l
x,
值 （特
图
unTn(t)Xn(x)
征值） 问题
uu(x,t)
u Tn(t)Xn(x)
偏微分方程 分 离 变量
常 微 分 方 程 （ 关 于 X ） + 边 界 条 件 故 有 ( 值 ) 函 数 常 微 分 方 程 （ 关 于 T ） + 初 始 条 件 叠 加 系 数
由第一式可得 C 2 0 ，因此第二式变为
C1sinl 0 而 C1 0 只有
sin l 0
l n
于是有固有值和固有函数
n n2 2 l2
Xn(x)
cos
n
l
x
n1,2,
综上所述，该问题的固有值和固有函数分别为
n n2 2 l2
Xn(x)
cos
n
l
x
n0,1,2,
现在需要求解 Ta2n2T0
当 n 0 时有解 T0(t)A0B0t
当 n 0 时有解
n a t n a t
T n (t) A nc o sl B nsinl (n 1 ,2 , ) 其中 A0 B0 An Bn 均为独立的任意常数。
由叠加原理，一般解为
u(x,t)＝ Tn(t)Xn(x)
＝ A n=0 0 B 0 t n 1 (A nc o sn la t B ns in n la t)c o sn lx
度分布 u(x, y) .
【解】先写出定解问题定解问题
uxx uyy 0
方程齐次
这组边界条件齐次
u 0u 0(0 y b )
x 0
x a
uy 0u 0
u U(0xa ) y b
用分离变量法
u |y0 u0 u U
yb
uxx uyy 0 u|x0u|xa=0
分 离 变 量
YY0
uX(x)Y(y) X(0)X(a)0
为确定叠加系数，将 u ( x , y ) 代入非齐次边界条件
n1
(An
Bn
)sin
nx
a
u0
n1
n
(Ane a
b
n
Bne a
b
)sin
nx
a
U
将等式右边展开为傅里叶正弦级数,并两边比较系数，得
A n B na 20 au 0sinn a xd x2 u 0(1 n ( 1 )n)
AnenabBnenab2U(1n (1)n)
联立求解得
An
(1(1)n)(U
nb
u0e a
nsh(nb)
)
a
nb
Bn
(1(1)n)(u0e a
nsh(nb)
U)
a
故原问题的解为
ny
u(x,y) (Anea
n1
Bnenay)sinnax
2(1(1)n) ny
n(by) nx
n1
nshnb [Ush
a
u0sh
a
]sin a
a
小结：对矩形域上拉普拉斯方程，只要一组边界条件
于是 Cn2 l 0l(x)sin(2n 2l1)xdx
如取 ( x ) A x ，则 l
Cn2 l 0l A lxsin(2n 2l1)xdx2l2 A[x((2n2l1)cos(2n 2l1)x)1 0
(2n2l1)0lcos(2n 2l1)xdx](1)n(2n 8A 1)22
从而下列问题 的解为
由初始条件得
A0
B0
n1
n1
An Bn
cosnx(x)
l
nacosnx(x)
l
l
(0xl)
把右边的函数展成傅里叶余弦级数, 比较两边的系数，得
A0
1 l
l()d
0
B0
1 l
l()d
0
An2 l0l()conlsd Bnn 2a0 l()conlsd
【例题2】研究细杆导热问题,初始时刻杆的一端温度为零度, 另一端跟外界绝热，杆上初始温度为 ( x ) ,试求无热源时细杆上 温度的变化。
XX 0
得到固有值问题 X(0)0 X(a)0
和常微分方程
YY0
得固有值：
n na222
(n1,2,.....)
固有函数：
Xn(x)
sin
nx
a
(n1,2,.....)
ny
ny
而
Yn(y)Anea Bnea
于是有
un(x,y)(A nena yBnena y)sinnax
叠加得
u(x,y)n 1(A nen a yB nen a y)sinnax
是齐次的，则可使用分离变量法求解。
图形如下： （程序：my2）
u u
6
5
4
6
3
5
2
4
1
3
2
0 3
2
2
1 2
1.5
0
1.5 1
1 Y
00
1 0.5
X
3 2.5
2 1.5
1 0.5
Y
0.5 X
00
a=2 b=3 U0=1 U=5 N=150
(a) 精确解图
(b) 瀑布图
【例2】求解下列问题 uxx uyy 0
X"X0
X (0 )0X (l)0
经讨论，当 0 时有解
XC1xC2
由定解条件得 C1 0,C2任意，于是有固有值和固有函数
00 , X 0(x)C 2 或X0(x)1
当 0 时有解
X (x ) C 1cox s C 2sin x
现确定积分常数，由条件知
C2 0 (C1sinlC2cosl)0
只有sin l 0才能保证 C 2 0 ，方程有非零解
于是 l n 或 n 此时 Xn(x)Cnsinnl x
n 2 l2
2
n1,2,
n 称为固有值，X n ( x ) 称为固有函数
再看关于T 的方程
T"
a2
n22
l2
T
0
分这离个变方量程的的形解式解T n(t)A n cosnlatB n sinnlat
25
25