复旦数学物理方法讲义14格林函数--偏微分方程解的积分表示

1第十四章格林函数--偏微分方程解的积分表示解偏微分方程主要有两种方法：A:数理方法中的分离变量法：正交的多项式或无穷级数解，但需要齐次边界条件。B:理论物理中的Green函数方法：既是简单的有理形式解，又允许任意的边界条件！1，Green函数（GF）的意义：物理上：点源产生的场（函数）在时空中的分布。特别是它在空间是源函数；在时空是传播函数。(Seebelow)数学上:具有点源的偏微分方程在齐次边界条件或者无界区域、初值条件下的解。2，GF的分类：边界值GF：(,')Grr即源函数；初始值GF：(,;',')Grtrt即传播函数。3，Green函数的性质：1）对称性：(,')(',)GrrGrr，它与定解问题相关，即与厄米性相关。(See4below)2）时间传播函数没有对称性：(,;',')(',';,)GrtrtGrtrt.（因果律引起）3）存在的必要条件：设方程2()(,')(')Grrrr，若λ是对应齐次方程的本征值，即2和附加齐次边界条件，则(,')Grr不存在。这是因为既有点源：(')rr矛盾于又无流：|0.nG本征值问题存在，但是没有激发，物理上自相矛盾！平面波(),ikxatAe球面波1()ikxatAre和柱面波1/2()ikatAe均是LaplaceEquation的解，但不是PossionEquation的解。球、柱面波分别来自于1x时（散射问题）渐近行为：(1,2)[(1/2)/2](1,2)[(1)/2]21(),().ixmixlmlxxHxehxe4，Green函数的边值条件：选取边值条件具有人为性，但要求简单并保证算子的厄米性。1）齐次边值条件：()|0GGn矛盾于上述3,3)，详见下面14.1.2末。2）有解性（解收敛）：|0rG--基本解。5，Green函数的用途：偏微分方程的积分解法：1）求(,')Grr,2）利用迭加原理给出物理问题()ur的积分形式解；Green函数的奇点与元激发的能量和寿命有关。6，Green函数的求法：1)特殊方法：2(')1/|'|.GrrGrr无界2）本征函数展开法：相应算子在同一边界下的本征函数作为基矢。3）方程齐次化方法：将非齐次项变成边值条件和初值条件。4）积分变换法：LaplaceTransforms，FourierTransforms.5)形式解：算子运算。214.1Green函数与偏微分方程1，定义：Green函数(源函数，影响函数，传播函数，传播子)数学上，含点源的偏微分方程在一定的边界条件或者初始条件下的解；物理上，点源在一定物理条件下产生的场。这种解（场）在时空中的分布与传播。例1，无界空间PossionEquation:224(),(,')4('),|0.|0.1(,')=,()(,')(')d'.|'|urGrrrruGGrrurGrrrrrr无界基本解---无界空间Green函数的叠加（seebelowforthesolution）.例2，无界空间HelmholtzEquation：22()4(),()(,')4('),|0.|0.()(,')(')d'(seebelowforthesolution.(,'):Field;('):Source).urGrrrruGurGrrrrGrrr例3，无界空间波动方程:22222222(,),|0,|0,|0.(,;',')(')('),|0,|0,|0.ttttttaurtuuutaGrtrtrrttGGGt在含时Green函数(,;',')Grtrt中，为方便计，我们将它简记为(,').Grr2,HelmholtzEquationandLaplaceEquation解的积分形式（在定解问题中求G）设定解问题:22()4(),(1)()|(),(2)()(,')4('),(3)()|0.((4)nnuruufGrrrrGG对应于齐次边界条件)假设(,')Grr已经求出，方法见14.3.作算符运算：[(,')Eq.(1)()Eq.(3)]dGrrurr得22[(,')()()(,')]d4[(,')()()(,')]d.GrrururGrrrGrrrurrrr对上式左边利用Gauss公式ˆdd,FrFn其中,,FGuuG可得如下的Green恒等式：22[(,')()()(,')]d()d.nnGrrururGrrrGuuGSee童裕孙等，高等数学下册，2nd版，PP154-156&PP161-162.故有形式解1()(,')(')(,')()d(,')()d.4urGrrurGrrrrGrrurnn（5）3在(,')Grr已知的前提下，解（5）也不是()ur用(,')Grr表示的最后形式。这是因为(),()nurur在边界上还未知(最多知道他们在的线性组合—BCs).幸好(,')Grr的边界条件还可以以多种形式提出；只要选定G的边界条件为齐次，则(),()nurur或者其线性组合就可为已知的边界条件，从而最后确定().ur例如：1）第一类边界条件：0,1,|(),|0ufG1()(')(,')()d()d.4GrurGrrrrurn2）第二类边界条件：1,0,|(),|0.nnufG1()(')(,')()d(,')d.4ururGrrrrGrrn3)第三类边界条件：(,)0,takingEq.(2)Eq.(4),()|().nnGuGuuGGf得到1(')(,')()d(,')()d.4urGrrrrGrrf这些形式解的前提是G要已知或者可求出。实际问题中，G可能就不存在。例如在第二类边界条件中，0的问题[PossionEquation,点源存在，但边界“流”为零--|0.nG物理上不通(既产生又绝缘--矛盾），数学上无解]。再例如，构成本征值问题的2()0,()|0nuuu，不存在相对应的方程（3）和（4）.这是因为即使在边界条件完全相同的情况下，0（无源）而方程（3）的源为(').rr在此情况下G无解（或适当修改Green函数的意义），其实分离变量法已解决了问题，GF就多此一举了。GF要解决更为复杂的物理问题！3.Green函数的物理意义以为边界的区域，既无论原方程是否齐次（即内有、无源），又无论原边界条件是否齐次（即上有、无源），GF总是定义在内除一点(点源)以外方程的非齐次项处处为零（问题本身总要有非齐次项，如源于边界条件）但要化成齐次边界条件的定解问题的解。因此GF是“点源影响函数”或者“作用的传播函数”。对于所讨论的线性方程而言，一旦知道了相应问题的Green函数，只要再做两个积分，把原方程的非齐次项所反应的连续源分布对各点所产生的影响线性迭加起来，便给出原问题的解。这是线性迭加原理的最成功应用。齐次方程的本征值问题的本征值解可用于表示相应非齐次方程定解问题的GF(求法见最后一节)。14.2Green函数的性质41.Green函数由线性算子L和边界条件和初始条件决定：(,;',')(')(')LGrtrtrrtt，加上齐次边界条件和初始条件。2.Green函数的叠加性1).01:GGG非齐次方程特解（）+齐次方程通解（=0）.例如：0101('),0,(,')('),|().|0;|().LGrrLGLGrrrrGfGGf2).()0.LGG3.Green函数的对称性：若算子L是厄米的，则由L产生的G有*(,')(',).GrrGrr特别地，对于实变Green函数，(,')(',).GrrGrr1）Helmholtz方程:2'()(,')('),...........(1)(')|0...(':source')......(2)nGrrrrGG2''()(,'')(''),.........(3)('')|0....(:source'')...(4)nGrrrrGG作[(,'')Eq.(1)(,')Eq.(3)]dGrrGrrr，则方程右端变为('',')(','')GrrGrr，而左端22[(,'')(,')(,')(,'')]d.GrrGrrGrrGrrr利用Guass公式和Green恒等式，上式变为'''(''')d0nnGGGG[这是因为Eqs.（2）和（4）的行列式为零]。故：(','')('',')GrrGrr.2)GreenEquations:'''(,')('),...................(1)(')|0............................(2)[(,'')][('')](''),.....(3)('')|0......................................nnLGrrrrGGLGrrrrrrGG.....(4)作[(,'')Eq.(1)(,')Eq.(3)]d,GrrGrrr则方程右端变为('',')(','')GrrGrr，而左端{(,'')(,')(,')[(,')]}d0,GrrLGrrGrrLGrrr(因为L的厄米性)。故(','')('','),or.GrrGrrGG5推论，*();(,')(',)GGGrrGrr，如果G为实变函数。总之，1(')(,')()(,')(')d'(,')(')d.4urGrrurGrrrrGrrurnn物理意义：点源'r处产生的势传播到r处：(,')Grr,再对整体求和，加上边界条件产生的势，总和为r处的势分布。4.Green函数的奇异性及其随空间维数降低的减弱性。1).（有）无界区域Green函数----基本解0(,')Grr.a).无界区域：(').LGrr设解为0(,')Grr（详解见下），则基本解0(,')Grr为有限形式，且中心对称（与|'|rr有关）以及关于'rr发散（奇异性与维度D有关）。b).有界区域：01GGG,0(')LGrr和11000,||(,')LGGGGr.因为基本解0G有奇异性，所以直接影响分布。点源'r通过边界，i.e.，0(,')Gr对r处的分布有直接影响。2).HelmholtzEquation22()(')kGrr的解(,')Grr在'rr附近的奇异性（具体求解见下节）：3D:11,|'|Grrr发散;2D:1ln|'|G，对数发散;1D:G连续，00||1xxxxGG，导数发散。由此可见：Green函数的奇异性随空间维数的降低而减弱。14.3Green函数的求法1.方程齐次化法：将非齐次项(')rr变为边界条件。1)HelmholtzEquation:a).1D:2''(')GkGxx22'''0,(')('')|(').xxGkGxxGkGxx|'|'2|'|'''[(,')(',)]/(2):1('')d'|'|()/(2).ikxxxikxxxxxGAeG