自在学

我们与你共同进步

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

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

探索

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

网站信息

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

加入社区

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

微信扫码,交流学习

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

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

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

高级物理三

  1. 01矢量分析与曲线坐标系
  2. 02质点系与刚体力学
  3. 03拉格朗日力学
  4. 04哈密顿力学
  5. 05振动与波的深入分析
  6. 06热力学深化
  7. 07统计物理基础
  8. 08量子统计简介
  9. 09连续介质力学入门
  10. 10非线性动力学与混沌入门
  11. 11对称性与守恒律
  12. 12连续介质的波动与声学
正在加载课程章节内容
课程物理高级物理三矢量分析与曲线坐标系

矢量分析与曲线坐标系

许多物理量带有方向,速度、力、电场强度都是如此。流体的温度、气体的压强则只需一个数就能描述。

把这两种量区分清楚,再用统一的记号写出梯度、散度、旋度,就能在直角坐标、柱坐标和球坐标之间切换,并把高斯定理、斯托克斯定理写成紧凑的矢量形式。下面从场与矢量的基本约定出发,逐步过渡到这些工具及其在简单物理问题中的用法。

welearn-20441038.png

书写习惯上,位置用 r⃗=xx^+yy^+zz^\vec{r} = x\hat{x} + y\hat{y} + z\hat{z}r=xx^+yy^​+zz^ 表示,单位矢量 x^\hat{x}x^、y^\hat{y}y^​、z^\hat{z}z^ 相互垂直且长度为 111。场量在一点处的取值随该点坐标变化而变化,因此求导时始终把坐标当作自变量,把场量当作因变量,后续梯度、散度、旋度的公式都建立在这一套约定之上。

同一物理对象可以在不同坐标下写出不同形式的方程,但几何关系与守恒律不变。熟练之后,往往先在直角坐标里核对公式与量纲,再按对称性换到柱坐标或球坐标做积分或解微分方程,既减少分量个数,又突出对称性带来的简化。

记号上始终区分“对坐标的偏导”与“对时间的导数”;下面出现的 ∂/∂x\partial/\partial x∂/∂x 等均在空间固定、对坐标求导。时间依赖场在后面振动与波、电磁感应等内容中再引入对 ttt 的偏导,与空间导数配合成完整的场方程。


标量场、矢量场与常用运算

空间中每一点对应一个数的映射称为标量场,记为 ϕ(x,y,z)\phi(x,y,z)ϕ(x,y,z) 或 ϕ(r⃗)\phi(\vec{r})ϕ(r)。每一点对应一个矢量的映射称为矢量场,记为 A⃗(x,y,z)\vec{A}(x,y,z)A(x,y,z)。电场强度 E⃗\vec{E}E、磁感应强度 B⃗\vec{B}B、流速 v⃗\vec{v}v 都是矢量场;温度 TTT、电势 φ\varphiφ 是标量场。

标量场在一点附近的变化用偏导数描述即可;矢量场有三个分量,每个分量都可以对 xxx、yyy、zzz 求偏导,因此需要更紧凑的符号把多种导数信息打包成梯度、散度、旋度。在此之前,点积与叉积是构造这些算符的基础。

两个矢量 a⃗\vec{a}a、b⃗\vec{b}b 的点积 a⃗⋅b⃗=∣a⃗∣∣b⃗∣cos⁡θ\vec{a} \cdot \vec{b} = |\vec{a}||\vec{b}|\cos\thetaa⋅b=∣a∣∣b∣cosθ 给出投影关系,结果为标量。叉积 a⃗×b⃗\vec{a} \times \vec{b}a×b 的大小为 ∣a⃗∣∣b⃗∣sin⁡θ|\vec{a}||\vec{b}|\sin\theta∣a∣∣b∣sinθ,方向按右手螺旋由 a⃗\vec{a}a 转向 b⃗\vec{b}b 确定,结果垂直于两者所在平面。点积为零表示垂直,叉积为零表示平行(或有一矢量为零)。

矢量模长 ∣a⃗∣=ax2+ay2+az2|\vec{a}| = \sqrt{a_x^2 + a_y^2 + a_z^2}∣a∣=ax2​+ay2​+az2​​ 与分量选取无关,换坐标系时模长不变。点积与叉积的公式在直角坐标中最易计算;换到曲线坐标时,基矢量可能随位置变化,需要区分对分量的导数与对几何方向的导数,因此后续专门给出 ∇\nabla∇、∇⋅\nabla \cdot∇⋅、∇×\nabla \times∇× 在柱、球坐标中的标准表,避免在弯曲坐标下误用直角分量公式。

