測地線

測地学
基本
  • 測地学
  • 地理力学(英語版)
  • ジオマティクス(英語版)
  • 地図学
  • 測地学の歴史(英語版)
概念
技術
基準(歴史(英語版)
NGVD29(英語版) 海面測地系1929年
OSGB36(英語版) イギリス陸上測量1936年
SK-42(英語版) Systema Koordinat 1942 goda
ED50(英語版) 欧州座標系1950年
SAD69(英語版) 南米測地系1969年
GRS80GRS80地球楕円体1980年
NAD83 北米測地系1983年
WGS84 世界測地系1984年
NAVD88(英語版) 北米垂直測地系1988年
ETRS89(英語版) 欧州地球基準座標系システム1989年
GCJ-02 中国の暗号化された座標系2002年
Geo URI(英語版) 地点へのインターネットリンク 2010年
或る測地線の球面においてはその大円が在る。赤い線でなぞった測地線の跡は、点PとQでの採り得る最短距離を表す。対蹠点 uとvでは、幾つかの測地線が最短距離となる。

微分幾何学において測地線(そくちせん、: geodesic)とは、曲面(より一般的にはリーマン多様体)上の曲線であって、その上の十分近い2つの離れた点が最短線で結ばれた曲線を言う。ユークリッド空間における直線の概念を、曲がった空間において一般化したものである。「測地線」という用語は、地球の大きさと形状を測定する学問である測地学に由来する。本来の意味では、測地線は地表の2点間の最短ルートであり、球体形状の地球の場合、大円の一部となる。測地線の中でその長さが最小のものは最短測地線という。

リーマン空間において、ある曲線が曲面上の測地線となるための必要十分条件は、曲線の主法線と曲面の接平面の法線とが曲線に沿って常に一致することである[1]

この概念は、数学的な空間にも拡張され[注釈 1]、例えばグラフ理論ではグラフ上の2つの頂点 (vertex) や結節点 (node) 間の測地線が定義されている。一般相対性理論では、光は曲がった空間での測地線を進むという原理に基づいて構築されている。

概要

1697年ヨハン・ベルヌーイは、曲面上の2点をその曲面上で結び、その長さを最小にする問題を考え、このような長さの最も短い曲線を、この曲面上の測地線(geodesic)と呼んだ[2]。そして、測地線上の点での接平面の法線がその点で曲面に垂直であることを発見した。1698年ヤコブ・ベルヌーイ円筒円錐回転面上の測地線を求めた。1728年レオンハルト・オイラーは自身の開発した変分法を用いて、曲面上の測地線が満たす微分方程式を導出した。

典型的な測地線は、測地学の対象でもある地球上の2点を結ぶ最短曲線である。地球を単純に球面であるとする。例えば、東京ニューヨークの間を最短距離で移動するためには、東京とニューヨークを通る大円に沿った移動を行えばよく、この大円の一部こそ、測地線と呼ばれるものになる[注釈 2]

2点間の最短距離を示す曲線は測地線となるので、2点を結ぶ測地線の中で最短のものが2点の最短距離を示すと考えてよい。その意味で、測地線というのは、2点間の最短距離を測るための曲線の候補の集まりであるともいえる。ちなみに、2点を北極南極のような対極の位置に取れば、この2点を結ぶ最短測地線は無数にあることにも注意されたい。

球面では測地線は閉曲線となるが、回転楕円体面上など一般には測地線は閉曲線とならない。

測地線の方程式

計量テンソル g i j {\displaystyle g_{ij}} を持つリーマン多様体上の微分可能な曲線 x ( t ) {\displaystyle x(t)_{}^{}} のある一点 x ( t 1 ) {\displaystyle x(t_{1})} から他の一点 x ( t 2 ) {\displaystyle x(t_{2})} までの長さ S は積分

S ( = s 1 s 2 d s = t 1 t 2 d s d t d t ) = t 1 t 2 g i j d x i d t d x j d t d t = t 1 t 2 F ( x , x ˙ ) d t {\displaystyle S\left(=\int _{s_{1}}^{s_{2}}\mathrm {d} s=\int _{t_{1}}^{t_{2}}{\frac {\mathrm {d} s}{\mathrm {d} t}}\mathrm {d} t\right)=\int _{t_{1}}^{t_{2}}{\sqrt {g_{ij}{\frac {\mathrm {d} x^{i}}{\mathrm {d} t}}{\frac {\mathrm {d} x^{j}}{\mathrm {d} t}}}}\mathrm {d} t=\int _{t_{1}}^{t_{2}}F(x,{\dot {x}})\mathrm {d} t}

で与えられる。この変分 δ S {\displaystyle \delta S} について δ S = 0 {\displaystyle \delta S=0} となる曲線 x ( t ) {\displaystyle x(t)} をそのリーマン多様体の測地線と呼ぶ。この曲線 x ( t ) {\displaystyle x(t)} について δ S = 0 {\displaystyle \delta S=0} となるための必要十分条件は曲線 x ( t ) {\displaystyle x(t)} オイラー=ラグランジュ方程式

d d t ( F x ˙ a ) F x a = 0 {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial F}{\partial {\dot {x}}^{a}}}\right)-{\frac {\partial F}{\partial x^{a}}}=0}  ただし、 F = g i j x ˙ i x ˙ j {\displaystyle F={\sqrt {g_{ij}{\dot {x}}^{i}{\dot {x}}^{j}}}}

