重积分（微积分）

由一元积分演化而来

• 定积分
  • 二重积分：直角坐标、极坐标
  • 三重积分：柱面坐标、球坐标、直角坐标

1 二重积分的概念与性质

一、概念

引入：曲顶柱体的体积、平面薄片的质量

• 定义1（分割求和取极限）：$\underset{D}{\iint }f(x,y)d\sigma =\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f(x_i,y_i)\Delta {\sigma }_i$（和式的极限），其中$D$为积分区域、$x,y$为积分变量、$f(x,y)$为被积函数、$d\sigma$为面积元素、$f(x,y)d\sigma$为被积表达式、$\sum f(x_i,y_i)\Delta {\sigma }_i$为积分和、 λ = max ⁡ { d i } 且 d i = max ⁡ ( x 1 , y 1 ) , ( x 2 , y 2 ) → Δ τ i { ∣ P 1 P 2 ∣ } \lambda=\max\{d_i\}且d_i=\max\limits_{(x_1,y_1),(x_2,y_2)\rightarrow\Delta\tau_i}\{|P_1P_2|\} λ=max{di​}且di​=(x1​,y1​),(x2​,y2​)→Δτi​max​{∣P1​P2​∣}.
  类比于定积分：${\int }_a^{b}f(x)dx=\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f(x_i)\Delta x_i$.
• 说明：
  • 连续必可积：$f(x,y)$在闭区间连续，二重积分必存在
  • 面积元素：$d\sigma =dxdy$
  • 几何意义：$f(x,y)$在$D$上的二重积分等于柱体体积的代数和（有正有负）

二、性质

• 可加，可数乘
• 单调，特别的$|\underset{D}{\iint }f(x,y)d\sigma |\le \underset{D}{\iint }|f(x,y)|d\sigma$.
• 估值不等式：$min\sigma \le \underset{D}{\iint }f(x,y)d\sigma \le MAX\sigma$.
• 中值定理：存在点$(x_0,y_0)$，使得$\underset{D}{\iint }f(x,y)d\sigma =f(x_0,y_0)\sigma$.

2 二重积分的计算

一、直角坐标

• X型：$V=$${\int }_a^b\left[{\int }_{{\phi }_1(x)}^{{\phi }_2(x)}f(x,y)dy\right]dx={\int }_a^{b}dx{\int }_{{\phi }_1(x)}^{{\phi }_2(x)}f(x,y)dy$.
  Y型：$V=$${\int }_c^d\left[{\int }_{{\psi }_1(y)}^{{\psi }_2(y)}f(x,y)dx\right]dy={\int }_c^{d}dy{\int }_{{\psi }_1(y)}^{{\psi }_2(y)}f(x,y)dx$.
• 分割，选择次序
• 做题时：
  • 只要见到积分区域（有时需要拆分区域）具有对称性，就要想到考察被积函数的奇偶性。
  • 考虑x、y能否轮换对称（关于y=x对称）。
  • 直接积分图形复杂时，考虑可加（减）性 。

例1：$求\underset{D}{\iint }{x}^2{e}^{-y^2}dxdy,其中D是以(0，0),(1,1),(0,1)为顶点的三角形$
解：
$$\begin{gathered} \because \int e^{-y^2}dy无法解出\therefore 先积分x \\ 原式={\int }_0^{1}dy{\int }_0^y{x}^2{e}^{-y^2}dx={\int }_0^1{e}^{-y^2}\cdot \frac{y^3}{3}dy={\int }_0^1{e}^{-y^2}\cdot \frac{y^2}{6}d{y}^2=-\frac{1}{6}{\int }_0^{1}ud{e}^{-u}=-\frac{1}{3e} \end{gathered}$$
例2：$求\underset{D}{\iint }y\sqrt{1+x^2-y^2}d\sigma ,其中D是由y=1,x=-1及y=x所围成的区域$
解：
$$\begin{gathered} 不妨以X型区域计算：原式={\int }_{-1}^{1}dx{\int }_x^{1}y\sqrt{1+x^2-y^2}dy \\ =...=-\frac{1}{3}{\int }_{-1}^1{\left[(1+x^2-y^2{)}^{\frac{3}{2}}\right]}_x^{1}dx=-\frac{1}{3}{\int }_{-1}^1(|x{|}^3-1)dx=-\frac{2}{3}{\int }_0^1(x^3-1)dx=\frac{1}{2}. \\ 不知x正负,平方开根号要绝对值，不然答案为奇函数错误 \end{gathered}$$
例3：$求两个底圆半径都等于R的直交援助面积所围成的立体的体积$
解：
$$\begin{gathered} 设两圆柱面方程为x^2+y^2=R^2和x^2+z^2=R^2 \\ V=8{\int }_0^{R}dx{\int }_0^{\sqrt{R^2-x^2}}\sqrt{R^2-x^2}dy=8{\int }_0^R(R^2-x^2)dx=\frac{16}{3}{R}^3. \end{gathered}$$

