数理方程习题综合

例 1.1.1 设v=v(线x,y),二阶性偏微分方程v xy =xy 的通解。 解 原方程可以写成 e／ex（ev／ey） =xy 两边对x 积分，得 v y =￠（y ）+1/2 x 2Y ,

其中￠（y ）是任意一阶可微函数。进一步地，两边对y 积分，得方程得通解为

v （x,y ）=∫v y dy+f （x ）=∫￠（y ）dy+f （x ）+1/4 x 2y 2

=f （x ）+g （y ）+1/4 x 2y 2

其中f （x ）,g （y ）是任意两个二阶可微函数。 例1.1.2

即 u(ξ,η) = F(ξ) + G(η),

其中F(ξ),G(η)是任意两个可微函数。

例1.2.1设有一根长为L 的均匀柔软富有弹性的细弦，平衡时沿直线拉紧，在受到初始小扰动下，作微小横振动。试确定该弦的运动方程。

取定弦的运动平面坐标系是O XU ,弦的平衡位置为x 轴，弦的长度为L ，两端固定在O,L 两点。用u(x,t)表示弦上横坐标为x 点在时刻t 的位移。由于弦做微小横振动，故u x ≈0.因此α≈0，cos α≈1,sin α≈tan α=u x ≈0,其中α表示在x 处切线方向同x 轴的夹角。下面用微元法建立u 所满足的偏微分方程。

在弦上任取一段弧'MM ,考虑作用在这段弧上的力。作用在这段弧上的力有张力和外力。可以证明，张力T 是一个常数，即T 与位置x 和时间t 的变化无关。

事实上，因为弧振动微小，则弧段'MM 的弧长

dx u x

x x

x ?

?++=?2

1s ≈x ?。

这说明该段弧在整个振动过程中始终未发生伸长变化。于是由Hooke 定律，张力T 与

时间t 无关。

因为弦只作横振动，在x 轴方向没有位移，故合力在x 方向上的分量为零，即

T(x+x ?)cos α’-T(x)cos α=0.

由于co's α’≈1，cos α≈1,所以T(X+?x)=T(x),故张力T 与x 无关。于是，张力是一个与位置x 和时间t 无关的常数，仍记为T.

作用于小弧段'MM 的张力沿u 轴方向的分量为 Tsin α’-T sin α≈T(u x (x+x ?,t)-u x (x,t)).

设作用在该段弧上的外力密度函数为F （x,t ）那么弧段'MM 在时刻t 所受沿u 轴方向的外力近似的等于F(x,t)x ?.由牛顿第二定律得

T （u x (x+x ?,t)-u x (x,t)+F(x,t)x ?=ρx ?tt u ,

其中ρ是线密度，由于弦是均匀的，故ρ为常数。这里tt u 是加速度tt u 在弧段'MM 上的平均值。设u=u(x,t)二次连续可微。由微分中值定理得

Tu zz （x+θx ?，t)x ?+F(x,t)x ?=ρtt u x ?, 0<θ<1. 消去x ?，并取极限x ?→0得 Tu xx （x,t ）+F(x,t)=ρu tt , 即

u tt =ɑ2

u xx +?(x,t), 00,

其中常数ɑ2

=T/ρ，函数?（x,t ）=F(x,t)/ρ表示在x 处单位质量上所受的外力。 上式表示在外力作用下弦的振动规律，称为弦的强迫横振动方程，又称一维非齐次波动方程。当外力作用为零时，即?=0时，方程称为弦的自由横振动方程。

类似地，有二维波动方程

u tt =ɑ2

（u xx +u y y ）+?（x.y.t ）, (x,y)Ω∈,t>0, 电场E 和磁场H 满足三维波动方程

E c E 2222t ?=??和H c H 2

22

2t

?=??, 其中c 是光速和

22

22222

x z

y ??+??+??=?=???=?。

