数学物理方程与定解条件
数学物理方程_定解问题
根据边界条件确定任意函数 f：
令 故
规定，当
时
4、定解问题是一个整体
达朗贝尔公式的求解过程，与大家熟知的常 微分方程的求结果成完全类似。
但遗憾的是，绝大多数偏微分方程很难求出 通解；即是求出通解，用定解条件确定其中待 定函数往往更为困难。这说明，达朗贝尔公式 不适用于普遍的数学物理定解问题的求解？
（7.1.8）
称式（7.1.8）为弦的自由振动方程。
(2) 如果在弦的单位长度上还有横向外力 作用，则式(7.1.8)应该改写为
(7.1.9)
式中
称为力密度，为 时刻作用于
处单位质量上的横向外力
式（7.1.9）称为弦的受迫振动方程.
2、 均匀杆的纵振动
一根杆，只要其中任一小段做纵向移动，必然使 它的邻段压缩或伸长，这邻段的压缩或伸长又使 它自己的邻段压缩或伸长。这样，任一小段的纵 振动必然传播到整个杆,这种振动的传播是纵波.
泊松方程和拉普拉斯方程的定解条件不包含初始条件, 而只有边界条件. 边界条件分为三类:
1、在边界上直接给定未知函数 , 即为第一类边界条件．
2、在边界上给定未知函数导数的值,即为第二类边界条件．
3、在边界上给定未知函数和它的导数的某种线性组合, 即第三类边界条件.
第一、二、三类边界条件可以统一地写成
第二类边界条件 规定了所研究的物理量在边界外法线方向上方向导数 的数值
u n
x0 , y0 ,z0
f (x0, y0, z0,t)
（7.2.3）
第三类边界条件 规定了所研究的物理量及其外法向导数的线性组合在 边界上的数值
（7.2.4）
其中 是时间 的已知函数， 为常系数．
7.2.2 泊松方程和拉普拉斯方程的定解条件
数学物理方程及其定解问题
3.定解问题的整体性(除上述两种类型外的 数学物理方程)
4.定解问题的适定性
4
一. 无界弦的自由振动
1. 无界弦的自由振动 (1)无界弦的含义：无界弦不是指无限长的弦，是指所关 心的那一段弦远离两端，在所讨论的时间内，弦两端的影响来 不及传到这段弦上，因而认为弦的两端在无限远，就不必给弦 的两端提出边界条件。 定解问题 初值问题
x at, x at
得方程的通解
u f1 ( x at) f 2 ( x at)
通解的物理意义: f2 ( x at ) 正行波, f2 ( x at ) 反行波
6
⑵ 利用定解条件来确定函数 f1 ( x), f 2 ( x)
由初始条件得
u ( x, 0) f1 ( x) f 2 ( x) ( x) ut ( x, 0) af1 ( x) af 2 ( x) ( x)
数学物理方程及其定解问题数学物理方程习题解答数学物理方程数学物理方程谷超豪数学物理方程pdf数学物理方程试卷数学物理方程视频数学物理方程答案数学物理方程第三版数学物理方程讲义
第七章 数学物理方程及其定解问题
1.数学物理方程的导出 2.定解条件 3.数学物理方程的分类 4.达朗贝尔公式 定解问题
1
3.数学物理方程的分类
15
三. 一般情况下的数学物理方程
一般情况下，无法像对无限长弦那样，先求通解，然后用定解条件 求特解。
定解问题的整体性
物理问题
数学问题
定解问题是一个整体
四 . 定解问题的适定性
如定解问题满足 (1) 有解 (2) 解是唯一的 (3) 解是稳定的 则称此定解问题是适定的。 因为定解问题来自实际。
习题11数学物理方程和定解条件
ρ + ε1Δ ρ ϕ +ε 2Δϕ
（ 0 < ε1 < 1 ， 0 < ε 2 < 1 ） ，
即
1
( ρ + Δρ )
Δu Δρ
−ρ
ρ + Δρ
Δu Δρ
ρ
ρ
Δρ
1 ∂ ρ ∂ρ
+
Δu 1 Δϕ
−
ϕ + Δϕ
Δu Δϕ
ϕ
ρ2
Δϕ
−
ρ m ∂ 2u
T ∂t 2
ρ + ε1Δρ ϕ + ε 2 Δϕ
=0
205．在铀块中，除了中子的扩散运动外，还进行着中子的吸收和增殖过程。设在单位时间 内单位体积中，吸收和增殖的中子数均正比于该时刻该处的中子浓度 u ( r , t ) ，因而净增中 子数可表为 α u ( r , t ) ， α 为比例常数。试导出 u ( r , t ) 所满足的方程。 用 q 表示单位时间流过某单位面积的中子数，有 q = − D∇u 。取一个六面体
− sin θ
θ +Δθ
∂u ⎤ 1 ⎛ ∂u k r + Δ ⎜ ⎥ ∂θ θ ⎦ Δϕ ⎜ ⎝ ∂ϕ
−
ϕ +Δϕ
= ρ ca 2 sin 2 θΔr
令 Δr , Δθ , Δϕ , Δt → 0 ，因为
Δu 。 Δt ⎡ ∂u ⎢sin (θ + Δθ ) ∂θ ⎣ − sin θ
θ +Δθ
1 Δθ
∂u ∂x
= 0 。由于左端点固定，故有 u
x=l
x=0
=0。
令（a）式中 t = 0 有 F − E S
∂u ∂x
第一章 三类典型方程和定解条件
a 其中，ij (x), bi (x), c x , f (x)都只是 x1 , x2, , xm 的已知 函数，与未知函数无关。
若一个函数具有某偏微分方程中所需 要的各阶连续偏导数，并且代入该方程中 能使它变成恒等式，则此函数称为该方程 的解（古典解）。 初始条件和边界条件都称为定解条件。 把某个偏微分方程和相应的定解条件 结合在一起，就构成了一个定解问题。 只有初始条件，没有边界条件的定解问题 称为始值问题（或柯西问题）。反之，只 有边界条件，没有初始条件的定解问题称 为边值问题。既有初始条件又有边界条件 的定解问题，称为混合问题。
数学物理方程
第一章 三类典型方程和定解条件 第二章 分离变量法 第三章 Laplace方程的格林函数法
第四章 贝塞尔函数及勒让德多项式
第一章 三类典型方程和定解条件
数学物理方程的研究对象——定解问题。 一个定解问题是由偏微分方程和相应的定解 条件组成。我们先来介绍三类典型的方程：
三类典型方程
一、波动方程 二、热传导方程
用以说明初始状态的条件称为初始条件。 用以说明边界上的约束情况的条件称为边 界条件。
一、初始条件
比如说波动方程（1.3）其初始条件有两 个，一个是参数u，一个是u的一阶导数。 即： u u t 0 及 都已知。 t
t 0
而热传导方程（1.7）其初始条件只有一 个，就是参数u。即：
Байду номын сангаасu t 0 是已知。
一个定解问题提的是否符合实际情况，从 数学角度来看，有三方面可以加以检验：
1、解的存在性，看定解问题是否有解。
2、解的唯一性，看是否只有一个解。
3、解的稳定性，看当定解条件有微小
变动时，解是否相应地只有微小的变 动，若确实如此，则称此解是稳定的。
数学物理方程及定解问题
这个初始问题有解
u( x, y) n2 sinh ny sin nx
数学物理方法2015.02
D 为扩散系数
数学物理方法2015.02
第三节 位势方程
稳定的温度场
a2u f ( x, y, z)
膜平衡方程
2u 2u a 2 2 f ( x, y ) x y
2
数学物理方法2015.02
第三节 位势方程
定解条件与定解问题的提法
第一类边界条件： u( x, y, z) g( x, y, z)
b u t2 b t2 u dx dx dt f0dx T0 ux (b, t ) ux (a, t ) dt a a t1 a t1 t t t2 t t t1 b
dt
t1
t2
b
a
数学物理方法2015.02
物理模型 在三维空间中，考虑均匀的、各向同性的物体。 假定它的内部有热源或汇，并且与周围的介质 有热交换，来研究物体内部温度的分布规律。 均匀物体： 物体的密度为常数 各向同性： 物体的比热和热传导系数均为常数
数学物理方法2015.02
第二节 热传导方程与扩散方程
数学模型的建立 设： u(x,y,z,t)表示物体于时刻 t 在位置 x,y,z 处的温度 C 表示是比热 (焦耳/度· 千克) 表示密度 (千克/米3)， k 表示导热系数 f 0(x,y,z,t)表示热源强度(焦耳/千克· 秒)
数学物理方法2015.02
第一节 波动方程及定解条件
三维波动方程或声波方程
2 2u 2u 2u 2 u a 2 2 2 f ( x, y , z , t ) 2 t z x y
第1章 数学物理方程及定解问题
2
T
ρ
, f (x, t) =
F(x, t)
ρ
, 得 力 用 ,弦 动 程 外 作 下 振 方 为
一维非齐次波动方程
∂ 2 u( x , t ) ∂ 2 u( x , t ) − a2 = f ( x , t ). 2 2 ∂t ∂x
二维波动方程或膜振动方程
一块均匀的紧张的薄膜，离开静止水平位置作垂直 于水平位置的微小振动，其运动规律满足
2 ∂ 2u ∂ 2u 2∂ u = a 2 + 2 + f ( x, y , t ) 2 ∂t ∂y ∂x
在时刻t , 弦段[ x , x + ∆x ]的动量为 x + ∆x ∂u( x , t ) ∫x ρ ∂t dx;
x + ∆x x
在时刻t + ∆t , 弦段[ x , x + ∆x ]的动量为 x + ∆x ∂u( x , t + ∆t ) dx . ∫x ρ ∂t
∫
=∫
∂u( x , t + ∆ t ) ∂u( x , t ) − ρ dx . ∂t ∂t
第一节 波动方程及定解条件
1.一维波动方程或弦振动方程 一维波动方程或弦振动方程
物理模型
一长为 l 的柔软、均匀的细弦，拉紧以后，让它离 的柔软、均匀的细弦，拉紧以后， 开平衡位置在垂直于弦线的外力作用下作微小横振 求弦上个点的运动规律。 动，求弦上个点的运动规律。
张紧的、静止的弦是一直线，该直线是弦的 平衡位置，以此为 x 轴。振动总是传播到整 根弦，横振动就是弦中的质点离开平衡位置 的位移垂直于 x 轴， 可用 t 时刻弦上各质点 x 离开平衡位置的横向位移 u ( x, t ) 来描述弦的 状态， 某一时刻 u ( x, t ) 的分布代表弦的形状， 称为位形。由于弦中质点的位移不同导致弦 的形变，形变产生应力，为了便于应力的描 述，不妨假定所研究的弦为“柔软的”弦。
数学物理方程 第一章典型方程和定解条件
sin ' tan ' u(x dx,t)
x
则
T T'
u
M'
ds
T'
'
M
gds
T
x
x dx x
T
u(
x dx, x
t)
u ( x, x
t
)
gds
ma
T
u(x dx,t) x
u ( x, x
t)
gds
ma
m ds
其中：
a 2u(x,t) t 2
ds dx
T
u(x dx,t) x
微小： 振幅极小， 张力与水平方向的夹角很小。
u
M'
ds
T'
'
M
gds
T
x
x dx x
牛顿运动定律：
横向：T cos T 'cos ' 0
纵向：T sin T 'sin ' gds ma 其中： cos 1 2 4 1
2! 4!
cos ' 1
sin tan u(x,t)
数学物理方程与特殊函数
☆ 数学与物理的关系
数理不分家
☆ 数学物理方程： 用数学方程来描述一定的物理现象
数学物理方程（简称数理方程）是指自然科学和工程技术的各门 分支学科中出现的一些偏微分方程(有时也包括积分方程、微分方程等)， 它们反映了物理量关于时间的导数和关于空间变量的导数 之间的制约关系。例如声学、流体力学、电磁学、量子力学等等 方面的基本方程都属于数学物理方程的研究对象。
• 如图，取杆长方向为x轴方向，垂直于杆长 方向的各截面均用平行位置x标记；在任一 时刻t，此截面相对于平衡位置的位移为u( x, t )
数学物理方程(谷超豪)课后习题完整解答
即对任何 x, G(x) C 0 又 G（x）=
1 1 x C ( x) ( )d x 2 2a 0 2a
x 2 ]的影响区域以外不发生变化；
（2） 在 x 轴区间[ x1 , x 2 ]上所给的初始条件唯一地确定区间[ x1 , x 2 ]的决定区 域中解的数值。 证： （1） 非齐次方程初值问题的解为 u(x,t)= [ ( x at ) ( x at )]
运动方程为:
2u 2u 2u t 2 x 2 y 2
2
x s x x
2u t
2
u u u ES x x ES x b x s x x x t t
利用微分中值定理，消去 x ,再令 x 0 得
2 2u u 2 u . b a t x 2 t 2
E
, 则得方程
所以 为原方程的通解。 由初始条件得
u
F x at G x at h x
1.
§2 达朗贝尔公式、 波的传抪 证明方程
2 2 x u 1 x 2 u h 0常数 1 1 2 2 h x h x a t
同理
5 2u 2 2 2 2 2 t x y t x 2 2 y 2 2 y
t
2 5 2 2 2 x y t 2x 2 y 2
u sin tg x. 2u u u [l ( x x)] ∣ x x g [l x] ∣ x g 2 x x t
1 F x Gx hx 1 x aF / x aG / x hx
数理方程重点总结
X (0) A 0 B 1 0
断 言: B 0， 于 是 有
u
u
0,
0 (2)
x x0
x xl
X ( x) A sin x
又 由 边 界 条 件u
0， 得
x xl
sin l 0
于 是 ， 得 到 空 间 变 量 问题 的 本 征 值
l n
或
n
( n l
)2
(n 1,2,3,)
据此，解得H( y)
H ( y) cos y 1 y2 1 H (0) 6
(7)
将 （5） 、 （7） 代 入 （4） 式 ， 即 得 特 解
u( x, y) 1 x3 y2 cos y 1 y2 1 x2
6
6
再另附：直接积分法 求偏微分方程的通解
2u u
t
2 2xt
xt x
可 以 由 两 个 边 界 条 件 唯一 地 被 确 定 。
例如 f (x) x
W (x)
1 6a 2
x3
C1 x C2
W (0) M1
M1 C2
W (l) M2
l3 M2 6a2 C1l M1
据此，得到W ( x) 的解
C1
M2
M1
l3 6a 2
l
M2
l
M1
l2 6a 2
X X 0
(1)
u x
0 , u
x0
x
0
xl
(2)
(1) 式的通解为
X ( x) Acos x B sin x
(3)
对上式求导，得
X ( x) A sin x B cos x
X ( x) A sin x B cos x
数学物理方法课件：第7章 数学物理方程定解问题
第七章 数学物理方程定解问题
§7.1 三类数学物理方程的导出 §7.2 定解条件 §7.3 数学物理方程的分类(自学） §7.4 达朗贝公式、定解问题
§7.1 三类数学物理方程的导出
（一）、梯度矢量
i
j
k
x y z
(i
j
k
) (i
j
k
)
x y z x y z
2 x2
2 y 2
2 z 2
令
2 2 2 x2 y2 z2
2 2 2 x2 y2 z2
记
utt
2u t 2
ut
u t
有时记
2
2 x2
2 y 2
u xx
2u x 2
2 2 2 3 x2 y2 z2
（二）、三类数学物理方程的导出
1、弦的横振动
弦的横向位移为 u(x,t)
dm ds T2 cos2 T1 cos1 0
(qx xdx qx x )dydzdt
qx dxdydzdt x
z
dx
y
dz
dy
(x, y, z)
x
x 方向净流入量为
qx dxdydzdt x
(D u )dxdydzdt x x
y 方向净流入量为
(D u )dxdydzdt y y
z 方向净流入量为 (D u )dxdydzdt z z
y
F (x,t)
M2
M1
1
T2
2
T2 sin 2 T1 sin 1 dmutt
T1
x
x+x
x
T2 sin 2 T1 sin 1 dsutt
T2 cos2 T1 cos1 0
数学物理方程学习指导书第3章经典方程建立及定解条件
第3章经典方程的成立和定解条件在议论数学物理方程的解法从前，我们第一要弄清楚数学物理方程所研究的问题应当如何提，为此，我们从双方面来议论，一方面要将一个详细的物理、力学等自然科学识题化为数学识题，即成立描绘某种物理过程的微分方程——数学物理方程，称此方程为泛定方程；另一方面要把一个特定的物理现象自己所拥有的详细条件用数学形式表达出来，即列出相应的初始条件和界限条件，二者合称为定解条件.定解条件提出详细的物理问题，泛定方程提供解决问题的依照，作为一个整体称之为定解问题.3．1 经典方程的成立在本节，我们将经过几个不一样的物理模型推导出数学物理方程中三种典型的方程，这些方程组成我们的主要研究对象.经典方程的导出步骤：（1）确立出所要研究的是哪一个物理量u；（2）用数学的“微元法”从所研究的系统中切割出一小部分，再依据相应的物理（力学）规律剖析周边部分和这个小部分间的作用（抓住主要作用，略去次要要素，即高等数学中的抓主部，略去高阶无量小），这类互相作用在一个短的时间间隔是如何影响物理量u3）把这类关系用数学算式（方程）表达出来，经化简整理就是所需求的数学物理方程.例1弦的振动弦的振动问题，固然是一个古典问题，但对于初学者仍旧拥有必定的启迪性.设有一根平均柔嫩的细弦，均衡时沿直线拉紧，并且除受不随时间而变的张力作用及弦自己的重力外，不受外力影响，下边研究弦的细小横向振动，即假定所有运动出此刻一个平面上，并且弦上的点沿垂直于x轴的方向运动（图3-1）.图3-1设弦上拥有横坐标为x的点，在时辰t时的地点为M，位移NM记作u.明显，在振动过程中位移u是变量x与t的函数u(x,t).此刻来成立位移u知足的方程.我们把弦上点的运动先看作小弧段的运动，而后再考虑小弧段趋于零的极限状况.在弦上任取一弧段MM，其长为ds，设是弦的线密度，弧段MM两头所受的张力记作T，T，此刻考虑孤段MM在t 时辰的受力状况，用牛顿运动定律，作用于弧段上任一方向上的力的总和等于这段孤的质量乘以该方向上的加快度.在x 轴方向弧段受力的总和为TcosTcos ，因为弦只作横向振动，所以TcosTcos0.（）假如弦的振动很小，并且在振动过程中弦上的切线倾角也很小，即0,0，则由24cos14!2!可知，当为无量小量时，cos 与1的差量是的高阶无量小量，能够略去不计，所以当0, 0时cos1,cos 1代入(3.1)式，即可近似获得TT .在u 方向弧段受力的总和为Tsin Tsingds ，此中是单位弧段的质量，gds 是弧段MM 的重力.又因当0， 0时 sin1 tg tgu(x,t)，tg 2xsin ' tg 'u(xdx,t)，x2ds 1u(x,t)dxdx,x2且小弧段在时辰 t 沿u 方向运动的加快度为u(x,t)，小弧段的质量为t 2TsinTsin2u(x,t)gdsds2t或gds ，所以（）Tu(xdx,t)u(x,t)gds2u(x,t)dx, x xt 2上式左侧方括号内的部分是因为x 产生dx 的变化而惹起的u(x,t)的改变量，可用微x分取代，即u(xdx,t)u(x,t)xu(x,t)dx2u(x,t)dx, xxxx 2于是T2u(x,t) gdx2u(x,t)dxx 2x 2或T2u(x,t) 2u(x,t) g.x2t22一般说来，张力较大时弧振动速度变化很快，即u要比g 大得多，所以又能够把g 略去.t 2经过这样逐渐略去一些次要的量，抓住主要的量，最后得出 u(x,t)应近似地知足方程2ua 22u（）t 2x 2这里的a2T.式(3.3)称为一维颠簸方程.假如在振动过程中，弦上此外还遇到一个与弦的振动方向平行的外力，且假定单位长 度所受外力的 F(x,t)，明显，在这里(3.1)及(3.2)分别为TcosTcos 0,FdsTsinTsin2u gdsds2.t利用上边的推导方法并略去弦自己的重量，可得弦的逼迫振动方程为2u22uf(x,t),（）’t2x2此中f(x,t)1F(x,t).方程(3.3)与(3.3)’的差异在于(3.3)’的右端多了一个与未知函数u 没关的项f(x,t)，这个项称为自由项，包括有非零自由项的方程称为非齐次方程，自由项恒等于零的方程称为齐次方程.(3.3)为齐次一维颠簸方程， (3.3)’为非齐次一维颠簸方程 .例2 传输线方程对于直流电或低频的沟通电，电路的基尔霍夫定律指出同一支路中电流相等 .但对于较高频次的电流（指频次还没有高到能明显地幅射电磁波的状况） ，电路中导线的自感和电容的效应不行忽视，因此同一支路中电流未必相等.现考虑一来一往的高频传输线，它被看作拥有散布参数的导体（图3-2）.在拥有散布参数的导体中，电流经过的状况，能够用电流强度与电压v来描绘，此处i与v都是x,t的函数，记作i(x,t)与v(x,t)，以R，L，C，G分别表示以下参数：R——每一回路单位的趾串连电阻，L——每一回路单位的串连电感，C——每单位长度的分路电容，G——每单位长度的分路电导.依据基尔霍夫第二定律，在长度为x的传输线中，电压降应等于电动势之和，即v(v v)R xii Lx.t而vx, vx故上式可写成v Ri L i.x T此外，由基尔霍夫第必定律，流入节点x的电流应等于流出该节点的电流，即i(i i)C x iGxv, t或i C vGv.x t将方程(3.4)与(3.5)归并，即得i与v应近似地知足以下方程组i C v Gv0,x tv iRi0.Ltx （）（（（（（（（（）为了确立函数i与v，将方程(3.5)对x微分，同时在方程(3.4)两头乘以C后再对t微分，并把两个结果相减，即得2i G v LC2i RC i x2x t20,t将(3.4)中的v代入上式，得x2i2i(RC i（）x2LG2GL)GRi,t t这就是电流i近似知足的微分方程，采纳近似的方法从（）与（）中消去i可得电压v近似知足的方程2v LG2v(RC GL)v GRv,（）x2t2t方程(3.6)或(3.7)称为传输线方程.依据不一样的详细状况，对参数R,L,C,G作不一样的假定，就能够获得传输线方程的各种特别形式.比如，在高频传输的状况下，电导与电阻所产生的效应能够忽视不计，也就是说可令G R0，此时方程(3.6)与(3.7)可简化为2i12ix2LC t2,2v12vt2LC t2.这两个方程称为高频传输线方程.若令a21这两个方程与(3.3)完整同样.因而可知，同一个方程能够用来描绘不一样的LC物理现象，一维颠簸方程不过颠簸方程中最简单的状况，在流体力学、声学及电磁场理论中，还要研究高维的颠簸方程.例3电磁场方程从物理学我们知道，电磁场的特征能够用电场强度E与磁场强度H以及电感觉强度D 与磁感觉强度B来描绘，联系这些量的麦克斯韦(Maxwell)方程组为rotH J D（）,trotE B,(3.9) tdivB 0,（）divD.(3.11)此中J为传导电流的体密度，为电荷的体密度.这组方程还一定与下述场的物质方程D eE,（）B H,（）J E,（）相联立，此中是介质的介电常数，是导磁率，为导电率，我们假定介质是平均并且是各向同性的，此时，，均为常数.方程(3.8)与(3.9)都同时包括有E与H，从中消去一个变量，就能够获得对于另一个变量的微分方程，比如先消去H，在(3,8)式两头求旋度并利用(3.12)与（）得rotrotH rotE rotE,t将(3.9)与(3.13)代入得rot rotH2H Ht2,t而rotrotH grad div2H,且divH1divB0,所以最后获得H所知足的方程为2H2HH t2t;同理，若消去H即得E所知足的方程2E2E E.t2t假如介质不导电(0)，则上边两个方程简化为2H12H,（）t22E12E,（）2t(3.15)与(3.16)称为三维颠簸方程.若将三维颠簸方程以标量函数的形式表示出来，则可写成2ua 2ua22u2u2u（t2x2y2z2,）此中a21,u是E或H的随意一个重量.从方程(3.11)与(3.12)还能够推导出静电场的电位所知足的微分方程.事实上，以(3.12)代入(3.11)得divD div E divE,而电场强度E与电位u之间存在关系E gradu,所以可得div(gradu)或2u,（）这个非齐次方程称为泊松（Poisson）方程.假如静电场是无源的，即0,则(3.18)变为2u 0，（）这个方程称为拉普拉斯（Laplace）方程.例4热传导方程一块热的物体，假如体内每一点的温度不全同样，则在温度较高的点处的热量就要向温度较低的点处流动，这类现象就是热传导.在工程技术上有很多传热问题都要归纳为求物体内温度的散布，此刻我们来推导传热过程中温度所知足的微分方程，与上例近似，我们不是先议论一点处的温度，而应当先考虑一个地区的温度.为此，在物体中任取一闭曲面S，它所包围的地区记作V（图3-3）.假定在时辰t，地区V内点M(x,y,z)处的温度为u(x,y,z,t)，n为曲面元素S的外法向（从V内指向V外）.图3-3由传热学可知，在t,t t时间内，从S流入地区V的热量与时间t，面积S，以及沿曲面的法线方向的温度变化率三者的乘积成正比，即Q kuk(gradu)n St Stnk(gradu)S t.此中k称为物体的热传导系数，当物体为平均导热体时，k为常数.于是，从时辰t1到时辰t2，经过曲面S流入地区V的所有热量为Q1t2dSdt.kgradut1S流入的热量使V内温度发生了变化，在△t时间内地区V内各点温度从u(x,y,z,t)变化到u(x,y,z,t+△t)，则在△t内V内温度高升所需要的热量为c[u(x,y,z,t t)u(x,y,z,t)]dVVu(x,y,z,t)ct tdV.V进而从时辰t1到时辰t2，因为温度高升所汲取的热量为t2uQ2cdV dt,t1tV此中c为物体的比热，为物体的密度，对平均物体来说，它们都是常数.因为热量守恒，流入的热量应等于物体温度高升所需汲取的热量，即t2t2ukgradudSdt c dVdt.t1t1tS V此式左端的由面积分中S是关闭曲面，能够利用奥-高公式将它化为三重积分，即kgradudS kdiv(gradu)dV k2udV,S V V 所以有t22udVdt t2udVdt.（）k ct1t1tV V因为时间间隔t,tt及地区V都是随意取的，并且被积函数是连续的，所以(3.20)式左右恒等的条件是它们的被积函数恒等，即u a22ua22u2u2u ,（）tx 2y 2z 2此中a 2k .方程(3.21)称为三维热传导方程.c若物体内有热源，其强度为F(x,y,z)，则相应的热传导方程为ua 22u2u2uf(x,y,z,t),tx 2y 2z 2此中fF .c作为特例，假如所考虑的物体是一根细杆 （或一块薄板），或许即便不是细杆（或薄板）而此中的温度u 只与x,t （或x,y,t ）相关，则方程(3.21)就变为一维热传导方程2u a 2u2;tx或二维热传导方程u a 22u2utx 22.y假如我们考虑稳恒温度场，即在热传导方程中物体的温度趋于某种均衡状态，这时温度u 已与时间t 没关，所以u 0,此时方程(3.21)就变为拉普拉斯方程(3.19).因而可知稳恒t温度场内的温度 u 也知足拉普拉斯方程 .在研究气体或液体的扩散过程时，若扩散系数是常数，则所得的扩散方程与热传导方程完整同样.3．2 初始条件与界限条件上边所议论的是如何将过程的物理规律用数学式子表达出来.除此之外，我们还需要把 详细条件也用数学形式表达出来， 这是因为任何一个详细的物理现象都是处在特定条件之下 的.比如弦振动问题，上节所推导出来的方程是全部柔嫩平均的弦作细小横向振动的共同规 律，在推导这个方程时没有考虑到弦在初始时辰物状态以及弦所受的拘束状况.假如我们不 是平常地研究弦的振动，必然就要考虑到弦所拥有的特定条件.因为任何一个详细振动现象 老是在某时辰的振动状态和此时辰从前的状态相关，进而就与初始时辰的状态相关.此外， 弦的两头所受的拘束也会影响弦的振动，端点所处的物理条件不一样会产生不一样的影响，因此弦的振动也不一样 .所以对弦振动问题来说，除了成立振动方程之外，还需列出它的详细条件对热传导方程，拉普拉斯方程也是这样.提出的条件应当恰好能够说明某一详细物理现象的初始状态以及界限上的拘束状况，.用以说明系统的初始状态的条件称为初始条件.用以说明界限上的拘束状况的条件称为界限条件.下边详细说明初始条件和界限条件的表达形式，先谈初始条件，对于弦振动问题来说，初始条件就是弦在开始时辰的位移及速度，若以(x),(x)分别表示初位移和初速度，则初始条件能够表达为u t0(x)u（）t(x) t0而对热传导方程来说，初始条件是指在开始时辰物体温度的散布状况，若以(M)表示t 0时物体内任一点M处的温度，则热传导方程的初始条件就是u(M,t)t0(M).(3.23)泊松方程与拉普拉斯方程都是描绘稳恒状态的，与初始状态无头，所以不提初始条件.再谈界限条件.假如界限条件直接给出了未知函数u(M,t)在界限S上的值，以s表示界限S上的动点，则这样的界限条件可表为u(M,t)MS(s,t),或简写成u S.（）这类界限条件称为第一类界限条件，此中(s,t)表示在界限S上给定的已知函数.比如，在杆的导热问题中，若在端点x a处温度保持为常数u0，这时在端点x a的界限条件为u xa u0.若在端点x a处温度随时间的变化规律f(t)为已知，在这点的界限条件为uxaf(t).又如在弦振动问题中，若弦的某端点x a是固定的，则在该点的位移为零，即uxa0.以上都是第一类界限条件的例子.总之，第一类界限条件直接给出了未知函数u(M,t)在边界S上的值但在很多状况下，界限上的物理条件其实不可以用第一类界限条件来描绘.比如，在杆的导热问题中，若杆的一端xa绝热，那么绝热这个条件就不可以直接给出杆的端点处的温度变化.因为从杆外经过杆端流入杆内的热量为kuSt(此中t为时间间隔，S为杆nxa的截面积，n为杆在端点x a处的外法向，若x a是杆的左端点，n的正向与x轴正向相反，则u u,若x a是杆的右端点，则n的正向与x轴正向同样，则uu), n x n x所以绝热这个条件能够表达为k uSt0, nxa即u0.nxa若在单位时间内经过x a端单位面积流入杆内的热量是t的已知函数f(t)，则这个条件可表示为k u f(t).nxa弦在对于弦振动问题来说，假如弦在x a处沿位移方向的张力（参照x a处是自由的，即沿着位移方向不受外力，中例1的推导）为则此时Tu0,nx a即u0.xx a总之，有时界限条件一定表达为u（）(s,t).n S的形式，此中u.表示函数沿界限外法向的变化率，这类界限条件称为第二类界限条件n除了上述两类界限条件外，有时还会碰到其余形式的界限条件.比如在杆的导热热问题中，若杆在某个端点x a自由冷却，那么自由冷却这个条件就是K uH(u1u xa), nxa(此中u1为四周介质的温度)即uu1h kuh.n xa H这是因为在单位时间内从四周介质传到杆的x a端单位面积上的热量与介质和杆端的温度差成正比，而在单位时间内经过u xa端单位面积传向杆内的热量与n考取例4）.成正比（参xa对于有界杆(0 x l)，若两头都是自由冷却，则在x l处，上述条件可表为uu h u1;nx i在x0处，这个条件可表为uu h u1.nx 0一般地，这类界限条件的形式为u(s,t).（）uhn s这样的界限条件称为第三类界限条件.无论哪一种界限条件，假如它的数学表达式中的右端自由项恒为零，则这类界限条件称为齐次的.定解问题的提法前方两节我们推导了三种不一样种类的偏微分方程并议论了与它们相应的初始条件与边界条件的表达方式.因为这些方程中出现的未知函数的偏导数的最高阶都是二阶，并且它们对于未知函数及其各阶偏导数来说都是线性的，所以这类方程称为二阶线性偏微分方程*）1 .在工程技术上二介线性偏微分方程碰到最多.假如一个函数拥有所需要的各阶连续编导数，并且代入某偏微方程中能使该方程变为恒等式，则此函数称为该方程的解.因为每一个物理过程都处在特定的条件之下，所以我们的任务是要求出合适初始条件和界限条件的解.初始条件和界限条件都称为定解条件.求一个偏微方程知足定解条件的解的问题称为定解问题.只有初始条件，没有界限条件的定解问题称为始值问题（或柯西问题）；而没有初始条件，只有界限条件的定解问题称为边值问题；既有初始条件也有界限条件的定解问题称为混合问题.一个定解问题提得能否切合实质状况，自然一定靠实质来证明，但是从数学角度来看，能够从三方面加以查验. 1）解的存在性，即看所结出来的定解问题能否有解；2）解的独一性，即看能否只有一个解；3）解的稳固性，即看当定解条件有细小改动时，解能否相应地只有细小的改动，假如*）二阶线性编微分方程能够按它们的二阶导数的系数的代数性质进行分类，在§中所推导的颠簸方程属于双曲型,拉普拉斯（或泊松）方程属于椭圆型，热传导方程属于抛物型，对于二阶线性偏微分方程的分类方法,读者可参阅复旦大学数学系编《数学物理方程》（第二版，上海科学技术第一版社第一版）第一章§5.确立这样，此解便称定的，否所得的解就无用价 .因定解条件往常是利用方法得的，因此所获得的果， 有必定的差，假如所以而解的很大， 那么种解然不可以切合客的要求 .假如一个定解存在独一且定的解，此称适定的，在此后中我把着眼点放在定解的解法上， 而极少它的适定性， 是因定解的适定 性常常十分困，而本所的定解都是古典的，它的适定性都是了然的 .习题一1. l 的平均杆，面，一端温度零，另一端有恒定流 q 入（即位内通位截面流入的量q ），杆的初始温度散布是x(ix)，写出相的定解.22. l 的弦两头固定，开始在xc 遇到冲量的作用，写出相的定解.有一平均杆，只需杆中任一小段有向位移或速度，必致段的或伸，种仲开去，就有波沿着杆播，推杆的振方程 .4.一平均杆原l ，一端固定，另一端沿杆的方向被拉e 而静止，忽然松手任其振，成立振方程与定解条件.若F(z),G(z)是随意二可微函数，uF(xat)G(xat)足方程2ua2 2ut 2x 26.若函数u 1(x,t),u 2(x,t), ,u n (x,t),⋯均性次方程2up2uq u r ux 2t 2xt的解，此中p,q,r 不过x,t 的函数，并且数uu k (x,t)收，并x,t 能够行两次k 1逐微分，求数uuk(x,t)足原方程（个叫做性次方程的叠加原理）.k1。
数学物理方法课件：7-数学物理定解问题
，a是弦的振动传播速度，则
utt a2uxx 0 (7.1.6)
如果，弦受到线密度为F(x,t) 的横向 力作用，弦 y方向方程应为：
T2ux xdx T1ux x F(x,t)dx (dx)utt
则弦的受迫振动方程为：
u
T2
B
α2
α1
C
T1 A
o x x+dx x
utt
a2uxx
F ( x, t )
y,
z, t )dxdydzdt
所以三维热传导方程为
cut
x
k
u x
y
k
u y
z
k
u z
F ( x,
y,
z,t)
17
三维 热传导方程
cut
x
k
u x
y
k
u y
z
k
u z
F ( x,
y,
z,t)
➢对于均匀物体，k、c、 ρ是常数
ut a23u f (x, y, z,t)
k k
x u
y u
k ux k uy k uz
z
15
确定物理量：温度的空间和时间分布u(x, y, z, t)
确定研究微元：x, x dxy, y dyz, z dz dV
数学物理方程复习
数学物理方程复习一．三类方程及定解问题（一）方程1.波动方程（双曲型）Utt = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1（t）;U(l,t)= Φ2（t）;U(x,0)= Ψ1（x）;Ut (x,0)=Ψ2（x）。
2.热传导方程（抛物型）Ut = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1（t）;U(l,t)= Φ2（t）;U(x,0)= Ψ1（x）.3.稳态方程（椭圆型）Uxx +Uyy=f; 0<x<a;0<y<b;t>0.U(0,x)= Φ1（x）;U(b,x)= Φ2（x）;U(y,0)= Ψ1（y）;Ut (y,a)=Ψ2（y）。
（二）解题的步骤1.建立数学模型，写出方程及定解条件2.解方程3.解的实定性问题（检验）（三）写方程的定解条件1.微元法：物理定理2.定解条件：初始条件及边界条件（四）解方程的方法1.分离变量法（有界区域内）2.行波法（针对波动方程，无界区域内）3.积分变换法（Fourier变换Laplace变换）Fourier变换:针对整个空间奇：正弦变换偶：余弦变换Laplace变换：针对半空间4.Green函数及基本解法5.Bessel函数及Legendre函数法例一：在弦的横震动问题中，若弦受到一与速度成正比的阻尼，试导出弦阻尼振动方程。
解：建立如图所示的直角坐标系，设位移函数为U（x,t）,取任意一小段△x进行受力分析，由题设，单位弦所受阻力为b U t（b为常数），在振动过程中有△x所受纵向力为：（T2COSa2-T1COSa1）横向力为：（T2SINa2-T1SINa1-b U t(x+n△x））(0<n<1). T2,T1为△x弦两端所受的张力，又因为弦做横振动而无纵振动，由牛顿定律有T2COSa2-T1COSa1=0，T2SINa2-T1SINa1-b(x+n△x)U t=p U tt(x+n△x)△x在小的振动下SINa1≈TANa1=Ux(x,t), SINa2≈TANa2=Ux(x+△x,t),COSa2≈COSa1≈1,T=T1=T2.(ρ是密度)即(T/ρ)[ Ux (x+△x,t)- Ux(x,t)]/ △x-(b/ρ) U t(x+n△x，t)即令△x→0时有：U tt+ aU t=a2U xx例二：设扩散物质的源强（即单位时间内单位体积所产生的扩散物质）为F （x,y,z,t），试导出扩散方程。
数学物理方法第七章
x dx
相对伸长
u x
x
u x
x
x dx
F
x dx
x
由虎克定律，B两端的 张应力（单位横截面 的力）分别为
u ( x)
u ( x dx)
A
B
C
u u u B段运动方程为 YS ) 2 x dx YS x ( Sdx x x t Y u x x dx u x x utt dx
u(x) u+u u 0 1
F B
T2 受力分析和牛顿运动定律： 2
沿x-方向，不出现平移
T2 cos 2 T1 cos 1 0
T1 x
x+x
（ 1）
沿垂直于x-轴方向
T2 sin 2 T1 sin 1 F ( x, t )dx ( dx )utt
在微小振动近似下：
如立方体内无源和汇 dt时间内粒子增加数为
u dxdydzdt (u t dt u t )dxdydz du x , y , z dxdydz t u dxdydz t
u u u ( D )dxdydz ( D )dxdydz ( D )dxdydz x x y y z z u u u u { [ ( D ) ( D ) ( D )]}dxdydz 0 t x x y y z z
7.1 数学物理方程的导出
步骤：
1、明确要研究的物理量是什么？
从所研究的系统中划出任一微元，分析邻近部 分与它的相互作用。 2、研究物理量遵循哪些物理规律？ 3、按物理定律写出数理方程（泛定方程）。
波动方程的导出
（一）均匀弦微小横振动方程 设：均匀柔软的细弦沿x轴绷紧，在平衡位置附
第五章 数理方程的建立,定解条件,傅里叶级数和傅里叶变换(简介),代尔塔函数的简介
ut = a2uxx + f ( x,t ) （有热源） 或 ut = a2uxx （无热源）。
若用 u 代表物体内某种物质的浓度。则扩散方程与热传导方程是一样的。
（3）泊松方程和拉普拉斯方程
若温度达到了稳定分布，即温度分布不随时间变化， ut = 0 ，则由热传 导方程可得温度稳定分布满足的方程为
dv
（积分形式），
∫v ∇ ⋅
JK Edv
=
1 ε
∫v
ρ dv
⇒
∇⋅
JK E
=
ρ ε
（微分形式）。
( ) 又
v∫l
JK E
⋅
d
K l
=
0
⇒
∫∫s
JK ∇×E
⋅
K ds
=
0
⇒
∇
×
JK E
=
0
，
∇
×
JK E
=
0
⇒
JK E
=
−∇u
，
u
为静电势（无旋场必为梯度场）
72
∴∇2u = − ρ ， 静电势满足的方程为泊松方程。
74
初始条件所反映的必须是物体上各点的初始状态，而不是仅仅某一点。
边界条件
一共有三类边界条件：
1.给定要求解的函数 u 在边界上的值 u ( x, y, z,t ) 边界上的x0 , y0 , z0
=
f
( x0, y0, z0,t ) ，称为
第一类边界条件。
例如：若研究长为 l 、两端固定的弦的振动情况，既然弦的两端 x = 0 ， x = l
所满足的方程。
设弦的质量密度为 ρ ，现在研究位于 x 到 x + dx 这一段弦的运动状况。这
数学物理方程学习指导书第3章经典方程的建立和定解条件
第3章经典方程的建立和定解条件在讨论数学物理方程的解法以前，我们首先要弄清楚数学物理方程所研究的问题应该怎样提，为此，我们从两方面来讨论，一方面要将一个具体的物理、力学等自然科学问题化为数学问题，即建立描述某种物理过程的微分方程一一数学物理方程，称此方程为泛定方程；另一方面要把一个特定的物理现象本身所具有的具体条件用数学形式表达出来，即列出相应的初始条件和边界条件，两者合称为定解条件•定解条件提出具体的物理问题，泛定方程提供解决问题的依据，作为一个整体称之为定解问题3. 1经典方程的建立在本节，我们将通过几个不同的物理模型推导出数学物理方程中三种典型的方程，这些方程构成我们的主要研究对象•经典方程的导出步骤：(1)确定出所要研究的是哪一个物理量u ；(2)用数学的“微元法”从所研究的系统中分割出一小部分，再根据相应的物理(力学)规律分析邻近部分和这个小部分间的作用(抓住主要作用，略去次要因素，即高等数学中的抓主部，略去高阶无穷小)，这种相互作用在一个短的时间间隔是如何影响物理量u(3)把这种关系用数学算式(方程)表达出来，经化简整理就是所需求的数学物理方程•例1弦的振动弦的振动问题，虽然是一个古典问题，但对于初学者仍然具有一定的启发性设有一根均匀柔软的细弦，平衡时沿直线拉紧，而且除受不随时间而变的张力作用及弦本身的重力外，不受外力影响，下面研究弦的微小横向振动，即假定全部运动出现在一个平面上，而且弦上的点沿垂直于x轴的方向运动(图3-1).图3-1设弦上具有横坐标为X的点，在时刻t时的位置为M,位移NM记作u .显然，在振动过程中位移u 是变量X与t的函数u(x,t) •现在来建立位移u满足的方程我们把弦上点的运动先看作小弧段的运动，然后再考虑小弧段趋于零的极限情况•在弦上任取一弧段MM •，其长为ds,设是弦的线密度，弧段MM •两端所受的张力记作T , T ,现在考虑孤段MM在t时刻的受力情况，用牛顿运动定律，作用于弧段上任一方向上的力的总和等于这段孤的质量乘以该方向上的加速度•在x轴方向弧段受力的总和为_T cos T cos，由于弦只作横向振动，所以T cos - T cos = 0 .( 3.1)如果弦的振动很小，并且在振动过程中弦上的切线倾角也很小，即:-:-0^ 0，则由2 4“ a a _ _cos -1 -2! 4!可知，当:为无穷小量时，cos与1的差量是:的高阶无穷小量，可以略去不计，因此当在u方向弧段受力的总和为-T sin驀"T Sin:「：gds，其中「是单位弧段的质量,-gds是弧段MM •的重力又因当：一0，厂0时上式左边方括号内的部分是由于x产生dx的变化而引起的u(x,t)的改变量，可用微代入(3.1)式，便可近似得到COS：1,cos： 1T T .tg：sin：J i +tg2«tg，皿)，ex=tg：：：u(x dx,t)：u(x，t)lx dx.且小弧段在时刻t沿u方向运动的加速度为2:u(x,t)-t2小弧段的质量为'gds，所以-T sin 二"T sin : 一：gds」ds：：2u(x, t):t2(3.2)T/^-晋」gds -：2u(x,t)吐2dx,:x分代替，即:x£u(x+dx,t) £u(x,t) ◎「初匕,灯［故 $u(x,t )d x—'ex 1 —x,2 2T j u(x,t) j u(x,t)2 ------- --- --------- + 、x 2:xT”g dx 「'5.x 2t 2 g. 般说来，张力较大时弧振动速度变化很快，即u 要比g 大得多，所以又可以把 g 略去•:t经过这样逐步略去一些次要的量，抓住主要的量, 最后得出u( x, t)应近似地满足方程-2:U 2 a 2 一 t x■2u (3.3)T这里的a.式(3.3)称为一维波动方程.如果在振动过程中，弦上另外还受到一个与弦的振动方向平行的外力，且假定单位长 度所受外力的F(x,t)，显然，在这里(3.1)及(3.2)分别为T cos ； " -T cos ：二 0, o2 u Fds-Tsin 二■ T sin ： -，gds …ds —— ct利用上面的推导方法并略去弦本身的重量，可得弦的强迫振动方程为-2 -.2:U 2： U 亨=口 —-1: Xf(x,t),(3.3)'1其中 f(x,t) F(x,t).方程(3.3)与(3.3)'的差别在于(3.3)'的右端多了一个与未知函数 U 无关的项f(x,t),这个项称为自由项，包含有非零自由项的方程称为非齐次方程，自由项恒等于零的方程称为齐次方程.(3.3)为齐次一维波动方程，(3.3)'为非齐次一维波动方程. 例2传输线方程对于直流电或低频的交流电，电路的基尔霍夫定律指出同一支路中电流相等 .但对于较高频率的电流(指频率还没有高到能显著地幅射电磁波的情况)，电路中导线的自感和电容的效应不可忽略，因而同一支路中电流未必相等现考虑一来一往的高频传输线，它被当作具有分布参数的导体(图3-2).在具有分布参数的导体中，电流通过的情况，可以用与电压v来描述，此处i与v都是x,t的函数，记作i(x,t)与v(x,t)，以R，L，C，G分别表示下列参数：R――每一回路单位的趾串联电阻，L——每一回路单位的串联电感，C――每单位长度的分路电容，G――每单位长度的分路电导•根据基尔霍夫第二定律，在长度为x的传输线中，电压降应等于电动势之和，即v -(v ：-v) = R x i L x而v x,ex故上式可写成Ri - L -(3.4)x 汀另外，由基尔霍夫第一定律，流入节点x的电流应等于流出该节点的电流，即=(i ：=i) C x G x v,(3.5) 将方程(3.4)与(3.5)合并，即得i与v 应近似地满足如下方程组这两个方程称为高频传输线方程.X 2=2:V;t 2LC ft 2 , 1；：2V1若令a 2这两个方程与(3.3)完全相同.由此可见，同一个方程可以用来描述不同的LC物理现象，一维波动方程只是波动方程中最简单的情况， 在流体力学、声学及电磁场理论中，还要研究高维的波动方程. 例3电磁场方程从物理学我们知道，电磁场的特性可以用电场强度 E 与磁场强度H 以及电感应强度 D与磁感应强度B 来描述，联系这些量的麦克斯韦(Maxwell)方程组为rotH = J, ct(3.8)cB耐百,(3.9)divB 二 0,(3.10)为了确定函数i 与v ，将方程(3.5)对x 微分，同时在方程(3.4)两端乘以C 后再对t 微分, 并把两个结果相减，即得将(3.4)中的兰代入上式，得CX这就是电流i 近似满足的微分方程，采用类似的方法从( 3.4)与(3.5)中消去i 可得电压V近似满足的方程方程(3.6)或(3.7)称为传输线方程.R, L,C,G 作不同的假定，就可以得到传输线方程的各种特殊形式•例如，在高频传输的情况下，电导与电阻所产生的效应可以忽略不计，也就是 说可令G = R = 0，此时方程(3.6)与(3.7)可简化为；：2i 1；：2i2・X 2iG 兰 一 LC ；x【RC 5, ct :t 2：2i2・i —二 LG X ;t 2 (RC GL)厂 GRi, ct(3.6).2轨LGxV(RC GL)二 GRv, .t :t(3.7)根据不同的具体情况，对参数divD 「.(3.11)将(3.9)与(3.13)代入得而 rot rotH 程为rot rotHrot rotHrotE 「rotE, :t21二 grad div - ' H ,且 divH divB = 0,所以最后得到H 所满足的方同理，若消去 H 即得E 所满足的方程2E如果介质不导电(卞-0)，则上面两个方程简化为2H 1 k 2Hr H , .t 二(3.15)(3.16)(3.15)与(3.佝称为三维波动方程.若将三维波动方程以标量函数的形式表示出来，则可写成其中J 为传导电流的体密度，「为电荷的体密度.这组方程还必须与下述场的物质方程D =eE,(3.13) (1.14)相联立，其中；是介质的介电常数， 」是导磁率，匚为导电率，我们假定介质是均匀而且是各向同性的，此时；，」，；「均为常数.方程(3.8)与(3.9)都同时包含有 E 与H ，从中消去一个变量，就可以得到关于另一个变 H ，在(3,8)式两端求旋度并利用(3.12)与(3.14)得(3.12) 量的微分方程，例如先消去(3.17)2 2 22- 2『召u d u & u= aJ = a — +一 +一2 , ^c x cy 氐丿 其中a_W ，u 是E 或H 的任意一个分量从方程(3.11)与(3.12)还可以推导出静电场的电位所满足的微分方程 •事实上，以(3.12)代入(3.11)得divD = div ；E = ； div E = ■,而电场强度E 与电位u 之间存在关系E = -gradu,所以可得Pdiv(gradu)= -一z或'、u = - — ,( 3.18)z这个非齐次方程称为泊松(Poisson )方程.如果静电场是无源的，即?=0,则(3.18)变成(3.19)这个方程称为拉普拉斯(Lap lace )方程. 例4热传导方程一块热的物体，如果体内每一点的温度不全一样，则在温度较高的点处的热量就要向 温度较低的点处流动，这种现象就是热传导.在工程技术上有许多传热问题都要归结为求物体内温度的分布，现在我们来推导传热过程中温度所满足的微分方程， 与上例类似，我们不是先讨论一点处的温度，而应该先考虑一个区域的温度.为此，在物体中任取一闭曲面S ，它所包围的区域记作 V (图3-3).假设在时刻t ，区域V 内点M(x,y,z)处的温度为u(x,y,z,t)，n 为曲面元素 S 的外法向(从V 内指向V 夕卜).图3-3:t2由传热学可知，在lt,^ t ］时间内，从 S 流入区域V 的热量与时间At ，面积 S ,以及沿曲面的法线方向的温度变化率三者的乘积成正比，即-k(grad u ) S t.其中k 称为物体的热传导系数，当物体为均匀导热体时，k 为常数.于是，从时刻t i 到时刻t 2，通过曲面S 流入区域V 的全部热量为t 2_Q i = tkgrad u1一 S流入的热量使 V 内温度发生了变化，在厶t 时间内区域V 内各点温度从u(x,y,z,t)变化到 u(x,y,z,t+ △ t)，则在△ t 内V 内温度升高所需要的热量为I I I c [u(x, y,z,t t)-u(x,y,z,t)]dVV从而从时刻t i 到时刻t 2，由于温度升高所吸收的热量为tJ ■ u二 c dV dt,飞 _ v ：t其中C 为物体的比热， '为物体的密度，对均匀物体来说，它们都是常数•由于热量守恒，流入的热量应等于物体温度升高所需吸收的热量，即上2 I.t 2. .：uJ t IJJkgradudS dt= J t | JJJ c°〒dV dt.ti.S _t^V醴—此式左端的由面积分中 S 是封闭曲面，可以利用奥 -高公式将它化为三重积分，即11 kgradudS 二 kdiv(gradu)dVSV因此有由于时间间隔't,r 11及区域V 都是任意取的，并且被积函数是连续的，所以 (3.20)式左 右恒等的条件是它们的被积函数恒等，即= Q= k S = t =cnk(grad udS dt.fu(x, y,z,t)dttdV.Q 2二 k 2udV,V：k 2udVti_ Vt2u IF t 」川 cP^dV dt. V ：t(3.20)3(3.21 ).u a :t或二维热传导方程如果我们考虑稳恒温度场，即在热传导方程中物体的温度趋于某种平衡状态，这时温 度u 已与时间t 无关，所以 —=0,此时方程(3.21)就变成拉普拉斯方程(3.19).由此可见稳恒Ct温度场内的温度 U 也满足拉普拉斯方程•在研究气体或液体的扩散过程时，若扩散系数是常数，则所得的扩散方程与热传导方 程完全相同•3. 2初始条件与边界条件上面所讨论的是如何将过程的物理规律用数学式子表达出来•除此以外，我们还需要把具体条件也用数学形式表达出来，这是因为任何一个具体的物理现象都是处在特定条件之下 的.例如弦振动问题，上节所推导出来的方程是一切柔软均匀的弦作微小横向振动的共同规 律，在推导这个方程时没有考虑到弦在初始时刻物状态以及弦所受的约束情况 •如果我们不是泛泛地研究弦的振动，势必就要考虑到弦所具有的特定条件 .因为任何一个具体振动现象总是在某时刻的振动状态和此时刻以前的状态有关，从而就与初始时刻的状态有关•另外，弦的两端所受的约束也会影响弦的振动， 端点所处的物理条件不同会产生不同的影响， 因而弦的振动也不同.所以对弦振动问题来说，除了建立振动方程以外，还需列出它的具体条件 对热传导方程，拉普拉斯方程也是如此.提出的条件应该恰恰能够说明某一具体物理现象的初始状态以及边界上的约束情况，其中a 2 = 2u=a 2；t2 2 2\ d u d u d u r + r + r cy cz厂■方程(3.21)称为三维热传导方程.若物体内有热源，其强度为F (x, y,z)，则相应的热传导方程为:u 2 a :t2 2 2u u u +r 2" 2x : yf(x, y,z,t),如果所考虑的物体是一根细杆 作为特例， 而其中的温度u 只与x,t (或x, y,t )有关,(或一块薄板)，或者即使不是细杆 (或薄板) 则方程(3.21)就变成一维热传导方程:U 2a：t-2；u-2:y 2用以说明系统的初始状态的条件称为初始条件 •用以说明边界上的约束情况的条件称为边界条件•下面具体说明初始条件和边界条件的表达形式，先谈初始条件，对于弦振动问题来说， 初始条件就是弦在开始时刻的位移及速度，若以 「(x),t (X )分别表示初位移和初速度，则 初始条件可以表达为t」= (x)-(X )t =0而对热传导方程来说，初始条件是指在开始时刻物体温度的分布情况，若以 「(M)表 示t=0时物体内任一点 M 处的温度，则热传导方程的初始条件就是u(M,t )tn 「(M).泊松方程与拉普拉斯方程都是描述稳恒状态的，与初始状态无头，所以不提初始条件 再谈边界条件•如果边界条件直接给出了未知函数 u(M ,t)在边界S 上的值，以s 表示边界S 上的动点，则这样的边界条件可表为u(M,t )M ：s 「(St),或简写成(3.24)这种边界条件称为第一类边界条件，其中(s,t)表示在边界S 上给定的已知函数•例如，在杆的导热问题中，若在端点 x = a 处温度保持为常数 u 0，这时在端点x = a 的边界条件为uxn =u 0・若在端点X =a 处温度随时间的变化规律 f (t)为已知，在这点的边界条件为ux 「f (t)・又如在弦振动问题中，若弦的某端点X = a 是固定的，则在该点的位移为零，即Uxn=0・以上都是第一类边界条件的例子 •总之，第一类边界条件直接给出了未知函数u(M ,t)在边界S 上的值但在许多情况下，边界上的物理条件并不能用第一类边界条件来描述 .例如，在杆的导热问题中，若杆的一端x 二a 绝热，那末绝热这个条件就不能直接给出杆的端点处的温度变cu 丄丄化.由于从杆外通过杆端流入杆内的热量为k —— 心S^t (其中A t 为时间间隔，AS 为杆cn(3.22)(3.23)的截面积，n 为杆在端点x=a 处的外法向，若 x=a 是杆的左端点，n 的正向与x 轴正向PuQuQu Qu相反，贝厂■,若x=a 是杆的右端点，则n 的正向与x 轴正向相同，贝厂 '), cnex cn ex所以绝热这个条件可以表达为0( Ik —也S 也 t = 0,cn xy即二 0.x =a若在单位时间内通过 x=a 端单位面积流入杆内的热量是 t 的已知函数f(t)，则这个条件可表示为对于弦振动问题来说， 如果弦在x=a 处是自由的，即沿着位移方向不受外力， 则此时弦在x=a 处沿位移方向的张力(参照 3.1中例1的推导) 为的形式，其中—表示函数沿边界外法向的变化率，这种边界条件称为第二类边界条件•cn除了上述两类边界条件外，有时还会遇到其他形式的边界条件 •例如在杆的导热热问题中，若杆在某个端点 x = a 自由冷却，那末自由冷却这个条件就是cuK 石=H (5 -ux =ax=a ),(其中U i 为周围介质的温度)即(加(k )u + h ———:u 1 h -1州」 X HI H 丿=f(t).x^a总之，有时边界条件必须表达为=0,(3.25)对于有界杆（0乞x 乞I ），若两端都是自由冷却，则在X = I 处，上述条件可表为-U 在x = 0处，这个条件可表为这是由于在单位时间内从周围介质传到杆的 X = a 端单位面积上的热量与介质和杆端的温度差成正比，而在单位时间内通过 考3.1中例4）.x = a 端单位面积传向杆内的热量与x=acn 丿X -i 成正比般地，这种边界条件的形式为(3.26)这样的边界条件称为第三类边界条件•不论哪一种边界条件，如果它的数学表达式中的右端自由项恒为零，则这种边界条件称为齐次的.3.3定解问题的提法前面两节我们推导了三种不同类型的偏微分方程并讨论了与它们相应的初始条件与边界条件的表达方式•由于这些方程中出现的未知函数的偏导数的最高阶都是二阶，而且它们对于未知函数及其各阶偏导数来说都是线性的，所以这种方程称为二阶线性偏微分方程*）1在工程技术上二介线性偏微分方程遇到最多•如果一个函数具有所需要的各阶连续编导数，并且代入某偏微方程中能使该方程变成恒等式，则此函数称为该方程的解•由于每一个物理过程都处在特定的条件之下，所以我们的任务是要求出适合初始条件和边界条件的解•初始条件和边界条件都称为定解条件•求一个偏微方程满足定解条件的解的问题称为定解问题只有初始条件，没有边界条件的定解问题称为始值问题（或柯西问题）；而没有初始条件，只有边界条件的定解问题称为边值问题；既有初始条件也有边界条件的定解问题称为混合问题•一个定解问题提得是否符合实际情况，当然必须靠实际来证实，然而从数学角度来看，可以从三方面加以检验•1）解的存在性，即看所结出来的定解问题是否有解；2）解的唯一性，即看是否只有一个解；3）解的稳定性，即看当定解条件有微小变动时，解是否相应地只有微小的变动，如果*）二阶线性编微分方程可以按它们的二阶导数的系数的代数性质进行分类，在§ 1・1中所推导的波动方程属于双曲型，拉普拉斯（或泊松）方程属于椭圆型，热传导方程属于抛物型，关于二阶线性偏微分方程的分类方法，读者可参阅复旦大学数学系编《数学物理方程》（第二版，上海科学技术岀版社岀版）第一章§5.确定如此，此解便称为稳定的，否则所得的解就无实用价值•因为定解条件通常总是利用实验方法获得的，因而所得到的结果，总有一定的误差，如果因此而解的变动很大，那末这种解显然不能符合客观实际的要求•如果一个定解问题存在唯一且稳定的解，则此问题称为适定的，在以后讨论中我们把着眼点放在讨论定解问题的解法上，而很少讨论它的适定性，这是因为讨论定解问题的适定性往往十分困难，而本书所讨论的定解问题都是古典的，它们的适定性都是经过证明了的.习题一1.长为I的均匀杆，侧面绝缘，一端温度为零，另一端有恒定热流q进入(即单位时间内通过单位截面积流入的热量为q ),杆的初始温度分布是x— x)，试写出相应的定解问题22.长为I的弦两端固定，开始时在x = C受到冲量的作用，试写出相应的定解问题3.有一均匀杆，只要杆中任一小段有纵向位移或速度，必导致邻段的压缩或伸长，这种仲缩传开去，就有纵波沿着杆传播，试推导杆的纵振动方程4.一均匀杆原长为I，一端固定，另一端沿杆的轴线方向被拉长e而静止，突然放手任其振动，试建立振动方程与定解条件.5.若F(z),G(z)是任意二阶可微函数，验证u = F (x at) G(x -at)满足方程-2 -2:-U 2： U2—a 2.t x6.若函数U i(x,t), U2(x, t)JM, U n(x, t)，…均为线性齐次方程-2 -2:u : u-~2 P 2~x :t的解，其中p,q, r只是x, t的函数，而且级数u k(x,t)收敛，并对x,t可以进行两次kz!逐项微分，求证级数u=：£u k(x,t)满足原方程(这个结论叫做线性齐次方程的叠加原理)k m。
数学物理方程第四章(调和)
1
4
(u ( 1 ) 1 u )dS S n rM0M rM0M n
调和函数在区域内任一点的值可以通过积分表达式用这个函数 在区域边界上的值和边界上的法向导数来表示。
数学物理方程
第4章 调和方程
三、调和函数的基本性质
1、调和方程的基本解
k
1
rM
0
M
ln
1
(x x0 )2 ln
P x
Q y
R z
d
Pdydz
Qdxdz
Rdxdy
Pcosn, x Qcosn, y Rcosn, zds 其中n cosn, x,cosn, y,cosn, z 是 在
点 x, y, z 处的外法向量
u
2u x2
2u y 2
2u z 2
0
调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。
二、 拉普拉斯方程边值问题的提法
1 第一边值问题（狄氏问题） 2 第二边值问题（牛曼问题）
u f
u f n
3、狄氏外问题
4、牛曼外问题
数学物理方程
三、泊松方程边值问题
2
z 2
x2
1 y2
z2
3z2 (x2 y2 z2 ) x2 y2 z2 5 2
三式相加，可得
2 1 0, r 0 r
数学物理方程
第4章 调和方程
② 当 r 0时，1 不可导，将 V 取为整个三维空间 r