调和函数.

若$\Delta u=0$，证明$u(\mathbf{x})=\underset{\partial B_r}{\int }-u(\mathbf{y})\text{d}S$。

$\Delta$为拉普拉斯算子，$\underset{\partial B_r}{\int }-$为以$\mathbf{x}$为球心、$r$为半径的球面上的平均积分。

---

由于$\Delta u=\text{div}(\nabla u)$，由高斯公式得

$$\underset{B_r}{\int }\Delta u(\mathbf{y})\text{d}V=\underset{B_r}{\int }\text{div}(\nabla u(\mathbf{y}))\text{d}V=\underset{\partial B_r}{\int }\nabla u(\mathbf{y})\cdot {\mathbf{e}}_n\text{d}S$$

由$\Delta u=0$，可得上式左端

$$\underset{B_r}{\int }\Delta u(\mathbf{y})\text{d}V=0$$

对于上式右端，由于$\mathbf{y}$位于以$\mathbf{x}$为球心、$r$为半径的球面上，可设$\mathbf{\omega }=\mathbf{y}-\mathbf{x}$，则$|\mathbf{\omega }|=r$，且$\mathbf{\omega }$沿径向方向。因此

$$\begin{array}{rl}\underset{\partial B_r}{\int }\nabla u(\mathbf{y})\cdot {\mathbf{e}}_n\text{d}S& =\underset{|\mathbf{\omega }|=r}{\int }\nabla u(\mathbf{x}+\mathbf{\omega })\cdot {\mathbf{e}}_n\text{d}{S}^{\prime }=\underset{|\mathbf{\omega }|=r}{\int }\nabla u(\mathbf{x}+\mathbf{\omega })\cdot \frac{\mathbf{\omega }}{|\mathbf{\omega }|}\text{d}{S}^{\prime }\\ & =r^2\underset{|\mathbf{\omega }|=1}{\int }\nabla u(\mathbf{x}+r\mathbf{\omega })\cdot \mathbf{\omega }\text{d}{S}^{\prime }=r^2\underset{|\mathbf{\omega }|=1}{\int }\nabla u(\mathbf{x}+r\mathbf{\omega })\cdot \mathbf{\omega }\text{d}{S}^{\prime }\end{array}$$

记$u$对第$i$维变量的偏导数为$u_i$，$\mathbf{\omega }=[{\omega }_1,{\omega }_2,\cdots ,{\omega }_n{]}^{\text{T}}$，则上式中

$$\nabla u(\mathbf{x}+r\mathbf{\omega })\cdot \mathbf{\omega }=\sum _{i=1}^n{u}_i(\mathbf{x}+r\mathbf{\omega })\cdot {\omega }_i=\frac{\partial }{\partial r}u(\mathbf{x}+r\mathbf{\omega })$$

因此之前的式子右端进一步化为

$$\underset{\partial B_r}{\int }\nabla u(\mathbf{y})\cdot {\mathbf{e}}_n\text{d}S=r^2\underset{|\mathbf{\omega }|=1}{\int }\frac{\partial }{\partial r}u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }=r^2\frac{\partial }{\partial r}\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }$$

由于左端等于$0$，所以$\frac{\partial }{\partial r}\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }=0$，即$\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }$关于$r$是常值函数。令$r\to 0$，得

$$\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }=\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x})\text{d}{S}^{\prime }=u(\mathbf{x})\underset{|\mathbf{\omega }|=1}{\int }\text{d}{S}^{\prime }=u(\mathbf{x})\cdot \alpha$$

其中$\alpha$为半径为$1$的球面面积。因此

$$\begin{array}{rl}u(\mathbf{x})& =\frac{1}{\alpha }\underset{|\mathbf{\omega }|=1}{\int }u(\mathbf{x}+r\mathbf{\omega })\text{d}{S}^{\prime }=\frac{1}{r^2\alpha }\underset{|\mathbf{\omega }|=r}{\int }u(\mathbf{x}+\mathbf{\omega })\text{d}{S}^{\prime }\\ & =\frac{1}{r^2\alpha }\underset{\partial B_r}{\int }u(\mathbf{y})\text{d}S=\frac{1}{|\partial B_r|}\underset{\partial B_r}{\int }u(\mathbf{y})\text{d}S=\underset{\partial B_r}{\int }-u(\mathbf{y})\text{d}S\end{array}$$