例1.2.2设物体Ω在内无热源。在Ω中任取一闭曲面S （图1.2）。以函数u(x,y,z,t)表示

物体在t 时刻，M=M(x,y,z)处的温度。根据Fourier 热传导定律，在无穷小时段dt 内流过物体的一个无穷小面积dS 的热量dQ 与时间dt ，曲面面积dS 以及物体温度u 沿曲面的外法线n 的方向导数三者成正比，即

dSdt n u

k

-??,

其中k=k(x,y,z)是在物体M(x,y,z)处的热传导系数，取正值。我们规定外法线n 方向所指的那一侧为正侧。上式中负号的出现是由于热量由温度高的地方流向温度低得地方。故当

0n

u

>??时，热量实际上是向-n 方向流去。

对于Ω内任一封闭曲面S ，设其所包围的空间区域为V ，那从时刻t 1

到时刻t

2

经曲面

流出的热量为

1Q =dSdt n

u

k

S

?????2

1t t - 设物体的比热容为c(x,y,z)，密度为ρ(x,y,z)，则在区域V 内，温度由u(x,y,z,1t )到u(x,y,z)所需的热量为

[]

dvdt t

u

c dv t z y x u t z y x u c t t V

V

??=-=???????2

1),,,(),,,(Q 122ρ

ρ. 根据热量守恒定律，有

12Q Q -=

即

[]dSst n

u

k

dv t z y x u t z y x t t S

????????=-2

1),,,(),,,u c 12V

（ρ 假设函数u(x,y,z,t)关于x,y,z 具有二阶连续偏导数，关于t 具有一阶连续偏导数，那么由高斯公式得

0][2

1t =??

? ??????-???? ??????-???? ??????-??????dvdt z u k z y u k y y u k x t u c t V

ρ

. 由于时间间隔[]21t ,t 及区域V 是任意的，且被积函数是连续的，因此在任何时刻t ，在Ω内任意一点都有

??

? ??????+???? ??????+???? ??????=??z u k z y u k y y u k y x u ρ

c

(1.2.6)

方程称为非均匀的各向同性体的热传导方程。如果物体是均匀的，此时k,c 及ρ均为常数，令2

a =

ρ

c k

，则方程(1.2.6)化为 u a z u y u x u ?=???

? ????+??+??=??2

222222

2a t u , (1.2.7) 它称为三维热传导方程

若物体内有热源，其热源密度函数为，则有热源的热传导方程为

),,,(a u 2t t z y x f u +?= (1.2.8)

其中ρ

c F f =

类似地，当考虑的物体是一根均匀细杆时如果它的侧面绝热且在同一截面上的温度分布相同，那么温度只与有关，方程变成一维热传导方程

xx u 2t a u = （1.2.9）

同样，如果考虑一块薄板的热传导，并且薄板的侧面绝热，则可得二维热传导方程 ）

yy 2t u (u +=xx u a （1.2.10）

（P16)例1.3.1一长为L 的弹性杆，一端固定，另一端被拉离平衡位置b 而静止，放手任其振动。试写出杆振动的定解问题。

解 取如图1.3所示的坐标系。

O L L+b x

泛定方程就是一维波动方程（杆的纵振动方程） u tt =a 2

u xx , 0

在初始时刻（即放手之时），杆振动的速度为零，即u t (x,0)=0,0≤x ≤L.

而在x=L 端拉离平衡位置，使整个弹性杆伸长了b 。这个b 是来自整个杆各部分伸长后的贡献，而不是x=L 一端伸长的贡献，故整个弹性杆的初始位移为

u|0=t =

L

b

x, 0≤x ≤L. 再看边界条件。一端x=0固定，即该端位移为零，故有u(0,t)=0,0≤x ≤L.另一端由于放手任其振动时未受外力，故有u x (L,t)=0,t ≥0.所以，所求杆振动的定解问题为

设弦在x 0点受到横向力T 作用后发生的位移为h,则弦的初始位移为 hx, 0≤x ≤x 0, u(x,0)= x 0

h(L-x), x 0≤x ≤L, L-x 0

其中h 待求。由牛顿第二定律得

F-Tsin α1-Tsin α2=0， 在微小振动的情况下，

Sin α1≈tan α1= h , sin α2≈tan α2= h ,

x

0 L-x 0

所以 F=Th +Th x 0 L-x 0

因此 h=Fx 0(L-x 0) .

TL F(L-x 0) , 0≤x ≤x 0, 从而初始位移为u(x,0)= TL Fx 0(L-x) , x 0≤x ≤L. TL 而初始速度u t (x,0)=0.

u tt =a 2

u xx , 00, u(x,0)=

L

b

x, u t （x,0)=0, 0≤x ≤L, u(0,t)=0, u x (L,t)=0, t ≥0.

