キャッチボール

「投げられたボールが 秒後に、どの位置 にあるか」を計算したい。このサイトでは、関数を表す際、 ではなく と表すことが多い（式()などを見やすくするため）。 の形式も場合によっては使う。

ボールの運動 は、初期値()と運動方程式から計算できる

キャッチボールをしている場面を想像してみよう。投げられたボールの運動は、投げるたびに異なったものとなるだろう。しかし経験的に分かるように、「同じ投げ方」をすればボールの運動は毎回同じになる（だからこそボールの軌道を予測して捕球できる）。これは非常に興味深いことである。というのも、初期値（＝投げ方）からボールの運動 （＝任意の時刻 におけるボールの位置）を一意的に決定する物理法則が存在することを意味しているからである。この物理法則を、（キャッチボールの）運動方程式という。

運動方程式を見つけて を計算したいのだが、その前にまず、「同じ投げ方」とは何なのかを具体的にしておく必要がある。これは経験的に分かるだろう。即ち、ボールの運動を一意的に決めるためには、以下の3つを固定すればよい： ボールをどの位置から投げるかどの方向に投げるかどれくらいの速さで投げるか とはいえ、これには不満がある。「どの位置から」とか「どの方向に」といった表現は感覚的すぎるのである。今後、何らかの計算を行うためには、これらを や といった数字で表す（＝定量化する）必要がある。細かいことを言えば、風やボールの回転による影響があるため、初期値()だけではボールの運動が確定しない。しかしここでは、それらの影響は無視できるとしている（例えば十分に重いボールを考える）。

以上をまとめると、ボールの運動 を計算するためには、3つのステップ：「初期値()の定量化」、「運動方程式の導出」そして「ボールの運動の計算」の順に考えるのが良さそうである。そこでこの章では、これらをそれぞれの節に分けて議論する： 初期値の定量化運動方程式の導出ボールの運動の計算

1.1初期値()の定量化

この節では、初期値()の定量化を行う。「どこから投げるか」は位置ベクトル （式()）で表され、残りの2つ「方向と速さ」は速度ベクトル （式()）で表される。

「どこから投げるか」は、デカルト座標で表す：式()

まず、初期値()の1つ目「どこから投げるか」である。これは単に、 でのボールの位置 である。従って、ボールの位置を定量化する方法を与えればよい。

これは特に問題ないだろう。即ち、慣れ親しんだデカルト座標を使えばよい。具体的には右図のように、ボールの位置 を、「原点 から 軸方向に 、 軸方向に 、 軸方向に 」という三つ組のベクトル： によって表すわけである。（あらかじめ「原点 」および「互いに直交する 軸」を適当に取っておく。）この のように、太字の変数はベクトルを表すものとする。

結局、「ボールをどこから投げるか」は、デカルト座標により と表すことができる。この章の冒頭でも、投げてから 秒後のボールの位置を と書いていたが、これについてもデカルト座標で表されているとする。もちろん、デカルト座標以外の座標系を使ってもよい。例えば、振り子の場合、振れ角によって重りの位置を表すしてもよい。そのような例は第9章で扱う。

「投げる方向と速さ」は、速度ベクトル()で表される

次は、初期値()の残りの2つ：「投げる方向」と「投げる速さ」を定量化する方法である。投げる「方向」については、その方向を向くベクトルで表現できそうだが、「速さ」については自明ではない。注目するのは、ボールの運動方向や速さを知るには、ごく短い間だけでも、ボールの運動を観察しなければならないという事実である。例えば、ボールの「運動方向や速さ」を、撮影された写真から判断しようとする時、1枚の写真だけでは駄目で、連続する複数の写真が必要となることからも分かるだろう。連続する写真（右図）における位置の「変化の方向」から運動の方向が、「とびの大きさ」から速さが分かる。 そこで、右図のように、微小時間 秒の測定で得られたボール位置 をもとに、「ボールの方向と速さ」を定量的に得る方法を探ろう。

まず、 は非常に小さいので、この間の運動は直線的（＝ が の1次関数）であるとみなせる。即ち、 の係数を とおいて、 は以下の形で近似できる： （詳しくは以下の【1.1-注1】の式()参照。）左辺 は、右辺にある2つの量 と から決まるわけだが、 は測定開始時刻 でのボールの位置なので、 を測定することで新たに得られる情報は、もう一方の ということになる。よって、「運動の方向と速さ」は、 に含まれているはずである。

そこで、 がどのような量なのか、詳しく見てみよう。まず、上式()の赤字部分は、 秒間にボールが からどれだけ移動したかを表すベクトルである（右図の赤矢印部分）。よって、 の向きは、運動方向に等しい。また、赤字部分の大きさ は、「 秒間での移動距離」なので、公式「移動距離=速さ・時間」と見比べれば、ベクトルの大きさ は、時刻 でのボールの速さに等しいことが分かる。つまり、 は、「ボールの運動方向を向き、大きさが運動の速さに一致する」ベクトルであり、まさに「運動の方向と速さ」を過不足なく持っている。この を速度ベクトル（あるいは単に速度）という。

以上により、「投げる方向と速さ」は、 における速度ベクトル： で表されることになる。非常にきれいな結果が得られたので、運動方程式についても簡潔に書けるのではないかという期待が持てる。なお、実際に を測定するには、 近傍での微小時間の観測を行って を特定し、式()を用いればよい。

1.2運動方程式の導出