例1: 设 a⃗=(1,0,0)\vec{a} = (1,0,0)a=(1,0,0),b⃗=(1,1,0)\vec{b} = (1,1,0)b=(1,1,0)(单位均取 m\text{m}m 仅表示分量无量纲化书写)。点积用分量形式计算:a⃗⋅b⃗=1×1+0×1+0×0=1\vec{a} \cdot \vec{b} = 1 \times 1 + 0 \times 1 + 0 \times 0 = 1a⋅b=1×1+0×1+0×0=1,与公式 ∣a⃗∣∣b⃗∣cos⁡θ|\vec{a}||\vec{b}|\cos\theta∣a∣∣b∣cosθ 一致,此时 ∣a⃗∣=1|\vec{a}|=1∣a∣=1,∣b⃗∣=2|\vec{b}|=\sqrt{2}∣b∣=2​,cos⁡θ=1/2\cos\theta = 1/\sqrt{2}cosθ=1/2​。叉积用行列式或分量公式得到 a⃗×b⃗=(0,0,1)\vec{a} \times \vec{b} = (0,0,1)a×b=(0,0,1)。点积刻画 b⃗\vec{b}b 在 xxx 轴上的投影长度,叉积给出垂直于 xyxyxy 平面的单位方向,在描述力矩、洛伦兹力方向时常用。

分量公式 (ax,ay,az)×(bx,by,bz)=(aybz−azby, azbx−axbz, axby−aybx)(a_x,a_y,a_z) \times (b_x,b_y,b_z) = (a_y b_z - a_z b_y,\ a_z b_x - a_x b_z,\ a_x b_y - a_y b_x)(ax​,ay​,az​)×(bx​,by​,bz​)=(ay​bz​−az​by​, az​bx​−ax​bz​, ax​by​−ay​bx​) 与行列式展开等价,计算时任选其一即可;熟练后对称情形可直接读出非零分量。

记号 ∇\nabla∇ 读作“德尔”,在直角坐标中写成 ∇=x^ ∂∂x+y^ ∂∂y+z^ ∂∂z\nabla = \hat{x}\,\dfrac{\partial}{\partial x} + \hat{y}\,\dfrac{\partial}{\partial y} + \hat{z}\,\dfrac{\partial}{\partial z}∇=x^∂x∂​+y^​∂y∂​+z^∂z∂​。它本身与后面的场量结合后,才得到梯度、散度或旋度。


梯度:标量场变化最快的方向

标量场 ϕ\phiϕ 的梯度定义为 ∇ϕ\nabla\phi∇ϕ。在直角坐标中:

∇ϕ=x^ ∂ϕ∂x+y^ ∂ϕ∂y+z^ ∂ϕ∂z\nabla\phi = \hat{x}\,\frac{\partial\phi}{\partial x} + \hat{y}\,\frac{\partial\phi}{\partial y} + \hat{z}\,\frac{\partial\phi}{\partial z}∇ϕ=x^∂x∂ϕ​+y^​∂y∂ϕ​+z^∂z∂ϕ​

梯度矢量指向 ϕ\phiϕ 增长最快的方向,其大小等于该方向上的方向导数最大值。方向导数 dϕds\dfrac{d\phi}{ds}dsdϕ​ 沿单位矢量 u^\hat{u}u^ 的方向可写成 u^⋅∇ϕ\hat{u} \cdot \nabla\phiu^⋅∇ϕ,因此沿梯度方向时该点乘取最大值 ∣∇ϕ∣|\nabla\phi|∣∇ϕ∣。在势场问题里,电场与电势满足 E⃗=−∇φ\vec{E} = -\nabla\varphiE=−∇φ,即电场指向电势下降最陡的方向,负号与电荷受力方向及电势能变化配套。

等势面是 φ=常数\varphi = \text{常数}φ=常数 的曲面;电场线与等势面处处正交,因为沿等势面切向移动时 φ\varphiφ 不变,方向导数为零,即切向与 ∇φ\nabla\varphi∇φ 垂直。画图时先画等势线或等势面,再按垂直关系勾电场线,静电问题的几何图像就完整了。

例2: 二维温度场 T(x,y)=T0+α(x2+y2)T(x,y) = T_0 + \alpha(x^2 + y^2)T(x,y)=T0​+α(x2+y2),其中 T0T_0T0​、α\alphaα 为常数,α>0\alpha > 0α>0。则

∂T∂x=2αx,∂T∂y=2αy\frac{\partial T}{\partial x} = 2\alpha x,\quad \frac{\partial T}{\partial y} = 2\alpha y∂x∂T​=2αx,∂y∂T​=2αy

故 ∇T=2α(x x^+y y^)\nabla T = 2\alpha(x\,\hat{x} + y\,\hat{y})∇T=2α(xx^+yy^​)。在点 (x,y)(x,y)(x,y) 处,梯度沿径向向外,热流若近似满足傅里叶定律 q⃗∝−∇T\vec{q} \propto -\nabla Tq​∝−∇T,则热量由高温区指向低温区,与梯度反向。等温线是 x2+y2=常数x^2+y^2 = \text{常数}x2+y2=常数,即同心圆;梯度沿半径方向,与圆周垂直,这与“梯度垂直于等值线”的一般结论一致。

