偏微分方程(Partial Differential Equation III)

---

偏微分方程(Partial Differential Equation I)
偏微分方程(Partial Differential Equation II)
偏微分方程(Partial Differential Equation III)
偏微分方程(Partial Differential Equation IV)

---

参考文献：

《数学物理方程》| 季孝达
《数学物理方法》| 吴崇试
《数学物理方法》| 梁昆淼
MOOC北京大学《数学物理方法》| 吴崇试 、高春媛

积分变换法

在积分变换中我们曾用拉普拉斯变换方法求解常微分方程。经过变换，常微分方程变为代数方程，解出代数方程，再进行反演就得常微分方程的解。积分变换在数学物理方程中亦有广泛的应用。

傅里叶变换法

傅里叶变换法常用于求解无界空间（含一维半无界空间）定解问题。本节通过几个例子给出几个重要的解的公式。

Fourier 积分定理：若 $f(x)$ 在 $(-\infty ,+\infty )$上满足：
(1) 在任一有限区间上满足狄利克雷(Dirichlet)条件¹；
(2) 在无限区间 $(-\infty ,+\infty )$上绝对可积，即 ${\int }_{-\infty }^{+\infty }|f(x)|dx$ 收敛
那么，对任意 $x\in (-\infty ,+\infty )$
$$f(x)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }[{\int }_{-\infty }^{+\infty }f(\tau )e^{-ik\tau }\text{d}\tau ]e^{ikx}\text{d}k\text{\hspace{1em}(1.1)}$$
在间断点处，上式左端为 $\frac{1}{2}[f(x^{-})+f(x^{+})]$

Fourier 变换：如果函数 $f(x)$ 满足Fourier 积分定理，由式 (1.1) 知，令
$$F(k)={\int }_{-\infty }^{+\infty }f(x)e^{-ikx}\text{d}x\text{\hspace{1em}(1.2)}$$
则有
$$f(x)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }F(k)e^{ikx}\text{d}k\text{\hspace{1em}(1.3)}$$
从上面两式可以看出，$f(x)$ 和 $F(k)$ 通过确定的积分运算可以互相转换。 $F(k)$称为 $f(x)$ Fourier 变换(Fourier transform)，或象函数(image function)，记为$F(k)=\mathcal{F}[f(x)]$ ；$f(x)$称为 $F(k)$ Fourier 逆变换(inverse Fourier transform)，或象原函数(original image function)，记为$f(x)={\mathcal{F}}^{-1}[F(k)]$ ；通常称$f(x)$与$F(k)$构成一个Fourier 变换对(transform pair)，记作 $f(x)\leftrightarrow F(k)$ 。

