查看“︁多重积分”︁的源代码

{{NoteTA
|1=zh-cn:变量;zh-tw:變數;
}}

{{微積分學}}
[[File:Areabetweentwographs.svg|thumb|该积分是两条曲线之间的面积。]]
[[File:Volume under surface.png|thumb|作为曲面下的体积的双重积分。该体积底部的矩形区域是积分的域，而曲面是被积的双变量函数的图像。]]
'''多重积分'''（{{lang-en|Multiple integral}}）是定积分的一类，它将[[定积分]]扩展到多元函数（多[[变量]]的[[函数]]），例如求<math>f(x,y)</math>或者<math>f(x,y,z)</math>类型的多元函数的积分。

== 简介 ==
正如单变量的正函数的定积分代表函数图像和''x''轴之间区域的[[面积]]一样，正的双变量函数的'''双重积分'''代表函数所定义的曲面和包含函数[[定义域]]的平面之间所夹的区域的[[体积]]。（注意同样的体积也可以通过三变量常函数''f''(''x'', ''y'', ''z'') = 1在上述曲面和平面之间的区域中的'''三重积分'''得到。若有更多变量，则[[多元函数]]的多重积分给出[[超体积]]。

''n''元函数''f''(''x''<sub>1</sub>, ''x''<sub>2</sub>,…, ''x''<sub>''n''</sub>)在定义域''D''上的多重积分通常用嵌套的积分号按照演算的逆序标识（最左边的积分号最后计算），后面跟着[[被积函数]]和正常次序的[[积分变量]]（最右边的变量最后使用）。积分域或者对每个积分变量在每个积分号下标识，或者用一个变量标在最右边的积分号下：

:<math>\int \ldots \int_\mathbf{D}\;f(x_1,x_2,\ldots,x_n) \;\mathrm{d}x_1 \!\ldots\mathrm{d}x_n</math>

因为不可能计算[[多元函数|多于一个自变量的函数]]的[[不定积分]]，“不定”多重积分是不存在的。因此所有多重积分都是“定”积分。

通常在坐标系中，多重积分都利用嵌套的[[累次积分]]计算。而累次积分为了简便可记为：

:<math>\int_{\varphi_1}^{\psi_1}\,\mathrm{d}x_1\,\int_{\varphi_2(x_1)}^{\psi_2(x_1)}\,\mathrm{d}x_2\!\dots\int_{\varphi_n(x_1,x_2,\dots,x_{n-1})}^{\psi_n(x_1,x_2,\dots,x_{n-1})}\,f(x_1,x_2,\dots,x_n)\mathrm{d}x_n</math>

其中积分域为：
:<math>D=\left\{(x_1,x_2,\dots,x_n)|\varphi_1\le x_1\le\psi_1,\varphi_2(x_1)\le x_2\le\psi_2(x_1),\dots,\varphi_n(x_1,x_2,\dots,x_{n-1})\le x_n\le\psi_n(x_1,x_2,\dots,x_{n-1})\right\}</math>
注意的是，该式一般情况下并不表示多个定积分的积，在实际计算中从最右侧积分变量开始积分，其结果会作为外一层积分的被积函数。

== 范例 ==
譬如，边长为4 × 6 × 5的[[长方体]]的体积可以通过两种方法得到：
* 通过函数''f''(''x'', ''y'') = 5在''xy''平面中的区域''D''，也就是长方体的底上的双重积分
:: <math>\iint_\mathrm{D} 5 \ \mathrm{d}x\, \mathrm{d}y</math>
* 或者是常函数1在长方体上的三重积分
::<math>\iiint 1\mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z</math> (1可以省略不写)

== 数学定义 ==

令''n''为大于1的自然数。考虑所谓的''半开'''''n维矩形'''（下面简称'''矩形'''）。对于[[平面 (数学)|平面]]，<math>n=2</math>。

: <math>T=[a_1, b_1)\times [a_2, b_2)\times\cdots \times [a_n, b_n)\subset \mathbb R^n</math>。

将每个区间[''a''<sub>''i''</sub>, ''b''<sub>''i''</sub>)分成有限个不重叠的子区间，每个都是左闭右开。将子区间记为''I''<sub>''i''</sub>。则，所有所有如下形式的子矩形的族

: <math>C=I_1\times I_2\times \cdots\times I_n</math>

是''T''的一个[[划分]]，也即，子矩形''C''是互不重叠的，而且它们聯集为''T''。''C''中的子矩形的'''直径'''按照定义是''C''中最大的边长，而''T''的划分的直径定義成划分中的所有子矩形的最大直径。

令''f'' : ''T'' → '''R'''为定义在''T''上的函数。考虑如上定义的''T''的划分

:<math>T=C_1\cup C_2\cup \cdots \cup C_m</math>

其中''m''是正整数。如下形式的和称为'''[[黎曼积分|黎曼和]]'''