在原点附近梯度为零,标量 TTT 在该点取极小值(在 α>0\alpha>0α>0 时);离原点越远,∣∇T∣|\nabla T|∣∇T∣ 越大,温度沿径向上升越快。同一形式的 x2+y2x^2+y^2x2+y2 依赖也常见于轴对称问题中只随到轴距离变化的标量场,梯度仍沿径向,等值线仍为同心圆族。

梯度把标量场升级为矢量场:知道等值面的形状,就能用梯度判断场点处标量变化的主方向,不必依赖数值搜索。


散度与旋度:源与涡

散度 ∇⋅A⃗\nabla \cdot \vec{A}∇⋅A 度量矢量场在某点附近是“向外发散”还是“向内汇聚”。直角坐标下:

∇⋅A⃗=∂Ax∂x+∂Ay∂y+∂Az∂z\nabla \cdot \vec{A} = \frac{\partial A_x}{\partial x} + \frac{\partial A_y}{\partial y} + \frac{\partial A_z}{\partial z}∇⋅A=∂x∂Ax​​+∂y∂Ay​​+∂z∂Az​​

三项分别是 xxx、yyy、zzz 方向上场分量沿各自方向的变化率之和;整体为正表示该点附近净“流出”,为负表示净“流入”。静电场高斯定律的微分形式 ∇⋅E⃗=ρ/ε0\nabla \cdot \vec{E} = \rho / \varepsilon_0∇⋅E=ρ/ε0​ 表明:电荷密度 ρ\rhoρ 正是电场的散度源,无电荷处散度为零,电场线在此既不产生也不消失。

旋度 ∇×A⃗\nabla \times \vec{A}∇×A 描述场的旋转性。直角坐标下可写成行列式便于记忆:

∇×A⃗=∣x^y^z^∂∂x∂∂y∂∂zAxAyAz∣\nabla \times \vec{A} = \begin{vmatrix} \hat{x} & \hat{y} & \hat{z} \\ \dfrac{\partial}{\partial x} & \dfrac{\partial}{\partial y} & \dfrac{\partial}{\partial z} \\ A_x & A_y & A_z \end{vmatrix}∇×A=​x^∂x∂​Ax​​y^​∂y∂​Ay​​z^∂z∂​Az​​​

按第一行展开,得到三个分量,每个分量都是两个偏导数之差。稳恒磁场的 ∇⋅B⃗=0\nabla \cdot \vec{B} = 0∇⋅B=0 说明磁场线无起点无终点;安培环路定理与法拉第定律则把 ∇×B⃗\nabla \times \vec{B}∇×B、∇×E⃗\nabla \times \vec{E}∇×E 与电流、变化磁场联系起来。

例3: A⃗=(x,y,0)\vec{A} = (x, y, 0)A=(x,y,0),三个分量为 Ax=xA_x = xAx​=x,Ay=yA_y = yAy​=y,Az=0A_z = 0Az​=0。于是 ∂Ax/∂x=1\partial A_x/\partial x = 1∂Ax​/∂x=1,∂Ay/∂y=1\partial A_y/\partial y = 1∂Ay​/∂y=1,∂Az/∂z=0\partial A_z/\partial z = 0∂Az​/∂z=0,故 ∇⋅A⃗=1+1+0=2\nabla \cdot \vec{A} = 1 + 1 + 0 = 2∇⋅A=1+1+0=2,处处为常正散度,场在 xyxyxy 平面内呈均匀扩张状。

例4: v⃗=(−ωy,ωx,0)\vec{v} = (-\omega y, \omega x, 0)v=(−ωy,ωx,0) 描写绕 zzz 轴的刚体式旋转,ω\omegaω 为常量。按行列式展开:(∇×v⃗)x=∂vz∂y−∂vy∂z=0(\nabla \times \vec{v})_x = \dfrac{\partial v_z}{\partial y} - \dfrac{\partial v_y}{\partial z} = 0(∇×v)x​=∂y∂vz​​−∂z∂vy​​=0,(∇×v⃗)y=∂vx∂z−∂vz∂x=0(\nabla \times \vec{v})_y = \dfrac{\partial v_x}{\partial z} - \dfrac{\partial v_z}{\partial x} = 0(∇×v)y​=∂z∂vx​​−∂x∂vz​​=0,(∇×v⃗)z=∂vy∂x−∂vx∂y=ω−(−ω)=2ω(\nabla \times \vec{v})_z = \dfrac{\partial v_y}{\partial x} - \dfrac{\partial v_x}{\partial y} = \omega - (-\omega) = 2\omega(∇×v)z​=∂x∂vy​​−∂y∂vx​​=ω−(−ω)=2ω。故 ∇×v⃗=(0,0,2ω)\nabla \times \vec{v} = (0, 0, 2\omega)∇×v=(0,0,2ω),旋度沿轴且大小与角速度同量级,体现涡旋运动。