を満たすことである。

d d t ( F x ˙ a ) F x a = d d t ( g a i x ˙ i d s d t ) 1 2 g i j x a x ˙ i x ˙ j d s d t = g a i x ¨ i s ˙ + g a i x j x ˙ i x ˙ j s ˙ g a i x ˙ i s ¨ ( s ˙ ) 2 1 2 g i j x a x ˙ i x ˙ j s ˙ = 0 {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial F}{\partial {\dot {x}}^{a}}}\right)-{\frac {\partial F}{\partial x^{a}}}={\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {g_{ai}{\dot {x}}^{i}}{\frac {\mathrm {d} s}{\mathrm {d} t}}}\right)-{\frac {1}{2}}{\frac {{\frac {\partial g_{ij}}{\partial x^{a}}}{\dot {x}}^{i}{\dot {x}}^{j}}{\frac {\mathrm {d} s}{\mathrm {d} t}}}={\frac {g_{ai}{\ddot {x}}^{i}}{\dot {s}}}+{\frac {{\frac {\partial g_{ai}}{\partial x^{j}}}{\dot {x}}^{i}{\dot {x}}^{j}}{\dot {s}}}-{\frac {g_{ai}{\dot {x}}^{i}{\ddot {s}}}{\left({\dot {s}}\right)^{2}}}-{\frac {1}{2}}{\frac {{\frac {\partial g_{ij}}{\partial x^{a}}}{\dot {x}}^{i}{\dot {x}}^{j}}{\dot {s}}}=0}

整理すると、

g a i ( x ¨ i x ˙ i s ¨ s ˙ ) + g a i x j x ˙ i x ˙ j 1 2 g i j x a x ˙ i x ˙ j = g a i ( x ¨ i x ˙ i s ¨ s ˙ ) + 1 2 ( g a i x j + g a j x i g i j x a ) x ˙ i x ˙ j = 0 {\displaystyle g_{ai}\left({\ddot {x}}^{i}-{\dot {x}}^{i}{\frac {\ddot {s}}{\dot {s}}}\right)+{\frac {\partial g_{ai}}{\partial x^{j}}}{\dot {x}}^{i}{\dot {x}}^{j}-{\frac {1}{2}}{\frac {\partial g_{ij}}{\partial x^{a}}}{\dot {x}}^{i}{\dot {x}}^{j}=g_{ai}\left({\ddot {x}}^{i}-{\dot {x}}^{i}{\frac {\ddot {s}}{\dot {s}}}\right)+{\frac {1}{2}}\left({\frac {\partial g_{ai}}{\partial x^{j}}}+{\frac {\partial g_{aj}}{\partial x^{i}}}-{\frac {\partial g_{ij}}{\partial x^{a}}}\right){\dot {x}}^{i}{\dot {x}}^{j}=0}

