第11章 多元函数的积分

Tip
在单变量积分中,我们通过定积分研究了一个变量在区间上的累积问题,例如曲线下的面积、物体的位移等。然而,现实世界中的许多问题涉及多个变量的相互作用,例如:三维空间中不均匀物体的质量(质量密度随位置变化),平面区域上的温度分布平均值等。这些问题要求我们将积分的概念从一维推广到高维空间。重积分(Multiple Integrals)正是解决这类问题的核心工具。本章将系统学习二重积分(Double Integrals)与三重积分(Triple Integrals)的定义、计算方法及其应用。

11.1 二重积分

Tip
跟单变量积分一样, 二重积分的本质仍然是“分割-近似-求和”的极限过程.

11.1.1 二重积分的概念

Tip

回顾: 单变量积分 我们要计算曲线 y=f(x)在区间 [a,b]下的面积,可以通过以下步骤实现: 将区间 [a,b]等分为 N个子区间,每个子区间的宽度为 Δx=baN。通过分割得到 N个窄矩形(小柱体),每个矩形的左端点或右端点记为 xi(例如 xi=a+iΔx).每个窄矩形的面积可近似为底边长度 Δx乘以高度 f(xi),即第 i个矩形的面积为 Δxf(xi).将所有窄矩形的面积相加,得到总面积 SN的近似值:

SN=i=1NΔxf(xi)

N(即子区间宽度 Δx0)时,近似值 SN趋近于精确面积 S

S=limNSN=limΔx0SN

最终,面积 S可表示为定积分:

S=abf(x)dx

Important

多变量积分

[图略,待补充] 如图所示,对于定义在平面区域 D 上的函数 z=f(x,y),若要计算其对应的体积或物理量(如质量、电荷量等),其思想是单变量积分在二维空间的自然推广.首先将积分区域 D 分割为 N 个面积微元 ΔA,每个微元可近似视为一个小矩形区域,并选取其代表点 (xi,yi)(例如取中心点或顶点坐标).每个微元的体积可近似为底面积 ΔA 乘以高度 f(xi,yi),即 ΔAf(xi,yi).将所有微元的贡献叠加,得到近似总体积:

VN=i=1NΔAf(xi,yi)

当分割无限细化(即 N,且每个微元面积 ΔA0)时,近似值趋于精确值:

V=limNVN=limΔA0VN

这一极限过程即为二重积分的定义,可表示为:

V=Df(x,y)dA

其中 D 是积分区域,f(x,y) 为被积函数,dA 为面积元(直角坐标系中通常写作 dxdy).

如何计算? [图略,待补充]

M=i=1Nρ(xi,yi)dAi=Dρ(x,y)dA

11.1.2 二重积分的计算

Note

例1:已知f(x,y)=1x2y2,D={(x,y)|0x1,0y1},计算Df(x,y)dA.

解:上式可展开为

Df(x,y)dA=Df(x,y)dxdy=0101(1x2y2)dxdy=01[01(1x2y2)dx]dy

Inner:

01(1x2y2)dx=[x13x3y2x]01=23y2

Outer:

01(23y2)dy=[23y13y3]01=13

同理

01[01(1x2y2)dx]dy=01[01(1x2y2)dy]dx=01([y13y3x2y]01)dx=01(23x2)dx=[23x13x3]01=13

例2:已知f(x,y)=1x2y2,D={(x,y)|x2+y21,x>0,y>0},计算Df(x,y)dA.

解:由图 [图略,待补充]

Df(x,y)dA=01[01x2(1x2y2)dy]dx

Inner:

01x2(1x2y2)dy=[yx2y13y3]01x2=1x2x21x213(1x2)32=1x2(1x213(1x2))=1x2(2323x2)=23(1x2)32

Outer:0123(1x2)32dx x=sinθ,dx=cosθdθ

0123(1x2)32dx=0π223cos4θdθ=320π2(1+cos2θ2)2dθ=320π2(14+12cos2θ+14cos22θ)dθ=32([θ4]0π2+[sin2θ4]0π2+140π2(1+cos4θ2)2dθ)=32(π8+π16+1814[sin4θ]0π2)=π8

Note

  • 重积分可交换顺序? 一般可以,注意上下限

例1:计算Dxydxdy,其中积分区域 D 由直线 y=1x=2y=x 围成.

解法一:由图 [图略,待补充]

Dxydxdy=12[1xxydy]dx=12([12xy2]1x)dx=12(x32x2)dx=[x48x24]12=1+18=98

解法二:由图 [图略,待补充]

Dxydxdy=12[y2xydx]dy=12[12x2y]y2dy=12(2yy32)dy=[y2y48]12=2(118)=98

