浅谈Green函数法

word格上传真是灾难
好在当时留了markdown的副本

调整了一下文本的结构

1 Green函数的相关概念

Green函数又称点源函数或影响函数，表示在$M_0$点的点源对任意场点$M$的影响，记为$G(M,M_0)$. 如果Green函数$G(M,M_0)$已知，则分布在任意区域$\Omega$内的“源”对$M$点的影响可以写成
$$u(M)={\int }_{\Omega }^{}G(M,M_0)\rho (M_0)\mathrm{d}{M}_0$$
其中$\rho$为源密度，$\mathrm{d}{M}_0$为包含$M_0$的微元.
当所要讨论的物理场带有边界条件和（或）初始条件时，Green函数还要与这些边界条件和（或）初始条件有关，因此一般说来，Green函数是一个点源在一定的边界条件和（或）初始条件下所产生的场，利用Green函数，可以求出任意分布的源所产生的场.

1.1 Poisson方程边值问题的Green函数法

三维Poisson方程
$$\{\begin{array}{l}\Delta u(M)=-h(M),M\in \Omega \\ \left[\alpha \frac{\partial u}{\partial n}+\beta u\right]{|}_{\partial \Omega }^{}=g(M),M\in \Omega \end{array}$$
其中$\partial \Omega$为$\Omega$的边界.
设$G(M,M_0)$是满足方程
$$\Delta G(M,M_0)=-\delta (M,M_0),M\in \Omega$$
的Green函数，则三维Poisson方程的解可以表示为
$$u(M)={\iiint }_{\Omega }^{}G(M,M_0)h(M_0)\mathrm{d}{\Omega }_0+{\iint }_{\partial \Omega }^{}G(M,M_0)\frac{\partial u}{\partial n_0}\mathrm{d}{S}_0-{\iint }_{\partial \Omega }^{}u(M_0)\frac{\partial G(M,M_0)}{\partial n_0}\mathrm{d}{S}_0.$$

1. 第一类边界条件$\alpha =0$，即$$u{|}_{\partial \Omega }^{}=\frac{1}{\beta }g(M)=f(M)$$
   这时要求$G(M,M_0)$满足第一类齐次边界条件$G(M,M_0)|\partial {\Omega }^{}=0$ ，方程的解为
   $$u(M)={\iiint }_{\Omega }^{}G(M,M_0)h(M_0)\mathrm{d}{\Omega }_0-{\iint }_{\partial \Omega }^{}u(M_0)\frac{\partial G(M,M_0)}{\partial n_0}\mathrm{d}{S}_0.$$
2. 第二类边界条件$\beta =0$，即$$\frac{\partial u}{\partial n}{|}_{\partial \Omega }^{}=\frac{1}{\alpha }g(M)=f(M)$$
   这时要求$G(M,M_0)$满足边界条件
   $$\frac{\partial G}{\partial n}{|}_{\partial \Omega }^{}=-\frac{1}{\Sigma }$$
   其中$\Sigma$为$\partial \Omega$的面积，方程的解为
   $$u(M)=c+{\iiint }_{\Omega }^{}G(M,M_0)h(M_0)\mathrm{d}{\Omega }_0+{\iint }_{\partial \Omega }^{}G(M,M_0)f(M_0)\mathrm{d}{S}_0$$
   其中$c$为待定常数.
3. 第三类边界条件，定解问题为三维Poisson方程的形式
   这是要求$G(M,M_0)$满足第三类其次边界条件
   $$\left[\alpha \frac{\partial G}{\partial n}+\beta G\right]{|}_{\partial \Omega }^{}=0$$
   方程的解为
   $$u(M)={\iiint }_{\Omega }^{}G(M,M_0)h(M_0)\mathrm{d}{\Omega }_0+\frac{1}{\alpha }{\iint }_{\partial \Omega }^{}G(M,M_0)g(M_0)\mathrm{d}{S}_0$$

1.2 Green函数的对易性质

前两天被问了Green函数如何证明对易性，这是一个很有意思的练习题。
大致思路是作差，然后用极值必定在边界上取得，在边界上都是0，所以恒为零即恒相等。
$$G(M,M_0)=G(M_0,M)$$

1.3 Green函数的一般求法

1. 无界区域的Green函数
   三维无界区域Green函数，即方程$\Delta G=-\delta (x-x_0,y-y_0,z-z_0)$的解是
   $$G(M,M_0)=\frac{1}{4\pi r_{M{M}_0}}=\frac{1}{4\pi \sqrt{(x-x_0{)}^2+(y-y_0{)}^2+(z-z_0{)}^2}}$$
   二维无界区域Green函数,即方程$\Delta G=-\delta (x-x_0,y-y_0)$的解是
   $$G(M,M_0)=\frac{1}{2\pi }\ln \frac{1}{r_{M{M}_0}}=\frac{1}{2\pi }\ln \left[\frac{1}{\sqrt{(x-x_0{)}^2+(y-y_0{)}^2}}\right]$$
2. 本征函数求下列第一类齐次边值问题
   $$\{\begin{array}{l}\Delta G(M,M_0)+\lambda G(M,M_0)=-\delta G(M,M_0),M\in \Omega \\ G(M,M_0){|}_{\partial \Omega }^{}=0\end{array}$$
   Green函数为
   $$G(M,M_0)=\sum _n\frac{1}{{\lambda }_n-\lambda }\overline{{\psi }_n}(M_0)\psi (M)$$
   其中$\lambda$和$\psi$为本征值问题
   $$\{\begin{array}{l}\Delta \psi (M)+\lambda \psi (M)=0,M\in \Omega \\ \psi (M){|}_{\partial \Omega }^{}=0\end{array}$$
   的本征值和本征函数.

1.4 电像法与某些特殊区域的第一类边界条件问题的Green函数

把带有边界问题的Green函数分成无界情况和边界的影响两部分，将边界的影响看成源点关于边界的对称点（成为像点）的影响，再通过物理的和几何的知识求出这一影响的方法称为电像法.

电像法的正确性是唯一性定理保证的，所以请放心大胆的使用。
电像法给了延拓的灵感。

1. 球域$(\rho \le a)$的Green函数
   Laplace方程第一类边值问题
   $$\{\begin{array}{l}\Delta u=0,\rho <a\\ u{|}_{\rho =a}^{}=f(M)\end{array}$$
   的Green函数为
   $$G(M,M_0)=\frac{1}{4\pi r_{M{M}_0}}-\frac{a/{\rho }_0}{4\pi r_{M{M}_1}}$$
   定解问题的解为
   $$u(\rho ,\theta ,\phi )=\frac{a}{4\pi }{\int }_0^{2\pi }{\int }_0^{\pi }f({\theta }_0,{\phi }_0)\frac{a^2-{\rho }^2}{(a^2+{\rho }^2-2a\rho \cos \gamma {)}^{3/2}}\sin {\theta }_0\mathrm{d}{\theta }_0\mathrm{d}{\phi }_0$$
   其中$a$为球面半径，其余量$r^{}$