--- title: '正交曲线坐标系与矢量算符' description: '从坐标映射与拉梅系数出发，推导线元、面元、体积元、梯度、散度与旋度，并以球坐标系为例。' pubDate: '2026-10-20' credit: 'Zewei Deng' tags: ['数学物理方法', '矢量分析', '正交曲线坐标系'] download: '/downloads/orthogonal-curvilinear-coordinates.md' --- 原文题名：**正交曲线坐标系基础** 原文日期：2026 年 3 月 21 日 本文由《正交曲线坐标系与矢量算符.pdf》转换为 Markdown，保留原文的章节层级与公式编号（1）至（22）。 ## 目录 - [1 正交曲线坐标系基础](#basics) - 1.1 正交曲线坐标系与拉梅系数 - 1.1.1 坐标映射与推移矢量 - 1.1.2 拉梅系数的定义与物理意义 - 1.1.3 局部正交基矢的构建 - 1.2 微元的几何由来 - 1.2.1 单位线元：全微分的结果 - 1.2.2 面元：矢量叉乘的几何意义 - 1.2.3 体积元：标量三重积 - 1.3 算符的导出：通量与环流密度定义 - 1.3.1 梯度（Gradient） - [2 矢量算符的几何导出](#operators) - 2.1 散度的推导 - 2.1.1 计算第一个坐标方向的净通量 - 2.1.2 汇总与最终公式 - 2.2 旋度的推导 - 2.2.1 封闭路径积分 - 2.2.2 环流密度计算 - 2.2.3 算符通用行列式 - [3 应用实例：球坐标系](#spherical) - 3.1 基矢与拉梅系数计算 - 3.2 微元表示 - 3.3 矢量算符公式 ## 1 正交曲线坐标系基础 在电动力学中，为了利用物理系统的对称性（如球对称或轴对称），我们通常需要在非直角坐标系下处理矢量场。假设存在一组正交曲线坐标 $(q_1,q_2,q_3)$。 ### 1.1 正交曲线坐标系与拉梅系数 在直角坐标系中，坐标的变化量（如 $\mathrm{d}x$）直接等同于空间中的位移长度。然而，在处理具有球对称或轴对称的物理系统（如点电荷或长直导线）时，我们引入了广义的正交曲线坐标 $(q_1,q_2,q_3)$。 在这些坐标系中，坐标分量本身的量纲可能不再是长度（例如球坐标中的角度 $\theta$ 和 $\phi$）。为了在物理定律中保持欧几里得空间的度量属性，我们需要一套系数来将“坐标的变化”缩放为“真实的物理距离”，这便是拉梅系数的几何本质。 #### 1.1.1 坐标映射与推移矢量 假设空间中任一点的位置矢量为 $\mathbf{r}$，它是广义坐标的函数： $$ \begin{aligned} \mathbf{r} &=\mathbf{r}(q_1,q_2,q_3)\\ &=x(q_1,q_2,q_3)\hat{\mathbf{i}} +y(q_1,q_2,q_3)\hat{\mathbf{j}} +z(q_1,q_2,q_3)\hat{\mathbf{k}}. \end{aligned} \tag{1} $$ 当我们让其中一个坐标 $q_i$ 发生微小变动 $\mathrm{d}q_i$，而保持其他坐标不变时，位置矢量的改变方向即为该点的局部坐标轴方向。我们定义**推移矢量**（Coordinate Tangent Vector）为： $$ \mathbf{v}_i=\frac{\partial\mathbf{r}}{\partial q_i}. \tag{2} $$ 这个矢量指向 $q_i$ 增加的方向，其模长代表了该坐标方向上的“拉伸程度”。 #### 1.1.2 拉梅系数的定义与物理意义 **定义 1.1（拉梅系数）**：推移矢量 $\mathbf{v}_i$ 的模长称为该坐标方向的**拉梅系数**（Lame Coefficients），通常记作 $h_i$： $$ h_i =\left|\frac{\partial\mathbf{r}}{\partial q_i}\right| =\sqrt{ \left(\frac{\partial x}{\partial q_i}\right)^2 +\left(\frac{\partial y}{\partial q_i}\right)^2 +\left(\frac{\partial z}{\partial q_i}\right)^2 }. \tag{3} $$ **物理理解**：拉梅系数本质上是坐标系的“比例尺”。在物理位移 $\mathrm{d}l_i=h_i\,\mathrm{d}q_i$ 中： - 若 $q_i$ 本身具有长度量纲（如 $r$），则 $h_i$ 为无量纲的纯数（通常为 $1$）。 - 若 $q_i$ 为无量纲量（如角度 $\theta$），则 $h_i$ 必然具有长度量纲，以保证 $\mathrm{d}l_i$ 的单位为米。 #### 1.1.3 局部正交基矢的构建 在正交坐标系中，任何一点的三个推移矢量 $\mathbf{v}_1,\mathbf{v}_2,\mathbf{v}_3$ 始终两两垂直。为了方便矢量运算，我们将这些矢量归一化，得到该点处的**局部单位正交基**： $$ \hat{\mathbf{e}}_i =\frac{\mathbf{v}_i}{|\mathbf{v}_i|} =\frac{1}{h_i}\frac{\partial\mathbf{r}}{\partial q_i}. \tag{4} $$ **注意**：与直角坐标系的固定基矢 $\{\hat{\mathbf{i}},\hat{\mathbf{j}},\hat{\mathbf{k}}\}$ 不同，曲线坐标系的基矢 $\{\hat{\mathbf{e}}_1,\hat{\mathbf{e}}_2,\hat{\mathbf{e}}_3\}$ 的方向通常会随着空间位置的改变而改变。在后续推导导数算符时，必须考虑到基矢本身的导数项。 ### 1.2 微元的几何由来 #### 1.2.1 单位线元：全微分的结果 根据多元函数全微分的定义，空间位移矢量 $\mathrm{d}\mathbf{r}$ 可以表示为： $$ \begin{aligned} \mathrm{d}\mathbf{l} =\mathrm{d}\mathbf{r} &=\frac{\partial\mathbf{r}}{\partial q_1}\mathrm{d}q_1 +\frac{\partial\mathbf{r}}{\partial q_2}\mathrm{d}q_2 +\frac{\partial\mathbf{r}}{\partial q_3}\mathrm{d}q_3\\ &=h_1\mathrm{d}q_1\hat{\mathbf{e}}_1 +h_2\mathrm{d}q_2\hat{\mathbf{e}}_2 +h_3\mathrm{d}q_3\hat{\mathbf{e}}_3. \end{aligned} \tag{5} $$ 这说明 $h_i$ 的作用是将“坐标增量”转换为“物理长度增量”。 #### 1.2.2 面元：矢量叉乘的几何意义 在正交坐标系中，由 $\mathrm{d}q_i$ 和 $\mathrm{d}q_j$ 围成的微小面元是一个矩形。利用矢量叉乘的几何意义（两个矢量围成的平行四边形面积及其法向），面元矢量 $\mathrm{d}\mathbf{a}$ 表示为： $$ \begin{aligned} \mathrm{d}\mathbf{a}_{ij} &=(h_i\mathrm{d}q_i\hat{\mathbf{e}}_i) \times(h_j\mathrm{d}q_j\hat{\mathbf{e}}_j)\\ &=(h_i h_j\,\mathrm{d}q_i\,\mathrm{d}q_j)\hat{\mathbf{e}}_k. \end{aligned} \tag{6} $$ 这里 $(i,j,k)$ 按右手正交基的循环顺序取值。例如，常数 $q_3$ 表面的面元为 $\mathrm{d}\mathbf{a}_3=h_1h_2\,\mathrm{d}q_1\,\mathrm{d}q_2\,\hat{\mathbf{e}}_3$。 #### 1.2.3 体积元：标量三重积 由三个坐标微元围成的微小体元的体积，等于三个边矢量的标量三重积： $$ \begin{aligned} \mathrm{d}\tau &=\left| (h_1\mathrm{d}q_1\hat{\mathbf{e}}_1) \cdot\left( h_2\mathrm{d}q_2\hat{\mathbf{e}}_2 \times h_3\mathrm{d}q_3\hat{\mathbf{e}}_3 \right) \right|\\ &=h_1h_2h_3\,\mathrm{d}q_1\,\mathrm{d}q_2\,\mathrm{d}q_3. \end{aligned} \tag{7} $$ ### 1.3 算符的导出：通量与环流密度定义 #### 1.3.1 梯度（Gradient） 根据方向导数的定义，标量场 $f$ 的全微分 $\mathrm{d}f=\nabla f\cdot\mathrm{d}\mathbf{l}$。对比全微分公式 $\mathrm{d}f=\sum_i\frac{\partial f}{\partial q_i}\mathrm{d}q_i$，可得： $$ \nabla f =\frac{1}{h_1}\frac{\partial f}{\partial q_1}\hat{\mathbf{e}}_1 +\frac{1}{h_2}\frac{\partial f}{\partial q_2}\hat{\mathbf{e}}_2 +\frac{1}{h_3}\frac{\partial f}{\partial q_3}\hat{\mathbf{e}}_3. \tag{8} $$ ## 2 矢量算符的几何导出 在正交曲线坐标系 $(q_1,q_2,q_3)$ 中，关键在于所有涉及“长度”的物理量都必须乘上对应的拉梅系数 $h_i$。我们考虑一个体积元 $\mathrm{d}\tau=(h_1\mathrm{d}q_1)(h_2\mathrm{d}q_2)(h_3\mathrm{d}q_3)$。 ### 2.1 散度（$\nabla\cdot\mathbf{A}$）的推导 散度定义为单位体积内的净通量。我们计算穿过微元六面体表面的总通量。 #### 2.1.1 计算 $q_1$ 方向的净通量 考虑一对垂直于 $q_1$ 方向的表面（面积为 $\mathrm{d}S_1=h_2h_3\,\mathrm{d}q_2\,\mathrm{d}q_3$）： **后表面（位于 $q_1$ 处）**：流入通量为该处的场强分量乘以面积。 $$ \begin{aligned} \Phi_{\mathrm{back}} ={}&-A_1(q_1,q_2,q_3)\\ &{}\cdot h_2(q_1,q_2,q_3)h_3(q_1,q_2,q_3) \,\mathrm{d}q_2\,\mathrm{d}q_3. \end{aligned} \tag{9} $$ 负号表示流入。 **前表面（位于 $q_1+\mathrm{d}q_1$ 处）**：利用一阶泰勒展开，该处的通量函数 $(A_1h_2h_3)$ 变为： $$ \Phi_{\mathrm{front}} =\left[ A_1h_2h_3 +\frac{\partial}{\partial q_1}(A_1h_2h_3)\,\mathrm{d}q_1 \right]\mathrm{d}q_2\,\mathrm{d}q_3. \tag{10} $$ **$q_1$ 方向的净通量**： $$ \begin{aligned} \mathrm{d}\Phi_1 &=\Phi_{\mathrm{front}}+\Phi_{\mathrm{back}}\\ &=\frac{\partial}{\partial q_1}(A_1h_2h_3) \,\mathrm{d}q_1\,\mathrm{d}q_2\,\mathrm{d}q_3. \end{aligned} \tag{11} $$ #### 2.1.2 汇总与最终公式 同理可得 $q_2,q_3$ 方向的贡献。总净通量 $\mathrm{d}\Phi=\mathrm{d}\Phi_1+\mathrm{d}\Phi_2+\mathrm{d}\Phi_3$。根据散度定义 $\nabla\cdot\mathbf{A}=\frac{\mathrm{d}\Phi}{\mathrm{d}\tau}$，其中 $\mathrm{d}\tau=h_1h_2h_3\,\mathrm{d}q_1\,\mathrm{d}q_2\,\mathrm{d}q_3$： $$ \begin{aligned} \nabla\cdot\mathbf{A} =\frac{1}{h_1h_2h_3}\Bigg[ &\frac{\partial}{\partial q_1}(h_2h_3A_1)\\ +&\frac{\partial}{\partial q_2}(h_1h_3A_2)\\ +&\frac{\partial}{\partial q_3}(h_1h_2A_3) \Bigg]. \end{aligned} \tag{12} $$ ### 2.2 旋度（$\nabla\times\mathbf{A}$）的推导 旋度的分量定义为单位面积的环流。我们推导其第 $3$ 分量 $(\nabla\times\mathbf{A})_3$，即绕 $q_1-q_2$ 面回路的环流。 #### 2.2.1 封闭路径积分 考虑由 $$ \begin{aligned} (q_1,q_2) &\longrightarrow(q_1+\mathrm{d}q_1,q_2)\\ &\longrightarrow(q_1+\mathrm{d}q_1,q_2+\mathrm{d}q_2)\\ &\longrightarrow(q_1,q_2+\mathrm{d}q_2) \end{aligned} $$ 构成的矩形路径，最后回到起点。 **路径 1（沿 $q_1$）**： $$ \mathrm{d}\Gamma_1=[A_1h_1]_{q_2}\,\mathrm{d}q_1. $$ **路径 2（沿 $q_2$）**： $$ \begin{aligned} \mathrm{d}\Gamma_2 &=[A_2h_2]_{q_1+\mathrm{d}q_1}\,\mathrm{d}q_2\\ &\approx\left[ A_2h_2+\frac{\partial}{\partial q_1}(A_2h_2)\,\mathrm{d}q_1 \right]\mathrm{d}q_2. \end{aligned} $$ **路径 3（沿 $-q_1$）**： $$ \begin{aligned} \mathrm{d}\Gamma_3 &=-[A_1h_1]_{q_2+\mathrm{d}q_2}\,\mathrm{d}q_1\\ &\approx-\left[ A_1h_1+\frac{\partial}{\partial q_2}(A_1h_1)\,\mathrm{d}q_2 \right]\mathrm{d}q_1. \end{aligned} $$ **路径 4（沿 $-q_2$）**： $$ \mathrm{d}\Gamma_4=-[A_2h_2]_{q_1}\,\mathrm{d}q_2. $$ #### 2.2.2 环流密度计算 叠加四段路径环流，抵消零阶项后得： $$ \mathrm{d}\Gamma_{12} =\left[ \frac{\partial}{\partial q_1}(A_2h_2) -\frac{\partial}{\partial q_2}(A_1h_1) \right]\mathrm{d}q_1\,\mathrm{d}q_2. \tag{13} $$ 该回路包围的面积为 $\mathrm{d}a_3=h_1h_2\,\mathrm{d}q_1\,\mathrm{d}q_2$。由定义 $(\nabla\times\mathbf{A})_3=\frac{\mathrm{d}\Gamma_{12}}{\mathrm{d}a_3}$ 得： $$ (\nabla\times\mathbf{A})_3 =\frac{1}{h_1h_2} \left[ \frac{\partial}{\partial q_1}(h_2A_2) -\frac{\partial}{\partial q_2}(h_1A_1) \right]. \tag{14} $$ #### 2.2.3 算符通用行列式 通过轮换下标，旋度算符可统一表示为： $$ \nabla\times\mathbf{A} =\frac{1}{h_1h_2h_3} \begin{vmatrix} h_1\hat{\mathbf{e}}_1 & h_2\hat{\mathbf{e}}_2 & h_3\hat{\mathbf{e}}_3\\ \dfrac{\partial}{\partial q_1} & \dfrac{\partial}{\partial q_2} & \dfrac{\partial}{\partial q_3}\\ h_1A_1 & h_2A_2 & h_3A_3 \end{vmatrix}. \tag{15} $$ ## 3 应用实例：球坐标系 $(r,\theta,\phi)$ 球坐标与直角坐标的关系为： $$ x=r\sin\theta\cos\phi,\qquad y=r\sin\theta\sin\phi,\qquad z=r\cos\theta. $$ ### 3.1 基矢与拉梅系数计算 根据定义计算拉梅系数： - $h_r=\left|\frac{\partial\mathbf{r}}{\partial r}\right|=1$； - $h_\theta=\left|\frac{\partial\mathbf{r}}{\partial\theta}\right|=r$； - $h_\phi=\left|\frac{\partial\mathbf{r}}{\partial\phi}\right|=r\sin\theta$。 直角坐标表示下的基矢： $$ \hat{\mathbf{e}}_r =\sin\theta\cos\phi\,\hat{\mathbf{i}} +\sin\theta\sin\phi\,\hat{\mathbf{j}} +\cos\theta\,\hat{\mathbf{k}}. \tag{16} $$ $$ \hat{\mathbf{e}}_\theta =\cos\theta\cos\phi\,\hat{\mathbf{i}} +\cos\theta\sin\phi\,\hat{\mathbf{j}} -\sin\theta\,\hat{\mathbf{k}}. \tag{17} $$ $$ \hat{\mathbf{e}}_\phi =-\sin\phi\,\hat{\mathbf{i}} +\cos\phi\,\hat{\mathbf{j}}. \tag{18} $$ ### 3.2 微元表示 **线元**： $$ \mathrm{d}\mathbf{l} =\mathrm{d}r\,\hat{\mathbf{e}}_r +r\,\mathrm{d}\theta\,\hat{\mathbf{e}}_\theta +r\sin\theta\,\mathrm{d}\phi\,\hat{\mathbf{e}}_\phi. $$ **面元（常数 $r$ 表面）**： $$ \mathrm{d}\mathbf{a}_r =r^2\sin\theta\,\mathrm{d}\theta\,\mathrm{d}\phi\,\hat{\mathbf{e}}_r. $$ **体积元**： $$ \mathrm{d}\tau=r^2\sin\theta\,\mathrm{d}r\,\mathrm{d}\theta\,\mathrm{d}\phi. $$ ### 3.3 矢量算符公式 将 $h_i$ 代入通用算符公式。 **梯度**： $$ \nabla f =\frac{\partial f}{\partial r}\hat{\mathbf{e}}_r +\frac{1}{r}\frac{\partial f}{\partial\theta}\hat{\mathbf{e}}_\theta +\frac{1}{r\sin\theta}\frac{\partial f}{\partial\phi}\hat{\mathbf{e}}_\phi. \tag{19} $$ **散度**： $$ \begin{aligned} \nabla\cdot\mathbf{A} &=\frac{1}{(1)(r)(r\sin\theta)} \Bigg[ \frac{\partial}{\partial r}(r^2\sin\theta\,A_r)\\ &\qquad\qquad +\frac{\partial}{\partial\theta}(r\sin\theta\,A_\theta) +\frac{\partial}{\partial\phi}(rA_\phi) \Bigg]\\ &=\frac{1}{r^2}\frac{\partial}{\partial r}(r^2A_r) +\frac{1}{r\sin\theta}\frac{\partial}{\partial\theta}(\sin\theta\,A_\theta) +\frac{1}{r\sin\theta}\frac{\partial A_\phi}{\partial\phi}. \end{aligned} \tag{20} $$ **旋度（以 $r$ 分量为例）**： $$ \begin{aligned} (\nabla\times\mathbf{A})_r &=\frac{1}{h_\theta h_\phi} \left[ \frac{\partial}{\partial\theta}(h_\phi A_\phi) -\frac{\partial}{\partial\phi}(h_\theta A_\theta) \right]\\ &=\frac{1}{r^2\sin\theta} \left[ \frac{\partial}{\partial\theta}(r\sin\theta\,A_\phi) -\frac{\partial}{\partial\phi}(rA_\theta) \right]. \end{aligned} \tag{21} $$ 简化常数 $r$ 后得： $$ (\nabla\times\mathbf{A})_r =\frac{1}{r\sin\theta} \left[ \frac{\partial}{\partial\theta}(\sin\theta\,A_\phi) -\frac{\partial A_\theta}{\partial\phi} \right]. \tag{22} $$