在圆柱和球坐标系中的del

下面是常用于正交曲线坐标系（英语：Curvilinear coordinates）中的一些向量微积分公式。 本文对球坐标使用标准符号ISO 80000-2，它取代了ISO 31-11，（部分其他来源可能有着颠倒θ和φ的定义）： 极角表示为θ：它是在z轴与连接原点和目标点的径向向量之间的角度。 方位角表示为φ：它是在x轴与径向向量在xy面上的投影之间的角度。 函数atan2(y, x)可以用于替代数学函数arctan(y/x)。这是由于它的定义域和像的缘故，经典arctan函数的像为(−π/2, +π/2)，而atan2定义的像为(−π, π]。 ^α 本页对极角采用 θ {\displaystyle \theta } 对方位角采用 φ {\displaystyle \varphi } ，这是在物理学中常用的符号。某些来源在这些公式中对方位角采用 θ {\displaystyle \theta } 对极角采用 φ {\displaystyle \varphi } ，这是常用数学符号，如果需要这种数学公式，可对换上表公式中的 θ {\displaystyle \theta } 和 φ {\displaystyle \varphi } 。 div grad ⁡ f ≡ ∇ ⋅ ∇ f ≡ ∇ 2 f {\displaystyle \operatorname {div} \,\operatorname {grad} f\equiv \nabla \cdot \nabla f\equiv \nabla ^{2}f} curl grad ⁡ f ≡ ∇ × ∇ f = 0 {\displaystyle \operatorname {curl} \,\operatorname {grad} f\equiv \nabla \times \nabla f=\mathbf {0} } div curl ⁡ A ≡ ∇ ⋅ ( ∇ × A ) = 0 {\displaystyle \operatorname {div} \,\operatorname {curl} \mathbf {A} \equiv \nabla \cdot (\nabla \times \mathbf {A} )=0} curl curl ⁡ A ≡ ∇ × ( ∇ × A ) = ∇ ( ∇ ⋅ A ) − ∇ 2 A {\displaystyle \operatorname {curl} \,\operatorname {curl} \mathbf {A} \equiv \nabla \times (\nabla \times \mathbf {A} )=\nabla (\nabla \cdot \mathbf {A} )-\nabla ^{2}\mathbf {A} } （del的拉格朗日公式） ∇ 2 ( f g ) = f ∇ 2 g + 2 ∇ f ⋅ ∇ g + g ∇ 2 f {\displaystyle \nabla ^{2}(fg)=f\nabla ^{2}g+2\nabla f\cdot \nabla g+g\nabla ^{2}f} div ⁡ A = lim V → 0 ∬ ∂ V A ⋅ d S ∭ V d V = A x ( x + d x ) d y d z − A x ( x ) d y d z + A y ( y + d y ) d x d z − A y ( y ) d x d z + A z ( z + d z ) d x d y − A z ( z ) d x d y d x d y d z = ∂ A x ∂ x + ∂ A y ∂ y + ∂ A z ∂ z {\displaystyle {\begin{aligned}\operatorname {div} \mathbf {A} =\lim _{V\to 0}{\frac {\iint _{\partial V}\mathbf {A} \cdot d\mathbf {S} }{\iiint _{V}dV}}&={\frac {A_{x}(x+dx)dydz-A_{x}(x)dydz+A_{y}(y+dy)dxdz-A_{y}(y)dxdz+A_{z}(z+dz)dxdy-A_{z}(z)dxdy}{dxdydz}}\\&={\frac {\partial A_{x}}{\partial x}}+{\frac {\partial A_{y}}{\partial y}}+{\frac {\partial A_{z}}{\partial z}}\end{aligned}}} ( curl ⁡ A ) x = lim S ⊥ x ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A z ( y + d y ) d z − A z ( y ) d z + A y ( z ) d y − A y ( z + d z ) d y d y d z = ∂ A z ∂ y − ∂ A y ∂ z {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{x}=\lim _{S^{\perp \mathbf {\hat {x}} }\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}&={\frac {A_{z}(y+dy)dz-A_{z}(y)dz+A_{y}(z)dy-A_{y}(z+dz)dy}{dydz}}\\&={\frac {\partial A_{z}}{\partial y}}-{\frac {\partial A_{y}}{\partial z}}\end{aligned}}} ( curl ⁡ A ) y {\displaystyle (\operatorname {curl} \mathbf {A} )_{y}} 和 ( curl ⁡ A ) z {\displaystyle (\operatorname {curl} \mathbf {A} )_{z}} 的表达式可以同理得出。 註：第一式中的 A x ( x + d x ) {\displaystyle A_{x}(x+dx)} 是 A x {\displaystyle A_{x}} 在 x + d x {\displaystyle x+dx} 時的量值，並非 A x {\displaystyle A_{x}} 值乘上 x + d x {\displaystyle x+dx} 。以下圓柱座標、球座標的推導中亦然。 div ⁡ A = lim V → 0 ∬ ∂ V A ⋅ d S ∭ V d V = A ρ ( ρ + d ρ ) ( ρ + d ρ ) d ϕ d z − A ρ ( ρ ) ρ d ϕ d z + A ϕ ( ϕ + d ϕ ) d ρ d z − A ϕ ( ϕ ) d ρ d z + A z ( z + d z ) d ρ ( ρ + d ρ / 2 ) d ϕ − A z ( z ) d ρ ( ρ + d ρ / 2 ) d ϕ ρ d ϕ d ρ d z = 1 ρ ∂ ( ρ A ρ ) ∂ ρ + 1 ρ ∂ A ϕ ∂ ϕ + ∂ A z ∂ z {\displaystyle {\begin{aligned}\operatorname {div} \mathbf {A} &=\lim _{V\to 0}{\frac {\iint _{\partial V}\mathbf {A} \cdot d\mathbf {S} }{\iiint _{V}dV}}\\&={\frac {A_{\rho }(\rho +d\rho )(\rho +d\rho )d\phi dz-A_{\rho }(\rho )\rho d\phi dz+A_{\phi }(\phi +d\phi )d\rho dz-A_{\phi }(\phi )d\rho dz+A_{z}(z+dz)d\rho (\rho +d\rho /2)d\phi -A_{z}(z)d\rho (\rho +d\rho /2)d\phi }{\rho d\phi d\rho dz}}\\&={\frac {1}{\rho }}{\frac {\partial (\rho A_{\rho })}{\partial \rho }}+{\frac {1}{\rho }}{\frac {\partial A_{\phi }}{\partial \phi }}+{\frac {\partial A_{z}}{\partial z}}\end{aligned}}} ( curl ⁡ A ) ρ = lim S ⊥ ρ ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A ϕ ( z ) ( ρ + d ρ ) d ϕ − A ϕ ( z + d z ) ( ρ + d ρ ) d ϕ + A z ( ϕ + d ϕ ) d z − A z ( ϕ ) d z ( ρ + d ρ ) d ϕ d z = − ∂ A ϕ ∂ z + 1 ρ ∂ A z ∂ ϕ {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{\rho }&=\lim _{S^{\perp {\boldsymbol {\hat {\rho }}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}\\&={\frac {A_{\phi }(z)(\rho +d\rho )d\phi -A_{\phi }(z+dz)(\rho +d\rho )d\phi +A_{z}(\phi +d\phi )dz-A_{z}(\phi )dz}{(\rho +d\rho )d\phi dz}}\\&=-{\frac {\partial A_{\phi }}{\partial z}}+{\frac {1}{\rho }}{\frac {\partial A_{z}}{\partial \phi }}\end{aligned}}} ( curl ⁡ A ) ϕ = lim S ⊥ ϕ ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A z ( ρ ) d z − A z ( ρ + d ρ ) d z + A ρ ( z + d z ) d ρ − A ρ ( z ) d ρ d ρ d z = − ∂ A z ∂ ρ + ∂ A ρ ∂ z {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{\phi }&=\lim _{S^{\perp {\boldsymbol {\hat {\phi }}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}\\&={\frac {A_{z}(\rho )dz-A_{z}(\rho +d\rho )dz+A_{\rho }(z+dz)d\rho -A_{\rho }(z)d\rho }{d\rho dz}}\\&=-{\frac {\partial A_{z}}{\partial \rho }}+{\frac {\partial A_{\rho }}{\partial z}}\end{aligned}}} ( curl ⁡ A ) z = lim S ⊥ z ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A ρ ( ϕ ) d ρ − A ρ ( ϕ + d ϕ ) d ρ + A ϕ ( ρ + d ρ ) ( ρ + d ρ ) d ϕ − A ϕ ( ρ ) ρ d ϕ ρ d ρ d ϕ = − 1 ρ ∂ A ρ ∂ ϕ + 1 ρ ∂ ( ρ A ϕ ) ∂ ρ {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{z}&=\lim _{S^{\perp {\boldsymbol {\hat {z}}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}\\&={\frac {A_{\rho }(\phi )d\rho -A_{\rho }(\phi +d\phi )d\rho +A_{\phi }(\rho +d\rho )(\rho +d\rho )d\phi -A_{\phi }(\rho )\rho d\phi }{\rho d\rho d\phi }}\\&=-{\frac {1}{\rho }}{\frac {\partial A_{\rho }}{\partial \phi }}+{\frac {1}{\rho }}{\frac {\partial (\rho A_{\phi })}{\partial \rho }}\end{aligned}}} curl ⁡ A = ( curl ⁡ A ) ρ ρ ^ + ( curl ⁡ A ) ϕ ϕ ^ + ( curl ⁡ A ) z z ^ = ( 1 ρ ∂ A z ∂ ϕ − ∂ A ϕ ∂ z ) ρ ^ + ( ∂ A ρ ∂ z − ∂ A z ∂ ρ ) ϕ ^ + 1 ρ ( ∂ ( ρ A ϕ ) ∂ ρ − ∂ A ρ ∂ ϕ ) z ^ {\displaystyle {\begin{aligned}\operatorname {curl} \mathbf {A} &=(\operatorname {curl} \mathbf {A} )_{\rho }{\hat {\boldsymbol {\rho }}}+(\operatorname {curl} \mathbf {A} )_{\phi }{\hat {\boldsymbol {\phi }}}+(\operatorname {curl} \mathbf {A} )_{z}{\hat {\boldsymbol {z}}}\\&=\left({\frac {1}{\rho }}{\frac {\partial A_{z}}{\partial \phi }}-{\frac {\partial A_{\phi }}{\partial z}}\right){\hat {\boldsymbol {\rho }}}+\left({\frac {\partial A_{\rho }}{\partial z}}-{\frac {\partial A_{z}}{\partial \rho }}\right){\hat {\boldsymbol {\phi }}}+{\frac {1}{\rho }}\left({\frac {\partial (\rho A_{\phi })}{\partial \rho }}-{\frac {\partial A_{\rho }}{\partial \phi }}\right){\hat {\boldsymbol {z}}}\end{aligned}}} div ⁡ A = lim V → 0 ∬ ∂ V A ⋅ d S ∭ V d V = A r ( r + d r ) ( r + d r ) d θ ( r + d r ) sin ⁡ θ d ϕ − A r ( r ) r d θ r sin ⁡ θ d ϕ + A θ ( θ + d θ ) sin ⁡ ( θ + d θ ) r d r d ϕ − A θ ( θ ) sin ⁡ ( θ ) r d r d ϕ + A ϕ ( ϕ + d ϕ ) ( r + d r / 2 ) d r d θ − A ϕ ( ϕ ) ( r + d r / 2 ) d r d θ d r r d θ r sin ⁡ θ d ϕ = 1 r 2 ∂ ( r 2 A r ) ∂ r + 1 r sin ⁡ θ ∂ ( A θ sin ⁡ θ ) ∂ θ + 1 r sin ⁡ θ ∂ A ϕ ∂ ϕ {\displaystyle {\begin{aligned}\operatorname {div} \mathbf {A} &=\lim _{V\to 0}{\frac {\iint _{\partial V}\mathbf {A} \cdot d\mathbf {S} }{\iiint _{V}dV}}\\&={\frac {A_{r}(r+dr)(r+dr)d\theta \,(r+dr)\sin \theta d\phi -A_{r}(r)rd\theta \,r\sin \theta d\phi +A_{\theta }(\theta +d\theta )\sin(\theta +d\theta )\,rdrd\phi -A_{\theta }(\theta )\sin(\theta )\,rdrd\phi +A_{\phi }(\phi +d\phi )(r+dr/2)drd\theta -A_{\phi }(\phi )(r+dr/2)drd\theta }{dr\,rd\theta \,r\sin \theta d\phi }}\\&={\frac {1}{r^{2}}}{\frac {\partial (r^{2}A_{r})}{\partial r}}+{\frac {1}{r\sin \theta }}{\frac {\partial (A_{\theta }\sin \theta )}{\partial \theta }}+{\frac {1}{r\sin \theta }}{\frac {\partial A_{\phi }}{\partial \phi }}\end{aligned}}} ( curl ⁡ A ) r = lim S ⊥ r ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A θ ( ϕ ) r d θ + A ϕ ( θ + d θ ) r sin ⁡ ( θ + d θ ) d ϕ − A θ ( ϕ + d ϕ ) r d θ − A ϕ ( θ ) r sin ⁡ ( θ ) d ϕ r d θ r sin ⁡ θ d ϕ = 1 r sin ⁡ θ ∂ ( A ϕ sin ⁡ θ ) ∂ θ − 1 r sin ⁡ θ ∂ A θ ∂ ϕ {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{r}=\lim _{S^{\perp {\boldsymbol {\hat {r}}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}&={\frac {A_{\theta }(\phi )\,rd\theta +A_{\phi }(\theta +d\theta )\,r\sin(\theta +d\theta )d\phi -A_{\theta }(\phi +d\phi )\,rd\theta -A_{\phi }(\theta )\,r\sin(\theta )d\phi }{rd\theta \,r\sin \theta d\phi }}\\&={\frac {1}{r\sin \theta }}{\frac {\partial (A_{\phi }\sin \theta )}{\partial \theta }}-{\frac {1}{r\sin \theta }}{\frac {\partial A_{\theta }}{\partial \phi }}\end{aligned}}} ( curl ⁡ A ) θ = lim S ⊥ θ ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A ϕ ( r ) r sin ⁡ θ d ϕ + A r ( ϕ + d ϕ ) d r − A ϕ ( r + d r ) ( r + d r ) sin ⁡ θ d ϕ − A r ( ϕ ) d r d r r sin ⁡ θ d ϕ = 1 r sin ⁡ θ ∂ A r ∂ ϕ − 1 r ∂ ( r A ϕ ) ∂ r {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{\theta }=\lim _{S^{\perp {\boldsymbol {\hat {\theta }}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}&={\frac {A_{\phi }(r)\,r\sin \theta d\phi +A_{r}(\phi +d\phi )dr-A_{\phi }(r+dr)(r+dr)\sin \theta d\phi -A_{r}(\phi )dr}{dr\,r\sin \theta d\phi }}\\&={\frac {1}{r\sin \theta }}{\frac {\partial A_{r}}{\partial \phi }}-{\frac {1}{r}}{\frac {\partial (rA_{\phi })}{\partial r}}\end{aligned}}} ( curl ⁡ A ) ϕ = lim S ⊥ ϕ ^ → 0 ∫ ∂ S A ⋅ d ℓ ∬ S d S = A r ( θ ) d r + A θ ( r + d r ) ( r + d r ) d θ − A r ( θ + d θ ) d r − A θ ( r ) r d θ ( r + d r / 2 ) d r d θ = 1 r ∂ ( r A θ ) ∂ r − 1 r ∂ A r ∂ θ {\displaystyle {\begin{aligned}(\operatorname {curl} \mathbf {A} )_{\phi }=\lim _{S^{\perp {\boldsymbol {\hat {\phi }}}}\to 0}{\frac {\int _{\partial S}\mathbf {A} \cdot d\mathbf {\ell } }{\iint _{S}dS}}&={\frac {A_{r}(\theta )dr+A_{\theta }(r+dr)(r+dr)d\theta -A_{r}(\theta +d\theta )dr-A_{\theta }(r)\,rd\theta }{(r+dr/2)drd\theta }}\\&={\frac {1}{r}}{\frac {\partial (rA_{\theta })}{\partial r}}-{\frac {1}{r}}{\frac {\partial A_{r}}{\partial \theta }}\end{aligned}}} curl ⁡ A = ( curl ⁡ A ) r r ^ + ( curl ⁡ A ) θ θ ^ + ( curl ⁡ A ) ϕ ϕ ^ = 1 r sin ⁡ θ ( ∂ ( A ϕ sin ⁡ θ ) ∂ θ − ∂ A θ ∂ ϕ ) r ^ + 1 r ( 1 sin ⁡ θ ∂ A r ∂ ϕ − ∂ ( r A ϕ ) ∂ r ) θ ^ + 1 r ( ∂ ( r A θ ) ∂ r − ∂ A r ∂ θ ) ϕ ^ {\displaystyle \operatorname {curl} \mathbf {A} =(\operatorname {curl} \mathbf {A} )_{r}\,{\hat {\boldsymbol {r}}}+(\operatorname {curl} \mathbf {A} )_{\theta }\,{\hat {\boldsymbol {\theta }}}+(\operatorname {curl} \mathbf {A} )_{\phi }\,{\hat {\boldsymbol {\phi }}}={\frac {1}{r\sin \theta }}\left({\frac {\partial (A_{\phi }\sin \theta )}{\partial \theta }}-{\frac {\partial A_{\theta }}{\partial \phi }}\right){\hat {\boldsymbol {r}}}+{\frac {1}{r}}\left({\frac {1}{\sin \theta }}{\frac {\partial A_{r}}{\partial \phi }}-{\frac {\partial (rA_{\phi })}{\partial r}}\right){\hat {\boldsymbol {\theta }}}+{\frac {1}{r}}\left({\frac {\partial (rA_{\theta })}{\partial r}}-{\frac {\partial A_{r}}{\partial \theta }}\right){\hat {\boldsymbol {\phi }}}} 坐标参数u的单位向量以如下方式定义，u的小的正值改变导致位置向量 r → {\displaystyle {\boldsymbol {\vec {r}}}} 在 u ^ {\displaystyle {\boldsymbol {\hat {u}}}} 方向上的改变。因此： ∂ r → ∂ u = ∂ s ∂ u u ^ {\displaystyle {\partial {\boldsymbol {\vec {r}}} \over \partial u}={\partial {s} \over \partial u}{\boldsymbol {\hat {u}}}} 这里的s是弧长参数。 对于两组坐标系 u i {\displaystyle u_{i}} 和 v j {\displaystyle v_{j}} ，依据链式法则： d r → = ∑ i ∂ r → ∂ u i d u i = ∑ i ∂ s ∂ u i u i ^ d u i = ∑ j ∂ s ∂ v j v j ^ d v j = ∑ j ∂ s ∂ v j v j ^ ∑ i ∂ v j ∂ u i d u i = ∑ i ∑ j ∂ s ∂ v j ∂ v j ∂ u i v j ^ d u i {\displaystyle d{\boldsymbol {\vec {r}}}=\sum _{i}{\partial {\boldsymbol {\vec {r}}} \over \partial u_{i}}du_{i}=\sum _{i}{\partial {s} \over \partial u_{i}}{\boldsymbol {\hat {u_{i}}}}du_{i}=\sum _{j}{\partial {s} \over \partial v_{j}}{\boldsymbol {\hat {v_{j}}}}dv_{j}=\sum _{j}{\partial {s} \over \partial v_{j}}{\boldsymbol {\hat {v_{j}}}}\sum _{i}{\partial {v_{j}} \over \partial u_{i}}du_{i}=\sum _{i}\sum _{j}{\partial {s} \over \partial v_{j}}{\partial {v_{j}} \over \partial u_{i}}{\boldsymbol {\hat {v_{j}}}}du_{i}} 现在，使除了一个之外的所有 d u i = 0 {\displaystyle du_{i}=0} 并在两边除以对应的坐标参数的微分，得到： ∂ s ∂ u i u i ^ = ∑ j ∂ s ∂ v j ∂ v j ∂ u i v j ^ {\displaystyle {\partial {s} \over \partial u_{i}}{\boldsymbol {\hat {u_{i}}}}=\sum _{j}{\partial {s} \over \partial v_{j}}{\partial {v_{j}} \over \partial u_{i}}{\boldsymbol {\hat {v_{j}}}}} Del算子 正交坐标系 曲线坐标系 在圆柱和球坐标中的向量场 Griffiths, David J. Introduction to Electrodynamics. Pearson. 2012. ISBN 978-0-321-85656-2.  Weisstein, Eric W. (编). Convective Operator. at MathWorld--A Wolfram Web Resource. Wolfram Research, Inc. [23 March 2011]. （原始内容存档于2016-03-03） （英语）.  Maxima Computer Algebra system scripts（页面存档备份，存于互联网档案馆） to generate some of these operators in cylindrical and spherical coordinates. Curl, Divergence, and Gradient in Cylindrical and Spherical Coordinate Systems.