最速降下曲線 †変分法の夜明けとなった重要な曲線である。この解法がレオンハルト・オイラーによって改良され後に変分法と呼ばれるものになった。 変分法は、 A Brachistochrone curve (brachistos - the shortest, chronos - time), or curve of fastest descent, is the curve between two points that is covered in the least time by a body that starts at the first point with zero speed and is constrained to move along the curve to the second point, under the action of constant gravity and assuming no friction. 最速降下曲線はサイクロイド †最速降下曲線というのは,曲線に沿って物を転がしたときに,物が一番速く転がり落ちる曲線 ベルヌーイの提起した問題とは 『決まった二点の間を,始点から終点まで玉が一番速く転がることが出来るような曲線を求める』 というものです. Given two points A and B, with A not lower than B, there is just one upside down cycloid that passes through A with infinite slope, passes also through B and does not have maximum points between A and B. This particular inverted cycloid is a brachistochrone curve. The curve does not depend on the body's mass or on the strength of the gravitational constant. The problem can be solved with the tools from the calculus of variations. Note that if the body is given an initial velocity at A, or if friction is taken into account, the curve that minimizes time will differ from the one described above. 質点の運動 †質点が原点A から点B= (x0; y0) まで曲線(x; y(x)) に沿って降下すると,速さv は で与えられる. 質量をm とすると位置エネルギーと運動エネルギーの和は一定なので,原点から初速0 で降下することにより, 1/2mv2 + mgy = 0 となり v =2gy^1/2 である。 原点から点Bまでの所要時間は となる。 最速降下線の問題は次の汎関数を最少にする変分法の問題である。 解法 †汎関数を使って最適化する。変分法の知識が必要です。
汎関数、最少作用の原理とオイラー・ラクランジェの方程式 †
ことです。 最小作用の原理: 力学系の運動は関数:F[t,y(t),y'(t)]によって特徴づけられ,実際に実現する運動の軌跡は 作用積分,S[y(t)]= ∫a->b F[t,y(t),y'(t)]dt を最少にするy*(t) である。 このような関数Fは,ラグランジェアンとよび Lで表わすことが多い。
φ'(ε)=∫ (∂F/∂y・∂y/∂ε+∂F/∂y'・∂y'/∂ε)dt y(t,ε)=y*(t)+εη(t), y'(t,ε)=y*'(t)+εη'(t)より φ'(ε)=∫ (∂F/∂y・η(t) +∂F/∂y'・εη'(t))dt =∫ (Fy・η +F/y'・η')dt この条件式を第一変分といいます 汎関数に極値を与える y*(t) は方程式 ∫ (Fy・η +F/y'・η')dt = 0 を満たすことがわかる。 さらに第2項の部分積分(t の関数として)を実行すると ∫a->b Fy'η'dt=[Fy'η] + ∫ dFy'/dt・ηdt 上の第1項はη(a)=η(b)=0 なので0 です。 結局,φ'(0)=0 の条件は ∫{Fy- dFy'/dt]・ηdt となる。上の式で、η(t)が任意に変化しても成立するためには Fy- dFy'/dt =0 ・・・・・・・オイラー方程式 解析力学でなじみのあるラグランジェアンLを用いて書けば ∂L/∂q-d(∂L/∂q')/dt=0 ・・・・・・ラグランジェの方程式 を得ることができた。 y(t)=y*(t)が汎関数を最少化するための必要条件は、y*(t)がオイラーの方程式を満たすことである オイラー方程式から解を求める †F=f(t,y,y')=(1+y'^2)^1/2・(2gy)^-1/2 を使って∂F/∂y'と∂F/∂yを計算する。 ∂F/∂y'= y'・(2gy)^-1/2・(1+y'^2)^-1/2 ∂F/∂y= -1/2・(1+y'^2)^1/2・(2g)^-1/2・y^-3/2 オイラーの方程式 Fy- dFy'/dt =0 を整理して、微分方程式を得る。 y(1 + y'^2) = C1 で y = k1+k2cosθ とおいてこれに代入して解く y = C1/2・(1 − cos θ) x = C1/2・(θ − sin θ) となり、サイクロイド曲線である。 変分法の例題:垂らした紐の形 †線密度ρ、長さl のひもの両端を固定して垂らしたときのひもの形を求めます。 ひもは、位置エネルギーが最小になるような形をとります。ひもの微小部分の質量は ρ(1 + y'^2)^1/2 なので 制約:∫(1 + y'^2)^1/2 dx=Length のもとで ∫gρy(1 + y'^2)^1/2 dx の極値を求めることになります。ただしgは重力係数 こたえは、懸垂線 懸垂線y +λ = C1 cosh[(x − C2)/C1 ] 歴史的記述 †
|