（P17）例1.3.2 ：长为L 的均匀弦，两端x=0和x=L 固定，弦中张力为T ，在x=x0处以横向力F 拉弦，达到稳定后放手任其振动。试写出初始条件。

解：建立如图坐标系。

(P18)例1.3.3考虑长为L 的均匀细杆的热传导问题。若（1）杆的两端保持零度；（2）杆的两端绝热；（3）杆的一端为恒温零度，另一端绝热。试写出该绝热传导问题在以上三种情况下的边界条件。

解：设杆的温度为u(x,t),则

（1） u （x,t ）=0，u(L,t)=0. （2） 当沿杆长方向有热量流动时，由Fourier 实验定律得

L

x x x

u k

q x

u

k

q ==??-=??=20

1,'

其中q1,q2分别为x=0和x=L 处的热流强度。而杆的两端绝热，这就意味着杆的两端与外界没有热交换，亦没有热量的流动，故有q1=q2=0和

,

0),0(=t x u 0),(=t L u x .

(3)显然，此时有

0),(,0),0(==t L u t u x .

例1.5.1求Poisson 方程Uxx +Uyy =X^2 +XY+Y^2的通解 解:先求出方程的一个特解V=V （x ，y),使其满足

Vxx +Vyy=X^2 +XY+Y^2 由于方程右端是一个二元二次齐次多项式，可设V （x ，y) 具有形式

V(x,y)=aX^4 +bX^3 Y+cY^4,其中a,b,c 是待定常数 Vx=4aX^3+3bX^2 Y Vy=bX^3+4cY^3 Vxx=12aX^2+6bXY Vyy=12cY^2 得Vxx+Vyy=12aX^2 +6bXY+12cY^2=X^2 +XY+Y^2 比较两边系数，可得 a=1/12,b=1/6,c=1/12

