《电动力学》第24讲 §4.4 分离变量法
电动力学二四(镜象法)
25
物理结果讨论: 物理结果讨论:
Q(Q0 − Q′) QQ′ 4πε0F = + 2 2 a (a − b) QQ Q R 2a − R 0 = 2 − 3 2 a a a −R
2
(
3 0
(
2
2 0 2 2 0
)
)
过渡到点 电荷相互 作用模型
R0 →0
吸引力, 吸引力, 趋于消失
26
QQ Q R 2a − R 0 4πε0F = 2 − 3 2 a a a −R
2
(
3 0
(
2
2 0 2 2 0
)
)
吸引力起主要作用 数值大于第一项) (数值大于第一项) 即使Q 即使Q和Q0同号 只要Q ,只要Q距球面足 够近, 够近,就受到导体 的吸引力。 的吸引力。
a→ R0
原因: 原因:虽然整个导 体的电荷与Q 体的电荷与Q同号 但在靠近Q ,但在靠近Q的球 面部分出现异号电 荷。从而相互吸引 起主要作用。 起主要作用。
可以看出,引入象电荷取代感应电荷, 可以看出,引入象电荷取代感应电荷,的确是 一种求解泊松方程的简洁方法。 一种求解泊松方程的简洁方法。
13
真空中有一半径为R 例2 真空中有一半径为R0的接 地导体球,距球心为a 地导体球,距球心为a(a>R0) 处有一点电荷Q,求空间各点的电 处有一点电荷Q 势(如图)。 如图)。
8
解
电荷: 电荷:一个点电荷 界面: 界面:接地无穷大导体 区域:上半空间(下半空间电势为零) 区域:上半空间(下半空间电势为零)
已知界面电势为零, 已知界面电势为零,满足唯一性定理 的要求,可以确定电势。 的要求,可以确定电势。
9
分离变量法常数变易法、行波法和积分变换法达朗贝尔
分离变量法常数变易法、行波法和积分变换法达朗贝尔设w=u+iV及z=x+iy分别是两个复平面上的点,复函数w=f(z)确定了这两个复平面之间的一个映射,当w=f(z))是一个目数不为零的解析函数时,所对应的映射称为保角映射。
保角映射这种映射必定是一对一的,且具有:(l)伸缩率的不变性,即在某一点Z0上沿不同的方向的曲线微元ds与映射后所得的象ds′的比值都是f′(z0);(2)旋转角的不变性并且保持角的定向,即若把z平面与w平面迭放在一起,且使ZO与W0=f(z0)重合,则过Z0的任一条曲线C到它的象C′的转角为定值。
如果X轴与U轴及y轴与V轴方向相同,这个转角就是Argf'(z0),因此交手Z0的任意两条曲线C1,C2的夹角与它们的象C1,C2的夹角相等且转向不变。
保角变换方法(conformaltransformationmethod)保角变换是利用复变量解析函数实部和虚部都满足拉普拉斯(Laplace)方程的特点,及通过复平面变换以简化求解二维拉普拉斯方程边值问题的一种方法。
由于在没有电荷分布的空间中静电势满足拉普拉斯方程,故此法可用来求解二维的静电势问题。
通过一适当的解析复变函数f(z),将复变数平面z=x+iy变换成另一复变数平面z′=f(z)=x′+iy′或z=g(z′)将z平面上位形复杂的边值问题,变换至z′平面上位形简单的相应边值问题,以便容易求出静电势的解φ′(x′,y′)。
由此在z′平面中构成解析的复变函数W′(z′)=φ′+i Ψ′。
最后再由z′平面换回z平面W(z)=W′(f(z))=φ(x,y)+iΨ(x,y),从而得到欲求的二维拉普拉斯方程边值问题的解。
由于通过解析函数变换时,分别在二复平面中任意二曲线元之间的夹角不变,故此种变换称为保角变换。
保角映射英文术语名:conformaltransformation【保角映射的定义】设f(z)是区域D到G的双射(既是单射又是满射),且在D内的每一点都具有保角性质,则称f(z)是区域D到G的保角映射,也称为保角变换或者共形映射。
分离变量法的精神和解题要领PPT教案
答案:
ut a2uxx (0 x l,t 0)
u(0,t) 0
ux
(l,t)
q k
u(x,0)
x(l 2
x)
第2页/共17页
解题方法:
1.建立方程, 2.定解条件:边界条件,初始条件 3.定解问题
建立方程——解题思路: • 由能量守恒定律
特征值问题同热导相同
第7页/共17页
u(x,t) uk (x,t) X k (x)Tk (t)
k 1
k 1
ka
ka k
k 1
(Ak cos
l
t Bk sin
l
t ) sin
l
x
Ak
2 l
l
( ) sin
0
n l
d
形式不变
Bk
2 ka
l
(
)
sin
n
0
l
d
第8页/共17页
有界杆上的热传导
u(x,t) X (x)T(t)
形式不变
特征值问题同振动方程相同
第9页/共17页
第10页/共17页
u(utx,0a) 2x2u(2x,),
x(0,l),t 0 x [0, l]
u(0,t) ux(l,t) 0, t 0
本征值 和本征
函数
n
n
1 2
l
2
,
X n (x)
sin
n
l
流沿x轴正向,强度为 c ρut = k uxx
u
q(x,t),温度分布为
ut = a2 uxx
u(x,t),则
第3页/共17页
第24讲分离变量法第4章介质中的电动力学4§4拉普拉斯方程分离变量法
第24讲 分离变量法 第4章 介质中的电动力学(4) §4.4 拉普拉斯方程 分离变量法以上两节给出静电问题的一般公式,并说明静电学的基本问题式求解满足给定边界条件的泊松方程的解。
只有在界面形状是比较简单的几何曲面时,这类问题的解才能以解析形式给出,而且视具体情况不同而有不同的解法。
在许多实际问题中,静电场是由带电导体决定的。
例如电容器内部的电场是由作为电极的两个导体板上所带电荷决定的;又如电子光学系统的静电透镜内部,电场是由于分布于电极上的自由电荷决定的。
这些问题的特点是自由电荷只出现在一些导体的表面上,在空间中没有其它自由电荷分布。
因此,如果我们选择这些导体表面作为区域V 的边界,则在V 内部自由电荷密度 ρ = 0 ,因而泊松方程化为比较简单的拉普拉斯(Laplace )方程20ϕ∇= (4.4---1) 产生这电场的电荷都分布于区域V 的边界上,它们的作用通过边界条件反映出来。
因此,这类问题的解法是求拉普拉斯方程的满足边界条件的解。
(4.4---1)式的通解可以用分离变量法求出。
先根据界面形状选择适当的坐标系,然后在该坐标系中由分离变量法解拉普拉斯方程。
最常用的坐标系有球坐标系和柱坐标系。
这里我们写出用球坐标系得出的通解形式(见附录Ⅱ)。
球坐标用(R ,θ,φ)表示,R 为半径,θ为极角,φ为方位角。
拉氏方程在球坐标系中的通解为1.(,,)()(cos )cos n mnm nm n n n mb R a R P m R ϕθφθφ+=+∑ 1,()(cos )sin n mnm nm n n n md c R P m Rθφ+++∑ (4.4---2) 式中 a n m ,b n m ,c n m 和 d n m 为任意常数,在具体问题中有边界条件定出。
P m n (cos θ) 为缔和勒让德(Legendre )函数。
若该问题中具有对称轴,取此轴为极轴,则电势φ不依赖于方位角φ,这情形下通解为 1()(cos ),n nn n n nb a R P Rϕθ+=+∑ (4.4---3) P n (cos θ)为勒让德函数,a n 和b n 由边界条件确定。
格里菲斯 电动力学 分离变量法
格里菲斯电动力学分离变量法下载提示:该文档是本店铺精心编制而成的,希望大家下载后,能够帮助大家解决实际问题。
文档下载后可定制修改,请根据实际需要进行调整和使用,谢谢!本店铺为大家提供各种类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you! In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!一、简介电动力学是物理学中的一个重要分支,研究电荷和电磁场之间的相互作用。
数学物理方程经典教案分离变量法
数学物理方程经典教案分离变量法分离变量法是数学物理中常用的求解偏微分方程的方法之一、它适用于一类特殊的二元函数方程,即能够通过变量分离的方式将方程化为两个只依赖于一个变量的常微分方程。
为了更好地理解和应用分离变量法,下面将从理论和实践两个方面进行介绍和解释。
一、理论介绍1.分离变量法的基本思想:对于具有特定形式的二元函数方程,我们可以通过合适的变量变换,将方程化为两个只依赖于一个变量的常微分方程。
从而达到求解方程的目的。
2.分离变量法的基本步骤:(1)假设原方程的解具有特定的形式,例如f(x,y)=X(x)Y(y)。
(2)将f(x,y)代入原方程,化简得到两个只依赖于一个变量的常微分方程。
一个关于X(x),一个关于Y(y)。
(3)解决两个常微分方程,得到X(x)和Y(y)的解。
(4)组合求解得到原方程的解。
3.分离变量法的适用范围:分离变量法适用于线性偏微分方程和一些非线性偏微分方程的特殊情况。
它在物理学、工程学和应用数学中有广泛的应用。
二、实践应用为了更好地说明分离变量法的应用,我们以热传导方程为例进行介绍。
热传导方程是一个描述物质内部热传导过程的重要方程,在热力学、材料科学等领域有广泛应用。
其方程形式为:∂u/∂t=α(∂²u/∂x²+∂²u/∂y²+∂²u/∂z²)其中u(x,y,z,t)表示温度分布随时间的变化,α为热扩散系数。
分离变量法的具体求解步骤如下:1.假设温度分布函数可以表示为u(x,y,z,t)=X(x)Y(y)Z(z)T(t)。
2.将u(x,y,z,t)代入热传导方程中,得到四个只依赖于一个变量的常微分方程:X''(x)/X(x)=Y''(y)/Y(y)=Z''(z)/Z(z)=T'(t)/αT(t)。
3.解决这四个常微分方程,得到X(x)、Y(y)、Z(z)和T(t)的解。
分离变量法在静电场问题中的应用
2 分 离变 量法 求解 的方程 及通 解
在许 多实际问题中, 静电场是由带电导体决定的。例如: 电容器 内部的电场是由作为电极的两个导 体板上所带电荷决定的。而这些问题 的特点是 : 空间中有 自由电荷分布 , 自由电荷只出现在一些导体 而
的表面上 , 在空间中没有其它 自由电荷分布 。因此 , J 如果选择这些导体表面作为区域 V的边界 , 则在
再次 , 从实 际情 况 出发 , 依次 找 出边 界条 件和 边值关 系来确 定通 解 中的待定 系数 。
概括起来 , 大致有 以下几种类型的边界条件 : J ( )两种绝缘介质界面上, 1 边值关系为
‘l P ‘2 P
ap ‘1
8
ap ‘2
应用这条件可以把界面两边的电势衔接起来 。
V内部 自由电荷密度 p 0 因而泊松方程 ‘=一 : , p 旦化为比较简单的拉普拉斯方程 :
‘= P 0 u () 1
上 式 称为拉 普拉 斯方程 ( 氏方 程 ) 拉 。产 生 电场的 电荷 都分 布 于 区域 V 的边 界上 , 们 的作 用 通过 它
边界条件反应出来。因此 , 这类问题的解法是求拉普拉斯方程的满足边界条件的解, 1 式 的通解可以 ()
9 3
昌吉学 院学 报
21 0 1年第 4期
‘=∑( “ p aR +
) C5 ) P (00
() 3
P(o0 为勒让德函数, 和 b 是任意常数 , c ) s a 由边界条件或边值关系确定。所以, 在每一个没有电
荷分布的区域 内, 满足拉普拉斯方程 , ‘ p 其通解 已有( ) 3 式给出, 2 或( ) 剩下的问题就是由具体 的边界条
电动力学第24讲44分离变量法
1
4 0
PgR R3
P
4 0 R 2
cos
山东大学物理学院 宗福建
25
§2.4 分离变量法
n=2时,
n
(an Rn
bn R n 1
)
Pn
(cos
)
(a2R 2
b2 R3
)P2 (cos
)=
1 2
(a2R 2
b2 R3
)(3cos2 (x)
1)
电四极子的势(参考教材第66-69页)
山东大学物理学院 宗福建
19
§2.4 分离变量法
拉氏方程在球坐标系中的通解为
( R,
,)
n.m
(anm R n
bnm Rn1
)Pnm
(cos
)
cos
m
n,m
(cnm Rn
dnm Rn1
)Pnm
(cos
) sin
m
式中 a n m ,b n m ,c n m 和 d n m 为任意常数,在 具体问题中由边界条件定出。Pnm(cosθ) 为缔和
山东大学物理学院 宗福建
7
上一讲习题解答
解:无限空间时
1 Q Q 4 r2 E j E I 4 r2 j 4 r
I 4 r2 j 4 r2 E Q 1 I
4 r
I (
1
4 (x x0 )2 ( y y0 )2 (z z0 )2
(2)
1
12 0
D33
1 R3
分离变量法与有限差分法在电场的应用
结论
使用有限差分法求解电磁场边值问题是 可行的,只要将网格取得足够小,是可 以将离散的点看成连续的。
但节点划分越多,迭代次数就越多,计 算量就越大,所以计算时应根据实场的应用
问题引出
y
100
b
50
a x
有一长方形的导体 槽,a = 20,b = 5, 设槽的长度为无限 长,槽上有两块与 槽绝缘的盖板,电 位分别为100V和50V, 其它板电位为0V, 求槽内的电位分布。
有限差分法基本思想
用有限差分方程组近似原微分方程和定解条件
把连续的定解区域用有限个离散点构
成的网格来代替; 把连续定解区域上的连续变量的函数 用在网格上定义的离散变量函数来近 似; 把原方程和定解条件中的微商用差商 来近似,积分用积分和来近似。
有限差分法计算机仿真
x=21; y=6; v1=ones(y,x); v1(y,:)=ones(1,x)*100; v1(1,:)=ones(1,x)*50; fori=1:y; v1(i,1)=0; v1(i,x)=0; end max=1; t=0; v2=v1; n=0; while(max>1e-6) n=n+1; max=0;
(m为奇数)
分离变量法计算解析解
(n为奇数)
分离变量法计算解析解
通过(2)得
分离变量法计算解析解
而
(s为奇数)
分离变量法计算解析解
x=0:1:20; y=0:0.25:5; [X,Y]= meshgrid(x,y); z=ones(21); for n=1:2:110; j=sin(n*pi.*X/a); c=(200*a)/(n*pi); d=1exp(2*n*b*pi/a); f=exp(n*pi.*Y/a);
电动力学二三分离变量法-文档资料
2R 1R R R
2
3
Q 1 2 2 2 R d R d R R 0 R R R R 3 2
8
将通解代入边界条件
0 2R 1R R
1
a 0
d c 0 R1
2R 1R R R
比较P1的系数得
b 1 E R c 0 0 1R 0 2 R 0
2 b 1 E c 0 3 1 R 0 0
可解出
3 0 b E R , 1 00 2 0
3 0 c E 1 0 2 0
c 0 , 其他Pn项的系数可解出为 b n n
其中
1 R 3 Q Q 1 1 1 1 R R R 1 2 3
利用这些 值得电势 的解 导体球上 的感应电 荷为
QQ 1 , 1 4 0R 1 1 Q 1 2 . 4 0 R R 1
(R R 3) (R 2 R R 1)
球内区域的电势
n d n c R n P cos 2 n 1 n R n
12
边界条件:
(1)无穷远处,
因而
1 0 01
E R cos E R P co
a E , 1 0 a 0 ( n 1 ) n
(2)R=0处,2为有限值,因此 dn 0 (3)在介质球面上,有 1 2
15
例3 半径为R0的导体球置于均 匀外电场E0中,求电势和导体上 的电荷面密度。
16
解
用导体表面边 界条件,照上 例方法可解出 导体球外电势 导体面上 电荷面密 度为
分离变量法 电动力学
四.应用举例
1、两无限大平行导体板,相距为l,两板间电势差
为V (与 x,y,z 无关),一板接地,求两板间的电势
和电场。
请大家先写出答案
解:(1)边界为平面,故
Z
应选直角坐标系
V
下板 S1 0 ,设为参考点
l
O
(2)定性分析:
y
由对称性知 与 x, y 无关
x
(z)
机动 目录 上页 下页 返回 结束
X ( x)
c eh1x 1
c eh1x 2
机动 目录 上页 下页 返回 结束
(2)若 (x, y)
d2X dx2
k2X
d
2Y
dy2
k 2Y
若 k2 0
X ( x) Aekx Bekx Y ( y) C sin ky Dcos ky
注意:在(1)和(2)两种情况中若考虑了某些边
界条件,k1,k2,k3,k 将与某些正整数有关,它们可
取1,2,3,… ,只有对它们取和后才得到通解。
(3)若 ( x)
d 2 dx 2
0,
Ax B
机动 目录 上页 下页 返回 结束
2.柱坐标
讨论:
2
1 r
r
(r
) r
1 r2
2 2
2 z 2
0
(1) (r, ), 令 (r,) f (r)g()
2.一对接地半无限大平板,相距为 b ,左端有一极
板电势为 V(常数),求两平行板之间的电势。
解:(1)边界为平面,
选直角坐标系;上、下两
y
平板接地,取为参考点;
且
0
x
V
(2) z 轴平行于平板,
偏微分(3)分离变量法PPT课件
(
)
sin
ka
L
(t
)sin
k
L
x
d
其中
Bk
(
)
2
ka
L f ( , ) sin k d
0
L
(2.10)
2021/3/22
21
分离变量法:
令
k
u(x, t) Ck(t)sin
k 1
L
x
(2.11)
是混合问题的解。 显见上述函数满足(2.2)。
(2.1)
Ck(t)sin k 1
k
L
x k
从而定出 X (x) 所适合的常微分方程齐次边值问题,以及
T (t) 适合的常微分方程。
本征
求解该常微分方程齐次边值问题,
第二步 求出全部本征值和本征函数,并求
值问 题
出相应的 T (t) 的表达式。
第三步 将所有变量分离形式的特解叠加起来,并
利用初始条件定出所有待定系数。
2021/3/22
12
物理意义
正弦展开的Fourier级数的系数,即
Ak
2 L
L
( ) sin
k
d
0
L
(k 1,2,)
(1.17)
Bk
2
ka
L
( ) sin
k
d
0
L
(1.18)
这样,我们就给出了混合问题(1.1)-(1.4)的形 式解(1.14),其中系数由公式(1.17)和(1.18)给 出。
2021/3/22
10
sin x, sin 2 x, sin k x, 是[0, L]上的正交函数列
[
2024版《电动力学》ppt课件
电势分布及等势面描绘方法
电势定义
单位正电荷在电场中某点所具有的电势能。
电势差与电势分布
描述电场中两点间电势的差值,电势分布可通过求解泊松方程或 拉普拉斯方程得到。
等势面描绘
电势相等的点构成的曲面,其描绘方法包括解析法、图解法等。
2024/1/24
10
导体在静电场中特性研究
导体静电平衡条件
导体内部电场强度为零,电荷只分布在导体表面。
物理意义
揭示了电磁现象的基本规律,是电磁学的基础理 论。
方程组包括
高斯定律、高斯磁定律、麦克斯韦-安培定律和法 拉第感应定律。
2024/1/24
5
电磁波传播特性及波动方程
2024/1/24
电磁波
01
电场和磁场相互激发并在空间中传播形成的波动现象。
传播特性
02
电磁波在真空中以光速传播,具有能量和动量。
铁磁材料在恒定磁场中表现出非线性、磁饱和、磁滞等特性。
2024/1/24
03
应用举例
利用铁磁材料的特性制作电感器、变压器、电机等电气设备,以及用于
磁记录、磁放大等领域。
16
恒定磁场能量储存与转换
2024/1/24
恒定磁场能量密度 恒定磁场中储存的能量与磁场强度的平方成正比,能量密 度w=(1/2)BH。
26
无线通信系统基本原理简介
无线通信系统组成
包括发射机、信道、接收机等部分,实现信息 的传输和接收。
2024/1/24
无线通信基本原理
利用电磁波作为信息载体,通过调制将信息加载到载 波上,经过信道传输后,在接收端进行解调还原出原 始信息。
无线通信关键技术
包括调制与解调、信道编码与解码、多址接入、 抗干扰等技术,保证通信系统的可靠性和有效 性。
新版数学物理方程学习指导书第4章 分离变量法
第4章 分离变量法物理学、力学和工程技术等方面的许多问题都可归结为偏微分方程的定解问题,上一章我们已初步看到怎样把具体的物理问题表达为定解问题.下面一个重要任务是怎样去解决这些定解问题,也就是说在已经列出的方程与定解条件之后,怎样去求既满足方程又满足定解条件的解.从微积分学得知,在计算诸如多元函数的微分及重积分时总是把它们转化为单元函数的相应问题来解决.与此类似,求解偏微分方程的定解问题也是要设法把它们转化为常微分方程问题,分离变量法就是常用的一种转化手法.本章我们将通过实例来说明分离变量法的步骤和实质.在4.2我们讨论了如何处理第三类齐次边界条件(当然也包括第二类边界条件).在4.3说明如何在极坐标系下使用分离变量法.在4.4及4.5我们讨论了如何处理非齐次方程及非齐次边界条件的问题,本章的最后还安排了两个较为综合性的例子作为总结.4.1 有界弦的自由振动为了使读者了解什么是分离变量法以及使用分离变量法应该具备什么条件,我们选取两端固定的弦的自由振动问题为例,通过具体地求解逐步回答这些问题.根据第3章所得的结论,讨论两端固定的弦的自由振动,就归结为求解下列定解问题22222000,0,0; (4.1)0,0;(4.2)(),().(4.3)x x l t t uu a x l t t x u u u u x x t ϕψ====⎧∂∂⎪=<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩这个定解问题的特点是:偏微分方程是线性齐次的,边界条件也是齐次的,求解这样的问题,可以运用叠加原理.我们知道.在求解常系数线性齐次常微分方程的初值问题时,是先求出足够多个特解(它们能构成通解),再利用叠加原理作这些特解的线性组合,使满足初始条件.这就启发我们,要解问题(4.1),(4.2),(4.3),先寻求齐次方程(4.1)的满足齐次边界条件(4.2)的足够多个具有简单形式(变量被分离的形式)的特解,再利用它们作线性组合使满足初始条件(4.3).现在我们试求方程(4.1)的变量分离形式(,)()()u x t X x T t =的非零解,并要求它满足齐次边界条件(4.2),式中(),()X x T t 分别表示仅与x 有关及仅与t 有关的待定函数. 由(,)()()u x t X x T t =得2222''()(),()''(),u uX x T t X x T t x t∂∂==∂∂代入方程(4.1)得2()()()()X x T t a X x T t ''''=或2()()()()X x T t X x a T t ''''=这个式子左端仅是x 的函数,右端仅是t 的函数,一般情况下二者不可能相等,只有当它们均为常数时才能相等.令此常数为λ,则有2()()()()X x T t X x a T t λ''''==.这样我们得到两个常微分方程:2()()0,T t a T t λ''-= (4.4)()()0.X x X x λ''-= (4.5)再利用边界条件(4.2),由于(,)()()u x t X x T t =),故有(0)()0,()()0.X T t X l T t ==但()0T t ≠,因为如果()0T t ≡,则(,)0u x t ≡,这种解显然不是我们所要求的,所以(0)0,()0.X X l == (4.6)因此,要求方程(4.1)满足条件(4.2)的分离变量形式的解,就先要从方程''()()0,(0)()0X x X x X X l λ-=⎧⎨==⎩中解出()X x .现在我们就来求非零解()X x ,但要求出()X x 并不是一个简单的问题,因为方程(4.5)中含有一个待定常数λ,所以我们的任务既要确定λ取何值时方程(4.5)才有满足条件(4.6)的非零解,又要求出这个非零数()X x .这种常微分方程问题称为固有值问题,λ称为特征值(固有值,本征值),函数()X x 称为特征函数(固有函数,本征函数).下面根据第1章所介绍的方法,我们对λ分三种情况来讨论.1°λ>0,此时方程(4.5)的通解为().X x Be =+由条件(4.6)得0A B +=,0.Be +=解出,A B 得0A B ==,即()0X x ≡,不符合非零解的要求,因此λ不能大于零.2°设λ=0,此时方程(4.5)的通解为()X x Ax B =+,由条件(4.6)还是得0A B ==,所以λ也不能等于零.3°设λ<0,并令ββλ,2-=为非零实数.此时方程(4.5)的通解为()cos sin ,X x A x B x ββ=+由条件(4.6)得0,A = Bsin 0.l β=由于B 不能为零(否则()0X x ≡),所以sin 0,l β=即),,3,2,1(. ==n ln πβ(n 为负整数可以不必考虑,因为例如21,sin n B x l π-=-实际上还是2sin B x lπ'的形式)从而222,n lπλ=- (4.7)这样,我们就求出了一系列固有值及相应的固有函数:222.(1,2,3,),n n n lπλ=-=()sin(1,2,3,),n n n X x B x n lπ== (4.8)限定了λ的值后,现在再来求函数()T t ,以(4.7)式中的λ值代入方程(4.4)中得2222()()0,n a n T t T t lπ''+= 显然,其通解为''()cossin (1,2,3,).n n n n at n at T t C D n l lππ=+= (4.9)于是由(4.8),(4.9)得到满足方程(4.1)及边界条件(4.2)的一组变量被分离的特解(,)cos sin sin(1,2,3,),n n n n at n at n x u x t C D n l l l πππ⎛⎫=+= ⎪⎝⎭(4.10)其中,n n n n n n C B C D B D ''==是任意常数,至此,我们的第一步工作已经完成了,求出了既满足方程(4.1)又满足边界条件(4.2)的无穷多个特解.为了求原定解问题的解,还需要满足条件(4.3).由(4.10)式所确定的一组函数虽然已经满足方程(4.1)及条件(4.2),但不一定满足初始条件(4.3).为了求出原问题的解,首先我们将(4.10)中所有函数(,)n u x t 叠加起来1(,)(,)n n u x t u x t ∞==∑1c o s s i n s i n ,n nn n a n a n C t D x l l lπππ∞=⎛⎫=+ ⎪⎝⎭∑ (4.11)如果(4.11)右端的无穷级数是收敛的,而且关于,x t 都能逐项微分两次,则它的和(,)u x t 也满足方程(4.1)和条件(4.2)(参考习题三第6题).现在我们要适当选择,n n C D ,使函数(,)u x t 也满足初始条件(4.3),为此必须有01(,)(,0)sin(),n t n n u x t u x C x x lπϕ∞=====∑ 10sin (),n n t u n a n D x x t l lππψ∞==∂==∂∑ 因为(),()x x ϕψ是定义在[0,]l 上的函数,所在只要选取n C 为()x ϕ的傅氏正弦级数展开式的系数,n n aD lπ为()x ψ的傅氏正弦级数展开式的系数,也就是 002()sin ,2()sin .l n l nn C x xdx l ln D x xdx n a l πϕπψπ⎧=⎪⎪⎨⎪=⎪⎩⎰⎰ (4.12) 初始条件(4.3)就能满足,以(4.12)所确定的,n n C D 代入(4.11)式,即得原定解解问题的解.当然,如上所述,要使(4.11)式所确定的函数u(x,t)确定是问题(4.1),(4.2),(2.3)的解,除了其中的系数,n n C D 必须由(4.12)确定以外,还要求只要对函数()x ϕ及()x ψ加一些条件就能满足,可以证明(参阅复旦大学数学系编《数学物理方程》第二章§1),如果()x ϕ三次连续可微,)(x ψ二次连续可微,且(0)()(0)()(0)()0l l l ϕϕϕϕψψ''''======,则问题(4.1),(4.2),(4.3)的解存在.并且这个解可以用(4.11)给出,其中,n n C D 由(4.12)式确定*).从上面的运算过程可以看出,用分离变量法求解定解问题的关键步骤是确定固有函数与运用叠加原理,这些运算之所以能够进行,就是因为偏微分方程与边界条件都是齐次的,这一点希望读者一定要注意.例1 设有一根长为10个单位的弦,两端固定,初速为零,初位移为1000)10()(x x x -=ϕ,求弦作微小横向振动时的位移.解 设位移函数为(,)u x t ,它是下列定解解问题2222201000,010,0;0,0;(10),01000x x l t uu a x t t x u u x x u u t ====⎧∂∂⎪=<<>∂∂⎪⎪==⎨⎪-∂⎪==⎪∂⎩的解,这时10l =,并给定210000a =(这个数字与弦的材料、张力有关).显然,这个问题的傅氏级数形式解可由(4.11)给出,其系数按(4.12)式为10033330,1(10)sin 5000102(1cos )50,4,5n n D n C x x xdx n n n n n ππππ==-=-⎧⎪=⎨⎪⎩⎰当为偶数 当为奇数 因此,所求的解为3341(21)(,)sin cos10(21).5(21)10n n u x t x n t n πππ∞=+=++∑为了加深理解,下面我们扼要地分析一下级数形式解(4.11)的物理意义,先分析一下级数中每一项(,)cos sin sin n n n a n a n u x t C t D t x l l l πππ⎛⎫=+ ⎪⎝⎭*)这里所讲的结论经适当修改即可用于下面几节将要讨论的定解问题,所以,本书中凡是用分离变量法求解的定解问题都假定它的定解条件具备一定的条件,保证定解问题的解可以表示成级数的形式,或者说可以运用叠加原理.的物理意义,分析的方法是:先固定时间t ,看看在任一指定时刻波是什么形状;再固定弦上一点,看看该点的振动规律.把括号内的式子改变一下形式,可得(,)cos()sin,n n n n u x t A t x lπωθ=-其中,.n n n n nD n aA arctg l C πωθ=== 当时间t 取定值0t 时,得(,)sin,n n u x t A x lπ'= 其中0cos()n n n n A A t ωθ'=-是一个定值,这表示在任一时刻,波0(,)n u x t 的形状都是一些正弦曲线,只是它的振幅随着时间的改变而改变.当弦上点的横坐标x 取定值0x 时,得0(,)cos(),n n n n u x t B t ωθ=-其中0sinn n n B A x l π=是一个定值.这说明弦上以0x 为横坐标的点作简谐振动,其振幅为n B ,角频率为n ω,初位相为n θ.若x 取另外一个定值时,情况也一样,只是振幅n B 不同罢了,所以(,)n u x t 表示这样一个振动波:在考察的弦上各点以同样的角频率n ω作简谐振动,各点处的初位相也相同,而各点的振幅则随点的位置改变而改变;此振动波在任一时刻的图形是一正弦曲线.这种振动波还有一个特点,即在[0,]l 范围内还有1n +个点(包括两个端点)永远保持不动,这是因为在(0,1,2,,)m mlx m n n==那些点上,sinsin 0mn x m l ππ==的缘故,这些点在物理上称为节点.这就说明(,)n u x t 的振动是在[0,]l 上的分段振动,其中有1n +个节点,人们把这种包含节点的振动波叫做驻波.另外驻波还在n 个点处振幅达到最大值(读者可自己讨论),这种使振幅达到最大值的点叫做波腹.图4-1画出了在某一时刻1,2,3n =的驻波形状.综合上述,可知(,),(1,2,3,)n u x t n =是一系列驻波,它们的频率、位相与振幅图4-1都随n 不同而不同,因此我们可以说,一维波动方程用分离变量法解出的结果(,)u x t )是由一系列驻波叠加而成的,而每一个驻波的波形由固有函数确定,它的频率由固有值确定.这完全符合实际情况,因为人们在考察弦的振动时,就发现许多驻波,它们的叠加又可以构成各种各样的波形,因此很自然地会想到用驻波的叠加表示弦振动方程的解.这就是分离变量法的物理背景.所以分离变量法也称为驻波法.4.2 有限长杆上的热传导设有一均匀细杆,长为l ,两端点的坐标为0x =与x l =,杆的侧面是绝热的,且在端点0x =处温度是零度,而在另一端x l =处杆的热量自由发散到周围温度是零度的介质中去(参考第3章1.2中第三类边界条件),已知初始温度分布为().x ϕ求杆上的温度变化规律,也就是要考虑下列定解问题:222,0,0; 2.13(,)(0,)0,(,)0;(2.14)(,0)().(2.15)u u a x l t t x u l t u t hu l t x u x x ϕ⎧∂∂=<<>⎪∂∂⎪∂⎪=+=⎨∂⎪=⎪⎪⎩()我们仍用分离变量法来解这个问题,首先求出满足边界条件而且是变量被分离形式的特解.设(,)()()u x t X x T t =,2()().()()T t X x a T t X x '''=上式左端不含有x ,右端不含有t ,所以只有当两端均为常数时才可能相等.令此常数为2β-,(读者可以从方程()(),X x X x λ''=结合边界条件按λ取值的三种不同情况像4.1那样讨论后得出),则有22()(),()()T t X x a T t X x β'''==- 从而得到两个线性常数微分方程222()()0,()()0.T t a T t X x X x ββ'+=''+= (4.16)解后一个方程得()cos sin ,X x A x B x ββ=+由边界条件(4.14)可知(0)0,'()()0.X X l hX l =+=从(0)0X =得 0A =,从()()0.X l hX l '+=得cos sin 0.l h l βββ+= (4.17)为了求出β,方程(4.17)可改写成tg a γγ=, (4.18)其中1,.l a hlγβ==-方程(4.18)的根可以看作是曲线1y tg γ=与直线2y a γ=交点的横坐标(图4-1),显然它们的交点有无穷多个,于是方程(4.18)有无穷多个根,由这些根可以确定出固有值β.设方程(4.18)的无穷多个正根(不取负根是因为负根与正根只差一个符号(图4-2),再根据4.1中所述的同样理由)为123,,,,,n γγγγ于是得到无穷多个固有值图4-21212,,,,nn lllγγγβββ===及相应的固有函数()sin .n n n X x B x β= (4.19)再由(4.16)中第一个方程解得22().n a tn n T t A eβ-= (4.20)由(4.19),(4.20)两式,我们得到方程(4.13)满足边界条件(4.14)的一组特解22(,)()()sin (1,2,3,),n a t n n n n n u x t X x T t C e x n ββ-=== (4.21)其中 .n n n C A B =由于方程(4.13)与边界条件(4.14)都是齐次的,所以2211(,)(,)sin n a t n n n n n u x t u x t C e x ββ∞∞-====∑∑ (4.22)仍满足方程与边界条件,最后考虑(,)u x t 是否能满足初始条件(4.15),从(4.22)式得1(,0)sin .n n n u x C x β∞==∑现在希望它等于已知函数()x ϕ,那么首先要问[0,]l 上定义的函数()x ϕ是否能展开为1sin nn n Cx β∞=∑级数形式,其次要问系数n C 如何确定,关于前者,只要()x ϕ在[0,]l 上满足狄氏条件就可以了,现在主要谈求系数的问题.回忆傅氏系数公式的得来是根据函数系的正交性,所以现在也要考察函数系{}sin n x β在[0,]l 上的正交性,可以证明(参阅§5.3关于固有函数正交性的证明方法)sin sin 0,.lmn x xdx m n ββ=≠⎰令 20sin ,ln n L xdx β=⎰于是把1()sin n n n x C x ϕβ∞==∑ (4.23)的两端乘上sin k x β,然后在[0,l ]上积分得()sin lkk k x xdx L C φβ=⎰即 01()sin .Lk k kC x xdx L φβ=⎰(4.24)把(4.24)代入(4.22)式即得原定解问题的解.通过上面两节的讨论,我们对分离变量法已经有了一个初步的了解,它的主要步骤大体为:一、首先将偏微分方程的定解问题通过分离变量转化为常微分方程的定解问题,这对线性齐次偏微分方程来说是可以做到的.二、确定固有值与固有函数.由于固有函数是要经过叠加的,所以用来确定固有函数的方程与条件,当函数经过叠加之后仍旧要满足,当边界条件是齐次时,求固有函数就是求一个常微分方程满足零边界条件的非零解.三、定出固有值、固有函数后,再解其他的常微分方程,把得到的解与固有函数乘起来成为(,)n u x t ,这时(,)n u x t 中还包含着任意常数.四、最后为了使解满足其余的定解条件,需要把所有的(,)n u x t 叠加起来成为级数形式,这时级数中的一系列任意常数就由其余的条件确定,在这最后的一步工作中,需要把已知函数展开为固有函数项的级数,所以,必须考虑固有函数的正交性.由本节的例子还可以看出,用分离变量法求解第三类边界条件(第二类边界条件也一样)的定解问题时,只要边界条件都是齐次的,其过程与解第一类边界条件的定解问题是相同的,但在确定固有值时,一般说来是比较复杂的.4.3 圆城内的二维拉普拉斯方程的定解问题一个半径为0ρ的薄圆盘,上下两面绝热,圆周边缘温度分布为已知,求达到稳恒状态时圆盘内的温度分布.在第3章讲过,热传导问题达到稳恒状态时温度分布与时间无关,应满足拉普拉斯方程20.u ∇=因为边界形状是个圆周,它在极坐标下的方程为0ρρ=,所以在极坐标系下边界条件可表为(),u f ρρθ==既然边界条件使用了极坐标系,所以我们将方程也采用极坐标的形式,于是有220220110,;(4.25)(,)().(4.26)u uu u f ρρρρρρρθρθθ⎧⎛⎫∂∂∂∇=+=<⎪ ⎪∂∂∂⎨⎝⎭⎪=⎩此外,因为自变量,ρθ的取值范围分别是[0,0ρ]与[0,2]π,圆盘中心点的温度决不可能是无穷的,并且(,ρθ)与(,2ρθπ+)实际上表示同一点,温度应该相同,即应该有lim u ρ→<+∞ (4.27)(,)(,2).u u ρθρθπ=+ (4.28)现在来求满足方程(4.25 )及条件(4.26),(4.27),(4.28)的解.先令(,)()(),u R ρθρθ=Φ代入方程(4.25)得2110R R R ρρ'''''Φ+Φ+Φ= 即2,R R Rρρ'''''+Φ=-Φ令比值为常数λ即得两个常微分方程0,λ''Φ+Φ=20.R R R ρρλ'''+-=再由条件(4.27)及(4.28)可得(0),R <∞(2)().θπθΦ+=Φ (4.29)这样一来,我们得到了两个常微分方程的定解问题0,(2)().λθπθ''Φ+Φ=⎧⎨Φ+=Φ⎩(4.30) 与20,(0).R R R R ρρλ'''⎧+-=⎨<∞⎩(4.31) 先解哪能一个呢?要看哪一个可以定出固有值,由于条件(4.29)满足可加性(即所有满足(4.29)的函数叠加起来仍旧满足(4.29),所以只能先解问题(4.30).采用与4.1中同样的方法可以得到 当0λ<时,问题(4.30)无解;当0λ=时,它的解为00()a θ'Φ=(常数);当0λ>时,它的解为0(),n n a b θ''Φ=+且λ必须是整数n ,取1,2,3,n =(只取正整数的理由与4.1相同),则0()cos sin .n n a n b n θθθ''Φ=+至此,我们已经定出了固有值与固有函数,接下去是解问题(4.31),其中的方程是欧拉(Euler)方程,它的通解为000ln ,R c d ρ=+当λ=0;,n n n n n R c d ρρ-=+ 当2(1,2,3,).n n λ==为了保证(0)R <∞,只有0n d = (0,1,2,),n =即 0(0,1,2,)nn R c n ρ==,因此利用叠加原理,方程(4.25)满足条件(4.27),(4.28)的解可以表示为级数01(,)(cos sin )2nn n n n a u a n b ρθρθθ∞==++∑ (2.32)此式中的2a 就是00;,n n a c ab '分别是.n n n n ac b c '',最后为了确定系数,n n a b ,我们利用边界条件(4.26)得001()(cos sin )2nn n n n a f a n b θρθθ∞==++∑,(4.33) 因此,000,,n nn n a a b ρρ就是()f θ展开为傅氏级数时的系数,即有2002002001(),1()cos ,1()sin .n n n n a f d a f n d b f n d πππθθπθθθρπθθθρπ⎧⎪=⎪⎪⎪=⎨⎪⎪=⎪⎪⎩⎰⎰⎰ (4.34) 将这些系数代入(4.32)式即得所求的解.为了以后应用起来方便,我们还可以将解(4.32)写成另一种形式.为此,将(4.34)式所确定的系数代入(4.32)式经过简化后可得201011(,)()cos ().2n n u f t n t dt πρρθθπρ∞=⎡⎤⎛⎫⎢⎥=+- ⎪⎢⎥⎝⎭⎣⎦∑⎰(4.35)利用下面已知的恒等式221111cos ()2212cos()n n k k n t k t k θθ∞=-+-=--+∑*)(1),k <*)这个恒等式的证明如下:[]∑∑∞=∞=---++=-+11)()(2121)(cos 21n n t in t in n n e e k t n k θθθ [][]∑∑∞=∞=---++=11)()(212121n n nt i n t i ke ke θθ)()()((1211)2121t i t i t i t i ke ke ke ke -------+-+=θθθθ(|K|<1) ⎥⎦⎤⎢⎣⎡-+-----+-----+-+=)sin()cos(1)sin()cos()sin()cos(1)sin()cos(121t ik t k t ik t k t ik t k t ik t k θθθθθθθθ可将解(,)u ρθ(4.35)表达为()2220022001(,)(),.22cos()u f t dt t πρρρθρρπρρρρθ-=<+--⎰(4.36) 公式(4.36)称为圆域上的泊松公式,它的作用在于把解写成了积分形式,这样便于作理论上的研究.4.4 有界弦的强迫振动前面所讨论的偏微分方程都限于齐次的,现在要讨论非齐次方程的解法,为方便起见,以弦的强迫振动为例,所用的方法对其他类型的方程也适用.我们研究的问题是一根弦在两端固定的情况下,受强迫力作用所产生的振动现象.即要考虑下列定解问题22222000(,),0,0;(2.37)0;(2.38)(),().(2.39)x x l t t uu a f x t x l t t x u u u u x x t ϕψ====⎧∂∂⎪=+<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩在现在的情况,弦的振动是由两部分干扰引起的,一是强迫力,一是初始状态,所以由物理意义可知,此时振动可以分解为仅由强近力引起的振动和仅由初始状态引起的振动的合成.由此得到启发,我们可设解为(,)(,)(,),U x t V x t W x t =+ (4.40)其中(,)V x t 表示仅由强迫力引起弦振动的位移,它满足22222000(,),0,0;0;(2.41)0.x x l t t VV a f x t x l t t x V V u V t ====⎧∂∂⎪=+<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩而(,)W x t 表示仅由初始状态引起弦振动的位移,它满足.)cos(21121)cos(212)cos(21212222kt k k k t k k t k +---=⎥⎦⎤⎢⎣⎡+----+=θθθ22222000,0,0;0;(2.42)(),().x x l t t WW a x l t t x W W W W x x t φψ====⎧∂∂⎪=<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩读者不难验证,若V 是(4.41)的解,W 是(4.42)的解,则U V W =+一定就是原定解问题的解.问题(4.42)可直接用分离变量法求解,因此现在的问题只要讨论如何解问题(4.41)就行了.关于问题(4.41),我们可以采用类似于线性非齐次常微分方程中所常用的参数变易法,并保持如下的设想,即这个定解问题的解可分解为无穷多个驻波的叠加,而每个驻波的波形仍然是由该振动体的固有函数所决定,这就是说,我们假设(4.41)的解具有如下的形式1(,)()sin,n n n V x t v t x lπ∞==∑ (4.43) 其中()n v t 是待定的函数,为了确定()n v t ,将自由项(,)f x t 展成的傅氏正弦级数,即1(,)()sin,n n n f x t f t x lπ∞==∑ (4.44) 其中 02()(,)sin .l n n f t f x t xdx l lπ=⎰ 将(4.43)及(4.44)代入(4.41)的第一个式了,得到222''21()()()sin 0,n n n n a n n v t v t f t x l l ππ∞=⎡⎤+-=⎢⎥⎣⎦∑ (4.45) 由此可见得222''2()()().nn n a n v t v t f t lπ+= 再将(4.41)中的初始条件代入(4.43)得'(0)0,(0)0.n n v v ==这样一来,确定函数()n v t 只需解下列定解问题:2222()()()(0)0,(0)0n n n nn a n v t v t f t l v v π⎧''+=⎪⎨⎪'==⎩ 1,2,n =. (4.46) 用拉氏变换法解出(4.46),得到()()()sintn n ln a t v t f d n a lπτττπ-=⎰*),*)在方程(4.46)的两端取关于t 的拉氏变换,得所以,01()(,)()sin sin .t n n ln a t n v x t f d x n al l πτπττπ∞=-⎡⎤=⎢⎥⎣⎦∑⎰ 将这个解与问题(4.42)的解加起来,就得到原定解问题(4.37),(4.38),(4.39)的解.这里所给的求解问题(4.41)的方法,其实质是将方程的自由项及解都按齐次方程所对应的一族固有函数展开.随着方程与边界条件的不同,固有函数族也就不同,但总是把非齐次方程的解按相应的固有函数展开.所以这种方法也叫固有函数法.4.5 非齐次边界条件的处理前面所讨论的定解问题的解法,不论方程是齐次的还是非齐次的,边界条件都是齐次的.如果遇到非齐次边界条件的情况,应该如何处理?总的原则是设法将边界条件化成齐次的.具体地说,就是取一个适当的未知函数之间的代换,使对新的未知函数,边界条件是齐次的.现在以下列定解问题为例,说明选取代换的方法.设有定解问题2222212000(,),0,0;(4.47)(),();(4.48)(),().(4.x x l t t uu a f x t x l t t x u u t u u t u u x x t ϕψ====⎧∂∂⎪=+<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩49) 我们设法作一代换将边界条件化为齐次的,为此令(,)(,)(,).u x t V x t W x t =+ (4.50) 选取W(x,t)使V(x,t)的边界条件为齐次的,即0,0x x lVV==== (4.51)120(),()x x lWu t Wu t ==== (4.52)也就是说,只要所选取的W 满足(4.52)就能达到我们的目的.而满足(4.52)的函数是容易找到的,例如取W 为x 的一次式,即设(,)()(),W x t A t x B t =+ 用条件(4.52)确定(),()A t B t 得2111()()(),()()A t u t u t B t u t l=-=⎡⎤⎣⎦.),()()(22222p F p U ln a p U p n n n =+π 其中U n (p),F n (p)分别为v n (t)与f n (t)的拉氏变换,解出),(1)(22222p F l n a p p U nn π+= 由于222221l n a p π+的逆拉氏变换为lat n an l ππsin,利用拉氏变换的卷积性质,即得v n (t).显然,函数211()()(,)()u t u t W x t u t x l-=+就满足(4.52)式,因而只要作代换211,u u u V u x l -⎡⎤=++⎢⎥⎣⎦(4.53)就能使新的未知函数满足齐次的边界条件.经过这个代换后,得到V 的定解问题为22212201100(,),0,0;0,0;(2.54)(),(),x x l t t VV a f x t x l t t x V V V V x x t ϕψ====⎧∂∂⎪=+<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩其中211121112111()()(,)(,)();(0)(0)()()(0);(0)(0)()()(0).u t u t f x t f x t u t x l u u x x u x l u u x x u x l ϕϕψψ⎧⎡⎤''''-''⎪=-+⎢⎥⎪⎢⎥⎣⎦⎪-⎡⎤⎪=-+⎨⎢⎥⎣⎦⎪⎪⎡⎤''-⎪'=-+⎢⎥⎪⎢⎥⎣⎦⎩(4.55)问题(4.54)可用上一节的方法解出,将(4.54)的解代入(4.53)即是原问题的解.上面由(4.52)式定(,)W x t 时,取(,)W x t 为x 的一次式是为了使(4.55)中几个式子简单一些,并且(,)W x t 也容易定,但若12,,f u u 都与t 无关,则可适当的选择()W x (也与t 无关),使(,)V x t 的方程与边界条件同时都化成齐次的,这样做就可以省掉下面对(,)V x t 要进行解非齐次方程的繁重工作.这种()W x 究竟怎么找,将在后面的例题中说明.若边界条件不全是第一类的,本节的方法仍然适用,不同的只是函数(,)W x t 的形式,读者可就下列几种边界条件的情况写出相应的(,)W x t 来:1201201201(),();2(),();3(),().x x lx l x x x lu u u t u t xuu t u u t x u u u t u t xx======∂==∂∂==∂∂∂==∂∂ 以上各节我们说明了如何用分离变量法来解定解问题,为便于读者掌握此方法,现将解定解问题的主要步骤简略小结如下:一、根据边界的形状选取适当的坐标系,选取的原则是使在此坐标系中边界条件的表达式最为简单.圆、圆环、扇形等域用极坐标系较方便,圆柱形域与球域分别用柱坐标系与球坐标系较方便.二、若边界条件是非齐次的,又没有其他条件可以用来定固有函数,则不论方程是否为齐次,必须先作代换使化为具有齐次边界条件的问题,然后再求解.三、非齐次方程、齐次边界条件的定解问题(不论初始条件如何)可以分为两个定解问题,其一是具有原来初始条件的齐次方程,其二是具有齐次定解条件的非齐次方程.第一个问题用分离变量法求解,第二个问题按固有函数法求解.在结束本章之前,我们再举两个综合性的例子,目的是帮助读者掌握分离变量法的全过程.例2 求下列定解问题5)22222000,0,0;(2.56)0,;(2.57)0(2.8x x l t t uu a A x l t t x u u B u u t ====⎧∂∂⎪=+<<>∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩的解,其中,A B 均为常数.解 这个定解问题的特点是:方程及边界条件都是非齐次的.根据上述原则,首先应将边界条件化成齐次的,由于方程(2.56)的自由项及边界条件都与t 无关,所以我们有可能通过一次代换将方程与边界条件都变成齐次的,具体做法如下: 令 (,)(,)(),u x t V x t W x =+ 代入方程(4.56)得22222''().V V a W x A t x ⎡⎤∂∂=++⎢⎥∂∂⎣⎦为了使这个方程及边界条件同时化成齐次的,选()W x 满足⎪⎩⎪⎨⎧===+==BW W A x W a l x X ,0,0)(''02(4.59) (4.59)是一个二阶常系数线性非齐次方程的边值问题,它的解可以通过两次积分求得222().22A Al B W x x x a al ⎛⎫=-++ ⎪⎝⎭ 求出函数()W x 之后,再由(4.58)可知函数(,)V x t 为下列定解问题的解.22222000,0,0;(4.60)00;(4.61)(),0.(4.62)x x l t t VV a x l t tx V V V V W x t ====⎧∂∂⎪=<<>∂∂⎪⎪==⎨⎪∂⎪=-=⎪∂⎩采用分离变量法,令(,)()()V x t X x T t =,代入(4.60)得2,XT a X T ''''=或22.X T X a Tβ''''==- 由此得到2220,(4.63)0.(4.64)X X T a T ββ''+=''+=由(4.63)及(4.61)可以确定固有值与固有函数为2222.n l πβ=()sin .n n X x x lπ=再由(4.64)得()cossin ,n n n n a n aT t C t D t l lππ=+ 利用(4.62)中第二个条件可得0n D =.于是定解问题(4.60),(4.61),(4.62)的解可表示为1(,)cossin .n n n a n V x t C t x l lππ∞==∑ 代入(4.62)中第一个条件得1()sin,n n n W x C x lπ∞=-=∑ 即2221sin .22nn A Al B n x x C x a a l l π∞=⎛⎫-+= ⎪⎝⎭∑ 由傅氏级数的系数公式可得22202sin .22l n A Al B n C x x xdx l a a l l π⎡⎤⎛⎫=-+ ⎪⎢⎥⎝⎭⎣⎦⎰ 2222002sinsin llAn A B n x xdx x xdx a ll a l l ππ⎛⎫=-+ ⎪⎝⎭⎰⎰2223322222cos Al Al B n a n n a n ππππ⎛⎫=-++ ⎪⎝⎭(4.65)因此,原定解解问题的解为2221(,)cos sin ,22n n A Al B n a n u x t x x C t x a a l l l ππ∞=⎛⎫=-+++ ⎪⎝⎭∑其中n C 由(4.65)确定. 例3在环形域a b ≤≤内求解下列定解问题22222212();0,0.a b u u x y x y ⎧∂∂+=-⎪∂∂⎪⎨⎪==⎪⎩解 由于求解区域是环形区域,所以我们选用平面极坐标系,利用直角坐标与极坐标系之间的关系cos ,sin ,x y ρθρθ=⎧⎨=⎩ 可将上述定解问题用极坐标,ρθ表示出来22211()122cos 2,;(4.66)0,0.(4.67)a b u ua b u u ρρρρθρρρρρθρ==⎧∂∂∂+=<<⎪∂∂∂⎪⎨∂⎪==⎪∂⎩这是一个非齐次方程附有齐次边界条件的定解问题,采用4.4所述的固有函数法,并注意在4.3中得到的关于圆域内拉普拉斯方程所对应的固有函数,可令问题(4.66),(4.67)的解为(,)()cos ()sin ,n n n u A n B n ρθρθρθ∞==+⎡⎤⎣⎦∑代入(4.66)并整理得到22''''''222011()()()cos ()()()sin 12cos ,n n n n n n n n n A A A n B B B n ρρρθρρρθρθρρρρ∞=⎫⎧⎡⎤⎡⎤⎪⎪+-++-=⎨⎬⎢⎥⎢⎥⎪⎣⎦⎣⎦⎪⎩⎭∑比较两端关于cos ,sin n n θθ的系数可得'''2222214()()()12,A A A ρρρρρρ+-= (4.68)2''''21()()()0(2),nnn n A A A n ρρρρρ+-=≠ (4.69)2'''21()()()0.nnn n B B B ρρρρρ+-= (4.70)再由条件(4.67)得()()0,n n A a B a == (4.71)()()0.n n A b B b ''== (4.72)由(4.69),(4.70),(4.71),(4.72)可知()0(2),()0.n n A n B ρρ≡≠≡下面的任务就是要确定2().A ρ方程(4.68)是一个非齐次的欧拉方程,它的通解为224212(),A C C ρρρρ-=++由条件(4.71),(4.72)确定12,C C 得661442,a b C a b +=+4422244(2).a b a b C a b-=-+ 因此664422224244442(2)(),a b a b a b A a b a b ρρρρ-+-=--+++原定解问题的解为()66244222444441(,)2(2)()cos 2.u a b a b a b a b a bρθρρρθ-⎡⎤=-++--+⎣⎦+习 题 四1、设弦的两端固定于0x =及x l =,弦的初始位移如图所示,初速度为零,又没有外力作用,求弦作横向振动时的位移函数(,)u x t .2、就下列初始条件及边界条件解弦振动方程(,0)0,u x =(,0)(),(0,)(,)0.u x x l x tu t u l t ∂=-∂== 3、就下列初始条件及边界条件解弦振动方程01,0211,1;2t x x u x x =⎧<≤⎪⎪=⎨⎪-<<⎪⎩01(1);0.t x x u x x tu u ===∂=-∂==4、解出习题三中第2题.5、试求适合于下列初始条件及边界条件的一维热传导方程的解00(),0.t x x t u x l x u u ====-==6、解一维热传导方程,其初始条件及边界条件为00,0,0.t x x tu x u u xx====∂∂==∂∂7、一根长为l 的细杆表面绝缘,其初始温度分布如图所示,由0t =开始两端温度保持于零度,求杆上温度分布.8、试解出具有放射衰变的热传导方程2220,xu u a Ae x tα-∂∂-+=∂∂ 已知边界条件为00,0,x x l u u ====初始条件为0t u T ==,T 为常数9、求下列定解问题22200;0;0x x l t u u a A t x u u u ===⎧∂∂=+⎪∂⎪⎪==⎨⎪=⎪⎪⎩的解.10、求满足下列定解条件的一维热传导方程的解010,5,x x l u u ==== 0,t u ks k ==为常数11、试确定下列定解问题22200();,;()x x l t u u a f x t x u A u B u g x ===⎧∂∂=+⎪∂∂⎪⎪==⎨⎪=⎪⎪⎩解的一般形式.12、求稳恒状态下,由直线120,,0,x x l y y l ====所围矩形板内各点的温度,假设在10,,0x x l y ===三边上温度保持为零度,2y l =这边上各点温度为(),x ϕ其中1(0)()0.l ϕϕ== 13、一半径为a 的半圆形平板,其圆周边界上的温度保持(,)()u a T θθπθ=-,而直径边界上温度保持为零度,板的侧面绿缘,试求稳恒状态下的温度分布规律(,).u ρθ14、一圆环形平板,内半径为1r ,外半径为2r ,侧面绝缘,如内圆温度保持零度,外圆温度保持1度,试求稳恒状态下的温度分布规律(,).u r θ 15、如何解下列定解问题2222212000();,;(),().x x l t t uu a f x t x u M u M u u x x t ϕψ====⎧∂∂⎪=+∂∂⎪⎪==⎨⎪∂⎪==⎪∂⎩16、在矩形域内求下列定解问题2120120(,);(),();(),().x x a y y b u f x y u y u y u x u x ϕϕψψ====⎧∇=⎪⎪==⎨⎪==⎪⎩ 的解.17、在扇形区域内求下列定解问题200;0;()a a u u u u f θθρθ===⎧∇=⎪⎪==⎨⎪=⎪⎩ 的解.18、在矩形域0,0x a y b ≤≤≤≤内求拉普拉斯方程的解,使满足边界的条件:000,;0,0.x x a y y b u u Ay uu y y ====⎧==⎪∂∂⎨==⎪∂∂⎩ 19、把高频输电线充电到具有电压E ,然后一端短路封闭,另一端仍保持断开,求以后的电压分布.20、求矩形膜振动的位移,即解下列定解问题22222220000;0;0;()(),0.x x a y y b y t u uu a tx y uu u u u u xy x a y b t ======⎧⎛⎫∂∂∂=+⎪ ⎪∂∂∂⎝⎭⎪⎪==⎪⎨==⎪⎪∂⎪=--=⎪∂⎩。
拉普拉斯方程的解——分离变量法
∵
∫ sin
0
b
0 mπy nπy sin dy = b b b / 2
b
n ≠ m b = δ mn sin x (正交归一性) n=m 2
∴
∞ mπy b dy = C n δ mn = C m b / 2 ∑ ∫0 b 2 n =1 b m πy 2 2V b mπ C m = ∫ V sin dy = ⋅ sin y ′dy ′ b 0 b b mπ ∫0 4V (m = 奇数) 2V mπ =[ [ − cos y ′] 0 = mπ mπ (m = 偶数) 0 ∞ 4V 1 mπy − mπx / b sin e ϕ ( x, y ) = ∑ b π m =1,3,5L m
S
∂ϕ 2 ∂n
表面无自由电荷。 V
S
z=l
O y
四.应用实例(习题课)
1. 两无限大平行导体板,相距为 l ,两板间电势差为 V (与 x, y , z 无关),一板接地,求两板间的电势 ϕ 和 E 解: (1)边界为平面,故应选直角坐标系 下板接地 ϕ
S1
r
x
= 0 ,为参考点
(2)定性分析:由于在 z = l 处, ϕ = V 常数,可考虑 ϕ 与 x, y 无关。 (3) 列出方程并给出解:在 0 < z < l 区域, ( 4) ∇ ϕ = 0
(ϕ (r = a ) ≡ 0) ) 。
选柱坐标系: 对称性分析: ① 导体为圆柱,柱上电荷均匀分布, ϕ 一定与 θ 无关。 ② 柱外无电荷,电力线从面上发出后,不会终止到面上,只能 终止到无穷远,且在导体面上电场只沿 er 方向,可认为 ϕ 与 z 无关,
y r θ o z x
r
静电场中的镜像法与分离变量法
静电场中的镜像法与分离变量法摘要:静电场的基本问题是求解给定边界条件下的泊松方程或拉普拉斯方程的解,本文分别阐述在求解区域内有和没有自由电荷分布的情况下,应用镜像法和分离变量法求解;同时,举例来演示应用镜像法和分离变量法的解题思路、步骤和结果讨论以及一些注意点,并在相同情况下分别应用镜像法和分离变量法进行对比讨论;深入理解镜像法和分离变量法及其特征。
关键词:静电场、镜像法、分离变量法。
The Method of Mirror Image and the Separate Variational Method in the Electrostatic FieldAbstract: The basic problem of electrostatic field is to explore the solution of Poisson equation or Laplace equation under its given boundary condition. This article respectively explains the approaches to explore the solution using mirror image and separate variational methods under the to-be-explore solution area situation which has and which lacks freedom electric charge distribution .Meanwhile, it takes some instances to demonstrate the problem-solving thoughts and steps applying mirror image and separate variational methods. It also provides some discussions about the result and the points needing to be noted in the process of this demonstration.This writer also tries to help the readers todeeply understand the methods of mirror image and separate variational methods and their characteristics by doing contrast discussion under the same condition. Keywords: the electrostatic field, the method of mirror image, the separate variational method.1、引言:静电场和电源外恒定电场的边值问题的求解,可归纳为在给定边界条件下,对拉普拉斯方程或泊松方程的求解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
其特解之一为:
2
4
V
1
( x ')
r
dV
有限区域分布电荷,选无限远处电势为零时的解。
山东大学物理学院 宗福建
17
§2.4 分离变量法
对一般情况,设泊松方程的解为:
则,
即:
1 ( x ') ' dV 4 V r '0
2
山东大学物理学院 宗福建 23
§2.4 分离变量法
Pn(cosθ )为勒让德函数
山东大学物理学院 宗福建
24
§2.4 分离变量法
n=0时,
n
bn (an R n 1 ) Pn (cos ) R n b0 a0 R 点电荷的势
山东大学物理学院 宗福建
25
§2.4 分离变量法
n
Pn(cosθ )为勒让德函数,an和bn由边界条件确定。
山东大学物理学院 宗福建
22
§2.4 分离变量法
Pn(cosθ )为勒让德函数
P0 ( x) 1 P0 (cos( x)) 1 P1 ( x) x P1 (cos( x)) cos( x) 1 2 1 P2 ( x) (3x 1) P2 (cos( x)) (3cos 2 ( x) 1) 2 2 1 3 1 P3 ( x) (5 x 3x) P3 (cos( x)) (5cos3 ( x) 3cos( x)) 2 2 1 dl 2 l Pl ( x) l ( x 1) l 2 l ! dx
1 Q 4 r
2
Q 4 r E
2
j E
I 4 r 2 j
I 4 r j 4 r E Q
2
1 I 4 r
1 ( 4 ( x x0 )2 ( y y0 )2 ( z z0 )2 1 ) ( x x0 )2 ( y y0 )2 ( z z0 )2
山东大学物理学院 宗福建
19
§2.4 分离变量法
拉普拉斯(Laplace)方程的通解可以用分离变量法
求出。先根据界面形状选择适当的坐标系,然后在
该坐标系中由分离变量法解拉普拉斯方程。最常用 的坐标系有球坐标系和柱坐标系。这里我们写出用 球坐标系得出的通解形式(见附录Ⅱ)。球坐标用 (R,θ,φ)表示,R为半径,θ为极角,φ为方位角。
I
山东大学物理学院 宗福建 9
上一讲习题解答
1. 设容器壁为良导体金属,则容器壁为等势 体 。为满足器壁(边界上) 0 0 ,应有 三对像电极(即三个像电流)。
山东大学物理学院 宗福建
10
上一讲习题解答
1. 设容器壁为良导体金属,则容器壁为等势体 。为 满足器壁(边界上)0 0 ,应有三对像电极(即 三个像电流)。
14
山东大学物理学院 宗福建
上一讲习题解答
2. 设容器壁为绝缘体,则流向容器壁的电流 为零, 为满足器壁(边界上)Jn=0 ,应有三 对像电极(即三个像电流)
山东大学物理学院 宗福建
15
本讲主要内容
分离变量法求解泊松方程
山东大学物理学院 宗福建
16
§2.4 分离变量法
我们知道静电场标势所满足的泊松方程为:
1 ( x x0 )2 ( y y0 )2 ( z z0 )2
1 ( x x0 )2 ( y y0 )2 ( z z0 )2 1 ( x x0 )2 ( y y0 )2 ( z z0 )2
1 ( x x0 )2 ( y y0 )2 ( z z0 )2
' 0
泊松方程的解为拉普拉斯方程的通解+泊松方程特解
山东大学物理学院 宗福建
18
§2.4 分离变量法
在许多实际问题中,静电场是由带电导体决定的。
例如电容器内部的电场是由作为电极的两个导体板
上所带电荷决定的;又如电子光学系统的静电透镜 内部,电场是由于分布于电极上的自由电荷决定的。 这些问题的特点是自由电荷只出现在一些导体的表 面上,在空间中没有其它自由电荷分布。因此,如 果我们选择这些导体表面作为区域V的边界,则在V 内部自由电荷密度 ρ = 0 ,因而泊松方程化为比较 简单的拉普拉斯(Laplace)方程 。
山东大学物理学院 宗福建
( R R0 ) ( R R0 )
30
§2.4 分离变量法
边界条件包括:
bn 1 (an R n1 ) Pn (cos ) R n 0
n
(1)无穷远处, E → E0
1 |R E0 R cos E0 RP 1 (cos )
11
山东大学物理学院 宗福建
上一讲习题解答
1. 设容器壁为良导体金属,则容器壁为等势 体 。为满足器壁(边界上)0 0 ,应有 三对像电极(即三个像电流)。
山东大学物理学院 宗福建
12
上一讲习题解答
2. 设容器壁为绝缘体,则流向容器壁的电流为零, 为满足器壁(边界上)Jn=0 ,应有三对像电极(即 三个像电流)
n
(2)
1 1 (3cos 2 ( x) 1) D33 3 12 0 R 2
山东大学物理学院 宗福建 27
§2.4 分离变量法
例1
电容率为 ε
的介质球置于均匀外电场 E0中,
求电势。 (取介质球球心处电势为零。)
山东大学物理学院 宗福建
28
§2.4 分离变量法
解
介质球在外电场中极化,在它表面上产生束缚
设有两平面围成的直角形无穷容器,其内充 满电导率为 的液体,取该两平面为 xz面和 yz 面,在 ( x0 , y0 , z0 ) 和 ( x0 , y0 , z0 ) 两点分别置 正负电极并通以电流 I ,求导电液体中的电 势。
山东大学物理学院 宗福建
8
上一讲习题解答
解:无限空间时
山东大学物理学院 宗福建
20
§2.4 分离变量法
拉氏方程在球坐标系中的通解为
bnm m ( R, , ) (anm R n 1 ) Pn (cos ) cos m R n.m d nm m n (cnm R n 1 ) Pn (cos )sin m R n,m
山东大学物理学院 宗福建
13
上一讲习题解答
2. 设容器壁为绝缘体,则流向容器壁的电流为零, 为满足器壁(边界上)Jn=0 ,应有三对像电极(即 三个像电流)
I 4 ( 1 ( x x )2 ( y y )2 ( z z ) 2 0 0 0 1 ( x x )2 ( y y )2 ( z z )2 0 0 0
n
式中 a
n m
,b
n m
,c
n m
和 d
n m
为任意常数,在
具体问题中由边界条件定出。Pnm(cosθ ) 为缔和 勒让德(Legendre)函数。
山东大学物理学院 宗福建
21
§2.4 分离变量法
若该问题中具有对称轴,取此轴为极轴,则电势不
依赖于方位角φ ,这情形下通解为
bn (an R n1 ) Pn (cos ), R n
n
(2)取介质球心处电势为零。
2 |R 0 0
c0 0
n 1
d n 0 (n 0,1, 2,...)
n
2 cn R Pn (cos )
山东大学物理学院 宗福建 32
§2.4 分离变量法
边界条件包括:
(3)在介质面上: (R=R0)
bn P (cos ) ( R R0 ) 1 (cos ) n 1 n 1 a0 E0 RP n 0 R 2 cn R n Pn (cos ) ( R R0 ) n 1 1 2 1 |R R0 2 |R R0 0 R R R0 R R R
《电动力学》第12讲
第二章 静电场(4)
§2.4 分离变量法
教师姓名: 宗福建 单位: 山东大学物理学院 2015年10月23日
上一讲复习
1、电像法的适用条件
我们设想,导体面上的感应电荷对空间中电场的影 响用导体内部某个或某几个假想电荷来代替。注意 我们在作这种代换时并没有改变空间中的电荷分布 (在求解电场的区域,即导体外部空间中仍然是只 有一个点电荷Q),因而并不影响泊松方程,问题的 关键在于能否满足边界条件。如果用这代换确实能 够满足边界条件,则我们所设想的假想电荷就可以 用来代替导体面上的感应电荷分布,从而问题的解 可以简单地表示出来。
0
山东大学物理学院 宗福建
33
§2.4 分离变量法
bn n a E R P (cos ) P (cos ) c R 0 0 1 n 0 P n (cos ) n 1 n 0 n 0 R0 n 1 b n 1 E P (cos ) (n 1) n P (cos ) nc R Pn (cos ) 0 1 n 0 n2 n R0 0 n 1 n 0 b0 b0 P0的系数: a 0 + 0 2 0 R0 R0
1 ( x x0 )2 ( y y0 )2 ( z z0 )2 1 ( x x0 )2 ( y y0 )2 ( z z0 )2 1 ( x x0 )2 ( y y0 )2 ( z z0 )2