正交座標系

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA數學裏,一個正交坐標系定義為一組正交坐標𝐪=(q1, q2, q3,  qn),其坐標曲面都以直角相交(注意:很多作者采用爱因斯坦记号对坐标标号使用上标并非表示指数)。坐標曲面定義為特定坐標qi等值曲面,即qi為常数的曲线曲面超曲面。例如,三維直角坐標(x, y, z)是一種正交坐標系,它的x為常數,y為常數,z為常數的坐標曲面,都是互相以直角相交的平面,都互相垂直。正交坐标系是曲线坐标系的特殊的但极其常见的形式。

动机

正交座標時常用來解析一些出現於量子力學流體動力學電動力學熱力學等等的偏微分方程。舉例而言,選擇一個恰當的的正交座標來解析氫離子H2波函數或消防水管的噴水,也許會比用直角座標方便的多。這主要是因為恰當的正交座標能夠與一個問題的對稱性相配合,從而促使應用分離變數法來成功的解析關於這問題的方程式。分離變數法是一種數學技巧,專門用來將一個複雜的n維問題變為n個一維問題。很多問題都可以簡化為拉普拉斯方程亥姆霍茲方程,這些方程式可以用很多種正交座標來分離。拉普拉斯方程可以在13个正交坐标系中分离(本文列出的14个中圆环坐标系除外),而亥姆霍茲方程可以在11个正交坐标系中分离[1][2]

概述

共形映射作用于矩形网格。注意,弯曲的网格的正交性被保留。

正交坐標的度規張量絕對沒有非對角項目。換句話說,無窮小距離的平方ds2,可以寫為無窮小坐標位移的平方和:

ds2=i=1n(hidqi)2

其中,n是維數,標度因子hi是度規張量的對角元素gii的平方根:

hi(𝐪) =def gii(𝐪)

這些標度因子可以用來計算一個正交坐標系的微分算子。例如,梯度拉普拉斯算子散度、或旋度

在數學裏,存在有各種各樣的正交座標系。應用二維直角座標系(x, y)共形映射方法,可以簡易的生成這些正交座標系。一個複數z=x+iy的任何全純函數w=f(z),其複值的導數,如果不等於零,則會造成一個共形映射。如果答案可以表達為w=u+iv,則uv的等值曲線以直角相交,就如同原本的xy的等值曲線以直角相交。

三維與更高維的正交座標系可以由一個二維正交座標系生成,只要將二維正交座標往一個新的座標軸投射(形成類似圓柱座標系的座標系),或者將二維正交座標繞著其對稱軸旋轉。可是,也有一些三維正交座標系,例如橢球座標系,則不能夠用上述方法得到。更一般的正交坐标可以从一些必要的坐标曲面/曲线起步并通过考虑它们的Template:En-link而得到。

向量代数

在正交坐標系裏,內積的公式仍舊不變:

𝐀𝐁=i=1nAiBi

向量微積分

從前面的距離公式,可以觀察出,一個正交坐標qi的無窮小改變dqi,其相伴的長度是dsi=hidqi。因此,一個位移向量的全微分d𝐫等於

d𝐫=i=1nhidqi𝐞i

其中,𝐞i是垂直於qi等值曲面的單位向量,指向著qi增值最快的方向,這些單位向量形成了一個局部直角坐標系的坐標軸。

因此,向量𝐅沿著周線的線積分等於

𝐅d𝐫=i=1nFihidqi

其中,Fi是向量𝐅在單位向量𝐞i方向的分量:

Fi =def 𝐞i𝐅

類似地,一個無窮小面積元素是

dA=dsidsj=hihjdqidqj,ij

一個無窮小體積元素是

dV=dsidsjdsk=hihjhkdqidqjdqk,ijk

例如,向量𝐅對於一個曲面𝕊的曲面積分是

𝕊𝐅d𝐀=𝕊F1h2h3dq2dq3+𝕊F2h3h1dq3dq1+𝕊F3h1h2dq1dq2

球坐標系實例

直角坐標(x, y, z)與球坐標(r, θ,ϕ)的變換方程式為

x=rsinθcosϕ
y=rsinθsinϕ
z=rcosθ

直角坐標的全微分是

dx=sinθcosϕdr+rcosθcosϕdθrsinθsinϕdϕ
dy=sinθsinϕdr+rcosθsinϕdθ+rsinθcosϕdϕ
dz=cosθdrrsinθdθ

所以,無窮小距離的平方是

