第八章分离变量法
§8.1有界弦的自由震动
8.1.1定解问题
⎧ ⎩ ⎨ ⎪ ⎪ u tt =a 2 u xx ,0<x<l(1)u| x=0 =0,u| x=l =0(2)u| t=0 =φ(x),u t | t=0 =ψ(x)(3)
8.1.2求解
思路:两端固定的弦形成驻波,故可用驻波法(即分离变量法)求解。
1.分离变量:
令u(x,t)=X(x)T(t)(4)
则(1)→{T ′′ −a 2 μT=0(5)X ′′ −μX=0(6)
(2)→X(0)=0,X(l)=0(7)
2.解本征值问题:
{X ′′ −μX=0(6)X(0)=0,X(l)=0(7)
定义:
本征值:(6)中μ不能任意取值只能依据边界条件
(7)取某些特定值称为本征值
本征函数:对于不同的μ值方程(6)所对应的解
本征值问题:求齐次方程带有齐次边界条件的本征值和本征函数的问题
二阶常微分方程的解:
对于y ′′ (x)+py ′ (x)+qy(x)=0
其特征方程为r 2 +pr+q=0
特征方程的解为:r=−p±p 2 −4q − − − − − − √ 2 ={r 1 r 2
①若r 1 ≠r 2 (实数),则y=c 1 e r 1 x +c 2 e r 2 x ;
②若r 1 =r 2 =a(实数),则y=(c 1 x+c 2 )e ax ;
③若r 1 =a+ib,r 2 =a−ib,则y=e ax (c 1 cosbx+c 2 sinbx)
本征值:μ=−n 2 π 2 l 2 ,n=1,2,3,⋯(8)
本征函数:X n (x)=C n sinnπl x(9)
3.求解T n (t)的方程
T ′′ n (t)+a 2 n 2 π 2 l 2 T n (t)=0
T n (t)=A ′ n cosnπal t+B ′ n sinnπal t(10)
4.叠加,用初始条件定系数
u(x,t)=∑ n=1 ∞ (A n cosnπal t+B n sinnπal t)sinnπl x(11)
A n =2l ∫ l 0 φ(α)sinnπl αdα,B n =2nπα ∫ l 0 ψ(α)sinnπl αdα(12)
周期函数的傅里叶展开
设以2l为周期的函数在[−l,l]上可以展开为三角函数:
f(x)=a 0 2 +∑ n=l ∞ (a n cosnπl x+b n sinnπl x)
则⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ a n =1l ∫ l −l f(x)cosnπl xdx,n=0,1,⋯b n =1l ∫ l −l f(x)sinnπl xdx,n=1,2,⋯
若f(x)为奇函数,则
f(x)=∑ n=1 ∞ b n sinnπxl ,b n =2l ∫ l 0 f(x)sinnπl xdx,n=1,2,⋯
若f(x)为偶函数,则
f(x)=a 0 2 +∑ n=1 ∞ a n cosnπxl ,a n =2l ∫ l 0 f(x)cosnπl xdx,n=0,1,⋯
8.1.3分析解答
1.存在性:
当φ∈c 2 ,ψ∈c l ,且满足边界条件(2),则(11),(12)给出的解存在.
2.物理意义:A n =N n cosδ n ,B n =N n sinδ n
即ω n =nπal ,u(x,t)=∑ n=l ∞ u n =∑ ∞ N n cos(ω n t−δ n )sinnπl x−−驻波叠加
振幅:N n sinnπl x
频率:ω n
初相位:δ n
波节:x m =mn l,m=0,1,⋯,n(共n+1个)
波腹:x k =2k−12n l,k=1,2,⋯,n(n个)
8.1.4小结
1.(1)−(3)式的解,由(11),(12)式给出;
⎧ ⎩ ⎨ ⎪ ⎪ u tt =a 2 u xx ,0<x<l(1)u| x=0 =0,u| x=l =0(2)u| t=0 =φ(x),u t | t=0 =ψ(x)(3)
u(x,t)=∑ n=1 ∞ (A n cosnπal t+B n sinnπal t)sinnπl x(11)
A n =2l ∫ l 0 φ(α)sinnπl αdα,B n =2nπa ∫ l 0 ψ(α)sinnπl αdα(12)
2.驻波法也适于求其它的齐次有界问题,又叫分离变量法。
3.分离变量法要领是,令u(x,y,z,⋯,t)=X(x)Y(y)Z(z)⋯T(t),从而将偏微分方程变成常微分方程求解。
4.分离变量法的解题步骤为:
①对齐次方程和齐次边界条件分离变量
②解常微分方程的本征值问题
③解其它变量的常微分方程
④叠加,用初始条件(或非齐次边界条件)定系数
8.1.5例题
考察在矩形薄板内稳定状态的温度分布,板的两对边绝热,而其余的两边一边温度保持为零,另一边的温度由f(x)规定.
解:其定解问题可表示为:
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ u xx +u yy =0,0<x<a,0<y<b(13)u x | x=0 =0u x | x=a =0 }(14)u| y=b =0(15)u| y=0 =f(x)(16)
1.分离变量:令u(x,y)=X(x)Y(y)
则(13)→{X ′′ −μX=0(17)Y ′′ +μY=0(18)
(14)→{X ′ (0)=0X ′ (a)=0 (19)
(15)→Y(b)=0(20)
2.解本征值问题(17),(19):
本征值:μ=−n 2 π 2 a 2 ,n=0,1,2,⋯
本征函数:X n (x)=A n cosnπa x
3.求解关于Y的方程:Y ′′ −n 2 π 2 a 2 Y=0(18) ′
Y n (y)=⎧ ⎩ ⎨ C 0 y+D 0 ,(n=0)C n coshnπa y+D n sinhnπa y=E n sinhnπa (y+F n ),(n≠0)
E n =D 2 n −C 2 n − − − − − − − √ ,F n =anπ th −1 C n D n
由Y n (b)=0有:C 0 =−D 0 b ,F n =−b(E n ≠0)
Y n (y)=⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ D 0 b−yb (n=0)E n shnπa (y−b)(n≠0)
4.叠加用非齐次边界条件定系数
u(x,y)=∑ n=0 ∞ X n (x)Y n (y)
=a 0 2 b−yb +∑ n=1 ∞ a n cosnπa xshnπa (y−b)
u| y=0 =f(x)(16)→f(x)=a 0 2 +∑ n=1 ∞ a n cosnπa xsh(−nπba )
记a ∗ n =2a ∫ a 0 f(x)cosnπxa dx
则u(x,y)=a ∗ 0 2 b−yb +∑ n=1 ∞ a ∗ n shnπa (b−y)shnπba cosnπa x
§8.2非齐次方程−−纯强迫震动
8.2.1定解问题
⎧ ⎩ ⎨ ⎪ ⎪ u tt =a 2 u xx +f(x,t),0<x<l,t>0(1)u| x=0 =0,u| x=l =0(2)u| t=0 =0,u t | t=0 =0(3)
8.2.2求解
思路1:⎧ ⎩ ⎨ ⎪ ⎪ u tt =a 2 u xx +f(x,t)u| t=0 =0u t | t=0 =0 →⎧ ⎩ ⎨ ⎪ ⎪ v tt −a 2 v xx =0v| t=τ =0v t | t=τ =f(x,τ)
u(x,t)=∫ t 0 v(x,t;τ)dτ
思路2:考虑二阶非齐次的常微分方程的求解:
对于y ′′ (x)+p(x)y ′ +Q(x)y=f(x)(A)
考虑齐次y ′′ (x)+p(x)y ′ +Q(x)y=0(B)
若(B)有通解:y g (x)=C 1 y 1 (x)+C 2 y 2 (x)
则由:常数变易法可令(A)有特解
y s (x)=C 1 (x)y 1 (x)+C 2 (x)y 2 (x)(C)
将(C)式代入(A)并补充条件:
C ′ 1 (x)y 1 (x)+C ′ 2 (x)y 2 (x)=0
则有:C ′ 1 (x)y ′ 1 (x)+C ′ 2 (x)y ′ 2 (x)=f(x)
于是可求得:C 1 (x)和C 2 (x),→y s (x)
1.对应的齐次问题:
{u tt =a 2 u xx u| x=0 =0,u| x=l =0
令u(x,t)=X(x)T(t)则可得:
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ X ′′ −μX=0μ=−(nπl ) 2 ,n=1,2,⋯X(0)=0X(l)=0X n (x)=C n sinnπxl
2.求对应的T n (t)方程的解
令u(x,t)=∑ n=1 ∞ T n (t)sinnπxl
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∑ n=1 ∞ [T ′′ n (t)+(anπl ) 2 T n (t)]sinnπxl =f(x,t)⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∑ n=1 ∞ T n (0)sinnπxl =0∑ n=1 ∞ T ′ n (0)sinnπxl =0 ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ T ′′ n (t)+(anπl ) 2 T n (t)=f n (t){T n (0)=0T ′ n (0)=0
f n (t)=2l ∫ l 0 f(α,t)sinnπαl dα
T n (t)=lnπa ∫ t 0 f n (τ)sinnπal (t−τ)dτ(4)
3.有界弦(杆)的纯强迫震动的解:
u(x,t)=∑ n=1 ∞ [lnπa ∫ t 0 f n (τ)sinnπal (t−τ)dτ]sinnπl x(5)
8.2.3小结
1.定解问题(1) (3)的解由(5)式给出。
2.本征函数法:以上求解非齐次方程的方法,显然也使用于求解带有其它齐次边界条件的各类方程。其中主要步骤为:
①用分离变量法求得对应的齐次问题的本征函数。
②将位置函数按求得的本征函数展开,其展开系数为另一变量的函数,代入非齐次方程和初始条件(或另一变量的边界条件),得到另一单元函数的非齐次常微分方程的定解问题。
③用常数变易法或拉氏变换法解非齐次常微分方程的定解问题,从而可求得原定解问题的解。此即本征函数法。
3.对于一般的两端固定的弦的强迫震动:
⎧ ⎩ ⎨ ⎪ ⎪ u tt =a 2 u xx +f(x,t)u| x=0 =0,u| x=l =0u| t=0 =φ(x),u t | t=0 =ψ(x)
令u=u i +u ii ,使:
§8.1:⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ u i tt =a 2 u i xx u i | x=0 =0,u i | x=l =0u i | t=0 =φ(x),u i t | t=0 =ψ(x) §8.2⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ u ii tt =a 2 u ii xx +f(x,t)u ii | x=0 =0,u ii | x=l =0u ii | t=0 =0,u ii t | t=0 =0
8.2.4例题
求解定解问题:⎧ ⎩ ⎨ ⎪ ⎪ u t −a 2 u xx =Asinωtu x | x=0 =0,u x | x=l =0u| t=0 =0
解:①对应的齐次方程的本征值问题为
{X ′′ −μX=0X ′ (0)=0,X ′ (l)=0 →μ=−n 2 π 2 l 2 ,
X n (x)=C n cosnπxl ,n=0,1,2,⋯
②u(x,t)=∑ n=0 ∞ T n (t)cosnπxl ,n=0,1,2,⋯
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∑ n=0 ∞ [T ′ n (t)+(anπl ) 2 T n (t)]cosnπxl =Asinωt∑ n=0 ∞ T n (0)cosnπxl =0
{T ′ 0 (t)=AsinωtT 0 (0)=0 ,n=0⎧ ⎩ ⎨ T ′ n (t)+(nπal ) 2 T n (t)=0T n (0)=0 ,n≠0
T 0 (t)=Aω (1−cosωt);T n (t)=0(n=1,2,3,⋯)
u(x,t)=Aω (1−cosωt)
§8.3非齐次边界条件的处理
8.3.1定解问题
⎧ ⎩ ⎨ ⎪ ⎪ u tt −a 2 u xx =0,0<x<l,t>0(1)u| x=0 =g(t),u| x=l =h(t)(2)u| t=0 =φ(x),u t | t=0 =ψ(x)(3)
8.3.2求解
1.思路:若令u(x,t)=X(x)T(t)
(2)式→{X(0)T(t)=g(t)X(l)T(t)=h(t) →{X(0)=g(t)/T(t)X(l)=h(t)/T(t)
2.求解:
(1)边界条件齐次化:
令u(x,t)=v(x,t)+w(x,t)(4)
使{w| x=0 =u| x=0 =g(t)(5)w| x=l =u| x=l =h(t)(6)
(2)确定辅助函数w(x,t)
令w(x,t)=A(t)x+B(t)
则w(x,t)=h(t)−g(t)l x+g(t)(7)
(3)求解v(x,t)的定解问题
(1)−(3)式→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ v tt −a 2 v xx =−(w tt −a 2 w xx )(8)v| x=0 =0,v| x=l =0(9)v| t=0 =φ(x)−w(x,0)v t | t=0 =ψ(x)−w t (x,0) }(10)
(4)定解问题(1)−(3)式的解
u(x,t)=v(x,t)+h(t)−g(t)l x+g(t)
8.3.3小结
1.以上介绍的方法也适用于带有其它非齐次边界条件的定解问题,其基本做法是:
①作变换:令u(x,t)=v(x,t)+w(x,t),
选择w(x,t)=A(t)x+B(t)
或w(x,t)=A(t)x 2 +B(t)x(当两端均为第2类非齐次边界条件时)
确定A(t),B(t)使关于v(x,t)的边界条件齐次化。
②解关于v(x,t)的定解问题,从而最后求得:u(x,t)=v(x,t)+w(x,t)
2.边界条件的齐次化,一般将导致方程的非齐次化
8.3.4例题
试研究一端固定,一端做周期运动sinωt的弦运动。
⎧ ⎩ ⎨ ⎪ ⎪ u tt −a 2 u xx =0,0<x<l(1)u(0,t)=0,u(l,t)=sinωt(2)u(x,0)=0,u t (x,0)=0,0≤x≤l(3)
①令u(x,t)=v(x,t)+w(x,t)(4)
选w(x,t)=sinωtl x+0=xl sinωt(5)
→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ v tt −a 2 v xx =ω 2 l xsinωt(6)v(0,t)=v(l,t)=0(7)v(x,0)=0,v t (x,0)=−ωl x(8)
②令v(x,t)=v i (x,t)+v ii (x,t)(9)
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ v i tt −a 2 v i xx =0v i (0,t)=v i (l,t)=0v i (x,0)=0,v i t (x,0)=−ωl x (10)⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ v ii tt −a 2 v ii xx =ω 2 l xsinωtv ii (0,t)=v ii (l,t)=0v ii (x,0)=v ii t (x,0)=0 (11)
③求解⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ v i tt −a 2 v i xx =0v i (0,t)=v i (l,t)=0v i (x,0)=0,v i t (x,0)=−ωl x (10)
v i (x,t)=∑ n=1 ∞ [A n cosnπal t+B n sinnπal t]sinnπl x
A n =0,B n =2nπa ∫ l 0 −ωl αsinnπαl dα=2ωl(−1) n (nπ) 2 a
v i (x,t)=2ωlπ 2 a ∑ n=1 ∞ (−1) n n 2 sinnπal tsinnπl x(12)
④令v ii (x,t)=∑ n=1 ∞ T n (t)sinnπl x
(11)→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ∑ n=1 ∞ [T ′′ n (t)+(anπ) 2 l 2 T n (t)]sinnπxl =ω 2 l xsinωt∑ n=1 ∞ T n (0)sinnπxl =0∑ n=1 ∞ T ′ n (0)sinnπxl =0
→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ T ′′ n (t)+a 2 n 2 π 2 l 2 T n (t)=f n (t)T n (0)=0T ′ n (0)=0 (13)
f n (t)=2l ∫ l 0 ω 2 l αsinωtsinnπαl dα=2ω 2 nπ sinωt(−1) n (14)
T n (t)=lnπa ∫ l 0 f n (τ)sinnπal (t−τ)dτ
ω n =nπal
T n (t)=lnπa ∫ l 0 2ω 2 (−1) n+1 nπ sinωτsinnπa(t−τ)l dτ
=ω 2 l(−1) n+1 a(nπ) 2 [sinω n t+sinωtω n +ω −sinω n t−sinωtω n −ω ](15)
v ii (x,t)=∑ n=1 ∞ T n (t)sinnπl xw(x,t)=xl sinωt(5)
v i (x,t)=2ωlπ 2 a ∑ n=1 ∞ (−1) n n 2 sinnπal tsinnπl x(12)
u(x,t)=v i (x,t)+v ii (x,t)+w(x,t)
§8.4正交曲线坐标系
8.4.1正交曲线坐标系
由三族互相正交的曲面而定义的坐标系。
1.柱坐标(ρ,φ,z):0≤ρ<∞,−∞<φ<∞,−∞<z<∞
⎧ ⎩ ⎨ ⎪ ⎪ x=ρcosφy=ρsinφz=z ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ r=x 2 +y 2 − − − − − − √ φ=tg −1 yx z=z
2.极坐标(ρ,φ):0≤ρ<∞,−∞<φ<∞
{ρcosφy=ρsinφ ⎧ ⎩ ⎨ ρ=x 2 +y 2 − − − − − − √ φ=tg −1 yx
3.球坐标系(r,θ,φ):0≤r<∞,0<θ<π,−∞<φ<∞
⎧ ⎩ ⎨ ⎪ ⎪ x=rsinθcosφy=rsinθcosφz=rcosθ ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ r=x 2 +y 2 +z 2 − − − − − − − − − − √ θ=tg −1 x 2 +y 2 − − − − − − √ z φ=tg −1 yx
8.4.2坐标系的选择
选择的坐标系使所研究问题的边界面和一个做表面重合。
边界 | 坐标 |
---|---|
长方形 | 直 |
球 | 球 |
圆柱 | 柱 |
圆锥 | 球 |
8.4.3正交曲线坐标系中的Δu
1.在柱坐标系中
∂u∂ρ =∂u∂x cosφ+∂u∂y sinφ(1)
∂ 2 u∂ρ 2 =∂ 2 u∂x 2 cos 2 φ+2∂ 2 u∂x∂y sinφcosφ+∂ 2 u∂y 2 sin 2 φ(2)
∂u∂φ =−∂u∂x ρsinφ+∂u∂y ρcosφ(3)
∂ 2 u∂φ 2 =∂ 2 u∂x 2 ρ 2 sin 2 φ−2∂ 2 u∂x∂y ρ 2 sinφcosφ+∂ 2 u∂y 2 ρ 2 cos 2 φ−(∂u∂x cosφ+∂u∂y sinφ)ρ(4)
Δu=1ρ ∂∂ρ (ρ∂u∂ρ )+1ρ 2 ∂ 2 u∂φ 2 +∂ 2 u∂z 2
2.在极坐标系中
Δu=1ρ ∂∂ρ (ρ∂u∂ρ )+1ρ 2 ∂ 2 u∂φ 2
3.在球坐标系中
Δu=1ρ ∂∂r (r 2 ∂u∂r )+1r 2 sinθ ∂∂θ (sinθ∂u∂θ )+1r 2 sin 2 θ ∂ 2 u∂φ 2
§8.5正交曲线坐标系中的分离变量
8.5.1Δu+λu=0(Δu=0)在物理学中的重要地位
令u(x,y,z;t)=T(t)v(x,y,z)
u tt =a 2 Δu→{T ′′ +a 2 λT=0Δv+λv=0
u t =DΔu→{T ′ +λDT=0Δv+λv=0
Δu=0
8.5.2柱坐标系中亥姆霍兹方程的分离变量
1.Δu+λu=0
1ρ (ρ∂u∂ρ )+1ρ 2 ∂ 2 u∂φ 2 +∂ 2 u∂z 2 +λu=0
令u(ρ,φ,z)=R(ρ)Φ(φ)Z(z)
→⎧ ⎩ ⎨ ⎪ ⎪ Z ′′ +μZ=0Φ ′′ +n 2 Φ=0ρ 2 R ′′ +ρR ′ +(k 2 ρ 2 −n 2 )R=0
本征值:μ,n 2 ,k 2
本征函数:Z,Φ,R
设λ−μ≥0x=kρ,y(x)=R(ρ),记λ−μ=k 2
x 2 y ′′ +xy ′ +(x 2 −n 2 )y=0n阶Bessel
2.Δu=0
Δu+λu=0
令u(ρ,φ,z)=R(ρ)Φ(φ)Z(z)
当λ=0时,Δu=0
→⎧ ⎩ ⎨ ⎪ ⎪ Z ′′ +μZ=0Φ ′′ +n 2 Φ=0ρ 2 R ′′ +ρR ′ +(k 2 ρ 2 −n 2 )R=0(−μ=k 2 )
x=kρ,y(x)=R(ρ)
x 2 y ′′ +xy ′ +(x 2 −n 2 )y=0
3.例:一个半径为a薄圆盘,上下两面绝热。若已知圆盘边温度,求圆盘上的稳定温度分布。
{Δu=0,ρ<au| ρ=a =f(φ)
(1)令u(ρ,φ)=R(ρ)Φ(φ)→{Φ ′′ +n 2 Φ=0ρ 2 R ′′ (ρ)+ρR ′ (ρ)−n 2 R=0
(2)解{Φ ′′ +n 2 Φ=0,n=0,1,2,⋯Φ(φ+2π)=Φ(φ)
Φ n (φ)=A ′ n cosnφ+B ′ n sinnφ
(3)解{ρ 2 R ′′ (ρ)+ρR ′ (ρ)−n 2 R=0R(ρ)| ρ=0 →有限 →R n (ρ)=C n ρ n
(4)u(ρ,φ)=∑ n=0 ∞ ρ n (A n cosnφ+B n sinnφ)
∑ n=0 ∞ (A n a n cosnφ+B n a n sinnφ)=f(φ)
α 0 =A 0 a 0 =12π ∫ π −π f(φ)dφ,α n =A n a n =1π ∫ π −π f(φ)cosφdφ
β n =B n a n =1π ∫ π −π f(φ)sinφdφ
u(ρ,φ)=∑ n=0 ∞ (ρa ) n (α n cosnφ+β n sinnφ)
8.5.3柱坐标系中亥姆霍兹方程的分离变量
1.Δu+λu=0
1r 2 ∂∂r (r 2 ∂u∂r )+1r 2 sinθ ∂∂θ (sinθ∂u∂θ )+1r 2 sin 2 θ ∂ 2 u∂φ 2 +λu=0(1)
令u(r,θ,φ)=R(r)y(θ,φ)
→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ r 2 d 2 Rdr 2 +2rdRdr +[k 2 r 2 −l(l+1)]R=0(2)1sinθ ∂∂θ (sinθ∂y∂θ )+1sin 2 θ ∂ 2 y∂φ 2 +l(l+1)y=0(3)
(3)→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ y(θ,φ)=Θ(θ)Φ(φ)Φ ′′ +m 2 Φ=0,m=0,1,2,⋯(4)1sinθ ddθ (sinθdΘdθ )+[l(l+1)−m 2 sin 2 θ ]Θ=0(5)
u(r,θ,φ)=R(r)Θ(θ)Φ(φ),(λ=k 2 )
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ r 2 d 2 Rdr 2 +2rdRdr +[k 2 r 2 −l(l+1)]R=0(2)Φ ′′ +m 2 Φ=0,m=0,1,2,⋯(4)1sinθ ddθ (sinθdΘθ )+[l(l+1)−m 2 sin 2 θ ]Θ=0(5)
本征值:l(l+1);m 2 ;l(l+1),m 2
本征函数:R(r),Φ,Θ
令:x=kr,y(x)/x √ =R(r)
(2)→x 2 y ′′ +xy ′ +[x 2 −(l+12 ) 2 ]y=0(2) ′
x=cosθ,y(x)=Θ(θ)
(5)→(1−x 2 )y ′′ −2xy ′ +[l(l+1)]−m 2 1−x 2 ]y=0(5) ′
(2),(2) ′ −球Bessel;(5),(5) ′ −缔合Legender
2.Δu=0
Δu+λu=0,令u(r,θ,φ)=R(r)Θ(θ)Φ(φ)
当λ=0时,Δu=0
→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ Φ ′′ +m 2 Φ=0,m=0,1,2,⋯(4)1sinθ ddθ (sinθdΘdθ )+[l(l+1)−m 2 sin 2 θ ]Θ=0(5)r 2 d 2 Rdr 2 +2rdRdr −l(l+1)R=0(6) −Euler方程
x=cosθ,y(x)=Θ(θ)
(5)→(1−x 2 )y ′′ −2xy ′ +[l(l+1)−m 2 1−x 2 ]y=0(5) ′ −缔合Legendre方程
8.5.4小结
1.在柱坐标系中
Δu=1ρ ∂∂ρ (ρ∂u∂ρ )+1ρ 2 ∂ 2 u∂φ 2 +∂ 2 u∂z 2
2.在极坐标系中
Δu=1ρ ∂∂ρ (ρ∂u∂ρ )+1ρ 2 ∂ 2 u∂φ 2
3.在球坐标系中
Δu=1r 2 ∂∂r (r 2 ∂u∂r )+1r 2 sinθ ∂∂θ (sinθ∂u∂θ )+1r 2 sin 2 θ ∂ 2 u∂φ 2
4.在柱坐标系中
令u(ρ,φ,z)=R(ρ)Φ(φ)Z(z)
Δu+λu=0Δu=0 →⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ Z ′′ +μZ=0Φ ′′ +n 2 Φ=00≤k 2 ={λ−μ−μ ρ 2 R ′′ +ρR ′ +(k 2 ρ 2 −n 2 )R=0
x=kρ,y(x)=R(ρ)
x 2 y ′′ +xy ′ +(x 2 −n 2 )y=0
本征值:μ,n,k 2
本征函数:Z,Φ,R
5.在球坐标系中
令u(r,θ,φ)=R(r)Θ(θ)Φ(φ)
Δu+λu=0→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ r 2 d 2 Rdr 2 +2rdRdr +[k 2 r 2 −l(l+1)]R=0Φ ′′ +m 2 Φ=0,m=0,1,2,⋯1sinθ ddθ (sinθdΘdθ )+[l(l+1)−m 2 sin 2 θ ]Θ=0
Δu=0→⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ Φ ′′ +m 2 Φ=0,m=0,1,2,⋯1sinθ ddθ (sinθdΘdθ )+[l(l+1)−m 2 sin 2 θ ]Θ=0r 2 d 2 Rdr 2 +2rdRdr −l(l+1)R=0
本征值:l(l+1),m 2 ,l(l+1),m 2
本征函数:R(r),Φ,Θ
(1−x 2 )y ′′ −2xy ′ +[l(l+1)−m 2 1−x 2 ]y=0(5) ′