于是V （x,y)=1/12(X^4 +2X^3 Y+Y^4）

下面求函数W=W(x,y),使其满足Wxx+Wyy=0.作变量代换e=x,n=iy(以下的偏导的符号记为d)

Ue=du/de=du/dx=Ux Un=du/dn=du/dy *dy/dn=-iy Uee=dUe/de=Uxx Unn=-Uyy 可得Wee-Wnn=0 再作变量代换

s=e+n,t=e-nUe=du/de(s,t)=Us+Ut Un=du/dn=Us-Ut Uee=dUe/de=d(Us+Ut)/de=Uss+Utt+2Ust Unn=dUn/dn=d(Us-Ut)/dn=Uss+Utt-2Ust 那么方程进一步化为Wst=0

其通解为W=f(s)+g(t)=f(e+n)+g(e-n)=f(x+iy)+g(x-iy),其中f,g 是任意两个二阶可微函数。那么根据叠加原理，方程的通解为u(x,y)=V+W=f(x+iy)+g(x-iy)+1/12(X^4+2X^3 Y+Y^4)

(P32)例2.1.1 判断方程U xx +2U xy -3U yy +2U x +6U y =0（2.1.22）的类型，并化简。

解： 因为a 11= 1,a 12= 1,a 22= -3,所以 =a 212-a 11a 22=4>0,故方程为双曲型方程。对应的特征方程组为

,31122

1112212=-+

=a a a a a d d x

y .111

22

1112212-=--=a a a a a d d x y

该方程组的特征曲线（即通解）为

.,321c x y c x y =+=-作自变量变换

x y x y +=-=ηξ,3则

;3ηξu u u xx +-= ,ηξu u u y +=

,69ηηξηξξu u u u xx +-= ,23ηηξηξξu u u u xy +--=

.2ηηξηξξu u u u yy ++=

将上述各式带入方程（2.1.22），得第一种标准形式

.02

1

=-

ηξηu u （2.1.23） 若令,2

,2

η

ξη

ξ-=

+=

t s 则得到第二种标准形式

.0=+--t s tt ss u u u u （2.1.24）

下面对式（2.1.24）进一步化简。令,t s Ve u

μλ+=则

.

22)2(,)2(,

)()(t s t tt tt t s s ss ss t s t t t s s s e V V V u e V V V u e V V u e V V u μλμλμλμλμμλλμλ++++++=++=+=+=，

代入方程，得

.

0)()21()12(22=-+-+-+-+-V V V V V t s tt ss λμμλμλ

我们取，2

1=

=μλ则式（2.1.24）化简为 ,0=-tt ss V V （2.1.25）

该方程不含一阶偏导数项。

例2.1.2

例 2.1.4 求值问题

4y2v xx+2(1-y2)v xy-v yy-2y/(1+y2) (2v x-v y)=0,xεR1,Y>0

V(X,0)=φ（X）,V Y（X,0）=￠（X）,XεR1

的解，其中φ（x）是已知任意二阶可微函数，￠（x）是任意一阶可微函数。

解先把所给方程化为标准型。特征方程组为

dy/dx =-1/2,dy/dx=1/2y^2.

其通解为

x+2y=C1,x-2y^3/3=C

做自变量变换

ξ=x+2y,п=x-2y^3/3,

这样给定的方程化为标准型

V

ξп=0

依次关于п和ξ积分两次，得通解v=F（ξ）+G（п）.代回原自变量x,y得原方程得通解

v？（x,y）=F（x+2y）+G（x-2y^2/3）

其中F,G是任意两个可微函数。进一步，由初始条件得

φ（x）=v（x,0）=F（x）+G（x）,￠（x）=V Y（x,0）=2F’（x）

从而求出

F（x）=F（0）+1/2∫x0￠（t）dt,G（x）=φ（x）-F（0）-1/2∫x0￠（t）dt.

所以原定解问题的解为

v（x,y）=φ（x-2y^3/3）+1/2∫x+2y x-2y^3/3￠（t）dt.

例2.1.3 设常数A,B,C满足B^2-4A C≠0,m1,m2是方程

Am^2+Bm+C=0 ①

的两个根。证明二阶线性偏微分方程

Au xx+Bu xy+Cu yy=0 ②

的通解具有如下形式：

u=u(x,y)=f(m1x+y)+g(m2x+y), ③

其中f,g是任意两个二阶可微函数。

证不失一般性，设A≠0和B^2-4AC>0.其它情况可以类似的处理。

令ξ=m1x+y,η=m2x+y.则

U x=m1uξ +m2uη，u y=uξ+uη，U xx=m1^2uξξ+2m1m2uξη+m2^2uηη

u yy=uξξ+2uξη+uηη ,u xy=m1uξξ+(m1+m2)uξη+uηη

上述式代入②得：

（Am1^2+Bm1+C）uξξ+（Am2^2+Bm2+C）uηη+(2Am1m2+B(m1+m2)+2C)uξη=0④

由题意得

Am1^2+Bm1+C=0，Am2^2+Bm2+C=0，m1+m2=B/A, m1m2=C/A

上述式代入④得

(1/A)(4AC-B^2)uξη=0

又由题意得4AC-B^2≠0

故uξη=0

对该方程两边分别关于和积分，得通解u=f(ξ)+g（η），代回自变量x,y,得方程②的通解是

u=u(x,y)=f(m1x+y)+g(m2x+y), ③

其中f,g是任意两个二阶可微函数。证毕。

端点自由的半无限长的均匀弦振动的定解问题

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

0,0,0,0,0,0,,0,,2≥+∞<≤>+∞<

?

??===+=t x t x t u x x u x x u t x f u a u x t xx tt ?φ (3.1.22) 因为()0,0=t u x ,我们对函数?φ,,f 关于x 做偶延拓。定义()()x t x F φ,,和()x ?如下：

()()().0,

0,,<≥??

?-=x x x x x φφφ ()()()??

?<-≥=.

0,,

0,x x x x x ??? ()()()?

?

?≥<-≥≥=.0,0,,,

0,0,,,t x t x f t x t x f t x F 函数()()()x x t x F ?φ,,,在+∞<<∞-x 上是偶函数。由推论 3.1.1，()t x U ,是关于x 的偶函数，且()().0,0,0==t U t u x x 这样得到定解问题（3.1.22）的解

()()).0,0(,,≥≥=t x t x U t x u 所以，

当at x ≥时，

()()()()()()()

()

???-+--+-++-++=t t a x t a x at

x at x d d f a d a at x at x t x u 0,21

2121,τττξτξξξ?φφ（3.1.23） 当at x <≤0时，

()()()()()()()()()()()()()???????--+------+-++??

????????????????

++??

????

++-++=t a x t t a x t a x a x t x t a t a x x at at

x d d f a d d f d f a d d a x at at x t x u τττττξτξτξτξξτξξξ?ξξ?φφ.,21,,212121,0000

0 （3.1.24） 例4.2.3

端点固定的半无限长的均匀弦振动的定解问题

考虑定解问题

求解上述问题的基本思路是以某种方式延拓函数使其在上也有定义，，这样把半无界区域上的问题转变成上的初值问题。然后利用达朗贝尔公式（3.1.15），求出在上的解u（x,t）。同时使此解u（0,t)满足u（0，t)=0.这样当x限制在上就是我们所要求的半无界区域上的解。