となる。上式に g k a {\displaystyle g^{ka}} を掛け合わせ a に関して足しあわせを行うと、

δ i k ( x ¨ i x ˙ i s ¨ s ˙ ) + 1 2 g k a ( g a i x j + g a j x i g i j x a ) x ˙ i x ˙ j = x ¨ k x ˙ k s ¨ s ˙ + 1 2 g k a ( g a i x j + g a j x i g i j x a ) x ˙ i x ˙ j = 0 {\displaystyle \delta _{i}^{k}\left({\ddot {x}}^{i}-{\dot {x}}^{i}{\frac {\ddot {s}}{\dot {s}}}\right)+{\frac {1}{2}}g^{ka}\left({\frac {\partial g_{ai}}{\partial x^{j}}}+{\frac {\partial g_{aj}}{\partial x^{i}}}-{\frac {\partial g_{ij}}{\partial x^{a}}}\right){\dot {x}}^{i}{\dot {x}}^{j}={\ddot {x}}^{k}-{\dot {x}}^{k}{\frac {\ddot {s}}{\dot {s}}}+{\frac {1}{2}}g^{ka}\left({\frac {\partial g_{ai}}{\partial x^{j}}}+{\frac {\partial g_{aj}}{\partial x^{i}}}-{\frac {\partial g_{ij}}{\partial x^{a}}}\right){\dot {x}}^{i}{\dot {x}}^{j}=0}

となるが、ここで弧長 s を媒介変数一次関数と置き直せば s ¨ = 0 {\displaystyle {\ddot {s}}=0} となる。特に媒介変数 t を弧長 s と置き直せばより簡単になり、

d 2 x k d s 2 + 1 2 g k a ( g a i x j + g a j x i g i j x a ) d x i d s d x j d s = 0 {\displaystyle {\frac {\mathrm {d} ^{2}x^{k}}{\mathrm {d} s^{2}}}+{\frac {1}{2}}g^{ka}\left({\frac {\partial g_{ai}}{\partial x^{j}}}+{\frac {\partial g_{aj}}{\partial x^{i}}}-{\frac {\partial g_{ij}}{\partial x^{a}}}\right){\frac {\mathrm {d} x^{i}}{\mathrm {d} s}}{\frac {\mathrm {d} x^{j}}{\mathrm {d} s}}=0}

を得る[3]。最後にクリストッフェル記号

{ k i j } = 1 2 g k a ( g a i x j + g a j x i g i j x a ) {\displaystyle \left\{{{k} \atop {ij}}\right\}={\frac {1}{2}}g^{ka}\left({\frac {\partial g_{ai}}{\partial x^{j}}}+{\frac {\partial g_{aj}}{\partial x^{i}}}-{\frac {\partial g_{ij}}{\partial x^{a}}}\right)}

で置き直せば、上式は、

d 2 x k d s 2 + { k i j } d x i d s d x j d s = 0 {\displaystyle {\frac {\mathrm {d} ^{2}x^{k}}{\mathrm {d} s^{2}}}+\left\{{{k} \atop {ij}}\right\}{\frac {\mathrm {d} x^{i}}{\mathrm {d} s}}{\frac {\mathrm {d} x^{j}}{\mathrm {d} s}}=0}

と表されることになる。これを測地線の方程式(geodesic equation)という。 s {\displaystyle s} x ( s ) {\displaystyle x(s)} の始点からの長さを表す弧長パラメータである。たとえば、3次元の空間が平坦である場合は、 g μ ν = d i a g ( 1 , 1 , 1 ) {\displaystyle g_{\mu \nu }^{}={\rm {diag}}(1,1,1)} であり、接続はすべて0となる為、測地線の方程式は単に d 2 x i d s 2 = 0 {\displaystyle {\frac {d^{2}x^{i}}{ds^{2}}}=0} となる。つまり、 x i {\displaystyle x^{i}} s {\displaystyle s} の1次式であり 通常の直線の方程式を表すこととなる。

