南邮数理方程4格林函数-上课用

NanjingUniversityofPostsandTelecommunicationsNanjingUniversityofPostsandTelecommunications2第四章格林函数法分离变量法主要适用于求解各种有界问题，而行波法则主要适用于求解各种无界问题，这两种方法所得到的解一般分别为无穷级数和无穷积分的形式。格林函数法给出的解则是有限的积分形式，十分便于理论分析和研究。NanjingUniversityofPostsandTelecommunications3格林函数又称为点源函数或影响函数。顾名思义，它表示一个点源在一定的边界条件和(或)初值条件下所产生的场或影响。由于任意分布的源所产生的场均可看成许许多多点源产生的场的叠加，因此格林函数一旦求出，就可算出任意源的场。格林函数法以统一的方式处理各类数学物理方程，既可以研究常微分方程，又可以研究偏微分方程；既可以研究齐次方程又可以研究非齐次方程；既可以研究有界问题，又可以研究无界问题。它的内容十分丰富，应用极其广泛。这一章，我们主要介绍用格林函数求解拉普拉斯方程的边值问题。NanjingUniversityofPostsandTelecommunications4.1拉普拉斯方程边值问题的提法zyxuu,,,0222222zuyuxu02u设满足拉普拉斯方程描述稳恒状态下的物理过程（如无源静电场的电位分布，稳恒温度场的温度分布等）。通常表示成或不存在初始条件.拉普拉斯方程的解称为调和函数.0uNanjingUniversityofPostsandTelecommunications1）第一边值问题狄利克雷(Direchlet)问题边界条件:|.uf)(0u2）第二边值问题)(0ufnu纽曼(Neumann)问题3）第三边值问题)(0u()uaubfnNanjingUniversityofPostsandTelecommunications4.2格林公式（曲面积分中高斯公式的推论）高斯公式：设是以光滑曲面为边界的有界区域，，，在闭域上连续,在内有一阶连续偏导数，则zyxP,,zyxQ,,zyxR,,cos,cos,cos,PQRdVPnxQnyRnzdSxyz其中为的外法向量。n高斯公式可简记为：SdnadVaNanjingUniversityofPostsandTelecommunications设(,,)uxyz和(,,)vxyz在区域内及其边界上具有连续一阶导数，而在中具有连续二阶导数。xvuzyxP,,yvuzyxQ,,zvuzyxR,,令RQP,,将的表达式代入高斯公式，等式右端dSznzvynyvxnxvu,cos,cos,cosdSnvu而等式左端为：()()PQRuvuvuvdVuvdVdVxyzxxyyzzNanjingUniversityofPostsandTelecommunications()()uvdVuvdV因此，高斯定理可改写为：()()PQRuvuvuvdVuvdVdVxyzxxyyzz()()vuvdVuvdVudSn()-()vuvdVudSuvdVn或----第一格林公式NanjingUniversityofPostsandTelecommunications()-()vuvdVudSuvdVn把中u,v交换位置，可得：vudVdSnuvvudV两式相减可得：()vuuvdSnn()uvvudV----第二格林公式NanjingUniversityofPostsandTelecommunications调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。1)拉普拉斯方程的钮曼问题有解的必要条件设是在以为边界的区域内的调和函数,在上有一阶连续偏导数,则在第二格林公式中取为上述调和函数,,则有.所以钮曼内问题()有解的必要条件为函数满足uu1v0dSnufnu|f0fdS事实上,这也是钮曼内问题有解的充分条件.()()SVvuuvdSuvvudVnnNanjingUniversityofPostsandTelecommunications2)拉普拉斯方程解的唯一性问题设是定解问题的两个解，则它们的差必是原问题满足零边界条件的解。对于狄利克雷问题，满足21,uu21uuvv0,|0vvv对于钮曼问题,满足v0,|0vvvnNanjingUniversityofPostsandTelecommunications在第一格林公式中取,由是调和函数，可得21uuvuvdSnvv0vvdV0vvdSn在两种边界条件下，都有20.vdV所以故在内必有,即：0v0zvyvxvCvC可得，其中为常数.()-()vuvdVudSuvdVnNanjingUniversityofPostsandTelecommunications对于狄利克雷问题,由于,故从而.0|v0C0v结论狄利克雷问题在原定解问题中的解是唯一确定的；钮曼问题的解在相差一个常数下也是唯一确定的.3)调和函数的积分表达式所谓调和函数的积分表达式,是指用调和函数在区域边界上的值及其在区域边界上的法向导数沿的积分来表达调和函数在内任一点的值.NanjingUniversityofPostsandTelecommunications设是内一固定点,下面求调和函数在这一点的值.0000,,zyxM为此构造一个辅助函数20202011zzyyxxrv可以证明函数除点外处处满足拉普拉斯方程,它称为三维拉普拉斯方程的基本解.（课后作业）1vr0M为了利用格林公式，我们在内挖去的球形邻域（球半径为），是其球面。0MK,0NanjingUniversityofPostsandTelecommunications在区域内及其边界上，是任意可导的。Krv1在第二格林公式中,取为调和函数,并假定它在上有一阶连续偏导数,而取,在区域上应用第二格林公式得：urv1K11KuudVrr11urudSnrn221/1/11rrnrr在球面上,NanjingUniversityofPostsandTelecommunications2221/1144rudSudSuun因此同理可得nudSnudSnur41104411nuudSnurrnu我们可得令,则：,000limMuu04lim0nu0001114MMMMuMuMuMdSnrrn于是NanjingUniversityofPostsandTelecommunications调和函数在区域内任一点的值可以通过积分表达式用这个函数在区域边界上的值和边界上的法向导数来表示。0001114MMMMuMuMuMdSnrrn0222000()()()MMrxxyyzzNanjingUniversityofPostsandTelecommunications(4)设),,(zyxu是区域内的调和函数，它在上有一阶连续偏导数，则,0dSnu其中n,的外法线方向。是证明只要在第二格林公式中取即证。1v注：此性质表明调和函数的法向导数沿区域边界的积分为零。对稳定的温度场，流入和流出物体界面的热量相等，否则就不能保持热的动态平衡，而使温度场不稳定。（5）(极值原理)设函数),,(zyxu在区域内调和，它在上连续且不为常数，则它的最大值与最小值只能在边界上达到。NanjingUniversityofPostsandTelecommunications6)平均值公式设函数在某区域内是调和函数,是内任一点,表示以为中心,为半径且完全落在内的球面,则有：MuaK0Ma0214aKuMudSa0M将调和函数的积分表达式应用于球面aK则在球面上，aK11,ra2111()()nrrra110aaKKuudSdSrnan0214aKuMudSa则有：NanjingUniversityofPostsandTelecommunications4.2格林函数NanjingUniversityofPostsandTelecommunications由于调和函数有积分表示:dSnMurrnMuMuMMMM)(1)1()(41)(000又因为Dirichlet边值问题fuxu,,0的解唯一,故希望将此问题的解用积分表示出来。但由于在积分表达示中，nuu在边界上的值虽然已知,而在边界上的值却不知道.那么只有想办法去掉|un为此，引入格林函数的概念。(4.2.1)NanjingUniversityofPostsandTelecommunications224.2.1格林函数的定义设在内有vuvu,;0,0在上有一阶连续偏导数，则由第二格林公式有:0)(dSnuvnvu(4.2.2)将（4.2.1）和(4.2.2)两式加起来：dSnuvrrnnvuMuMMMM)41()]1(41[)(000(4.2.3)选择调和函数v满足041MMrv,于是有：dSvrnuMuMM)41()(00(4.2.4)NanjingUniversityofPostsandTelecommunications23记vrMMGMM041),(0(4.2.5)则有dsnGuMu)(0(4.2.6)称),(0MMG为Laplace方程的格林函数。若),(0MMG上有一阶连续偏导数，则当Dirichlet问题且在),,(),,(,0zyxfuzyxu上具有一阶连续偏导数的解存在时，解可以表示为在dSnGzyxfMu),,()(0(4.2.7)存在NanjingUniversityofPostsandTelecommunications24对Poisson方程的Dirichlet问题上存在具有一阶连续偏导数的解，则解可以如果在fuzyxFu),,(,表示为FGdVdSnGfMu)(0由此可见,求解Dirichlet问题,关键是求Green函数(4.2.5),其中v满足一个特殊的Dir