多元函數積分的幾種特殊解法

法; 重積分的等值面 (線) 法; 重積分及曲面積分的微元法。

關鍵詞. 多元函數積分、“先二後一” 法、 等值面 (線) 法、 微元法、 面積、 體積。

多元函數積分一般是指曲線積分、 曲面積分、 二重及二重以上的積分。 這一類積分的計算 往往是根據積分的定義、 積分的幾何意義或物理意義、 積分的變量替換等方法轉化成計算一個 或若干個定積分來完成的。

眾所周知, 一個積分量是由積分區域及積分的被積函數決定的。 因而當積分的積分區域或 被積函數具有某種特殊性時, 我們也可採用特殊的方法去求解。 本文的目的就是介紹幾種較特 殊的方法來求解一類具有特殊性的多元函數數積分。 需要說明的是, 文中所引的例子均有其他 解法, 讀者不妨將其他解法與本文所給解法作一比較。

一 . 三重積分的 “先二後一” 法

三重積分可以通過轉化為三個累次積分進行計算。 但當用垂直於某一坐標軸 (如 z 軸) 的 平面去截積分區域所得的截面有某種規律時, 則可先在截面上積分, 然後再關於 z 積分, 這種 先計算某兩個變量的二重積分再計算另一個變量積分的方法通常稱為 “先二後一” 法或 “先重 後單” 法, 也稱 “坐標軸投影法”, 而這種方法實質上是定積分中用截面法求體積方法的推廣。

設 Ω = {(x, y, z) | (x, y) ∈ D(z), z1 ≤ z ≤ z2}, 其中 D(z) 是 x, y 平面上隨 z 連 續變化的有界閉區域 (如圖 1 所示)。 如果 f (x, y, z) 在 Ω 上有界可積, 對任意的 z ∈ [z¹, z₂], f(x, y, z) 作為 x, y 的函數在 D(z) 上可積, 則

Z Z Z

Ω

f(x, y, z)dxdydz =

Z

z2

z1

dz

Z Z

D(z)

f(x, y, z)dxdy, (1)

其中 D(z) 是用垂直於 z 軸的平面 z = z 與 Ω 相交的截面。

52

53

圖 1

註: (i) 這個公式從物理上可以給出這樣的解釋: 把 Ω 看作是一個空間物體, f (x, y, z) 為 物體在 Ω 上的分布密度, (1) 式左端的三重積分即物體的質量, 而 (1) 式右端則表示先把物體 切成薄片, 再把所有薄片的質量累積起來, 故這種方法也稱為 “切片法”。

(ii) 這種方法也適用於垂直 x 軸或 y 軸的平面與 Ω 相截而得到的積分公式。

Z Z Z

Ω

f(x, y, z)dxdydz =

Z

x2

x1

dx

Z Z

D(x)

f(x, y, z)dydz, (2) 其中 Ω = {(x, y, z) | (y, z) ∈ D(x), x1 ≤ x ≤ x²}。

Z Z Z

Ω

f(x, y, z)dxdydz =

Z

y2

y1

dy

Z Z

D(y)

f(x, y, z)dxdz, (3)

其中 Ω = {(x, y, z) | (x, z) ∈ D(y), y1 ≤ y ≤ y2}。

例1. 計算三重積分: I =

^RRR

_Ω(x²+ y²+ z²)dxdydz, 其中 Ω = {(x, y, z) | x²+ y² ≤ 2az, x²+ y²+ z² ≤ 3a²} (a > 0)。

解_{: 解方程組}







x²+ y² = 2az, x²+ y²+ z² = 3a², 得 z = a, z = −3a (捨去)。

平面 z = a 把閉區域 Ω 分成兩部分, 記上半部分為 Ω₁, 下半部分為 Ω2 (圖 2), 在 Ω1, Ω₂ 上分別運用 “先二後一” 法, 則有

I=

Z Z Z

Ω¹

(x²+ y²+ z²)dxdydz +

Z Z Z

Ω²

(x²+ y²+ z²)dxdydz

=

Z

