数学物理方法课件Chp16.pptVIP

保角变换法求解定解问题 第十四章 王建东 沙河校区计算机楼东206 jdwang@uestc.edu.cn 在许多物理问题中（如电学、热学、光学、流体力学和弹性力学等）经常会遇到解平面场的拉普拉斯方程或泊松方程的问题。尽管可用前几章的理论方法如分离变量法或格林函数法等来解决，但当边值问题中的边界形状变得十分复杂时，分离变量法和格林函数法却显得十分困难，甚至不能解决。对于复杂的边界形状，拉普拉斯方程定解问题常采用保角变换法求解。 保角变换法解定解问题的基本思想是：通过解析函数的变换（或映射，这部分知识在复变函数论中已经学习过）将z平面上具有复杂边界形状的边值问题变换为平面上具有简单形状（通常是圆、上半平面或带形域）的边值问题，而后一问题的解易于求得。于是再通过逆变换就求得了原始定解问题的解。 14.1 保角变换与拉普拉斯方程边值问题 的关系 在复变函数论中我们已经知道，由解析函数ω=f(z)实现的从z平面到ω平面的变换在f?(z)≠0的点具有保角性质，因此这种变换称为保角变换。下面我们主要讨论一一对应的保角变换，即假定ω=f(z)和它的反函数都是单值函数，或者如果它们之中有多值函数就规定取它的黎曼面的一叶。 证明：利用复合函数求导法则 定理14.1.1 如果将由z=x+iy到ω=u+iv的保角变换看成为二元（实变）函数Φ(x,y)的变量由x、y到u、v的变量代换，则z平面上的边界变成了ω平面上的边界。我们能证明，如果Φ(x,y)满足拉普拉斯方程，则经过保角变换后得到的Φ(u,v)也满足拉普拉斯方程 同理： 两式相加 解析函数ω=f(z)=u+iv的C—R条件： 解析函数的实部和虚部分别满足拉普拉斯方程： 解析函数的导数： 因此，可得 对于保角变换，ω?(z)=f?(z)≠0， 因而只要Φ(x,y)满足拉普拉斯方程，则Φ(u,v)也满足拉普拉斯方程，即 这样我们就有结论：如果在z=x+iy平面上给定了Φ(x,y)的拉普拉斯方程边值问题，则利用保角变换ω=f(z)，可以将它转化为ω=u+iv平面上Φ(u,v)的拉普拉斯方程边值问题。 同理可得，在单叶解析函数ω=f(z)变换下，泊松方程： 仍然变为泊松方程： 即在保角变换 下，泊松方程中的电荷密度发生了变化。 经保角变换后仍然为亥姆霍兹方程： 对于亥姆霍兹方程： 14.2 保角变换法求解定解问题 例14.2.1 设有半无限平板y0，在边界y=0上，|x|a (a0)范围内保持温度u=u0， |x|a范围内保持温度u=0。求平板上的稳定温度分布。 解：根据题意可得描述上述问题的定解问题 分两步作保角变换： (1) 作分式线性变换 z平面上实轴(z=x, y=0)的对应关系: (i) 若|x|a (-axa) ，则 (ii) 若|x|a (x-a或xa) 对于x-a (a0)，则x+a0、 x-a-2a0， 因此 ② 对于xa (a0)，则x+a2a0、 x-a0，因此 如图所示，原定解问题中的边界条件中对应于|x|a范围温度为u0，变换后对应到ζ1平面的负实轴(ξ10)温度为u0；而|x|a温度为0则对应于变换后的ζ1平面的正实轴温度保持为0。 (2) 作变换 把ζ1平面的上半平面变成ζ平面上平行于实轴，宽为π的一个带型区域， ζ1平面的正实轴变换为ζ平面的实轴(正实轴辐角为0，故对应于η=0，温度u|y=0=0)， ζ1平面的负实轴变换为ζ平面的平行于实轴的直线(负实轴辐角为π，故对应于η=π，温度u|y=0=u0)。 原定解问题变为： 上述定解问题表明：等温线是与实轴平行的直线η=常数，热流线则是与虚轴平行的直线ξ=常数，即沿虚轴没有温度梯度，因此，在(ξ, η)坐标系中，拉普拉斯方程的解与ξ无关。这样，定解问题又简化为： 易解得： 回到变换前的z平面，则： 例14.2.2 试求平面静电场的电势分布φ(x,y) ： 解：作保角变换ω=lnz，使上半z平面变成ω平面上的带形区域： 原定解问题变为： 上述定解问题表明：等势线平行于虚轴，即电势φ与u无关。这样，定解问题又简化为： 易解得： 回到变换前的z平面： 可得原定解问题的解为： 例14.2.3 把柱面充电到： 试用保角变换法求解一半径为a的无限长导体圆柱壳内的电场分布。 解：根据题意可得描述上述问题的定解问题： 分三步作保角变换： (1) 作变换 将任意的圆周变换为单位圆，即把原图缩小为原来的1/a； (2) 作变换 把|z1|1变换为Imz20，即将原边界的下半圆周变换为负半实轴，上半圆周变为正半实轴； (3) 作变换 把z2平面的