自在学

我们与你共同进步

  • 分类课程
  • 文章
  • 工作台
  • 订阅

  • 关于我们
  • 隐私政策
  • 使用条款

探索

  • 分类课程
  • 文章
  • 工作台
  • 订阅

网站信息

  • 关于我们
  • 隐私政策
  • 使用条款

加入社区

自在学学习社区微信二维码

微信扫码,交流学习

株洲市自在学教育科技有限公司© 2025 - 2026 版权所有

© 2025 - 2026 株洲市自在学教育科技有限公司 版权所有

湘公网安备43020302000292号|湘ICP备2025148919号-1
分类课程工作台文章订阅
分类课程工作台文章价格

多变量微积分 I

  1. 01多变量微积分研究什么
  2. 02三维坐标、向量与空间几何
  3. 03曲面、二次曲面与坐标变换初识
  4. 04向量值函数与空间曲线
  5. 05多变量函数、定义域、图像与等值线
  6. 06多变量极限与连续
  7. 07偏导数:沿坐标方向看变化
  8. 08切平面、线性近似与全微分
  9. 09多变量链式法则与隐函数求导
  10. 10方向导数与梯度
  11. 11多变量极值与优化
  12. 12约束优化与拉格朗日乘子
  13. 13二重积分:从小矩形累加到曲面下体积
  14. 14一般区域上的二重积分与积分次序
  15. 15极坐标中的二重积分
  16. 16三重积分:空间区域中的累积
  17. 17柱坐标、球坐标与空间坐标变换
  18. 18变量替换与 Jacobian
  19. 19向量场:空间中每一点都有方向
  20. 20线积分:沿路径累加标量和向量
  21. 21保守场、势函数与路径无关
  22. 22Green 定理:平面区域的边界与内部
  23. 23曲面参数化与曲面积分
  24. 24散度、旋度与三维向量场
  25. 25Stokes 定理与 Divergence Theorem
正在加载课程章节内容
课程数学多变量微积分 I变量替换与 Jacobian

变量替换与 Jacobian:把难区域换成好坐标

上一章里,极坐标的面积微元是 r dr dθr\,dr\,d\thetardrdθ,柱坐标和球坐标的体积微元里又分别出现了 rrr 与 ρ2sin⁡ϕ\rho^2\sin\phiρ2sinϕ。第一次见到这些因子时,很容易把它们当成需要单独背诵的规则:换成极坐标就补一个 rrr,换成球坐标就补一个 ρ2sin⁡ϕ\rho^2\sin\phiρ2sinϕ。

可这样记很快就会卡住。只要题目换成一组没见过的新变量,口诀立刻失效。真正通用的问题其实只有一个:新坐标里的一小块,被送回原来的平面或空间后,究竟被拉伸了多少?

Jacobian 回答的就是这个问题。它先把非线性变换在一个点附近看成线性变换,再用行列式计算小面积或小体积的缩放倍数。理解这条主线以后,极坐标、球坐标和一般变量替换就不再是三套互不相干的公式。


坐标一换,小面积也跟着变

设

T(u,v)=(x(u,v),y(u,v))T(u,v)=\bigl(x(u,v),y(u,v)\bigr)T(u,v)=(x(u,v),y(u,v))

把 uvuvuv 平面中的点送到 xyxyxy 平面。我们在 uvuvuv 平面里画一个边长为 dududu、dvdvdv 的小矩形,它的面积是 du dvdu\,dvdudv。老办法最容易犯的错,是把它送到 xyxyxy 平面后,仍然把面积写成 du dvdu\,dvdudv。

问题在于,dududu 和 dvdvdv 只是参数变化量,不一定是目标平面里的实际边长。沿 uuu 方向走一点,像点可能向右走、向上走,也可能两个方向同时走;沿 vvv 方向也是一样。原来的直角小矩形通常会变成一个斜着的、大小也不同的小块。

先看线性变换

(xy)=(abcd)(uv)\begin{pmatrix} x\\ y \end{pmatrix} = \begin{pmatrix} a & b\\ c & d \end{pmatrix} \begin{pmatrix} u\\ v \end{pmatrix}(xy​)=(ac​bd​)(uv​)

当 uuu 增加 dududu 而 vvv 不动时,目标平面里的位移向量是

(ac)du\begin{pmatrix} a\\ c \end{pmatrix}du(ac​)du

当 vvv 增加 dvdvdv 而 uuu 不动时,位移向量是

(bd)dv\begin{pmatrix} b\\ d \end{pmatrix}dv(bd​)dv

这两个向量撑起一个平行四边形。它的面积是

∣ad−bc∣ du dv|ad-bc|\,du\,dv∣ad−bc∣dudv

所以 ∣ad−bc∣|ad-bc|∣ad−bc∣ 就是面积缩放倍数。若它等于 333,新坐标中的每单位面积会在目标平面里铺出面积 333;若它等于 1/41/41/4,面积会被压到原来的四分之一。