由微积分知识可知，如果一个连续可微函数g（x)在上是奇函数，则必有g（0）=0.因此要使解u=u（x,t)满足u（0,t)=0，只要u(x,t)是x的奇函数便可。而由推论3.1.1，只

要

f(x,t),是x 的奇函数。因此对函数

f

和关于x 作奇延拓。我们定义F

（x,t),

和

如下：

显然函数F

和

在

上是奇函数。然后考虑初值问题

（3.1.17）

由（3,1,15)，问题（3.1.17)的解是

（3.1.18）

所以问题（3.1.16)的解u （x,t)在

上的限制。于是当

时，

（3.1.19)

当

时，

（3.1.20）

例2.2.1确定下列方程标准型

(1)u xx +2u xy -2u xz +2u yy +6u zz =0 (2).0244=++--z y yz xy xx u u u u u 解:(1)方程对应的系数矩阵是 .60102

1111???

?

?

?

?--=A ,06222=++-+zz yy xz xy xx u u u u u 利用线性代数中把对称矩阵化为对角型的方法,我们可选取

????

?

? ?

?--=212

11011

001B , 则令为三阶矩阵这里,E ,E BAB T

=

.

0 .22=++?

??????

?

?+--=????? ??=????? ????ηηξξ?ηξu u u z y x x y x z y x B 则给定的方程化简为

(2)方程对应的系数矩阵是 .01

0102

024

????

? ?

?----=A 因为

,10

0010

001

????

?

?

?-=T BAB 其中

'11

2

10121

0021

???????

?

??--=B 所以取

.

22

2

??????

??

?

?+--+=????

? ??=????? ??z y x

y x

x z y x B ?ηξ