:<math>\sum_{k=1}^m f(P_k)\, \operatorname{m}(C_k)</math>

其中，对于每个''k''，点''P''<sub>''k''</sub>在''C''<sub>''k''</sub>中，而m(''C''<sub>''k''</sub>)是笛卡尔积为''C''<sub>''k''</sub>的区间的边长之积。

函数''f''称为'''黎曼可积'''，如果如下[[极限]]存在

:<math>S = \lim_{\delta \to 0} \sum_{k=1}^m f(P_k)\, \operatorname{m}\, (C_k)</math>

其中极限取遍所有直径小於δ的''T''的划分。若''f''[[黎曼积分|黎曼可积]]，''S''称为''f''在''T''上的'''[[黎曼积分]]'''。记为

: <math>\int_T\!f(x)\,\mathrm{d}x</math>。

定义在任意有界''n''维集合上的函数的黎曼和可以通过将函数延拓到一个半开半闭矩形上来求出，其取值在原来的定义域之外为0。然后，原来的函数的积分就定义为延展的函数在矩形区域中的积分（如果存在的话）。

下文中''n''维黎曼积分简称'''多重积分'''。

=== 性质 ===
多重积分具有很多与单变量函数的积分一样的性质（线性，可加性，单调性，等等）。而且，和单变量情况一样，可以用多重积分找出函数在给定集合上的积分。具体来讲，给定集合''D'' ⊆ '''R'''<sup>''n''</sup>和''D''上的可积函数''f''，''f''在定义域上的平均值为

:<math>\bar{f} = \frac{1}{m(D)} \int_D f(x)\, \mathrm{d}x,</math>

其中''m''（''D''）是''D''的[[测度]]。

=== 特例 ===

''T'' ⊆ '''R'''<sup>2</sup>时，积分

:<math>\ell = \iint_T f(x,y)\, \mathrm{d}x\, \mathrm{d}y </math>

是''f''在''T''上的'''双重积分'''，而若''T'' ⊆ '''R'''<sup>3</sup>，积分

:<math>\ell = \iiint_T f(x,y,z)\, \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z</math>

是''f''在''T''上的'''三重积分'''。

注意，按常规，双重积分用两个积分号，而三重积分有三个；这只是记法上方便，也是为了通过重复积分来计算多重积分（参看本条目后文）。

== 积分方法 ==
多重积分问题的解决在多数情况下依赖于将多重积分转化为一系列单变量积分，而其中每个单变量积分都是直接可解的。

=== 直接检验 ===
有时可以直接获得积分的结果，而无需任何直接计算。

==== 常数 ====
在[[常函数]]的情况中，结果很直接：只要将常函数''c''乘以测度就可以了。如果''c'' = 1，而且是在'''R'''<sup>2</sup>的子集中积分，则乘积就是区域面积，而在'''R'''<sup>3</sup>中，它就是区域的体积。

* 例如：

::<math>D = \{ (x,y) \in \mathbb{R}^2 \ : \ 2 \le x \le 4 \ ; \ 3 \le y \le 6 \}</math> and <math>f(x,y) = 2\,\!</math>

:在''D''上积分''f''：

::<math>\int_3^6 \int_2^4 \ 2 \ \mathrm{d}x\, \mathrm{d}y = \mbox{area}(D) \cdot 2 = (2 \cdot 3) \cdot 2 = 12</math>。

==== 利用可能的对称性 ====
如果定义域存在沿着某条轴的对称性而且函数对于那个变量是奇函数，则积分为0（因为相反的两部分加起来为0）。

对于'''R'''<sup>''n''</sup>中的函数，只要相关变量对于形成对称的轴是奇变量就可以了。
<blockquote>例一:
:给定<math>f(x,y) = 2 \sin(x)-3y^3+5</math>以及<math>T=\left \{ ( x,y) \in \mathbf{R}^2 \ : \ x^2+y^2\le 1 \right \}</math>为积分区域（[[半径]]为1的圆盘，包含[[边界]]）。

:利用线性性质，积分可以分解为三部分：

:<math>\iint_T (2\sin x - 3y^3 + 5) \, \mathrm{d}x \, \mathrm{d}y = \iint_T 2 \sin x \, \mathrm{d}x \, \mathrm{d}y - \iint_T 3y^3 \, \mathrm{d}x \, \mathrm{d}y + \iint_T 5 \, \mathrm{d}x \, \mathrm{d}y</math>

:<math display="inline">2\sin(x)</math> 和 <math display="inline">3y^3</math> 都是奇函数，而且显然''T''对于''x''和''y''轴都是对称的；因此唯一有贡献的部分是常函数''5''因为其它两个都贡献0.</blockquote>
<blockquote>例二：
:考虑函数<math display="inline">f(x,y,z)=x \exp(y^2+z^2)</math>以及圆心在原点的半径为2的[[球]]
:<math>T = \left \{ ( x,y, z) \in \mathbf{R}^3 \ : \ x^2+y^2+z^2 \le 4 \right \}.</math>
:该球显然是对于三条轴都对称，但是只要对于''x''轴积分就可以看出结果是0，因为''f''对于该变量是奇函数。</blockquote>