この方程式は、最短測地線の満たすべき「2点間の最短距離を示す」という性質、或いは、「測地線 x ( τ ) {\displaystyle x(\tau )} 接ベクトル場レヴィ-チビタ接続(リーマン接続)に関して平行である」という性質により導かれる。

繰り返しになるが、微分方程式は局所的な情報を与えるものなので、大域的な曲線の長さなどを表すものではなく、したがってこの方程式で定義される測地線が必ずしも最短測地線とはならないことに注意されたい。

応用例

一般相対性理論では時空を4次元の擬リーマン多様体として記述する。時空上の試験粒子(時空への重力的な反作用を与えない仮想的な質点である。電荷やスピンなどの性質は通常持たないと考える)や光の経路は測地線で記述されると考えられている。いわゆる自由落下している物体の軌跡は測地線で表されると考えるのである。たとえば、地上でボールを放り投げたときに描く放物線も4次元の時空の中でその軌跡を捉えれば測地線である。一般相対性理論では測地線は時空の因果構造を定義するときに重要な役割を果たす。ブラックホールの定義や特異点定理、そのほか数学的な時空の定式化には欠かせない道具である。

回転楕円体面上の測地線

詳細は「:en:Geodesics on an ellipsoid」を参照

回転楕円体面上の測地線は、地球の場合に大圏コースに対応する。経線に沿う測地線は子午線弧

脚注

[脚注の使い方]

注釈

  1. ^ 測地線や極小曲面の概念をM次元の幾何学的対象に一般化するにはリーマン多様体で考える必要があろう(測地線は1次元リーマン多様体であり、極小曲面は2次元リーマン多様体である)。その際これら測地線の(両)端点や極小曲面の縁の曲線(あるいは端点)は、それら対象となっている多様体から変形運動するそれら多様体が置かれるN次元の空間である多様体(たとえば球面上の測地線で考えるならばその次元は2である)との写像に関する、変分問題の境界条件として捉え直される(西川 2006, pp. 89–124、p. 105 図3.1。'多様体間の調和写像'の項をも見よ。) 。
  2. ^ しかしながら、一般に、大円をその上の2点で分けると優弧と劣弧に分かれる。東京からニューヨークへ大円に沿った移動をしても、東京からニューヨークに行くには大円の周り方によって遠い移動と近い移動とある。この場合、劣弧に沿って移動すれば最短距離、優弧に沿えば直線的な移動としては最も遠回りになるわけである。大円の一部である弧は測地線となるが、必ずしも2点間の最短距離を示す曲線とはならない。

出典

  1. ^ 矢野 1949, p. 58
  2. ^ 近藤 基吉、井関 清志『現代数学の黎明期 近代数学[上]』日本評論社、1982年。  p.275
  3. ^ 矢野 1949, pp. 120–121

参考文献

  • 矢野, 健太郎『微分幾何学』朝倉書店、1949年。 
  • リーマン、リッチ、レビ=チビタ、アインシュタイン、マイヤー 著、矢野 健太郎(訳) 編『リーマン幾何とその応用』 10巻、共立出版株式会社〈現代数学の系譜〉、1971年。 
  • 西川青季『幾何学的変分問題』岩波書店、2006年4月5日。ISBN 4-00-005243-8。 

関連項目

Glossary of tensor theory(英語版)
範囲 (Scope)
数学
物理学 • 工学
表記法 (Notation)
テンソルの定義
算法
関連事項
有名なテンソル
数学
物理学
数学者
カテゴリカテゴリ
スタブアイコン

この項目は、微分幾何学に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています(プロジェクト:数学/Portal:数学)。

  • 表示
  • 編集
典拠管理データベース: 国立図書館 ウィキデータを編集
  • ドイツ
  • イスラエル
  • アメリカ