例4：$更换积分次序，I={\int }_0^{4}dy{\int }_{-\sqrt{4-y}}^{\sqrt{4y-y^2}}f(x,y)dx$.
解：
$$I={\int }_{-2}^{0}dx{\int }_0^{4-x^2}f(x,y)dy+{\int }_0^{2}dx{\int }_{2-\sqrt{4-x^2}}^{2+\sqrt{4-x^2}}f(x,y)dy.$$

例5：$计算z=xy,x+y+z=1,z=0所围闭合区域的体积$.
解：法1
$$\begin{gathered} 分为两部分，一部分z=xy为顶，另一部分z=1-x-y为顶 \\ \because xy=1-x-y⟺y=\frac{1-x}{1+x} \\ \therefore V={\int }_0^{1}dx{\int }_0^{\frac{1-x}{1+x}}xydy+{\int }_0^{1}dx{\int }_{\frac{1-x}{1+x}}^{1-x}1-x-ydy. \end{gathered}$$

​ 法2：参考例9，截面法
例6：$平面区域D由曲线\{\begin{array}{l}x=t-\sin t,\\ y=1-\cos t\end{array}(0\le t\le 2\pi )与x轴围成，计算\underset{D}{\iint }(x+2y)dxdy$.
$$\begin{gathered} 将参数方程选做换元积分的变量替换 \\ 设D:0\le x\le 2\pi ,0\le y\le g(x),则 \\ \underset{D}{\iint }(x+2y)dxdy={\int }_0^{2\pi }dx{\int }_0^{g(x)}(x+2y)dy={\int }_0^{2\pi }[xg(x)+g^2(x)]dx \\ ={\int }_0^{2\pi }[(t-\sin t)(1-\cos t)+(1-\cos t{)}^2](1-\cos t)dt=...=3{\pi }^2+5\pi . \end{gathered}$$

二、极坐标

• 由于$\Delta {\sigma }_i\approx {\rho }_i\Delta {\rho }_i\Delta {\theta }_i$，元素面积可以表示为$d\sigma =\rho d\rho d\theta$.
  因此$\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f(x_i,y_i)\Delta {\sigma }_i=\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f({\rho }_i\cos {\theta }_i,{\rho }_i\sin {\theta }_i){\rho }_i\Delta {\rho }_i\Delta {\theta }_i$，
  即==$\underset{D}{\iint }f(x,y)d\sigma =\underset{D}{\iint }f(\rho \cos \theta ,\rho \sin \theta )\rho d\rho d\theta$==.
• 其中按顺序：$\underset{D}{\iint }f(\rho \cos \theta ,\rho \sin \theta )\rho d\rho d\theta ={\int }_{\alpha }^{\beta }d\theta {\int }_{{\phi }_1(\theta )}^{{\phi }_2(\theta )}f(\rho \cos \theta ,\rho \sin \theta )\rho d\rho$.

例7：$转换成直角坐标系I={\int }_0^{\frac{\pi }{2}}d\theta {\int }_{a\cos \theta }^{2a\cos \theta }\sqrt{4a^2-{\rho }^2}\rho d\rho$.
解：
$$难点：画图I=I={\int }_0^{a}dx{\int }_{\sqrt{ax-x^2}}^{\sqrt{2ax-x^2}}\sqrt{4a^2-x^2-y^2}dy+{\int }_0^{2a}dx{\int }_0^{\sqrt{2ax-x^2}}\sqrt{4a^2-x^2-y^2}dy.$$

例8：$闭区域D:x^2+y^2\le y(1),x\ge 0(2),f(x,y)为D上连续函数，且f(x,y)=\sqrt{1-x^2-y^2}-\frac{8}{\pi }\underset{D}{\iint }f(u,v)dudv(\ast ),求f(x,y)$.
解：
$$\begin{gathered} 设A=\underset{D}{\iint }f(u,v)dudv,又\underset{D}{\iint }dxdy=\frac{\pi }{8}, \\ 则在(\ast )两边求二重积分得到：A=\underset{D}{\iint }\sqrt{1-x^2-y^2}dxdy-\underset{D}{\iint }\frac{8}{\pi }Adxdy. \\ \therefore A=\frac{1}{2}\underset{D}{\iint }\sqrt{1-x^2-y^2}dxdy=\frac{1}{2}{\int }_0^{\frac{\pi }{2}}d\theta {\int }_0^{\sin \theta }\rho \sqrt{1-{\rho }^2}d\rho =\frac{1}{2}{\int }_0^{\frac{\pi }{2}}{\left[-\frac{1}{3}(1-{\rho }^2{)}^{\frac{3}{2}}\right]}_0^{\sin \theta }d\theta =\frac{\pi }{12}-\frac{1}{9}. \\ \therefore f(x,y)=\sqrt{1-x^2-y^2}-\frac{8}{\pi }(\frac{\pi }{12}-\frac{1}{9}). \end{gathered}$$

三、坐标变换

3 三重积分

一、概念

定义1：$\underset{\Omega }{\iiint }f(x,y,z)d\sigma =\underset{\lambda \to 0}{\lim }\sum _{i=1}^{n}f(x_i,y_i,z_1)\Delta {\sigma }_i$.
三重积分也记作：$\underset{\Omega }{\iiint }f(x,y,z)d\sigma =\underset{\Omega }{\iiint }f(x,y,z)dxdydz$.

二、计算

• 先一后二（投影法”压薄“）：
  $$\underset{\Omega }{\iiint }f(x,y,z)dxdydz=\underset{D}{\iint }F(x,y)dxdy=\underset{D}{\iint }\left[{\int }_{z_1(x,y)}^{z_2(x,y)}f(x,y,z)dz\right]dxdy$$

• 先二后一（截面法”集中“）：
  $$\underset{\Omega }{\iiint }f(x,y,z)dxdydz={\int }_{C_1}^{C_2}F(x)dz={\int }_{C_1}^{C_2}\left[\underset{D_z}{\iint }f(x,y,z)dxdy\right]dz$$
  例9：$计算I=\underset{\Omega }{\iiint }{e}^{|z|}dxdydz，其中\Omega 为x^2+y^2+z^2\le 1所围成的立体$.
  解：
  $$\begin{gathered} \because e^{|z|}为偶函数，采用先二后一的思想， \\ I=2{\int }_0^1{e}^{|z|}dz\underset{D_{xy}}{\iint }dxdy.其中\underset{D_{xy}}{\iint }dxdy=\pi (r(z){)}^2,r(z)=\sqrt{1-z^2}. \\ \therefore I=2{\int }_0^1{e}^z(1-z^2)dz=2\pi {\left[-(z-1{)}^2{e}^z\right]}_0^1=2\pi . \end{gathered}$$
  例10：$计算I=\underset{\Omega }{\iiint }xdxdydz，其中\Omega 由z=xy,x+y+z=1及z=0所围成的立体$.
  解：
  $$\begin{gathered} 发现当x为一定值时,平行于O_{yz}面的切面切立体的截面为一过原点的三角形,则原式表示成I={\int }_0^{1}xdx\underset{D_{yz}}{\iint }S(x)dydz. \\ 对于S(x),当x=k时，则z=ky,z=1-k-y,z=0，\sin \theta =\frac{k}{\sqrt{1+k^2}} \\ 有S=\frac{1}{2}(1-k)(\sqrt{1+k^2}\frac{1-k}{1+k})\sin \theta =\frac{k(1-k{)}^2}{2(1+k)}. \\ \therefore I=\frac{1}{2}{\int }_0^1\frac{x(1-x{)}^2}{1+x}dx+1=\frac{1}{2}{\int }_1^2\frac{(t-2{)}^2(t-1{)}^2}{t}dt=\frac{1}{2}{\int }_0^1{t}^3-6t^2+13t-12+\frac{4}{t}dt=2ln2-\frac{11}{8}. \end{gathered}$$
  

• 利用柱面坐标（case：被积函数具有$f(\sqrt{x^2+y^2})$的形式，积分区域是直柱体，且投影面用极坐标表示方便）：
  $$\{\begin{array}{l}x=\rho \cos \theta .\\ y=\rho \sin \theta .\\ z=z.\end{array}$$
  $\underset{\Omega }{\iiint }f(x,y,z)dxdydz=\underset{\Omega }{\iiint }f(\rho \cos \theta ,\rho \sin \theta ,z)\rho d\rho d\theta dz={\int }_{\alpha }^{\beta }d\theta {\int }_{{\phi }_1}^{{\phi }_2}\rho d\rho {\int }_{z_1}^{z_2}f(\rho \cos \theta ,\rho \sin \theta ,z)dz$

  例11：$将I=\underset{\Omega }{\iiint }f(x,y,z)dxdydz表示成柱面坐标系的累次积分，其中\Omega 是由2z=x^2+y^2,z=2,z=1所围成的立体$.
  解：
  $$I={\int }_0^{2\pi }d\theta {\int }_0^{\sqrt{2}}\rho d\rho {\int }_1^{2}f(\rho \cos \theta ,\rho \sin \theta ,z)dz+{\int }_0^{2\pi }d\theta {\int }_{\sqrt{2}}^2\rho d\rho {\int }_{\frac{{\rho }^2}{2}}^{2}f(\rho \cos \theta ,\rho \sin \theta ,z)dz.$$
  

• 利用球面坐标（case：被积函数具有$f(\sqrt{x^2+y^2+z^3})$的形式，积分区域是球面或锥面）：
  $$\{\begin{array}{l}x=r\sin \phi \cos \theta .\\ y=r\sin \phi \sin \theta .\\ z=r\cos \phi .\end{array}$$
  $\underset{\Omega }{\iiint }f(x,y,z)dxdydz=\underset{\Omega }{\iiint }f(r\sin \phi \cos \theta ,r\sin \phi \sin \theta ,r\cos \phi )r^2\sin \phi drd\phi d\theta$.

  例12：$将I=\underset{\Omega }{\iiint }f(x,y,z)dxdydz表示成球面坐标系的累次积分，其中\Omega 是由2z=x^2+y^2,z=2,z=1所围成的立体$.
  解：
  $$\begin{gathered} I={\int }_0^{2\pi }d\theta {\int }_0^{\frac{\pi }{4}}\sin \phi d\phi {\int }_{\frac{1}{\cos \phi }}^{\frac{2}{\cos \phi }}f(r\sin \phi \cos \theta ,r\sin \phi \sin \theta ,r\cos \phi )r^{2}dr \\ +{\int }_0^{2\pi }d\theta {\int }_{\frac{\pi }{4}}^{\arctan \sqrt{2}}\sin \phi d\phi {\int }_{\frac{1}{\cos \phi }}^{\frac{2\cos \phi }{{\sin }^2\phi }}f(r\sin \phi \cos \theta ,r\sin \phi \sin \theta ,r\cos \phi )r^{2}dr. \end{gathered}$$
  

4 重积分的应用

一、曲面面积

• 由定义可知：
  • $\underset{D}{\iint }1d\sigma =$平面区域面积
  • $\underset{\Omega }{\iiint }1dV=\underset{D_{xy}}{\iint }zd\sigma =$空间区域体积
• 二重积分计算曲面区域$z=f(x,y)$的面积:
  由于$\vec{n}=(f_x,f_y,-1)$且$|\cos \alpha |=\vec{n}\times (0,0,1)$，故$|\cos \alpha |=\frac{1}{\sqrt{1+f_x^2(x,y)+f_y^2(x,y)}}$
  得到 $S=\underset{D}{\iint }\sqrt{1+f_x^2(x,y)+f_y^2(x,y)}d\sigma$.

二、质心

• 质心：使物体上的点$P$对空间任意一点$P_0$力矩积分为0，即$\underset{\Omega }{\iiint }(\overrightarrow{OP}-\overrightarrow{O{P}_0})\rho (x,y,z)dV=0$。
• 故质心坐标$(x_0,y_0,z_0)$满足：$\{\begin{array}{l}\underset{\Omega }{\iiint }(x-x_0)\rho (x,y,z)dV=0\\ \underset{\Omega }{\iiint }(y-y_0)\rho (x,y,z)dV=0\\ \underset{\Omega }{\iiint }(z-z_0)\rho (x,y,z)dV=0\end{array}$，
  即
  $$\{\begin{array}{l}{x}_0=\frac{\underset{\Omega }{\iiint }x\rho (x,y,z)dV}{m}\\ y_0=\frac{\underset{\Omega }{\iiint }y\rho (x,y,z)dV}{m}\\ z_0=\frac{\underset{\Omega }{\iiint }z\rho (x,y,z)dV}{m}\end{array}$$
  其中$m=\underset{\Omega }{\iiint }\rho (x,y,z)dV$即物体的质量。
• 假如物体为一薄板，且水平放置在Oxy面上，则
  $$\{\begin{array}{l}{x}_0=\frac{\underset{D}{\iint }x\rho (x,y)d\sigma }{m}\\ y_0=\frac{\underset{D}{\iiint }y\rho (x,y)d\sigma }{m}\end{array}$$

三、转动惯量

• 整个物体对于坐标轴的转动变量：
  例如对于$z$轴：$J_z=\underset{\Omega }{\iiint }(x^2+y^2)\rho (x,y,z)dV$.
• 整个物体对于一点$P_0(x_0,y_0,z_0)$的转动变量：$J_{P_0}=\underset{\Omega }{\iiint }\left[(x-x_0{)}^2+(y-y_0{)}^2+(z-z_0{)}^2\right]\rho (x,y,z)dV$，
  特别的，物体关于原点的转动惯量为$J_0=\underset{\Omega }{\iiint }(x^2+y^2+z^2)\rho (x,y,z)dV$.
• 由于一个质点关于一个平面的转动惯量为质量乘以质点到平面最短距离之平方，故整个物体对于坐标面的转动惯量：
  例如对于$Oxy$面：$J_{xy}=\underset{\Omega }{\iiint }{z}^2\rho (x,y,z)dV$.

四、引力

• 物体对于质点$P_0$的总引力为
  $$F=G\underset{\Omega }{\iiint }\frac{M_{0}dm}{r^2}\cdot r_0=G\underset{\Omega }{\iiint }\frac{M_0\cdot \vec{r}(x,y,z)\cdot \rho (x,y,z)dV}{r^3}$$
  其中$r_0=\frac{\vec{r}(x,y,z)}{r}$为单位向量，$r=\sqrt{(x-x_0{)}^2+(y-y_0{)}^2+(z-z_0{)}^2}$，$\vec{r}(x,y,z)=(x-x_0,y-y_0,z-z_0)$，则坐标分量为
  $$\{\begin{array}{l}{F}_x=G{M}_0\underset{\Omega }{\iiint }\frac{x-x_0}{r^3}\rho (x,y,z)dV\\ F_y=G{M}_0\underset{\Omega }{\iiint }\frac{y-y_0}{r^3}\rho (x,y,z)dV\\ F_z=G{M}_0\underset{\Omega }{\iiint }\frac{z-z_0}{r^3}\rho (x,y,z)dV\end{array}$$