浅谈Green公式和外微分形式

数学分析大作业6，时隔一年半才想起来没发。
衷心感谢北航的杨小远老师。
希望学弟学妹们能够继续认真的完成杨老师的大作业。

经过一年半的继续学习和理解才发现
Green公式的作用，到用Green函数解偏微分方程时候才能真正体现出来。
Green函数法解偏微分方程，具体就有很多应用了。
比如，简单的静电场的边值问题，直接就是Green函数模型提出的背景
比如，量子力学计算散射问题，声子的引入就类似于Green函数提出的思想。

关于外微分形式，是出于如何理解stokes定理和Green定理诡异的相似性而衍生出来的讨论。这个讨论在龚老师的《简明微积分》里有比较全面的讨论，文中的部分也主要来自于《简明微积分》
强烈安利《简明微积分》，这本书看完以后，对于微积分的理解会更上升一个层次。

关于雅可比行列式的讨论，其实可以归结到度量张量底下。这一部分内容如果用坐标系的选择与换元，和nabla算子在不同坐标系下的形式联系起来会更有趣。继续静待后文。

1 二维平面上的Green公式

1.1 单连通区域上的Green公式

设$S$是由分段光滑的封闭曲线$C$围成的闭区域,满足(1)向量场$F=(P(x,y),Q(x,y))$在$S$上连续;(2)$P(x,y),Q(x,y)$在$S$上有连续的一阶偏导数,则成立
$$\begin{gathered} {\iint }_S\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dxdy={\int }_{C}Pdx+Qdy \\ {\iint }_S\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dxdy={\int }_C[P\cos (t,x)+Q\cos (t,y)]ds \end{gathered}$$
取$C$正向.$(\cos (t,x),\cos (t,y))$为$C$正向一致的切线方向.

设$\vec{n}$为$C$的外法线单位向量,则
$${\iint }_S\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}\right)dxdy={\int }_C[P\cos (\vec{n},x)+Q\cos (\vec{n},y)]ds$$

注意 用$(\vec{a},\vec{b})$表示向量$\vec{b}$转到向量$\vec{a}$走过的角度,则有
$$\begin{gathered} (t,x)=(\vec{n},x)+\frac{\pi }{2},(t,y)=(\vec{n},y)+\frac{\pi }{2},(t,x)=(t,y)+\frac{\pi }{2} \\ \cos (t,x)=-cos(\vec{n},y),\cos (t,y)=\cos (\vec{n},x) \end{gathered}$$

Green第一公式