例2:计算Dy1+x2y2dσ,其中D 由直线 y=xx=1y=1 围成的闭区域.

解法一:由图 [图略,待补充]

Dy1+x2y2dxdy=11dxx1y1+x2y2dy=11[[13(1+x2y2)32]x1]dx=1113(1|x|3)dx=2301(1x3)dx=12

解法二:由图 [图略,待补充]

Dy1+x2y2dσ=11y[1y1+x2y2dx]dy

Note

极坐标下的积分

例:已知f(x,y)=1x2y2,D={(x,y)|x2+y21,x0,y0},计算Df(x,y)dA.

解:由图 [图略,待补充] 由于

dA=dxdy=rdθdr
D(1x2y2)dxdy=D(1x2y2)rdrdθ=D(1r2)rdrdθ=0π2[01(1r2)rdr]dθ=π2[12r214r4]01=π8

11.2 多变量积分的换元法

Important

变量替换与雅可比行列式

考虑坐标变换 {u=u(x,y)v=v(x,y),其全微分形式可表示为:

{du=uxdx+uydydv=vxdx+vydy

该微分系统可写成矩阵形式:

[dudv]=[uxuyvxvy][dxdy]

[图略,待补充]

dudv=J=(u,v)(x,y)=|uxdxuydyvxdxvydy|=|uxvyuyvx|dxdy

 

Important

雅可比行列式与坐标变换的深层关系

在坐标变换 (x,y)(u,v) 中,雅可比行列式不仅描述面积元的缩放比例,更揭示了双向变换的对称性。由基本关系:

dudv=|J|dxdy=|(u,v)(x,y)|dxdy

通过反函数定理,可得逆变换关系:

dxdy=|(x,y)(u,v)|dudv

这建立了双向变换的精确对应:雅可比行列式互为倒数,即:

|(x,y)(u,v)|=|(u,v)(x,y)|1

经典示例:极坐标变换
对极坐标变换 {x=ρcosθy=ρsinθ,其雅可比矩阵为:

J(ρ,θ)=(x,y)(ρ,θ)=|xρxθyρyθ|=|cosθρsinθsinθρcosθ|=ρ

因此面积元变换为:

dxdy=|ρ|dρdθ=ρdρdθ

矩阵互逆性的几何诠释
变换矩阵的互逆性直接体现在坐标系转换中:

[xuxvyuyv][uxuyvxvy]

这两个雅可比矩阵的乘积为单位矩阵,其行列式乘积为1。例如当 Juv=5 时,逆变换的行列式为 15,此时:

dxdy=15dudv

这种精确的倒数关系确保了积分在不同坐标系下计算结果的一致性,是多重积分变量替换定理的核心数学基础。

 

Note

例:计算+ex2dx.

解: I=+ex2dx

I2=++e(x2+y2)dxdy=+ey2dy+ex2dx

x=rcosθ,y=rsinθ,则

I2=02π0+er2rdrdθ=02πdθ(120+er2dr2)=1202π[(er2)]0+dθ=π

+ex2dx=I=π

11.3 重积分的应用

Important

质量

对于分布在平面区域 D上的物体(如薄片),若其密度函数为 ρ(x,y),则总质量可通过二重积分计算:

M=Dρ(x,y)dσ

其中 dσ表示面积元(即直角坐标系中的 dxdy).当密度均匀时(ρ为常数),总质量退化为密度与面积的乘积 M=ρA,其中 A=Ddσ为区域 D的面积.

质心

质心是物体质量分布的加权平均位置.对于离散质点系,质心坐标 (x¯,y¯)定义为:

x¯=i=1nmixii=1nmi,y¯=i=1nmiyii=1nmi

将其推广到连续情形,则质心坐标通过二重积分表达。当密度函数为 μ(x,y)时:

x¯=Dxμ(x,y)dσDμ(x,y)dσ,y¯=Dyμ(x,y)dσDμ(x,y)dσ

Note

例1:如图所示,求I=Dr2(ρ,θ)ρdρdθ.

[图略,待补充] 解:

I=Dr2(ρ,θ)ρdρdθ=π2π202cosθρ3dρdθ=π2π2[14ρ4]02cosθdθ=0π24cos4θdθ=32π

例2:如图所示,求I=Dr2(ρ,θ)ρdρdθ.

[图略,待补充] 解法一:

I=Dr2(ρ,θ)ρdρdθ=π2π202cosθ(1+ρ22ρcosθ)ρdρdθ=π2

解法二:

I=02π01ρ2ρdρdθ=π2

Warning

转动惯量与质心的关系(不要求)

对于平面区域 D上的质量分布,若密度函数为 μ(x,y),则系统对某点 (x0,y0)的转动惯量为:

I(x0,y0)=Dμ(x,y)[(xx0)2+(yy0)2]dxdy

此式可视为质量元 μ(x,y)dxdy到参考点距离平方的加权积分.特别地,当参考点为原点时,转动惯量简化为:

I=Dμ(x,y)r2(x,y)dxdy(r2=x2+y2)

质心的极值特性:
质心具有重要物理性质:它是使系统转动惯量取极小值的位置.为证明这一点,对 I(x0,y0)分别关于 x0,y0求偏导并令其为零:

Ix0=D2μ(x,y)(x0x)dxdy=0Iy0=D2μ(x,y)(y0y)dxdy=0

由此可得方程组:

x0Dμdxdy=Dμxdxdyy0Dμdxdy=Dμydxdy

最终解得质心坐标为:

x0=DμxdxdyDμdxdy,y0=DμydxdyDμdxdy

这与离散质点系的质心公式一致:

{x0=miximi(质量加权平均)y0=miyimi(质量加权平均)

Warning

平行轴定理

平行轴定理揭示了物体绕不同轴转动惯量之间的关系:绕任意轴的转动惯量等于绕质心轴的转动惯量加上总质量与该轴到质心距离平方的乘积。其数学推导可通过二重积分展开:

设质心坐标为 (xc,yc),绕点 (x0,y0) 的转动惯量为:

I(x0,y0)=Dμ[(xx0)2+(yy0)2]dxdy=Dμ[(xxc+xcx0)2+(yyc+ycy0)2]dxdy

展开平方项并重组积分:

=Dμ[(xxc)2+(yyc)2]dxdyIc (绕质心的转动惯量)+Dμ[(xcx0)2+(ycy0)2]dxdyMr2 (总质量×距离平方)+2(xcx0)Dμ(xxc)dxdy0 (质心定义)+2(ycy0)Dμ(yyc)dxdy0 (质心定义)

关键消去项:根据质心坐标的定义 xc=1MDμxdxdy,交叉项中的积分 Dμ(xxc)dxdyDμ(yyc)dxdy 必然为零。最终得到简洁形式:

I(x0,y0)=Ic+Mr2

其中 r=(xcx0)2+(ycy0)2 为两轴间距离,M=Dμdxdy 为总质量.

 

 

Note

例3:求球体 x2+y2+z24a2 被圆柱面 x2+y2=2ax (a>0) 截得的含在圆柱面内的立体体积.

[图略,待补充] 解:

V=4D4a2x2y2dxdy=40π202acosθ4a2ρ2ρdρdθ=40π283a3(1sin3θ)dθ=163a3π8383a3=323a3(π223)

Note
换元法

(x,y)(u,v)

例1(椭圆面积):已知x2a2+y2b2=1,求Ddxdy.

解法一: [图略,待补充]

Ddxdy=40a0b1x2a2dydx=40ab1x2a2dy

y=asinθ,dy=acosθdθ,则

=40π2bcosθacosθdθ=4ab0π21+cos2θ2dθ=πab

解法二: x=au,y=bv,1u1,1v1,u2+v2=1,则

dx=adudy=bdvdxdy=adubdv
Ddxdy=Dabdudv=πab

 

Note

例2:已知{u=3x2yv=x+y

[图略,待补充] 解: 方法一:

|ω1×ω2|=|ω1||ω2|sinθ=|3211|=5

方法二:

J=|uxuyvxvy|=|3211|=5

dudv=(u,v)(x,y)dxdy=5dxdy
Ddxdy=D15dudv

11.4 三重积分

Note

例1:计算Ωxdxdydz,其中 Ω 由平面 x+2y+z=1 与三个坐标平面 (x=0, y=0, z=0) 围成.

[图略,待补充] 解法一:

Ωxdxdydz=01dz012(1z)dy012yzxdx=01dz012(1z)12(12yz)2dy=1201[(12yz)36]012(1z)dz=11201(1z)3dz=112[(1z)44]01=148.

解法二:

Ωxdxdydz=01dx012(1x)dy01x2yxdz=01xdx012(1x)(1x2y)dy=0114x(1x)2dx=148

例2:利用柱面坐标计算三重积分 Ωzdxdydz,其中Ω 是由曲面 z=x2+y2 与平面 z=4 围成的闭区域.

[图略,待补充] 解:

{x=ρcosθy=ρsinθz=z

dV=dxdydz=ρdρdθdz
Ωzdxdydz=Ωzρdρdθdz=04dz02πdθ0zzρdρ=04dz02π12z2dθ=04z2πdz=13π[z3]04=643π