示例：求函数 $f(x)=e^{-b{x}^2}$ 的傅里叶变换，其中常数 $b>0$
解：由定义和分部积分得
$$\begin{array}{rl}F(k)& ={\int }_{-\infty }^{+\infty }f(x)e^{-ikx}dx={\int }_{-\infty }^{+\infty }{e}^{-ikx-b{x}^2}dx\\ & =-\frac{1}{\mathrm{i}k}{e}^{-ikx-b{x}^2}{|}_{-\infty }^{+\infty }-\frac{1}{\mathrm{i}k}{\int }_{-\infty }^{+\infty }2bx{e}^{-ikx-b{x}^2}dx\\ & =\frac{2b\mathrm{i}}{k}\mathcal{F}[xf(x)]=-\frac{2b}{k}\frac{\mathrm{d}}{\mathrm{d}k}F(k)\end{array}$$
取 $k=0$ 可得
$$F(0)={\int }_{-\infty }^{+\infty }{e}^{-b{x}^2}dx=\sqrt{\frac{\pi }{b}}$$
转化为解常微分方程初值问题
$$\{\begin{array}{l}{F}^{\prime }(k)+\frac{k}{2b}F(k)=0\\ F(0)=\sqrt{\frac{\pi }{b}}\end{array}$$
其解为
$$F(k)=\sqrt{\frac{\pi }{b}}\exp (-\frac{k^2}{4b})$$
即
$$\mathcal{F}[e^{-b{x}^2}]=\sqrt{\frac{\pi }{b}}\exp (-\frac{k^2}{4b})(b>0)\text{\hspace{1em}(1.4)}$$
特别的，取 $b=\frac{1}{4c^{2}t}$ 时，有
$$\mathcal{F}[\exp (-\frac{x^2}{4c^{2}t})]=2c\sqrt{\pi t}{e}^{-c^2{k}^{2}t}\text{\hspace{1em}(1.5)}$$
逆变换为
$${\mathcal{F}}^{-1}[e^{-c^2{k}^{2}t}]=\frac{1}{2c\sqrt{\pi t}}\exp (-\frac{x^2}{4c^{2}t})\text{\hspace{1em}(1.6)}$$
由式 (1.4) 还可以得到
$${\int }_0^{+\infty }{e}^{-b{x}^2}\cos kx\mathrm{d}x=\frac{1}{2}\sqrt{\frac{\pi }{b}}\exp (-\frac{k^2}{4b})\text{\hspace{1em}(1.7)}$$
多重傅里叶变换：设 $\mathbf{x}=(x_1,x_2,\cdots ,x_n)\in {\mathbb{R}}^n,\mathbf{k}=(k_1,k_2,\cdots ,k_n),\mathrm{d}\mathbf{x}=d{x}_{1}d{x}_2\cdots d{x}_n$。若 $f(\mathbf{x})$ 在 ${\mathbb{R}}^n$ 上连续，分片光滑且连续可积，令
$$F(\mathbf{k})={\int }_{{\mathbb{R}}^n}f(\mathbf{x})e^{-\mathrm{i}\mathbf{k}\cdot \mathbf{x}}d\mathbf{x}\text{\hspace{1em}(1.8)}$$
则有
$$f(\mathbf{x})=\frac{1}{(2\pi {)}^n}{\int }_{{\mathbb{R}}^n}F(\mathbf{k})e^{\mathrm{i}\mathbf{k}\cdot \mathbf{x}}d\mathbf{k}\text{\hspace{1em}(1.9)}$$
其中 $F(\mathbf{k})$ 称为 $f(\mathbf{x})$ 的 多重傅里叶变换，记为 $F(\mathbf{k})=\mathcal{F}[f(\mathbf{x})]$；$f(\mathbf{x})$ 称为 $F(\mathbf{k})$ 的 多重傅里叶逆变换，记为 $f(\mathbf{x})={\mathcal{F}}^{-1}[F(\mathbf{k})]$

求无限长弦的初值问题
$$\{\begin{array}{l}{u}_{tt}-a^2{u}_{xx}=0\\ u{|}_{t=0}=\varphi (x),\frac{\partial u}{\partial t}{|}_{t=0}=\psi (x)\end{array}$$
其中 $\varphi (x),\psi (x)$ 分别表示初始位移和初始速度。
解：应用傅里叶变换，即方程及初始条件两边同乘以 $e^{-\mathrm{i}kx}$ ，并对空间变量 $x$ 进行积分（时间变量 $t$ 视作参数）。
记 $U(k,t)=\mathcal{F}[u(x,t)]={\int }_{-\infty }^{+\infty }u(x,t)e^{-\mathrm{i}kx}dx$ ，运用用含参变量的积分及傅里叶变换的微分性质

则定解问题变为关于 $t$ 的常微分方程及初值条件
$$\{\begin{array}{l}{U}^{\prime \prime }+k^2{a}^{2}U=0\\ U(0)=\Phi (k),U^{\prime }(0)=\Psi (k)\end{array}$$
其中 $\Phi (k)=\mathcal{F}[\varphi (x)],\Psi (k)=\mathcal{F}[\psi (x)]$ 分别是 $\varphi (x),\psi (x)$ 关于 $x$ 的傅里叶变换。其解为
$$U(k,t)=\Phi (k)\cos kat+\frac{1}{ka}\Psi (k)\sin kat$$

最后，对 $U(k,t)$ 做傅里叶逆变换，用延迟性质和积分性质，结果是
$$u(x,t)=\frac{1}{2}[\varphi (x+at)+\varphi (x-at)]+\frac{1}{2a}{\int }_{x-at}^{x+at}\psi (\xi )d\xi$$
这个公式正是达朗贝尔(d’Alembert)公式。

求无限长杆的有源热传导问题
$$\{\begin{array}{l}{u}_t-a^2{u}_{xx}=f(x,t)\\ u{|}_{t=0}=\varphi (x)\end{array}$$
解：做傅里叶变换，定解问题变为
$$\{\begin{array}{l}{U}^{\prime }+k^2{a}^{2}U=F(k,t)\\ U(0)=\Phi (k)\end{array}$$
其中 $U(k,t)=\mathcal{F}[u(x,t)],F(k,t)=\mathcal{F}[f(x,t)],\Phi (k)=\mathcal{F}[\varphi (x)]$ 分别是 $u(x,t),f(x,t),\varphi (x)$ 关于 $x$ 的傅里叶变换。这个常微分方程初值问题的解为
$$U(k,t)=\Phi (k)e^{-k^2{a}^{2}t}+{\int }_0^{t}F(k,\tau )e^{-k^2{a}^2(t-\tau )}d\tau$$

最后，对 $U(k,t)$ 做傅里叶逆变换，用卷积定理，结果是
$$u(x,t)={\int }_{-\infty }^{+\infty }\varphi (\xi )\frac{1}{2a\sqrt{\pi t}}\exp [-\frac{(x-\xi {)}^2}{4a^{2}t}]d\xi +{\int }_0^t{\int }_{-\infty }^{+\infty }\frac{f(\xi ,\tau )}{2a\sqrt{\pi (t-\tau )}}\exp [-\frac{(x-\xi {)}^2}{4a^2(t-\tau )}]d\xi d\tau$$

傅里叶正弦变换或余弦变换：如果 $f(x)$ 定义在半无界区间 $[0,+\infty )$ 上，满足狄利克雷(Dirichlet)条件¹且绝对可积，则有傅里叶正弦变换
$$\begin{gathered} F(k)={\int }_0^{+\infty }f(x)\sin kx\text{d}x \\ f(x)=\frac{1}{2\pi }{\int }_0^{+\infty }F(k)\sin kx\text{d}k \end{gathered}$$
或余弦变换
$$\begin{gathered} G(k)={\int }_0^{+\infty }g(x)\cos kx\text{d}x \\ g(x)=\frac{1}{2\pi }{\int }_0^{+\infty }G(k)\cos kx\text{d}k \end{gathered}$$
对于半无界空间，存在自然边界条件
$$\underset{x\to \infty }{\lim }f(x)=0,\underset{x\to \infty }{\lim }{f}^{\prime }(x)=0$$
可以采用正弦变换或余弦变换，对于正弦变换
导数的正弦变换为
$$\begin{array}{rl}& {\int }_0^{+\infty }{f}^{\prime }(x)\sin kx\text{d}x\\ =& f(x)\sin kx{|}_0^{+\infty }-k{\int }_0^{+\infty }f(x)\cos kx\text{d}x\\ =& -k{\int }_0^{+\infty }f(x)\cos kx\text{d}x\end{array}$$
二阶导正弦变换为
$$\begin{array}{rl}& {\int }_0^{+\infty }{f}^{\prime \prime }(x)\sin kx\text{d}x\\ =& -k{\int }_0^{+\infty }{f}^{\prime }(x)\cos kx\text{d}x\\ =& -k[f(x)\cos kx{|}_0^{+\infty }+k{\int }_0^{+\infty }f(x)\sin kx\text{d}x]\\ =& kf(0)-k^{2}F(k)\end{array}$$

由此可见，对于二阶偏微分方程的定解问题，只有在半无界空间的 $x=0$ 端给出第一类边界条件时，才可以采用正弦变换。

同样对于余弦变换，也有
$$\begin{gathered} {\int }_0^{+\infty }{g}^{\prime }(x)\cos kx\text{d}x=-g(0)+k{\int }_0^{+\infty }g(x)\sin kx\text{d}x \\ {\int }_0^{+\infty }{g}^{\prime \prime }(x)\cos kx\text{d}x=-g^{\prime }(0)-k^{2}G(k) \end{gathered}$$
由此可见，对于二阶偏微分方程的定解问题，只有在半无界空间的 $x=0$ 端给出第二类边界条件时，才可以采用余弦变换。

求半无界杆的热传导问题
$$\{\begin{array}{ll}{u}_t-a^2{u}_{xx}=0& (x>0)\\ u{|}_{t=0}=0& \\ u{|}_{x=0}=u_0& \end{array}$$
解：采用傅里叶正弦变换，定解问题变为
$$\{\begin{array}{l}{U}^{\prime }+k^2{a}^{2}U=k{a}^2{u}_0\\ U(0)=0\end{array}$$
其中 $U(k,t)$ 是 $u(x,t)$ 关于 $x$ 的傅里叶正弦变换。这个常微分方程初值问题的解为
$$U(k,t)=\frac{u_0}{k}(1-e^{-k^2{a}^{2}t})$$

最后，对 $U(k,t)$ 做傅里叶逆变换，结果是
$$u(x,t)=\frac{2u_0}{\pi }{\int }_0^{+\infty }\frac{1}{k}(1-e^{-k^2{a}^{2}t})\sin kxdk$$
通常记误差函数 (error function)
$$\mathrm{e}\mathrm{r}\mathrm{f}(x)=\frac{2}{\sqrt{\pi }}{\int }_0^x{e}^{-s^2}ds$$
和余误差函数 (error function complement)
$$\mathrm{e}\mathrm{r}\mathrm{f}\mathrm{c}(x)=1-\mathrm{e}\mathrm{r}\mathrm{f}(x)=\frac{2}{\sqrt{\pi }}{\int }_x^{\infty }{e}^{-s^2}ds$$
故 $u(x,t)$ 可进一步变换为²
$$u(x,t)=\frac{2u_0}{\pi }[\frac{\pi }{2}-\frac{\pi }{2}\mathrm{e}\mathrm{r}\mathrm{f}(\frac{x}{2a\sqrt{t}})]=u_0\mathrm{e}\mathrm{r}\mathrm{f}\mathrm{c}(\frac{x}{2a\sqrt{t}})$$

求三维无界空间中的波动问题
$$\{\begin{array}{l}{u}_{tt}-a^2\Delta u=0\\ u{|}_{t=0}=\varphi (x,y,z),\frac{\partial u}{\partial t}{|}_{t=0}=\psi (x,y,z)\end{array}$$
解：作三重傅里叶变换，记 $\mathbf{r}=(x,y,z),\mathbf{k}=(k_1,k_2,k_3)$，定解问题变为
$$\{\begin{array}{l}{U}^{\prime \prime }+{\mathbf{k}}^2{a}^{2}U=0\\ U(0)=\Phi (\mathbf{k}),U^{\prime }(0)=\Psi (\mathbf{k})\end{array}$$
其中 $U(\mathbf{k},t),\Phi (\mathbf{k}),\Psi (\mathbf{k})$ 分别是 $u(\mathbf{r},t),\varphi (\mathbf{r}),\psi (\mathbf{r})$ 关于 $\mathbf{r}$ 的三维傅里叶变换。这个常微分方程初值问题的解为
$$U(\mathbf{k},t)=\Phi (\mathbf{k})\cos kat+\Psi (\mathbf{k})\frac{\sin kat}{ka}$$

其中 $k=|\mathbf{k}|=\sqrt{k_1^2+k_2^2+k_3^2}$ ，再进行傅里叶逆变换
$$\begin{array}{rl}u(\mathbf{r},t)& =\frac{1}{(2\pi {)}^3}\underset{-\infty }{\overset{+\infty }{\iiint }}[\Phi (\mathbf{k})\cos kat+\Psi (\mathbf{k})\frac{\sin kat}{ka}]e^{\mathrm{i}\mathbf{k}\cdot \mathbf{r}}d\mathbf{k}\\ & =\frac{1}{4\pi a}\frac{\partial }{\partial t}\underset{S_{at}^{\mathbf{r}}}{\iint }\frac{\varphi (\mathbf{r})}{at}dS+\frac{1}{4\pi a}\frac{\partial }{\partial t}\underset{S_{at}^{\mathbf{r}}}{\iint }\frac{\psi (\mathbf{r})}{at}dS\end{array}$$
上式称为泊松公式。式中 $S_{at}^{\mathbf{r}}$ 表示以 $\mathbf{r}$ 为圆心，以 $at$ 为半径的球面，$dS$ 表示 $S_{at}^{\mathbf{r}}$ 的面积元。

拉普拉斯变换法

拉普拉斯变换法适合求解初值问题，不管方程和边界条件是否为齐次的。

Laplace变换：设函数$f(t)$ 在 $t⩾0$ 时有定义，且积分 ${\int }_0^{+\infty }f(t)e^{-st}dt$ 收敛，则此积分所确定的函数
$$F(s)={\int }_0^{+\infty }f(t)e^{-st}\text{d}t\text{\hspace{1em}(2.1)}$$
称为函数 $f(t)$ 的 Laplace 变换，记为 $F(s)=\mathcal{L}[f(t)]$，函数 $F(s)$ 也可称为 $f(t)$的象函数。
Laplace逆变换：令 $s=\beta +i\omega$ ，则有
$$f(t)=\frac{1}{2\pi i}{\int }_{\beta -i\omega }^{\beta +i\omega }F(s)e^{st}\text{d}s(t>0)$$
称为 Laplace 逆变换，记为 $f(t)={\mathcal{L}}^{-1}[F(s)]$ 。在Laplace 变换中，只要求$f(t)$在 $[0,+\infty )$ 内有定义即可。为了研究方便，以后总假定在$(-\infty ,0)$ 内，$f(t)\equiv 0$
还可用留数就算拉普拉斯逆变换：设在复平面内只有有限个孤立奇点 $s_1,s_2,\cdots ,s_n$ ，实数 $\beta$使这些奇点全在半平面 $\text{Re}(s)<\beta$ 内，且 $\underset{s\to \infty }{\lim }F(s)=0$ ，则有
$$f(t)=\sum _{k=1}^n\text{Res}[F(s)e^{st},s_k](t>0)$$

求半无界杆的热传导问题
$$\{\begin{array}{ll}{u}_t-a^2{u}_{xx}=0& (x>0)\\ u{|}_{t=0}=0& (x>0)\\ u{|}_{x=0}=u_0& \end{array}$$
解：对方程和边界条件关于 $t$ 进行拉普拉斯变换，采用微分性质，变换结果为
$$\{\begin{array}{l}sU-a^2{U}^{\prime \prime }=0\\ U(0)=\frac{1}{s}{u}_0\end{array}$$
其中 $U(s,x)$ 是 $u(x,t)$ 关于 $t$ 的傅里叶正弦变换。这个常微分方程通解为
$$U(s,x)=A\exp (-\frac{\sqrt{sx}}{a})+B\exp (\frac{\sqrt{sx}}{a})$$
考虑到自然边界条件 $\underset{x\to \infty }{\lim }U$应为有限值，带入初值条件可求得
$$U(s,x)=\frac{1}{s}{u}_0\exp (-\frac{\sqrt{sx}}{a})$$
最后，对 $U(s,x)$ 进行拉普拉斯逆变换[^L1]
$$u(x,t)=u_0\mathrm{e}\mathrm{r}\mathrm{f}\mathrm{c}(\frac{x}{2a\sqrt{t}})$$

求长 $l$ 均匀细杆的热传导问题
$$\{\begin{array}{ll}{u}_t-a^2{u}_{xx}=0& (0<x<l)\\ u{|}_{t=0}=0& (0<x<l)\\ u{|}_{x=0}=u_0,\frac{\partial u}{\partial x}{|}_{x=l}=0& \end{array}$$
解：对方程和边界条件关于 $t$ 进行拉普拉斯变换，采用微分性质，变换结果为
$$\{\begin{array}{l}sU-a^2{U}^{\prime \prime }=0\\ U(0)=\frac{1}{s}{u}_0,U^{\prime }(l)=0\end{array}$$
其中 $U(s,x)$ 是 $u(x,t)$ 关于 $t$ 的傅里叶正弦变换。这个二阶常微分方程的解为
$$U(s,x)=\frac{1}{s}{u}_0\frac{\cosh \frac{(l-x)\sqrt{s}}{a}}{\cosh \frac{l\sqrt{s}}{a}}$$
最后，利用留数定理对 $U(s,x)$ 进行拉普拉斯逆变换
$$u(x,t)=u_0-\frac{\pi u_0}{4}\sum _{n=0}^{\infty }\frac{1}{2n+1}\sin \left(\frac{2n+1}{2l}\pi x\right)\exp \left[-{\left(\frac{2n+1}{2l}\pi \right)}^2{a}^{2}t\right]$$

基本解和格林函数

格林函数，又称点源影响函数，代表一个点源在一定的边界条件和（或）初始条件下所产生的场。均匀分布的函数可看做点源的叠加，该想法来源于静电场叠加原理。实际上，这种做法只不过是利用了偏微分方程的积分叠加原理。

泊松方程的基本解

引例：先举一个静电场的例子，设无界空间中电荷密度为 $\rho (\mathbf{r})$ ，这样在坐标 ${\mathbf{r}}_0=(x_0,y_0,z_0)$ 的体积元 $d{V}_0$ 内的电荷量为 $\rho ({\mathbf{r}}_0)d{V}_0$ ，它在空间点 $\mathbf{r}=(x,y,z)$ 产生的电势为
$$\frac{1}{4\pi {\epsilon }_0}\frac{\rho ({\mathbf{r}}_0)}{|\mathbf{r}-{\mathbf{r}}_0|}$$
根据电势叠加原理，可叠加求得任意密度分布引起的总电势分布
$$\phi (\mathbf{r})=\frac{1}{4\pi {\epsilon }_0}\iiint \frac{\rho ({\mathbf{r}}_0)}{|\mathbf{r}-{\mathbf{r}}_0|}d{V}_0$$
泊松方程的基本解：对于无界空间的泊松方程
$$\Delta u=f(\mathbf{r})\text{\hspace{1em}(1.1)}$$
在物理上可看做电荷密度分布 $-{\epsilon }_{0}f(\mathbf{r})$ 在无界空间的电势方程。为了研究点源产生的场， $\delta$ 函数恰是一个表示点源密度的函数，由 $\delta$ 函数的性质知
$$\iiint f({\mathbf{r}}_0)\delta (\mathbf{r}-{\mathbf{r}}_0)\mathrm{d}{V}_0=f(\mathbf{r})$$
这说明一般源 $f(\mathbf{r})$ 可看成 ${\mathbf{r}}_0$ 点的点源 $f({\mathbf{r}}_0)\delta (\mathbf{r}-{\mathbf{r}}_0)$ 的积分叠加，由第一章积分叠加原理知道，只需求出方程
$$\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\text{\hspace{1em}(1.2)}$$
的解 $G({\mathbf{r},\mathbf{r}}_0)$ ，便可得无界空间泊松方程的解
$$u(\mathbf{r})=\iiint G({\mathbf{r},\mathbf{r}}_0)f({\mathbf{r}}_0)\mathrm{d}{V}_0\text{\hspace{1em}(1.3)}$$
其中 $G({\mathbf{r},\mathbf{r}}_0)$ 称为格林函数，又称点源影响函数。
由于是在无界空间，不妨先做平移变换，求方程（拉普拉斯算符平移不变性）
$$\Delta U(\mathbf{r})=\delta (\mathbf{r})\text{\hspace{1em}(1.4)}$$
的解 $U(\mathbf{r})$，称为基本解，代表置于原点的点源引起的场。则
$$G(\mathbf{r},{\mathbf{r}}_0)=U(\mathbf{r}-{\mathbf{r}}_0)\text{\hspace{1em}(1.5)}$$
进而有
$$u(\mathbf{r})=\iiint U(\mathbf{r}-{\mathbf{r}}_0)f({\mathbf{r}}_0)\mathrm{d}{V}_0=U(\mathbf{r})\ast f(\mathbf{r})\text{\hspace{1em}(1.6)}$$
基本解的求法：在物理上， 基本解 $U(\mathbf{r})$ 描述了位于原点电荷量为 $-{\epsilon }_0$ 的电荷在无界空间 $\mathbf{r}$ 处的电势。这里基本解没有定解条件限制，因此不是惟一的，通常根据问题的物理意义和数学的需要选定其中一个。
下面介绍数学上的一种求解方法。将方程用球坐标表示，当 $r\ne 0$ 时，由于函数关于 $r$ 对称，所以方程化为
$$\frac{1}{r^2}\frac{d}{dr}(r^2\frac{dU}{dr})=0$$
其解为
$$U=-\frac{c_1}{r}+c_2$$
一般令无穷远处 $U=0$，则 $c_2=0$ 。为了求出 $c_1$ ，将方程包含原点的区域（不妨取以 $ϵ$ 为半径的球 $V_{ϵ}$）进行体积分，利用 $\delta$ 函数的性质
$$\underset{V_{ϵ}}{\iiint }\Delta UdV=1$$
利用格林公式，将上述体积分化为面积分
$$\underset{V_{ϵ}}{\iiint }\Delta UdV=\underset{{\Sigma }_{ϵ}}{\iint }\frac{\partial U}{\partial r}dS={\int }_0^{2\pi }{\int }_0^{\pi }{r}^2\sin \theta d\theta d\varphi =4\pi c_1$$
于是 $c_1=\frac{1}{4\pi }$ ，从而三维泊松方程基本解
$$U(\mathbf{r})=-\frac{1}{4\pi r}\text{\hspace{1em}(1.7)}$$
类似的，用平面极坐标可求得二维泊松方程的基本解
$$U(\mathbf{r})=-\frac{1}{2\pi }\ln \frac{1}{r}=\frac{1}{2\pi }\ln r\text{\hspace{1em}(1.8)}$$

泊松方程的格林函数

泊松方程的边值问题：泊松方程第一、第二、第三类边值问题可统一表示为
$$\{\begin{array}{ll}\Delta u=f(\mathbf{r})& (\mathbf{r}\in V)\\ (\alpha \frac{\partial u}{\partial n}+\beta u){|}_{\Sigma }=\varphi (\mathbf{r})& (\mathbf{r}\in \Sigma )\end{array}\text{\hspace{1em}(2.1)}$$
若 $\alpha =0,\beta \ne 0$ 为第一类边值问题；若 $\alpha \ne 0,\beta =0$ 为第二类边值问题；若 $\alpha \ne 0,\beta \ne 0$ 为第三类边值问题。由上节知道格林函数满足方程
$$\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\text{\hspace{1em}(2.2)}$$
从物理上看，格林函数就是位于 ${\mathbf{r}}_0$点电荷量为 $-{\epsilon }_0$ 的电荷在区域 $V$ 内 $\mathbf{r}$ 点产生的电势。
现在，我们开始使用格林公式，叠加出泊松方程边值问题的解。为此，我们将 (2.1) 中泊松方程和方程 (2.2) 分别乘上 $G(\mathbf{r},{\mathbf{r}}_0)$ 和 $u(\mathbf{r})$ ，相减，然后在区域 $V$ 内积分，得到
$$\underset{V}{\iiint }(G\Delta u-u\Delta G)dV=\underset{V}{\iiint }GfdV-\underset{V}{\iiint }u(\mathbf{r})\delta (\mathbf{r}-{\mathbf{r}}_0)dV$$
根据格林公式，可以将上式左端化为面积分
$$\underset{V}{\iiint }(G\Delta u-u\Delta G)dV=\underset{\Sigma }{\iint }(G\frac{\partial u}{\partial n}-u\frac{\partial G}{\partial n})\mathrm{d}S$$
右端第二项根据 $\delta$ 函数的性质可以得到
$$\underset{V}{\iiint }u(\mathbf{r})\delta (\mathbf{r}-{\mathbf{r}}_0)dV=u({\mathbf{r}}_0)$$
于是，可以得到
$$u({\mathbf{r}}_0)=\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f(\mathbf{r})dV-\underset{\Sigma }{\iint }[G(\mathbf{r},{\mathbf{r}}_0)\frac{\partial u(\mathbf{r})}{\partial n}-u(\mathbf{r})\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}]\mathrm{d}S\text{\hspace{1em}(2.3)}$$
上式称为泊松方程的基本积分公式。
需要注意的是，$G(\mathbf{r},{\mathbf{r}}_0)$ 在 $\mathbf{r}={\mathbf{r}}_0$ 是不连续的，格林公式并不适用。严格的证明是，先在区域 $V$ 内奇点 ${\mathbf{r}}_0$ 处挖去半径为 $\epsilon$ 的球形区域 $V_{\epsilon }$ ，应用格林公式，再令 $\epsilon \to 0$ 取极限求得。今后类似使用时将不再加以说明。

基本积分公式将泊松方程的解 $u$ 用体积分和边界上的面积分表示了出来。对于面积分，$u$ 的边界条件是已知的，如果我们对 $G(\mathbf{r},{\mathbf{r}}_0)$ 提出适当的边界条件，从而获得格林函数确切的解，就可以将 $u$ 确切的表示出来。

(1) 如果泊松方程满足第一类边界条件
$$u{|}_{\Sigma }=\varphi (\mathbf{r})(\mathbf{r}\in \Sigma )$$
同时要求 $G(\mathbf{r},{\mathbf{r}}_0)$ 满足第一类齐次边界条件，即解决边值问题
$$\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\\ G{|}_{\Sigma }=0\end{array}$$
则积分公式 (2.3) 含 $\frac{\partial u}{\partial n}$ 的一项为零，所以不需要知道$\frac{\partial u}{\partial n}$ 在边界上的值，上述边值问题的解称为泊松方程第一边值问题的格林函数。在物理上，可看做边界接地条件下，区域 $V$ 内 ${\mathbf{r}}_0$ 点（点源）电荷为 $-{\epsilon }_0$ 的点源在 $V$ 内 $\mathbf{r}$ 点的电势。基本积分公式
$$u({\mathbf{r}}_0)=\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f(\mathbf{r})dV+\underset{\Sigma }{\iint }\varphi (\mathbf{r})\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}S$$
(2) 如果泊松方程满足第三类边界条件
$$(\alpha \frac{\partial u}{\partial n}+\beta u){|}_{\Sigma }=\varphi (\mathbf{r})(\mathbf{r}\in \Sigma )$$
令 $G(\mathbf{r},{\mathbf{r}}_0)$ 满足第三类齐次边界条件，即解决边界问题
$$\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\\ (\alpha \frac{\partial G}{\partial n}+\beta G){|}_{\Sigma }=0\end{array}$$
分别用 $G,u$ 或 $\frac{\partial G}{\partial n},\frac{\partial u}{\partial n}$ 交叉相乘上述两方程，并相减，可以得到
$$(G\frac{\partial u}{\partial n}-u\frac{\partial G}{\partial n}){|}_{\Sigma }=\frac{1}{\alpha }G\varphi =-\frac{1}{\beta }\frac{\partial G}{\partial n}\varphi$$
上述边值问题的解称为泊松方程第三边值问题的格林函数，此时
$$\begin{array}{rl}u({\mathbf{r}}_0)& =\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f(\mathbf{r})dV-\frac{1}{\alpha }\underset{\Sigma }{\iint }\varphi (\mathbf{r})G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}S\\ & =\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f(\mathbf{r})dV+\frac{1}{\beta }\underset{\Sigma }{\iint }\varphi (\mathbf{r})\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}S\end{array}$$
(3) 至于第二类边界条件
$$\frac{\partial u}{\partial n}{|}_{\Sigma }=\varphi (\mathbf{r})(\mathbf{r}\in \Sigma )$$
似乎可以按照上面的方法，即解决边值问题
$$\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\\ \frac{\partial G}{\partial n}{|}_{\Sigma }=0\end{array}$$
在格林公式中，令 $u(\mathbf{r})=1,v(\mathbf{r})=G(\mathbf{r},{\mathbf{r}}_0)$ ，则有
$$\underset{V}{\iiint }\Delta G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}V=\underset{\Sigma }{\iint }\nabla G(\mathbf{r},{\mathbf{r}}_0)\cdot \mathrm{d}\mathbf{S}=\underset{\Sigma }{\iint }\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}S$$
对边值问题中方程进行积分，根据 $\delta$ 函数的性质又得到
$$\underset{V}{\iiint }\Delta G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}V=1$$
于是格林函数在边界上的积分必须满足
$$\underset{\Sigma }{\iint }\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}S=1\ne 0$$
这显然和上述边值问题中边界条件是矛盾的，第二类齐次边界问题一定无解。此时，需要引进广义的格林函数
$$\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)-\frac{1}{v}\\ \frac{\partial G}{\partial n}{|}_{\Sigma }=0\end{array}$$
其中 $v$ 是区域 $V$ 的体积，并且当且仅当
$$\underset{V}{\iiint }f(\mathbf{r})dV=-\underset{\Sigma }{\iint }\varphi (\mathbf{r})\mathrm{d}S$$
相应边值问题解的积分公式为
$$u({\mathbf{r}}_0)=\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f(\mathbf{r})dV-\underset{\Sigma }{\iint }\varphi (\mathbf{r})G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}S$$
格林函数的对称性：上述积分公式，似乎没有明确的物理意义。接下来我们先讨论格林函数的一个重要性质
$$G({\mathbf{r}}_1,{\mathbf{r}}_2)=G({\mathbf{r}}_2,{\mathbf{r}}_1)\text{\hspace{1em}(2.4)}$$
引入两个格林函数 $G(\mathbf{r},{\mathbf{r}}_1),G(\mathbf{r},{\mathbf{r}}_2)$，以第一类边界问题为例
$$\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_1)=\delta (\mathbf{r}-{\mathbf{r}}_1)\\ G(\mathbf{r},{\mathbf{r}}_1){|}_{\Sigma }=0\end{array},\{\begin{array}{l}\Delta G(\mathbf{r},{\mathbf{r}}_2)=\delta (\mathbf{r}-{\mathbf{r}}_2)\\ G(\mathbf{r},{\mathbf{r}}_2){|}_{\Sigma }=0\end{array}(\mathbf{r},{\mathbf{r}}_1,{\mathbf{r}}_2\in V)$$
利用 $\delta$ 函数的性质和格林公式，上述方程可以得到
$$\begin{array}{rl}& G({\mathbf{r}}_1,{\mathbf{r}}_2)-G({\mathbf{r}}_2,{\mathbf{r}}_1)\\ =& \underset{V}{\iiint }[G(\mathbf{r},{\mathbf{r}}_2)\delta (\mathbf{r}-{\mathbf{r}}_1)-G(\mathbf{r},{\mathbf{r}}_1)\delta (\mathbf{r}-{\mathbf{r}}_2)]\mathrm{d}V\\ =& \underset{V}{\iiint }[G(\mathbf{r},{\mathbf{r}}_2)\Delta G(\mathbf{r},{\mathbf{r}}_1)-G(\mathbf{r},{\mathbf{r}}_1)\Delta G(\mathbf{r},{\mathbf{r}}_2)]\mathrm{d}V\\ =& \underset{\Sigma }{\iint }[G(\mathbf{r},{\mathbf{r}}_2)\nabla G(\mathbf{r},{\mathbf{r}}_1)-G(\mathbf{r},{\mathbf{r}}_1)\nabla G(\mathbf{r},{\mathbf{r}}_2)]\cdot \mathrm{d}\mathbf{S}\end{array}$$
带入边界条件，可得出面积分等于零，于是
$$G({\mathbf{r}}_1,{\mathbf{r}}_2)=G({\mathbf{r}}_2,{\mathbf{r}}_1)$$
对于第二类、第三类边界条件也可以得到同样的结果。

综上所述：对于泊松方程积分公式中的 $\mathbf{r}$ 和 ${\mathbf{r}}_0$ 互换下位置，并利用格林函数的对称性可得
第一边值问题解的积分表达式为
$$u(\mathbf{r})=\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f({\mathbf{r}}_0)\mathrm{d}{V}_0+\underset{\Sigma }{\iint }\varphi ({\mathbf{r}}_0)\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}{S}_0$$
第二边值问题解的积分表达式为
$$u(\mathbf{r})=\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f({\mathbf{r}}_0)\mathrm{d}{V}_0-\underset{\Sigma }{\iint }\varphi ({\mathbf{r}}_0)G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}{S}_0$$
第三边值问题解的积分表达式为
$$\begin{array}{rl}u(\mathbf{r})& =\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f({\mathbf{r}}_0)\mathrm{d}{V}_0-\frac{1}{\alpha }\underset{\Sigma }{\iint }\varphi ({\mathbf{r}}_0)G(\mathbf{r},{\mathbf{r}}_0)\mathrm{d}{S}_0\\ & =\underset{V}{\iiint }G(\mathbf{r},{\mathbf{r}}_0)f({\mathbf{r}}_0)d{V}_0+\frac{1}{\beta }\underset{\Sigma }{\iint }\varphi ({\mathbf{r}}_0)\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}{S}_0\end{array}$$
此时，积分公式有明确的物理意义，右边第一个积分表示在区域 $V$ 内分布的点源在 $\mathbf{r}$ 处产生的场的总和，而第二项则代表边界面上感生场对 $\mathbf{r}$ 处场的影响的总和。

用镜像法求格林函数

泊松方程第一、第三类边值问题对应的格林函数满足
$$\{\begin{array}{ll}\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)& (\mathbf{r}\in V)\\ (\alpha \frac{\partial G}{\partial n}+\beta G){|}_{\Sigma }=0& (\mathbf{r}\in \Sigma )\end{array}\text{\hspace{1em}(3.1)}$$
下面介绍格林函数的两种解法

• 第一种是按相应齐次问题及边界条件的本征函数展开，用分离变量法求得，但这样得到的解往往是无穷级数。
• 格林函数的物理意义启发我们，对于某些特殊区域，格林函数可以通过镜像法求得，可以取得有限形式的解。

镜像法：例如泊松方程第一边值问题的格林函数
$$\{\begin{array}{ll}\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)& (\mathbf{r}\in V)\\ G{|}_{\Sigma }=0& (\mathbf{r}\in \Sigma )\end{array}$$
在物理上可理解为，一接地导体 $V$ 内 ${\mathbf{r}}_0$ 点电荷量为 $-{\epsilon }_0$ 的点电荷在 $V$ 内 $\mathbf{r}$ 点的电势。
由叠加原理，通常将格林函数 $G$ 分成两部分
$$G(\mathbf{r},{\mathbf{r}}_0)=G_0(\mathbf{r},{\mathbf{r}}_0)+G_1(\mathbf{r},{\mathbf{r}}_0)$$
其中 $U$ 满足
$$\Delta G_0=\delta (\mathbf{r}-{\mathbf{r}}_0)$$
是 ${\mathbf{r}}_0$ 点的点电荷产生的场，$G_1$ 满足
$$\{\begin{array}{l}\Delta G_1=0\\ G_1{|}_{\Sigma }=-G_0{|}_{\Sigma }\end{array}$$
是导体内 ${\mathbf{r}}_0$ 点的点电荷在边界上的感应电荷产生的场。
利用第一节中的基本解可知，在三维情形下
$$G_0(\mathbf{r},{\mathbf{r}}_0)=-\frac{1}{4\pi }\frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}\text{\hspace{1em}(3.2)}$$
类似的，二维情形下
$$G_0(\mathbf{r},{\mathbf{r}}_0)=-\frac{1}{2\pi }\ln \frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}=\frac{1}{2\pi }\ln |\mathbf{r}-{\mathbf{r}}_0|\text{\hspace{1em}(3.3)}$$
由于区域 $V$ 外的电源在 $V$ 内产生的场满足拉普拉斯方程，镜像法的中心思想是把边界上的感生电荷用一个等价的点电荷（像电荷）代替，困难在于 $V$ 内点电荷的电场在边界上必须和像电荷的电场相抵消，只有在某些特殊区域（例如，球形，半无界空间，等等）才能实现。

