等時曲線

トートクロンカーブを表すオブジェクト

サイクロイド曲線を異なる地点から出発して滑り降りる4つの質点は同時に最下点へ到達する。青い矢印は曲線に沿った質点の加速度を示す。右上には時間-位置グラフを示す。

等時曲線（とうじきょくせん）または等時降下曲線（とうじこうかきょくせん）とは、物体が一様重力場の下でその曲線に沿って摩擦なく滑り降りるとき、最下点に達するまでの時間が出発点に依存しなくなるような曲線をいう。英語では テンプレート:Lang または テンプレート:Lang（ギリシャ語の接頭辞 テンプレート:ラテン翻字 「同じ」、 テンプレート:ラテン翻字「等しい」、テンプレート:ラテン翻字「時間の」から）と呼ばれる。この曲線はサイクロイドであり、また降下時間はサイクロイドの動円の半径を重力加速度で除したものの平方根に テンプレート:Mvar を乗じたものとなる。また、等時曲線は全ての出発点について最速降下曲線と同一である。

等時曲線問題

等時曲線問題とは、等時曲線がどのような曲線であるかを同定する問題である。この問題は1659年にクリスティアーン・ホイヘンスによって解かれた。彼は1673年初出の自著 テンプレート:Lang 中で、等時曲線がサイクロイドであることを幾何学的に証明している。

軸を鉛直とし、頂点を最下点としたサイクロイドの降下時間、すなわち物体が最下点に到達するまでの時間は、サイクロイド上のどの点から出発したとしても互いに等しい...

ホイヘンスは同時に、降下時間は物体がサイクロイドの動円半径と等しい距離を自由落下するのにかかる時間に テンプレート:Math を乗じたものに等しいことを証明した。現代的記法を用いれば、降下時間は テンプレート:Mvar を動円半径、テンプレート:Mvar を重力加速度として テンプレート:Math と表わされる。

振幅の異なる5つの等時サイクロイド振り子

その後、この解は最速降下曲線問題への取り組みに活用された。 ヤコブ・ベルヌーイは最速降下曲線問題を解析学を用いて解き、初めて積分の用語が用いられた論文 (テンプレート:Lang, 1690) において発表した。

サイクロイド振り子の図

ホイヘンスが等時曲線問題を詳しく研究する過程で、円軌道を描く振り子は厳密には等時性をもっておらず、したがって当時の振り子時計は振幅によって異なる時間を刻むことが明らかになった。正しい軌道を得たホイヘンスは、おもりを糸で吊し、糸の最上点の傍に障害物を設置しておもりに等時曲線を描かせることを試みた。これらの試みは多くの理由から実用的でなかった。まず、糸が曲がることにより摩擦が生じる。さらに、等時曲線から外れることよりも大きな誤差要因はいくつもあり、あまり改善に役立たない。その上、振り子の「円軌道誤差」は振幅が小さくなるにつれて減少するため、テンプレート:仮リンクの改善によりこの誤差は大きく低減できる。

その後、ジョゼフ＝ルイ・ラグランジュとレオンハルト・オイラーがこの問題を解析的に解いた。

ラグランジュの解

質点の位置を最下点からの弧長 テンプレート:Math を媒介変数として表示すると、 運動エネルギーは テンプレート:Math に比例する。ポテンシャルエネルギーは高さ テンプレート:Math に比例する。曲線が等時曲線となるのは、ラグランジアンが単純調和振動子のものとなる条件、すなわち、高さが弧長の2乗に等しい時である。長さの単位を比例係数を1となるようにとると

$y(s)=s^2$

と書ける。この関係を微分形式で書き下し、媒介変数を消去すると以下を得る。

$\mathrm{d}y=2s\mathrm{d}s$

$\mathrm{d}{y}^2=4s^2\mathrm{d}{s}^2=4y(\mathrm{d}{x}^2+\mathrm{d}{y}^2)$

これを解くため、 テンプレート:Math を テンプレート:Math で積分する。

$\frac{\mathrm{d}x}{\mathrm{d}y}=\frac{\sqrt{1-4y}}{2\sqrt{y}}$

$x=\int \sqrt{1-4u^2}\mathrm{d}u$

ここで、 $u=\sqrt{y}$ と置いた。この積分は円よりも下の面積であり、三角形と扇形に切り分けられる。

$x=\frac{1}{2}u\sqrt{1-4u^2}+\frac{1}{4}\arcsin 2u$

$y=u^2$

これが見慣れない媒介変数表示のサイクロイドであることを示すため、テンプレート:Math と置いて代数的に整理すると以下を得る。

$8x=2\sin \theta \cos \theta +2\theta =\sin 2\theta +2\theta$