则给定的方程化简为

.0=++-η??ηηξξu u u u

例3.1.1 求解下列初值问题

,0,,9>+∞<<-∞-=--t x e e u u x x xx tt

()().,sin 0,,0,+∞<<-∞==x x x u x x u t

解：利用达朗贝尔公式（3.3.15）得

()()()()()()()

()

??

?-+--+-++

-++=t

t a x t a x at

x at

x d d f a

d a

at x at x t x u 0

,21

21

21,τττ

ζτζζζψφφ

()()

()

()τζζζττζζ

d d

e e d t x t x t t x t x t

x t

x ???

-+---+--+

+-++=

0333361sin 6

13321

,3cosh sinh 9

2

sinh 923sin sin 31t x x t x x +-+

= 易见，解

()t x u ,关于x 是奇函数。

4.2.1波动方程的初边值问题

例4.2.1 设边长为L 的弦，两端固定，作微小横振动。已知初位移为φ（x ）,初始速度为ψ（x ）,试求弦的运动规律。

解： 该物理问题可归为下列定解问题：

?

????

?

?????><<=====0,0,2

0),(),0()()0,(),()0,(t L x xx u a tt u t L u t u x x t u x x u ψφ [1] 设上述问题有非零变量分离解 u(x,t)=X(x)T(t).代入上述问题[1]中得：

X(x)T ``(t)= a 2 X ``

(x) T(t)，

由此设： T ··(t)∕a 2 T(t) = X ··

(x) ∕X(x) =-λ（记-λ为比值常数），并得：

T ··(t)+ λa 2

T(t) =0 [2] X ··

(x)+ λX(x) = 0, [3]

再根据边界条件 u(0,t)=u(L ,t)=0,得：X(0) T(t)= X(L)T(t)=0 , T(t)≠0,则 X(0)= X(L)=0，

由上分析，得：

[4]

{

)()"(0

)()0(=+==x X x X L X X λ （1）λ=-β2

＜0时，方程组[4]的通解为：X(x)=C 1e βx

+C 2 e -βx

,代入X(0)= X(L)=0，

解得常数C1=C2=0,即得零解X(x)=0（u=0），不合初设u 为非零解，舍去； （2）λ=0时，方程组[4]的通解为： X(x)= C 1 x+ C 2, 代入 X(0)= X(L)=0，解得零解X(x)=0（u=0），舍去；

（3）λ=β2

＞0时，方程组［4］的通解为：X(x)= C 1cos βx+C 2sin βx.代入X(0)= X(L)=0，解得C 1=0, C 2sin βL=0

则 λ=λn =βn 2=（n π/L ）2

,n=1,2,... 对应λn 的特征函数为： X n （x ）=C n sin

L

x

n π , n=1,2,..［

5］

将特征值λn 代入［2］得： T ``(t)+ λn a 2

T(t) =0通解为

T n (t)=A n cos L at

n π+B n sin L

at n π ［6］

综上可

知

定

解

问

题

的

变

量

分

离

特

解

为

：

u n （x ，t ）=（a n cos

L at n π+b n sin L

at n π ）sin L x

n π ［7］

其中，a n = A n C n ,b n =B n C n 为任意常

数

，

n=1,2…

根据线性叠加原理，将特解u n （x ，t ）叠加起来，得到通解：

u(x,t)=

∑∞

=1

),(n n t x u =

L

x n L at n b L at n a n n n πππsin )sin cos

(1

+∑∞

=. [8]

由原定解问题： Φ（x ）=u(x,0)=

∑∞

=1

sin

n n L

x

n a π Ψ (x) =u t (x,0)=

L

x

n L a n b n n

ππsin

1

∑∞

=， 可将Φ（x ），Ψ (x)看作是[0,L]上的傅里叶级数，则有：

??

?????==L dx L

x

n x a dx

L x

n x a n b L

n L n πφπ?π00

sin )(2sin )(2

把上面得到的a n ,b n 代入[8]中，得级数通解 u(x,t)=∑∞

=1

),(n n t x u =

L

x

n L at n b L at n a n n n πππsin )sin cos