前節では、初期位置 と初期速度 を決めれば、その後のボールの運動 が一意的に決まることを見た。この状況を模式的に表すと以下のようになる： 「時間を秒進める法則」 この節では、「時間を 秒進める法則」、即ち、運動方程式()を書き下す。

ボールの運動は、加速度ベクトル()によって決まる

式()に肝である「時間を 秒進める法則」を求めたい。そのためにまず、初期時刻を として、非常に小さな値 だけ時間を進めることを考える： 「時間を秒進める法則」 両辺を対称にするため、右辺にも速度 を加えた。これは、「 の関数形が決まっている→その時間微分 も決まる」ことから正当化できる。これを漸化式の要領で連続して使用すれば、 から始めて ずつ、任意の時刻まで進めることができる。よって、「時間を 秒進める法則」が分かれば、式()についても理解したことになる。

さて、なぜ微小時間をｓ考えているかというと、前節で有効だった1次近似がここでも上手くいくことを期待しているからである。実際、式()の右辺を1次近似すると、式()の要領で以下のようになる： この式の右辺の未知数は のみである（それ以外の は式()の左辺として与えられている）。 は、速度ベクトル の時間微分、即ち、速度の変化率＝右上図のように、微小時間間隔 を掛けると、 秒間での速度の変化量になるような量である。 を、加速度ベクトル（または単に加速度）と言う。

よって結局、ボールの運動を決めるのに必要な量は、任意の時刻 における加速度ベクトル ということになる。そうすれば、漸化式()によって1ステップずつ、任意の時刻まで時間を進めていくことができる（「積分する」という）。

なお、式()は近似式なので、等式と見なして時間を進めていくと誤差が生じる。よって、刻み幅 を小さくしていった時に、任意の時刻 での誤差 がゼロになることを言う必要がある。直感的に言うと、 は ステップ数ステップあたりの誤差・より、・以下の【注】より、 程度の量となる。よって、 の極限で となる。ただしこの説明では、1ステップでの誤差が次のステップで増幅されていく可能性が考慮されていないので不十分である。厳密に議論しても同じ結論になるのだが、それは第3章で行う。

キャッチボールの運動方程式：式()

後は、加速度 を決定するだけである。単純な例として、重力のない宇宙空間では である。なぜなら、ボールは一定の速度で進み続けるため、式()の第2式が となるからである。 がどのような法則に従うかを見出すには、実際に様々な初期値・時刻のもとで を測定してみればよい。測定によって の値を求めるには、短い時間間隔で位置 をで3回測定して、以下の【1.2-注2】の式()を使えばよい。

実際にそのような測定が行われ、その結果、 は、重力加速度と呼ばれるある一定の値 となることが知られている（以下の【1.2-注3】）： （ボールは下向きに加速するのだから、 は下方向を向くベクトル。）これが求めたかった運動方程式である。

以上により、初期値 が与えられれば、運動方程式()によって、任意の時刻 でのボールの位置 が決まることになる。ボールの素材は無数に考えられるが、それら運動は全てこの1つの単純な運動方程式に支配されているのである（空気抵抗は無視している）。ボールの運動の違いは、初期値が異なることに起因する。

1.3運動方程式()の解

ここまでの議論により、初期位置 ・初期速度 が与えられれば、運動方程式()：（再掲） を解くことで、任意の時刻 での位置 が分かる。

この節では、式()の解として、数値解と解析解について述べる： 「数値解」＝コンピュータを用いて計算する近似的な解「解析解」＝数学的に厳密な解 「数値解」は、様々な場面で一般的に使える手法だが、近似値なので、誤差が許容範囲内か考慮する必要がある。（数値界に限ったことではないが、誤差はこれ以外にも、実験における「初期値の測定誤差」や空気抵抗を無視するといった「単純化による誤差」も考慮すべきである。）「解析解」は、誤差を含まない厳密な解だが、問題に応じて解法が異なり、そもそも計算できないことのほうが多い。

数値解

数値的に運動方程式()を解くには、上述の1次近似式()を使って逐次的に計算を進めていけばよい。即ち、以下の漸化式：（ を1つのベクトルにまとめた） を用いて、 から始めて ずつ進めていけばよい。 を小さくしていけば精度が良くなる。（前節でも述べたが、数学的な正当化は第3章で行う。）このように、「いきなり大きな時間幅を考えるのは難しいので、微小な時間幅 を考えることで簡単な関係式を探す」というのが微分の考えであり、「そうして得られた関係式()を繋いでいって大きな時間幅での結果を得る」というのが積分の考え方である。このように、微分と積分はペアとなっており、今後も時間発展を計算する際には、微分形の方程式→積分というロジックを考えることになる。

実際に数値計算を行うと、右図のようになる。グリッドの間隔は である。

解析解：式()

運動方程式()の解 は、解析的に示すこともできる。実際、初期値 を用いて となる。この式は確かに、 における初期条件「 および 」が成立し、任意の において運動方程式 を満たしている。ボールの運動は時刻 の2次関数であることが、式()から分かる。これは上述の数値解法では分からなかった性質である。このように、解析解のほうが得られる情報は多い。ただし、一部の単純な場合を除いて、解析解の導出は困難である。逆に、数値解法は汎用性がある。【式()の導出の概要】まず、運動方程式 より、 を2階微分すると定数 になるのだから、 は の2次関数になるはずである：（ は定数） 後は、初期値条件 および と、運動方程式 を満たすように を求めると、一意的に解()が得られる。