量主要图像典型方程中出现的情形
散度源或汇高斯定律、不可压缩流体 ∇⋅v⃗=0\nabla \cdot \vec{v} = 0∇⋅v=0
旋度小涡旋、环流安培定律、法拉第定律、涡旋速度场

散度与旋度都是局部量:在某一点取值,由该点邻域内场的变化率决定。整体性质则通过高斯定理、斯托克斯定理转化为大范围的通量与环量。先算微分形式,再积分,是理论力学与电动力学里反复使用的路线。

不可压缩流体的速度场满足 ∇⋅v⃗=0\nabla \cdot \vec{v} = 0∇⋅v=0,表示任意微小体积内净流量为零,流体既不堆积也不抽空,只从一侧流入、从另一侧流出。有源或有汇的流动则对应散度非零的区域,与静电场中有电荷处 ∇⋅E⃗≠0\nabla \cdot \vec{E} \neq 0∇⋅E=0 的图像类似。


直角、柱坐标与球坐标

柱坐标 (ρ,φ,z)(\rho, \varphi, z)(ρ,φ,z) 与直角坐标的关系为 x=ρcos⁡φx = \rho\cos\varphix=ρcosφ,y=ρsin⁡φy = \rho\sin\varphiy=ρsinφ,z=zz = zz=z,其中 ρ≥0\rho \ge 0ρ≥0,0≤φ<2π0 \le \varphi < 2\pi0≤φ<2π。由前两式可得 ρ2=x2+y2\rho^2 = x^2 + y^2ρ2=x2+y2,tan⁡φ=y/x\tan\varphi = y/xtanφ=y/x,但由 tan⁡φ\tan\varphitanφ 反求 φ\varphiφ 时必须看 (x,y)(x,y)(x,y) 所在象限,单独用反正切函数容易丢掉 π\piπ 的修正。

球坐标 (r,θ,φ)(r, \theta, \varphi)(r,θ,φ) 常取 x=rsin⁡θcos⁡φx = r\sin\theta\cos\varphix=rsinθcosφ,y=rsin⁡θsin⁡φy = r\sin\theta\sin\varphiy=rsinθsinφ,z=rcos⁡θz = r\cos\thetaz=rcosθ,其中 r≥0r \ge 0r≥0,θ\thetaθ 为与 zzz 轴夹角,φ\varphiφ 为方位角。径向距离 r=x2+y2+z2r = \sqrt{x^2+y^2+z^2}r=x2+y2+z2​ 与柱坐标下 ρ\rhoρ 的关系为 r2=ρ2+z2r^2 = \rho^2 + z^2r2=ρ2+z2。

柱坐标体积元中的因子 ρ\rhoρ 来自 ρ dφ\rho\,d\varphiρdφ 与 dρd\rhodρ 所张成的面积元;球坐标中的 r2sin⁡θr^2\sin\thetar2sinθ 把微小盒子的三条边长乘在一起,在积分算体积或质量时直接代入即可。

从直角到柱、到球的变换,本质是同一空间点的三种参数化;画图时把 ρ\rhoρ 看成 xyxyxy 平面内的极径,把 rrr 看成到原点的直线距离,把 θ\thetaθ 看成与 zzz 轴的夹角,就不易混淆。线元长度在柱坐标中为 ds2=dρ2+ρ2dφ2+dz2ds^2 = d\rho^2 + \rho^2 d\varphi^2 + dz^2ds2=dρ2+ρ2dφ2+dz2,在球坐标中为 ds2=dr2+r2dθ2+r2sin⁡2θ dφ2ds^2 = dr^2 + r^2 d\theta^2 + r^2\sin^2\theta\, d\varphi^2ds2=dr2+r2dθ2+r2sin2θdφ2,需要弧长或动能表达式时可直接引用。

例5: 点 PPP 的直角坐标为 (1,1,0)(1, 1, 0)(1,1,0)。则 ρ=x2+y2=2\rho = \sqrt{x^2+y^2} = \sqrt{2}ρ=x2+y2​=2​,φ=arctan⁡(y/x)=π/4\varphi = \arctan(y/x) = \pi/4φ=arctan(y/x)=π/4,z=0z = 0z=0,柱坐标为 (2,π/4,0)(\sqrt{2}, \pi/4, 0)(2​,π/4,0)。又 r=2r = \sqrt{2}r=2​,θ=π/2\theta = \pi/2θ=π/2(在 xyxyxy 平面内),球坐标为 (2,π/2,π/4)(\sqrt{2}, \pi/2, \pi/4)(2​,π/2,π/4)。轴对称问题常用柱坐标,点电荷、球壳问题常用球坐标,可少算一个或两个分量。