下面的交互把这件事直接画出来。可以依次试“放大”“剪切”“翻转”和“压扁”:剪切会改变形状却不一定改变面积,压扁则会让二维小块失去面积。

变量替换不是只换掉公式里的字母。它会同时改变点的坐标、区域的边界以及每个小块的面积。漏掉其中任何一项,得到的都不是原来的积分。


非线性变换为什么仍能用一个矩阵

一般的 T(u,v)T(u,v)T(u,v) 会把直线网格弯成曲线网格。整个平面不可能再由同一个矩阵控制,可如果只盯住点 (u0,v0)(u_0,v_0)(u0​,v0​) 附近极小的一块,弯曲还来不及明显出现,一阶变化已经占了主导。

在这个点附近,两个坐标的增量满足

(ΔxΔy)≈(xuxvyuyv)(u0,v0)(ΔuΔv)\begin{pmatrix} \Delta x\\ \Delta y \end{pmatrix} \approx \begin{pmatrix} x_u & x_v\\ y_u & y_v \end{pmatrix}_{(u_0,v_0)} \begin{pmatrix} \Delta u\\ \Delta v \end{pmatrix}(ΔxΔy​)≈(xu​yu​​xv​yv​​)(u0​,v0​)​(ΔuΔv​)

中间这张矩阵叫 Jacobian 矩阵,也就是变换 TTT 在该点的导数矩阵:

DT(u,v)=(xuxvyuyv)DT(u,v)= \begin{pmatrix} x_u & x_v\\ y_u & y_v \end{pmatrix}DT(u,v)=(xu​yu​​xv​yv​​)

它的第一列 Tu=⟨xu,yu⟩\mathbf T_u=\langle x_u,y_u\rangleTu​=⟨xu​,yu​⟩ 告诉我们“只让 uuu 动”时,像点朝哪里、走多快;第二列 Tv=⟨xv,yv⟩\mathbf T_v=\langle x_v,y_v\rangleTv​=⟨xv​,yv​⟩ 则记录“只让 vvv 动”的效果。于是参数平面中的小矩形,局部近似为由 Tu du\mathbf T_u\,duTu​du 和 Tv dv\mathbf T_v\,dvTv​dv 撑起的平行四边形。

它的有向面积因子是

J(u,v)=det⁡DT(u,v)=det⁡(xuxvyuyv)=xuyv−xvyuJ(u,v)= \det DT(u,v) = \det \begin{pmatrix} x_u & x_v\\ y_u & y_v \end{pmatrix} =x_u y_v-x_v y_uJ(u,v)=detDT(u,v)=det(xu​yu​​xv​yv​​)=xu​yv​−xv​yu​

因此真正的面积是

dA≈∣J(u,v)∣ du dvdA\approx |J(u,v)|\,du\,dvdA≈∣J(u,v)∣dudv

规则小网格经非线性变换后成为弯曲网格,并在一个采样点附近用平行四边形作局部线性近似

图里远处的网格弯曲程度和疏密都在变化;但把某个小块放大以后,它仍可由两条切向量撑起的平行四边形近似。

下面的交互使用

x=u+qsin⁡v,y=v+psin⁡ux=u+q\sin v,\qquad y=v+p\sin ux=u+qsinv,y=v+psinu

拖动采样点时,局部面积因子会随位置变化。调大弯曲参数以后,还能看到 JJJ 变成负数或接近 000 时,局部取向与网格形状怎样改变。

非线性变量替换的核心节奏是:把区域切得足够细,在每个小块上使用各自的线性近似,再把所有局部贡献加起来。Jacobian 可以随位置变化,正是因为不同位置的拉伸程度本来就可能不同。


行列式的正负与绝对值

行列式不只记录面积大小,还记录两个坐标方向的先后关系。

  • J>0J>0J>0 时,uuu 方向转向 vvv 方向的取向在映射后保持不变。
  • J<0J<0J<0 时,映射像镜面反射那样翻转了取向。
  • J=0J=0J=0 时,两条局部方向被压到同一直线上,甚至都被压成零向量,二维面积因此消失。

积分里的普通面积不能为负,所以变量替换公式使用 ∣J∣|J|∣J∣,而不是 JJJ。负号提供的是取向信息,不是“负面积”。

这个区别可以用一个最小例子看清。设

x=u,y=−2vx=u,\qquad y=-2vx=u,y=−2v

则

J=det⁡(100−2)=−2J= \det \begin{pmatrix} 1 & 0\\ 0 & -2 \end{pmatrix} =-2J=det(10​0−2​)=−2

这个变换把图形上下翻转,同时把竖直方向拉长为原来的 222 倍。面积缩放因子显然是 222,所以 dA=2 du dvdA=2\,du\,dvdA=2dudv。

看到 J<0J<0J<0 后把积分整体乘成负数,是把“取向”与“面积”混在了一起。计算面积、质量和普通二重积分时始终使用 ∣J∣|J|∣J∣;以后处理带方向的曲线或曲面时,符号才会以别的方式重新出现。


二重积分的变量替换公式