$8y=2{\sin }^2\theta =1-\cos 2\theta$

これは、係数を除いてサイクロイド曲線の通常の媒介変数表示となっている。

「仮想重力」解

等時曲線問題の最も単純な解は、傾斜角とその傾斜角で質点が感じる重力の関係式を直接書き下すものである。傾斜角 90° の鉛直線上の質点は重力の影響を完全に受け、水平線上の質点は重力の影響を全く受けない。中間の傾斜角では質点の感じる「仮想重力」 テンプレート:Math を受けると考えることができる。まず、望ましい振る舞いを産み出す「仮想重力」はどんなものかを調べる。

等時降下を実現するのに必要な「仮想重力」は、残りの距離に単純に比例するものであるから、以下を得る。

$\frac{{\mathrm{d}}^{2}s}{{\mathrm{d}t}^2}=-k^{2}s$

$s=A\cos kt$

下式が上の微分方程式の解になっていることは容易に確かめられる。また、どんな高さ テンプレート:Mvar からも時間 テンプレート:Math で テンプレート:Math に達することも明らかである。問題は、このような「仮想重力」を産み出す曲線をつきとめることに帰着する。

$g\sin \theta =-k^{2}s$

残り距離を明示的に表現するのは困難だが、微分することにより問題を単純化できる。

$g\cos \theta \mathrm{d}\theta =-k^2\mathrm{d}s$

$\mathrm{d}s=-\frac{g}{k^2}\cos \theta \mathrm{d}\theta$

この方程式は曲線の傾斜角の変化と曲線に沿った距離の変化とを関係づけている。ここでピタゴラスの定理と曲線の傾斜は傾斜角の正接に等しいことおよび三角恒等式をいくつか用いると、テンプレート:Math を テンプレート:Math を用いて表現することができる。

$\begin{array}{cc}\mathrm{d}{s}^2& =\mathrm{d}{x}^2+\mathrm{d}{y}^2\\ & =(1+{\left(\frac{\mathrm{d}y}{\mathrm{d}x}\right)}^2)\mathrm{d}{x}^2\\ & =(1+{\tan }^2\theta )\mathrm{d}{x}^2\\ & ={\sec }^2\theta \mathrm{d}{x}^2\\ \mathrm{d}s& =\sec \theta \mathrm{d}x\end{array}$

これに最初の微分方程式を代入すると、テンプレート:Mvar を テンプレート:Mvar について解くことができる。

$\begin{array}{cc}\mathrm{d}s& =-\frac{g}{k^2}\cos \theta \mathrm{d}\theta \\ \sec \theta \mathrm{d}x& =-\frac{g}{k^2}\cos \theta \mathrm{d}\theta \\ \mathrm{d}x& =-\frac{g}{k^2}{\cos }^2\theta \mathrm{d}\theta \\ & =-\frac{g}{2k^2}(\cos 2\theta +1)\mathrm{d}\theta \\ x& =-\frac{g}{4k^2}(\sin 2\theta +2\theta )+C_x\end{array}$

同様に、テンプレート:Math を テンプレート:Math で表わし、テンプレート:Mvar を テンプレート:Mvar について解くことができる。

$\begin{array}{cc}\frac{\mathrm{d}y}{\mathrm{d}x}& =\tan \theta \\ \mathrm{d}x& =\cot \theta \mathrm{d}y\\ \cot \theta \mathrm{d}y& =-\frac{g}{k^2}{\cos }^2\theta \mathrm{d}\theta \\ \mathrm{d}y& =-\frac{g}{k^2}\sin \theta \cos \theta \mathrm{d}\theta \\ & =-\frac{g}{2k^2}\sin 2\theta \mathrm{d}\theta \\ y& =\frac{g}{4k^2}\cos 2\theta +C_y\end{array}$

テンプレート:Math および テンプレート:Math を代入すると、これら テンプレート:Mvar および テンプレート:Mvar についての方程式は水平線を転がる円上の一点の軌跡、すなわちサイクロイドであることがわかる。

$\begin{array}{cc}x& =r(\sin \varphi +\varphi )+C_x\\ y& =r(\cos \varphi )+C_y\end{array}$

テンプレート:Mvar を解き、降下時間が テンプレート:Math であることを用いると、降下時間を動円半径 テンプレート:Mvar で表わすことができる。

$\begin{array}{cc}r& =\frac{g}{4k^2}\\ k& =\frac{1}{2}\sqrt{\frac{g}{r}}\\ T& =\pi \sqrt{\frac{r}{g}}\end{array}$

（おおよそ テンプレート:Harvtxt に基く）

アーベルの解