代回验证:x=rsin⁡θcos⁡φ=2⋅1⋅22=1x = r\sin\theta\cos\varphi = \sqrt{2} \cdot 1 \cdot \dfrac{\sqrt{2}}{2} = 1x=rsinθcosφ=2​⋅1⋅22​​=1,y=2⋅1⋅22=1y = \sqrt{2} \cdot 1 \cdot \dfrac{\sqrt{2}}{2} = 1y=2​⋅1⋅22​​=1,z=2⋅0=0z = \sqrt{2} \cdot 0 = 0z=2​⋅0=0,与直角坐标一致。θ=π/2\theta = \pi/2θ=π/2 时 cos⁡θ=0\cos\theta = 0cosθ=0,z=rcos⁡θz = r\cos\thetaz=rcosθ 自然为零,说明点落在赤道平面上。

同一符号 φ\varphiφ 在文献中有时表示方位角,有时表示电势,阅读时需结合上下文区分,避免把角度代进势函数。


高斯定理与斯托克斯定理(矢量形式)

高斯定理(散度定理)把闭曲面积分与体积分联系起来:

∯SA⃗⋅dS⃗=∭V(∇⋅A⃗) dV\oiint_S \vec{A} \cdot d\vec{S} = \iiint_V (\nabla \cdot \vec{A})\,dV∬​S​A⋅dS=∭V​(∇⋅A)dV

左侧是矢量场穿过闭曲面 SSS 的通量,面积元 dS⃗d\vec{S}dS 的方向取外法向;右侧是散度在 SSS 所围体积 VVV 内的积累。静电学中用它由对称电荷分布求 E⃗\vec{E}E 十分有效:先由对称性判断 E⃗\vec{E}E 的方向与大小在面上的取值,再把 E⃗⋅dS⃗\vec{E} \cdot d\vec{S}E⋅dS 对整个闭面积分,与体内总电荷联系起来。

welearn-32091812.webp

斯托克斯定理把曲面积分与边界曲线积分联系起来:

∮CA⃗⋅dl⃗=∬S(∇×A⃗)⋅dS⃗\oint_C \vec{A} \cdot d\vec{l} = \iint_S (\nabla \times \vec{A}) \cdot d\vec{S}∮C​A⋅dl=∬S​(∇×A)⋅dS

左侧是 A⃗\vec{A}A 沿闭曲线 CCC 的环量,线元 dl⃗d\vec{l}dl 的方向与 CCC 的走向一致;右侧是旋度穿过以 CCC 为边界的曲面 SSS 的通量。曲面 SSS 的法向与 CCC 的走向仍由右手定则配合:四指沿 CCC 方向弯曲时,大拇指指向 dS⃗d\vec{S}dS 的正方向。同一边界 CCC 可以张成不同曲面,定理保证只要 SSS 以 CCC 为边界,积分值相同。

例6: 半径为 RRR 的球面,球心有点电荷 qqq,球面上各点电场大小均为 E=14πε0qR2E = \dfrac{1}{4\pi\varepsilon_0}\dfrac{q}{R^2}E=4πε0​1​R2q​,方向沿径向与外法向一致。于是 E⃗⋅dS⃗=E dS\vec{E} \cdot d\vec{S} = E\,dSE⋅dS=EdS,通量 ΦE=∯SE dS=E⋅4πR2\Phi_E = \displaystyle\oiint_S E\,dS = E \cdot 4\pi R^2ΦE​=∬​S​EdS=E⋅4πR2。代入 EEE 得 ΦE=q/ε0\Phi_E = q/\varepsilon_0ΦE​=q/ε0​,与高斯定理给出的体积分结果一致;体积分中散度仅在点电荷位置奇异,整体仍给出同一通量。

该结果与球半径 RRR 无关:只要闭曲面把点电荷包在内部,电通量恒为 q/ε0q/\varepsilon_0q/ε0​。电荷在曲面外时,穿入与穿出的通量抵消,总通量为零。斯托克斯定理一侧的典型用法是:已知 ∇×A⃗\nabla \times \vec{A}∇×A 在曲面上的分布,通过选简单曲面(如平面圆盘)把环量化为面积分;另一侧则是沿边界逐段计算线积分作检验。

定理几何对象联系的两类积分
高斯定理闭曲面 SSS 与所围体积 VVV面积分 ⇔\Leftrightarrow⇔ 体积分
斯托克斯定理曲面 SSS 与其边界 CCC线积分 ⇔\Leftrightarrow⇔ 面积分