设$C$为封闭光滑曲线,$S$为$C$围成的区域,$u$有连续的二阶偏导数,$\omega$有连续的偏导数,$\vec{n}$为$C$的外法线单位向量,取$C$的正向,则
$$\begin{gathered} {\iint }_S\omega \frac{\partial u}{\partial x}dxdy+{\iint }_{S}u\frac{\partial \omega }{\partial x}dxdy={\int }_{C}uwdy \\ {\iint }_S\omega \Delta udxdy+{\iint }_S\left(\frac{\partial u}{\partial x}\frac{\partial \omega }{\partial x}+\frac{\partial u}{\partial y}\frac{\partial \omega }{\partial y}\right)dxdy={\int }_C\omega \frac{\partial u}{\partial \vec{n}}ds \end{gathered}$$
特别地,当$\omega =u$时,有
$${\iint }_S\left[{\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2\right]dxdy+{\iint }_{S}u\Delta udxdy={\oint }_{C}u\frac{\partial u}{\partial \vec{n}}ds$$
证明 注意到
$$\omega \frac{\partial u}{\partial \vec{n}}=\omega \frac{\partial u}{\partial x}\cos (\vec{n},x)+\omega \frac{\partial u}{\partial y}\cos (\vec{n},y)$$
即可得到结论.

Green第二公式

设$C$为封闭光滑曲线,$S$为$C$围成的区域,$u,v$有连续的二阶偏导数,$\vec{n}$为$C$的外法线单位向量,取$C$的正向,则
$${\iint }_S\left|\begin{array}{cc}\Delta u& \Delta v\\ u& v\end{array}\right|dxdy={\oint }_c\left|\begin{array}{cc}\frac{\partial u}{\partial \vec{n}}& \frac{\partial v}{\partial \vec{n}}\\ u& v\end{array}\right|ds$$
证明 在Green第一公式中分别取$u,v$并交换位置,得
$$\begin{gathered} {\iint }_{S}v\Delta udxdy+{\iint }_S\left(\frac{\partial u}{\partial x}\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\frac{\partial v}{\partial y}\right)dxdy={\int }_{C}v\frac{\partial u}{\partial \vec{n}}ds \\ {\iint }_{S}u\Delta vdxdy+{\iint }_S\left(\frac{\partial v}{\partial x}\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\frac{\partial u}{\partial y}\right)dxdy={\int }_{C}u\frac{\partial v}{\partial \vec{n}}ds \end{gathered}$$
两等式作差,即得结论.

Green第三公式

若$u$为有界闭区域$S$内的调和函数,$C$为$S$的边界,$$\frac{\partial u}{\partial \vec{n}}$$为$u$沿着$C$的单位外法线的方向导数,$$\vec{r}=(a-x,b-y)$$, $$r=\sqrt{(a-x{)}^2+(b-y{)}^2}$$为$(x,y)$与$C$上动点$(a,b)$之间的距离,则
$$u(x,y)=\frac{1}{2\pi }{\oint }_C\left(u\frac{\cos (\vec{r},\vec{n})}{r}-\ln r\frac{\partial u}{\partial \vec{n}}\right)ds=\frac{1}{2\pi }{\oint }_C\left(u\frac{\partial \ln r}{\partial \vec{n}}-\ln r\frac{\partial u}{\partial \vec{n}}\right)ds$$
证明 在结论等式中有一点需要额外说明
$$\frac{\cos (\vec{r},\vec{n})}{r}=\frac{1}{r}\frac{a-x}{r}\cos (\vec{n},x)+\frac{1}{r}\frac{y-b}{r}\cos (\vec{n},y)=\frac{\partial \ln r}{\partial \vec{n}}$$

$$\frac{{\partial }^2\ln r}{\partial a^2}+\frac{{\partial }^2\ln r}{\partial b^2}=\frac{r^2-2(a-x{)}^2}{r^4}+\frac{r^2-2(b-y{)}^2}{r^4}=\frac{2r^2-2r^2}{r^4}=0,r̸=0$$

即$\ln r,r̸=0$为调和函数.设 S 0 = { ( a , b ) ∣ ( a − x ) 2 + ( b − y ) 2 ≤ ρ 2 } ⊂ S S_0=\{(a,b)|(a-x)^2+(b-y)^2\le \rho^2\}\subset S S0​={(a,b)∣(a−x)2+(b−y)2≤ρ2}⊂S,$C_0$为$S_0$的边界,则
$${\oint }_{C-C_0}\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds={\iint }_{S-S_0}\left|\begin{array}{cc}\Delta \ln r& \Delta u\\ \ln r& u\end{array}\right|dadb=0$$

$$\begin{array}{rl}{\oint }_C\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds& ={\oint }_{C_0}\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds\\ & =\frac{1}{\rho }{\oint }_{C_0}uds-\ln \rho {\oint }_{C_0}\frac{\partial u}{\partial \vec{n}}ds\\ & =2\pi u(x+{\theta }_1\rho ,y+{\theta }_2\rho ),0<{\theta }_1<1,0<{\theta }_2<1\end{array}$$

对上式取极限$\rho \to 0^{+}$,注意到$u(x,y)$的连续性,则证毕.

调和函数的中值定理

设 Δ u = 0 , C ∗ = { ( x , y ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 = R 2 } , S ∗ = { ( x , y ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 ≤ R 2 } \Delta u=0,C^*=\{(x,y)|(x-x_0)^2+(y-y_0)^2=R^2\},S^*=\{(x,y)|(x-x_0)^2+(y-y_0)^2\le R^2\} Δu=0,C∗={(x,y)∣(x−x0​)2+(y−y0​)2=R2},S∗={(x,y)∣(x−x0​)2+(y−y0​)2≤R2},则
$$u(x_0,y_0)=\frac{1}{2\pi R}{\oint }_{C^{\ast }}u(x,y)ds$$
证明 用Green第三公式,取$C=C^{\ast }$,即
$$u(x_0,y_0)=\frac{1}{2\pi }{\oint }_{C^{\ast }}\left(u\frac{\partial \ln r}{\partial \vec{n}}-\ln r\frac{\partial u}{\partial \vec{n}}\right)ds$$
在$C^{\ast }$上满足
$${\frac{\partial \ln r}{\partial \vec{n}}|}_{r=R}={\frac{\partial \ln r}{\partial r}|}_{r=R}={\frac{1}{r}|}_{r=R}=\frac{!}{R}$$
且
$${\oint }_{C^{\ast }}\frac{\partial u}{\partial \vec{n}}ds=0$$
该结论在之前证明Green第三公式时也使用过,则可得到
$$u(x_0,y_0)=\frac{1}{2\pi }{\oint }_{C^{\ast }}\left(\frac{u}{R}-\ln R\frac{\partial u}{\partial \vec{n}}\right)ds=\frac{1}{2\pi R}{\oint }_{C^{\ast }}u(x,y)ds$$

极大值原理

连通的有界闭区域上的非常数调和函数$u(x,y)$在此区域内部的点不能达到其最大值或最小值.

证明 设该有界闭区域为$\overline{S}$,由开区域$S$和边界$C=\partial S$构成.用反证法证明,假设在$S$上能取到$u(x,y)$的最大值或最小值,不妨设某一个最大值点为$(x_0,y_0)$,设圆域 S R = { ( x , y ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 ≤ R 2 } S_R=\{(x,y)|(x-x_0)^2+(y-y_0)^2\le R^2\} SR​={(x,y)∣(x−x0​)2+(y−y0​)2≤R2},圆周 C R = { ( x , y ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 = R 2 } C_R=\{(x,y)|(x-x_0)^2+(y-y_0)^2=R^2\} CR​={(x,y)∣(x−x0​)2+(y−y0​)2=R2}.

1. 若$S_R\subset S$,对$\forall r\in (0,R]$,由调和函数的中值定理有
   $$u(x_0,y_0)=\frac{1}{2\pi r}{\oint }_{C_r}u(x,y)ds$$
   故
   $$\frac{1}{2\pi r}{\oint }_{C_r}[u(x_0,y_0)-u(x,y)]ds=0$$
   因$u(x_0,y_0)$为最大值,则有$u(x_0,y_0)-u(x,y)\ge 0$,因而在$C_r$上满足
   $$u(x_0,y_0)-u(x,y)\equiv 0$$
   根据$r$的任意性,可知对于任何$(x,y)\in S_R$都有$u(x_0,y_0)-u(x,y)\equiv 0$成立,则$u(x,y)$在$S_R$上是一个常数.

2. 由于该有界区域是连通的,则不同于(x_0,y_0)的任意内点$(x,y)\in S$,存在一条完全含于$S$的曲线将$(x,y)$和$(x_0,y_0)$连接起来,并且可以用有限个完全含于$S$的圆$S_1^{\ast },S_2^{\ast },\cdots .S_k^{\ast }$完全覆盖这条曲线,并且满足
   $$(x,y)\in S_1^{\ast },(x_0,y_0)\in S_k^{\ast },S_1^{\ast }\cap S_2^{\ast }̸=\varnothing ,S_2^{\ast }\cap S_3^{\ast }̸=\varnothing ,\cdots ,S_{k-1}^{\ast }\cap S_k^{\ast }̸=\varnothing$$
   根据1知$S_k^{\ast }$上$u$为常数,则$S_{k-1}^{\ast },\cdots ,S_2^{\ast },S_1^{\ast }$上$u$均为常数,特别地,$u(x,y)$也为常数.根据$(x,y)$的任意性,则$u(x,y)$在$S$上是一个常数.

3. 根据$u(x,y)$在$\overline{S}$上的连续性,通过$S$的点趋向$C=\partial S$上的点取极限,则$u(x,y)$在$\overline{S}$上为常数.

综上所述,$u$为常数,这与条件"非常数调和函数"矛盾,因而假设有误.假设最小值点同理可证,则$u(x,y)$的最值点不能在区域内部取得.

1.2 复连通区域上的Green公式

对于复连通区域,边界的正向按照沿边界正向增加时研究区域始终在左手侧来定义,曲线的外法线方向始终由研究区域内部出发指向研究区域的外侧.例如,对于满足
S = { ( x , y ) ∣ a 2 ≤ x 2 + y 2 ≤ b 2 } , C 1 = { ( x , y ) ∣ x 2 + y 2 = a 2 } , C 2 = { ( x , y ) ∣ x 2 + y 2 = b 2 } , a &lt; b , ∂ S = C 1 ∪ C 2 S=\{(x,y)|a^2\le x^2+y^2\le b^2\},C_1=\{(x,y)|x^2+y^2=a^2\},C_2=\{(x,y)|x^2+y^2=b^2\},a&lt;b,\partial S=C_1\cup C_2 S={(x,y)∣a2≤x2+y2≤b2},C1​={(x,y)∣x2+y2=a2},C2​={(x,y)∣x2+y2=b2},a<b,∂S=C1​∪C2​
的区域$S$,边界$C_1$上的外法线单位向量为$$\vec{n}=(-\frac{x}{a},-\frac{y}{a})$$,而边界$C_2$上的外法线单位向量为$$\vec{n}=(\frac{x}{b},\frac{y}{b})$$.圆环内侧由圆环指向原点,圆环外侧背离原点方向.因而通过外法线单位向量来表述Green公式,可以得到叙述形式不变的结论:

$S$是由分段光滑的封闭曲线$C_1,C_2,\cdots ,C_k$围成的闭区域,满足(1)向量场$F=(P(x,y),Q(x,y))$在$S$上连续;(2)$P(x,y),Q(x,y)$在$S$上有连续的一阶偏导数,设$\vec{n}$为$C_1,C_2,\cdots ,C_k$的外法线单位向量,用"+"表示取正向,则
$${\iint }_S\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}\right)dxdy={\int }_{C_1+C_2+\cdots +C_k}[P\cos (\vec{n},x)+Q\cos (\vec{n},y)]ds$$
其余公式的证明没有较大的区别,就不给出了.

Green第一公式

设$C_1,C_2,\cdots ,C_k$为封闭光滑曲线,$S$为$C_1,C_2,\cdots ,C_k$围成的区域,$u$有连续的二阶偏导数,$\omega$有连续的偏导数,$\vec{n}$为$C_1,C_2,\cdots ,C_k$的外法线单位向量,取$C_1,C_2,\cdots ,C_k$的正向,则
$$\begin{gathered} {\iint }_S\omega \frac{\partial u}{\partial x}dxdy+{\iint }_{S}u\frac{\partial \omega }{\partial x}dxdy={\int }_{C_1+C_2+\cdots +C_k}uwdy \\ {\iint }_S\omega \Delta udxdy+{\iint }_S\left(\frac{\partial u}{\partial x}\frac{\partial \omega }{\partial x}+\frac{\partial u}{\partial y}\frac{\partial \omega }{\partial y}\right)dxdy={\int }_{C_1+C_2+\cdots +C_k}\omega \frac{\partial u}{\partial \vec{n}}ds \end{gathered}$$
特别地,当$\omega =u$时,有
$${\iint }_S\left[{\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2\right]dxdy+{\iint }_{S}u\Delta udxdy={\oint }_{C_1+C_2+\cdots +C_k}u\frac{\partial u}{\partial \vec{n}}ds$$

Green第二公式

设$C_1,C_2,\cdots ,C_k$为封闭光滑曲线,$S$为$C_1,C_2,\cdots ,C_k$围成的区域,$u,v$有连续的二阶偏导数,$\vec{n}$为$C_1,C_2,\cdots ,C_k$的外法线单位向量,取$C_1,C_2,\cdots ,C_k$的正向,则
$${\iint }_S\left|\begin{array}{cc}\Delta u& \Delta v\\ u& v\end{array}\right|dxdy={\oint }_{C_1+C_2+\cdots +C_k}\left|\begin{array}{cc}\frac{\partial u}{\partial \vec{n}}& \frac{\partial v}{\partial \vec{n}}\\ u& v\end{array}\right|ds$$

Green第三公式

若$u$为有界闭区域$S$内的调和函数,$C_1,C_2,\cdots ,C_k$为$S$的边界,$$\frac{\partial u}{\partial \vec{n}}$$为$u$沿着$C_1,C_2,\cdots ,C_k$的单位外法线的方向导数,$$\vec{r}=(a-x,b-y)$$, $$r=\sqrt{(a-x{)}^2+(b-y{)}^2}$$为$(x,y)$与$C_1,C_2,\cdots ,C_k$上动点$(a,b)$之间的距离,则
$$u(x,y)=\frac{1}{2\pi }{\oint }_{C_1+C_2+\cdots +C_k}\left(u\frac{\cos (\vec{r},\vec{n})}{r}-\ln r\frac{\partial u}{\partial \vec{n}}\right)ds=\frac{1}{2\pi }{\oint }_{C_1+C_2+\cdots +C_k}\left(u\frac{\partial \ln r}{\partial \vec{n}}-\ln r\frac{\partial u}{\partial \vec{n}}\right)ds$$
注意
$${\oint }_{C_1+C_2+\cdots +C_k-C_0}\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds={\iint }_{S-S_0}\left|\begin{array}{cc}\Delta \ln r& \Delta u\\ \ln r& u\end{array}\right|dadb=0$$

$${\oint }_{C_1+C_2+\cdots +C_k}\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds={\oint }_{C_0}\left|\begin{array}{cc}\frac{\partial \ln r}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \ln r& u\end{array}\right|ds$$

中值定理和极大值原理同样满足,这里不再重复.

2 三维空间上的Gauss公式

2.1 单连通区域上的Gauss公式

设$V$是由分片光滑的双侧封闭曲面$S$围成的闭区域,满足(1)向量场$F=(P(x,y,z),Q(x,y,z),R(x,y,z))$在$V$上连续;(2)$P(x,y,z),Q(x,y,z),R(x,y,z)$在$V$上有连续的一阶偏导数,则成立
$${\iiint }_V\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dxdydz={\iint }_{S}Pdydz+Qdzdx+Rdxdy$$
取$S$的外侧.

设$\vec{n}$为$S$的外法线单位向量,则
$${\iiint }_V\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dxdydz={\iint }_S[P\cos (\vec{n},x)+Q\cos (\vec{n},y)+R\cos (\vec{n},z)]dS$$

空间Green第一公式

设$S$为封闭光滑曲面,$V$为$S$围成的区域,$u$在$S$和$V$上有二阶连续的偏导数,$\omega$有连续的偏导数,$\vec{n}$为$S$的外法线单位向量,取$S$的外侧,则
$$\begin{gathered} {\iiint }_V\omega \frac{\partial u}{\partial x}dxdydz+{\iiint }_{V}u\frac{\partial \omega }{\partial x}dxdydz={\iint }_{S}uwdydz \\ {\iiint }_V\omega \Delta udxdydz+{\iiint }_V\left(\frac{\partial u}{\partial x}\frac{\partial \omega }{\partial x}+\frac{\partial u}{\partial y}\frac{\partial \omega }{\partial y}+\frac{\partial u}{\partial z}\frac{\partial \omega }{\partial z}\right)dxdydz={\iint }_S\omega \frac{\partial u}{\partial \vec{n}}dS \end{gathered}$$
特别地,当$\omega =u$时,有
$${\iiint }_V\left[{\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2+{\left(\frac{\partial u}{\partial z}\right)}^2\right]dxdydz+{\iiint }_{V}u\Delta udxdydz={\iint }_{S}u\frac{\partial u}{\partial \vec{n}}dS$$

证明 证明方法同上,即有
$$\omega \frac{\partial u}{\partial \vec{n}}=\omega \frac{\partial u}{\partial x}\cos (\vec{n},x)+\omega \frac{\partial u}{\partial y}\cos (\vec{n},y)+\omega \frac{\partial u}{\partial y}\cos (\vec{n},z)$$
即可证得结论.

空间Green第二公式

设$S$为封闭光滑曲面,$V$为$S$围成的区域,$u,v$在$S$和$V$上有二阶连续的偏导数,$\vec{n}$为$S$的外法线单位向量,取$S$的外侧,则
$${\iiint }_V\left|\begin{array}{cc}\Delta u& \Delta v\\ u& v\end{array}\right|dxdydz={\iint }_S\left|\begin{array}{cc}\frac{\partial u}{\partial \vec{n}}& \frac{\partial v}{\partial \vec{n}}\\ u& v\end{array}\right|dS$$

证明 在空间格林第一公式里面写出关于$u,v$的等式并交换位置,有
$$\begin{gathered} {\iiint }_{V}v\Delta udxdydz+{\iiint }_V\left(\frac{\partial u}{\partial x}\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\frac{\partial v}{\partial y}+\frac{\partial u}{\partial z}\frac{\partial v}{\partial z}\right)dxdydz={\iint }_{S}v\frac{\partial u}{\partial \vec{n}}ds \\ {\iiint }_{V}u\Delta vdxdydz+{\iiint }_V\left(\frac{\partial v}{\partial x}\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\frac{\partial u}{\partial y}+\frac{\partial v}{\partial z}\frac{\partial u}{\partial z}\right)dxdydz={\iint }_{S}u\frac{\partial v}{\partial \vec{n}}ds \end{gathered}$$
二式相减即得结论.

空间Green第三公式

若$u$为有界闭区域$V$内的调和函数,$S$为$V$的边界,$$\frac{\partial u}{\partial \vec{n}}$$为$u$沿着$S$的单位外法线的方向导数,$$\vec{r}=(a-x,b-y,c-z)$$,$$r=\sqrt{(a-x{)}^2+(b-y{)}^2+(c-z{)}^2}$$为$(x,y,z)$与$S$上动点$(a,b,c)$之间的距离,则
$$u(x,y,z)=\frac{1}{4\pi }{\iint }_S\left(u\frac{\cos (\vec{r},\vec{n})}{r^2}+\frac{1}{r}\frac{\partial u}{\partial \vec{n}}\right)dS=\frac{1}{4\pi }{\iint }_S\left(u\frac{\partial \frac{-1}{r}}{\partial \vec{n}}-\frac{-1}{r}\frac{\partial u}{\partial \vec{n}}\right)dS$$

证明 在结论等式中有一点需要额外说明
$$\frac{\cos (\vec{r},\vec{n})}{r^2}=\frac{1}{r^2}\frac{a-x}{r}\cos (\vec{n},x)+\frac{1}{r^2}\frac{y-b}{r}\cos (\vec{n},y)+\frac{1}{r^2}\frac{z-c}{r}\cos (\vec{n},y)=\frac{\partial \frac{-1}{r}}{\partial \vec{n}}$$

$$\frac{{\partial }^2\frac{-1}{r}}{\partial a^2}+\frac{{\partial }^2\frac{-1}{r}}{\partial b^2}+\frac{{\partial }^2\frac{-1}{r}}{\partial c^2}=\frac{r^3-3r(a-x{)}^2}{r^6}+\frac{r^3-3r(b-y{)}^2}{r^6}+\frac{r^3-3r(c-z{)}^2}{r^6}=\frac{3r^3-3r^3}{r^6}=0,r̸=0$$

即$$\frac{-1}{r},r̸=0$$为调和函数.设 V 0 = { ( a , b , c ) ∣ ( a − x ) 2 + ( b − y ) 2 + ( c − z ) 2 ≤ ρ 2 } ⊂ V V_0=\{(a,b,c)|(a-x)^2+(b-y)^2+(c-z)^2\le \rho^2\}\subset V V0​={(a,b,c)∣(a−x)2+(b−y)2+(c−z)2≤ρ2}⊂V,$S_0$为$V_0$的边界,则
$${\iint }_{S-S_0}\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS={\iiint }_{V-V_0}\left|\begin{array}{cc}\Delta \frac{-1}{r}& \Delta u\\ \frac{-1}{r}& u\end{array}\right|dadbdc=0$$

$$\begin{array}{rl}{\iint }_S\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS& ={\iint }_{S_0}\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS\\ & =\frac{1}{{\rho }^2}{\iint }_{S_0}udS-\frac{-1}{\rho }{\iint }_{S_0}\frac{\partial u}{\partial \vec{n}}dS\\ & =4\pi u(x+{\theta }_1\rho ,y+{\theta }_2\rho ,z+{\theta }_3\rho ),0<{\theta }_1<1,0<{\theta }_2<1,0<{\theta }_3<1\end{array}$$

对上式取极限$\rho \to 0^{+}$,则证毕.

调和函数的中值定理

设 Δ u = 0 , S ∗ = { ( x , y , z ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 = R 2 } \Delta u=0,S^*=\{(x,y,z)|(x-x_0)^2+(y-y_0)^2+(z-z_0)^2=R^2\} Δu=0,S∗={(x,y,z)∣(x−x0​)2+(y−y0​)2+(z−z0​)2=R2}, V ∗ = { ( x , y , z ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 ≤ R 2 } V^*=\{(x,y,z)|(x-x_0)^2+(y-y_0)^2+(z-z_0)^2\le R^2\} V∗={(x,y,z)∣(x−x0​)2+(y−y0​)2+(z−z0​)2≤R2},则
$$u(x_0,y_0,z_0)=\frac{1}{4\pi R^2}{\iint }_{S^{\ast }}u(x,y,z)dS$$
证明 用空间Green第三公式,取$S=S^{\ast }$,即
$$u(x_0,y_0,z_0)=\frac{1}{4\pi }{\iint }_{S^{\ast }}\left(u\frac{\partial \frac{-1}{r}}{\partial \vec{n}}-\frac{-1}{r}\frac{\partial u}{\partial \vec{n}}\right)dS$$
在$C^{\ast }$上满足
$${\frac{\partial \frac{-1}{r}}{\partial \vec{n}}|}_{r=R}={\frac{\partial \frac{-1}{r}}{\partial r}|}_{r=R}={\frac{1}{r^2}|}_{r=R}=\frac{1}{R^2}$$
且
$${\iint }_{S^{\ast }}\frac{\partial u}{\partial \vec{n}}dS=0$$
该结论在之前证明空间Green第三公式时也使用过,则可得到
$$u(x_0,y_0,z_0)=\frac{1}{4\pi }{\iint }_{S^{\ast }}\left(\frac{u}{R^2}-\frac{-1}{R}\frac{\partial u}{\partial \vec{n}}\right)dS=\frac{1}{4\pi R^2}{\iint }_{S^{\ast }}u(x,y,z)dS$$

极大值原理

连通的有界闭区域上的非常数调和函数$u(x,y,z)$在此区域内部的点不能达到其最大值或最小值.

证明 设该有界闭区域为$\overline{V}$,由开区域$V$和边界$S=\partial V$构成.用反证法证明,假设在$V$上能取到$u(x,y,z)$的最大值或最小值,不妨设某一个最大值点为$(x_0,y_0,z_0)$,设球 V R = { ( x , y , z ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 ≤ R 2 } V_R=\{(x,y,z)|(x-x_0)^2+(y-y_0)^2+(z-z_0)^2\le R^2\} VR​={(x,y,z)∣(x−x0​)2+(y−y0​)2+(z−z0​)2≤R2},球面 S R = { ( x , y , z ) ∣ ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 = R 2 } S_R=\{(x,y,z)|(x-x_0)^2+(y-y_0)^2+(z-z_0)^2=R^2\} SR​={(x,y,z)∣(x−x0​)2+(y−y0​)2+(z−z0​)2=R2}.

1. 若$V_R\subset V$,对$\forall r\in (0,R]$,由调和函数的中值定理有
   $$u(x_0,y_0,z_0)=\frac{1}{4\pi r^2}{\iint }_{S_r}u(x,y,z)dS$$
   故
   $$\frac{1}{4\pi r^2}{\iint }_{S_r}[u(x_0,y_0,z_0)-u(x,y,z)]dS=0$$
   因$u(x_0,y_0,z_0)$为最大值,则有$u(x_0,y_0,z_0)-u(x,y,z)\ge 0$,因而在$C_r$上满足
   $$u(x_0,y_0,z_0)-u(x,y,z)\equiv 0$$
   根据$r$的任意性,可知对于任何$(x,y,z)\in V_R$都有$u(x_0,y_0,z_0)-u(x,y,z)\equiv 0$成立,则$u(x,y,z)$在$V_R$上是一个常数.

2. 由于该有界区域是连通的,则不同于$(x_0,y_0,z_0)$的任意内点$(x,y,z)\in V$,存在一条完全含于$V$的曲线将$(x,y,z)$和$(x_0,y_0,z_0)$连接起来,并且可以用有限个完全含于$V$的球$V_1^{\ast },V_2^{\ast },\cdots .V_k^{\ast }$完全覆盖这条曲线,并且满足
   $$(x,y,z)\in V_1^{\ast },(x_0,y_0,z_0)\in V_k^{\ast },V_1^{\ast }\cap V_2^{\ast }̸=\varnothing ,V_2^{\ast }\cap V_3^{\ast }̸=\varnothing ,\cdots ,V_{k-1}^{\ast }\cap V_k^{\ast }̸=\varnothing$$
   根据1知$V_k^{\ast }$上$u$为常数,则$V_{k-1}^{\ast },\cdots ,V_2^{\ast },V_1^{\ast }$上$u$均为常数,特别地,$u(x,y,z)$也为常数.根据$(x,y,z)$的任意性,则$u(x,y,z)$在$V$上是一个常数.

3. 根据$u(x,y,z)$在$\overline{V}$上的连续性,通过$V$的点趋向$S=\partial V$上的点取极限,则$u(x,y,z)$在$\overline{V}$上为常数.

综上所述,$u$为常数,这与条件"非常数调和函数"矛盾,因而假设有误.假设最小值点同理可证,则$u(x,y,z)$的最值点不能在区域内部取得.

2.2 复连通区域上的Gauss公式

对于复连通区域,曲线的外法线方向始终由研究区域内部出发指向研究区域的外部,外法线方向定义为外侧.例如,对于满足
V = { ( x , y , z ) ∣ a 2 ≤ x 2 + y 2 + z 2 ≤ b 2 } S 1 = { ( x , y , z ) ∣ x 2 + y 2 + z 2 = a 2 } , S 2 = { ( x , y , z ) ∣ x 2 + y 2 + z 2 = b 2 } , a &lt; b , ∂ V = S 1 ∪ S 2 V=\{(x,y,z)|a^2\le x^2+y^2+z^2\le b^2\}\\S_1=\{(x,y,z)|x^2+y^2+z^2=a^2\},S_2=\{(x,y,z)|x^2+y^2+z^2=b^2\},a&lt;b,\partial V=S_1\cup S_2 V={(x,y,z)∣a2≤x2+y2+z2≤b2}S1​={(x,y,z)∣x2+y2+z2=a2},S2​={(x,y,z)∣x2+y2+z2=b2},a<b,∂V=S1​∪S2​
的区域$V$,边界$S_1$上的外法线单位向量为$$\vec{n}=(-\frac{x}{a},-\frac{y}{a},-\frac{z}{a})$$,而边界$S_2$上的外法线单位向量为$$\vec{n}=(\frac{x}{b},\frac{y}{b},\frac{z}{b})$$,在区域的内表面上指向原点,区域的外表面上背离原点方向.因而类似地有:

设$V$是由分片光滑的双侧封闭曲面$S_1,S_2,\cdots ,S_k$围成的闭区域,满足(1)向量场$F=(P(x,y,z),Q(x,y,z),R(x,y,z))$在$V$上连续;(2)$P(x,y,z),Q(x,y,z),R(x,y,z)$在$V$上有连续的一阶偏导数,设$\vec{n}$为$S_1,S_2,\cdots ,S_k$的外法线单位向量,用"+"表示取曲面的外侧,则
$${\iiint }_V\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dxdydz={\iint }_{S_1+S_2+\cdots +S_k}[P\cos (\vec{n},x)+Q\cos (\vec{n},y)+R\cos (\vec{n},z)]dS$$
证明同样省略.

空间Green第一公式

设$S_1,S_2,\cdots ,S_k$为封闭光滑曲面,$V$为$S_1,S_2,\cdots ,S_k$围成的区域,$u$在$S$和$V$上有二阶连续的偏导数,$\omega$有连续的偏导数,$\vec{n}$为$S_1,S_2,\cdots ,S_k$的外法线单位向量,则
$$\begin{gathered} {\iiint }_V\omega \frac{\partial u}{\partial x}dxdydz+{\iiint }_{V}u\frac{\partial \omega }{\partial x}dxdydz={\iint }_{S_1+S_2+\cdots +S_k}uwdydz \\ {\iiint }_V\omega \Delta udxdydz+{\iiint }_V\left(\frac{\partial u}{\partial x}\frac{\partial \omega }{\partial x}+\frac{\partial u}{\partial y}\frac{\partial \omega }{\partial y}+\frac{\partial u}{\partial z}\frac{\partial \omega }{\partial z}\right)dxdydz={\iint }_{S_1+S_2+\cdots +S_k}\omega \frac{\partial u}{\partial \vec{n}}dS \end{gathered}$$
特别地,当$\omega =u$时,有
$${\iiint }_V\left[{\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2+{\left(\frac{\partial u}{\partial z}\right)}^2\right]dxdydz+{\iiint }_{V}u\Delta udxdydz={\iint }_{S_1+S_2+\cdots +S_k}u\frac{\partial u}{\partial \vec{n}}dS$$
空间Green第二公式

空间Green第二公式

设$S_1,S_2,\cdots ,S_k$为封闭光滑曲面,$V$为$S_1,S_2,\cdots ,S_k$围成的区域,$u,v$在$S_1,S_2,\cdots ,S_k$和$V$上有二阶连续的偏导数,$\vec{n}$为$S_1,S_2,\cdots ,S_k$的外法线单位向量,则
$${\iiint }_V\left|\begin{array}{cc}\Delta u& \Delta v\\ u& v\end{array}\right|dxdydz={\iint }_{S_1+S_2+\cdots +S_k}\left|\begin{array}{cc}\frac{\partial u}{\partial \vec{n}}& \frac{\partial v}{\partial \vec{n}}\\ u& v\end{array}\right|dS$$

空间Green第三公式

若$u$为有界闭区域$V$内的调和函数,$S_1,S_2,\cdots ,S_k$为$V$的边界,$$\frac{\partial u}{\partial \vec{n}}$$为$u$沿着$S_1,S_2,\cdots ,S_k$的单位外法线的方向导数,$$\vec{r}=(a-x,b-y,c-z)$$,$$r=\sqrt{(a-x{)}^2+(b-y{)}^2+(c-z{)}^2}$$为$(x,y,z)$与$S$上动点$(a,b,c)$之间的距离,则
$$u(x,y,z)=\frac{1}{4\pi }{\iint }_{S_1+S_2+\cdots +S_k}\left(u\frac{\cos (\vec{r},\vec{n})}{r^2}+\frac{1}{r}\frac{\partial u}{\partial \vec{n}}\right)dS=\frac{1}{4\pi }{\iint }_{S_1+S_2+\cdots +S_k}\left(u\frac{\partial \frac{-1}{r}}{\partial \vec{n}}-\frac{-1}{r}\frac{\partial u}{\partial \vec{n}}\right)dS$$
注意
$${\iint }_{S_1+S_2+\cdots +S_k-S_0}\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS={\iiint }_{V-V_0}\left|\begin{array}{cc}\Delta \frac{-1}{r}& \Delta u\\ \frac{-1}{r}& u\end{array}\right|dadbdc=0$$

$${\iint }_{S_1+S_2+\cdots +S_k}\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS={\iint }_{S_0}\left|\begin{array}{cc}\frac{\partial \frac{-1}{r}}{\partial \vec{n}}& \frac{\partial u}{\partial \vec{n}}\\ \frac{-1}{r}& u\end{array}\right|dS$$

中值定理和极大值原理同样满足.

3 黎曼公式证明

$$\begin{gathered} {\iint }_S\left|\begin{array}{cc}L(u)& M(v)\\ u& v\end{array}\right|dxdy={\int }_{C}Pdx+Qdy \\ L(u)=\frac{{\partial }^{2}u}{\partial x\partial y}+a\frac{\partial u}{\partial x}+b\frac{\partial u}{\partial y}+cu \\ M(v)=\frac{{\partial }^{2}v}{\partial x\partial y}-a\frac{\partial v}{\partial x}-b\frac{\partial v}{\partial y}+cv \end{gathered}$$

其中$a,b,c$为常数,$C$围成的区域为$S$,$u,v$在$S$上连续并且有连续的二阶偏导数.$P,Q$为确定的函数.

证明 这个问题的证明需要用到格林公式.对于左边等式先凑形式.
$$\begin{array}{rl}& {\iint }_S\left|\begin{array}{cc}L(u)& M(v)\\ u& v\end{array}\right|dxdy\\ =& {\iint }_S\left[vL(u)-uM(v)\right]dxdy\\ =& {\iint }_S\left(v\frac{{\partial }^{2}u}{\partial x\partial y}+av\frac{\partial u}{\partial x}+bv\frac{\partial u}{\partial y}+cuv-u\frac{{\partial }^{2}v}{\partial x\partial y}+au\frac{\partial v}{\partial x}+bu\frac{\partial v}{\partial y}-cuv\right)dxdy\\ =& {\iint }_S\left[\left(v\frac{{\partial }^{2}u}{\partial x\partial y}+\frac{\partial v}{\partial x}\frac{\partial u}{\partial y}\right)+a\left(v\frac{\partial u}{\partial x}+u\frac{\partial v}{\partial x}\right)+b\left(v\frac{\partial u}{\partial y}+u\frac{\partial v}{\partial y}\right)-\left(u\frac{{\partial }^{2}v}{\partial x\partial y}+\frac{\partial u}{\partial y}\frac{\partial v}{\partial x}\right)\right]dxdy\\ =& {\iint }_S\left[\frac{\partial }{\partial x}\left(v\frac{\partial u}{\partial y}\right)+a\frac{\partial vu}{\partial x}+b\frac{\partial vu}{\partial y}-\frac{\partial }{\partial y}\left(u\frac{\partial v}{\partial x}\right)\right]dxdy\\ =& {\int }_C\left(u\frac{\partial v}{\partial x}-buv\right)dx+\left(v\frac{\partial u}{\partial y}+auv\right)dx\end{array}$$
则令
$$\begin{gathered} P=u\frac{\partial v}{\partial x}-buv \\ Q=v\frac{\partial u}{\partial y}+auv \end{gathered}$$
即证毕.

4 外微分形式和三大公式的统一形式

4.1 外积和外微分形式

以二重积分为例,如果认为面积元有向,在应用换元公式的时候不对雅可比行列式取绝对值,会得到这样一个结果
$$dxdy=\frac{\partial (x,y)}{\partial (u,v)}=\left|\begin{array}{cc}\frac{\partial x}{\partial u}& \frac{\partial x}{\partial v}\\ \frac{\partial y}{\partial u}& \frac{\partial y}{\partial v}\end{array}\right|dudv$$
如果交换$dx$和$dy$,则有
$$dydx=\frac{\partial (y,x)}{\partial (u,v)}=-\left|\begin{array}{cc}\frac{\partial x}{\partial u}& \frac{\partial x}{\partial v}\\ \frac{\partial y}{\partial u}& \frac{\partial y}{\partial v}\end{array}\right|dudv$$
默认以上等式成立的条件都满足,那么就有
$$dxdy=-dydx$$
在当时学习二重积分时,我们强制取了雅可比行列式的绝对值,得到了一个看起来很对称的结果.那么如果不取绝对值,可以考虑以下内容.

如果认为在$dx$和$dy$之间存在一种运算,使得其满足反对称性,即$dxdy=-dydx$,立刻可以得到$dxdx=0$.可以想到这种运算应该是和向量外积类似的运算,那么

在微分$dx,dy,dz,\cdots$之间定义一种外积运算,用记号来表示: 满足以下规定的微分乘积叫做外乘积
$$dx\wedge dy=-dy\wedge dx\Rightarrow dx\wedge dx=0$$
更严格的数学定义描述如下

1. $(adx)\wedge dy=a(dx\wedge dy),a\in R$
2. 外乘积运算对加法的分配律: $dx\wedge (dy+dz)=dx\wedge dy+dx\wedge dz$
3. 反交换律: $dx\wedge dy=-dy\wedge dx$
4. 结合律: $dx\wedge (dy\wedge dz)=(dx\wedge dy)\wedge dz$

这些运算定律和向量的外积运算是对应的.并且认为普通的两个不相同微分之间的乘积交换顺序是不变的,即$dxdy=dydx$.

由微分的外乘积和函数组成的线性组合称为外微分形式,具体地说,设$P,Q,R,A,B,C,H$都为$x,y,z$的函数,则

1. $Pdx+Qdy+Rdz$称为一次外微分形式(一次没有乘积,与普通的微分形式是一样的);
2. $Ady\wedge dz+Bdz\wedge dx+Cdx\wedge dy$称为二次外微分形式;
3. $Hdx\wedge dy\wedge dz$称为三次外微分形式;
4. 特别地,函数$f$成为零次外微分形式,$P,Q,R,A,B,C,H$称为微分形式的系数.

对于两个外微分形式$\lambda ,\mu$,也可以定义外积$\lambda \wedge \mu$,只需要将对应的各项外微分进行外积就可以了.如果设$P,Q,R,A,B,C,E,F,G,H$都为$x,y,z$的函数,那么
$$\begin{array}{rl}\lambda =& Adx+Bdy+Cdz\\ \mu =& Edx+Fdy+Gdz\\ \nu =& Pdx\wedge dy+Qdz\wedge dx+Rdx\wedge dy\\ \lambda \wedge \mu =& (Adx+Bdy+Cdz)\wedge (Edx+Fdy+Gdz)\\ =& (BG-CF)dx\wedge dy+(CE-AG)dz\wedge dx+(AF-BE)dx\wedge dy\\ \lambda \wedge \nu =& (Adx+Bdy+Cdz)\wedge (Pdx\wedge dy+Qdz\wedge dx+Rdx\wedge dy)\\ =& (AP+BQ+CR)dx\wedge dy\wedge dz\end{array}$$
进一步,可以给出如下运算律,设$\lambda ,\mu ,\nu$是三个任意的外微分形式,则有

1. 对加法的分配律: $(\lambda +\mu )\wedge \nu =\lambda \wedge \nu +\mu \wedge \nu ,\lambda \wedge (\mu +\nu )=\lambda \wedge \mu +\lambda \wedge \nu$
2. 结合律: $\lambda \wedge (\mu \wedge \nu )=(\lambda \wedge \mu )\wedge \nu$
3. 外微分形式的外积不满足交换律.但若$\lambda$是$p$次外微分形式,$\mu$是$q$次外微分形式,则$\lambda \wedge \mu =(-1{)}^{pq}\mu \wedge \lambda$

4.2 外微分运算与Poincaré引理

对外微分形式$\omega$引进外微分运算,用$d$表示,称$d$为外微分算子.在全文中出现的$d$格式都有些问题,为了方便起见,就不修改了.定义运算如下

1. 对零次外微分形式$\omega =f$ ,就是通常的全微分运算,
   $$d\omega =df=\frac{\partial f}{\partial x}dx+\frac{\partial f}{\partial y}dy+\frac{\partial f}{\partial z}dz$$

2. 对一次外微分形式$\omega =Pdx+Qdy+Rdz$,
   $$\begin{gathered} d\omega =dP\wedge dx+dQ\wedge dy+dR\wedge dz \\ d\omega =\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)dy\wedge dz+\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)dz\wedge dx+\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\wedge dy \end{gathered}$$

3. 对二次外微分形式$\omega =Ady\wedge dz+Bdz\wedge dx+Cdx\wedge dy$,
   $$\begin{gathered} d\omega =dA\wedge dy\wedge dz+dB\wedge dz\wedge dx+dC\wedge dx\wedge dy \\ d\omega =\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dx\wedge dy\wedge dz \end{gathered}$$

4. 对三次外微分形式$Hdx\wedge dy\wedge dz$,
   $$d\omega =dH\wedge dx\wedge dy\wedge dz=0$$

在这些规定下,外微分算子$d$与普通微分算子$d$就是一样的了,即对每一项进行运算,在每一项中分别对每一个因子运算,保持其余因子不动,将得出的各项相加.区别在于普通的微分将各项得出后进行普通的乘法,而外微分需要进行外积.

Poincaré引理叙述如下:

若$\omega$为一外微分形式,其微分形式的系数具有二阶连续的偏微商,则$d^2\omega =d(d\omega )=0$.

证明 仅在三维空间里讨论,而三维空间的外微分形式只有四种,不妨逐个验证.
$$\begin{array}{rl}dd\omega =ddf=& d\left(\frac{\partial f}{\partial x}\right)\wedge dx+d\left(\frac{\partial f}{\partial y}\right)\wedge dy+d\left(\frac{\partial f}{\partial z}\right)\wedge dz\\ =& \left(\frac{{\partial }^{2}f}{\partial y\partial z}-\frac{{\partial }^{2}f}{\partial z\partial y}\right)dy\wedge dz+\left(\frac{{\partial }^{2}f}{\partial z\partial x}-\frac{{\partial }^{2}f}{\partial x\partial z}\right)dz\wedge dx+\left(\frac{{\partial }^{2}f}{\partial x\partial y}-\frac{{\partial }^{2}f}{\partial y\partial x}\right)dx\wedge dy\\ =& 0\end{array}$$

$$\begin{array}{rl}dd\omega =& d\left[\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)dy\wedge dz+\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)dz\wedge dx+\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\wedge dy\right]\\ =& \left(\frac{{\partial }^{2}R}{\partial y\partial x}-\frac{{\partial }^{2}Q}{\partial z\partial x}+\frac{{\partial }^{2}P}{\partial z\partial y}-\frac{{\partial }^{2}R}{\partial x\partial y}+\frac{{\partial }^{2}Q}{\partial x\partial z}-\frac{{\partial }^{2}P}{\partial y\partial z}\right)dx\wedge dy\wedge dz\\ =& 0\end{array}$$

$$dd\omega =d(dA\wedge dy\wedge dz+dB\wedge dz\wedge dx+dC\wedge dx\wedge dy)=0$$

$$dd\omega =dd(Hdx\wedge dy\wedge dz)=d(0)=0$$

当然Poincaré引理还有逆定理:若$\omega$是一个$p,p=1,2,3$次外微分形式,且$d\omega =0$,则一定存在一个$p-1$次外微分形式$\alpha$,使得
$$\omega =d\alpha$$
证明 还是仅在三维空间内讨论.

若$p=1$,则有$\omega =Pdx+Qdy+Rdz$,并且
$$\begin{gathered} d\omega =\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)dy\wedge dz+\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)dz\wedge dx+\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\wedge dy=0 \\ \Rightarrow \frac{\partial R}{\partial y}=\frac{\partial Q}{\partial z},\frac{\partial P}{\partial z}=\frac{\partial R}{\partial x},\frac{\partial Q}{\partial x}=\frac{\partial P}{\partial y} \end{gathered}$$
取$$f(x,y,z)={\int }_0^{x}P(t,y,z)dt$$,则
$$\begin{array}{rl}df=& P(x,y,z)dx+{\int }_0^x\frac{\partial P(t,y,z)}{\partial y}dtdy+{\int }_0^x\frac{\partial P(t,y,z)}{\partial z}dtdz\\ =& P(x,y,z)dx+{\int }_0^x\frac{\partial Q(t,y,z)}{\partial t}dtdy+{\int }_0^x\frac{\partial R(t,y,z)}{\partial t}dtdz\\ =& P(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz-Q(0,y,z)dy-R(0,y,z)dz\end{array}$$
因而可得到,令
$$\alpha ={\int }_0^{x}P(t,y,z)dt+{\int }_0^{y}Q(0,t,z)dt+{\int }_0^{z}R(0,y,t)dt$$
即可.

若$p=2$,则有$\omega =Ady\wedge dz+Bdz\wedge dx+Cdx\wedge dy$,并且,
$$d\omega =\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dx\wedge dy\wedge dz=0\Rightarrow \frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}=0$$

$$\begin{gathered} d\left({\int }_0^{y}A(x,t,z)dtdz\right)=A(x,y,z)dy\wedge dz+{\int }_0^y\frac{\partial A(x,t,z)}{\partial x}dtdx\wedge dz \\ d\left({\int }_0^y-C(x,t,z)dtdx\right)=-C(x,y,z)dy\wedge dx-{\int }_0^y\frac{\partial C(x,t,z)}{\partial z}dtdz\wedge dx \end{gathered}$$

$$\begin{array}{rl}& d\left({\int }_0^{y}A(x,t,z)dtdz-{\int }_0^{y}C(x,t,z)dtdx\right)\\ =& A(x,y,z)dy\wedge dz+C(x,y,z)dx\wedge dy-{\int }_0^y\left(\frac{\partial A(x,t,z)}{\partial x}+\frac{\partial C(x,t,z)}{\partial z}\right)dtdz\wedge dx\\ =& A(x,y,z)dy\wedge dz+C(x,y,z)dx\wedge dy+{\int }_0^y\frac{\partial B(x,t,z)}{\partial t}dtdz\wedge dx\\ =& A(x,y,z)dy\wedge dz+B(x,y,z)dz\wedge dx+C(x,y,z)dx\wedge dy-B(x,0,z)dz\wedge dx\end{array}$$

令
$$\alpha ={\int }_0^{y}A(x,t,z)dtdz-{\int }_0^{y}C(x,t,z)dtdx+{\int }_0^{z}R(x,0,t)dtdx$$
若$p=3$,则有$\omega =Hdx\wedge dy\wedge dz$,令
$$\alpha ={\int }_0^{x}H(t,y,z)dtdy\wedge dz$$
即可.

当然以上定理需要满足一定的条件,当$p=1$要求定义域是曲面单连通的,当$p=2$要求定义域是空间单连通的,这是为了保证上述所有积分是有意义的.

4.3 梯度、旋度和散度的数学意义

对于零次外微分形式$\omega =f$,梯度和零次外微分形式的外微分相当,即
$$gradf=\frac{\partial f}{\partial x}i+\frac{\partial f}{\partial y}j+\frac{\partial f}{\partial z}k,d\omega =df=\frac{\partial f}{\partial x}dx+\frac{\partial f}{\partial y}dy+\frac{\partial f}{\partial z}dz$$
对一次外微分形式$\omega =Pdx+Qdy+Rdz$,旋度和一次外微分形式的外微分相当,即
$$\begin{gathered} rotv=\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)i+\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)j+\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)k \\ d\omega =\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z}\right)dy\wedge dz+\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x}\right)dz\wedge dx+\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\wedge dy \end{gathered}$$
对二次外微分形式$\omega =Ady\wedge dz+Bdz\wedge dx+Cdx\wedge dy$,散度和二次外微分形式的外微分相当,即
$$divv=\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z},d\omega =\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dx\wedge dy\wedge dz$$
由这个角度可以看出,三维空间不会再出现其他的"度".