(1

+∑∞

= ,其中 ?

?

?????==L dx L

x

n x a dx L x

n x a n b L

n L n

πφπ?π00sin )(2sin )(2

经检验，得到的通解u （x,t ）满足关于x 和t 逐项微分二次后一致收敛，因而满足定解问题[1]中方程和相应条件，即通解u （x,t ）存在，是定解问题的解

例4.2.2设长为L ，且两端自由的均匀细杆，作纵振动，且初始位移为φ（x ），初始速度为ψ（x ）。试求杆做自由纵振动的位移规律。

解：

()()()()()()2tt xx t =a ,0,t 0=,0,00,,0,t 0x x x L x x x x L t L μμμφμψμμ??<<>?

?

=≤≤???

?=≥??

x ，0，

令()()(),x t X x T t μ=，代入上式得：

()()()

()

''''2

T t X x a T t X x λ==- 得到两个独立的常微分方程

()()()()''2''

0X 0

T t a T t x X x λλ+=+=

又由边界条件，得()()'

'

00.X X L ==所以特征值问题为

()()()()''''

0,000.X x X x x L X X L λ??==<

当λ<0时，上述问题只有零解。当λ=0时，可得非零的常数解()000X x A =≠.当

2,0λββ=>时，边值问题中方程的通解为

()cos sin X x A x B x ββ=+

由

边

界

条

件

()()''00,B 0Asin L 0

X X L β====得和.因为

0=,0,1,2,.....

n A n L

π

β

≠=，所以。因此得到一系列特征值和对应的特征函数列 ()2

,cos ,n 0,1,2,...n n n n n x X x A L L ππλ??

=== ???

将n λ代入前式，得到相应的

()00,0cos sin ,n 1,2,...n n n C D t n T t n at n at

C D L L ππ+=??

??

=??+=????

因此函数

()()()00+,n 0,cos sin cos ,1,2,...n n n n n a b t x t T t X x n at n at n x a b n L L L μπππ=??

??==??

??+= ???

????

设所求的形式解为

()001,cos

sin cos ,n n n n at n at n x

x t a b t a b L L L

πππμ∞

=?

?=+++ ??

?∑ 其中系数中的初始条件确定，即 ()()01

,0cos

,n

n n x x x a a L

πφμ∞

===+∑

()()01

,0cos t n

n n a n a x x b b L L

ππψμ∞

===+∑

从而得（n=1,2,…）

()()()()00000012,cos

,12,cos L L

n L L

n n x

a x dx a x dx L L L n x

b x dx b x dx L n a L πφφπψψπ??==????????==????

????

例 4.2.4 设有一均匀细杆，长为L ，两端点坐标分别为X=0和X=L.杆的表面绝热，再X=0端保持零度，在X=L 端热量自由发散到温度为零度的介质中去，已知初始温度为(x),Φ求杆上温度分布规律

解： 设U(x,t)表示在x 处，时刻t 时的温度，那么由第1章可知，所给物理问题可以归结为求解下列定解问题：

2x , 0,0,

(x,0)(x),0,u(0,t)0,u (L,t)hu(L,t)0,t 0,

t xx u a u x L t u x L φ= <<>= ≤≤ (4.2.38)=+= ≥

其中常数h>0,函数(x)φ 在[0,L] 满足狄利克雷条件。我们仍然用分离变量法解这个问题。设(x,t)X(x)T(t).u = 将其代入到（4.2.38）方程中，得

2

'(t)''(x)

.(4.2.39)(t)(x)

T X a T X λ==- 从而得到关于(t),X(x)T 的常微分方程

2'(t)a (t)0,(4.2.40)X''(x)X(x)0,x 4.2.41)

T T λλ+= += 0<< (L.由（4.2.38）

中的边界条件，得

(0)0,X'(L)h (L)0.(4.2.42)X X = +=

下面求解由方程（4.2.41）和边界条件（4.2.42）组成的特征值问题当0λ≤ 时，边值问

题（4.2.41）和（4.2.42）只有零解。当

2

,0λββ=> 时，方程（4.2.41）的通解为 (x)Acos x Bsin x X ββ=+

由边界条件（4.2.42），得

0,(cos L hsin L)0.A B βββ= +=

为求特征值和特征函数，设0.B ≠ 所以

cos sin 0,L h L βββ+=

记,L λβ= 则上式可表示为

1tan ,.hL

γαγα= =-

方程（4.2.43）的根可以看作切曲线1tan y γ= 与直线2y αγ= 的交点的横坐标，见图

4.6.

由此可见，他们交点有无穷多个，他们关于原点对称，设方程（4.2.43）的无穷多个正根依次为

120,n γγγ<<<<<

于是边值问题（4.2.41）和（4.2.42）的特征值n λ和相应的特征函数

(x)n X 为

22

n n 2

,(x)B sin x,n 1,2,(4.2.44)n

n n

n X L γλβ

β==

= =

现在证明特征函数系{}12sin ,sin ,,sin ,

n x x x βββ 在[0,L] 是正交

系.记n (x)sin x,n 1,2,

n φβ= = 则m

(x),(x)n φφ 分别满足

n n m m ''(x)(x)0,0x L,(4.2.45)

(x)(x)0,0x L (4.2.46)

n m φλφφλφ+= << += <<

和边界条件（4.2.42）.用m φ 乘以（4.2.46），然后（4.2.46），然后两式相减，并且在[0,L] 上积分，得

n m n m m n 0

n m m n n m m n 0n m n m n m n ()(x)(x)((x)''(x)(x)''(x))((x)'(x)(x)'(x))'dx ((x)'(x)(x)'(x))|

0'(L)(L)(L)'(L)h (L)(L)h (L)(L)0

L

L

n m L

m dx dx

L λλφφφφφφφφφφφφφφφφφφφφφφ-=-=-=- =-=-+=???因为当

n m ≠ 时，n m λλ≠ .所以

m 0

(x)(x)dx 0,n m.(4.2.47)L

n

φ

φ= ≠ ?

即特征函数系}{sin |1n

x n β∞

= 是[0,L] 上的正交函数系.

下面将n λ

λ= 代入到方程（4.2.40），得解

2n (t)A e ,1,2,

(4.2.48)n n a t

T n λ-= =

由此得到满足方程（4.2.38）中的方程和边界条件的一组特解

2

n n n (x,t)X (x)T (t)C e sin ,1,2,

,

n n n a t

x n u λβ-

== = 其中任意常数,1,2,n n n C A B n = = ,由于方程和边界条件是其次的利用叠

加原理，可设定解问题（4.2.38）的形式解为

2n 1

1

(x,t)(x,t)C e sin (4.2.50)n n n n n a t

u u x λβ∞

∞

-===

=

∑∑ 用sin

m x β 乘以式（4.2.50），并且利用}{sin |1

n x n β∞

= 是[0,L] 上的

正交函数性，我们得到

m

1

(x)sin xdx,1,2,(4.2.51)

L

m m

C m L φβ

=

= ?这里

2

sin ,L

m m L xdx β=

?

将是（4.2.51）代入（4.2.49），即得原定解问题（4.2.38）的形式解.

例4.4.1解下列非齐次边界的定解问题

)(u 2tt x f u a xx += 00,

)0,(x u =)(x φ,)0,(x u t =)(x ψ, 0≤x ≤L,

,),(,),0(B t L u A t u == t ≥0。

其中A ，B 是常数。

解：设u （x,t ）=),(t x υ +)(x ω,将其带入到上述的方程中，得

tt u =[]).()(2x f x a xx +''+ωυ