两式都把“局部微分信息”与“整体积分信息”连成一体:高斯定理看闭曲面内外通量与体内源强;斯托克斯定理看曲面边缘环量与面上旋度分布。实际解题时常先判断对称性,再选高斯面或安培回路,最后把面积分或线积分化为代数运算。


拉普拉斯算符在不同坐标中的形式(查阅用)

拉普拉斯算符 ∇2\nabla^2∇2 对标量场定义为散度的梯度:∇2ϕ=∇⋅(∇ϕ)\nabla^2\phi = \nabla \cdot (\nabla\phi)∇2ϕ=∇⋅(∇ϕ)。直角坐标:

∇2ϕ=∂2ϕ∂x2+∂2ϕ∂y2+∂2ϕ∂z2\nabla^2\phi = \frac{\partial^2\phi}{\partial x^2} + \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2}∇2ϕ=∂x2∂2ϕ​+∂y2∂2ϕ​+∂z2∂2ϕ​

把 ∇ϕ\nabla\phi∇ϕ 的三个分量依次代入散度公式,二阶混合导数在 ϕ\phiϕ 足够光滑时会相互抵消或合并,最终只剩下三个纯二阶偏导之和。光滑性指 ϕ\phiϕ 在所讨论区域内连续且所需阶偏导存在并连续;物理上常见的势函数在电荷、质量分布之外一般都满足这一条件。柱坐标与球坐标下表达式较长,但在解拉普拉斯方程 ∇2ϕ=0\nabla^2\phi = 0∇2ϕ=0 或泊松方程时必不可少。下面给出标准形式供直接代入计算。

坐标系∇2ϕ\nabla^2\phi∇2ϕ 的展开(标准写法)
柱 (ρ,φ,z)(\rho,\varphi,z)(ρ,φ,z)1ρ∂∂ρ ⁣(ρ∂ϕ∂ρ)+1ρ2∂2ϕ∂φ2+∂2ϕ∂z2\dfrac{1}{\rho}\dfrac{\partial}{\partial\rho}\!\left(\rho\dfrac{\partial\phi}{\partial\rho}\right) + \dfrac{1}{\rho^2}\dfrac{\partial^2\phi}{\partial\varphi^2} + \dfrac{\partial^2\phi}{\partial z^2}ρ1​∂ρ∂​(ρ∂ρ∂ϕ​)+ρ21​∂φ2∂2ϕ​+∂z2∂2ϕ​
球 (r,θ,φ)(r,\theta,\varphi)(r,θ,φ)1r2∂∂r ⁣(r2∂ϕ∂r)+1r2sin⁡θ∂∂θ ⁣(sin⁡θ∂ϕ∂θ)+1r2sin⁡2θ∂2ϕ∂φ2\dfrac{1}{r^2}\dfrac{\partial}{\partial r}\!\left(r^2\dfrac{\partial\phi}{\partial r}\right) + \dfrac{1}{r^2\sin\theta}\dfrac{\partial}{\partial\theta}\!\left(\sin\theta\dfrac{\partial\phi}{\partial\theta}\right) + \dfrac{1}{r^2\sin^2\theta}\dfrac{\partial^2\phi}{\partial\varphi^2}r21​∂r∂​(r2∂r∂ϕ​)+r2sinθ1​∂θ∂​(sinθ∂θ∂ϕ​)+r2sin2θ1​∂φ2∂2ϕ​

轴对称且与 zzz 无关的柱坐标问题中,∂/∂φ=0\partial/\partial\varphi = 0∂/∂φ=0,式子明显简化;球对称时 ϕ\phiϕ 仅依赖 rrr,则 ∇2ϕ=1r2ddr ⁣(r2dϕdr)\nabla^2\phi = \dfrac{1}{r^2}\dfrac{d}{dr}\!\left(r^2\dfrac{d\phi}{dr}\right)∇2ϕ=r21​drd​(r2drdϕ​)。展开导数还可写成 d2ϕdr2+2rdϕdr\dfrac{d^2\phi}{dr^2} + \dfrac{2}{r}\dfrac{d\phi}{dr}dr2d2ϕ​+r2​drdϕ​,在求点电荷电势或引力势时经常遇到。

泊松方程 ∇2ϕ=−4πGρm\nabla^2\phi = -4\pi G\rho_m∇2ϕ=−4πGρm​(引力)或静电中的 ∇2φ=−ρ/ε0\nabla^2\varphi = -\rho/\varepsilon_0∇2φ=−ρ/ε0​ 在球对称情形下都退化为对 rrr 的常微分方程,先积出通解,再用边界条件定常数。柱坐标下若仅与 ρ\rhoρ 有关,拉普拉斯算符化为 1ρddρ ⁣(ρdϕdρ)\dfrac{1}{\rho}\dfrac{d}{d\rho}\!\left(\rho\dfrac{d\phi}{d\rho}\right)ρ1​dρd​(ρdρdϕ​),与均匀带电长直线的电势问题对应。