设参数平面中的区域 SSS 经 TTT 映到 xyxyxy 平面中的区域 RRR。在 TTT 连续可微、区域内部基本上一一对应,并且除去面积为 000 的少量点外 JJJ 不为 000 时,变量替换公式是

∬Rf(x,y) dA=∬Sf(x(u,v),y(u,v))∣∂(x,y)∂(u,v)∣du dv\iint_R f(x,y)\,dA = \iint_S f\bigl(x(u,v),y(u,v)\bigr) \left| \frac{\partial(x,y)}{\partial(u,v)} \right| du\,dv∬R​f(x,y)dA=∬S​f(x(u,v),y(u,v))​∂(u,v)∂(x,y)​​dudv

这条公式看起来长,实际只要求完成三笔“对账”:

把被积函数中的 x,yx,yx,y 全部写成 u,vu,vu,v。若还剩下无法解释的 xxx 或 yyy,换元就没有完成。

把原区域 RRR 的边界翻译到 uvuvuv 平面,得到新区域 SSS。积分限应描述 SSS,不能继续沿用原区域的数值。

用 ∣∂(x,y)/∂(u,v)∣du dv\left|\partial(x,y)/\partial(u,v)\right|du\,dv∣∂(x,y)/∂(u,v)∣dudv 替换 dAdAdA。这里的分子必须是原积分中的坐标,分母才是新变量。

正向 Jacobian 与反向 Jacobian

题目有时直接给出 x=x(u,v)x=x(u,v)x=x(u,v)、y=y(u,v)y=y(u,v)y=y(u,v)。另一些题目却先给 u=u(x,y)u=u(x,y)u=u(x,y)、v=v(x,y)v=v(x,y)v=v(x,y)。如果反变换存在,两个 Jacobian 矩阵互为逆矩阵:

D(T−1)(T(u,v))=(DT(u,v))−1D(T^{-1})\bigl(T(u,v)\bigr)=\bigl(DT(u,v)\bigr)^{-1}D(T−1)(T(u,v))=(DT(u,v))−1

取行列式后得到

∂(x,y)∂(u,v)=1∂(u,v)∂(x,y)\frac{\partial(x,y)}{\partial(u,v)} = \frac{1}{\dfrac{\partial(u,v)}{\partial(x,y)}}∂(u,v)∂(x,y)​=∂(x,y)∂(u,v)​1​

这里默认两边在对应点计算,而且分母不为 000。因此

dA=1∣∂(u,v)∂(x,y)∣du dvdA = \frac{1}{\left|\dfrac{\partial(u,v)}{\partial(x,y)}\right|} du\,dvdA=​∂(x,y)∂(u,v)​​1​dudv

多元情形里,不能把某一个偏导数单独倒过来。例如 ∂r/∂x\partial r/\partial x∂r/∂x 通常不等于 1/(∂x/∂r)1/(\partial x/\partial r)1/(∂x/∂r),因为两边保持不动的变量不同。可靠的倒数关系发生在完整的 Jacobian 矩阵及其行列式之间。


先把边界变简单,再决定积分限

变量替换真正省力的地方,往往不在被积函数,而在区域。一个好的新变量通常直接来自边界方程:若两组边分别写成 g(x,y)=cg(x,y)=cg(x,y)=c 与 h(x,y)=dh(x,y)=dh(x,y)=d,就优先试

u=g(x,y),v=h(x,y)u=g(x,y),\qquad v=h(x,y)u=g(x,y),v=h(x,y)

这样原来斜着或弯着的四条边,到了 uvuvuv 平面里可能只是 u=u0u=u_0u=u0​、u=u1u=u_1u=u1​、v=v0v=v_0v=v0​、v=v1v=v_1v=v1​,新区域自然就是矩形。

定限时可以按下面的顺序走。

先列出原区域的每一条边界,并观察哪些表达式在同一组边上保持常数。不要急着算 Jacobian;变量选得不好,后面算得再熟也省不了事。

把边界代入新变量,写出它们对应的是 uuu 的边还是 vvv 的边。若有交点,把交点也映过去,检查新区域的角点是否接得上。

选原区域内部的一个测试点,确认不等号方向。只变边界方程、不检查区域落在哪一侧,很容易把内部和外部写反。

确认新区域没有重叠计数或漏掉一块,再选择积分顺序。若 SSS 是矩形,两个顺序通常都可以;若仍有曲边,就选上下界更简单的顺序。

例题:把平移椭圆拉回单位圆

设

R:(x−1)29+(y+2)24≤1R:\quad \frac{(x-1)^2}{9}+\frac{(y+2)^2}{4}\le 1R:9(x−1)2​+4(y+2)2​≤1

求

∬R(x−2y) dA\iint_R (x-2y)\,dA∬R​(x−2y)dA

椭圆难写积分限,是因为它相对单位圆做了平移和不等比例拉伸。把这些动作倒过来,令

x=1+3u,y=−2+2vx=1+3u,\qquad y=-2+2vx=1+3u,y=−2+2v

原区域立即变成

S:u2+v2≤1S:\quad u^2+v^2\le 1S:u2+v2≤1