=== 简化公式 ===
简化公式基于简单积分区域来将多重积分转化为单变量积分的序列。它们必须从右至左计算，过程中将其它变量暂时视为常数（和[[偏导数]]的计算类似）。

==== R<sup>2</sup>中的常规区域 ====
{{seealso|积分次序}}
此種方法適用於滿足下述條件的任何定義域 ''D'':
* ''D'' 投影到 ''x''軸或 ''y''軸任一軸，形成一個有邊界的範圍, 以 ''a'', ''b'' 代表邊界值。
* 通過 ''a'', ''b'' 兩點並與<math>\overline{ab}</math> 垂直的直線與 ''D'' 相交後的兩個端點，可以用 2 個函數<math>\alpha</math>, <math>\beta</math> 定義。
===== ''x''轴 =====
將 ''D'' 對 ''x''軸做垂直投影，函數<math>f: D \longrightarrow \mathbb{R}</math>是连续函数，并且''D''可以视为（定义在[''a'',''b'']区间上的）α(''x'')和β(''x'')之间的区域。则

:<math>\iint_D f(x,y)\ \mathrm{d}x\, \mathrm{d}y = \int_a^b \mathrm{d}x \int_{ \alpha (x)}^{ \beta (x)} f(x,y)\, \mathrm{d}y</math>

===== ''y''轴 =====
將 ''D''對''y''軸做垂直投影，函數<math>f: D \longrightarrow \mathbb{R}</math>是连续函数，并且''D''可以视为（定义在[''a'',''b'']区间上的）α(''y'')和β(''y'')之间的区域。则

:<math>\iint_D f(x,y)\ \mathrm{d}x\, \mathrm{d}y = \int_a^b \mathrm{d}y \int_{ \alpha (y)}^{ \beta (y)} f(x,y)\, \mathrm{d}x</math>

===== 范例 =====
[[File:Esempio-formulediriduzione-r2.svg|thumb|160px|例：可以采用简化公式的D区域。]]

考虑区域：<math>D = \{ (x,y) \ : \ x \ge 0, y \le 1, y \ge x^2 \}</math>（参看附图）。计算

:<math>\iint_D (x+y) \, \mathrm{d}x \, \mathrm{d}y</math>

该区域可以沿x或者y轴分解。要采用公式，必须先找到限制''D''的两个函数和定义区间。
这个例子中，这两个函数为：

:<math>\alpha (x) = x^2\,\!</math> 和 <math>\beta (x) = 1\,\!</math>

而区间为<math>[a,b] = [0,1]\,\!</math>（这里为了直观起见采用沿x轴分解）。

应用简化公式，得到：

:<math>\iint_D (x+y) \, \mathrm{d}x \, \mathrm{d}y = \int_0^1 \mathrm{d}x \int_{x^2}^1 (x+y) \, \mathrm{d}y = \int_0^1 \mathrm{d}x \ \left[xy \ + \ \frac{y^2}{2} \ \right]^1_{x^2}</math>

（首先，第二个积分将''x''作为常数）。然后就是用积分的基本技术：

:<math>\int_0^1 \left[xy \ + \ \frac{y^2}{2} \ \right]^1_{x^2} \, \mathrm{d}x = \int_0^1 \left(x + \frac{1}{2} - x^3 - \frac{x^4}{2} \right) \mathrm{d}x = \cdots = \frac{13}{20}</math>。

如果沿着y轴分解，可以计算

:<math>\int_0^1 \mathrm{d}y \int_0^{\sqrt{y}} (x+y) \, \mathrm{d}x</math>

并得到同样的结果。

