正如单变量的正函数的定积分代表函数图像和x 轴之间区域的面积 一样,正的双变量函数的双重积分 代表函数所定义的曲面和包含函数定义域 的平面之间所夹的区域的体积 。(注意同样的体积也可以通过三变量常函数f (x , y , z ) = 1在上述曲面和平面之间的区域中的三重积分 得到。若有更多变量,则多元函数 的多重积分给出超体积 。
n 元函数f (x 1 , x 2 ,…, x n )在定义域D 上的多重积分通常用嵌套的积分号按照演算的逆序标识(最左边的积分号最后计算),后面跟着被积函数 和正常次序的积分变量 (最右边的变量最后使用)。积分域或者对每个积分变量在每个积分号下标识,或者用一个变量标在最右边的积分号下:
∫ … ∫ D f ( x 1 , x 2 , … , x n ) d x 1 … d x n {\displaystyle \int \ldots \int _{\mathbf {D} }\;f(x_{1},x_{2},\ldots ,x_{n})\;\mathrm {d} x_{1}\!\ldots \mathrm {d} x_{n}} 因为不可能计算多于一个自变量的函数 的不定积分 ,“不定”多重积分是不存在的。因此所有多重积分都是“定”积分。
通常在坐标系中,多重积分都利用嵌套的累次积分 计算。而累次积分为了简便可记为:
∫ φ 1 ψ 1 d x 1 ∫ φ 2 ( x 1 ) ψ 2 ( x 1 ) d x 2 … ∫ φ n ( x 1 , x 2 , … , x n − 1 ) ψ n ( x 1 , x 2 , … , x n − 1 ) f ( x 1 , x 2 , … , x n ) d x n {\displaystyle \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}} 其中积分域为:
D = { ( x 1 , x 2 , … , x n ) | φ 1 ≤ x 1 ≤ ψ 1 , φ 2 ( x 1 ) ≤ x 2 ≤ ψ 2 ( x 1 ) , … , φ n ( x 1 , x 2 , … , x n − 1 ) ≤ x n ≤ ψ n ( x 1 , x 2 , … , x n − 1 ) } {\displaystyle D=\left\{(x_{1},x_{2},\dots ,x_{n})|\varphi _{1}\leq x_{1}\leq \psi _{1},\varphi _{2}(x_{1})\leq x_{2}\leq \psi _{2}(x_{1}),\dots ,\varphi _{n}(x_{1},x_{2},\dots ,x_{n-1})\leq x_{n}\leq \psi _{n}(x_{1},x_{2},\dots ,x_{n-1})\right\}} 注意的是,该式一般情况下并不表示多个定积分的积,在实际计算中从最右侧积分变量开始积分,其结果会作为外一层积分的被积函数。
譬如,边长为4 × 6 × 5的长方体 的体积可以通过两种方法得到:
通过函数f (x , y ) = 5在xy 平面中的区域D ,也就是长方体的底上的双重积分 ∬ D 5 d x d y {\displaystyle \iint _{\mathrm {D} }5\ \mathrm {d} x\,\mathrm {d} y} ∭ 1 d x d y d z {\displaystyle \iiint 1\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z} (1可以省略不写) 数学定义
编辑
令n 为大于1的自然数。考虑所谓的半开 n维矩形 (下面简称矩形 )。对于平面 ,n = 2 {\displaystyle n=2} 。
T = [ a 1 , b 1 ) × [ a 2 , b 2 ) × ⋯ × [ a n , b n ) ⊂ R n {\displaystyle T=[a_{1},b_{1})\times [a_{2},b_{2})\times \cdots \times [a_{n},b_{n})\subset \mathbb {R} ^{n}} 。将每个区间[a i , b i )分成有限个不重叠的子区间,每个都是左闭右开。将子区间记为I i 。则,所有所有如下形式的子矩形的族
C = I 1 × I 2 × ⋯ × I n {\displaystyle C=I_{1}\times I_{2}\times \cdots \times I_{n}} 是T 的一个划分 ,也即,子矩形C 是互不重叠的,而且它们聯集为T 。C 中的子矩形的直径 按照定义是C 中最大的边长,而T 的划分的直径定義成划分中的所有子矩形的最大直径。
令f : T → R 为定义在T 上的函数。考虑如上定义的T 的划分
T = C 1 ∪ C 2 ∪ ⋯ ∪ C m {\displaystyle T=C_{1}\cup C_{2}\cup \cdots \cup C_{m}} 其中m 是正整数。如下形式的和称为黎曼和
∑ k = 1 m f ( P k ) m ( C k ) {\displaystyle \sum _{k=1}^{m}f(P_{k})\,\operatorname {m} (C_{k})} 其中,对于每个k ,点P k 在C k 中,而m(C k )是笛卡尔积为C k 的区间的边长之积。
函数f 称为黎曼可积 ,如果如下极限 存在
S = lim δ → 0 ∑ k = 1 m f ( P k ) m ( C k ) {\displaystyle S=\lim _{\delta \to 0}\sum _{k=1}^{m}f(P_{k})\,\operatorname {m} \,(C_{k})} 其中极限取遍所有直径小於δ的T 的划分。若f 黎曼可积 ,S 称为f 在T 上的黎曼积分 。记为
∫ T f ( x ) d x {\displaystyle \int _{T}\!f(x)\,\mathrm {d} x} 。定义在任意有界n 维集合上的函数的黎曼和可以通过将函数延拓到一个半开半闭矩形上来求出,其取值在原来的定义域之外为0。然后,原来的函数的积分就定义为延展的函数在矩形区域中的积分(如果存在的话)。
下文中n 维黎曼积分简称多重积分 。
多重积分具有很多与单变量函数的积分一样的性质(线性,可加性,单调性,等等)。而且,和单变量情况一样,可以用多重积分找出函数在给定集合上的积分。具体来讲,给定集合D ⊆ R n 和D 上的可积函数f ,f 在定义域上的平均值为
f ¯ = 1 m ( D ) ∫ D f ( x ) d x , {\displaystyle {\bar {f}}={\frac {1}{m(D)}}\int _{D}f(x)\,\mathrm {d} x,} 其中m (D )是D 的测度 。
T ⊆ R 2 时,积分
ℓ = ∬ T f ( x , y ) d x d y {\displaystyle \ell =\iint _{T}f(x,y)\,\mathrm {d} x\,\mathrm {d} y} 是f 在T 上的双重积分 ,而若T ⊆ R 3 ,积分
ℓ = ∭ T f ( x , y , z ) d x d y d z {\displaystyle \ell =\iiint _{T}f(x,y,z)\,\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z} 是f 在T 上的三重积分 。
注意,按常规,双重积分用两个积分号,而三重积分有三个;这只是记法上方便,也是为了通过重复积分来计算多重积分(参看本条目后文)。
积分方法
编辑
多重积分问题的解决在多数情况下依赖于将多重积分转化为一系列单变量积分,而其中每个单变量积分都是直接可解的。
直接检验
编辑
有时可以直接获得积分的结果,而无需任何直接计算。
在常函数 的情况中,结果很直接:只要将常函数c 乘以测度就可以了。如果c = 1,而且是在R 2 的子集中积分,则乘积就是区域面积,而在R 3 中,它就是区域的体积。
D = { ( x , y ) ∈ R 2 : 2 ≤ x ≤ 4 ; 3 ≤ y ≤ 6 } {\displaystyle D=\{(x,y)\in \mathbb {R} ^{2}\ :\ 2\leq x\leq 4\ ;\ 3\leq y\leq 6\}} and f ( x , y ) = 2 {\displaystyle f(x,y)=2\,\!} 在D 上积分f : ∫ 3 6 ∫ 2 4 2 d x d y = area ( D ) ⋅ 2 = ( 2 ⋅ 3 ) ⋅ 2 = 12 {\displaystyle \int _{3}^{6}\int _{2}^{4}\ 2\ \mathrm {d} x\,\mathrm {d} y={\mbox{area}}(D)\cdot 2=(2\cdot 3)\cdot 2=12} 。利用可能的对称性
编辑
如果定义域存在沿着某条轴的对称性而且函数对于那个变量是奇函数,则积分为0(因为相反的两部分加起来为0)。
对于R n 中的函数,只要相关变量对于形成对称的轴是奇变量就可以了。
例一:
给定f ( x , y ) = 2 sin ( x ) − 3 y 3 + 5 {\displaystyle f(x,y)=2\sin(x)-3y^{3}+5} 以及T = { ( x , y ) ∈ R 2 : x 2 + y 2 ≤ 1 } {\displaystyle T=\left\{(x,y)\in \mathbf {R} ^{2}\ :\ x^{2}+y^{2}\leq 1\right\}} 为积分区域(半径 为1的圆盘,包含边界 )。 利用线性性质,积分可以分解为三部分: ∬ T ( 2 sin x − 3 y 3 + 5 ) d x d y = ∬ T 2 sin x d x d y − ∬ T 3 y 3 d x d y + ∬ T 5 d x d y {\displaystyle \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} 2 sin ( x ) {\textstyle 2\sin(x)} 和 3 y 3 {\textstyle 3y^{3}} 都是奇函数,而且显然T 对于x 和y 轴都是对称的;因此唯一有贡献的部分是常函数5 因为其它两个都贡献0.
例二:
考虑函数f ( x , y , z ) = x exp ( y 2 + z 2 ) {\textstyle f(x,y,z)=x\exp(y^{2}+z^{2})} 以及圆心在原点的半径为2的球
T = { ( x , y , z ) ∈ R 3 : x 2 + y 2 + z 2 ≤ 4 } . {\displaystyle T=\left\{(x,y,z)\in \mathbf {R} ^{3}\ :\ x^{2}+y^{2}+z^{2}\leq 4\right\}.}
该球显然是对于三条轴都对称,但是只要对于x 轴积分就可以看出结果是0,因为f 对于该变量是奇函数。
简化公式
编辑
简化公式基于简单积分区域来将多重积分转化为单变量积分的序列。它们必须从右至左计算,过程中将其它变量暂时视为常数(和偏导数 的计算类似)。
R2 中的常规区域
编辑
此種方法適用於滿足下述條件的任何定義域 D :
D 投影到 x 軸或 y 軸任一軸,形成一個有邊界的範圍, 以 a , b 代表邊界值。
通過 a , b 兩點並與a b ¯ {\displaystyle {\overline {ab}}} 垂直的直線與 D 相交後的兩個端點,可以用 2 個函數α {\displaystyle \alpha } , β {\displaystyle \beta } 定義。
將 D 對 x 軸做垂直投影,函數f : D ⟶ R {\displaystyle f:D\longrightarrow \mathbb {R} } 是连续函数,并且D 可以视为(定义在[a ,b ]区间上的)α(x )和β(x )之间的区域。则
∬ D f ( x , y ) d x d y = ∫ a b d x ∫ α ( x ) β ( x ) f ( x , y ) d y {\displaystyle \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}
將 D 對y 軸做垂直投影,函數f : D ⟶ R {\displaystyle f:D\longrightarrow \mathbb {R} } 是连续函数,并且D 可以视为(定义在[a ,b ]区间上的)α(y )和β(y )之间的区域。则
∬ D f ( x , y ) d x d y = ∫ a b d y ∫ α ( y ) β ( y ) f ( x , y ) d x {\displaystyle \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}
例:可以采用简化公式的D区域。 考虑区域:D = { ( x , y ) : x ≥ 0 , y ≤ 1 , y ≥ x 2 } {\displaystyle D=\{(x,y)\ :\ x\geq 0,y\leq 1,y\geq x^{2}\}} (参看附图)。计算
∬ D ( x + y ) d x d y {\displaystyle \iint _{D}(x+y)\,\mathrm {d} x\,\mathrm {d} y} 该区域可以沿x或者y轴分解。要采用公式,必须先找到限制D 的两个函数和定义区间。
这个例子中,这两个函数为:
α ( x ) = x 2 {\displaystyle \alpha (x)=x^{2}\,\!} 和 β ( x ) = 1 {\displaystyle \beta (x)=1\,\!} 而区间为[ a , b ] = [ 0 , 1 ] {\displaystyle [a,b]=[0,1]\,\!} (这里为了直观起见采用沿x轴分解)。
应用简化公式,得到:
∬ D ( x + y ) d x d y = ∫ 0 1 d x ∫ x 2 1 ( x + y ) d y = ∫ 0 1 d x [ x y + y 2 2 ] x 2 1 {\displaystyle \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]_{x^{2}}^{1}} (首先,第二个积分将x 作为常数)。然后就是用积分的基本技术:
∫ 0 1 [ x y + y 2 2 ] x 2 1 d x = ∫ 0 1 ( x + 1 2 − x 3 − x 4 2 ) d x = ⋯ = 13 20 {\displaystyle \int _{0}^{1}\left[xy\ +\ {\frac {y^{2}}{2}}\ \right]_{x^{2}}^{1}\,\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}}} 。如果沿着y轴分解,可以计算
∫ 0 1 d y ∫ 0 y ( x + y ) d x {\displaystyle \int _{0}^{1}\mathrm {d} y\int _{0}^{\sqrt {y}}(x+y)\,\mathrm {d} x} 并得到同样的结果。
R3中可(沿xy 平面)分解区域的例子。 R3 中的分解
编辑
这些公式可以推广到三重积分:
T 是一个可以投影到xy 平面的体,它夹在α (x ,y )和β(x ,y )两个函数之间。那么:
∭ T f ( x , y , z ) d x d y d z = ∬ D d x d y ∫ α ( x , y ) β ( x , y ) f ( x , y , z ) d z {\displaystyle \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} (此定义和其它R 3 中的分解类似)。
变量替换
编辑
积分的极限常常不易交换(区域无法分解或者公式很复杂),这时可以采用变量替换 来重写积分,令区域更加简易,从而可以用更简单的公式表达。为此,函数必须变换到新坐标系下。
例(1-a) :
函数为f ( x , y ) = ( x − 1 ) 2 + y {\displaystyle f(x,y)=(x-1)^{2}+{\sqrt {y}}} ;
若采用替换x ′ = x − 1 , y ′ = y {\displaystyle x'=x-1,\ y'=y\,\!} 则x = x ′ + 1 , y = y ′ {\displaystyle x=x'+1,\ y=y'\,\!}
可以得到新函数f 2 ( x , y ) = ( x ′ ) 2 + y {\displaystyle f_{2}(x,y)=(x')^{2}+{\sqrt {y}}} . 对于定义域要进行类似处理,因为原来是采用变换前的变量表达的(本例中的x 和y )。
微分dx 和dy 要通过包含被替换的变量对于新变量的偏微分的雅可比行列式 来变换。(譬如,极坐标的微分变换)。 常用的变量替换有三种(R 2 中一种,R 3 中两种);但是,更普遍的变换可以用同样的原理来发现。
从笛卡尔坐标系变换到极坐标系。 在R 2 中,若定义域有某种圆形对称性而函数也有某种特征,则可以采用极坐标变换(参看图中的例子),也就是说将点P(x,y) 从笛卡尔坐标变换到相应的极坐标中。这使得定义域的形状改变,从而简化运算。
该变换的基本关系如下:
f ( x , y ) → f ( ρ cos ϕ , ρ sin ϕ ) {\displaystyle f(x,y)\rightarrow f(\rho \ \cos \phi ,\rho \ \sin \phi )} 。
例(2-a) :
函数为f ( x , y ) = x + y {\displaystyle f(x,y)=x+y\,\!} 应用该变换得到 f ( ρ , ϕ ) = ρ cos ϕ + ρ sin ϕ = ρ ( cos ϕ + sin ϕ ) {\displaystyle f(\rho ,\phi )=\rho \cos \phi +\rho \sin \phi =\rho \ (\cos \phi +\sin \phi )} 。例(2-b) :
函数为f ( x , y ) = x 2 + y 2 {\displaystyle f(x,y)=x^{2}+y^{2}\,\!}
这裡有: f ( ρ , ϕ ) = ρ 2 ( cos 2 ϕ + sin 2 ϕ ) = ρ 2 {\displaystyle f(\rho ,\phi )=\rho ^{2}(\cos ^{2}\phi +\sin ^{2}\phi )=\rho ^{2}\,\!} 这里使用了勾股定理 (在简化操作时很有用)。
定义域的变换是根据x 和y 通过环厚和角度的幅度来限定ρ, φ的区间。
从笛卡尔到极坐标的区域变换。
例(2-c) :
区域为D = x 2 + y 2 ≤ 4 {\displaystyle D=x^{2}+y^{2}\leq 4\,\!} ,圆周半径2;很明显,这个区域所覆盖的角度是整个圆周角,所以φ从0变化到2π,而环半径从0变化到2(内环为0的环形就是圆)。 例(2-d) :
区域为D = { x 2 + y 2 ≤ 9 , x 2 + y 2 ≥ 4 , y ≥ 0 } {\displaystyle D=\{x^{2}+y^{2}\leq 9,\ x^{2}+y^{2}\geq 4,\ y\geq 0\}} ,这是在正y 半平面中的圆环(参看示意图);注意φ表示平面角而ρ从2变化到3。因此变换出来的区域为矩形 : T = { 2 ≤ ρ ≤ 3 , 0 ≤ ϕ ≤ π } {\displaystyle T=\{2\leq \rho \leq 3,\ 0\leq \phi \leq \pi \}} .该变换的雅可比行列式 为:
∂ ( x , y ) ∂ ( ρ , ϕ ) = | cos ϕ − ρ sin ϕ sin ϕ ρ cos ϕ | = ρ {\displaystyle {\frac {\partial (x,y)}{\partial (\rho ,\phi )}}={\begin{vmatrix}\cos \phi &-\rho \sin \phi \\\sin \phi &\rho \cos \phi \end{vmatrix}}=\rho } 这可以通过将x = ρ cos(φ), y = ρ sin(φ)代入关于ρ的第一行和关于φ的第二行的偏微分中得到,所以微分dx dy 变换为ρ d ρ d φ.
一旦函数和区域的变换完成后,可以定义极坐标中的变量变换公式:
∬ D f ( x , y ) d x d y = ∬ T f ( ρ cos ϕ , ρ sin ϕ ) ρ d ρ d ϕ {\displaystyle \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 } 。注意φ在[0, 2π]区间中有效,而ρ测量长度,因此只能取非负值。
此外,应用变量变换公式的前提是,雅可比行列式的值在变换后的积分变量(如此例中的ρ和φ)组成的有界闭区域(如此例中φ和ρ构成的二维域)上恒不为零。但是在极坐标中当且仅当ρ为零时,才有雅可比行列式为零,故可证明该变量变换公式成立。
例 (2-e) :
函数为f ( x , y ) = x {\displaystyle f(x,y)=x\,\!} 区域和例2-d相同。
从前面对D 的分析,我们知道ρ的区间为[2,3],而φ的为[0,π].函数变换为: f ( x , y ) = x ⟶ f ( ρ , ϕ ) = ρ cos ϕ {\displaystyle f(x,y)=x\longrightarrow f(\rho ,\phi )=\rho \ \cos \phi } 。最后,应用积分公式: ∬ D x d x d y = ∬ T ρ cos ϕ ρ d ρ d ϕ {\displaystyle \iint _{D}x\,\mathrm {d} x\,\mathrm {d} y=\iint _{T}\rho \cos \phi \ \rho \,\mathrm {d} \rho \,\mathrm {d} \phi } 。一旦区间给定,就可以得到 ∫ 0 π ∫ 2 3 ρ 2 cos ϕ d ρ d ϕ = ∫ 0 π cos ϕ d ϕ [ ρ 3 3 ] 2 3 = [ sin ϕ ] 0 π ( 9 − 8 3 ) = 0 {\displaystyle \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} 。
柱极坐标
编辑
柱极坐标。 R 3 中,在有圆形底面的定义域上的积分可以通过变换到柱极坐标系 来完成;函数的变换用如下的关系进行:
f ( x , y , z ) → f ( ρ cos ϕ , ρ sin ϕ , z ) {\displaystyle f(x,y,z)\rightarrow f(\rho \cos \phi ,\rho \sin \phi ,z)}
区域的变换可以从图形中得到,因为底面的形状可能不同,而高遵循初始区域的形状。
例(3-a) :
区域为D = { x 2 + y 2 ≤ 9 , x 2 + y 2 ≥ 4 , 0 ≤ z ≤ 5 } {\displaystyle D=\{x^{2}+y^{2}\leq 9,\ x^{2}+y^{2}\geq 4,\ 0\leq z\leq 5\}} (也即底面为例2-d中的圆环的高度为5的"管道");如果采用变换,可以得到区域T = { 2 ≤ ρ ≤ 3 , 0 ≤ ϕ ≤ π , 0 ≤ z ≤ 5 } {\displaystyle T=\{2\leq \rho \leq 3,\ 0\leq \phi \leq \pi ,\ 0\leq z\leq 5\}} (这是一个底面为例2-d中的矩形而高为5的长方体)。 因为z 分量没有变化,dx dy dz 和在极坐标中一样变化:变为ρ dρ dφ dz 。
最后,变换到柱极坐标的最后公式为:
∭ D f ( x , y , z ) d x d y d z = ∭ T f ( ρ cos ϕ , ρ sin ϕ , z ) ρ d ρ d ϕ d z {\displaystyle \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} 。这个方法在柱形或者锥形区域的情况较为适用,也适用于容易分辨z 区间和变换圆形底面和函数的其它情况。
例(3-b) :
函数为f ( x , y , z ) = x 2 + y 2 + z {\displaystyle f(x,y,z)=x^{2}+y^{2}+z\,\!} 而积分区域为圆柱 : D = { x 2 + y 2 ≤ 9 , − 5 ≤ z ≤ 5 } {\displaystyle D=\{x^{2}+y^{2}\leq 9,\ -5\leq z\leq 5\}} .
将D 变换到柱极坐标如下:
T = { 0 ≤ ρ ≤ 3 , 0 ≤ ϕ ≤ 2 π , − 5 ≤ z ≤ 5 } {\displaystyle T=\{0\leq \rho \leq 3,\ 0\leq \phi \leq 2\pi ,\ -5\leq z\leq 5\}} 。函数变为
f ( ρ cos ϕ , ρ sin ϕ , z ) = ρ 2 + z {\displaystyle f(\rho \cos \phi ,\rho \sin \phi ,z)=\rho ^{2}+z} 最有应用积分公式:
∭ D ( x 2 + y 2 + z ) d x d y d z = ∭ T ( ρ 2 + z ) ρ d ρ d ϕ d z ; {\displaystyle \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;} 推演一下公式,得到
∫ − 5 5 d z ∫ 0 2 π d ϕ ∫ 0 3 ( ρ 3 + ρ z ) d ρ = 2 π ∫ − 5 5 [ ρ 4 4 + ρ 2 z 2 ] 0 3 d z = 2 π ∫ − 5 5 ( 81 4 + 9 2 z ) d z = ⋯ = 405 π . {\displaystyle \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 .}
球极坐标
编辑
球极坐标。【注意某些地区(如北美)角度标识相反】 R 3 中,有些区域有球形对称性,所以将积分区域的每点用两个角度和一个距离标识较为合适。因此可以采用变换到球极坐标系 ;函数变换由如下关系产生:
f ( x , y , z ) ⟶ f ( ρ cos θ sin ϕ , ρ sin θ sin ϕ , ρ cos ϕ ) {\displaystyle f(x,y,z)\longrightarrow f(\rho \cos \theta \sin \phi ,\rho \sin \theta \sin \phi ,\rho \cos \phi )\,\!}
注意z 轴上的点没有唯一表示,θ {\displaystyle \theta } 可以在0到2π间变化。
这个方法最为适用的区域显然是球。
例(4-a) :
区域为D = x 2 + y 2 + z 2 ≤ 16 {\displaystyle D=x^{2}+y^{2}+z^{2}\leq 16} (球心在原点半径为4的球);应用变换后得到:T = { 0 ≤ ρ ≤ 4 , 0 ≤ ϕ ≤ π , 0 ≤ θ ≤ 2 π } {\displaystyle T=\{0\leq \rho \leq 4,\ 0\leq \phi \leq \pi ,\ 0\leq \theta \leq 2\pi \}} 。 坐标变换的雅可比行列式为: ∂ ( x , y , z ) ∂ ( ρ , θ , ϕ ) = | cos θ sin ϕ − ρ sin θ sin ϕ ρ cos θ cos ϕ sin θ sin ϕ ρ cos θ sin ϕ ρ sin θ cos ϕ cos ϕ 0 − ρ sin ϕ | = − ρ 2 sin ϕ {\displaystyle {\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 } 因此d x d y d z {\displaystyle \mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z} 变换为ρ 2 sin ϕ d ρ d θ d ϕ {\displaystyle \rho ^{2}\sin \phi \,\mathrm {d} \rho \,\mathrm {d} \theta \,\mathrm {d} \phi } . 得到最后公式: ∭ D f ( x , y , z ) d x d y d z = ∭ T f ( ρ sin θ cos ϕ , ρ sin θ sin ϕ , ρ cos θ ) ρ 2 sin ϕ d ρ d θ d ϕ {\displaystyle \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 } 。
应当在积分区域为球形对称并且 函数很容易通过基本三角公式简化的时候才使用这个方法。(参看例4-b);其它情况下,可能使用柱极坐标更为合适(参看例4-c)。
∭ T f ( a , b , c ) ρ 2 sin ϕ d ρ d θ d ϕ {\displaystyle \iiint _{T}f(a,b,c)\rho ^{2}\sin \phi \,\mathrm {d} \rho \,\mathrm {d} \theta \,\mathrm {d} \phi } 。注意从雅可比行列式来的ρ 2 {\displaystyle \rho ^{2}} 和sin ϕ {\displaystyle \sin \phi } 因子。
注意下面例子中,φ和θ的作用反过来了。
例(4-b) :
D 和例4-a相同,而f ( x , y , z ) = x 2 + y 2 + z 2 {\displaystyle f(x,y,z)=x^{2}+y^{2}+z^{2}\,\!} 是被积函数。很容易变换为: f ( ρ sin θ cos ϕ , ρ sin θ sin ϕ , ρ cos θ ) = ρ 2 , {\displaystyle f(\rho \sin \theta \cos \phi ,\rho \sin \theta \sin \phi ,\rho \cos \theta )=\rho ^{2},} ,而从D 到T 的变换是已知的: ( 0 ≤ ρ ≤ 4 , 0 ≤ ϕ ≤ 2 π , 0 ≤ θ ≤ π ) {\displaystyle (0\leq \rho \leq 4,\ 0\leq \phi \leq 2\pi ,\ 0\leq \theta \leq \pi )} 。应用积分公式: ∭ D ( x 2 + y 2 + z 2 ) d x d y d z = ∭ T ρ 2 ρ 2 sin θ d ρ d θ d ϕ , {\displaystyle \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 ,} 并展开: ∭ T ρ 4 sin θ d ρ d θ d ϕ = ∫ 0 π sin θ d θ ∫ 0 4 ρ 4 d ρ ∫ 0 2 π d ϕ = 2 π ∫ 0 π sin θ [ ρ 5 5 ] 0 4 d θ {\displaystyle \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 } = 2 π [ ρ 5 5 ] 0 4 [ − cos θ ] 0 π = 4 π ⋅ 1024 5 = 4096 π 5 {\displaystyle =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}}} 。例(4-c) :
区域D 是球心在原点半径为3a 的球(D = x 2 + y 2 + z 2 ≤ 9 a 2 {\displaystyle D=x^{2}+y^{2}+z^{2}\leq 9a^{2}\,\!} )而f ( x , y , z ) = x 2 + y 2 {\displaystyle f(x,y,z)=x^{2}+y^{2}\,\!} 是被积函数。 看起来采用球极坐标变换较为合适,但是事实上,限定新区域T 的变量很明显应该是: 0 ≤ ρ ≤ 3 a , 0 ≤ ϕ ≤ 2 π , 0 ≤ θ ≤ π {\displaystyle 0\leq \rho \leq 3a,\ 0\leq \phi \leq 2\pi ,\ 0\leq \theta \leq \pi } 。但是采用这个变换就有 f ( x , y , z ) = x 2 + y 2 ⟶ ρ 2 sin 2 θ cos 2 ϕ + ρ 2 sin 2 θ sin 2 ϕ = ρ 2 sin 2 θ {\displaystyle 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 } .应用积分公式得到: ∭ T ρ 2 sin 2 θ ρ 2 sin θ d ρ d θ d ϕ = ∭ T ρ 4 sin 3 θ d ρ d θ d ϕ {\displaystyle \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 } 这很难求解。而如果采用柱极坐标,新的T 区间为: 0 ≤ ρ ≤ 3 a , 0 ≤ ϕ ≤ 2 π , − 9 a 2 − ρ 2 ≤ z ≤ 9 a 2 − ρ 2 ; {\displaystyle 0\leq \rho \leq 3a,\ 0\leq \phi \leq 2\pi ,\ -{\sqrt {9a^{2}-\rho ^{2}}}\leq z\leq {\sqrt {9a^{2}-\rho ^{2}}};} z 区间可以通过将球切成两个半球并求解从D 的公式来的不等式 得到(然后直接变换x2 + y2 到ρ2 )。新函数就是ρ2 .采用积分公式∭ T ρ 2 ρ d ρ d ϕ d z {\displaystyle \iiint _{T}\rho ^{2}\rho \ \mathrm {d} \rho \mathrm {d} \phi \mathrm {d} z} .得到 ∫ 0 2 π d ϕ ∫ 0 3 a ρ 3 d ρ ∫ − 9 a 2 − ρ 2 9 a 2 − ρ 2 d z = 2 π ∫ 0 3 a 2 ρ 3 9 a 2 − ρ 2 d ρ {\displaystyle \int _{0}^{2\pi }\mathrm {d} \phi \int _{0}^{3a}\rho ^{3}\mathrm {d} \rho \int _{-{\sqrt {9a^{2}-\rho ^{2}}}}^{\sqrt {9a^{2}-\rho ^{2}}}\,\mathrm {d} z=2\pi \int _{0}^{3a}2\rho ^{3}{\sqrt {9a^{2}-\rho ^{2}}}\,\mathrm {d} \rho } 。然后应用变换 9 a 2 − ρ 2 = t ⟶ d t = − 2 ρ d ρ ⟶ d ρ = d t − 2 ρ {\displaystyle 9a^{2}-\rho ^{2}=t\,\!\longrightarrow \mathrm {d} t=-2\rho \,\mathrm {d} \rho \longrightarrow \mathrm {d} \rho ={\frac {\mathrm {d} t}{-2\rho }}\,\!} (新区间变为0 , 3 a ⟶ 9 a 2 , 0 {\displaystyle 0,3a\longrightarrow 9a^{2},0} )。得到 − 2 π ∫ 9 a 2 0 ρ 2 t d t {\displaystyle -2\pi \int _{9a^{2}}^{0}\rho ^{2}{\sqrt {t}}\,\mathrm {d} t} 因为ρ 2 = 9 a 2 − t {\displaystyle \rho ^{2}=9a^{2}-t\,\!} ,所以 − 2 π ∫ 9 a 2 0 ( 9 a 2 − t ) t d t , {\displaystyle -2\pi \int _{9a^{2}}^{0}(9a^{2}-t){\sqrt {t}}\,\mathrm {d} t,} 将积分限反过来,然后分配括号中的项,很容易将积分分解为可以直接积分的两部分: 2 π [ ∫ 0 9 a 2 9 a 2 t d t − ∫ 0 9 a 2 t t d t ] = 2 π [ 9 a 2 2 3 t 3 2 − 2 5 t 5 2 ] 0 9 a 2 {\displaystyle 2\pi \left[\int _{0}^{9a^{2}}9a^{2}{\sqrt {t}}\,\mathrm {d} t-\int _{0}^{9a^{2}}t{\sqrt {t}}\,\mathrm {d} t\right]=2\pi \left[9a^{2}{\frac {2}{3}}t^{\frac {3}{2}}-{\frac {2}{5}}t^{\frac {5}{2}}\right]_{0}^{9a^{2}}} = 2 ⋅ 27 π a 5 ( 6 − 18 5 ) = 648 π 5 a 5 . {\displaystyle =2\cdot 27\pi a^{5}\left(6-{\frac {18}{5}}\right)={\frac {648\pi }{5}}a^{5}.} 由于采用柱极坐标,很容易就将这个三重积分变换为简单的单变量积分。
参看柱极和球极坐标下的∇ 中讨论的不同的体积元。
数学应用范例-体积计算
编辑
多重广义积分
编辑
定义域无界或者函数值在靠近定义域边界时无界的情况下,可以引入二重广义积分 或者三重广义积分 。
多重积分和累次积分
编辑
富比尼定理 断言若
∫ A × B | f ( x , y ) | d ( x , y ) < ∞ , {\displaystyle \int _{A\times B}|f(x,y)|\,\mathrm {d} (x,y)<\infty ,} 也即,积分绝对收敛,则多重积分和累次积分 给出同样的结果,
∫ A × B f ( x , y ) d ( x , y ) = ∫ A ( ∫ B f ( x , y ) d y ) d x = ∫ B ( ∫ A f ( x , y ) d x ) d y {\displaystyle \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} 。一个特例是如果| f ( x , y ) | {\displaystyle |f(x,y)|} 是有界函数 而A 和B 为有界集 时。
如果积分不是绝对收敛,必须小心,不要混淆多重积分 和累次积分 的概念,特别是当它们采用形式上相同的记法的时候。记法
∫ 0 1 ∫ 0 1 f ( x , y ) d y d x {\displaystyle \int _{0}^{1}\int _{0}^{1}f(x,y)\,\mathrm {d} y\,\mathrm {d} x} 在某些情况下表示累次积分而非真正的双重积分。累次积分中,外围的积分
∫ 0 1 ⋯ d x {\displaystyle \int _{0}^{1}\cdots \,\mathrm {d} x} 是对于如下x 的函数关于x 的积分
g ( x ) = ∫ 0 1 f ( x , y ) d y {\displaystyle g(x)=\int _{0}^{1}f(x,y)\,\mathrm {d} y} 。双重积分却是定义在xy 平面的区域上。若双重积分存在,则它等于两个累次积分中的任何一个(或者"d y d x {\displaystyle \mathrm {d} y\,\mathrm {d} x} "或者"d x d y {\displaystyle \mathrm {d} x\,\mathrm {d} y} "),它也就是通过其中之一来计算的。但是有时这两个累次积分存在,而双重积分不存在。这种情况下,有时两个累次积分不相等,也即,
∫ 0 1 ∫ 0 1 f ( x , y ) d y d x ≠ ∫ 0 1 ∫ 0 1 f ( x , y ) d x d y {\displaystyle \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} 。这是条件收敛 的积分的重排序的一个例子。
如果要强调使用双重积分而非累次积分时,可以采用如下记法
∫ [ 0 , 1 ] × [ 0 , 1 ] f ( x , y ) d x d y {\displaystyle \int _{[0,1]\times [0,1]}f(x,y)\,\mathrm {d} x\,\mathrm {d} y} 一些实际应用
编辑
很普遍地,像单变量一样,我们通过多重积分可以找到给定集合上的函数的平均值。给定一个集合D ⊆ R n 和一个在D 上可积的函数f ,f 在区域上的平均值是
f ¯ = 1 m ( D ) ∫ D f ( x ) d x , {\displaystyle {\bar {f}}={\frac {1}{m(D)}}\int _{D}f(x)\,\mathrm {d} x,} 其中m (D )是D 的测度。
此外,这些积分在物理 中有大量应用。
力学 中,转动惯量 可以作为密度 乘以刚体和转轴的距离的平方的体积分(三重积分)计算:
I z = ∭ V ρ r 2 d V {\displaystyle I_{z}=\iiint _{V}\rho r^{2}\,\mathrm {d} V} 。与三维欧氏空间 R 3 中的质量测度 d m {\displaystyle \mathrm {d} m} 表示的质量分布 关联的引力势 :
V ( x ) = − ∫ R 3 G | x − y | d m ( y ) . {\displaystyle V(\mathbf {x} )=-\int _{\mathbf {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {y} |}}\,\mathrm {d} m(\mathbf {y} ).} 如果存在一个连续函数 ρ ( x ) {\displaystyle \rho (\mathbf {x} )} 表示 x 处的密度分布, 那么 d m ( x ) = ρ ( x ) d 3 x {\displaystyle \mathrm {d} m(\mathbf {x} )=\rho (\mathbf {x} )\,\mathrm {d} ^{3}\mathbf {x} } , 其中d 3 x {\displaystyle \mathrm {d} ^{3}\mathbf {x} } 是欧几里得体积元 , 那么引力势就是
V ( x ) = − ∫ R 3 G | x − y | ρ ( y ) d 3 y . {\displaystyle V(\mathbf {x} )=-\int _{\mathbf {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {y} |}}\,\rho (\mathbf {y} )\,\mathrm {d} ^{3}\mathbf {y} .} 在电磁学 中,麦克斯韦方程组 可以写作多重积分,用以计算总磁场和电场。下例中,由电荷分布 产生的电场 通过向量函数的三重积分 得到:
E → = 1 4 π ϵ 0 ∭ r → − r → ′ ‖ r → − r → ′ ‖ 3 ρ ( r → ′ ) d 3 r ′ {\displaystyle {\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'} 。 外部链接
编辑