中心在一逗号负二、半轴长为三和二的椭圆经平移缩放后对应单位圆

先算面积因子:

J=det⁡(3002)=6J= \det \begin{pmatrix} 3 & 0\\ 0 & 2 \end{pmatrix} =6J=det(30​02​)=6

所以 dA=6 du dvdA=6\,du\,dvdA=6dudv。

再换被积函数:

x−2y=(1+3u)−2(−2+2v)=5+3u−4vx-2y=(1+3u)-2(-2+2v)=5+3u-4vx−2y=(1+3u)−2(−2+2v)=5+3u−4v

写出新积分:

∬R(x−2y) dA=6∬S(5+3u−4v) du dv\iint_R (x-2y)\,dA = 6\iint_S(5+3u-4v)\,du\,dv∬R​(x−2y)dA=6∬S​(5+3u−4v)dudv

单位圆关于两条坐标轴对称,所以 uuu 与 vvv 的积分都为 000。只剩常数项:

6⋅5⋅area⁡(S)=30π6\cdot 5\cdot \operatorname{area}(S)=30\pi6⋅5⋅area(S)=30π

这个结果还能反过来检查。椭圆面积是 6π6\pi6π,而线性函数 x−2yx-2yx−2y 在关于中心对称的区域上的平均值,等于它在中心 (1,−2)(1,-2)(1,−2) 的值 555。所以积分应为 5⋅6π=30π5\cdot 6\pi=30\pi5⋅6π=30π。

例题:把斜平行四边形摊成矩形

设 RRR 由

x−y=0,x−y=2,x+y=1,x+y=4x-y=0,\qquad x-y=2,\qquad x+y=1,\qquad x+y=4x−y=0,x−y=2,x+y=1,x+y=4

围成。求

∬R(x+y) dA\iint_R(x+y)\,dA∬R​(x+y)dA

四条边已经把新变量写在脸上了。令

u=x−y,v=x+yu=x-y,\qquad v=x+yu=x−y,v=x+y

那么新区域直接变成

0≤u≤2,1≤v≤40\le u\le 2,\qquad 1\le v\le 40≤u≤2,1≤v≤4

由 u=x−yu=x-yu=x−y、v=x+yv=x+yv=x+y 解出反变换:

x=u+v2,y=v−u2x=\frac{u+v}{2},\qquad y=\frac{v-u}{2}x=2u+v​,y=2v−u​

计算原坐标对新变量的 Jacobian:

∂(x,y)∂(u,v)=det⁡(1212−1212)=12\frac{\partial(x,y)}{\partial(u,v)} = \det \begin{pmatrix} \dfrac12 & \dfrac12\\ -\dfrac12 & \dfrac12 \end{pmatrix} =\frac12∂(u,v)∂(x,y)​=det​21​−21​​21​21​​​=21​

因此 dA=12 du dvdA=\dfrac12\,du\,dvdA=21​dudv。

被积函数 x+yx+yx+y 正好就是 vvv,所以

∬R(x+y) dA=∫02∫14v2 dv du\iint_R(x+y)\,dA = \int_0^2\int_1^4 \frac{v}{2}\,dv\,du∬R​(x+y)dA=∫02​∫14​2v​dvdu

完成计算:

∫02∫14v2 dv du=∫02154 du=152\int_0^2\int_1^4 \frac{v}{2}\,dv\,du = \int_0^2\frac{15}{4}\,du =\frac{15}{2}∫02​∫14​2v​dvdu=∫02​415​du=215​

例题:非线性边界也能变成矩形

设 RRR 位于第一象限,由两条双曲线 xy=1xy=1xy=1、xy=4xy=4xy=4 以及两条射线 y=xy=xy=x、y=9xy=9xy=9x 围成。求 RRR 的面积。

令

u=xy,v=yxu=xy,\qquad v=\frac{y}{x}u=xy,v=xy​

新区域就是

1≤u≤4,1≤v≤91\le u\le 4,\qquad 1\le v\le 91≤u≤4,1≤v≤9

在第一象限内,反变换唯一写成

x=uv,y=uvx=\sqrt{\frac{u}{v}},\qquad y=\sqrt{uv}x=vu​​,y=uv​

直接对反变换求导当然可以,但这里算正向 Jacobian 更快:

∂(u,v)∂(x,y)=det⁡(yx−yx21x)=2yx=2v\frac{\partial(u,v)}{\partial(x,y)} = \det \begin{pmatrix} y & x\\ -\dfrac{y}{x^2} & \dfrac1x \end{pmatrix} =2\frac{y}{x}=2v∂(x,y)∂(u,v)​=det(y−x2y​​xx1​​)=2xy​=2v

于是

dA=12v du dvdA=\frac{1}{2v}\,du\,dvdA=2v1​dudv

面积为

area⁡(R)=∫14∫1912v dv du=3ln⁡3\operatorname{area}(R) = \int_1^4\int_1^9\frac{1}{2v}\,dv\,du =3\ln 3area(R)=∫14​∫19​2v1​dvdu=3ln3