求球内泊松方程第一边值问题格林函数：
$$\{\begin{array}{ll}\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)& (0<r,r_0<a)\\ G{|}_{r=a}=0& \end{array}$$
(1) 像电荷如果存在的话，一定在球外。这是由于感应电荷的电势在球内是处处连续的，在球内的任何电荷都不能产生同样的效果。
(2) 考虑到对称性，这个像电荷一定存在于真实电荷所在半径的延长线上。
记球内电荷位于点 $Q_0({\mathbf{r}}_0)$ ，像电荷位于点 $Q_1({\mathbf{r}}_1)$ 电量为 $q$ ，如图

在球内任取一点 $P(\mathbf{r})$ ，由叠加原理知道，总电势是球内电荷产生的电势和像电荷产生的电势叠加
$$G=-\frac{1}{4\pi }\frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}+\frac{q}{4\pi {\epsilon }_0}\frac{1}{|\mathbf{r}-{\mathbf{r}}_1|}$$
引入球坐标系（原点在球心），由于 ${\mathbf{r}}_0,{\mathbf{r}}_1$ 共线，设 ${\mathbf{r}}_1=\lambda {\mathbf{r}}_0$ 则
$$\begin{gathered} |\mathbf{r}-{\mathbf{r}}_0|=P{Q}_0=\sqrt{r^2-2r{r}_0\cos \gamma +r_0^2} \\ |\mathbf{r}-{\mathbf{r}}_1|=P{Q}_1=\sqrt{r^2-2\lambda r{r}_0\cos \gamma +{\lambda }^2{r}_0^2} \end{gathered}$$
其中 $\gamma$ 是矢径 $\mathbf{r}$ 和 ${\mathbf{r}}_0$ $({\mathbf{r}}_1)$ 之间的夹角
$$\cos \gamma =\cos \theta \cos {\theta }_0+\sin \theta \sin {\theta }_0\cos (\varphi -{\varphi }_0)$$
当观察点 $P$ 位于球面上时，考虑边界条件
$$G{|}_{r=a}=0$$
可得到球面上
$$-\frac{1}{4\pi }\frac{1}{\sqrt{a^2-2a{r}_0\cos \gamma +r_0^2}}+\frac{q}{4\pi {\epsilon }_0}\frac{1}{\sqrt{a^2-2\lambda a{r}_0\cos \gamma +{\lambda }^2{r}_0^2}}=0$$
整理移项得
$$\frac{q}{{\epsilon }_0}\sqrt{a^2-2a{r}_0\cos \gamma +r_0^2}-\sqrt{a^2-2\lambda a{r}_0\cos \gamma +{\lambda }^2{r}_0^2}=0$$
为使上式在球面上恒成立（与球坐标 $\theta ,\varphi$ 无关），可以得到
$$\{\begin{array}{l}\frac{q}{{\epsilon }_0}>0\\ \lambda =(\frac{q}{{\epsilon }_0}{)}^2\\ a^2+{\lambda }^2{r}_0^2=(\frac{q}{{\epsilon }_0}{)}^2(a^2+r_0^2)\end{array}$$
于是我们可得到
$$\{\begin{array}{l}q=\frac{{\epsilon }_{0}a}{r_0}\\ \lambda =(\frac{a}{r_0}{)}^2\end{array}$$
这个设想的等效电荷 $q$ 称为球内点电荷的点像。这样，球内任意一点总电势为
$$G(\mathbf{r},{\mathbf{r}}_0)=-\frac{1}{4\pi }(\frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}-\frac{a}{r_0}\frac{1}{|\mathbf{r}-{\mathbf{r}}_1|})$$
其中 $a$ 为球半径，点像位置 ${\mathbf{r}}_1=(\frac{a}{r_0}{)}^2{\mathbf{r}}_0$
球坐标表示式为
$$G(r,\theta ,\varphi )=-\frac{1}{4\pi }(\frac{1}{\sqrt{r^2-2r{r}_0\cos \gamma +r_0^2}}-\frac{a}{\sqrt{r_0^2{r}^2-2a^{2}r{r}_0\cos \gamma +a^4}})$$
类似的，圆内泊松方程第一边值问题的格林函数
$$\{\begin{array}{l}{\Delta }_{2}G=\delta (\mathbf{r}-{\mathbf{r}}_0)\\ G{|}_{r=a}=0\end{array}$$
用电像法求得其解为
$$G(\mathbf{r},{\mathbf{r}}_0)=-\frac{1}{2\pi }(\ln \frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}-\ln \frac{1}{|\mathbf{r}-{\mathbf{r}}_1|}-\ln \frac{a}{r_0})$$
其中 $a$ 为圆半径，电像位置 ${\mathbf{r}}_1=(\frac{a}{r_0}{)}^2{\mathbf{r}}_0$
极坐标表示式为
$$G(r,\theta )=-\frac{1}{4\pi }[-\ln (r^2-2r{r}_0\cos (\theta -{\theta }_0)+r_0^2)+\ln (r^2-2r\frac{a^2}{r_0}\cos (\theta -{\theta }_0)+\frac{a^4}{r_0^2})+2\ln \frac{a}{r_0}]$$
示例 1：在球内求解拉普拉斯方程的第一边值问题
$$\{\begin{array}{ll}\Delta u=0& (r<a)\\ u{|}_{r=a}=f& \end{array}$$
解：其用格林函数表示的解为
$$u(\mathbf{r})=\underset{\Sigma }{\iint }f({\mathbf{r}}_0)\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}{S}_0$$
其中球内泊松方程第一边值问题的格林函数前面已用电像法求得，球坐标表示式如下
$$G(r,\theta ,\varphi )=-\frac{1}{4\pi }(\frac{1}{\sqrt{r^2-2r{r}_0\cos \gamma +r_0^2}}-\frac{a}{\sqrt{r_0^2{r}^2-2a^{2}r{r}_0\cos \gamma +a^4}})$$
注意到，在球面上外法线方向与 $r_0$ 所在半径的方向一致，因此
$$\frac{\partial G}{\partial n_0}{|}_{r_0=a}=\frac{\partial G}{\partial r_0}{|}_{r_0=a}=\frac{1}{4\pi a}\frac{a^2-r^2}{(a^2-2ar\cos \gamma +r^2{)}^{3/2}}$$
带入第一边值问题解的积分公式
$$u(r,\theta ,\varphi )=\frac{a^2-r^2}{4\pi a}{\int }_0^{2\pi }\mathrm{d}{\varphi }_0{\int }_0^{\pi }\frac{f({\theta }_0,{\varphi }_0)}{(a^2-2ar\cos \gamma +r^2{)}^{3/2}}\sin {\theta }_0\mathrm{d}{\theta }_0$$
上式称为球域上的泊松公式。

示例 2：在圆内解拉普拉斯方程的第一边值问题
$$\{\begin{array}{ll}{u}_{xx}+u_{yy}=0& (r<a)\\ u{|}_{r=a}=f& \end{array}$$
解：和上例用类似的方法可求得
$$u(r,\theta )=\frac{a^2-r^2}{2\pi }{\int }_0^{2\pi }\frac{f({\theta }_0)}{a^2-2ar\cos (\theta -{\theta }_0)+r^2}\mathrm{d}{\theta }_0$$
示例 3：在上半空间 $z>0$ 内求解拉普拉斯方程的第一边值问题
$$\{\begin{array}{ll}\Delta u=0& (z>0)\\ u{|}_{z=0}=f& \end{array}$$
解：其用格林函数表示的解为
$$u(\mathbf{r})=\underset{\Sigma }{\iint }f({\mathbf{r}}_0)\frac{\partial G(\mathbf{r},{\mathbf{r}}_0)}{\partial n}\mathrm{d}{S}_0$$
格林函数 $G(\mathbf{r},{\mathbf{r}}_0)$ 满足的方程为
$$\{\begin{array}{l}\Delta G=\delta (x-x_0)\delta (y-y_0)\delta (z-z_0)\\ G{|}_{z=0}=0\end{array}$$
这相当于接地导体平面 $z=0$ 上方的电势，如图，在点 $Q_0(x_0,y_0,z_0)$ 处放置电荷量为 $-{\epsilon }_0$ 的点电荷。电势可用电像法求得，设想在 $Q_1(x_0,y_0,-z_0)$ 放置电量为 $+{\epsilon }_0$ 的点电荷，不难验证，在 $z=0$ 上电势处处为零，$Q_1$ 即为 $Q_0$ 的电像。

格林函数
$$\begin{array}{rl}G(\mathbf{r},{\mathbf{r}}_0)& =-\frac{1}{4\pi }\frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}+\frac{1}{4\pi }\frac{1}{|\mathbf{r}-{\mathbf{r}}_1|}\\ & =-\frac{1}{4\pi }\frac{1}{\sqrt{(x-x_0{)}^2+(y-y_0{)}^2+(z-z_0{)}^2}}\\ & +\frac{1}{4\pi }\frac{1}{\sqrt{(x-x_0{)}^2+(y-y_0{)}^2+(z+z_0{)}^2}}\end{array}$$
外法向方向与 $z_0$ 正方向一致，因此
$$\frac{\partial G}{\partial n_0}{|}_{z_0=0}=-\frac{\partial G}{\partial z_0}{|}_{z_0=0}=\frac{1}{2\pi }\frac{z}{[(x-x_0{)}^2+(y-y_0{)}^2+z^2{]}^{3/2}}$$
于是可得上半空间的泊松积分
$$u(x,y,z)=\frac{z}{2\pi }\underset{-\infty }{\overset{+\infty }{\iint }}\frac{f(x_0,y_0)}{[(x-x_0{)}^2+(y-y_0{)}^2+z^2{]}^{3/2}}\mathrm{d}{x}_0\mathrm{d}{y}_0$$
示例 4：在圆内解拉普拉斯方程的第一边值问题
$$\{\begin{array}{ll}{u}_{xx}+u_{yy}=0& (y>0)\\ u{|}_{y=0}=f& \end{array}$$
解：和上例用类似的方法可求得
$$u(x,y,z)=\frac{y}{\pi }{\int }_{-\infty }^{+\infty }\frac{f(x_0)}{(x-x_0{)}^2+y^2}\mathrm{d}{x}_0$$

演化问题的基本解

冲量原理：考虑一维受迫振动定解问题
$$\{\begin{array}{ll}\frac{{\partial }^{2}u}{\partial t^2}-a^2\frac{{\partial }^{2}u}{\partial x^2}=f(x,t)& (t>0,0<x<l)\\ u{|}_{x=0}=u{|}_{x=l}=0& \\ u{|}_{t=0}=0,\frac{\partial u}{\partial t}{|}_{t=0}=0& \end{array}$$

这里 $f(x,t)=\frac{F(x,t)}{\rho }$ 是作用在弦单位长度单位质量上的外力。

考虑 $[\tau ,\tau +\Delta \tau )$ 时间段内的位移变化，令 $w{|}_{t=\tau }=0,w{|}_{t=\tau +\Delta \tau }=\Delta u$
$f(x,\tau )\Delta \tau$ 表示 $\Delta \tau$ 内冲量，这个冲量使得系统的速度有一定的增量。现在，我们把 $\Delta \tau$ 时间内的速度增量看成是 $t=\tau$ 瞬时得到的，而在 $\Delta \tau$ 的其余时间内没有冲量的作用，即在这段时间内没有力的作用，故方程是齐次的。$t=\tau$ 时的集中速度可置于“初始”条件中，得到的关于瞬时力引起的振动的定解问题是
$$\{\begin{array}{l}\frac{{\partial }^{2}w}{\partial t^2}-a^2\frac{{\partial }^{2}w}{\partial x^2}=0\\ w{|}_{x=0}=w{|}_{x=l}=0\\ w{|}_{t=\tau }=0,\frac{\partial u}{\partial t}{|}_{t=\tau }=f(x,\tau )\Delta \tau \end{array}(\tau <t<\tau +\Delta \tau ,0<x<l)$$
令 $w(x,t;\tau ,\Delta \tau )=v(x,t;\tau )\Delta \tau$ 则
$$\{\begin{array}{l}\frac{{\partial }^{2}v}{\partial t^2}-a^2\frac{{\partial }^{2}v}{\partial x^2}=0\\ v{|}_{x=0}=v{|}_{x=l}=0\\ v{|}_{t=\tau }=0,\frac{\partial v}{\partial t}{|}_{t=\tau }=f(x,\tau )\end{array}(t>\tau ,0<x<l)$$
于是
$$u(x,t)=\underset{\Delta \tau \to 0}{\lim }\sum _{\tau =0}^{t}w(x,t;\tau )={\int }_0^{t}v(x,t;\tau )\mathrm{d}\tau$$

冲量原理：（齐次化原理）设 $L$ 是关于 $\mathbf{r}=(x_1,x_2,\cdots ,x_n)$ 的线性偏微分算子，若 $v(\mathbf{r},t;\tau )$ 满足齐次方程定解问题
$$\{\begin{array}{l}\frac{{\partial }^{m}v}{\partial t^m}-L[v]=0\\ (\alpha \frac{\partial v}{\partial t}+\beta v){|}_{\Sigma }=0\\ v{|}_{t=\tau }=\frac{\partial v}{\partial t}{|}_{t=\tau }=\cdots =\frac{{\partial }^{m-2}v}{\partial t^{m-2}}{|}_{t=\tau }=0\\ \frac{{\partial }^{m-1}v}{\partial t^{m-1}}{|}_{t=\tau }=f(\mathbf{r},\tau )\end{array}(t>\tau )$$
则
$$u(\mathbf{r},t)={\int }_0^{t}v(\mathbf{r},t;\tau )\mathrm{d}\tau$$
是以下非齐次方程定解问题的解
$$\{\begin{array}{l}\frac{{\partial }^{m}u}{\partial t^m}-L[u]=f(\mathbf{r},t)\\ (\alpha \frac{\partial u}{\partial t}+\beta u){|}_{\Sigma }=0\\ u{|}_{t=0}=\frac{\partial u}{\partial t}{|}_{t=0}=\cdots =\frac{{\partial }^{m-1}u}{\partial t^{m-1}}{|}_{t=0}=0\end{array}(t>0)$$
以无界热传导初值问题为例，若 $v(\mathbf{r},t;\tau )$ 满足齐次方程初值问题
$$\{\begin{array}{ll}\frac{\partial v}{\partial t}-a^2\Delta v=0& (t>\tau )\\ v{|}_{t=\tau }=f(\mathbf{r},\tau )& \end{array}$$
则积分 $u(\mathbf{r},t)={\int }_0^{t}v(\mathbf{r},t;\tau )\mathrm{d}\tau$ 满足非齐次方程初值问题
$$\{\begin{array}{ll}\frac{\partial u}{\partial t}-a^2\Delta u=f(\mathbf{r},t)& (t>0)\\ u{|}_{t=0}=0& \end{array}$$

热传导方程初值问题：对于无界空间的热传导方程初值问题
$$\{\begin{array}{ll}\frac{\partial u}{\partial t}-a^2\Delta u=f(\mathbf{r},t)& (t>0)\\ u{|}_{t=0}=\varphi (\mathbf{r})& \end{array}\text{\hspace{1em}(4.1)}$$
由叠加原理知道，$u=u_1+u_2$ ，其中 $u_1$ 满足初值问题
$$\{\begin{array}{l}\frac{\partial u_1}{\partial t}-a^2\Delta u_1=0(t>0)\\ u_1{|}_{t=0}=\varphi (\mathbf{r})\end{array}$$
代表在初始时刻瞬时给予热量的传导问题，$u_2$ 满足初值问题
$$\{\begin{array}{l}\frac{\partial u_2}{\partial t}-a^2\Delta u_2=f(\mathbf{r},t)(t>0)\\ u_2{|}_{t=0}=0\end{array},$$
代表持续热源下的热传导。由 $\delta$ 函数的性质知道，空间持续热源可看做瞬时点源的叠加
$$f(\mathbf{r},t)=\iiint \mathrm{d}{V}_0\int f({\mathbf{r}}_0,t_0)\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\mathrm{d}{t}_0$$
将空间点 ${\mathbf{r}}_0$ 在 $t_0$ 时刻的热源在 $\mathbf{r},t$ 引起的温度记作 $G(\mathbf{r},t;{\mathbf{r}}_0,t_0)$ 称为热传导方程的格林函数。 $G$ 满足的初值问题
$$\{\begin{array}{l}{G}_t-a^2\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\\ G{|}_{t=t_0}=0\end{array}\text{\hspace{1em}(4.2)}$$
方程中 $t=t_0$ 时刻即代表初始时刻，之前无任何热源作用，$t<t_0$ 时刻 $G$ 均为零。
取 $\tau =t-t_0$ 由冲量原理，若 $v(\mathbf{r},\tau ;{\tau }_0)$ 满足齐次方程初值问题
$$\{\begin{array}{l}{v}_{\tau }-a^2\Delta v=0\\ v{|}_{\tau ={\tau }_0}=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta ({\tau }_0)\end{array}$$

则格林函数 $G={\int }_0^{\tau }v(\mathbf{r},\tau ;{\tau }_0)\mathrm{d}{\tau }_0$
对 $v(\mathbf{r},\tau ;{\tau }_0)$ 满足的齐次方程和初值条件积分可以得到
$$\{\begin{array}{l}{G}_{\tau }-a^2\Delta G=0\\ G{|}_{\tau ={\tau }_0}=\delta (\mathbf{r}-{\mathbf{r}}_0)\end{array}$$
上式取 ${\tau }_0=0$ ，并将 $\tau =t-t_0$ 代回，可知格林函数满足的定解问题 (4.2) 等价于以下定解问题
$$\{\begin{array}{l}{G}_t-a^2\Delta G=0\\ G{|}_{t=t_0}=\delta (\mathbf{r}-{\mathbf{r}}_0)\end{array}$$
由冲量原理和积分叠加原理，进一步可求得
$$u_2(\mathbf{r},t)=\iiint \mathrm{d}{V}_0{\int }_0^{t}G(\mathbf{r},t;{\mathbf{r}}_0,t_0)f({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0$$
为求得 $G(\mathbf{r},t;{\mathbf{r}}_0,t_0)$ ，不妨先求满足定解问题
$$\{\begin{array}{l}{U}_t-a^2\Delta U=0\\ U{|}_{t=0}=\delta (\mathbf{r})\end{array}\text{\hspace{1em}(4.3)}$$
的基本解 $U(\mathbf{r},t)$ ，代表初始瞬间原点给予热量 $Q=c\rho$ 后的温度分布。做变量变换有
$$G(\mathbf{r},t;{\mathbf{r}}_0,t_0)=U(\mathbf{r}-{\mathbf{r}}_0,t-t_0)\text{\hspace{1em}(4.4)}$$
进而有
$$u_2(\mathbf{r},t)=\iiint \mathrm{d}{V}_0{\int }_0^{t}U(\mathbf{r}-{\mathbf{r}}_0,t)f({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0={\int }_0^{t}U(\mathbf{r},t-t_0)\ast f(\mathbf{r},t_0)\mathrm{d}{t}_0$$
由于 $\varphi (\mathbf{r})=\delta (\mathbf{r})\ast \varphi (\mathbf{r})$ ，根据叠加原理有
$$u_1(\mathbf{r},t)=\iiint U(\mathbf{r}-{\mathbf{r}}_0,t)\varphi ({\mathbf{r}}_0)\mathrm{d}{V}_0=U(\mathbf{r},t)\ast \varphi (\mathbf{r})$$
所以
$$u(\mathbf{r},t)=U(\mathbf{r},t)\ast \varphi (\mathbf{r})+{\int }_0^{t}U(\mathbf{r},t-t_0)\ast f(\mathbf{r},t_0)\mathrm{d}{t}_0\text{\hspace{1em}(4.5)}$$
可带入 (4.1) 可直接验证上述结论。

波动方程初值问题：对于无界空间的波动方程初值问题
$$\{\begin{array}{ll}\frac{{\partial }^{2}u}{\partial t^2}-a^2\Delta u=f(\mathbf{r},t)& (t>0)\\ u{|}_{t=0}=\varphi (\mathbf{r}),\frac{\partial u}{\partial t}{|}_{t=0}=\psi (\mathbf{r})& \end{array}\text{\hspace{1em}(4.6)}$$
同样由叠加原理知道，$u=u_1+u_2+u_3$ ，分别满足初值问题
$$\{\begin{array}{l}\frac{{\partial }^2{u}_1}{\partial t^2}-a^2\Delta u_1=f(\mathbf{r},t)(t>0)\\ u_1{|}_{t=0}=\frac{\partial u_1}{\partial t}{|}_{t=0}=0\end{array}$$
代表持续作用的力引起的波动
$$\{\begin{array}{l}\frac{{\partial }^2{u}_2}{\partial t^2}-a^2\Delta u_2=0(t>0)\\ u_2{|}_{t=0}=0,\frac{\partial u_2}{\partial t}{|}_{t=0}=\psi (\mathbf{r})\end{array}$$
代表初始时刻瞬时冲量引起的波动
$$\{\begin{array}{l}\frac{{\partial }^2{u}_3}{\partial t^2}-a^2\Delta u_3=0(t>0)\\ u_3{|}_{t=0}=\varphi (\mathbf{r}),\frac{\partial u_3}{\partial t}{|}_{t=0}=0\end{array}$$
同样林函数 $G(\mathbf{r},t;{\mathbf{r}}_0,t_0)$ 满足的定解问题为
$$\{\begin{array}{l}{G}_{tt}-a^2\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\\ G{|}_{t=t_0}=0,\frac{\partial G}{\partial t}{|}_{t=t_0}=0\end{array}\text{\hspace{1em}(4.7)}$$
方程中 $t=t_0$ 时刻即代表初始时刻，之前无任何力的作用，$t<t_0$ 时刻 $G$ 均为零。
根据冲量原理和积分叠加原理定解问题 (4.7) 等价于
$$\{\begin{array}{l}{G}_{tt}-a^2\Delta G=0\\ G{|}_{t=t_0}=0,\frac{\partial G}{\partial t}{|}_{t=t_0}=\delta (\mathbf{r}-{\mathbf{r}}_0)\end{array}$$
不妨先求 基本解 $U(\mathbf{r},t)$ 满足的定解问题
$$\{\begin{array}{l}{U}_{tt}-a^2\Delta U=0\\ U{|}_{t=0}=0,\frac{\partial U}{\partial t}{|}_{t=0}=\delta (\mathbf{r})\end{array}\text{\hspace{1em}(4.8)}$$
代表初始时刻在原点处的瞬时冲量引起的波动。变量代换可求得
$$G(\mathbf{r},t;{\mathbf{r}}_0,t_0)=U(\mathbf{r}-{\mathbf{r}}_0,t-t_0)\text{\hspace{1em}(4.9)}$$
根据冲量原理和积分叠加原理，进而有
$$\begin{array}{rl}& u_1(\mathbf{r},t)=\iiint \mathrm{d}{V}_0{\int }_0^{t}U(\mathbf{r}-{\mathbf{r}}_0,t-t_0)f({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0={\int }_0^{t}U(\mathbf{r},t-t_0)\ast f(\mathbf{r},t_0)\mathrm{d}{t}_0\\ & u_2(\mathbf{r},t)=\iiint U(\mathbf{r}-{\mathbf{r}}_0,t)\varphi ({\mathbf{r}}_0)\mathrm{d}{V}_0=U(\mathbf{r},t)\ast \psi (\mathbf{r})\end{array}$$

设 $u_3=\frac{\partial v}{\partial t}$ ，则有
$$\{\begin{array}{l}\frac{{\partial }^{2}v}{\partial t^2}-a^2\Delta v=0(t>0)\\ v{|}_{t=0}=0,\frac{\partial v}{\partial t}{|}_{t=0}=\varphi (\mathbf{r})\end{array}$$
即 $v=U(\mathbf{r},t)\ast \varphi (\mathbf{r})$ ， $u_3=\frac{\partial v}{\partial t}$ 满足 $u_3$ 的初值问题，所以
$$u(\mathbf{r},t)=U(\mathbf{r},t)\ast \psi (\mathbf{r})+\frac{\partial }{\partial t}[U(\mathbf{r},t)\ast \varphi (\mathbf{r})]+{\int }_0^{t}U(\mathbf{r},t-t_0)\ast f(\mathbf{r},t_0)\mathrm{d}{t}_0\text{\hspace{1em}(4.10)}$$
基本解的求法： 基本解 $U(\mathbf{r},t)$ 没有边界条件限制，因此不是惟一的，适当选取即可。

三维无界热传导方程初值问题的基本解：（可用傅里叶变换法求得）
$$U(\mathbf{r},t)=(\frac{1}{2a\sqrt{\pi t}}{)}^3\exp (-\frac{|\mathbf{r}{|}^2}{4a^{2}t})$$

三维无界波动方程初值问题的基本解：（可用傅里叶变换法求得）
$$U(\mathbf{r},t)=(\frac{1}{2\pi }{)}^3\underset{-\infty }{\overset{-\infty }{\iiint }}\frac{\sin (r_{0}at)}{r_{0}a}\exp (\mathrm{i}{\mathbf{r}}_0\cdot \mathbf{r})\mathrm{d}{V}_0$$

其中 $r_0=|{\mathbf{r}}_0|,r=|\mathbf{r}|$ 。以 $\mathbf{r}$ 为极轴方向取球坐标，则
$$U(\mathbf{r},t)=\frac{\delta (r-at)}{4\pi ar}$$

一般演化问题的格林函数

波动方程定解问题：一般强迫振动波动方程的定解问题
$$\{\begin{array}{ll}{u}_{tt}-a^2\Delta u=f(\mathbf{r},t)& (t>0,\mathbf{r}\in V)\\ (\alpha \frac{\partial u}{\partial n}+\beta u){|}_{\Sigma }=\sigma (\mathbf{r},t)& (t>0,\mathbf{r}\in \Sigma )\\ u{|}_{t=0}=\varphi (\mathbf{r}),\frac{\partial u}{\partial t}{|}_{t=0}=\psi (\mathbf{r},)& (t>0,\mathbf{r}\in V)\end{array}\text{\hspace{1em}(5.1)}$$
由 $\delta$ 函数的性质知道，持续力 $f(\mathbf{r},t)$ 可表示为
$$f(\mathbf{r},t)=\underset{V}{\iiint }\mathrm{d}{V}_0{\int }_{t}f({\mathbf{r}}_0,t_0)\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\mathrm{d}{t}_0$$
波动方程的格林函数$G(\mathbf{r},t;{\mathbf{r}}_0,t_0)$ 满足的定解问题是
$$\{\begin{array}{ll}{G}_{tt}-a^2\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)& (t,t_0>0,\mathbf{r}\in V)\\ (\alpha \frac{\partial G}{\partial n}+\beta G){|}_{\Sigma }=0& (t,t_0>0,\mathbf{r}\in \Sigma )\\ G{|}_{t=t_0}=0,\frac{\partial G}{\partial t}{|}_{t=t_0}=0& (t,t_0>0,\mathbf{r}\in V)\end{array}\text{\hspace{1em}(5.2)}$$
方程中 $t=t_0$ 时刻即代表初始时刻，$G{|}_{t<t_0}\equiv 0$ 。我们可以用和解泊松方程类似的方法求解波动方程解的积分表达式，首先讨论格林函数的对称性。

格林函数的对称性
$$G({\mathbf{r}}_1,t_1;{\mathbf{r}}_2,t_2)=G({\mathbf{r}}_2,-t_2;{\mathbf{r}}_1,-t_1)\text{\hspace{1em}(5.3)}$$
引入两个格林函数 $G(\mathbf{r},t;{\mathbf{r}}_1,t_1),G(\mathbf{r},-t;{\mathbf{r}}_2,-t_2)$，简记为 $G_1,G_2$ ，分别是下面定解问题的解
$$\{\begin{array}{l}\frac{{\partial }^2{G}_1}{\partial t^2}-a^2\Delta G_1=\delta (\mathbf{r}-{\mathbf{r}}_1)\delta (t-t_1)\\ (\alpha \frac{\partial G_1}{\partial n}+\beta G_1){|}_{\Sigma }=0\\ G_1{|}_{t=t_1}=0,\frac{\partial G_1}{\partial t}{|}_{t=t_1}=0\end{array},\{\begin{array}{l}\frac{{\partial }^2{G}_2}{\partial (-t{)}^2}-a^2\Delta G_2=\delta (\mathbf{r}-{\mathbf{r}}_2)\delta (t+t_2)\\ (\alpha \frac{\partial G_2}{\partial n}+\beta G_2){|}_{\Sigma }=0\\ G_2{|}_{-t=-t_2}=0,\frac{\partial G_2}{\partial t}{|}_{-t=-t_2}=0\end{array}$$
其中 ${\mathbf{r}}_1,{\mathbf{r}}_2\in V,t>t_1,t_2>0$ 。
利用 $\delta$ 函数的性质和格林公式，在空间区域 $V$ 和 时间区间$[0,t]$ 上积分，上述方程可以得到
$$\begin{array}{rl}& G({\mathbf{r}}_1,t_1;{\mathbf{r}}_2,t_2)-G({\mathbf{r}}_2,-t_2;{\mathbf{r}}_1,-t_1)\\ =& \underset{V}{\iiint }\mathrm{d}V{\int }_0^t[G_2\delta (\mathbf{r}-{\mathbf{r}}_1)\delta (t-t_1)-G_1\delta (\mathbf{r}-{\mathbf{r}}_2)\delta (t+t_2)]\mathrm{d}t\\ =& \underset{V}{\iiint }\mathrm{d}V{\int }_0^t[G_2(\frac{{\partial }^2{G}_1}{\partial t^2}-a^2\Delta G_1)-G_1(\frac{{\partial }^2{G}_2}{\partial (-t{)}^2}-a^2\Delta G_2)]\mathrm{d}t\\ =& \underset{V}{\iiint }\mathrm{d}V{\int }_0^t[(G_2\frac{{\partial }^2{G}_1}{\partial t^2}-G_1\frac{{\partial }^2{G}_2}{\partial t^2})+a^2(G_1\Delta G_2-G_2\Delta G_1)]\mathrm{d}t\\ =& \underset{V}{\iiint }\mathrm{d}V[G_2\frac{\partial G_1}{\partial t}-G_1\frac{\partial G_2}{\partial t}]{|}_0^t+a^2{\int }_0^t\mathrm{d}t\underset{\Sigma }{\iint }(G_1\frac{\partial G_2}{\partial n}-G_2\frac{\partial G_1}{\partial n})\mathrm{d}S\end{array}$$
将上述两个边界条件分别乘以 $G_2,G_1$ ，相减可以得到
$$(G_1\frac{\partial G_2}{\partial n}-G_2\frac{\partial G_1}{\partial n}){|}_{\Sigma }=0$$
带入初始条件我们又可以得到
$$[G_2\frac{\partial G_1}{\partial t}-G_1\frac{\partial G_2}{\partial t}]{|}_0^t=0$$
于是有
$$G({\mathbf{r}}_1,t_1;{\mathbf{r}}_2,t_2)=G({\mathbf{r}}_2,-t_2;{\mathbf{r}}_1,-t_1)$$

tips：

1. $\frac{\partial G}{\partial t}=\frac{\partial G}{\partial (-t)}\frac{\mathrm{d}(-t)}{\mathrm{d}t}=-\frac{\partial G}{\partial (-t)}$
   将 $t_1,t_2$ 位置互换时出现的负号，正好保证了时间的先后次序不变，否则就会有悖于因果律的要求。
2. $\frac{{\partial }^{2}G}{\partial t^2}=\frac{{\partial }^{2}G}{\partial (-t{)}^2}$
3. 波动方程中重要的偏微分
   $\begin{array}{rl}& \frac{\partial }{\partial t}(G_2\frac{\partial G_1}{\partial t}-G_1\frac{\partial G_2}{\partial t})\\ =& (\frac{\partial G_2}{\partial t}\frac{\partial G_1}{\partial t}+G_2\frac{{\partial }^2{G}_1}{\partial t^2})-(\frac{\partial G_1}{\partial t}\frac{\partial G_2}{\partial t}+G_1\frac{{\partial }^2{G}_2}{\partial t^2})\\ =& G_2\frac{{\partial }^2{G}_1}{\partial t^2}-G_1\frac{{\partial }^2{G}_2}{\partial t^2}\end{array}$
4. 热传导方程中重要的偏微分
   $G_2\frac{\partial G_1}{\partial t}+G_1\frac{\partial G_2}{\partial t}=\frac{\partial (G_1{G}_2)}{\partial t}$

解的积分表达式：将波动方程定解问题中的 $\mathbf{r},t$ 改写成 ${\mathbf{r}}_0,t_0$
$$\{\begin{array}{l}\frac{{\partial }^{2}u({\mathbf{r}}_0,t_0)}{\partial t_0^2}-a^2{\Delta }_{0}u({\mathbf{r}}_0,t_0)=f({\mathbf{r}}_0,t_0)\\ [\alpha \frac{\partial u({\mathbf{r}}_0,t_0)}{\partial n_0}+\beta u({\mathbf{r}}_0,t_0)]{|}_{\Sigma }=\sigma ({\mathbf{r}}_0,t_0)\\ u({\mathbf{r}}_0,t_0){|}_{t_0=0}=\varphi ({\mathbf{r}}_0),\frac{\partial u({\mathbf{r}}_0,t_0)}{\partial t_0}{|}_{t_0=0}=\psi ({\mathbf{r}}_0)\end{array}$$
再将格林函数定解问题中的 $\mathbf{r},t$ 改换成 ${\mathbf{r}}_0,-t_0$ ，将 ${\mathbf{r}}_0,t_0$ 改换成 $\mathbf{r},-t$ 同时利用对称关系，得
$$\{\begin{array}{l}\frac{{\partial }^{2}G(\mathbf{r},t;{\mathbf{r}}_0,t_0)}{\partial t_0^2}-a^2{\Delta }_{0}G(\mathbf{r},t;{\mathbf{r}}_0,t_0)=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\\ [\alpha \frac{\partial G(\mathbf{r},t;{\mathbf{r}}_0,t_0)}{\partial n_0}+\beta G(\mathbf{r},t;{\mathbf{r}}_0,t_0)]{|}_{\Sigma }=0\\ G(\mathbf{r},t;{\mathbf{r}}_0,t_0){|}_{-t_0=-t}=0,\frac{\partial G(\mathbf{r},t;{\mathbf{r}}_0,t_0)}{\partial t_0}{|}_{-t_0=-t}=0\end{array}$$
两方程交叉相乘 $G(\mathbf{r},t;{\mathbf{r}}_0,t_0),u({\mathbf{r}}_0,t_0)$ 相减，再积分
$$\begin{array}{rl}& \underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t+0}[G\frac{{\partial }^{2}u({\mathbf{r}}_0,t_0)}{\partial t_0^2}-u({\mathbf{r}}_0,t_0)\frac{{\partial }^{2}G}{\partial t_0^2}]\mathrm{d}{t}_0\\ & -a^2\underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t+0}[G{\Delta }_{0}u({\mathbf{r}}_0,t_0)-u({\mathbf{r}}_0,t_0){\Delta }_{0}G]\mathrm{d}{t}_0\\ =& \underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t+0}Gf({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0\\ & -\underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t+0}u({\mathbf{r}}_0,t_0)\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)\mathrm{d}{t}_0\end{array}$$
利用 $\delta$ 函数的性质和格林公式可得到
$$\begin{array}{rl}u(\mathbf{r},t)=& \underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t+0}Gf({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0\\ & -\underset{V}{\iiint }[G\frac{\partial u({\mathbf{r}}_0,t_0)}{\partial t_0}-u({\mathbf{r}}_0,t_0)\frac{\partial G}{\partial t_0}]{|}_{t_0=0}^{t_0=t+0}\mathrm{d}{V}_0\\ & +a^2{\int }_0^{t+0}\mathrm{d}{t}_0\underset{\Sigma }{\iint }[G\frac{\partial u({\mathbf{r}}_0,t_0)}{\partial n_0}-u({\mathbf{r}}_0,t_0)\frac{\partial G}{\partial n_0}]\mathrm{d}{S}_0\end{array}$$
带入初始条件和边界条件
$$\begin{array}{rl}u(\mathbf{r},t)=& \underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t}Gf({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0\\ & +\underset{V}{\iiint }[G\psi ({\mathbf{r}}_0)-\varphi ({\mathbf{r}}_0)\frac{\partial G}{\partial t_0}]{|}_{t_0=0}\mathrm{d}{V}_0\\ & +a^2{\int }_0^t\mathrm{d}{t}_0\underset{\Sigma }{\iint }[G\frac{\partial u({\mathbf{r}}_0,t_0)}{\partial n_0}-u({\mathbf{r}}_0,t_0)\frac{\partial G}{\partial n_0}]\mathrm{d}{S}_0\end{array}\text{\hspace{1em}(5.4)}$$
对于不同类型的边界条件条件，可令 $G$ 满足相应的齐次边界条件，从而得到适用于不同边界条件的解以 $G$ 表示的解的积分表达式。

热传导方程定解问题
$$\{\begin{array}{ll}{u}_t-a^2\Delta u=f(\mathbf{r},t)& (t>0,\mathbf{r}\in V)\\ (\alpha \frac{\partial u}{\partial n}+\beta u){|}_{\Sigma }=\sigma (\mathbf{r},t)& (t>0,\mathbf{r}\in \Sigma )\\ u{|}_{t=0}=\varphi (\mathbf{r})& (t>0,\mathbf{r}\in V)\end{array}\text{\hspace{1em}(5.5)}$$
格林函数满足的定解问题是
$$\{\begin{array}{ll}{G}_t-a^2\Delta G=\delta (\mathbf{r}-{\mathbf{r}}_0)\delta (t-t_0)& (t,t_0>0,\mathbf{r}\in V)\\ (\alpha \frac{\partial G}{\partial n}+\beta G){|}_{\Sigma }=0& (t,t_0>0,\mathbf{r}\in \Sigma )\\ G{|}_{t=t_0}=0& (t,t_0>0,\mathbf{r}\in V)\end{array}\text{\hspace{1em}(5.6)}$$
类似上面的讨论，同样可以得到解的积分表达式
$$\begin{array}{rl}u(\mathbf{r},t)=& \underset{V}{\iiint }\mathrm{d}{V}_0{\int }_0^{t}Gf({\mathbf{r}}_0,t_0)\mathrm{d}{t}_0\\ & +\underset{V}{\iiint }[G\varphi ({\mathbf{r}}_0)]{|}_{t_0=0}\mathrm{d}{V}_0\\ & +a^2{\int }_0^t\mathrm{d}{t}_0\underset{\Sigma }{\iint }[G\frac{\partial u({\mathbf{r}}_0,t_0)}{\partial n_0}-u({\mathbf{r}}_0,t_0)\frac{\partial G}{\partial n_0}]\mathrm{d}{S}_0\end{array}\text{\hspace{1em}(5.7)}$$

格林函数的求法

示例 1：对于一维热传导问题的格林函数
$$\{\begin{array}{ll}{G}_t-a^2{G}_{xx}=\delta (x-x_0)\delta (t-t_0)& (t,t_0>0,0<x<l)\\ G{|}_{x=0}=G{|}_{x=l}=0& \\ G{|}_{t=t_0}=0& \end{array}$$
由特征函数展开，分离变量法可求得
$$G(x,t;x_0,t_0)=\frac{l}{2}\sum _{n=1}^{\infty }\exp [-(\frac{n\pi a}{l}{)}^2(t-t_0)]\sin \frac{n\pi x_0}{l}\sin \frac{n\pi x}{l}(t>t_0)$$
示例 2：对于一维受迫振动的格林函数
$$\{\begin{array}{ll}{G}_t-a^2{G}_{xx}=\delta (x-x_0)\delta (t-t_0)& (t,t_0>0,0<x<l)\\ G{|}_{x=0}=G{|}_{x=l}=0& \\ G{|}_{t=t_0}=0,\frac{\partial G}{\partial t}{|}_{t=t_0}=0& \end{array}$$
由特征函数展开，分离变量法可求得
$$G(x,t;x_0,t_0)=\frac{2}{\pi a}\sum _{n=1}^{\infty }\frac{1}{n}\sin \frac{n\pi x_0}{l}\sin \frac{n\pi a(t-t_0)}{l}\sin \frac{n\pi x}{l}(t>t_0)$$

附录

静电场理论

库伦定律
$$\mathbf{F}=\frac{1}{4\pi {\epsilon }_0}\frac{Qq}{r^2}{\mathbf{e}}_r\text{\hspace{1em}(1.1)}$$
电场强度：是用来表示电场的强弱和方向的物理量
$$\mathbf{E}=\frac{\mathbf{F}}{q}=\frac{1}{4\pi {\epsilon }_0}\frac{Q}{r^2}{\mathbf{e}}_r\text{\hspace{1em}(1.2)}$$
高斯定理：穿过闭合曲面 $\Sigma$ 向外的电场强度通量等于闭合曲面围成的空间 $V$ 内电量的 $1/{\epsilon }_0$，其中 ${\epsilon }_0$ 为真空介电常数。
$$\underset{\Sigma }{\iint }\mathbf{E}\cdot \mathrm{d}\mathbf{S}=\frac{1}{{\epsilon }_0}\underset{V}{\iiint }\rho \mathrm{d}V\text{\hspace{1em}(1.3)}$$
简要证明：先对点电荷的场证明，再推广到一半的电荷分布。

(1) 取包围点电荷 $q$ 的任意闭合曲面 $\Sigma$，在闭合曲面上任取面元 $\mathrm{d}\mathbf{S}$ ，由库伦定律计算电通量
$$\mathrm{d}\Phi =\mathbf{E}\cdot \mathrm{d}\mathbf{S}=\frac{1}{4\pi {\epsilon }_0}\frac{q}{r^2}\mathrm{d}S\cos \theta$$
面元 $\mathrm{d}\mathbf{S}$ 在球面上（以点电荷为球心，$r$ 为半径）的投影 $\mathrm{d}{S}_0=\mathrm{d}S\cos \theta$ ，在此，引入球面立体角³ $\mathrm{d}\Omega =\frac{\mathrm{d}{S}_0}{r^2}$ ，然后对上式积分
$$\underset{\Sigma }{\iint }\mathbf{E}\cdot \mathrm{d}\mathbf{S}=\frac{q}{4\pi {\epsilon }_0}\underset{\Sigma }{\iint }\mathrm{d}\Omega =\frac{q}{{\epsilon }_0}$$
(2) 根据场强叠加原理，上式可进一步扩展为
$$\underset{\Sigma }{\iint }\mathbf{E}\cdot \mathrm{d}\mathbf{S}=\frac{1}{{\epsilon }_0}\sum _{q_i\text{ in }V}{q}_i$$
若电荷为连续分布，电荷密度为 $\rho$ 则上式可改写为
$$\underset{\Sigma }{\iint }\mathbf{E}\cdot \mathrm{d}\mathbf{S}=\frac{1}{{\epsilon }_0}\underset{V}{\iiint }\rho \mathrm{d}V$$
通过数学上的高斯公式，上式左端可化为体积分
$$\underset{V}{\iiint }\nabla \cdot \mathbf{E}\mathrm{d}V=\frac{1}{{\epsilon }_0}\underset{V}{\iiint }\rho \mathrm{d}V$$
由于闭合曲面的随意性，取极限 $V\to M(x,y,z)$ 上式可化为微分形式
$$\nabla \cdot \mathbf{E}=\frac{\rho }{{\epsilon }_0}\text{\hspace{1em}(1.4)}$$
静电场环路定理： 在静电场中，场强沿任意闭合路径的线积分等于 0
$${\oint }_L\mathbf{E}\cdot \mathrm{d}l=0\text{\hspace{1em}(1.5)}$$
取点电荷 $q$ 的电场
$$\mathbf{E}\cdot \mathrm{d}l=\frac{1}{4\pi {\epsilon }_0}\frac{q}{r^2}\mathrm{d}l\cos =\frac{1}{4\pi {\epsilon }_0}\frac{q}{r^2}\mathrm{d}r$$
环路积分
$${\oint }_L\mathbf{E}\cdot \mathrm{d}l=\frac{q}{4\pi {\epsilon }_0}{\oint }_L\frac{1}{r^2}\mathrm{d}r=0$$
根据场强叠加原理，可推广至多个点电荷及电荷连续分布情形。
根据数学上的斯托克斯公式，上式左端可以化为面积分
$$\underset{S}{\iint }(\nabla \times \mathbf{E})\cdot \mathrm{d}\mathbf{S}=0$$
其中 $S$ 为以 $L$ 为环边界线的任意曲面。取极限 $S\to M(x,y,z)$ 上式可化为微分形式
$$\nabla \times \mathbf{E}=0\text{\hspace{1em}(1.6)}$$

电势能和电势：由静电场环路定理知道
$${\oint }_L\mathbf{F}\cdot \mathrm{d}l={\oint }_{L}q\mathbf{E}\cdot \mathrm{d}l=0$$
由此可见，电场力做功与路径无关，只和起始和终止位置有关，由此引入势能的概念。
定义：静电场力做功等于相应电势能的减小量
$${\int }_A^B\mathbf{F}\cdot \mathrm{d}l=W_A-W_B$$
做功的大小还和电荷量有关，在此引入电势 ${\phi }_P=\frac{W_P}{q}$ ，电势差与路径无关
$${\int }_A^B\mathbf{E}\cdot \mathrm{d}l={\phi }_A-{\phi }_B$$
取无穷远处为电势零点，可求得点电荷电势场方程为
$$\phi (\mathbf{r})={\int }_r^{\infty }\mathbf{E}\cdot \mathrm{d}l={\int }_r^{\infty }\frac{1}{4\pi {\epsilon }_0}\frac{Q}{r^2}\mathrm{d}r=\frac{Q}{4\pi {\epsilon }_{0}r}$$
根据电场叠加原理进一步推广到电荷连续分布
$$\phi (\mathbf{r})=\frac{1}{4\pi {\epsilon }_{0}r}\underset{V}{\iiint }\rho (\mathbf{r})\mathrm{d}V\text{\hspace{1em}(1.7)}$$
由电势的定义，相距为 $\mathrm{d}l$ 的两点的电势差为
$$\mathrm{d}\phi =-\mathbf{E}\cdot \mathrm{d}l$$
由于
$$\mathrm{d}\phi =\frac{\partial \phi }{\partial x}\mathrm{d}x+\frac{\partial \phi }{\partial y}\mathrm{d}y+\frac{\partial \phi }{\partial z}\mathrm{d}z=\nabla \phi \cdot \mathrm{d}l$$
所以电场强度和电势的关系为
$$\mathbf{E}=-\nabla \phi \text{\hspace{1em}(1.8)}$$
带入 $(1.4)$ 可得
$$\Delta \phi =-\frac{\rho }{{\epsilon }_0}\text{\hspace{1em}(1.9)}$$
这就是电势函数应当满足的静电场方程，$\mathbf{E}$ 是矢量，而 $\phi$ 是标量，求解标量方程相对简单些。

格林公式

高斯公式：(Gauss formula) 设空间闭区域 $V$ 由分片光滑的闭曲面 $\Sigma$ 所围成，函数$P(x,y,z),Q(x,y,z),R(x,y,z)$ 在 $V$ 上有连续的一阶偏导数，则有
$$\begin{array}{rl}\underset{V}{\iiint }(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z})\mathrm{d}V& =\underset{\Sigma }{\iint }P\mathrm{d}y\mathrm{d}z+Q\mathrm{d}x\mathrm{d}z+R\mathrm{d}x\mathrm{d}y\\ & =\underset{\Sigma }{\iint }(P\cos \alpha +Q\cos \beta +R\cos \gamma )\mathrm{d}S\end{array}$$
曲面 $\Sigma$ 的方向取外侧，$\cos \alpha ,\cos \beta ,\cos \gamma$ 为曲面 $\Sigma$ 在点 $(x,y,z)$ 处外法线的方向余弦⁴
高斯公式向量形式为
$$\underset{V}{\iiint }\nabla \cdot \mathbf{A}\mathrm{d}V=\underset{\Sigma }{\iint }\mathbf{A}\cdot \mathrm{d}\mathbf{S}=\underset{\Sigma }{\iint }\mathbf{A}\cdot \mathbf{n}\mathrm{d}S$$
其中 $\mathbf{A}=(P,Q,R)$，$\mathbf{n}=(\cos \alpha ,\cos \beta ,\cos \gamma )$ 为曲面 $\Sigma$ 单位法向量，$\mathrm{d}\mathbf{S}=\mathbf{n}\mathrm{d}S=\mathrm{d}y\mathrm{d}z\mathbf{i}+\mathrm{d}x\mathrm{d}z\mathbf{j}+\mathrm{d}x\mathrm{d}y\mathbf{k}$ 为单位元。

格林公式 ：设函数 $u(x,y,z),v(x,y,z)$ 在空间闭区域 $V$ 及边界 $\Sigma$ 上有一阶连续偏导数，在边界 $\Sigma$ 上有二阶连续偏导数 。利用高斯公式可得到
$$\underset{V}{\iiint }\nabla \cdot (u\nabla v)\mathrm{d}V=\underset{\Sigma }{\iint }u\nabla v\cdot \mathrm{d}\mathbf{S}$$
于是我们得到第一格林公式
$$\underset{\Sigma }{\iint }u\nabla v\cdot \mathrm{d}\mathbf{S}=\underset{V}{\iiint }u\Delta v\mathrm{d}V+\underset{V}{\iiint }\nabla u\cdot \nabla v\mathrm{d}V$$
同理我们可以得到
$$\underset{\Sigma }{\iint }v\nabla u\cdot \mathrm{d}\mathbf{S}=\underset{V}{\iiint }v\Delta u\mathrm{d}V+\underset{V}{\iiint }\nabla v\cdot \nabla u\mathrm{d}V$$
两式相减可得到
$$\underset{V}{\iiint }(u\Delta v-v\Delta u)\mathrm{d}V=\underset{\Sigma }{\iint }(u\nabla v-v\nabla u)\cdot \mathrm{d}\mathbf{S}=\underset{\Sigma }{\iint }(u\frac{\partial v}{\partial n}-v\frac{\partial u}{\partial n})\mathrm{d}S$$
其中 $\frac{\partial }{\partial n}$ 表示沿边界 $\Sigma$ 外法线的方向导数。上式称为第二格林公式，简称格林公式。

调和函数的边界性质：设函数 $u(\mathbf{r})$ 是区域 $V$ 内的调和函数，则有
$$\underset{\partial V}{\iint }\frac{\partial v}{\partial n}\mathrm{d}S=0$$
证明：在第二格林公式中取 $u(\mathbf{r})$ 为调和函数，即满足 $\Delta u=0$ ，取 $v=1$ ，即得上式。

调和函数的均值定理：设区域 $V$ 是以 ${\mathbf{r}}_0$ 为球心 $a$ 为半径的球，函数 $u(\mathbf{r})$ 是 $V$ 内的调和函数，则
$$u({\mathbf{r}}_0)=\frac{1}{4\pi a^2}\underset{\partial V}{\iint }u(\mathbf{r})\mathrm{d}S$$

证明：已知调和函数 $u(\mathbf{r})$ 满足
$$\Delta u(\mathbf{r})=0$$
取函数 $v(\mathbf{r})=\frac{1}{|\mathbf{r}-{\mathbf{r}}_0|}=\frac{1}{r}$ ，如图，利用多维 $\delta$ 函数的性质得到
$$\Delta v(\mathbf{r})=-4\pi \delta (\mathbf{r}-{\mathbf{r}}_0)$$
将上述两方程分别乘以 $v(\mathbf{r}),u(\mathbf{r})$ ，并相减，在体积 $V$ 内积分
$$\underset{V}{\iiint }(u\Delta v-v\Delta u)\mathrm{d}V=-4\pi \underset{V}{\iiint }u(\mathbf{r})\delta (\mathbf{r}-{\mathbf{r}}_0)=-4\pi u({\mathbf{r}}_0)$$
利用第二格林公式得到
$$u({\mathbf{r}}_0)=\frac{1}{4\pi }\underset{\partial V}{\iint }[v(\mathbf{r})\frac{\partial u(\mathbf{r})}{\partial n}-u(\mathbf{r})\frac{\partial v(\mathbf{r})}{\partial n}]\mathrm{d}S$$
因为在球面 $\partial V$ 上，外法线 $\mathbf{n}$ 的方向与 $r$ 所在半径的方向一致，所以球面上
$$\frac{\partial v(\mathbf{r})}{\partial n}=\frac{\partial }{\partial n}(\frac{1}{r})=\frac{\partial }{\partial r}(\frac{1}{r})=-\frac{1}{r^2}$$
又因为调和函数带入格林公式
$$\underset{V}{\iiint }\Delta u\mathrm{d}V=\underset{\partial V}{\iint }\frac{\partial u}{\partial n}\mathrm{d}S=0$$
于是最终得到
$$u({\mathbf{r}}_0)=\frac{1}{4\pi a^2}\underset{\partial V}{\iint }u(\mathbf{r})\mathrm{d}S$$
调和函数的极值原理：设函数 $u(\mathbf{r})$ 是区域 $V$ 内的调和函数，则 $u(\mathbf{r})$ 必在 $V$ 的边界面上取得最大值最小值。
证明：结合均值定理，可用反证法证明。

$\delta$ 函数简介

$\delta$ 函数起源于集中分布物理量的描述。
对于连续分布的物理量 $Q$ ，通常有两种描述方式，一种是局部性的，给出密度分布函数
$$\rho (\mathbf{r})=\frac{\mathrm{d}Q}{\mathrm{d}\mathbf{r}}$$
另一种是整体性的
$$Q={\int }_V\rho (\mathbf{r})\mathrm{d}\mathbf{r}$$
对于集中分布的物理量同样有两种方式描述。

$\delta$ 函数：(点电荷的线密度) 设在直线 $L$ 上，仅在 $x=0$ 处有一单位电荷，可以看成单位电荷分布在 $[-\epsilon ,\epsilon ]$ 上当 $\epsilon \to 0$ 的极限情况，后者密度可表示为
$${\rho }_{\epsilon }(x)=\{\begin{array}{ll}\frac{1}{2\epsilon }& (|x|⩽\epsilon )\\ 0& (|x|>\epsilon )\end{array}$$
且对任意 $\epsilon >0$ 直线上的电荷总量为
$$Q={\int }_{-\infty }^{+\infty }{\rho }_{\epsilon }(x)\mathrm{d}x=1$$
令 $\epsilon \to 0$ 可由 ${\rho }_{\epsilon }(x)$ 的极限推得单位点电荷的分布
$$\rho (x)=\{\begin{array}{ll}\infty & (x=0)\\ 0& (x\ne 0)\end{array}$$
且保持直线上的电荷总量为 1。
对于集中于 $x=0$ 点的单位物理量引起的密度函数叫做 $\delta$ 函数，$\delta (x)$ 满足条件
$$\delta (x)=\{\begin{array}{ll}\infty & (x=0)\\ 0& (x\ne 0)\end{array}$$
和
$${\int }_{-\infty }^{+\infty }\delta (x)\mathrm{d}x=1$$
注意：
(1) $\delta$ 函数并不是经典意义下的函数，因此通常称其为广义函数(或者奇异函数)。
(2) 它不能用常规意义下的值的对应关系来理解和使用，而总是通过它的定义和性质来使用它。
(3) $\delta$ 函数还有其他多种定义方式。

$\delta$ 函数的平移：对于集中于 $x=x_0$ 点的单位物理量引起的密度函数， $\delta$ 函数平移满足
$$\delta (x-x_0)=\{\begin{array}{ll}\infty & (x=x_0)\\ 0& (x\ne x_0)\end{array}$$
和
$${\int }_{-\infty }^{+\infty }\delta (x-x_0)\mathrm{d}x=1$$
筛选性质
$${\int }_a^b\delta (x)f(x)dx=\{\begin{array}{ll}f(0)& 0\in [a,b]\\ 0& 0\notin [a,b]\end{array}$$
特别的
$${\int }_{-\infty }^{+\infty }\delta (t)f(t)dt=f(0)$$
也可以把上述性质作为 $\delta$ 函数的另一种定义，此时我们对 $\delta$ 函数有了全新的认识，它实际上是一种映射，把元素 $f(x)$ 映射成了一个数 $f(0)$ 。

性质和运算

(1) $\delta (x)$ 和常数 $c$ 的乘积 $c\delta (x)$
$${\int }_{-\infty }^{+\infty }[c\delta (x)]f(x)dx={\int }_{-\infty }^{+\infty }\delta (x)[cf(x)]dx=cf(0)$$
(2) 筛选性质
$${\int }_{-\infty }^{+\infty }\delta (x-x_0)f(x)dx=f(x_0)$$
(3) 对称性
$$\delta (x-x_0)=\delta (x_0-x)$$
特别的
$$\delta (x)=\delta (-x)$$
(4) 与连续分布函数 $f(x)$ 的乘积
$${\int }_{-\infty }^{+\infty }f(x)\delta (x-x_0)={\int }_{-\infty }^{+\infty }f(x_0)\delta (x-x_0)$$
即
$$f(x)\delta (x-x_0)=f(x_0)\delta (x-x_0)$$
(5) $\delta$ 函数的导数 ${\delta }^{\prime }(x)$，对于在 $x=0$ 点连续并有连续导数的任意函数 $f(x)$ ，应用分部积分
$${\int }_{-\infty }^{+\infty }{\delta }^{\prime }(x)f(x)dx=\delta (x)f(x){|}_{-\infty }^{+\infty }-{\int }_{-\infty }^{+\infty }\delta (x)f^{\prime }(x)dx=-f^{\prime }(0)$$
对于 $\delta$ 函数的高阶导数 ${\delta }^{(n)}(x)$ ，对于在 $t=0$ 点连续并有连续导数的任意函数 $f(x)$ ，有
$${\int }_{-\infty }^{+\infty }{\delta }^{(n)}(x)f(x)dx=(-1{)}^n{f}^{(n)}(0)$$
(6) $\delta$ 函数是单位阶跃函数的导数
$$\frac{\mathrm{d}H(x)}{\mathrm{d}x}=\delta (x)$$
$\delta$ 函数的原函数为单位阶跃函数
$$H(x)=\{\begin{array}{ll}0& (x<0)\\ 1& (x>0)\end{array}={\int }_{-\infty }^x\delta (s)ds$$

(7) $\delta$ 函数的卷积
$$\begin{gathered} \delta (x)\ast f(x)=f(x) \\ {\delta }^{\prime }(x)\ast f(x)=\delta (x)\ast f^{\prime }(x)=f^{\prime }(x) \\ L[f\ast g]=L[f]\ast g=f\ast L[g] \end{gathered}$$

(8) 连续分布的质量、电荷或持续作用的力也可以用 $\delta$ 函数表示。现在以从 $t=a$ 持续作用到 $t=b$ 的作用力 $f(t)$ 为例说明。将时间段 $[a,b]$ 分成许多小段，在某个 $\tau$ 到 $\tau +d\tau$ 的短时间上，力 $f(t)$ 的冲量为 $f(\tau )d\tau$ ，既然 $d\tau$ 很短，不妨将这段时间的力看成是瞬时力，记作 $f(\tau )\delta (t-\tau )d\tau$ ，这样许许多多瞬时力的总和就是持续力 $f(t)$ ，即
$$f(t)=\sum f(\tau )\delta (t-\tau )d\tau ={\int }_a^{b}f(\tau )\delta (t-\tau )d\tau$$

δ函数的Fourier 变换
(1) 根据 $\delta$ 函数筛选性质可得
$$\begin{gathered} F(\omega )=\mathcal{F}[\delta (t)]={\int }_{-\infty }^{+\infty }\delta (t)e^{-i\omega t}\text{d}t=e^{-i\omega t}{|}_{t=0}=1 \\ \delta (t)={\mathcal{F}}^{-1}[1]=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{e}^{i\omega t}\text{d}\omega \end{gathered}$$
或写为
$$\delta (t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\cos \omega t\text{d}\omega =\frac{1}{\pi }{\int }_0^{+\infty }\cos \omega t\text{d}\omega$$
由此可见，δ函数包含所有频率成份，且它们具有相等的幅度，称此为均匀频谱或白色频谱。
我们可以得到 ：
$$\begin{array}{rl}& \delta (t)\leftrightarrow 1\\ & \delta (t-t_0)\leftrightarrow e^{-i\omega t_0}\\ & 1\leftrightarrow 2\pi \delta (\omega )\\ & e^{-i{\omega }_{0}t}\leftrightarrow 2\pi \delta (\omega -{\omega }_0)\\ & t\leftrightarrow 2\pi i{\delta }^{\prime }(\omega )\\ & \cos ({\omega }_{0}t)=\pi [\delta (\omega +{\omega }_0)+\delta (\omega -{\omega }_0)]\\ & \sin ({\omega }_{0}t)=i\pi [\delta (\omega +{\omega }_0)-\delta (\omega -{\omega }_0)]\end{array}$$
(2) 有许多重要的函数不满足Fourier 积分定理条件（绝对可积），例如常数、符号函数、单位阶跃函数、正弦函数和余弦函数等，但它们的广义Fourier 变换⁵也是存在的，利用δ函数及其Fourier 变换可以求出它们的Fourier 变换。

δ函数的Fourier 展开：当 $x,x_0\in (-\pi ,\pi )$ 时
$$\delta (x-x_0)=\frac{a_0}{2}+\sum _{n=1}^{\infty }(a_n\cos nx+b_n\sin nx)$$
其中傅里叶系数
$$\{\begin{array}{l}{a}_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }\delta (x-x_0)\cos nxdx=\frac{1}{\pi }\cos n{x}_0\\ b_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }\delta (x-x_0)\sin nxdx=\frac{1}{\pi }\sin n{x}_0\end{array}$$
多维 $\delta$ 函数：例如位于三维空间的坐标原点质量为 $m$ 的质点，其密度函数可表示为 $m\delta (\mathbf{r})$。 在三维空间中的 $\delta$ 函数定义如下：
$$\begin{gathered} \delta (\mathbf{r})=\{\begin{array}{ll}0& (\mathbf{r}\ne 0)\\ \infty & (\mathbf{r}=0)\end{array} \\ \underset{-\infty }{\overset{+\infty }{\iiint }}\delta (\mathbf{r})\mathrm{d}\mathbf{r}=1 \end{gathered}$$

(1) 三维 $\delta$ 函数可表示为三个一维 $\delta$ 函数乘积表示，在直角坐标系中
$$\delta (\mathbf{r})=\delta (x)\delta (y)\delta (z)$$
三维空间点 ${\mathbf{r}}_0=(x_0,y_0,z_0)$ 处密度分布函数就是
$$\delta (\mathbf{r}-{\mathbf{r}}_0)=\delta (x-x_0)\delta (y-y_0)\delta (z-z_0)$$
(2) 变量代换：当
$$\{\begin{array}{l}x=x(\xi ,\eta ,\zeta )\\ y=y(\xi ,\eta ,\zeta )\\ z=z(\xi ,\eta ,\zeta )\end{array}$$
时有
$$\delta (x-x_0,y-y_0,z-z_0)=\frac{1}{|J|}\delta (\xi -{\xi }_0,\eta -{\eta }_0,\zeta -{\zeta }_0)$$
其中 $|J|\ne 0$ 是 Jacobi 行列式的绝对值，$(x_0,y_0,z_0)$ 和 $({\xi }_0,{\eta }_0,{\zeta }_0)$ 相对应。
直角坐标系换算到柱坐标系 $\mathbf{r}=(r,\theta ,z)$
$$\delta (\mathbf{r}-{\mathbf{r}}_0)=\frac{1}{r_0}\delta (r-r_0)\delta (\theta -{\theta }_0)\delta (z-z_0)$$
直角坐标系换算到球坐标系 $\mathbf{r}=(r,\theta ,\varphi )$
$$\delta (\mathbf{r}-{\mathbf{r}}_0)=\frac{1}{r_0^2\sin {\theta }_0}\delta (r-r_0)\delta (\theta -{\theta }_0)\delta (\varphi -{\varphi }_0)$$
(3) 筛选性质
$$\begin{gathered} \underset{-\infty }{\overset{+\infty }{\iiint }}f(\mathbf{r})\delta (\mathbf{r}-{\mathbf{r}}_0)\mathrm{d}\mathbf{r}=f({\mathbf{r}}_0) \\ \underset{-\infty }{\overset{+\infty }{\iiint }}f(\mathbf{r})[\nabla \delta (\mathbf{r}-{\mathbf{r}}_0)]\mathrm{d}\mathbf{r}=-\nabla f(\mathbf{r}){|}_{\mathbf{r}={\mathbf{r}}_0} \end{gathered}$$
位矢的微分：
$$\Delta \frac{1}{r}=-4\pi \delta (\mathbf{r})$$
其中 $r=\sqrt{x^2+y^2+z^2}$

(4) 混合偏导：
$$\frac{{\partial }^{3}H(x,y,z)}{\partial x\partial y\partial z}=\delta (x,y,z)$$
其中 $H(x,y,z)=H(x)H(y)H(z)$ 为单位阶跃函数

(5) 多重傅里叶变换
$$\begin{gathered} \delta (x,y,z)\leftrightarrow 1 \\ 1\leftrightarrow (2\pi {)}^3\delta (\lambda ,\mu ,\nu ) \\ x\leftrightarrow (2\pi {)}^{3}i\frac{\partial \delta (\lambda ,\mu ,\nu )}{\partial \lambda } \\ {x}^2+y^2+z^2\leftrightarrow (2\pi {)}^3\delta (\lambda ,\mu ,\nu ) \\ {e}^{iax}\leftrightarrow (2\pi {)}^3\delta (\lambda -a,\mu ,\nu ) \end{gathered}$$
(6) 多重卷积定义
$$f\ast g=\underset{-\infty }{\overset{+\infty }{\iiint }}f(\mathbf{r})g(\mathbf{r}-{\mathbf{r}}_0)\mathrm{d}{\mathbf{r}}_0$$
性质如下
$$\begin{gathered} \delta \ast f=f \\ \frac{\partial \delta }{\partial x}\ast f=\delta \ast \frac{\partial f}{\partial x}=\frac{\partial f}{\partial x} \\ \frac{\partial }{\partial x}(f\ast g)=\frac{\partial f}{\partial x}\ast g=f\ast \frac{\partial g}{\partial x} \\ L[f\ast g]=L[f]\ast g=f\ast L[g] \\ \mathcal{F}(f\ast g)=\mathcal{F}(f)\cdot \mathcal{F}(g) \end{gathered}$$

---

1. 若函数 $f(x)$ 在区间 $D$ 上满足：
   (1) 连续或只有有限个第一类间断点；
   (2) 只有有限个极值点
   则称函数 $f(x)$ 在区间 $D$上满足狄利克雷(Dirichlet)条件 ↩︎ ↩︎

2. 这里要用到拉普拉斯变换得到的两个公式
   $$\begin{gathered} {\int }_0^{+\infty }\frac{\sin kx}{k}dk=\frac{\pi }{2} \\ {\int }_0^{+\infty }\frac{\sin kx}{k}{e}^{-k^2{a}^{2}t}dk=\sqrt{\pi }{\int }_0^{x/2a\sqrt{t}}{e}^{-{\xi }^2}d\xi \end{gathered}$$ ↩︎

3. 立体角：常用字母 $\Omega$ 表示。以观测点为球心，构造一个单位球面；任意物体投影到该单位球面上的投影面积，即为该物体相对于该观测点的立体角。 因此，立体角是单位球面上的一块面积，这和“平面角是单位圆上的一段弧长”类似。
   立体角的国际制单位是球面度 (steradian , sr)
   在球坐标系中，任意球面的极小面积为：
   $$\mathrm{d}A=(r\sin \theta \mathrm{d}\varphi )\cdot (r\mathrm{d}\theta )=r^2\sin \theta \mathrm{d}\theta \mathrm{d}\varphi$$
   因此，极小立体角（单位球面上的极小面积）为：
   $$\mathrm{d}\Omega =\frac{\mathrm{d}A}{r^2}=\sin \theta \mathrm{d}\theta \mathrm{d}\varphi$$
   所以，立体角是投影面积与球半径平方值的比，这和“平面角是圆的弧长与半径的比”类似。 对极小立体角做曲面积分即可得立体角
   $$\Omega =\underset{S}{\iint }\mathrm{d}\Omega =\underset{S}{\iint }\sin \theta \mathrm{d}\theta \mathrm{d}\varphi$$
   一个完整的球面对于球内任意一点的立体角为 $4\pi$， 这个定理对所有封闭曲面皆成立。
   $$\underset{S}{\iint }\mathrm{d}\Omega =\underset{S}{\iint }\sin \theta \mathrm{d}\theta \mathrm{d}\varphi ={\int }_0^{\pi }\sin \theta \mathrm{d}\theta {\int }_0^{2\pi }\mathrm{d}\varphi =4\pi$$ ↩︎

4. 方向向量与坐标轴的夹角 $\alpha ,\beta ,\gamma$ 称为方向角，$\cos \alpha ,\cos \beta ,\cos \gamma$ 称为方向向量的方向余弦。 ↩︎

5. 在δ函数的Fourier变换中，其广义积分是根据δ函数的性质直接给出的，而不是按通常的积分方式得到的，称这种方式的Fourier 变换为广义Fourier 变换。 ↩︎