概述
由一元积分演化而来
- 定积分
- 二重积分:直角坐标、极坐标
- 三重积分:柱面坐标、球坐标、直角坐标
1 二重积分的概念与性质
一、概念
引入:曲顶柱体的体积、平面薄片的质量
- 定义1(分割求和取极限):
∬
D
f
(
x
,
y
)
d
σ
=
lim
λ
→
0
∑
i
=
1
n
f
(
x
i
,
y
i
)
Δ
σ
i
iintlimits_{D}f(x,y),dsigma=limlimits_{lambdarightarrow 0}sumlimits_{i=1}^n f(x_i,y_i),Deltasigma_i
D∬f(x,y)dσ=λ→0limi=1∑nf(xi,yi)Δσi(和式的极限),其中
D
D
D为积分区域、
x
,
y
x,y
x,y为积分变量、
f
(
x
,
y
)
f(x,y)
f(x,y)为被积函数、
d
σ
dsigma
dσ为面积元素、
f
(
x
,
y
)
d
σ
f(x,y),dsigma
f(x,y)dσ为被积表达式、
∑
f
(
x
i
,
y
i
)
Δ
σ
i
sum f(x_i,y_i),Deltasigma_i
∑f(xi,yi)Δσ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=maxlimits_{(x_1,y_1),(x_2,y_2)rightarrowDeltatau_i}{|P_1P_2|}
λ=max{di}且di=(x1,y1),(x2,y2)→Δτimax{∣P1P2∣}.
类比于定积分: ∫ a b f ( x ) d x = lim λ → 0 ∑ i = 1 n f ( x i ) Δ x i int_a^bf(x),dx=limlimits_{lambdarightarrow 0}sumlimits_{i=1}^n f(x_i),Delta x_i ∫abf(x)dx=λ→0limi=1∑nf(xi)Δxi. - 说明:
- 连续必可积: f ( x , y ) f(x,y) f(x,y)在闭区间连续,二重积分必存在
- 面积元素: d σ = d x d y dsigma=dxdy dσ=dxdy
- 几何意义: f ( x , y ) f(x,y) f(x,y)在 D D D上的二重积分等于柱体体积的代数和(有正有负)
二、性质
- 可加,可数乘
- 单调,特别的 ∣ ∬ D f ( x , y ) d σ ∣ ≤ ∬ D ∣ f ( x , y ) ∣ d σ |iintlimits_{D}f(x,y),dsigma|leiintlimits_{D}|f(x,y)|,dsigma ∣D∬f(x,y)dσ∣≤D∬∣f(x,y)∣dσ.
- 估值不等式: m i n σ ≤ ∬ D f ( x , y ) d σ ≤ M A X σ minsigma leiintlimits_{D}f(x,y),dsigma le MAXsigma minσ≤D∬f(x,y)dσ≤MAXσ.
- 中值定理:存在点 ( x 0 , y 0 ) (x_0,y_0) (x0,y0),使得 ∬ D f ( x , y ) d σ = f ( x 0 , y 0 ) σ iintlimits_{D}f(x,y),dsigma=f(x_0,y_0)sigma D∬f(x,y)dσ=f(x0,y0)σ.
2 二重积分的计算
一、直角坐标
- X型:
V
=
V=
V=
∫
a
b
[
∫
φ
1
(
x
)
φ
2
(
x
)
f
(
x
,
y
)
d
y
]
d
x
=
∫
a
b
d
x
∫
φ
1
(
x
)
φ
2
(
x
)
f
(
x
,
y
)
d
y
int_a^bleft[int_{varphi_1(x)}^{varphi_2(x)}f(x,y),dyright],{color{red}dx} = int_a^b,{color{red}dx}int_{varphi_1(x)}^{varphi_2(x)}f(x,y),dy
∫ab[∫φ1(x)φ2(x)f(x,y)dy]dx=∫abdx∫φ1(x)φ2(x)f(x,y)dy.
Y型: V = V= V= ∫ c d [ ∫ ψ 1 ( y ) ψ 2 ( y ) f ( x , y ) d x ] d y = ∫ c d d y ∫ ψ 1 ( y ) ψ 2 ( y ) f ( x , y ) d x int_c^dleft[int_{psi_1(y)}^{psi_2(y)}f(x,y),dxright],{color{red}dy} = int_c^d,{color{red}dy}int_{psi_1(y)}^{psi_2(y)}f(x,y),dx ∫cd[∫ψ1(y)ψ2(y)f(x,y)dx]dy=∫cddy∫ψ1(y)ψ2(y)f(x,y)dx. - 分割,选择次序
- 做题时:
- 只要见到积分区域(有时需要拆分区域)具有对称性,就要想到考察被积函数的奇偶性。
- 考虑x、y能否轮换对称(关于y=x对称)。
- 直接积分图形复杂时,考虑可加(减)性 。
例1:
求
∬
D
x
2
e
−
y
2
d
x
d
y
,
其
中
D
是
以
(
0
,
0
)
,
(
1
,
1
)
,
(
0
,
1
)
为
顶
点
的
三
角
形
color{blue}求iintlimits_{D}x^2e^{-y^2}dxdy,其中D是以(0,0),(1,1),(0,1)为顶点的三角形
求D∬x2e−y2dxdy,其中D是以(0,0),(1,1),(0,1)为顶点的三角形
解:
∵
∫
e
−
y
2
d
y
无
法
解
出
∴
先
积
分
x
原
式
=
∫
0
1
d
y
∫
0
y
x
2
e
−
y
2
d
x
=
∫
0
1
e
−
y
2
⋅
y
3
3
d
y
=
∫
0
1
e
−
y
2
⋅
y
2
6
d
y
2
=
−
1
6
∫
0
1
u
d
e
−
u
=
−
1
3
e
{color{red}because int e^{-y^2}dy无法解出therefore 先积分x} \ 原式 = int_0^1 dyint_0^y x^2e^{-y^2}dx = int_0^1e^{-y^2}cdot dfrac{y^3}{3}dy {color{red}=} int_0^1e^{-y^2}cdot dfrac{y^2}{6}dy^2 = -dfrac{1}{6}int_0^1 u de^{-u} = -dfrac{1}{3e}
∵∫e−y2dy无法解出∴先积分x原式=∫01dy∫0yx2e−y2dx=∫01e−y2⋅3y3dy =∫01e−y2⋅6y2dy2=−61∫01ude−u=−3e1
例2:
求
∬
D
y
1
+
x
2
−
y
2
d
σ
,
其
中
D
是
由
y
=
1
,
x
=
−
1
及
y
=
x
所
围
成
的
区
域
color{blue}求iintlimits_{D}ysqrt{1+x^2-y^2}dsigma,其中D是由y=1,x=-1及y=x所围成的区域
求D∬y1+x2−y2dσ,其中D是由y=1,x=−1及y=x所围成的区域
解:
不
妨
以
X
型
区
域
计
算
:
原
式
=
∫
−
1
1
d
x
∫
x
1
y
1
+
x
2
−
y
2
d
y
=
.
.
.
=
−
1
3
∫
−
1
1
[
(
1
+
x
2
−
y
2
)
3
2
]
x
1
d
x
=
−
1
3
∫
−
1
1
(
∣
x
∣
3
−
1
)
d
x
=
−
2
3
∫
0
1
(
x
3
−
1
)
d
x
=
1
2
.
不
知
x
正
负
,
平
方
开
根
号
要
绝
对
值
,
不
然
答
案
为
奇
函
数
错
误
不妨以X型区域计算:原式 = int_{-1}^1dxint_x^1 ysqrt{1+x^2-y^2}dy \ = ...=-dfrac{1}{3}int_{-1}^1left[(1+x^2-y^2)^{dfrac{3}{2}}right]_x^1dx = -dfrac{1}{3}int_{-1}^1({color{red}|x|^3}-1 )dx = -dfrac{2}{3}int_0^1(x^3-1 )dx = dfrac{1}{2}. \ color{red}不知x正负,平方开根号要绝对值,不然答案为奇函数错误
不妨以X型区域计算:原式=∫−11dx∫x1y1+x2−y2dy=...=−31∫−11⎣⎡(1+x2−y2)23⎦⎤x1dx=−31∫−11(∣x∣3−1)dx=−32∫01(x3−1)dx=21.不知x正负,平方开根号要绝对值,不然答案为奇函数错误
例3:
求
两
个
底
圆
半
径
都
等
于
R
的
直
交
援
助
面
积
所
围
成
的
立
体
的
体
积
color{blue}求两个底圆半径都等于R的直交援助面积所围成的立体的体积
求两个底圆半径都等于R的直交援助面积所围成的立体的体积
解:
设
两
圆
柱
面
方
程
为
x
2
+
y
2
=
R
2
和
x
2
+
z
2
=
R
2
V
=
8
∫
0
R
d
x
∫
0
R
2
−
x
2
R
2
−
x
2
d
y
=
8
∫
0
R
(
R
2
−
x
2
)
d
x
=
16
3
R
3
.
设两圆柱面方程为x^2+y^2=R^2和x^2+z^2=R^2 \ V=8int_0^Rdxint_0^{sqrt{R^2-x^2}} sqrt{R^2-x^2},dy = 8int_0^R(R^2-x^2)dx = dfrac{16}{3}R^3.
设两圆柱面方程为x2+y2=R2和x2+z2=R2V=8∫0Rdx∫0R2−x2R2−x2dy=8∫0R(R2−x2)dx=316R3.
例4:
更
换
积
分
次
序
,
I
=
∫
0
4
d
y
∫
−
4
−
y
4
y
−
y
2
f
(
x
,
y
)
d
x
color{blue}更换积分次序,I=int_0^4,dyint_{-sqrt{4-y}}^{sqrt{4y-y^2}}f(x,y),dx
更换积分次序,I=∫04dy∫−4−y4y−y2f(x,y)dx.
解:
I
=
∫
−
2
0
d
x
∫
0
4
−
x
2
f
(
x
,
y
)
d
y
+
∫
0
2
d
x
∫
2
−
4
−
x
2
2
+
4
−
x
2
f
(
x
,
y
)
d
y
.
I=int_{-2}^0,dxint_0^{4-x^2}f(x,y),dy+int_0^2,dxint_{2{color{red}-}sqrt{4-x^2}}^{2{color{red}+}sqrt{4-x^2}}f(x,y),dy.
I=∫−20dx∫04−x2f(x,y)dy+∫02dx∫2−4−x22+4−x2f(x,y)dy.
例5:
计
算
z
=
x
y
,
x
+
y
+
z
=
1
,
z
=
0
所
围
闭
合
区
域
的
体
积
color{blue}计算z=xy,x+y+z=1,z=0所围闭合区域的体积
计算z=xy,x+y+z=1,z=0所围闭合区域的体积.
解:法1
分
为
两
部
分
,
一
部
分
z
=
x
y
为
顶
,
另
一
部
分
z
=
1
−
x
−
y
为
顶
∵
x
y
=
1
−
x
−
y
⟺
y
=
1
−
x
1
+
x
∴
V
=
∫
0
1
d
x
∫
0
1
−
x
1
+
x
x
y
d
y
+
∫
0
1
d
x
∫
1
−
x
1
+
x
1
−
x
1
−
x
−
y
d
y
.
分为两部分,一部分z=xy为顶,另一部分z=1-x-y为顶 \ because xy=1-x-y iff y=dfrac{1-x}{1+x} \ therefore V=int_0^1,dxint_0^{frac{1-x}{1+x}}xy,dy+int_0^1,dxint_{frac{1-x}{1+x}}^{1-x}1-x-y,dy.
分为两部分,一部分z=xy为顶,另一部分z=1−x−y为顶∵xy=1−x−y⟺y=1+x1−x∴V=∫01dx∫01+x1−xxydy+∫01dx∫1+x1−x1−x1−x−ydy.
法2:参考例9,截面法
例6:
平
面
区
域
D
由
曲
线
{
x
=
t
−
sin
t
,
y
=
1
−
cos
t
(
0
≤
t
≤
2
π
)
与
x
轴
围
成
,
计
算
∬
D
(
x
+
2
y
)
d
x
d
y
color{blue}平面区域D由曲线begin{cases}x=t-sin t, \ y=1-cos tend{cases}(0le tle 2pi)与x轴围成,计算iintlimits_D (x+2y),dxdy
平面区域D由曲线{x=t−sint,y=1−cost(0≤t≤2π)与x轴围成,计算D∬(x+2y)dxdy.
将
参
数
方
程
选
做
换
元
积
分
的
变
量
替
换
设
D
:
0
≤
x
≤
2
π
,
0
≤
y
≤
g
(
x
)
,
则
∬
D
(
x
+
2
y
)
d
x
d
y
=
∫
0
2
π
d
x
∫
0
g
(
x
)
(
x
+
2
y
)
d
y
=
∫
0
2
π
[
x
g
(
x
)
+
g
2
(
x
)
]
d
x
=
∫
0
2
π
[
(
t
−
sin
t
)
(
1
−
cos
t
)
+
(
1
−
cos
t
)
2
]
(
1
−
cos
t
)
d
t
=
.
.
.
=
3
π
2
+
5
π
.
color{red}将参数方程选做换元积分的变量替换 \ 设D:0le xle 2pi,0le yle g(x),则 \ iintlimits_D (x+2y),dxdy=int_0^{2pi},dxint_0^{g(x)}(x+2y),dy=int_0^{2pi}[xg(x)+g^2(x)],dx \ =int_0^{2pi}[(t-sin t)(1-cos t)+(1-cos t)^2](1-cos t),dt =,...,=3pi^2+5pi.
将参数方程选做换元积分的变量替换设D:0≤x≤2π,0≤y≤g(x),则D∬(x+2y)dxdy=∫02πdx∫0g(x)(x+2y)dy=∫02π[xg(x)+g2(x)]dx=∫02π[(t−sint)(1−cost)+(1−cost)2](1−cost)dt=...=3π2+5π.
二、极坐标
- 由于
Δ
σ
i
≈
ρ
i
Δ
ρ
i
Δ
θ
i
Deltasigma_iapprox rho_i Deltarho_i Deltatheta_i
Δσi≈ρiΔρiΔθi,元素面积可以表示为
d
σ
=
ρ
d
ρ
d
θ
dsigma = rho drho dtheta
dσ=ρdρdθ.
因此 lim λ → 0 ∑ i = 1 n f ( x i , y i ) Δ σ i = lim λ → 0 ∑ i = 1 n f ( ρ i cos θ i , ρ i sin θ i ) ρ i Δ ρ i Δ θ i limlimits_{lambdarightarrow 0}sumlimits_{i=1}^n f(x_i,y_i),Deltasigma_i = limlimits_{lambdarightarrow 0}sumlimits_{i=1}^n f(rho_icostheta_i,rho_isintheta_i),rho_iDeltarho_iDeltatheta_i λ→0limi=1∑nf(xi,yi)Δσi=λ→0limi=1∑nf(ρicosθi,ρisinθi)ρiΔρiΔθi,
即== ∬ D f ( x , y ) d σ = ∬ D f ( ρ cos θ , ρ sin θ ) ρ d ρ d θ iintlimits_{D}f(x,y),dsigma = iintlimits_{D}f(rhocostheta,rhosintheta),rho,drho,dtheta D∬f(x,y)dσ=D∬f(ρcosθ,ρsinθ)ρdρdθ==. - 其中按顺序: ∬ D f ( ρ cos θ , ρ sin θ ) ρ d ρ d θ = ∫ α β d θ ∫ φ 1 ( θ ) φ 2 ( θ ) f ( ρ cos θ , ρ sin θ ) ρ d ρ iintlimits_{D}f(rhocostheta,rhosintheta),rho,drho,dtheta = int_{alpha}^{beta},dtheta int_{varphi_1(theta)}^{varphi_2(theta)}f(rhocostheta,rhosintheta),rho,drho D∬f(ρcosθ,ρsinθ)ρdρdθ=∫αβdθ∫φ1(θ)φ2(θ)f(ρcosθ,ρsinθ)ρdρ.
例7:
转
换
成
直
角
坐
标
系
I
=
∫
0
π
2
d
θ
∫
a
cos
θ
2
a
cos
θ
4
a
2
−
ρ
2
ρ
d
ρ
color{blue}转换成直角坐标系I= int_0^{frac{pi}{2}},dtheta int_{acostheta}^{2acostheta}sqrt{4a^2-rho ^2},rho,drho
转换成直角坐标系I=∫02πdθ∫acosθ2acosθ4a2−ρ2ρdρ.
解:
难
点
:
画
图
I
=
I
=
∫
0
a
d
x
∫
a
x
−
x
2
2
a
x
−
x
2
4
a
2
−
x
2
−
y
2
d
y
+
∫
0
2
a
d
x
∫
0
2
a
x
−
x
2
4
a
2
−
x
2
−
y
2
d
y
.
难点:画图 I=I=int_0^a,dxint_{sqrt{ax-x^2}}^{sqrt{2ax-x^2}}sqrt{4a^2-x ^2-y^2},dy+int_0^{2a},dxint_0^{sqrt{2ax-x^2}}sqrt{4a^2-x ^2-y^2},dy.
难点:画图I=I=∫0adx∫ax−x22ax−x24a2−x2−y2dy+∫02adx∫02ax−x24a2−x2−y2dy.
例8:
闭
区
域
D
:
x
2
+
y
2
≤
y
(
1
)
,
x
≥
0
(
2
)
,
f
(
x
,
y
)
为
D
上
连
续
函
数
,
且
f
(
x
,
y
)
=
1
−
x
2
−
y
2
−
8
π
∬
D
f
(
u
,
v
)
d
u
d
v
(
∗
)
,
求
f
(
x
,
y
)
color{blue}闭区域D:x^2+y^2le y(1),xge 0(2),f(x,y)为D上连续函数,且f(x,y)=sqrt{1-x^2-y^2}-dfrac{8}{pi}iintlimits_ D f(u,v),dudv(*),求f(x,y)
闭区域D:x2+y2≤y(1),x≥0(2),f(x,y)为D上连续函数,且f(x,y)=1−x2−y2−π8D∬f(u,v)dudv(∗),求f(x,y).
解:
设
A
=
∬
D
f
(
u
,
v
)
d
u
d
v
,
又
∬
D
d
x
d
y
=
π
8
,
则
在
(
∗
)
两
边
求
二
重
积
分
得
到
:
A
=
∬
D
1
−
x
2
−
y
2
d
x
d
y
−
∬
D
8
π
A
d
x
d
y
.
∴
A
=
1
2
∬
D
1
−
x
2
−
y
2
d
x
d
y
=
1
2
∫
0
π
2
d
θ
∫
0
sin
θ
ρ
1
−
ρ
2
d
ρ
=
1
2
∫
0
π
2
[
−
1
3
(
1
−
ρ
2
)
3
2
]
0
sin
θ
d
θ
=
π
12
−
1
9
.
∴
f
(
x
,
y
)
=
1
−
x
2
−
y
2
−
8
π
(
π
12
−
1
9
)
.
设A=iintlimits_ D f(u,v),dudv,又iintlimits_ D ,dxdy=dfrac{pi}{8}, \则{color{red}在(*)两边求二重积分}得到:A=iintlimits_ Dsqrt{1-x^2-y^2},dxdy-iintlimits_ Ddfrac{8}{pi}A,dxdy. \ therefore A=dfrac{1}{2}iintlimits_ Dsqrt{1-x^2-y^2},dxdy = dfrac{1}{2}int_0^{piover2},dthetaint_0^{sintheta}rhosqrt{1-rho^2},drho = dfrac{1}{2}int_0^{piover2}left[-dfrac{1}{3}(1-rho^2)^{3over2}right]_0^{sintheta} ,dtheta = dfrac{pi}{12}-dfrac{1}{9}. \ therefore f(x,y) = sqrt{1-x^2-y^2}-dfrac{8}{pi}(dfrac{pi}{12}-dfrac{1}{9}).
设A=D∬f(u,v)dudv,又D∬dxdy=8π,则在(∗)两边求二重积分得到:A=D∬1−x2−y2dxdy−D∬π8Adxdy.∴A=21D∬1−x2−y2dxdy=21∫02πdθ∫0sinθρ1−ρ2dρ=21∫02π[−31(1−ρ2)23]0sinθdθ=12π−91.∴f(x,y)=1−x2−y2−π8(12π−91).
三、坐标变换
3 三重积分
一、概念
定义1:
∭
Ω
f
(
x
,
y
,
z
)
d
σ
=
lim
λ
→
0
∑
i
=
1
n
f
(
x
i
,
y
i
,
z
1
)
Δ
σ
i
iiintlimits_{Omega}f(x,y,z),dsigma=limlimits_{lambdarightarrow 0}sumlimits_{i=1}^n f(x_i,y_i,z_1),Deltasigma_i
Ω∭f(x,y,z)dσ=λ→0limi=1∑nf(xi,yi,z1)Δσi.
三重积分也记作:
∭
Ω
f
(
x
,
y
,
z
)
d
σ
=
∭
Ω
f
(
x
,
y
,
z
)
d
x
d
y
d
z
iiintlimits_{Omega}f(x,y,z),dsigma=iiintlimits_{Omega}f(x,y,z),dx,dy,dz
Ω∭f(x,y,z)dσ=Ω∭f(x,y,z)dxdydz.
二、计算
-
先一后二(投影法”压薄“):
∭ Ω f ( x , y , z ) d x d y d z = ∬ D F ( x , y ) d x d y = ∬ D [ ∫ z 1 ( x , y ) z 2 ( x , y ) f ( x , y , z ) d z ] d x d y iiintlimits_{Omega}f(x,y,z),dx,dy,dz = iintlimits_{D}F(x,y),dx,dy = iintlimits_{D} left[int_{z_1(x,y)}^{z_2(x,y)}f(x,y,z),dz right],dx,dy Ω∭f(x,y,z)dxdydz=D∬F(x,y)dxdy=D∬[∫z1(x,y)z2(x,y)f(x,y,z)dz]dxdy -
先二后一(截面法”集中“):
∭ Ω f ( x , y , z ) d x d y d z = ∫ C 1 C 2 F ( x ) d z = ∫ C 1 C 2 [ ∬ D z f ( x , y , z ) d x d y ] d z iiintlimits_{Omega}f(x,y,z),dx,dy,dz = int_{C_1}^{C_2} F(x),dz = int_{C_1}^{C_2} left[iintlimits_{D_z}f(x,y,z),dx,dyright] ,dz Ω∭f(x,y,z)dxdydz=∫C1C2F(x)dz=∫C1C2⎣⎡Dz∬f(x,y,z)dxdy⎦⎤dz
例9: 计 算 I = ∭ Ω e ∣ z ∣ d x d y d z , 其 中 Ω 为 x 2 + y 2 + z 2 ≤ 1 所 围 成 的 立 体 color{blue}计算I=iiintlimits_{Omega}e^{|z|},dxdydz,其中Omega为x^2+y^2+z^2le 1所围成的立体 计算I=Ω∭e∣z∣dxdydz,其中Ω为x2+y2+z2≤1所围成的立体.
解:
∵ e ∣ z ∣ 为 偶 函 数 , 采 用 先 二 后 一 的 思 想 , I = 2 ∫ 0 1 e ∣ z ∣ d z ∬ D x y d x d y . 其 中 ∬ D x y d x d y = π ( r ( z ) ) 2 , r ( z ) = 1 − z 2 . ∴ I = 2 ∫ 0 1 e z ( 1 − z 2 ) d z = 2 π [ − ( z − 1 ) 2 e z ] 0 1 = 2 π . because e^{|z|}为偶函数,采用先二后一的思想,\ I=2int_0^1e^{|z|},dz iintlimits_{D_{xy}},dxdy. 其中iintlimits_{D_{xy}},dxdy=pi (r(z))^2,r(z)=sqrt{1-z^2}. \ therefore I = 2int_0^1 e^z(1-z^2),dz = 2pileft[-(z-1)^2e^zright]_0^1 = 2pi. ∵e∣z∣为偶函数,采用先二后一的思想,I=2∫01e∣z∣dzDxy∬dxdy. 其中Dxy∬dxdy=π(r(z))2,r(z)=1−z2.∴I=2∫01ez(1−z2)dz=2π[−(z−1)2ez]01=2π.
例10: 计 算 I = ∭ Ω x d x d y d z , 其 中 Ω 由 z = x y , x + y + z = 1 及 z = 0 所 围 成 的 立 体 color{blue}计算I=iiintlimits_{Omega}x,dxdydz,其中Omega由z=xy,x+y+z=1及z=0所围成的立体 计算I=Ω∭xdxdydz,其中Ω由z=xy,x+y+z=1及z=0所围成的立体.
解:
发 现 当 x 为 一 定 值 时 , 平 行 于 O y z 面 的 切 面 切 立 体 的 截 面 为 一 过 原 点 的 三 角 形 , 则 原 式 表 示 成 I = ∫ 0 1 x d x ∬ D y z S ( x ) d y d z . 对 于 S ( x ) , 当 x = k 时 , 则 z = k y , z = 1 − k − y , z = 0 , sin θ = k 1 + k 2 有 S = 1 2 ( 1 − k ) ( 1 + k 2 1 − k 1 + k ) sin θ = k ( 1 − k ) 2 2 ( 1 + k ) . ∴ I = 1 2 ∫ 0 1 x ( 1 − x ) 2 1 + x d x + 1 = 1 2 ∫ 1 2 ( t − 2 ) 2 ( t − 1 ) 2 t d t = 1 2 ∫ 0 1 t 3 − 6 t 2 + 13 t − 12 + 4 t d t = 2 l n 2 − 11 8 . 发现当x为一定值时,平行于O_{yz}面的切面切立体的截面为一{color{red}过原点的三角形},则原式表示成I=int_0^1 x,dx iintlimits_{D_{yz}}S(x),dydz. \ 对于S(x),当x=k时,则z=ky,z=1-k-y,z=0,sintheta=dfrac{k}{sqrt{1+k^2}} \ 有S=dfrac{1}{2}(1-k)(sqrt{1+k^2}dfrac{1-k}{1+k})sintheta = dfrac{k(1-k)^2}{2(1+k)}. \ therefore I = dfrac{1}{2}int_0^1 dfrac{x(1-x)^2}{1+x},dx+1 = dfrac{1}{2}int_1^2 dfrac{(t-2)^2(t-1)^2}{t},dt = dfrac{1}{2}int_0^1 t^3-6t^2+13t-12+dfrac{4}{t},dt = 2ln2-dfrac{11}{8}. 发现当x为一定值时,平行于Oyz面的切面切立体的截面为一过原点的三角形,则原式表示成I=∫01xdxDyz∬S(x)dydz.对于S(x),当x=k时,则z=ky,z=1−k−y,z=0,sinθ=1+k2k有S=21(1−k)(1+k21+k1−k)sinθ=2(1+k)k(1−k)2.∴I=21∫011+xx(1−x)2dx+1=21∫12t(t−2)2(t−1)2dt=21∫01t3−6t2+13t−12+t4dt=2ln2−811.
-
利用柱面坐标(case:被积函数具有 f ( x 2 + y 2 ) f(sqrt{x^2+y^2}) f(x2+y2)的形式,积分区域是直柱体,且投影面用极坐标表示方便):
{ x = ρ cos θ . y = ρ sin θ . z = z . begin{cases}x=rhocostheta. \ y=rhosintheta. \z=z. end{cases} ⎩⎪⎨⎪⎧x=ρcosθ.y=ρsinθ.z=z.
∭ Ω f ( x , y , z ) d x d y d z = ∭ Ω f ( ρ cos θ , ρ sin θ , z ) ρ d ρ d θ d z = ∫ α β d θ ∫ φ 1 φ 2 ρ d ρ ∫ z 1 z 2 f ( ρ cos θ , ρ sin θ , z ) d z iiintlimits_{Omega}f(x,y,z),dx,dy,dz = iiintlimits_{Omega}f(rhocostheta, rhosintheta, z)rho,drho dtheta dz = int_{alpha}^{beta},dtheta int_{varphi_1}^{varphi_2},{color{red}rho},drhoint_{z_1}^{z_2}f(rhocostheta,rhosintheta,z),dz Ω∭f(x,y,z)dxdydz=Ω∭f(ρcosθ,ρsinθ,z)ρdρdθdz=∫αβdθ∫φ1φ2ρdρ∫z1z2f(ρcosθ,ρsinθ,z)dz例11: 将 I = ∭ Ω f ( x , y , z ) d x d y d z 表 示 成 柱 面 坐 标 系 的 累 次 积 分 , 其 中 Ω 是 由 2 z = x 2 + y 2 , z = 2 , z = 1 所 围 成 的 立 体 color{blue}将I=iiintlimits_{Omega}f(x,y,z),dxdydz表示成柱面坐标系的累次积分,其中Omega是由2z=x^2+y^2,z=2,z=1所围成的立体 将I=Ω∭f(x,y,z)dxdydz表示成柱面坐标系的累次积分,其中Ω是由2z=x2+y2,z=2,z=1所围成的立体.
解:
I = ∫ 0 2 π d θ ∫ 0 2 ρ d ρ ∫ 1 2 f ( ρ cos θ , ρ sin θ , z ) d z + ∫ 0 2 π d θ ∫ 2 2 ρ d ρ ∫ ρ 2 2 2 f ( ρ cos θ , ρ sin θ , z ) d z . I=int_0^{2pi},dtheta int_0^{sqrt{2}}rho,drho int_1^2f(rhocostheta,rhosintheta,z),dz+int_0^{2pi},dtheta int_{sqrt{2}}^2rho,drho int_{rho^2over2}^2f(rhocostheta,rhosintheta,z),dz. I=∫02πdθ∫02ρdρ∫12f(ρcosθ,ρsinθ,z)dz+∫02πdθ∫22ρdρ∫2ρ22f(ρcosθ,ρsinθ,z)dz.
-
利用球面坐标(case:被积函数具有 f ( x 2 + y 2 + z 3 ) f(sqrt{x^2+y^2+z^3}) f(x2+y2+z3)的形式,积分区域是球面或锥面):
{ x = r sin φ cos θ . y = r sin φ sin θ . z = r cos φ . begin{cases}x=rsinvarphicostheta. \ y=rsinvarphisintheta. \z=rcosvarphi. end{cases} ⎩⎪⎨⎪⎧x=rsinφcosθ.y=rsinφsinθ.z=rcosφ.
∭ Ω f ( x , y , z ) d x d y d z = ∭ Ω f ( r sin φ cos θ , r sin φ sin θ , r cos φ ) r 2 sin φ d r d φ d θ iiintlimits_{Omega}f(x,y,z),dx,dy,dz = iiintlimits_{Omega}f(rsinvarphicostheta, rsinvarphisintheta, rcosvarphi){color{red}r^2sinvarphi},dr,dvarphi,dtheta Ω∭f(x,y,z)dxdydz=Ω∭f(rsinφcosθ,rsinφsinθ,rcosφ)r2sinφdrdφdθ.例12: 将 I = ∭ Ω f ( x , y , z ) d x d y d z 表 示 成 球 面 坐 标 系 的 累 次 积 分 , 其 中 Ω 是 由 2 z = x 2 + y 2 , z = 2 , z = 1 所 围 成 的 立 体 color{blue}将I=iiintlimits_{Omega}f(x,y,z),dxdydz表示成球面坐标系的累次积分,其中Omega是由2z=x^2+y^2,z=2,z=1所围成的立体 将I=Ω∭f(x,y,z)dxdydz表示成球面坐标系的累次积分,其中Ω是由2z=x2+y2,z=2,z=1所围成的立体.
解:
I = ∫ 0 2 π d θ ∫ 0 π 4 sin φ d φ ∫ 1 cos φ 2 cos φ f ( r sin φ cos θ , r sin φ sin θ , r cos φ ) r 2 d r + ∫ 0 2 π d θ ∫ π 4 arctan 2 sin φ d φ ∫ 1 cos φ 2 cos φ sin 2 φ f ( r sin φ cos θ , r sin φ sin θ , r cos φ ) r 2 d r . I=int_0^{2pi},dtheta int_0^{piover4}sinvarphi,dvarphi int_{1over{cosvarphi}}^{2over{cosvarphi}}f(rsinvarphicostheta, rsinvarphisintheta, rcosvarphi)r^2,dr \+ int_0^{2pi},dtheta int_{piover4}^{arctan{sqrt2}}sinvarphi,dvarphi {color{red}int_{1over{cosvarphi}}^{2cosvarphi over sin^2varphi}} f(rsinvarphicostheta, rsinvarphisintheta, rcosvarphi)r^2,dr. I=∫02πdθ∫04πsinφdφ∫cosφ1cosφ2f(rsinφcosθ,rsinφsinθ,rcosφ)r2dr+∫02πdθ∫4πarctan2sinφdφ∫cosφ1sin2φ2cosφf(rsinφcosθ,rsinφsinθ,rcosφ)r2dr.
4 重积分的应用
一、曲面面积
- 由定义可知:
- ∬ D 1 d σ = iintlimits_D 1,dsigma = D∬1dσ=平面区域面积
- ∭ Ω 1 d V = ∬ D x y z d σ = iiintlimits_{Omega} 1,dV= iintlimits_{D_{xy}}z,dsigma = Ω∭1dV=Dxy∬zdσ=空间区域体积
- 二重积分计算曲面区域
z
=
f
(
x
,
y
)
z=f(x,y)
z=f(x,y)的面积:
由于 n ⃗ = ( f x , f y , − 1 ) vec n=(f_x,f_y,-1) n=(fx,fy,−1)且 ∣ cos α ∣ = n ⃗ × ( 0 , 0 , 1 ) |cos alpha|= vec ntimes(0,0,1) ∣cosα∣=n×(0,0,1),故 ∣ cos α ∣ = 1 1 + f x 2 ( x , y ) + f y 2 ( x , y ) |cos alpha|= dfrac{1}{sqrt{1+f_x^2(x,y)+f_y^2(x,y)}} ∣cosα∣=1+fx2(x,y)+fy2(x,y)1
得到 S = ∬ D 1 + f x 2 ( x , y ) + f y 2 ( x , y ) d σ S=iintlimits_ Dsqrt{1+f_x^2(x,y)+f_y^2(x,y)},dsigma S=D∬1+fx2(x,y)+fy2(x,y)dσ.
二、质心
- 质心:使物体上的点 P P P对空间任意一点 P 0 P_0 P0力矩积分为0,即 ∭ Ω ( O P → − O P 0 → ) ρ ( x , y , z ) d V = 0 iiintlimits_{Omega}(overrightarrow{OP}- overrightarrow{OP_0})rho(x,y,z),dV=0 Ω∭(OP−OP0)ρ(x,y,z)dV=0。
- 故质心坐标
(
x
0
,
y
0
,
z
0
)
(x_0,y_0,z_0)
(x0,y0,z0)满足:
{
∭
Ω
(
x
−
x
0
)
ρ
(
x
,
y
,
z
)
d
V
=
0
∭
Ω
(
y
−
y
0
)
ρ
(
x
,
y
,
z
)
d
V
=
0
∭
Ω
(
z
−
z
0
)
ρ
(
x
,
y
,
z
)
d
V
=
0
begin{cases}iiintlimits_{Omega}(x-x_0)rho(x,y,z),dV=0 \ iiintlimits_{Omega}(y-y_0)rho(x,y,z),dV=0 \ iiintlimits_{Omega}(z-z_0)rho(x,y,z),dV=0 end{cases}
⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧Ω∭(x−x0)ρ(x,y,z)dV=0Ω∭(y−y0)ρ(x,y,z)dV=0Ω∭(z−z0)ρ(x,y,z)dV=0,
即
{ x 0 = ∭ Ω x ρ ( x , y , z ) d V m y 0 = ∭ Ω y ρ ( x , y , z ) d V m z 0 = ∭ Ω z ρ ( x , y , z ) d V m color{red}begin{cases}x_0=dfrac{iiintlimits_{Omega}xrho(x,y,z),dV}{m} \ y_0=dfrac{iiintlimits_{Omega}yrho(x,y,z),dV}{m} \ z_0=dfrac{iiintlimits_{Omega}zrho(x,y,z),dV}{m} end{cases} ⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎪⎧x0=mΩ∭xρ(x,y,z)dVy0=mΩ∭yρ(x,y,z)dVz0=mΩ∭zρ(x,y,z)dV
其中 m = ∭ Ω ρ ( x , y , z ) d V m=iiintlimits_{Omega}rho(x,y,z),dV m=Ω∭ρ(x,y,z)dV即物体的质量。 - 假如物体为一薄板,且水平放置在Oxy面上,则
{ x 0 = ∬ D x ρ ( x , y ) d σ m y 0 = ∭ D y ρ ( x , y ) d σ m begin{cases}x_0=dfrac{iintlimits_ D xrho(x,y),dsigma}{m} \ y_0=dfrac{iiintlimits_ D yrho(x,y),dsigma}{m} end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧x0=mD∬xρ(x,y)dσy0=mD∭yρ(x,y)dσ
三、转动惯量
- 整个物体对于坐标轴的转动变量:
例如对于 z z z轴: J z = ∭ Ω ( x 2 + y 2 ) ρ ( x , y , z ) d V J_z=iiintlimits_{Omega}(x^2+y^2)rho(x,y,z),dV Jz=Ω∭(x2+y2)ρ(x,y,z)dV. - 整个物体对于一点
P
0
(
x
0
,
y
0
,
z
0
)
P_0(x_0,y_0,z_0)
P0(x0,y0,z0)的转动变量:
J
P
0
=
∭
Ω
[
(
x
−
x
0
)
2
+
(
y
−
y
0
)
2
+
(
z
−
z
0
)
2
]
ρ
(
x
,
y
,
z
)
d
V
J_{P_0}=iiintlimits_{Omega}left[(x-x_0)^2+(y-y_0)^2+(z-z_0)^2 right]rho(x,y,z),dV
JP0=Ω∭[(x−x0)2+(y−y0)2+(z−z0)2]ρ(x,y,z)dV,
特别的,物体关于原点的转动惯量为 J 0 = ∭ Ω ( x 2 + y 2 + z 2 ) ρ ( x , y , z ) d V J_0=iiintlimits_{Omega}(x^2+y^2+z^2)rho(x,y,z),dV J0=Ω∭(x2+y2+z2)ρ(x,y,z)dV. - 由于一个质点关于一个平面的转动惯量为质量乘以质点到平面最短距离之平方,故整个物体对于坐标面的转动惯量:
例如对于 O x y Oxy Oxy面: J x y = ∭ Ω z 2 ρ ( x , y , z ) d V J_{xy}=iiintlimits_{Omega}z^2rho(x,y,z),dV Jxy=Ω∭z2ρ(x,y,z)dV.
四、引力
- 物体对于质点
P
0
P_0
P0的总引力为
F = G ∭ Ω M 0 d m r 2 ⋅ r 0 = G ∭ Ω M 0 ⋅ r ⃗ ( x , y , z ) ⋅ ρ ( x , y , z ) d V r 3 F=Giiintlimits_{Omega}dfrac{M_0,dm}{r^2}cdot r_0 = Giiintlimits_{Omega}dfrac{M_0cdot vec r(x,y,z)cdot rho(x,y,z),dV}{r^3} F=GΩ∭r2M0dm⋅r0=GΩ∭r3M0⋅r(x,y,z)⋅ρ(x,y,z)dV
其中 r 0 = r ⃗ ( x , y , z ) r r_0=dfrac{vec r(x,y,z)}{r} r0=rr(x,y,z)为单位向量, r = ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 r=sqrt{(x-x_0)^2+(y-y_0)^2+(z-z_0)^2} r=(x−x0)2+(y−y0)2+(z−z0)2, r ⃗ ( x , y , z ) = ( x − x 0 , y − y 0 , z − z 0 ) vec r(x,y,z)=(x-x_0,y-y_0,z-z_0) r(x,y,z)=(x−x0,y−y0,z−z0),则坐标分量为
{ F x = G M 0 ∭ Ω x − x 0 r 3 ρ ( x , y , z ) d V F y = G M 0 ∭ Ω y − y 0 r 3 ρ ( x , y , z ) d V F z = G M 0 ∭ Ω z − z 0 r 3 ρ ( x , y , z ) d V color{red}begin{cases}F_x=GM_0iiintlimits_{Omega}dfrac{x-x_0}{r^3}rho(x,y,z),dV \ F_y=GM_0iiintlimits_{Omega}dfrac{y-y_0}{r^3}rho(x,y,z),dV \ F_z=GM_0iiintlimits_{Omega}dfrac{z-z_0}{r^3}rho(x,y,z),dV end{cases} ⎩⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎧Fx=GM0Ω∭r3x−x0ρ(x,y,z)dVFy=GM0Ω∭r3y−y0ρ(x,y,z)dVFz=GM0Ω∭r3z−z0ρ(x,y,z)dV
最后
以上就是心灵美玫瑰为你收集整理的重积分(微积分)的全部内容,希望文章能够帮你解决重积分(微积分)所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复