那么来考虑$dd\omega =0$,若$\omega =f,dd\omega =0$相当于$f$梯度场必为无旋场,即
$$rotgradf=0$$
若$\omega =Pdx+Qdy+Rdz,dd\omega =0$相当于$v$的旋度场必为无源场,即
$$divrotv=0$$
考虑Poincaré引理,若$\omega$为一次外微分形式,相当于无旋场为有势场,存在一个势函数;若$\omega$为二次外微分形式,相当于无源场有向量势,存在一个向量势.

4.4 多变量微积分的基本定理(一般的Stokes公式)

Green公式,Gauss公式,Stokes公式可以写成一个统一的形式
$${\int }_{\partial \Omega }\omega ={\int }_{\Omega }d\omega$$
这里的$\omega$为外微分形式,$d\omega$为外微分运算,$\Omega$为$d\omega$的积分区域,$\partial \Omega$表示$\Omega$的边界,$$\int$$表示重积分,其维数和区域的维数相等.这个公式可以推广到更高维的空间中去.更详细的理论可以参见Wikipedia.

附录:

这次大作业的前半部分主要讨论了平面和空间上Green公式的和调和函数的相关定理性质.在整个的证明体系中,Green和Gauss公式在课堂上已经证明过,而Green第一公式、Green第二公式、Green第三公式、调和函数的中值定理和极大值原理是给出证明的,并且在证明过程中使用了前一条定理的结论,使得整个证明形成一套体系,可以更好地认识和理解Green和Gauss公式.从中可以看到,平面和空间两种情况下的定理证明有极为相似的内在思想和外在形式.如果用外法线单位向量的形式来表述整个定理体系,可以认为,平面情况和空间情况基本上"没有什么不同".在单连通区域和复连通区域,如果也用外法线单位向量的形式来表述整个定理,得到的结果也极为相似.那么,如何定义这个"外法线单位向量"就成了很关键的问题.定义好外法线单位向量的含义,可以简明而系统地表述这些定理.此外,Green第三公式在平面和空间两种情况下有着相似的形式,相同的积分和行列式形式,不同的系数和关于$r$的函数,这其中是否有必然的联系?这是一个可以继续研究和尝试推广到$n$维空间下的问题.