^√3a a

dz

Z Z

D1(z)

(x²+ y²+ z²)dxdy +

Z

a 0

dz

Z Z

D2(z)

(x²+ y²+ z²)dxdy

圖 2

例2. 求曲面 (^x_a²² + ^y_b²² + ^z_c²²)² = ^x_h 所圍立體體積。

解: 顯然 x ≥ 0。 用垂直於 x 軸的平面去截所給立體, 截口曲線為一橢圓:

D(x) : y² b² +z²

c² =

r

x h− x²

a². 因為

^q

^x_h − ^xa²² ≥ 0, 所以 0 ≤ x ≤ (^ah⁴)¹³。

記 x₀ = (^a_h⁴)¹³, λ² =

^q

^x_h −^xa²² (x 6= 0), 則上面的橢圓方程即為:

y²

(λb)² + z²

(λc)² = 1.

此橢圓所圍的面積為 πλ²bc= πbc(

^q

^π_h − ^xa²²), 因此所求立體體積為 V =

Z

x0

0

dx

Z Z

D(x)

dydz = πbc

Z

x0

0

r

x h − x²

a²

!

dx= πbca² h.

若將計算三重積分的 “先二後一” 法推廣到計算四重或四重以上積分, 則可有 “先三後 一”、 “先二後二” 等方法。 例如, 求四維空間球體 x²+ y²+ z²+ t² ≤ a² 的體積, 可用 “先三 後一” 的方法。

55

二 . 等值面 (線) 法

我們以三重積分 I =

^RRR

_Ωf(x, y, z)dxdydz 的計算來說明這種方法。

設 f (x, y, z) = t 是函數 f (x, y, z) 的等值面, 又 {f(x, y, z) ≤ t} ∩ Ω 的體積為 V (t)。

設 f 在 Ω 上的值域為 [α, β], 在 [α, β] 中引進分劃:

α= t0 ≤ t¹ ≤ t² ≤ · · · ≤ tⁿ= β,

則和式 n

X

i=1

t_i[V (ti) − V (tⁱ−1)] (4) 在分劃最大長度 kλk → 0 時的極限為 I。 若 V (t) 為連續可微函數, 則 (4) 式以

^R

α^βtV^′(t)dt 為極限, 故由極限的唯一性知

Z Z Z

Ω

f(x, y, z)dxdydz =

Z

β α

tV^′(t)dt =

Z

β α

tdV(t), (5)

由此可將三重積分化為一重積分。

註: 對於二重積分及三重以上的積分有類似結論; 當被積函數是 f (x, y, z) = t 的函數時 也有類似結論。

例3. 計算積分 I =

^RR

_D(

^q

^x_a +

^q

^y_b)³dxdy, 其中 D = {(x, y) | x ≥ 0, y ≥ 0,

^q

^xa +

q

_y

b ≤ 1}。

解: 設

^q

^x_a +

^q

^y_b = t 為被積函數的等值線, 易知 0 ≤ t ≤ 1。

由

^q

^x_a +

^q

^y_b = t 知 y = b(t² + ^x_a − 2t

^q

^xa), 由 0 ≤

^q

^xa ≤ t 知 0 ≤ x ≤ at²。 等值線

q

_x

a +

^q

^y_b = t 與坐標軸所圍成圖形面積 S(t) = b

Z

at² 0



t² +x a − 2t

r

x a



dx= 1 6abt⁴, 故 dS(t) = ²₃abt³dt, 於是由等值線求積公式知

I = 2 3ab

Z

1

0 t⁶dt= 2 27ab.

例4. 求函數 f (x, y, z) = e

q

x2 a2+^y2

b2+^z2

c2

在橢球 ^x2

a² +^y_b²² + ^z_c²² ≤ 1 內的平均值。

解: 設 Ω = {(x, y, z) | ^xa²² +^y_b²² +^z_c²² ≤ 1}, 則 Ω 的體積 Ω = ⁴3πabc。

設 ^x_a²2 + ^y_b²2 + ^z_c²2 = t 為被積函數的等值面, 則 0 ≤ t ≤ 1。