ニールス・アーベルは一般化された等時曲線問題（アーベルの力学問題）、すなわちある関数 テンプレート:Math が与えられたとき、ある点からの総降下時間がこの関数となるような曲線を求める問題に取り組んだ。等時曲線問題は、アーベルの力学問題において テンプレート:Math が定数関数である特殊例に相当する。

アーベルの解は、エネルギー保存の原理を基礎とする。質点は摩擦を受けないので熱として失われるエネルギーはなく、運動エネルギーと位置エネルギーの和はつねに出発点における位置エネルギーに等しい。運動エネルギーは テンプレート:Math と表わされ、質点は曲線に沿って動くよう拘束されているので、その速さは単純に テンプレート:Math と書ける。ここで テンプレート:Mvar は曲線に沿って測った距離である。同様に、初期高さ テンプレート:Math から高さ テンプレート:Mvar まで落ちたときの重力による位置エネルギーの増加は テンプレート:Math であるから、以下を得る。

$\begin{array}{cc}\frac{1}{2}m{\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)}^2& =mg(y_0-y)\\ \frac{\mathrm{d}s}{\mathrm{d}t}& =\pm \sqrt{2g(y_0-y)}\\ \mathrm{d}t& =\pm \frac{\mathrm{d}s}{\sqrt{2g(y_0-y)}}\\ \mathrm{d}t& =-\frac{1}{\sqrt{2g(y_0-y)}}\frac{\mathrm{d}s}{\mathrm{d}y}\mathrm{d}y\end{array}$

最後の等式において、曲線に沿った残り距離は単調減少し（よって複号は負）、また高さの関数として表わせる (テンプレート:Math) ことを仮定してある。さらに、連鎖律 $\mathrm{d}s=\frac{\mathrm{d}s}{\mathrm{d}y}\mathrm{d}y$ を用いた。

ここで、テンプレート:Math から テンプレート:Math まで積分し、総降下時間を求めると以下を得る。

$T(y_0)=\underset{y=y_0}{\overset{y=0}{\int }}\mathrm{d}t=\frac{1}{\sqrt{2g}}\underset{0}{\overset{y_0}{\int }}\frac{1}{\sqrt{y_0-y}}\frac{\mathrm{d}s}{\mathrm{d}y}\mathrm{d}y$

これをアーベルの積分方程式と呼ぶ。これにより、曲線が与えられれば（かつ テンプレート:Math が計算できるならば）総降下時間を計算することができる。しかし、アーベルの力学問題では逆に、テンプレート:Math が与えられたときに テンプレート:Math を得たい。これもこの式から得ることができる。そのために、右辺を テンプレート:Math と テンプレート:Math の畳み込みとみなし、両辺をラプラス変換する。

$ℒ[T(y_0)]=\frac{1}{\sqrt{2g}}ℒ\left[\frac{1}{\sqrt{y}}\right]ℒ\left[\frac{\mathrm{d}s}{\mathrm{d}y}\right]$

$ℒ\left[\frac{1}{\sqrt{y}}\right]=\sqrt{\pi }{z}^{-\frac{1}{2}}$ であるから、テンプレート:Math のラプラス変換は テンプレート:Math のラプラス変換を用いて次のように書ける。

$ℒ\left[\frac{\mathrm{d}s}{\mathrm{d}y}\right]=\sqrt{\frac{2g}{\pi }}{z}^{\frac{1}{2}}ℒ[T(y_0)]$

一般の テンプレート:Math について言えることはここまでである。ひとたび テンプレート:Math が与えられれば、そのラプラス変換が計算でき、 テンプレート:Math のラプラス変換を得ることができ、逆変換を施すことにより テンプレート:Math が得られる。

等時曲線問題では、テンプレート:Math が定数関数である。テンプレート:Math のラプラス変換は テンプレート:Math であるから、次を得る。

$\begin{array}{cc}ℒ\left[\frac{\mathrm{d}s}{\mathrm{d}y}\right]& =\sqrt{\frac{2g}{\pi }}{z}^{\frac{1}{2}}ℒ[T_0]\\ & =\sqrt{\frac{2g}{\pi }}{T}_0{z}^{-\frac{1}{2}}\end{array}$

上述のラプラス変換を逆に利用することにより、次の結論を得る。

$\frac{\mathrm{d}s}{\mathrm{d}y}=T_0\frac{\sqrt{2g}}{\pi }\frac{1}{\sqrt{y}}$

サイクロイドがこの式に従うことは示すことができる。

関連項目

• 変分法
• テンプレート:仮リンク
• サイクロイド
• カテナリー曲線
• 等加速度運動
• 最速降下曲線

参考文献

• テンプレート:Cite book
• テンプレート:Cite book, posted by Cornell University Library.

外部リンク

• Mathworld