ds2=dx2+dy2+dz2=dr2+(rdθ)2+(rsinθdϕ)2

標度因子是

hr=1
hθ=r
hϕ=rsinθ

向量𝐅沿著周線的線積分等於

𝐅d𝐫=Fr dr+Fθ rdθ+Fϕ rsinθdϕ

向量𝐅對於一個曲面𝕊的曲面積分是

𝕊𝐅d𝐀=𝕊Fr r2sinθdθdϕ+𝕊Fθ rsinθdrdϕ+𝕊Fϕ rdrdθ

三維微分算子

Template:Main

算子 正交坐標公式
标量场梯度 Φ=𝐞^11h1Φq1+𝐞^21h2Φq2+𝐞^31h3Φq3
向量场散度 𝐅=1h1h2h3[q1(F1h2h3)+q2(F2h3h1)+q3(F3h1h2)]
向量场旋度 ×𝐅=𝐞1h2h3[q2(h3F3)q3(h2F2)]+𝐞2h3h1[q3(h1F1)q1(h3F3)]+𝐞3h1h2[q1(h2F2)q2(h1F1)]
标量场拉普拉斯算子 2Φ=1h1h2h3[q1(h2h3h1Φq1)+q2(h3h1h2Φq2)+q3(h1h2h3Φq3)]

上面表达式可以使用列維-奇維塔符號ϵ的更简洁形式书写,定义H=h1h2h3,并使用爱因斯坦记号,即在同时出现上标和下标的项目上求此项所有可能的总和:

算子 表达式
标量场梯度 ϕ=𝐞^khkϕqk
向量场散度 𝐅=1Hqk(HhkFk)
向量场(只3D)的旋度 (×𝐅)k=hk𝐞^kHϵijkqi(hjFj)
标量场的拉普拉斯算子 2ϕ=1Hqk(Hhk2ϕqk)

二维正交坐标系表格

坐标系 复数变换

x+iy=f(u+iv)

uv等值线的形状 注释
直角 u+iv 直线, 直线
Template:En-link exp(u+iv) 圆, 直线 u=lnr则为极坐标系
抛物线 12(u+iv)2 抛物线, 抛物线
点偶极 (u+iv)1 圆, 圆
椭圆 cosh(u+iv) 椭圆, 双曲线 对于大距离看似对数极
双极 coth(u+iv) 圆, 圆 对于大距离看似点偶极
u+iv 双曲线, 双曲线
u=x2+2y2, y=vx2 椭圆, 抛物线

Template:Multiple image

三维正交坐标系表格

除了直角坐标系之外,下表列出其他常见的正交坐标系[3],为了简明性在坐标列中使用了区间符号

曲线坐标 (q1, q2, q3) 从直角坐标(x, y, z)转换 缩放因子
球极坐标系

(r,θ,ϕ)[0,)×[0,π]×[0,2π)

x=rsinθcosϕy=rsinθsinϕz=rcosθ h1=1h2=rh3=rsinθ
圆柱坐标系

(ρ,ϕ,z)[0,)×[0,2π)×(,)

x=ρcosϕy=ρsinϕz=z h1=h3=1h2=ρ
抛物柱面坐标系

(u,v,z)(,)×[0,)×(,)

x=12(u2v2)y=uvz=z h1=h2=u2+v2h3=1
抛物线坐标系

(u,v,ϕ)[0,)×[0,)×[0,2π)

x=uvcosϕy=uvsinϕz=12(u2v2) h1=h2=u2+v2h3=uv
椭圆柱坐标系

(u,v,z)[0,)×[0,2π)×(,)

x=acoshucosvy=asinhusinvz=z h1=h2=asinh2u+sin2vh3=1
长球面坐标系

(ξ,η,ϕ)[0,)×[0,π]×[0,2π)

x=asinhξsinηcosϕy=asinhξsinηsinϕz=acoshξcosη h1=h2=asinh2ξ+sin2ηh3=asinhξsinη
扁球面坐标系

(ξ,η,ϕ)[0,)×[π2,π2]×[0,2π)

x=acoshξcosηcosϕy=acoshξcosηsinϕz=asinhξsinη h1=h2=asinh2ξ+sin2ηh3=acoshξcosη
双极圆柱坐标系

(u,v,z)[0,2π)×(,)×(,)