这个例子说明,变量替换不只会旋转或缩放直线边界。只要新变量顺着边界的“家族”来选,双曲线与射线也能被整理成横平竖直的矩形。

下面的交互把椭圆、斜平行四边形和极坐标扇形放在一起比较。观察重点是每组边界到了新坐标里分别变成了什么常数。


极坐标为什么多出一个 r

现在回头看最熟悉的变换:

x=rcos⁡θ,y=rsin⁡θx=r\cos\theta,\qquad y=r\sin\thetax=rcosθ,y=rsinθ

它的 Jacobian 矩阵是

DT(r,θ)=(cos⁡θ−rsin⁡θsin⁡θrcos⁡θ)DT(r,\theta)= \begin{pmatrix} \cos\theta & -r\sin\theta\\ \sin\theta & r\cos\theta \end{pmatrix}DT(r,θ)=(cosθsinθ​−rsinθrcosθ​)

两列分别代表径向小位移和角向小位移。第一列长度是 111,第二列长度是 rrr,而且两列互相垂直。因此

J=det⁡DT(r,θ)=rJ=\det DT(r,\theta)=rJ=detDT(r,θ)=r

从而

dA=r dr dθdA=r\,dr\,d\thetadA=rdrdθ

这正好与几何画面一致:drdrdr 是径向厚度,r dθr\,d\thetardθ 才是角向弧长。Jacobian 没有发明一个额外的 rrr,它只是把这两条局部边的面积用行列式记了下来。

极坐标也提醒我们注意两个条件。第一,r=0r=0r=0 时 J=0J=0J=0,所有角度都落到同一个原点;第二,θ\thetaθ 相差 2π2\pi2π 的参数表示同一点。做完整圆盘积分时,我们通常取一个长度为 2π2\pi2π 的角度区间,并允许角度区间的两条边在目标平面重合。重复只发生在边界或原点这类面积为 000 的集合上,不会改变普通二重积分。


一一对应与可逆性不能只看 J

J(u0,v0)≠0J(u_0,v_0)\ne 0J(u0​,v0​)=0 能保证变换在这个点附近没有把二维小块压扁,并且局部可以反解。这里的关键词是“附近”。它并不自动保证整个大区域上都一一对应。

例如

x=u2−v2,y=2uvx=u^2-v^2,\qquad y=2uvx=u2−v2,y=2uv

有

J=4(u2+v2)J=4(u^2+v^2)J=4(u2+v2)

除原点外,JJJ 都不为 000。可 (u,v)(u,v)(u,v) 与 (−u,−v)(-u,-v)(−u,−v) 总会映到同一个 (x,y)(x,y)(x,y),所以在包含这两个点的大区域上,变换仍然会重复覆盖。

使用变量替换公式前,要分开检查两件事:

  • 局部是否退化:看 JJJ 是否在区域内部变成 000。
  • 全局是否重复:看不同参数点是否会落到同一个目标点。

若只是边界上有少量重合,例如极坐标角度区间的两条边,通常不影响积分。若区域内部被覆盖两次,就必须缩小参数范围、把区域拆成几块,或者在明确知道覆盖次数时作相应修正。

“Jacobian 处处不为零”是局部可逆条件,不是全局一一对应的保证。变量替换题中,画出几条坐标曲线、检查角点与内部测试点,往往比只盯着行列式更能发现重复覆盖。


三重积分中的 3×3 Jacobian

空间里的逻辑没有变化,只是小平行四边形升级成了小平行六面体。设

T(u,v,w)=(x(u,v,w),y(u,v,w),z(u,v,w))T(u,v,w)=\bigl(x(u,v,w),y(u,v,w),z(u,v,w)\bigr)T(u,v,w)=(x(u,v,w),y(u,v,w),z(u,v,w))

它的 Jacobian 矩阵是

DT=(xuxvxwyuyvywzuzvzw)DT= \begin{pmatrix} x_u & x_v & x_w\\ y_u & y_v & y_w\\ z_u & z_v & z_w \end{pmatrix}DT=​xu​yu​zu​​xv​yv​zv​​xw​yw​zw​​​

三列向量 Tu\mathbf T_uTu​、Tv\mathbf T_vTv​、Tw\mathbf T_wTw​ 是参数小长方体三条边在真实空间里的局部像。它们的标量三重积

J(u,v,w)=det⁡(xuxvxwyuyvywzuzvzw)J(u,v,w)= \det \begin{pmatrix} x_u & x_v & x_w\\ y_u & y_v & y_w\\ z_u & z_v & z_w \end{pmatrix}J(u,v,w)=det​xu​yu​zu​​xv​yv​zv​​xw​yw​zw​​​

给出有向体积缩放,普通体积则取绝对值:

dV=∣J(u,v,w)∣ du dv dwdV=|J(u,v,w)|\,du\,dv\,dwdV=∣J(u,v,w)∣dudvdw

因此三重积分的换元公式是