從而 f (x, y, z) 在 Ω 的平均值為 f(x, y, z) = 1

Ω

Z Z Z

Ω

f(x, y, z)dxdydz = 3(e − 2).

三 . 微元法

從公式 (5) 可以看出, 利用等值面法將三重積分化為一重積分可理解為作變量替換: f (x, y, z) = t, 所以我們除了要知道 t 的取值範圍即新的積分變量的上、 下限外, 還要知道原來的 體積元素 dxdydz 在變量替換下新的體積元素 dV (t)。 但問題是求 dV (t) 並不那麼簡單, 例如 求例 3中的 dS(t)。 但由於 dV (t) (或 dS(t)) 是 V (t) (或 S(t)) 對 t 微分, 因而啟發我們可 直接用微元法求 dV (t) (或 dS(t)) 而不必先求出 V (t) (或 S(t)) 再去對 t 微分。 利用等值面 (線) 法將其他類型的積分 (如曲面積分) 化為一重積分時也可類似處理。

例5. 設 f (t) 為連續函數, 證明:

Z Z

D

f(x − y)dxdy =

Z

A

−A

f(t)(A − |t|)dt, (6)

其中 A > 0 為常數, D = {(x, y) | |x| ≤ ^A₂,|y| ≤ ^A₂}。

證: 令 x − y = t, 則在 D 上, −A ≤ t ≤ A。 如圖3, 直線 y = x 將 D 分成 D¹, D2

兩部分。 在D₁ 上, x − y > 0, 從而 t > 0。

圖 3

57

如圖 3, 設直線 A0B₀ 的方程為 x − y = t, C⁰D₀ 的方程為 x − y = t + dt (dt > 0)。

易知 A₀, B0 的坐標分別為 (^A₂,^A₂ − t), (−^A₂ + t, −^A₂), A0, B0 兩點間的距離為 |A⁰B₀| =

√2(A − t), 坐標原點到直線 A⁰B₀ 的距離為 √^|t|

2 = ^√t₂。

以 |A0B₀| 為長, ^√dt₂ 為寬的矩形面積作為面積元素 dS(t), 即 dS(t) =√

2(A − t) dt

√2 = (A − t)dt, 從而在 D1 上, 有

Z Z

D1

f(x − y)dxdy =

Z

A

0 f(t)(A − t)dt.

同理可知, 在 D2 上有

Z Z

D2

f(x − y)dxdy =

Z

0

−Af(t)(A + t)dt (t < 0),

於是

^RR

_Df(x−y)dxdy =

^RR

D1f(x−y)dxdy +

^RR

D2f(x−y)dxdy =

^R

_−A^A f(t)(A−|t|)dt。

由本題結論不難證明:

設 f (x) 為連續偶函數, 則有

Z Z

D

f(x − y)dxdy = 2

Z

2a

0 (2a − u)f(u)du, (7) 其中 a > 0 為常數, D = {(x, y) | |x| ≤ a, |y| ≤ a}。

例6. 求三重積分 I =

^RRR

_Ω(1+x+y+z)^{dxdydz 3}, 其中 Ω = {(x, y, z) | x + y + z ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0}。

解法1 (等值面法): 設 x + y + z = t 為被積函數的等值面, 易知 0 ≤ t ≤ 1。

取 V (t) = {(x, y, z) | x + y + z ≤ t, x ≥ 0, y ≥ 0, z ≥ 0}, 則 V (t) 的體積為 V(t) = ¹₆t³, 從而 dV (t) = ¹₂t²dt, 所以

I = 1 2

Z

₁

0

t²

(1 + t)³dt= 1 2



ln 2 − 5 8



.

解法2 (微元法): 設 x + y + z = t, 則 0 ≤ t ≤ 1。 被積函數的等值面 x + y + z = t 與 Ω 交成一個邊長為 √

2t 的等邊三角形, 其面積為 ^√₂³t²。 坐標原點到三角形所在平面的距離為

√t