x=asinhvcoshvcosuy=asinucoshvcosuz=z h1=h2=acoshvcosuh3=1
圆环坐标系

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinhvcosϕcoshvcosuy=asinhvsinϕcoshvcosuz=asinucoshvcosu h1=h2=acoshvcosuh3=asinhvcoshvcosu
双球坐标系

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinucosϕcoshvcosuy=asinusinϕcoshvcosuz=asinhvcoshvcosu h1=h2=acoshvcosuh3=asinucoshvcosu
圆锥坐标系

(λ,μ,ν)ν2<b2<μ2<a2λ[0,)

x=λμνaby=λa(μ2a2)(ν2a2)a2b2z=λb(μ2b2)(ν2b2)a2b2 h1=1h22=λ2(μ2ν2)(μ2a2)(b2μ2)h32=λ2(μ2ν2)(ν2a2)(ν2b2)
抛物面坐标系

(λ,μ,ν)λ<b2<μ<a2<ν

x2qia2+y2qib2=2z+qi

其中 (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)
椭球坐标系

(λ,μ,ν)λ<c2<b2<a2,c2<μ<b2<a2,c2<b2<ν<a2,

x2a2qi+y2b2qi+z2c2qi=1

其中 (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)(c2qi)

微分算子导引

梯度導引

一個函數ϕ的梯度朝某個方向𝐧^的分量,等於方向導數 dϕds𝐧^方向的值:

Φ𝐧^=dϕds

其中,ds是朝𝐧^方向的無窮小位移。

假若,這𝐧^與正交坐標軸𝐞^i同方向。那麼,ds=hidqi。所以,函數ϕ的梯度朝𝐞^i的分量是ϕhiqi;也就是說,

Φ=𝐞^11h1Φq1+𝐞^21h2Φq2+𝐞^31h3Φq3

散度導引

𝐅=(𝐞^1F1+𝐞^2F2+𝐞^3F3)

取右手邊第一個項目,

(𝐞^1F1)=[(𝐞^1h2h3)(h2h3F1)]

應用向量恆等式(𝐀ϕ)=ϕ𝐀+𝐀(ϕ)(ϕ1×ϕ2)=0,可以得到

(𝐞^1F1)=(h2h3F1)(𝐞^1h2h3)+(𝐞^1h2h3)(h2h3F1)=(h2h3F1)[(q2)×(q3)]+(𝐞^1h2h3)(h2h3F1)=(𝐞^1h2h3)(h2h3F1)=1h1h2h3q1(F1h2h3)

總合所有項目,

𝐅=1h1h2h3[q1(F1h2h3)+q2(F2h3h1)+q3(F3h1h2)]


旋度導引

×𝐅=×(𝐞^1F1+𝐞^2F2+𝐞^3F3)

取右手邊第一個項目,

×(𝐞^1F1)=×[(𝐞^1h1)(h1F1)]

應用向量恆等式×(𝐀ϕ)=ϕ×𝐀𝐀×(ϕ)

×(𝐞^1F1)=(h1F1)×(𝐞^1h1)(𝐞^1h1)×(h1F1)=(h1F1)×(q1)(𝐞^1h1)×(𝐞^2h2q2(h1F1)+𝐞^3h3q3(h1F1))

應用向量恆等式×(ϕ)=0

×(𝐞^1F1)=𝐞^2h1h3q3(h1F1)𝐞^3h1h2q2(h1F1)

總合所有項目,

×𝐅=𝐞1h2h3[q2(h3F3)q3(h2F2)]+𝐞2h3h1[q3(h1F1)q1(h3F3)]+𝐞3h1h2[q1(h2F2)q2(h1F1)]

拉普拉斯算子

2Φ=1h1h2h3[q1(h2h3h1Φq1)+q2(h3h1h2Φq2)+q3(h1h2h3Φq3)]

引用

Template:Reflist

参见

參考文獻

  • Korn GA and Korn TM. (1961) Mathematical Handbook for Scientists and Engineers, McGraw-Hill, pp. 164-182。
  • Morse PM and Feshbach H. (1953) Methods of Theoretical Physics, McGraw-Hill, pp. 494-523, 655-666。
  • Margenau H. and Murphy GM. (1956) The Mathematics of Physics and Chemistry, 2nd. ed., Van Nostrand, pp. 172-192。

Template:正交坐標系

  1. Template:Cite web
  2. Template:Harvnb
  3. Vector Analysis (2nd Edition), M.R. Spiegel, S. Lipschutz, D. Spellman, Schaum’s Outlines, McGraw Hill (USA), 2009, Template:Isbn