∭Ef(x,y,z) dV=∭Gf(x(u,v,w),y(u,v,w),z(u,v,w))∣J(u,v,w)∣ du dv dw\iiint_E f(x,y,z)\,dV = \iiint_G f\bigl(x(u,v,w),y(u,v,w),z(u,v,w)\bigr) |J(u,v,w)|\,du\,dv\,dw∭E​f(x,y,z)dV=∭G​f(x(u,v,w),y(u,v,w),z(u,v,w))∣J(u,v,w)∣dudvdw

球坐标因子也来自同一个行列式

球坐标变换为

x=ρsin⁡ϕcos⁡θ,y=ρsin⁡ϕsin⁡θ,z=ρcos⁡ϕx=\rho\sin\phi\cos\theta,\qquad y=\rho\sin\phi\sin\theta,\qquad z=\rho\cos\phix=ρsinϕcosθ,y=ρsinϕsinθ,z=ρcosϕ

对应的 Jacobian 矩阵是

(sin⁡ϕcos⁡θρcos⁡ϕcos⁡θ−ρsin⁡ϕsin⁡θsin⁡ϕsin⁡θρcos⁡ϕsin⁡θρsin⁡ϕcos⁡θcos⁡ϕ−ρsin⁡ϕ0)\begin{pmatrix} \sin\phi\cos\theta & \rho\cos\phi\cos\theta & -\rho\sin\phi\sin\theta\\ \sin\phi\sin\theta & \rho\cos\phi\sin\theta & \rho\sin\phi\cos\theta\\ \cos\phi & -\rho\sin\phi & 0 \end{pmatrix}​sinϕcosθsinϕsinθcosϕ​ρcosϕcosθρcosϕsinθ−ρsinϕ​−ρsinϕsinθρsinϕcosθ0​​

计算行列式得到

J=ρ2sin⁡ϕJ=\rho^2\sin\phiJ=ρ2sinϕ

在 0≤ϕ≤π0\le\phi\le\pi0≤ϕ≤π 的常用范围内,sin⁡ϕ≥0\sin\phi\ge 0sinϕ≥0,于是

dV=ρ2sin⁡ϕ dρ dϕ dθdV=\rho^2\sin\phi\,d\rho\,d\phi\,d\thetadV=ρ2sinϕdρdϕdθ

柱坐标同样是三维变换 x=rcos⁡θx=r\cos\thetax=rcosθ、y=rsin⁡θy=r\sin\thetay=rsinθ、z=zz=zz=z。第三个方向没有缩放,前两个方向仍贡献极坐标的 rrr,所以 3×3 行列式就是 rrr。

例题:线性空间变换下的积分

令单位立方体 G:0≤u,v,w≤1G:0\le u,v,w\le 1G:0≤u,v,w≤1 通过

x=u+v,y=v+w,z=u+wx=u+v,\qquad y=v+w,\qquad z=u+wx=u+v,y=v+w,z=u+w

映到空间区域 EEE。求

∭E(x+y+z) dV\iiint_E(x+y+z)\,dV∭E​(x+y+z)dV

先算 3×3 Jacobian:

J=det⁡(110011101)=2J= \det \begin{pmatrix} 1 & 1 & 0\\ 0 & 1 & 1\\ 1 & 0 & 1 \end{pmatrix} =2J=det​101​110​011​​=2

所以单位立方体被送成体积为 222 的平行六面体。

被积函数变成 x+y+z=2(u+v+w)x+y+z=2(u+v+w)x+y+z=2(u+v+w),再乘体积因子 222:

∭E(x+y+z) dV=∫01∫01∫014(u+v+w) du dv dw\iiint_E(x+y+z)\,dV = \int_0^1\int_0^1\int_0^1 4(u+v+w)\,du\,dv\,dw∭E​(x+y+z)dV=∫01​∫01​∫01​4(u+v+w)dudvdw

单位立方体上 u,v,wu,v,wu,v,w 的平均值都为 1/21/21/2,所以积分为

4(12+12+12)=64\left(\frac12+\frac12+\frac12\right)=64(21​+21​+21​)=6

验算也很直接:EEE 的体积是 222,而 x+y+zx+y+zx+y+z 在 EEE 上的平均值是 333,所以积分应为 2⋅3=62\cdot 3=62⋅3=6。


定限与计算后的验算

变量替换题的错误往往不会在最后一步自己暴露出来。最好在写完新积分后停一下,用几项成本很低的检查把它拦住。

写积分前的检查

  • 看边界:原区域的每条边是否都在新区域中找到对应边界?
  • 看内部:至少映射一个内部测试点,确认不等号方向和区域位置。
  • 看覆盖:新区域是否把原区域恰好覆盖一次?
  • 看方向:需要的是 ∂(x,y)/∂(u,v)\partial(x,y)/\partial(u,v)∂(x,y)/∂(u,v),还是它的倒数?
  • 看绝对值:面积与体积因子是否已经取绝对值?
  • 看变量:被积函数、微元和积分限是否都只含新变量?

算出结果后的检查