如果再加上Stokes公式,这三个公式类别具有着相似的特性,这就启发我们思考这几个公式是否可以统一?也就是说,这几个公式具有一定的内在联系.事实上,在本章的最后的确给出了一个统一的表述.那么高维空间是否也有相似的结论呢?这也是一个可以继续推广的问题.

参考文献

1. 杨小远,孙玉泉,杨卓琴,薛玉梅编著.工科数学分析分析教程(下册)[M].北京:科学出版社,2012.

2. 费定晖,周学圣编演.吉米多维奇数学分析习题集题解6[M].济南:山东科学技术出版社,2012.

3. 龚昇.简明微积分[M].北京:高等教育出版社,2006.

4. 中国科学技术大学数学科学学院编著.微积分学导论(下册)[M].合肥:中国科学技术大学出版社,2016.

5. 李雪平.关于外微分形式[J].韶关师专学报,1981,(02):14-20.

6. wikipedia.Master theorem[EB/OL].2017-06-15.

   https://en.wikipedia.org/wiki/Exterior_derivative

7. wikipedia.Green theorem[EB/OL].2017-06-15.

   https://en.wikipedia.org/wiki/Green's_theorem

8. wikipedia.Stokes theorem[EB/OL].2017-06-15.

   https://en.wikipedia.org/wiki/Stokes'_theorem