第四章格林函数法
合集下载
第四章拉氏方程的格林函数法.docx
第四章拉氏方程的格林函数法.docx第四章拉氏方程的格林函数法前面儿章，介绍了儿种求解PDE定解问题的方法：分离变量法、行波法、积分变换法。
?本章介绍令一种求解拉氏方程的格林函数法。
首先来看一下我们要研究的定解问题是怎么捉出的。
§4.1拉氏方程边值问题的提法在第一章中，我们知道，对于无源的稳恒热传导问题满足拉氏方程，它的边值问题一般有三种提法。
研究最多的就是前面两种。
1)第一边值问题边界条件为:心要求的解心C2(Q)AC°(Q),既比在区域Q上连续，在Q上有二阶连续导数,满足拉氏方程且在边界上与/吻合。
Q = Q + 「为边界；称第一-边值问题为狄利克莱(Dirichlet)问题，简称狄氏问题。
通常称拉氏方程的连续解，也就是说，具有二阶连续偏导数并满足拉氏方程的连续函数为调和函数。
2)第二边值问题边界条件为:単=/,on r要求的解ue C2(Q)nC*(Q),既u在区域豆上有一阶连续导数，在Q 上有二阶连续导数，满足拉氏方程且在边界上满足上边界条件。
称第二边值问题为牛曼(Neumarm)问题，简称牛氏问题。
前面两种边值问题都是在Q内求解拉氏方程，故称此类方程为内问题。
另外, 冇这样一类问题，如已知某区域边界上的温度，要求该区域Z外的温度分布情况, 这就归结为在区域Q外求解拉氏问题，称这样的问题为外问题。
注：对于外问题來说，求解通常都是在无界区域上，这时需不需要对解加些限制条件呢？看下面一例了。
Aw = 0, r > 1, r = Vx2 + r2 + z2易知u = Vu = \!r都是上定解问题的解，这就出现了解的不唯一性，为了保证解的唯一性，通常我们要加一些限制条件，三维问题时limw = 0厂T8二维问题通常假定解冇界。
3）狄氏外问题（略）4）牛氏外问题（略）§ 4.2格林公式及其应用一、格林公式的推导为建立拉氏方程解的积分公式,我们先推导格林函数，它由曲面积分的Guass 公式直接导出。
4第四章格林函数法
，于是有 除在 M 0 点外处处满足三维Laplace方程 u0 内调和,则 上有一阶连续偏导数,且在 定理:若函数 u 在 调和函数在区域内任一点的值可以通过积分表达式用这个函 数在区域边界上的值和边界上的法向导数来表示。
u ( M ) 1[ 1 1 u ( M ) u ( M ) ( ) ] dS 0 4 n r r n MM MM 0 0
3
P Q R { P , Q , R } n dS ( ) dV 由高斯公式 x y z v v v
2019/2/12
4.1.3 调和函数的积分表达式
由Green公式可导出调和函数的积分表示。由于函数：
1 1 2 2 2 r ( x x ) ( y y ) ( z z ) MM 0 0 0 0
为二维Laplace方程的基本解.
其通解为： 为任意常数)。 V ( r ) c ln r c , ( r 0 , c , c 1 2 1 2
4.1.2 格林公式
令P u , Q u ,R u ,则得到格林第一公式： x y z u v u v u v v u vdV ( ) dV u dS x x y y z z n u v u v u v u v udV ( ) dV vdS x x y y z z n 将以上两公式相减，得到格林第二公式： v u ( u v v u ) dV ( u v ) dS n n 调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。
1
2019/2/12
4.1.1 拉普拉斯方程的基本解 对拉普拉斯方程 , 其球坐标形式为： u u u u 0 xx yy zz
拉普拉斯方程的格林函数法
注：利用调和函数的极值原理，可证狄氏问题在 C2() C0()
内的解是唯一的。
整理课件
17
§3 格林函数
为 什 么 引 入 格 林 函 数 ？
调 和 函 数 的 积 分 表 达 式 为
1 1 1 u
u(M0)4u(M)n(rM M0)rM M0
dS n
对于狄氏问题或牛曼问题，利用上面公式都不能直接得到想要问
( P x Q y R z)d V P d y d z Q d z d x R d x d y
其 中 取 外 侧 位 正 向 .
由两类曲面积分之间的关系得高斯公式的另一种形式：
整理课件
7
(P xQ yR z)dV (Pcos(n,x)Qcos(n,y)Rcos(n,z))dS.
_ u n 1 r 1 r u n d S 4 u 4 u n 0
当 0 时 ， 有 l i m 0 u u ( M 0 ) ,(u 连 续 ) ____
lim 4u0( u一 阶 连 续 可 微 ， u有 界 ）
0 n
定理：若格林函数G(M,M0)存在，且G(M,M0)C1()，则狄氏
问题2uu0,f(M in)的解（存在的话）可表示为
这 时 需 不 需 要 对 解 加 些 限 制 条 件 呢 ？ 看 下 面 一 例 子 。
u0,r1,
u 1 r1
其 中 rx2y2z2
易 知
u1,
u1/r
都 是 上 述 定 解 问 题 的 解 ， 即 解 不 唯 一 .为 了 保 证 解 的 唯 一 性 ，
通 常 我 们 要 加 一 些 限 制 条 件 .
取 u 为 调 和 函 数 ， 并 假 定 其 在 上 有 一 阶 连 续 偏 导 数 ， 取 v 1/r
内的解是唯一的。
整理课件
17
§3 格林函数
为 什 么 引 入 格 林 函 数 ？
调 和 函 数 的 积 分 表 达 式 为
1 1 1 u
u(M0)4u(M)n(rM M0)rM M0
dS n
对于狄氏问题或牛曼问题，利用上面公式都不能直接得到想要问
( P x Q y R z)d V P d y d z Q d z d x R d x d y
其 中 取 外 侧 位 正 向 .
由两类曲面积分之间的关系得高斯公式的另一种形式：
整理课件
7
(P xQ yR z)dV (Pcos(n,x)Qcos(n,y)Rcos(n,z))dS.
_ u n 1 r 1 r u n d S 4 u 4 u n 0
当 0 时 ， 有 l i m 0 u u ( M 0 ) ,(u 连 续 ) ____
lim 4u0( u一 阶 连 续 可 微 ， u有 界 ）
0 n
定理：若格林函数G(M,M0)存在，且G(M,M0)C1()，则狄氏
问题2uu0,f(M in)的解（存在的话）可表示为
这 时 需 不 需 要 对 解 加 些 限 制 条 件 呢 ？ 看 下 面 一 例 子 。
u0,r1,
u 1 r1
其 中 rx2y2z2
易 知
u1,
u1/r
都 是 上 述 定 解 问 题 的 解 ， 即 解 不 唯 一 .为 了 保 证 解 的 唯 一 性 ，
通 常 我 们 要 加 一 些 限 制 条 件 .
取 u 为 调 和 函 数 ， 并 假 定 其 在 上 有 一 阶 连 续 偏 导 数 ， 取 v 1/r
第四章_拉普拉斯方程的格林函数法
Laplace方程的连续解，即具有二阶连续偏导数并满足Laplace 方程的连续函数. (1.u ∈ C 2 (Ω) I C 0 (Ω) 2. ∇ 2u = 0)
边值问题的提法： 二. 边值问题的提法：
1）第一边值问题(Dirichlet问题/狄氏问题) ∇ 2u = 0, in Ω uΓ = f 数学解释： 在Ω内寻求一个调和函数u, 它在边界Γ上与已知
连续函数f 吻合，即u Γ = f .
2）第二边值问题( Neumann 问题/牛氏问题 ) ∇ 2 u = 0, in Ω (其中 n是 Γ 的外法向量, f 是 连续函数 ) ∂u = f ∂n Γ
数学解释： Ω内寻求一个调和函数，它在闭区域Ω上有一阶 在 连续偏导数,即u ∈ C 2 (Ω) I C1 (Ω)，且在边界上满足边界条件。 ，且在边界上满足边界条件。
∂u ∫∫ u (M )dS + 4π a ∫∫ ∂n dS Γ Γ 1
∫∫ u (M )dS
Γ
（4）Laplace方程解的唯一性问题 ） 方程解的唯一性问题
定理：狄氏问题在C 2 (Ω) I C 1 (Ω)内解唯一，牛曼问题除相差一个
常数外解也是唯一确定的。
证明：
设u1 , u2为上述两类问题的解，则它们的差v = u1 − u2必是原问题的
第四章 拉普拉斯方程的格林函数法
第一节 第二节 第三节 第四节 拉普拉斯方程边值问题的提法 格林公式 格林函数 两种特殊区域的格林函数及狄氏 问题的解
格 林 函 数 法
格林函数:又称点源影响函数，是数学物理中的一个重要概念 格林函数 又称点源影响函数， 又称点源影响函数
格林函数代表一个点源在一定的边界条件和初始条件下所产生的 格林函数代表一个点源在一定的边界条件和初始条件下所产生的 知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 场,知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 知道了点源的场就可以用迭加的方法计算出任意源所产生的场
边值问题的提法： 二. 边值问题的提法：
1）第一边值问题(Dirichlet问题/狄氏问题) ∇ 2u = 0, in Ω uΓ = f 数学解释： 在Ω内寻求一个调和函数u, 它在边界Γ上与已知
连续函数f 吻合，即u Γ = f .
2）第二边值问题( Neumann 问题/牛氏问题 ) ∇ 2 u = 0, in Ω (其中 n是 Γ 的外法向量, f 是 连续函数 ) ∂u = f ∂n Γ
数学解释： Ω内寻求一个调和函数，它在闭区域Ω上有一阶 在 连续偏导数,即u ∈ C 2 (Ω) I C1 (Ω)，且在边界上满足边界条件。 ，且在边界上满足边界条件。
∂u ∫∫ u (M )dS + 4π a ∫∫ ∂n dS Γ Γ 1
∫∫ u (M )dS
Γ
（4）Laplace方程解的唯一性问题 ） 方程解的唯一性问题
定理：狄氏问题在C 2 (Ω) I C 1 (Ω)内解唯一，牛曼问题除相差一个
常数外解也是唯一确定的。
证明：
设u1 , u2为上述两类问题的解，则它们的差v = u1 − u2必是原问题的
第四章 拉普拉斯方程的格林函数法
第一节 第二节 第三节 第四节 拉普拉斯方程边值问题的提法 格林公式 格林函数 两种特殊区域的格林函数及狄氏 问题的解
格 林 函 数 法
格林函数:又称点源影响函数，是数学物理中的一个重要概念 格林函数 又称点源影响函数， 又称点源影响函数
格林函数代表一个点源在一定的边界条件和初始条件下所产生的 格林函数代表一个点源在一定的边界条件和初始条件下所产生的 知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 场,知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 知道了点源的场就可以用迭加的方法计算出任意源所产生的场
4格林函数法
u( r ) f ( r0 ) dr0 4 | r r0 |
那么，如何求解某边界条件下的泊松方程呢？
格林函数法求解拉普拉斯方程的边值问题
边值问题与无界空间的问题不同，要受到边界的影响， 边值问题的解 u(r) 与格林函数 G 的关系就更复杂了，需要 用到格林第二公式。 格林第二公式 设函数 u(r)， v(r) 在区域 Ω 直到其边界 Γ 上具有连续一阶 微商，而在 Ω 中有连续二阶微商，有格林第二公式
若知道一个点源在一定的边界条件和(或)初值 条件下所产生的场（称为格林函数），就可用叠加 的方法计算出任意源产生的场。这就是格林函数法 的基本思想。故也称点源法。
δ – 函数（狄拉克函数）
用来描述物理中集中分布的量，如点电荷、点热源、质点、单 位脉冲等，这是通常的函数概念不能描述的。
0 x x0 ( x x0 ) x x0 0 x0 ( a , b ) a ( x x0 )dx 1 x0 (a , b)
r r0 o
2
R (x0, y0, z0)
G G x x 0 , x R R
2G 2G x x0 G 1 ( x x0 )2 2 2 3 x R R R R R
2G 2G 同样计算可 2 和 ，于是 2 y z 2G R 2 G 3 R 2 2G 2 G 2G 3 2 2 R R R 2 R R 0 R R R
于是有积分形式的解
拉普拉斯方程的格林函数
根据格林函数法的基本思想，先求解方程
(x, y, z) r r0 o
G ( r , r0 ) ( r r0 )
2
那么，如何求解某边界条件下的泊松方程呢？
格林函数法求解拉普拉斯方程的边值问题
边值问题与无界空间的问题不同，要受到边界的影响， 边值问题的解 u(r) 与格林函数 G 的关系就更复杂了，需要 用到格林第二公式。 格林第二公式 设函数 u(r)， v(r) 在区域 Ω 直到其边界 Γ 上具有连续一阶 微商，而在 Ω 中有连续二阶微商，有格林第二公式
若知道一个点源在一定的边界条件和(或)初值 条件下所产生的场（称为格林函数），就可用叠加 的方法计算出任意源产生的场。这就是格林函数法 的基本思想。故也称点源法。
δ – 函数（狄拉克函数）
用来描述物理中集中分布的量，如点电荷、点热源、质点、单 位脉冲等，这是通常的函数概念不能描述的。
0 x x0 ( x x0 ) x x0 0 x0 ( a , b ) a ( x x0 )dx 1 x0 (a , b)
r r0 o
2
R (x0, y0, z0)
G G x x 0 , x R R
2G 2G x x0 G 1 ( x x0 )2 2 2 3 x R R R R R
2G 2G 同样计算可 2 和 ，于是 2 y z 2G R 2 G 3 R 2 2G 2 G 2G 3 2 2 R R R 2 R R 0 R R R
于是有积分形式的解
拉普拉斯方程的格林函数
根据格林函数法的基本思想，先求解方程
(x, y, z) r r0 o
G ( r , r0 ) ( r r0 )
2
数学物理方程课件第四章拉普拉斯方程的格林函数法
r M 0 M
M 1
1
4 xx02 y y02 zz02
解：
1
4 xx02 y y02 zz02
u(M 0)G (M n,M 0)f(M )dS G(M z,M0)|z0 f(x,y)dS
数学物理方程与特殊函数
第4章格林函数法
1
1
G ( M , M 0 ) 4 x x 0 2 y y 0 2 z z 0 2 4 x x 0 2 y y 0 2 z z 0 2
调和函数的积分表达式
k
拉 普l1r拉n 斯1
1 方x程2的基y本2 解z
ln 1
2
r
x2 y2
三维 二维
1 1 1 u
u (M 0)4 S(u n(r)r n)d S
调和函数在区域内任一点的值可以通过积分表达式用这个
函数在区域边界上的值和边界上的法向导数来表示。
2 牛曼内问题有解的必要条件
V (u 2 v v 2 u )d V S (u n v v u n )d S
一 拉普拉斯方程边值问
题 的 1提 第法一边值问题（狄氏问题）
第四章
拉普 u f
2 第二边值问题（牛曼问题）
拉斯方程的格 u f 林函数法 n
3 内问题与外问题
4 调和函数：具有二阶偏导数并且满足拉普拉斯方程 的连续函数。
二 格林公式及其结论
V (u 2 v )d V S u n vd S V u v d V 格V 林(u 公 2 式v 的v 结 2 论u )：d V S (u n v v u n )d S
半空间的格林函数
1 1 1
G(M,M0)4rM
r M 0 M
M 1
M0q d
格林函数法
两种边值问题: 两种边值问题:
第一边值问题
u |Γ = f .
这类问题也叫做狄利克雷 问题。 这类问题也叫做狄利克雷(Dirichlet)问题。 狄利克雷 问题
拉普拉斯方程的连续解，也叫调和函数， 拉普拉斯方程的连续解，也叫调和函数，所以 调和函数 狄利克雷问题也可以叙述为： 狄利克雷问题也可以叙述为：在区域 Ω 内找 一个调和函数, 上的值已知。 一个调和函数 它在边界 Γ 上的值已知。 第二边值问题 在光滑的闭曲面 Γ 上给出连续函数 f，寻找函数 ， u(x,y,z) ：在 Γ 的内部 Ω 是调和函数 在 ( Ω + Γ ) 是调和函数, 上连续， 上连续，在 Γ 上任一点法向导数存在并且等于 已知函数 f ，即： ∂u =f ∂n Γ 这类问题也叫做纽曼 纽曼(Neumann)问题。 问题。 这类问题也叫做纽曼 问题
在球面坐标下, 拉普拉斯方程为: 在球面坐标下, 拉普拉斯方程为:
1 ∂ 2 ∂u 1 ∂ ∂u 1 ∂2u =0 r + 2 sin θ + 2 2 2 2 r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ
球对称解 u=u(x,y,z)在以原点为中心的同一球面的 在以原点为中心的同一球面的 值为常数。 的函数:u=u(r)。 值为常数。u 仅为半径 r 的函数 。
(
)
两式相减, 两式相减 得
2 2
第二格林公式
∂v ∂u ∫∫∫ (u∇ v − v∇ u)dV = ∫∫ (u ∂n − v ∂n )dS Ω Γ
利用格林公式, 可以得到调和函数的一些性质： 利用格林公式, 可以得到调和函数的一些性质： 1) 纽曼内问题有解的必要条件 内的调和函数, 设 u 是以 Γ 为边界的区域 Ω 内的调和函数, 在 上有一阶连续偏导数, Ω + Γ上有一阶连续偏导数, 则在第二格林公式 为上述调和函数 调和函数, 中取 u 为上述调和函数,取 v ≡ 1, 有
数学物理方程第四章_格林函数
1 ⎧ ⎪∆G (r , r0 ) = − δ (r − r0 ) ε ⎨ ⎪G Γ = 0 ⎩
(4.3.7) (4.3.8)
以 G (r , r0 ) 乘式 (4.3.5), u (r ) 乘式 (4.3.7), 二式相减后在 Ω 上对 r 积分 ,以 dr 表示 r 点处的体积微元,有
∫
Ω
(G∆u − u∆G )dr = −
第 4 章 格林函数
在这一章里,我们介绍数学物理方程中另外一种常用的方法—格林函数法.从物理上看, 一个数学物理方程是表示一种特定的“场”和产生这种场的“源”之间的关系.例如,热传导 方程表示温度场和热源之间的关系,泊松方程表示静电场和电荷分布的关系,等等.这样,当源 被分解成很多点源的叠加时,如果能设法知道点源产生的场,利用叠加原理,我们可以求出同 样边界条件下任意源的场,这种求解数学物理方程的方法就叫格林函数法.而点源产生的场就 叫做格林函数. 4.1
⎧0, T ( x) = ⎨ ⎩∞,
x≠0 x=0
且
∫Байду номын сангаас
所以有
+∞
−∞
cρT ( x)dx = Q
T ( x) =
Q δ ( x) cρ
通过以上两个例题,我们对 δ ( x) 有了进一步的认识.如果将坐标平移 x0 ,即集中量 出现在点 x = x 0 处,则有
δ ( x − x0 ) = ⎨
且
⎧0, ⎩∞,
∫
= ∫ (u∆v)dΩ + ∫ gradu ⋅ gradvdΩ
Ω Ω
=∫u
Γ
∂v dS ∂n
或表示为
∫
Ω
(u∆v)dΩ = ∫ u
Γ
∂v dS − ∫ gradu ⋅ gradvdΩ Ω ∂n
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
HUST 数学物理方程与特殊函数
第4章格林函数法
思考：Laplace方程Neumann问题有解的必要条件是什么？
u u 0, | f . n
性质2 (平均值定理)
u n dS f dS 0.
设函数 u( M ) 在区域 内调和，
M0 是
的球面，此球完全落在区域 的内部，则有
4 r
2u(M ) r
0 F ()
u(M ) F ()
2
1 u(M ) F ()dV0 4 rMM0
HUST 数学物理方程与特殊函数
第4章格林函数法
4.2.1 格林函数的定义
设在 内有 u 0, v 0; u, v 在 上有一阶连续
1 1 u ( M ) u(M ) ( ) dS (4.2.1) n rMM 0 rMM 0 n
将此问题的解用积分表示出来。但由于在积分表达示中，u
u 在边界上的值虽然已知,而 在边界上的值却不知道.那么, n u | 的值呢？显然这是行不通的， 能否作为边界条件加上 n u | 因为,此时的解已经是唯一的了.那么只有想办法去掉 n 为此，引入格林函数的概念。
n 次后， 为中心 , 半径小于 d 的球
， …,
k n 内 , 因而 u ( N ) u ( M n ) u ( M 1 ) , 由 N 的 任意性 , 就得到整个 上有 u ( N ) u ( M 1 ) , 这与 u 不为
常数矛盾.
HUST 数学物理方程与特殊函数
第4章格林函数法
4.1.4 调和函数的性质
性质1. 设 u ( x, y, z ) 是区域 内的调和函数，它在 u 上有一阶连续偏导数，则 dS 0, 其中 , n n 是 的外法线方向。
证明 只要在Green公式中取 v 1 即证。 注：此性质表明调和函数的法向导数沿区域边界的积分为零。 对稳定的温度场，流入和流出物体界面的热量相等，否则就 不能保持热的动态平衡，而使温度场不稳定。
偏导数，则由格林第二公式有
v u (u n v n )dS 0
将（4.2.1）和 (4.2.2)两式加起来：
(4.2.2)
1 1 u v 1 u(M 0 ) u[ ( )] ( v) dS (4.2.3) 4 rMM 0 n n 4 n rMM 0
HUST 数学物理方程与特殊函数
如图4.1 , 以
u(M ) u(M 1 ) .
M1
第4章格林函数法
证明： （ 反证法）假设 u 在 内某点 M 1 达到最大值， 为中心，任意长 R 为半径作球 k R ,使 它完全落在区域 中，记 kR 的球面为 S R , 则在 SR 上有 这是因为，若 M ，使u ( M ) u ( M 1 ) ,则由函数的 连续性，必可找到此点在球面 S R上的一个邻域，在此 邻域中，也有 u ( M )
将以上两公式相减，得到格林第二公式：
v u (uv vu )dV (u n v n )dS
调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。
HUST 数学物理方程与特殊函数
第4章格林函数法
4.1.3 调和函数的积分表达式 由Green公式可导出调和函数的积分表示。由于函数： 1 1 rMM ( x x0 ) 2 ( y y 0 ) 2 ( z z 0 ) 2
若取 c1 1, c2 0,则得到特解 V0 ( r ) ,称此解为三维Laplace
r
c1 c2 , (r 0, c1 , c2 为任意常数)。 r 1
方程的基本解,它在研究三维拉普拉斯方程中起着重要的作用.
HUST 数学物理方程与特殊函数
第4章格林函数法
对二维拉普拉斯方程 u uxx u yy 0 ,其极坐标形式为:
若取 c1 1, c2
1 0 , 则得到特解 V0 (r ) ln , 称此解为二维 r
Laplace方程的基本解.
HUST 数学物理方程与特殊函数
第4章格林函数法
4.1.2 格林公式
由高斯公式 P Q R dV P cos(n, x) Q cos(n, y) R cos(n, z) d S x y z v v v 令P u , Q u , R u ,则在
a
a
上式称为调和函数的球面平均值公式。 性质3 (极值原理)
1 1 1 ，所以 上有 n ( r ) r 2 a 2 1 u(M 0 ) udS . 2 a 4a
1 u 1 dS r n a
a
u dS 0 n
调和函数在区域内任一点的值可以通过积分表达式用这个 函数在区域边界上的值和边界上的法向导数来表示。
HUST 数学物理方程与特殊函数
第4章格林函数法
若函数
内满足Poisson方程 u
1 u(M 0 ) 4
u 在 上有一阶连续偏导数，且在
F ，则同样有
1 1 u ( M ) 1 u(M ) ( ) dS n rMM 0 rMM 0 n 4 F (M ) rMM dV 0
设函数 u ( x, y, z ) 在区域 内调和，
它在 上连续且不为常数，则它的最大值与最小值
只能在边界上达到。 推论1 设在 内有 u 0, v 0; u, v 在
上连续且在边界 上有 u v，则在 内有 u v.
u 0, ( x, y, z ) 推论2 Dirichlet问题 的解是唯一的。 u f ( x, y, z )
且在 上有一阶连续偏导数，则当Dirichlet问题
G u ( M 0 ) u ds n
(4.2.6)
u 0 , ( x, y, z ) u f ( x, y, z ) 在 上具有一阶连续偏导数的解存在时，解可以表示为
1 u(M 0 ) udS 2 a 4a 证明： 由调和函数的积分表示：
内任意一点，若 是以 M 为中心，a为半径 a 0
1 u(M 0 ) 4
及由性质1，有
1 1 u u ( ) dS a r n n r
HUST 数学物理方程与特殊函数
M 2 K1 M1 K2 M3 S1 S2
Kn N Mn Sn
l
图4.1
HUST 数学物理方程与特殊函数
第4章格林函数法
4.2 格林函数
由于调和函数有积分表示:
1 u(M 0 ) 4
u 0, x , 又因为Dirichlet边值问题 的解唯一,故希望 u f
选择调和函数v满足 v
1 4 rMM 0
,于是有：
1 u(M 0 ) u ( v)dS n 4r MM 0
(4.2.4)
HUST 数学物理方程与特殊函数
第4章格林函数法
记
1 G( M , M 0 ) v 4 rMM 0
(4.2.5)
则有
称 G( M , M 0 ) 为Laplace方程的格林函数。若G( M , M 0 ) 存在
u ( M ) u ( M 1 ) 。设 M 2 是 K 1 的球面 S1 与折线 L 的交点，
则 u ( M 2 ) u ( M 1 ) 。以 M 2 为中心，以小于 d 的数为半径 在 内作球 k 2 ，在 k 2上 u ( M ) u ( M 2 ) u ( M 1 ) 点 N 一定包含在以某点 M n
0
除在 M 0 点外处处满足三维Laplace方程u
定理：若函数 内调和，则
0，于是有
u 在 上有一阶连续偏导数，且在
1 1 u ( M ) u(M ) ( ) dS n rMM 0 rMM 0 n
1 u(M 0 ) 4
HUST 数学物理方程与特殊函数
第4章格林函数法
格林函数又称为点源函数或影响函数。顾名思 义，它表示一个点源在一定的边界条件和(或)初值条 件下所产生的场或影响。由于任意分布的源所产生的 场均可看成许许多多点源产生的场的叠加，因此格林 函数一旦求出，就可算出任意源的场。格林函数法以 统一的方式处理各类数学物理方程，既可以研究常微
2u 1 u 1 2u 2 0 2 2 r r r r
(4.1.2)
求方程(4.1.2)的径向对称解 u V (r ) (即与 无关的解) ,则有：
其通解为： V (r ) c1 ln r c2 , (r
d 2V 1 dV 0 2 dr r dr
0, c1 , c2 为任意常数)。
分方程，又可以研究偏微分方程；既可以研究齐次方 程又可以研究非齐次方程；既可以研究有界问题，又 可以研究无界问题。它的内容十分丰富，应用极其广 泛。这一章，我们主要介绍用格林函数求解拉普拉斯
方程的边值问题。
HUST 数学物理方程与特殊函数
第4章格林函数法
4.1 格林公式及其应用
4.1.1 基本解 对拉普拉斯方程 u uxx u yy uzz 0 , 其球坐标形式为：
1 4R 2
u ( M 1 ) 。因此有
1 4R 2
S
u ( M ) ds
R
S u ( M 1 )ds u ( M 1 )
R
这与平均值性质矛盾。由于 R 的任意性，则在球 k R 中 恒有 u u ( M 1 ) .任取 N ,在 中作连结 M 1 ， N 两 点的折线 L ，记L 到 的边界 的最小距离为 d ，以M 1 为中心，小于d 的数为半径在 内作求 K1 ，则在 K1上