用直角坐标验算:ϕ=x2+y2+z2\phi = x^2 + y^2 + z^2ϕ=x2+y2+z2 时 ∂2ϕ/∂x2=2\partial^2\phi/\partial x^2 = 2∂2ϕ/∂x2=2,三项相加得 ∇2ϕ=6\nabla^2\phi = 6∇2ϕ=6。换到球对称写法 ϕ=r2\phi = r^2ϕ=r2,代入 ∇2ϕ=1r2ddr(r2dϕdr)\nabla^2\phi = \dfrac{1}{r^2}\dfrac{d}{dr}(r^2 \dfrac{d\phi}{dr})∇2ϕ=r21​drd​(r2drdϕ​),其中 dϕ/dr=2rd\phi/dr = 2rdϕ/dr=2r,r2⋅2r=2r3r^2 \cdot 2r = 2r^3r2⋅2r=2r3,再对 rrr 求导得 6r26r^26r2,除以 r2r^2r2 仍为 666,与直角坐标结果一致,可用来检查坐标公式是否抄错。


练习题

选择题

1. 下列关于梯度 ∇ϕ\nabla\phi∇ϕ 的说法,正确的是(  )

A. 梯度方向总是标量场 ϕ\phiϕ 下降最快的方向

B. 梯度的大小等于 ϕ\phiϕ 在空间中的最大值

C. 梯度方向是 ϕ\phiϕ 增长最快的方向,其模等于该方向的方向导数

D. 梯度与等值面平行

答案:C

梯度指向 ϕ\phiϕ 增加最快的方向,且该方向上的方向导数等于 ∣∇ϕ∣|\nabla\phi|∣∇ϕ∣。静电学中 E⃗=−∇φ\vec{E} = -\nabla\varphiE=−∇φ,电场沿电势降低方向,故选项 A 把符号说反。梯度大小是局部变化率,不是场的全局最大取值,B 错误。梯度垂直于等值面,D 错误。

方向导数沿任意 u^\hat{u}u^ 为 u^⋅∇ϕ\hat{u}\cdot\nabla\phiu^⋅∇ϕ,由柯西—施瓦茨不等式 ∣u^⋅∇ϕ∣≤∣∇ϕ∣|\hat{u}\cdot\nabla\phi| \le |\nabla\phi|∣u^⋅∇ϕ∣≤∣∇ϕ∣,等号当且仅当 u^\hat{u}u^ 与 ∇ϕ\nabla\phi∇ϕ 同向或反向时成立,故增长最快方向必沿 ∇ϕ\nabla\phi∇ϕ。


2. 直角坐标下,矢量场 A⃗=(x,y,z)\vec{A} = (x, y, z)A=(x,y,z) 的散度 ∇⋅A⃗\nabla \cdot \vec{A}∇⋅A 为(  )

A. 000

B. 111

C. 222

D. 333

答案:D

∂x∂x+∂y∂y+∂z∂z=1+1+1=3\dfrac{\partial x}{\partial x} + \dfrac{\partial y}{\partial y} + \dfrac{\partial z}{\partial z} = 1 + 1 + 1 = 3∂x∂x​+∂y∂y​+∂z∂z​=1+1+1=3。

该场 A⃗\vec{A}A 沿径向从原点向外增长,xxx、yyy、zzz 三个方向分量各自随坐标线性增大,三个方向偏导各贡献 111,合起来散度处处为常数 333,与位置无关。


3. 对稳恒磁场 B⃗\vec{B}B,下列恒成立的是(  )

A. ∇×B⃗=0\nabla \times \vec{B} = 0∇×B=0 处处成立

B. ∇⋅B⃗=0\nabla \cdot \vec{B} = 0∇⋅B=0 处处成立

C. ∇⋅B⃗=μ0ρm\nabla \cdot \vec{B} = \mu_0 \rho_m∇⋅B=μ0​ρm​,ρm\rho_mρm​ 为磁荷密度

D. B⃗\vec{B}B 必为保守场

答案:B

磁场无源,磁感应线闭合,故 ∇⋅B⃗=0\nabla \cdot \vec{B} = 0∇⋅B=0。有电流处 ∇×B⃗=μ0J⃗\nabla \times \vec{B} = \mu_0 \vec{J}∇×B=μ0​J,一般不为零,A 错误。不存在磁单极磁荷的标准表述,C 错误。有旋场一般不是保守场,D 错误。

稳恒情形麦克斯韦方程组中 ∇⋅B⃗=0\nabla\cdot\vec{B}=0∇⋅B=0 与 ∇×H⃗=J⃗\nabla\times\vec{H}=\vec{J}∇×H=J(真空下 B⃗=μ0H⃗\vec{B}=\mu_0\vec{H}B=μ0​H)配套;散度恒为零是对磁场源结构的总体约束,不因是否有电流而改变。


