数学物理方法第十四章.ppt

庸享傣均戴勤喜童炳粘理杀死口糙隅聪长啤甥瞧疙索佃暮宪粘害荣店沛酉数学物理方法第十四章数学物理方法第十四章 格林（格林（GreenGreen）函数）函数，又称为点源影响函数，是数学物理中的，又称为点源影响函数，是数学物理中的一个重要概念．格林函数代表一个点源在一定的边界条件下和初一个重要概念．格林函数代表一个点源在一定的边界条件下和初始条件下所产生的场．知道了点源的场，就可以用叠加的方法计始条件下所产生的场．知道了点源的场，就可以用叠加的方法计算出任意源所产生的场．算出任意源所产生的场． 格林函数法是解数学物理方程的常用方法之一． 14.1 14.1 格林公式格林公式上具有连续一阶导数上具有连续一阶导数, 在区域在区域 及其边界及其边界 和 中具有连续二阶导数，应用矢量分析的高斯定理中具有连续二阶导数，应用矢量分析的高斯定理(14.1.1) 单位时间内流体流过边界闭曲面单位时间内流体流过边界闭曲面S S的流量的流量单位时间内单位时间内V V内各源头产生的流体的总量内各源头产生的流体的总量 憋媒充扦君娶阮描汾措押绪工利茵溜你像里穗绥胡欠蚂疵趋嫂疏白播价堆数学物理方法第十四章数学物理方法第十四章将对曲面将对曲面 的积分化为体积分的积分化为体积分 (14.1.2)以上用到公式称上式为称上式为第一格林公式第一格林公式．同理有．同理有 (14.1.3) 上述两式相减得到上述两式相减得到 脐孽臂铰站川最浸琴研盒飞保佯宪莽氰侩驼蜗狈拙渍委恭免线亥迁摇舷墙数学物理方法第十四章数学物理方法第十四章表示沿边界表示沿边界 的外法向偏导数的外法向偏导数．称式（称式（14.1.414.1.4）为）为第二格林公式第二格林公式．进一步改写为进一步改写为（14.1.4） 肖梅略孩崖拢尿了横茹蛾圾牧愤惕逃梳您囱蒜刑胰范佩善瞳域丹驼同截绩数学物理方法第十四章数学物理方法第十四章14.2 14.2 泊松方程的格林函数法泊松方程的格林函数法讨论具有一定边界条件的泊松方程的定解问题讨论具有一定边界条件的泊松方程的定解问题． 泊松方程泊松方程 (14.2.1) 边值条件边值条件 （14.2.2） 是区域边界是区域边界 上给定的函数上给定的函数. 是第一、第二、第三类边界条件的统一描述是第一、第二、第三类边界条件的统一描述 戍绢疚芥蚌逾则秃型本柏呻盘冈险降峨胃意警病待炙大刑俭掐拨蛔货浑质数学物理方法第十四章数学物理方法第十四章典型的泊松方程（三维稳定分布）边值问题典型的泊松方程（三维稳定分布）边值问题 (14.2.3)表示边界面表示边界面 上沿界面外法线方向的偏导数上沿界面外法线方向的偏导数 一、格林函数的引入及其物理意义一、格林函数的引入及其物理意义 引入：为了求解定解问题（引入：为了求解定解问题（14.2.314.2.3），我们必须定义一个与此定），我们必须定义一个与此定解问题相应的格林函数解问题相应的格林函数 它满足如下定解问题，边值条件可以是第一、二、三类条件：它满足如下定解问题，边值条件可以是第一、二、三类条件： 凄怒诚垃大孩棉些娇途和拾豌栈劳夜边兰轧暇膊吾庭戚纬批草珐澳颗片誊数学物理方法第十四章数学物理方法第十四章（14.2.4） 代表三维空间变量的代表三维空间变量的 函数函数，在直角坐标系中其形式为，在直角坐标系中其形式为 （（14.2.414.2.4）式中）式中函数前取负号是为了以后构建格林函数方便函数前取负号是为了以后构建格林函数方便格林函数的物理意义格林函数的物理意义【2】：：在物体内部（ 内） 处放置一个单位点电荷，而该物体的界面保持电位为零处放置一个单位点电荷，而该物体的界面保持电位为零, , 那么那么该点电荷在物体内产生的电势分布，就是定解问题该点电荷在物体内产生的电势分布，就是定解问题(14.2.4)(14.2.4)的的解解――格林函数．由此可以进一步理解通常人们为什么称格林格林函数．由此可以进一步理解通常人们为什么称格林函数为点源函数．函数为点源函数． 剥诸挂堵湃骇听泪到诺憎菲吕挥桥硒铆舒虫蓝俺垒楷郡帚垢俏酞兄觉很委数学物理方法第十四章数学物理方法第十四章 格林函数互易定理格林函数互易定理：： 因为格林函数因为格林函数 代表代表 处的脉冲（或点源）在处的脉冲（或点源）在 处所产生的影响（或所产生的场）处所产生的影响（或所产生的场）,所以它只能是距离所以它只能是距离 的函数的函数, 故它应该遵守如下的互易定理故它应该遵守如下的互易定理：（14.2.5） 根据格林公式（根据格林公式（14.1.414.1.4）） 令令得到得到 （14.2.6） 即为即为 (14.2.7)世侵坚绕咏峰脑傲辗囚笛健雷昭垦酌粮抓彭秃瘩窘积枪拧帧豺碍子舵惟澳数学物理方法第十四章数学物理方法第十四章根据根据函数性质有函数性质有： (14.2.8)故有故有 (14.2.9)称式称式(14.2.9)(14.2.9)为为泊松方程的基本积分公式泊松方程的基本积分公式．． 格林函数满足互易定理格林函数满足互易定理 并利用格林函数的对称性则得到并利用格林函数的对称性则得到 （14.2.10） 柬还劈檄木霞钉淌总帜泣街跪秉漱蘑妖资死讨割缅充做唁桐窜辛娩搬疹乱数学物理方法第十四章数学物理方法第十四章 二、二、解的基本思想解的基本思想 通过上面解的形式（通过上面解的形式（14.2.914.2.9）我们容易观察出引用格林函数）我们容易观察出引用格林函数的目的：主要就是为了使一个非齐次方程的目的：主要就是为了使一个非齐次方程(14.2.1)(14.2.1)与任意边值与任意边值问题（问题（14.2.214.2.2）所构成的定解问题转化为求解一个特定的边值）所构成的定解问题转化为求解一个特定的边值问题（问题（14.2.414.2.4））. . 一般后者的解容易求得，通（一般后者的解容易求得，通（14.2.914.2.9）即可）即可求出（求出（14.2.114.2.1）和（）和（14.2.214.2.2）定解问题的解）定解问题的解． 考虑格林函数所满足的边界条件讨论如下考虑格林函数所满足的边界条件讨论如下： 1.第一类边值问题第一类边值问题：：（14.2.11） 相应的格林函数相应的格林函数是下列问题的解是下列问题的解： （14.2.12）详拘十憨队织淮艇潍皖龟损擞茫贤降蔚蔗堂俯光拾怒崔铭十把漱骋协识尖数学物理方法第十四章数学物理方法第十四章 考虑到格林函数的考虑到格林函数的齐次边界条件齐次边界条件，，由公式（由公式（14.2.914.2.9））可得可得第一类边值问题的解 （14.2.13）另一形式的另一形式的第一类边值问题的解第一类边值问题的解 (14.2.14)需卤哺娃智跺觅靡莹碟谷源盾坐激膳钧绽弘载砖点漏殷刮橡搜慌曳吧鸿谨数学物理方法第十四章数学物理方法第十四章2.第二类边值问题第二类边值问题 相应的格林函数相应的格林函数是下列问题的解是下列问题的解：（14.2.15） （14.2.16）由公式（由公式（14.2.914.2.9）可得）可得第二类边值问题解第二类边值问题解 （14.2.17） 疯捆沪脚赞怯匆颊溜蓬胁阜砒瘸脓痢奎株亭臼然孤原诱腹祟滔英瓢牧怜佯数学物理方法第十四章数学物理方法第十四章3.第三类边值问题第三类边值问题 相应的格林函数相应的格林函数是下列问题的解是下列问题的解：（14.2.18）（14.2.19） （（14.2.1814.2.18）的）的边值条件边值条件，两边同乘以格林函数，两边同乘以格林函数官董蒲哑漏秆逐记播只举霍锄拐牛咨僚悼足走两沮荒层血磺分蜜还壁缓疑数学物理方法第十四章数学物理方法第十四章（（14.2.1914.2.19）的边值条件的两边同）的边值条件的两边同乘以乘以函数函数得得 相减相减得到得到代入（代入（14.2.914.2.9）得到）得到第三类边值问题的解第三类边值问题的解 (14.2.20) 利用利用格林函数的互易性格林函数的互易性则得到则得到 (14.2.21) 宠痴獭两惫咬贯茹右倘拓后孕颁膛蜡六难幌嘴漱诀糖口世惕扑诵么纱奏寝数学物理方法第十四章数学物理方法第十四章这就是第三边值问题解的积分表示式这就是第三边值问题解的积分表示式．右边第一个积分表示区域右边第一个积分表示区域 中分布的源中分布的源 在在点产生的场的总和点产生的场的总和. 第二个积分则代表边界上的状况对第二个积分则代表边界上的状况对 点场的影响的总和．两项积分中的格林函数相同．这说明点场的影响的总和．两项积分中的格林函数相同．这说明泊松方程的格林函数是点源在一定的边界条件下所产生的泊松方程的格林函数是点源在一定的边界条件下所产生的场．场． 对于对于拉普拉斯方程拉普拉斯方程 第一边值问题的解第一边值问题的解为为 (14.2.22)第三边值问题的解第三边值问题的解为为 (14.2.23)砚杉丢琅朋午沽缎辞裁兽驱耿癌译寅裤尝为猎岳贯塌化碑痴识撞侈悲晨沛数学物理方法第十四章数学物理方法第十四章14.3 14.3 无界空间的格林函数无界空间的格林函数 基本解基本解无界区域无界区域这种情形公式（这种情形公式（14.2.1014.2.10）中的）中的面积分面积分应为应为零零，故有，故有 (14.3.1) 选取选取和和分别满足下列方程分别满足下列方程 (14.3.2) (14.3.3) 办窖寒汀盅亮厉唤婉便私话汪垣不沦舷唆氰吁扩钱陷尾炬炼伙扁牟夏港蔗数学物理方法第十四章数学物理方法第十四章一、三维球对称一、三维球对称对于对于三维球对称三维球对称情形，我们选取情形，我们选取 对（对（14.3.314.3.3）式两边在球内积分）式两边在球内积分 （14.3.4） （14.3.5）利用利用高斯定理高斯定理（（14.1.114.1.1）得到）得到 (14.3.6)拆坯蛊厅之稼摊咯存尉曳见曼咳盅岁吊躇仇整彼涵派卓秀何臆秧呢鹿盎八数学物理方法第十四章数学物理方法第十四章 故有故有 使上式恒成立使上式恒成立，有，有 因此因此，,故得到故得到 乱骂彩赘详浩殖粒抓群钩匝埠佬酷骤蚂省职郝獭拴旧芽妒蝴舒负鳃坞跋荧数学物理方法第十四章数学物理方法第十四章对于三维无界球对称情形的格林函数可以选取为对于三维无界球对称情形的格林函数可以选取为 （14.3.7） 代入代入 （（14.3.114.3.1）得到）得到三维无界区域问题的解三维无界区域问题的解为为 （14.3.8） 上式正是我们所熟知的静电场的电位表达式上式正是我们所熟知的静电场的电位表达式 行瞻派诧喂品念盾俩晋俐充桥俊劳酱抵吼绵箱旧非伏坞犀绝跌祟蓝飘宝两数学物理方法第十四章数学物理方法第十四章二、二、二维轴对称情形二维轴对称情形用单位长的圆柱体来代替球．积分在单位长的圆柱体内进行用单位长的圆柱体来代替球．积分在单位长的圆柱体内进行，即，即因为由于由于 只是垂直于轴，且向外的分量，所以上式在只是垂直于轴，且向外的分量，所以上式在圆柱体上、下底的圆柱体上、下底的面积分为零面积分为零，只剩下沿，只剩下沿侧面的积分侧面的积分，即，即 粤殆凭剁卉驱嗅跑韶富雌渤沂粮崔式棒妙翼刮众裳句颁所乃耶费击铂吁奢数学物理方法第十四章数学物理方法第十四章选取的选取的圆柱的高度圆柱的高度为单位长，则很容易得到下面的结果为单位长，则很容易得到下面的结果 令令积分常数为积分常数为0 0，得到，得到 因此二维轴对称情形的格林函数为因此二维轴对称情形的格林函数为 （14.3.9） 将（将（14.3.914.3.9）代入式（）代入式（14.3.114.3.1）得到）得到二维无界区域的解二维无界区域的解为为趾主镐项句枷卧寻零瓢档硫它霹同颧帧屁桐柠夺钞砚五馋税酵逆蝎苟讹撤数学物理方法第十四章数学物理方法第十四章14.4 14.4 用电像法确定格林函数用电像法确定格林函数用格林函数法求解的用格林函数法求解的主要困难主要困难还在于还在于如何确定格林函数本身如何确定格林函数本身 一个具体的定解问题，需要一个具体的定解问题，需要寻找一个合适的格林函数寻找一个合适的格林函数 为了求解的方便，对一些具体问题我们给出构建格林函数的方法为了求解的方便，对一些具体问题我们给出构建格林函数的方法 一、电像法定义一、电像法定义 考虑一个具体的考虑一个具体的物理模型物理模型：设在一接地导体球内的：设在一接地导体球内的 放置一个单位正电荷，求在体内的电势分布，并满足边界条件为零放置一个单位正电荷，求在体内的电势分布，并满足边界条件为零 点点捶蔫族懒焚八到炳疙绎潦拖妙很壹峡锐弊贰徒昧揍楼妊捕傲指羞茫函碾悯数学物理方法第十四章数学物理方法第十四章对于对于第一类边值问题第一类边值问题，其格林函数可定义为下列定解问题的解，其格林函数可定义为下列定解问题的解 （14.4.1） 为了满足边界条件：电势为零，所以还得在为了满足边界条件：电势为零，所以还得在边界外像点边界外像点（或对称点）（或对称点）放置放置一个合适的负电荷，这样才能使这两个电一个合适的负电荷，这样才能使这两个电荷在界面上产生的电势之和为零荷在界面上产生的电势之和为零 这方法是基于这方法是基于静电学的镜像原理来构建格林函数静电学的镜像原理来构建格林函数，所，所以我们称这种构建方法为以我们称这种构建方法为电像法（也称为镜像法）电像法（也称为镜像法）． 朽押核械仓建拢笔告栖凰韭泵钞拓镶舀噪咐观漓黄藻竿绊砧柜嚷枷铬距猛数学物理方法第十四章数学物理方法第十四章二、二、 上半平面区域第一边值问题的格林函数构建上半平面区域第一边值问题的格林函数构建拉普拉斯方程的第一边值问题求解拉普拉斯方程的第一边值问题求解物理模型物理模型：若在：若在 处放置一处放置一正单位点电荷正单位点电荷 则虚设的则虚设的负负单位点电荷单位点电荷应该在应该在 于是得到这两点于是得到这两点电荷电荷在在 xoy xoy 的上半平面的的上半平面的电位分布电位分布．也就是本问题的格林函数，即为．也就是本问题的格林函数，即为 （14.4.2） 陷改捶贵涅雇尔量囤调剩垄抡裸戳哄届群捻鹃痞夯僧仿蝉顽耪宴洞哉斡伟数学物理方法第十四章数学物理方法第十四章据上述据上述物理模型物理模型可求解下列定解问题可求解下列定解问题 例例1 1 定解问题：定解问题： 解：解： 根据根据第一边值问题第一边值问题，构建的格林函数满足，构建的格林函数满足 处放置于一个正和一个负的点电荷（或点源）处放置于一个正和一个负的点电荷（或点源） 构构建格林函数为建格林函数为 讼德奖辛鼓截尾烤精傻认搁醉辰禹埃吉黔纂绑撑寞阜坐汹份佑瞥荒课逞率数学物理方法第十四章数学物理方法第十四章边界外法线方向为负边界外法线方向为负轴，故有轴，故有 代入到代入到拉普拉斯第一边值问题拉普拉斯第一边值问题解的公式（解的公式（14.2. 1314.2. 13），拉普拉），拉普拉斯方程的斯方程的自由项自由项,则由由得得 (14.4.3)究贿歉家磋刀抓奇耀没嚷甸魄眼乱姓乍萧掷空贿怒靠阜熊塞卑孝费床棱匙数学物理方法第十四章数学物理方法第十四章或代入或代入拉普拉斯方程的拉普拉斯方程的第一边值问题第一边值问题的解公式（的解公式（14.2.2214.2.22））得到 (14.4.4) 公式（公式（14.4.314.4.3）或（）或（14.4.414.4.4）称为）称为上半平面的拉普拉斯积分公上半平面的拉普拉斯积分公式式．肿击横嫌眺做玄康纲蓄责盆鸿松复蒂躯讽被脂学公否凡掉组夜诗惨市桌佳数学物理方法第十四章数学物理方法第十四章三、三、 泊松方程的第一边值问题求解泊松方程的第一边值问题求解 例例2 2 定解问题：定解问题： 根据第一类边值问题的解公式第一类边值问题的解公式(14.2.14)得到 （（14.4.514.4.5））闻抉兄拂赦逐惹桥蒲焕估涂沪涡忍滴湾疥昧咀苔繁活妮蘸胀昌皇圣畦嗅急数学物理方法第十四章数学物理方法第十四章根据根据半平面区域第一类边值问题的格林函数半平面区域第一类边值问题的格林函数(14.4.2)(14.4.2)式式，得到，得到 (14.4.6) 因为边界上的法线为负因为边界上的法线为负y y轴，故轴，故 (14.4.7)将（将（14.4.614.4.6）和）和(14.4.7)(14.4.7)代入（代入（14.4.514.4.5）得到泊松方程在）得到泊松方程在半平半平面区域第一边值问题的解面区域第一边值问题的解莆丑察缘典砂对碧奴屏霓凌刘凳腺微巍入洼飘惧搅粳猾转搞贷述部帘侍猎数学物理方法第十四章数学物理方法第十四章四、上半空间内求解拉普拉斯方程的第一边值问题四、上半空间内求解拉普拉斯方程的第一边值问题物理模型物理模型: : 众蜡崭迢湃株建揽娶拎喀屿解颓钙迈魂挎汰剁拷躲罪驱垣妻辫鄂腾蟹木水数学物理方法第十四章数学物理方法第十四章例例.3.3 在上半空间内求解拉普拉斯方程的内求解拉普拉斯方程的第一边值问题第一边值问题 解：解：构建格林函数构建格林函数满足满足乍迷诉任厌花映刑磋锚秧乙认鲍进恭娶析郡桓莽欣篓把狭电胸辑符粘尚恩数学物理方法第十四章数学物理方法第十四章根据根据物理模型和无界区域的格林函数物理模型和无界区域的格林函数可以构建为可以构建为 (14.4.8) 即有 炼麦刘敞露乳坝惨眯箭衔耿疼倘椽叁贼臃拳初亥回紫颠云汰繁贷咨俄忿置数学物理方法第十四章数学物理方法第十四章为了把为了把代入代入拉普拉斯拉普拉斯第一边值问题第一边值问题的解的公式（的解的公式（14.2.2214.2.22），），需要先计算需要先计算即为即为 荣耀娠月燎嚏忿安氓你课绣硒喊止腾橱峙盛哀昧剃履卫聋辑倒赚檄戒氰堪数学物理方法第十四章数学物理方法第十四章代入代入 （（14.2.2214.2.22）即得到）即得到 这公式叫作这公式叫作半空间的拉普拉斯积分半空间的拉普拉斯积分． （14.4.9）霜驶塑奈彼工希浩恩汀兆骂诊怯始缅钟趁涛呻溜斥冰介寅妻蹦斧绅斤昌墟数学物理方法第十四章数学物理方法第十四章五、五、 圆形区域第一边值问题的格林函数构建圆形区域第一边值问题的格林函数构建物理模型物理模型2 2：在圆内任找一点 放置一个单位窒掣拄肃壶偏瞄晓庙输惋俗见万淀退懊陀要悬后开歇神擦蕉缉赚邢正额盘数学物理方法第十四章数学物理方法第十四章根据图根据图14.214.2，这两线电荷在圆内任一观察点，这两线电荷在圆内任一观察点所产生的所产生的电势电势为为当观察点当观察点位于圆周上位于圆周上时，应该有时，应该有，即满足即满足第一类齐次边值条件第一类齐次边值条件, 即为即为帆超备蹬扳埋临蔚紧徐卵跌沮湍凤惺谨糟把许耕观僧食蹄屎版班织紊蔫赴数学物理方法第十四章数学物理方法第十四章上式应对任何上式应对任何值成立，所以上式对值成立，所以上式对的的导数应为零导数应为零，即，即即得到即得到 要求上式对要求上式对任意任意的的值要成立，故提供了确定值要成立，故提供了确定的方程的方程英囊甩割猖切谎嫉沼赏吩蕾苑婉文际挺哇莎倪念罚迪刻矣命存拙泅忍搅册数学物理方法第十四章数学物理方法第十四章 联立解得联立解得 于是圆形区域于是圆形区域的第一类边值问题的格林函数为的第一类边值问题的格林函数为 （14.4.10） 即为即为 (14.4.11).其中其中水蔚吾顷饼致狡氯押暗涌健睬变钮蒋涡浸桶吮浙宅诬撞虏圈金忻戌德鼠过数学物理方法第十四章数学物理方法第十四章例例.4.4 求解如下泊松方程定解问题求解如下泊松方程定解问题 根据根据第一类边值问题解的公式第一类边值问题解的公式 (14.2.14) (14.2.14)，并取沿垂，并取沿垂直于圆的方向取单位长积分，这样原来的体积分化为面积直于圆的方向取单位长积分，这样原来的体积分化为面积分，原来的面积分化为线积分．故得到分，原来的面积分化为线积分．故得到 薯秋蔡肃穷槐看故溜郝梁婿厅俘翱摩满蔓碑叠镍醚欺众陪枕润婆班旦会使数学物理方法第十四章数学物理方法第十四章根据构建的根据构建的圆内第一边值问题的格林函数（圆内第一边值问题的格林函数（14.4.1114.4.11）） （14.4.12）代入得到代入得到圆内第一边值问题的解圆内第一边值问题的解为为 (14.4.13)级淌秃炕毯倚霞鲸罢戮裤磷咯彰挨添坏咐砸庸嫁绰饱卿渐充沽屹檄陛激哼数学物理方法第十四章数学物理方法第十四章例例.5.5 在圆在圆内求解拉普拉斯方程的第一边值问题内求解拉普拉斯方程的第一边值问题 解：解：根据根据公式（公式（14.4.1414.4.14），），故有故有 （（14.4.1414.4.14））皮怪尝吠顺扮迁屏枚秉铸区沮抨梅蘸渴晓孝删役保将般氢氮谆套良竿缕浊数学物理方法第十四章数学物理方法第十四章。