第一种检查是量纲。若 x,yx,yx,y 有长度单位,那么 dAdAdA 必须有长度平方单位。极坐标中的 drdrdr 有长度单位,dθd\thetadθ 没有,因此必须由 rrr 补出另一个长度。

第二种检查是先令 f=1f=1f=1。变量替换公式应退化成区域面积或体积。如果连常数函数都算不出合理的几何大小,通常是 Jacobian 或积分限出了问题。

第三种检查是估计数量级。若区域面积约为 101010,且被积函数在区域里始终介于 222 与 555,积分就应在 202020 与 505050 之间。精确结果跑出这个范围,无须重算全部步骤就知道有错。

第四种检查是利用对称性。关于坐标轴对称的区域上,奇函数项应当抵消;中心对称区域上线性函数的平均值等于中心处的函数值。

第五种检查是正反变换。若两个方向的 Jacobian 都算得出来,它们的乘积应为 111。线性变换尤其适合用这一点快速核对。


常见误区

只换微元,不换函数

把 dAdAdA 换成 ∣J∣du dv|J|du\,dv∣J∣dudv 后,被积函数中的 x,yx,yx,y 也必须换掉。新积分里同时混着 x,y,u,vx,y,u,vx,y,u,v,通常意味着工作只做了一半。

把 Jacobian 的方向写反

从 dx dydx\,dydxdy 换到 du dvdu\,dvdudv,需要的是

∣∂(x,y)∂(u,v)∣\left|\frac{\partial(x,y)}{\partial(u,v)}\right|​∂(u,v)∂(x,y)​​

若手里算的是 ∂(u,v)/∂(x,y)\partial(u,v)/\partial(x,y)∂(u,v)/∂(x,y),要在可逆处取倒数。

忘记绝对值

负行列式表示取向翻转。面积和体积仍为正,普通二重、三重积分使用绝对值。

只变边界方程,不检查区域内部

边界变成矩形的四条边,并不自动告诉你矩形内部对应原区域内部。代入一个测试点,才能确认不等号方向。

以为 J 不为零就一定一一对应

J≠0J\ne 0J=0 只给出局部可逆性。周期角度、平方映射和其他折叠变换都可能让远处的不同参数点落到同一点。

把 J 等于零一概判成不能积分

若 J=0J=0J=0 只发生在边界或孤立点上,例如极坐标的原点,常规积分通常仍可进行。真正危险的是区域内部出现整片折叠或退化,并造成无法控制的重复覆盖。


练习

练习 1:设 x=2u+vx=2u+vx=2u+v、y=u−vy=u-vy=u−v。求 Jacobian、面积缩放因子,并说明取向是否翻转。

Jacobian 为

J=det⁡(211−1)=−3J= \det \begin{pmatrix} 2 & 1\\ 1 & -1 \end{pmatrix} =-3J=det(21​1−1​)=−3

因此面积缩放因子是 333,变换把面积放大为原来的 333 倍。因为 J<0J<0J<0,取向发生翻转。

练习 2:用变量替换求椭圆 x2/4+y2/25≤1x^2/4+y^2/25\le 1x2/4+y2/25≤1 的面积。

令 x=2ux=2ux=2u、y=5vy=5vy=5v,椭圆变成单位圆 u2+v2≤1u^2+v^2\le 1u2+v2≤1。Jacobian 为

J=det⁡(2005)=10J= \det \begin{pmatrix} 2 & 0\\ 0 & 5 \end{pmatrix} =10J=det(20​05​)=10

所以

area⁡(R)=∬u2+v2≤110 du dv=10π\operatorname{area}(R) = \iint_{u^2+v^2\le 1}10\,du\,dv =10\piarea(R)=∬u2+v2≤1​10dudv=10π

长半轴乘短半轴再乘 π\piπ 也给出 5⋅2⋅π=10π5\cdot2\cdot\pi=10\pi5⋅2⋅π=10π。

练习 3:区域 RRR 由 x+y=0x+y=0x+y=0、x+y=3x+y=3x+y=3、x−2y=1x-2y=1x−2y=1、x−2y=5x-2y=5x−2y=5 围成。令 u=x+yu=x+yu=x+y、v=x−2yv=x-2yv=x−2y,写出新区域,并求 dAdAdA 与 du dvdu\,dvdudv 的关系。

新区域是

0≤u≤3,1≤v≤50\le u\le 3,\qquad 1\le v\le 50≤u≤3,1≤v≤5

反变换为

x=2u+v3,y=u−v3x=\frac{2u+v}{3},\qquad y=\frac{u-v}{3}x=32u+v​,y=3u−v​

因此

∂(x,y)∂(u,v)=det⁡(231313−13)=−13\frac{\partial(x,y)}{\partial(u,v)} = \det \begin{pmatrix} \dfrac23 & \dfrac13\\ \dfrac13 & -\dfrac13 \end{pmatrix} =-\frac13∂(u,v)∂(x,y)​=det​32​31​​31​−31​​​=−31​

面积微元取绝对值,所以 dA=13 du dvdA=\dfrac13\,du\,dvdA=31​dudv。