[[File:Dominio-normalità r3 esempio.svg|thumb|160px|R3中可(沿''xy''平面）分解区域的例子。]]

==== R<sup>3</sup>中的分解 ====
这些公式可以推广到三重积分：

''T''是一个可以投影到''xy''平面的体，它夹在α (''x'',''y'')和β(''x'',''y'')两个函数之间。那么：

:<math>\iiint_T f(x,y,z) \ \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z = \iint_D \mathrm{d}x\, \mathrm{d}y \int_{\alpha (x,y)}^{\beta (x,y)} f(x,y,z) \, \mathrm{d}z</math>

（此定义和其它'''R'''<sup>3</sup>中的分解类似）。

=== 变量替换 ===

积分的极限常常不易交换（区域无法分解或者公式很复杂），这时可以采用''变量替换''来重写积分，令区域更加简易，从而可以用更简单的公式表达。为此，函数必须变换到新坐标系下。

:<u>例（1-a）</u>:
::函数为<math>f(x, y) = (x-1)^2 +\sqrt y</math>;
::若采用替换<math>x' = x-1, \ y'= y \, \!</math>则<math>x = x' + 1, \ y=y' \,\!</math>
::可以得到新函数<math>f_2(x,y) = (x')^2 +\sqrt y</math>.
* 对于定义域要进行类似处理，因为原来是采用变换前的变量表达的（本例中的''x''和''y''）。
* 微分''dx''和''dy''要通过包含被替换的变量对于新变量的偏微分的[[雅可比行列式]]来变换。（譬如，极坐标的微分变换）。

常用的变量替换有三种（'''R'''<sup>2</sup>中一种，'''R'''<sup>3</sup>中两种）；但是，更普遍的变换可以用同样的原理来发现。

==== 极坐标 ====
{{see also|极坐标系}}
[[File:Passaggio in coordinate polari.svg|thumb|270px|从笛卡尔坐标系变换到极坐标系。]]
在'''R'''<sup>2</sup>中，若定义域有某种圆形对称性而函数也有某种特征，则可以采用极坐标变换（参看图中的例子），也就是说将点''P(x,y)''从笛卡尔坐标变换到相应的极坐标中。这使得定义域的形状改变，从而简化运算。

该变换的基本关系如下：

:<math>f(x,y) \rightarrow f(\rho \ \cos \phi,\rho \ \sin \phi )</math>。

<blockquote>
<u>例（2-a）</u>:<br />
:函数为<math>f(x,y) = x + y\,\!</math>

:应用该变换得到

::<math>f(\rho, \phi) = \rho \cos \phi + \rho \sin \phi = \rho \ (\cos \phi + \sin \phi )</math>。

<u>例（2-b）</u>:
:函数为<math>f(x,y) = x^2 + y^2\,\!</math>
:这裡有：

::<math>f(\rho, \phi) = \rho^2 (\cos^2 \phi + \sin^2 \phi) = \rho^2\,\!</math>

:这里使用了[[勾股定理]]（在简化操作时很有用）。</blockquote>

定义域的变换是根据''x''和''y''通过环厚和角度的幅度来限定ρ, φ的区间。
[[File:Esempio trasformazione dominio da cartesiano polare.svg|thumb|230px|从笛卡尔到极坐标的区域变换。]]

<blockquote>
<u>例（2-c）</u>:
:区域为<math>D = x^2 + y^2 \le 4\,\!</math>，圆周半径2；很明显，这个区域所覆盖的角度是整个圆周角，所以φ从0变化到2π，而环半径从0变化到2（内环为0的环形就是圆）。

<u>例（2-d）</u>:
:区域为<math>D = \{ x^2 + y^2 \le 9, \ x^2 + y^2 \ge 4, \ y \ge 0 \}</math>，这是在正''y''半平面中的圆环（参看示意图）；注意φ表示平面角而ρ从2变化到3。因此变换出来的区域为[[矩形]]：

::<math>T = \{ 2 \le \rho \le 3, \ 0 \le \phi \le \pi \}</math>.

该变换的[[雅可比行列式]]为：

:<math>\frac{\partial (x,y)}{\partial (\rho, \phi)} =
\begin{vmatrix}
\cos \phi & - \rho \sin \phi \\
\sin \phi & \rho \cos \phi
\end{vmatrix} = \rho</math>

这可以通过将''x'' = ρ cos(φ), ''y'' = ρ sin(φ)代入关于ρ的第一行和关于φ的第二行的偏微分中得到，所以微分''dx dy''变换为ρ ''d''ρ ''d''φ.

一旦函数和区域的变换完成后，可以定义极坐标中的变量变换公式：

:<math>\iint_D f(x,y) \ \mathrm{d}x\, \mathrm{d}y = \iint_T f(\rho \cos \phi, \rho \sin \phi) \rho \, \mathrm{d} \rho\, \mathrm{d} \phi</math>。

注意φ在[0, 2π]区间中有效，而ρ测量长度，因此只能取非负值。

此外，应用变量变换公式的前提是，雅可比行列式的值在变换后的积分变量(如此例中的ρ和φ)组成的有界闭区域(如此例中φ和ρ构成的二维域)上恒不为零。但是在极坐标中当且仅当ρ为零时，才有雅可比行列式为零，故可证明该变量变换公式成立。



<u>例 (2-e)</u>:
:函数为<math>f(x,y) = x\,\!</math>区域和例2-d相同。
:从前面对''D''的分析，我们知道ρ的区间为[2,3]，而φ的为[0,π].函数变换为：

::<math>f(x,y) = x \longrightarrow f(\rho,\phi) = \rho \ \cos \phi</math>。

:最后，应用积分公式：

::<math>\iint_D x \, \mathrm{d}x\, \mathrm{d}y = \iint_T \rho \cos \phi \ \rho \, \mathrm{d}\rho\, \mathrm{d}\phi</math>。

:一旦区间给定，就可以得到

::<math>\int_0^{\pi} \int_2^3 \rho^2 \cos \phi \ \mathrm{d} \rho \ \mathrm{d} \phi = \int_0^{\pi} \cos \phi \ \mathrm{d} \phi \left[ \frac{\rho^3}{3} \right]_2^3 = \left[ \sin \phi \right]_0^{\pi} \ \left(9 - \frac{8}{3} \right) = 0</math>。
</blockquote>

==== 柱极坐标 ====
[[File:Cylindrical Coordinates.svg|thumb|190px|柱极坐标。]]
'''R'''<sup>3</sup>中，在有圆形底面的定义域上的积分可以通过变换到[[柱极坐标系]]来完成；函数的变换用如下的关系进行：
<center><math>f(x,y,z) \rightarrow f(\rho \cos \phi, \rho \sin \phi, z)</math> </center>

区域的变换可以从图形中得到，因为底面的形状可能不同，而高遵循初始区域的形状。

<blockquote>
<u>例（3-a）</u>: <br />
:区域为<math>D = \{ x^2 + y^2 \le 9, \ x^2 + y^2 \ge 4, \ 0 \le z \le 5 \}</math>（也即底面为例2-d中的圆环的高度为5的"管道"）；如果采用变换，可以得到区域<math>T = \{ 2 \le \rho \le 3, \ 0 \le \phi \le \pi, \ 0 \le z \le 5 \}</math>（这是一个底面为例2-d中的矩形而高为5的长方体）。

因为''z''分量没有变化，''dx dy dz''和在极坐标中一样变化：变为''ρ dρ dφ dz''。

最后，变换到柱极坐标的最后公式为：

:<math>\iiint_D f(x,y,z) \, \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z = \iiint_T f(\rho \cos \phi, \rho \sin \phi, z) \rho \, \mathrm{d}\rho\, \mathrm{d}\phi\, \mathrm{d}z</math>。

这个方法在柱形或者锥形区域的情况较为适用，也适用于容易分辨''z''区间和变换圆形底面和函数的其它情况。

<u>例（3-b）</u>:
函数为<math>f(x,y,z) = x^2 + y^2 + z\,\!</math>而积分区域为[[圆柱]]: <math>D = \{ x^2 + y^2 \le 9, \ -5 \le z \le 5 \}</math>.
将''D''变换到柱极坐标如下：

:<math>T = \{ 0 \le \rho \le 3, \ 0 \le \phi \le 2 \pi, \ -5 \le z \le 5 \}</math>。

函数变为

:<math>f(\rho \cos \phi, \rho \sin \phi, z) = \rho^2 + z</math>

最有应用积分公式：

:<math>\iiint_D (x^2 + y^2 +z) \, \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z = \iiint_T ( \rho^2 + z) \rho \, \mathrm{d}\rho\, \mathrm{d}\phi\, \mathrm{d}z;</math>

推演一下公式，得到

:<math>\int_{-5}^5 \mathrm{d}z \int_0^{2 \pi} \mathrm{d}\phi \int_0^3 ( \rho^3 + \rho z )\, \mathrm{d}\rho = 2 \pi \int_{-5}^5 \left[ \frac{\rho^4}{4} + \frac{\rho^2 z}{2} \right]_0^3 \, \mathrm{d}z = 2 \pi \int_{-5}^5 \left( \frac{81}{4} + \frac{9}{2} z\right)\, \mathrm{d}z = \cdots = 405 \pi.</math></blockquote>

==== 球极坐标 ====
[[File:Spherical Coordinates (Colatitude, Longitude).svg|thumb|190px|球极坐标。【注意某些地区（如北美）角度标识相反】]]
'''R'''<sup>3</sup>中，有些区域有球形对称性，所以将积分区域的每点用两个角度和一个距离标识较为合适。因此可以采用变换到[[球极坐标系]]；函数变换由如下关系产生：
<center> <math>f(x,y,z) \longrightarrow f(\rho \cos \theta \sin \phi, \rho \sin \theta \sin \phi, \rho \cos \phi)\,\!</math> </center>

注意''z''轴上的点没有唯一表示，<math>\theta </math>可以在0到2π间变化。

这个方法最为适用的区域显然是球。

<blockquote>
<u>例（4-a）</u>: <br />
:区域为<math>D = x^2 + y^2 + z^2 \le 16</math>（球心在原点半径为4的球）；应用变换后得到：<math>T = \{ 0 \le \rho \le 4, \ 0 \le \phi \le \pi, \ 0 \le \theta \le 2 \pi \}</math>。

:坐标变换的雅可比行列式为：

::<math>\frac{\partial (x,y,z)}{\partial (\rho, \theta, \phi)} =
\begin{vmatrix}
\cos \theta \sin \phi & - \rho \sin \theta \sin \phi & \rho \cos \theta \cos \phi \\
\sin \theta \sin \phi & \rho \cos \theta \sin \phi & \rho \sin \theta \cos \phi \\
\cos \phi & 0 & - \rho \sin \phi
\end{vmatrix} = -\rho^2 \sin \phi</math>

:因此<math> \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z </math>变换为<math> \rho^2 \sin \phi \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi</math>.

:得到最后公式:

::<math>\iiint_D f(x,y,z) \, \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z = \iiint_T f(\rho \sin \theta \cos \phi, \rho \sin \theta \sin \phi, \rho \cos \theta) \rho^2 \sin \phi \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi</math>。</blockquote>

应当在积分区域为球形对称'''并且'''函数很容易通过基本三角公式简化的时候才使用这个方法。（参看例4-b）；其它情况下，可能使用柱极坐标更为合适（参看例4-c）。

:<math>\iiint_T f(a,b,c) \rho^2 \sin \phi \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi</math>。

注意从雅可比行列式来的<math>\rho^2</math>和<math>\sin \phi</math>因子。

注意下面例子中，φ和θ的作用反过来了。

<blockquote>
<u>例（4-b）</u>:<br />
:''D''和例4-a相同，而<math>f(x,y,z) = x^2 + y^2 + z^2\,\!</math>是被积函数。

:很容易变换为：

::<math>f(\rho \sin \theta \cos \phi, \rho \sin \theta \sin \phi, \rho \cos \theta) = \rho^2,</math>，

:而从''D''到''T''的变换是已知的：

::<math>(0 \le \rho \le 4, \ 0 \le \phi \le 2 \pi, \ 0 \le \theta \le \pi)</math>。

:应用积分公式：

::<math>\iiint_D (x^2 + y^2 +z^2) \, \mathrm{d}x\, \mathrm{d}y\, \mathrm{d}z = \iiint_T \rho^2 \ \rho^2 \sin \theta \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi,</math>

:并展开：

::<math>\iiint_T \rho^4 \sin \theta \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi = \int_0^{\pi} \sin \theta \,\mathrm{d}\theta \int_0^4 \rho^4 \mathrm{d} \rho \int_0^{2 \pi} \mathrm{d}\phi = 2 \pi \int_0^{\pi} \sin \theta \left[ \frac{\rho^5}{5} \right]_0^4 \, \mathrm{d} \theta</math>

::<math>= 2 \pi \left[ \frac{\rho^5}{5} \right]_0^4 \left[- \cos \theta \right]_0^{\pi} = 4 \pi \cdot \frac{1024}{5} = \frac{4096 \pi}{5}</math>。

<u>例（4-c）</u>:
:区域''D''是球心在原点半径为''3a''的球（<math>D = x^2 + y^2 + z^2 \le 9a^2 \,\!</math>）而<math>f(x,y,z) = x^2 + y^2\,\!</math>是被积函数。

:看起来采用球极坐标变换较为合适，但是事实上，限定新区域''T''的变量很明显应该是：

::<math>0 \le \rho \le 3a, \ 0 \le \phi \le 2 \pi, \ 0 \le \theta \le \pi</math>。

:但是采用这个变换就有

::<math>f(x,y,z) = x^2 + y^2 \longrightarrow \rho^2 \sin^2 \theta \cos^2 \phi + \rho^2 \sin^2 \theta \sin^2 \phi = \rho^2 \sin^2 \theta</math>.

:应用积分公式得到：

::<math>\iiint_T \rho^2 \sin^2 \theta \rho^2 \sin \theta \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi = \iiint_T \rho^4 \sin^3 \theta \, \mathrm{d}\rho\, \mathrm{d}\theta\, \mathrm{d}\phi</math>

:这很难求解。而如果采用柱极坐标，新的''T''区间为：

::<math>0 \le \rho \le 3a, \ 0 \le \phi \le 2 \pi, \ - \sqrt{9a^2 - \rho^2} \le z \le \sqrt{9a^2 - \rho^2};</math>

:''z''区间可以通过将球切成两个半球并求解从''D''的公式来的[[不等式]]得到（然后直接变换''x<sup>2</sup> + y<sup>2</sup>''到''ρ<sup>2</sup>''）。新函数就是''ρ<sup>2</sup>''.采用积分公式

::<math>\iiint_T \rho^2 \rho \ \mathrm{d} \rho \mathrm{d} \phi \mathrm{d}z</math>.

:得到

::<math>\int_0^{2 \pi} \mathrm{d}\phi \int_0^{3a} \rho^3 \mathrm{d}\rho \int_{- \sqrt{9a^2 - \rho^2} }^{\sqrt{9 a^2 - \rho^2} }\, \mathrm{d}z = 2 \pi \int_0^{3a} 2 \rho^3 \sqrt{9 a^2 - \rho^2} \, \mathrm{d}\rho</math>。

:然后应用变换

::<math>9 a^2 - \rho^2 = t\,\! \longrightarrow \mathrm{d}t = -2 \rho\, \mathrm{d}\rho \longrightarrow \mathrm{d}\rho = \frac{\mathrm{d} t}{- 2 \rho}\,\!</math>

:(新区间变为<math>0, 3a \longrightarrow 9 a^2, 0</math>)。得到

::<math>- 2 \pi \int_{9 a^2}^{0} \rho^2 \sqrt{t}\, \mathrm{d}t</math>

:因为<math>\rho^2 = 9 a^2 - t\,\!</math>，所以

::<math>-2 \pi \int_{9 a^2}^0 (9 a^2 - t) \sqrt{t}\, \mathrm{d}t,</math>

:将积分限反过来，然后分配括号中的项，很容易将积分分解为可以直接积分的两部分：

::<math>2 \pi \left[ \int_0^{9 a^2} 9 a^2 \sqrt{t} \, \mathrm{d}t - \int_0^{9 a^2} t \sqrt{t} \, \mathrm{d}t\right] = 2 \pi \left[9 a^2 \frac{2}{3} t^{ \frac{3}{2} } - \frac{2}{5} t^{ \frac{5}{2}} \right]_0^{9 a^2}</math>

::<math>= 2 \cdot 27 \pi a^5 \left ( 6 - \frac{18}{5} \right ) = \frac{648 \pi}{5} a^5.</math>

:由于采用柱极坐标，很容易就将这个三重积分变换为简单的单变量积分。</blockquote>

参看[[柱极和球极坐标下的∇]]中讨论的不同的体积元。

== 数学应用范例-体积计算 ==

利用上面描述的方法，很容易计算一些立体的体积。
* '''[[圆柱]]'''：半径为''R''的圆形底面作为定义域，将等于高度''h''的常函数作为积分对象。可以在极坐标中将体积写作：
::体积<math> = \int_0^{2 \pi } \mathrm{d} \phi \int_0^R h \rho \ \mathrm{d} \rho = h 2 \pi \left[\frac{\rho^2}{2 }\right]_0^R = \pi R^2 h</math>

: ''验证''：体积=底面积×高 = <math>\pi R^2 \cdot h</math>
* '''[[球]]''':可以作为常函数''1''在球极坐标下的半径为''R''的球中积分：

::体积<math> = \int_0^{2 \pi }\, \mathrm{d} \phi \int_0^{ \pi } \sin \theta\, \mathrm{d} \theta \int_0^R \rho^2\, \mathrm{d} \rho = 2 \pi \int_0^{ \pi } \sin \theta \frac{R^3}{3 }\, \mathrm{d} \theta = \frac{2}{3 } \pi R^3 [- \cos \theta]_0^{ \pi } = \frac{4}{3 } \pi R^3</math>。
* '''[[四面体]]''' (三[[棱锥]]或者说3维[[单纯形]])：顶点在原点，三条长度为''l''的边沿着各个笛卡尔坐标系轴向的四面体的体积可以通过简化公式计算，因为''xy''平面和'z'轴互相垂直，''x''和''y''垂直，被积函数是常数''1''。

::体积<math> = \int_0^\ell \mathrm{d}x \int_0^{\ell-x }\, \mathrm{d}y \int_0^{\ell-x-y }\, \mathrm{d}z = \int_0^\ell \mathrm{d}x \int_0^{\ell-x } (\ell - x - y)\, \mathrm{d}y</math>

:: <math>= \int_0^\ell (\ell^2 - 2\ell x + x^2 - \frac{ (\ell-x)^2 }{2 })\, \mathrm{d}x = \ell^3 - \ell \ell^2 + \frac{\ell^3}{3 } - \left[\frac{\ell^2}{2 } - \ell x + \frac{x^2}{2 }\right]_0^\ell = </math>

:: <math> = \frac{\ell^3}{3 } - \frac{\ell^3}{6 } = \frac{\ell^3}{6}</math>

: ''验证''：体积 = 底面积×高/3 = <math>\frac{\ell^2}{2 } \cdot \ell/3 = \frac{\ell^3}{6}</math>。

[[File:Dominio improprio.svg|thumb|140px|反常区间的例子。]]

== 多重广义积分 ==
定义域无界或者函数值在靠近定义域边界时无界的情况下，可以引入二重[[广义积分]]或者'''三重广义积分'''。

== 多重积分和累次积分 ==
{{参看|积分次序}}
[[富比尼定理]]断言若
:<math>\int_{A\times B} |f(x,y)|\,\mathrm{d}(x,y)<\infty,</math>
也即，积分绝对收敛，则多重积分和[[累次积分]]给出同样的结果，
:<math>\int_{A\times B} f(x,y)\,\mathrm{d}(x,y)=\int_A\left(\int_B f(x,y)\,dy\right)\,dx=\int_B\left(\int_A f(x,y)\,\mathrm{d}x\right)\,\mathrm{d}y</math>。
一个特例是如果<math>|f(x,y)|</math>是[[有界函数]]而''A''和''B''为[[有界集]]时。

如果积分不是绝对收敛，必须小心，不要混淆''多重积分''和''累次积分''的概念，特别是当它们采用形式上相同的记法的时候。记法

:<math>\int_0^1\int_0^1 f(x,y)\,\mathrm{d}y\,\mathrm{d}x</math>

在某些情况下表示累次积分而非真正的双重积分。累次积分中，外围的积分

:<math>\int_0^1 \cdots \, \mathrm{d}x</math>

是对于如下''x''的函数关于''x''的积分

:<math>g(x)=\int_0^1 f(x,y)\,\mathrm{d}y</math>。

双重积分却是定义在''xy''平面的区域上。若双重积分存在，则它等于两个累次积分中的任何一个（或者"<math>\mathrm{d}y\,\mathrm{d}x</math>"或者"<math>\mathrm{d}x\,\mathrm{d}y</math>"），它也就是通过其中之一来计算的。但是有时这两个累次积分存在，而双重积分不存在。这种情况下，有时两个累次积分不相等，也即，

:<math>\int_0^1\int_0^1 f(x,y)\,\mathrm{d}y\,\mathrm{d}x \neq \int_0^1\int_0^1 f(x,y)\,\mathrm{d}x\,\mathrm{d}y</math>。

这是[[条件收敛]]的积分的重排序的一个例子。

如果要强调使用双重积分而非累次积分时，可以采用如下记法
:<math>\int_{[0,1]\times[0,1]} f(x,y)\,\mathrm{d}x\,\mathrm{d}y</math>

== 一些实际应用 ==

很普遍地，像单变量一样，我们通过多重积分可以找到给定集合上的函数的平均值。给定一个集合''D'' ⊆ '''R'''<sup>''n''</sup>和一个在''D''上可积的函数''f''，''f''在区域上的平均值是

:<math>\bar{f} = \frac{1}{m(D)} \int_D f(x)\, \mathrm{d}x,</math>

其中''m''(''D'')是''D''的测度。

此外，这些积分在[[物理]]中有大量应用。

[[力学]]中，[[转动惯量]]可以作为[[密度]]乘以刚体和转轴的距离的平方的体积分（三重积分）计算：

:<math>I_z = \iiint_V \rho r^2\, \mathrm{d}V</math>。

与三维[[欧氏空间]]'''R'''<sup>3</sup>中的质量[[博雷尔测度|测度]]<math>\mathrm{d}m</math>表示的{{link-en|质量分布|Mass distribution}}关联的[[引力势]]：
:<math>V(\mathbf{x}) = -\int_{\mathbf{R}^3} \frac{G}{|\mathbf{x} - \mathbf{y}|}\,\mathrm{d}m(\mathbf{y}).</math>

如果存在一个连续函数 <math>\rho(\mathbf {x}) </math> 表示 '''x''' 处的密度分布, 那么 <math>\mathrm{d}m(\mathbf {x})=\rho(\mathbf {x}) \, \mathrm{d}^3\mathbf {x}</math>, 其中'''<math>\mathrm{d}^3\mathbf {x}</math>''' 是欧几里得[[体积元]], 那么引力势就是
:<math>V(\mathbf{x}) = -\int_{\mathbf{R}^3} \frac{G}{|\mathbf{x}-\mathbf{y}|}\,\rho(\mathbf{y})\,\mathrm{d}^3\mathbf{y}.</math>

在[[电磁学]]中，[[麦克斯韦方程组]]可以写作多重积分，用以计算总磁场和电场。下例中，由[[电荷分布]]产生的[[电场]]通过向量函数的''三重积分''得到：

:<math>\vec E = \frac {1}{4 \pi \epsilon_0} \iiint \frac {\vec r - \vec r'}{\left \| \vec r - \vec r' \right \|^3} \rho (\vec r')\, \operatorname{d}^3 r'</math>。

== 参看 ==
* 和多重积分相关的[[实分析]]定理：
** [[散度定理]]
** [[斯托克斯定理]]
** [[格林定理]]
*[[重复积分的柯西公式]]

== 参考 ==
* Robert A. Adams - Calculus: A Complete Course (5th Edition) ISBN 0-201-79131-5.

== 外部链接 ==
* {{MathWorld | urlname=MultipleIntegral | title=多重积分}}
* {{Springer | urlname=M/m065370 | title=Multiple integral | author=L.D. Kudryavtsev}}
* [http://user.mendelu.cz/marik/maw/index.php?lang=en&form=integral2 Mathematical Assistant on Web] {{Wayback|url=http://user.mendelu.cz/marik/maw/index.php?lang=en&form=integral2 |date=20140407155645 }} [[笛卡尔坐标系]]下和[[极坐标]]下的双重积分的在线计算 （包括求解过程中的中间步骤，由[[Maxima]](軟體)支持)

[[Category:积分学]]
[[Category:多变量微积分]]

[[de:Integralrechnung#Mehrdimensionale Integration]]