4. 高斯定理 ∯SA⃗⋅dS⃗=∭V(∇⋅A⃗) dV\displaystyle\oiint_S \vec{A} \cdot d\vec{S} = \iiint_V (\nabla \cdot \vec{A})\,dV∬​S​A⋅dS=∭V​(∇⋅A)dV 中,曲面 SSS 与体积 VVV 应满足的关系是(  )

A. SSS 为任意开曲面,VVV 为 SSS 一侧的任意区域

B. SSS 为闭曲面,VVV 为 SSS 所包围的体积,且法向通常取外法向

C. SSS 与 VVV 无必然关系

D. SSS 必须是平面

答案:B

散度定理要求 SSS 为闭曲面,VVV 为其内部区域,通量积分用外法向与体积分中的散度配套。开曲面适用于斯托克斯定理与曲面积分,不直接套此闭曲面积分形式。

若误把开曲面当作闭曲面使用,就没有唯一确定的“内部体积”,通量与体积分之间的恒等关系也失去意义。外法向约定保证“流出为正”,与体积分中散度的符号约定一致。


计算题

计算题侧重分步求偏导、代入数值与单位自检;选择题侧重概念辨析与定理条件。完成后可把梯度、散度、旋度的结果代回原场,检查量纲与对称性是否与题设一致。

书写答案时保持矢量分量顺序与题中坐标轴约定一致,避免符号差错。

5. 标量场 ϕ(x,y,z)=x2y+z\phi(x,y,z) = x^2 y + zϕ(x,y,z)=x2y+z,求 ∇ϕ\nabla\phi∇ϕ,并求点 (1,2,0)(1, 2, 0)(1,2,0) 处梯度的模。

解:

对 xxx 求偏导时把 yyy、zzz 看作常数,得 ∂ϕ∂x=2xy\dfrac{\partial\phi}{\partial x} = 2xy∂x∂ϕ​=2xy。对 yyy 求偏导时把 xxx、zzz 看作常数,得 ∂ϕ∂y=x2\dfrac{\partial\phi}{\partial y} = x^2∂y∂ϕ​=x2。对 zzz 求偏导得 ∂ϕ∂z=1\dfrac{\partial\phi}{\partial z} = 1∂z∂ϕ​=1。故

∇ϕ=2xy x^+x2 y^+z^\nabla\phi = 2xy\,\hat{x} + x^2\,\hat{y} + \hat{z}∇ϕ=2xyx^+x2y^​+z^

在 (1,2,0)(1, 2, 0)(1,2,0) 处代入 x=1x=1x=1、y=2y=2y=2、z=0z=0z=0:∇ϕ=4x^+y^+z^\nabla\phi = 4\hat{x} + \hat{y} + \hat{z}∇ϕ=4x^+y^​+z^。模长用勾股定理:

∣∇ϕ∣=42+12+12=18=32|\nabla\phi| = \sqrt{4^2 + 1^2 + 1^2} = \sqrt{18} = 3\sqrt{2}∣∇ϕ∣=42+12+12​=18​=32​

量纲上若 x,y,zx,y,zx,y,z 以 m\text{m}m 计,ϕ\phiϕ 以 V\text{V}V 计,则梯度单位为 V/m\text{V/m}V/m,与电场强度单位一致。


6. 用柱坐标表示,点 PPP 的直角坐标为 (0,3,4)(0, 3, 4)(0,3,4)。求 ρ\rhoρ、φ\varphiφ(取 0≤φ<2π0 \le \varphi < 2\pi0≤φ<2π)及 zzz;并求该点到原点距离 rrr(球坐标中的径向距离)。

解:

柱坐标径向距离在 xyxyxy 平面内量度:ρ=x2+y2=0+9=3 m \rho = \sqrt{x^2 + y^2} = \sqrt{0 + 9} = 3\ \text{m}ρ=x2+y2​=0+9​=3 m(长度单位仅示意)。

方位角 φ\varphiφ:x=0x=0x=0 且 y=3>0y=3>0y=3>0,点落在 yyy 轴正半轴,故 φ=π/2\varphi = \pi/2φ=π/2,不能用 arctan⁡(y/x)\arctan(y/x)arctan(y/x) 的单一主值公式机械代入,需按象限判定。

zzz 与直角坐标相同:z=4 m z = 4\ \text{m}z=4 m。

球坐标径向距离 r=ρ2+z2=9+16=5 m r = \sqrt{\rho^2 + z^2} = \sqrt{9 + 16} = 5\ \text{m}r=ρ2+z2​=9+16​=5 m,与 r=x2+y2+z2r = \sqrt{x^2+y^2+z^2}r=x2+y2+z2​ 一致,说明两套公式配套。

下一章质点系与刚体力学