练习 4:设第一象限区域由 xy=2xy=2xy=2、xy=6xy=6xy=6、y=2xy=2xy=2x、y=8xy=8xy=8x 围成。使用 u=xyu=xyu=xy、v=y/xv=y/xv=y/x 求区域面积。

新区域是 2≤u≤62\le u\le 62≤u≤6、2≤v≤82\le v\le 82≤v≤8。与正文例题相同,∂(u,v)/∂(x,y)=2v\partial(u,v)/\partial(x,y)=2v∂(u,v)/∂(x,y)=2v,所以

dA=12v du dvdA=\frac{1}{2v}\,du\,dvdA=2v1​dudv

区域面积为

∫26∫2812v dv du=4ln⁡2\int_2^6\int_2^8\frac{1}{2v}\,dv\,du =4\ln 2∫26​∫28​2v1​dvdu=4ln2

练习 5:直接计算极坐标变换的 Jacobian,并解释为什么 r=0r=0r=0 处与角度区间需要额外小心。

Jacobian 为

J=det⁡(cos⁡θ−rsin⁡θsin⁡θrcos⁡θ)=rJ= \det \begin{pmatrix} \cos\theta & -r\sin\theta\\ \sin\theta & r\cos\theta \end{pmatrix} =rJ=det(cosθsinθ​−rsinθrcosθ​)=r

当 r=0r=0r=0 时,所有角度都对应同一个原点,变换在那里不是一一对应,Jacobian 也为 000。另外,θ\thetaθ 与 θ+2π\theta+2\piθ+2π 表示同一方向,因此角度区间若取得过长,会重复覆盖。常规圆盘积分把角度限制在任意一个长度为 2π2\pi2π 的区间。

练习 6:设 x=u2−v2x=u^2-v^2x=u2−v2、y=2uvy=2uvy=2uv。求 Jacobian,并判断能否在圆环 1≤u2+v2≤41\le u^2+v^2\le 41≤u2+v2≤4 上直接把它当成一一对应变换。

Jacobian 为

J=det⁡(2u−2v2v2u)=4(u2+v2)J= \det \begin{pmatrix} 2u & -2v\\ 2v & 2u \end{pmatrix} =4(u^2+v^2)J=det(2u2v​−2v2u​)=4(u2+v2)

它在这个圆环上始终为正,但变换仍不是一一对应,因为 T(u,v)=T(−u,−v)T(u,v)=T(-u,-v)T(u,v)=T(−u,−v)。圆环内部的点会成对映到同一个目标点,直接套公式会重复计数。可以把参数角度限制在长度为 π\piπ 的区间,再处理边界重合。

练习 7:设单位立方体经 x=2ux=2ux=2u、y=u+vy=u+vy=u+v、z=v+3wz=v+3wz=v+3w 映到区域 EEE。求 EEE 的体积。

Jacobian 为

J=det⁡(200110013)=6J= \det \begin{pmatrix} 2 & 0 & 0\\ 1 & 1 & 0\\ 0 & 1 & 3 \end{pmatrix} =6J=det​210​011​003​​=6

参数区域的体积是 111,所以

volume⁡(E)=∭[0,1]36 du dv dw=6\operatorname{volume}(E) = \iiint_{[0,1]^3}6\,du\,dv\,dw =6volume(E)=∭[0,1]3​6dudvdw=6

练习 8:有人把 u=x+yu=x+yu=x+y、v=x−yv=x-yv=x−y 的换元写成

dx dy=∣∂(u,v)∂(x,y)∣du dvdx\,dy= \left| \frac{\partial(u,v)}{\partial(x,y)} \right|du\,dvdxdy=​∂(x,y)∂(u,v)​​dudv

指出错误并写出正确结果。

他把 Jacobian 的方向写反了。先算

∂(u,v)∂(x,y)=det⁡(111−1)=−2\frac{\partial(u,v)}{\partial(x,y)} = \det \begin{pmatrix} 1 & 1\\ 1 & -1 \end{pmatrix} =-2∂(x,y)∂(u,v)​=det(11​1−1​)=−2

所以 ∂(x,y)/∂(u,v)=−1/2\partial(x,y)/\partial(u,v)=-1/2∂(x,y)/∂(u,v)=−1/2,面积微元为

dx dy=12 du dvdx\,dy=\frac12\,du\,dvdxdy=21​dudv

从坐标变换走向向量场

到这里,Jacobian 已经不只是积分前面那个“补上的因子”了。它是一张局部说明书:给它一个很小的参数位移,它会告诉你真实平面或空间里的位移向量;取它的行列式,才进一步得到面积或体积怎样缩放。

下一章研究向量场时,我们会把注意力从“一个小块被坐标变换送到哪里”转向“空间中每一点本身带着怎样的向量”。不过这里的两列或三列偏导向量不会消失。它们已经让我们习惯了一件事:在多元世界里,局部变化不再由一个斜率说完,而要由一组方向共同描述。

上一章柱坐标、球坐标与空间坐标变换下一章向量场:空间中每一点都有方向