前回の記事「「2つのボールをぶつけると円周率がわかる」らしいのでシミュレーションしてみた」の続き。 この記事では、衝突回数が円周率(の適当に桁をズラした値)になることを解析的に導いてみます。
- 大人になってからの再学習 『2つのボールをぶつけると円周率がわかる』
目次
問題設定
質量がそれぞれ
回衝突後の速度 
回数をカウントする衝突には2種類あります:
: 質点同士の衝突
: 質点2と壁との衝突
前回見たように、この2種類の衝突は交互に起こります。 奇数回目の衝突は 、偶数回目の衝突は
です。
回衝突した後の
の速度をそれぞれ
とおきましょう(速度は2つの質点から壁の方向を正):

これらを成分とした2成分のベクトルを考えます(要は速度空間を考える):
の漸化式と初期条件
前回の結果から、2つの衝突
行列 は前回計算していて
で与えられます。 ちなみに初期条件は
となります。
の一般項(行列を用いた表式)
行列
で定義すると、上記の初期条件と漸化式より
となります。 は
に
を作用させても得られますが、上記の方が簡単。 よって、次に
を求めましょう。
を求める
この節では
の書き換え
質点同士の衝突での速度の変換を表す行列
のように常に になることを踏まえると、
を満たす と
を用いて、
は
と書けます( のときは
)。 もともと
は
だけで書けてたので、この
と
は独立ではありませんが、それは後ほど。
をこの形に書いたとき、行列
は
のように回転行列みたく書くことができます。
を求める
上記の
となります。 証明は数学的帰納法で。 のときは成り立つのは明らか(
も単位行列になるので成り立ちますが)。
のときに成り立つと仮定すると
よって のときも成り立ちます。
を求める
以上をまとめると、
これで 回衝突した後の速度が求まりました。 次は、
(これは
によって
と関係していて、
は
を使って定義されていた)と衝突回数
の関係を導きます。
と衝突回数
との関係
以下の手順によって、この関係を求めましょう:
- 衝突が起こらなくなる条件を求める
- ベクトル
を導入
- この条件から衝突回数
と角度
(後述)の間の関係を求める
- 衝突回数
と
との関係を求める
衝突が起こらなくなる条件
この条件は前回見ました。
衝突が だけ起こるとき、
が初めて上記の条件を満たせばいいことが分かります。
ベクトル 
ベクトル
で定義しましょう。 このとき
となり、これらは -空間で円周
上にあることが分かります。 また、衝突が起こらなくなる条件は
となります。 これらを図に描くと

となります。 角度
で関係しています()。
について、
が奇数のものを見ていくと、
から円周の上半分を反時計回りに
ずつ回転していきます。 また、
が偶数のものは、その1つ前の(
が奇数の)
の真下に位置します。 つまり、円周の下半分を時計回りに
ずつ回転していきます。 衝突が起こらなくなる条件は、この円周上の点が図で示した青色の領域(境界を含む)に入ることです。
補足
前節で
と関係しているので は互いに独立ではないと書きましたが、
の代わりに
と
との関係を見ておきましょう。 正接の倍角の公式(こちらを参照)より
となって、 であることがわかります。 したがって、行列
は以下のように
いずれかだけで書けます:
まぁ、この表式を使うわけではありませんが。
衝突回数
と
の関係を求める
まずは
偏角が に等しくなる場合、2回前の
回衝突したときの偏角が
となり、その次の
回衝突したときの偏角が
となるので、この時点で衝突が終わります。
次は 、つまり
が偶数の場合。
が円周の下半分にあることに注意して、これの偏角は
です。 このとき衝突が起こらなくなる条件は
偏角が に等しくなる場合、その1回前の
の衝突の時点で偏角が
になるので、それ以上衝突は起こりません。
以上をまとめると、衝突回数が偶数、奇数どちらの場合もまとめて書けて
もしくはガウス記号 [ ] を用いて
となります。
衝突回数
と
との関係を求める
が成り立つので、結局衝突回数 は
で与えられます。 前回行ったシミュレーションでは、 が充分大きいという条件は必要なく、0以上の自然数で成り立っていましたが、評価が面倒なのでこの辺で終了。
まとめ
思ったより長い導出になりましたが(最後の近似の評価は投げたけど)、一応高校の数学と物理の知識で大体出せましたね。 と言いつつ、行列は高校の数学からなくなったらしいですが。 まぁこういう話は、導出はある程度理解できても、そもそも「衝突回数が円周率を与える」ということに気付くなんてことがそう簡単にできねぇよ、という結論になってしまうんですが。 原論文とか読めばそのあたりのことも書かれてるのかもしれませんが、あんまり読む気しない・・・ ビリヤードっぽい基本的な力学系なら、量子ビットか何かの関連で量子系への拡張することを念頭においてまずは古典系でやってみた、みたいなもんじゃないかと推測してみたり。ちなみに
この記事は海外では1ヶ月前くらいに雑誌に掲載されたとか書いてましたが、3月といえば3月14日の円周率の日にネタとして載ってたのかな?
- 作者: Alfred S. Posamentier,Ingmar Lehmann,松浦俊輔
- 出版社/メーカー: 日経BP社
- 発売日: 2005/11/03
- メディア: 単行本
- 購入: 2人 クリック: 32回
- この商品を含むブログ (19件) を見る
*1:挟まれた数を超えない最大の整数を与える。