3。 用柱體體積近似代替平面 x + y + z = t, x + y + z = t + dt (dt > 0) 與第一卦限圍成 封閉區域 (臺體) 的體積作為體積元素 dV (t) (積分方向沿矢量 n = {1, 1, 1} 的方向), 即

dV(t) =

√3

2 t²d t

√3

!

= 1 2t²dt,

Schwarz 不等式知

(ax + by + cz)² ≤ (a²+ b²+ c²)(x²+ y²+ z²) = a² + b²+ c², 故 −1 ≤ u ≤ 1。

坐標原點到平面 √^ax+by+cz

a²+b²+c² = u 的距離為 h = |u|。 兩平面 ^√ax+by+cz_a²_+b²_+c² = u 與 ^√ax+by+cz_a2+b²+c² = u+ du (du > 0) 之間的距離為 dh = du, 取它們截

^P

所得球帶面積為面積元素 dS(u) (積 分方向沿矢量 n = ^√_a2 ¹

+b²+c²{a, b, c} 的方向), 即 dS(u) = 2πdu, 從而

Z Z

Σ

f(ax + by + cz)dS = 2π

Z

₁

−1

f(√

a²+ b²+ c²u)du. (8)

例8. 設函數 f (z) 在圓錐面

^P

= {(x, y, z) | z = a√

x²+ y²,0 ≤ z ≤ h} 上連續 (a > 0 為常數), 證明:

Z Z

Σ

f(z)dS = 2π a²

√1 + a²

Z

h 0

zf(z)dz. (9)

證: 取 z = t 為被積函數的等值面, 則 0 ≤ t ≤ h。 考慮圓錐面

^P

夾在平面 z = t 與 z = t + dt (dt > 0) 之間部分的面積。

曲面 z = a√

x²+ y² 可視為 xOz 平面上的直線 z = ax (z ≥ 0) 繞 z 軸旋轉一周而生 成。 由旋轉體側面積公式知所考慮的面積為

2π

Z

t+dt t

z a

s

1 + 1

a²dz = π a²

√1 + a²[2tdt + (dt)²],

故面積元素 dS(t) = ^2π_a²√

1 + a²tdt, 於是

Z Z

Σ

f(z)dS = 2π a²

√1 + a²

Z

h 0

tf(t)dt = 2π a²

√1 + a²

Z

h 0

zf(z)dz. (9) 類似可證:

59

設函數 f (z) 在旋轉拋物面

^P

= {(x, y, z) | z = a(x² + y²), 0 ≤ z ≤ h} 上連續 (a > 0 為常數), 則

Z Z

Σ

f(z)dS = π a

Z

h 0

f(z)√

1 + 4azdz.

參考文獻

1. 謝惠民, 惲自求, 易法槐, 錢定邊, 數學分析習題課講義 (下冊), 北京: 高等教育出版社, 2004。

2. 吳良森, 毛羽輝, 數學分析習題精解 (多變量部分), 北京: 科學出版社, 2004。

3. 同濟大學應用數學系, 微積分 (下冊), 北京: 高等教育出版社, 2000。

4. 沈運彤, 葉牡才, 苑金臣, “先重後單” 法用於求空間曲面所圍立體的體積, 數學學習, 1998 年第 1 期。

5. 馬娜蕊, 對三重積分“先二後一” 計算方法的討論, 高等數學研究, 2000 年第 1 期。

6. 姚允龍, 數學分析, 上海: 復旦大學出版社, 2002。

7. Γ. M. 菲赫金哥爾茨著, 路見可譯, 微積分學教程 (第三卷第一分冊), 北京: 人民教育出版社, 1978。

8. Γ. M. 菲赫金哥爾茨著, 吳親仁, 路見可譯, 微積分學教程 (第三卷第二分冊), 北京: 人民教育出版 社, 1978。

9. B.

^Q

. 吉米多維奇著, 李榮凍譯, 數學分析習題集, 北京: 人民教育出版社, 1978。

—本文作者任教於中國安徽省合肥工業大學數學與信息科學系_—