2

放物線上をはね続けるボールの研究③焦点を通らない軌道と包絡線の話

96
0

はじめに

 この記事は、先日私が公開した

【Mathlog】放物線上をはね続けるボールの研究①フィボナッチ数が生えた話

【Mathlog】放物線上をはね続けるボールの研究②はね続ける条件の話

の記事や、子葉さんが公開した

【Mathlog】放物線上を跳ねる物体とフィボナッチ数

と関連のある内容となっています。

 時間のある方はこの記事と合わせてお読みいただくとよりいっそう楽しめると思います。

 これまでの記事では、
「放物線上を、その放物線の焦点を通る軌道で発射した質点は永久にはね続ける」
という現象を取り扱っていましたが、この記事では主に
「放物線上を、その放物線の焦点を通らない軌道で発射した質点が描く軌道」
を取り扱います。

 まずは実際に、「下に凸な放物線の上を、その放物線の焦点を通らない軌道で発射した質点が描く軌道」を 2 パターンほどみてみましょう。

包絡線のパターンA 包絡線のパターンA

包絡線のパターンB 包絡線のパターンB

 質点がはね続けると、質点の通り得る領域を囲む包絡線が現れるのが見えますね。
 @aoki_taichi さんは現れる包絡線が、ベースとなる放物線と焦点を同じくする放物線になるのではないかと予想されましたが……

……結論から言うと、その予想が正しかったことを証明できました!

 実際に動かせるDesmosファイルを用意しましたので、よろしければ実際に動かして体感してみてください。

Desmosファイルへのリンク

 この記事ではこれらの包絡線について主に取り上げます。
 これまでの記事と同様、直感的に理解しやすいように、x 軸正方向を「右向き」、y 軸負方向を「下向き」などと表現しています。
 座標平面上に重力加速度が一様に下向きに大きさ 1 でかかっているものとします。衝突は完全弾性衝突(反発係数=1)とし、空気抵抗等は考慮しないものとします。また、反射のベースとなる放物線は原点に焦点を持ち、軸がy軸と重なるものとします。
 なお、これまでの記事とは違い、ベースとなる放物線は上に凸なものだけでなく、下に凸なものも含めて考えます。

包絡線のパターン

 まず、結果から書きます。

放物線上で反射する質点の描く放物線の焦点・準線・包絡線の関係

 質点の描く放物線の「焦点の原点からの距離」及び「準線の y 座標」は反射前後で変化しない不変量となる。それぞれ r1,Y1 とおく。

 原点焦点で (0,Y1±r12) を頂点とする放物線は、どちらも質点の軌跡と必ず接する。反射を繰り返すと包絡線として現れる。

 では実際の様子を、「①下に凸な放物線上、 Y1>r1 」「②下に凸な放物線上、 Y1<r1 」「③上に凸な放物線上、 Y1>r1 」「④上に凸な放物線上、 Y1<r1 」のパターン別に見ていきましょう。

①下に凸な放物線上、 Y1>r1 のとき

 
①下に凸な放物線上、 !FORMULA[13][-275474569][0] のとき ①下に凸な放物線上、 Y1>r1 のとき

 原点焦点で (0,Y1±r12) を頂点とする放物線はいずれも上に凸となり、間に挟まれた領域を通過する軌道となっています。

②下に凸な放物線上、 Y1<r1 のとき

②下に凸な放物線上、 !FORMULA[16][-277321611][0] のとき ②下に凸な放物線上、 Y1<r1 のとき

 原点焦点で (0,Y1+r12) を頂点とする放物線は上に凸、(0,Y1r12) を頂点とする放物線は下に凸となり、アリジゴクの巣のような形の領域を通過する軌道となっています。

③上に凸な放物線上、 Y1>r1 のとき

③上に凸な放物線上、 !FORMULA[20][-275474569][0] のとき ③上に凸な放物線上、 Y1>r1 のとき

 原点焦点で (0,Y1±r12) を頂点とする放物線はいずれも上に凸となり、間に挟まれた領域を通過する軌道となっています。

 なお、ベースとなる放物線が (0,Y1r12) を頂点とする放物線よりさらに下にある場合には、質点は一度も反射することなく奈落へ落ちていくことになります。

④上に凸な放物線上、 Y1<r1 のとき

④上に凸な放物線上、 !FORMULA[24][-277321611][0] のとき ④上に凸な放物線上、 Y1<r1 のとき

 原点焦点で (0,Y1+r12) を頂点とする放物線は上に凸、(0,Y1r12) を頂点とする放物線は下に凸となり、アリジゴクの巣のような形の領域を通過する軌道となっています。

 ベースとなる放物線の位置によっては、右向きに発射した質点が途中で左向きになったりします。

 なお、Y1=r1 のときは永久にはね続ける軌道になります。

包絡線

 それではここから、これらの性質が成り立つことを数式で確認していきたいと思います。

放物線族

 a0とします。原点に焦点を持ち、 (0,a) に頂点を持つ放物線の方程式は

 y=x24a+a

 と書くことができます。
 このとき、a>0 のときは上に凸な放物線となります。このような放物線をまとめて C+ と呼ぶことにします。
 逆に、a<0 のときは下に凸な放物線となります。このような放物線をまとめて C と呼ぶことにします。

 すると、C+C は必ず直交します。


証明はここをクリック / タップ

 y=x24a+a

a の方程式と見て解くと

 a=y±x2+y22

 復号のうち正の方を a1 、負の方を a2 とする。

 このとき、 a1>0,a2<0 となることに注意。

 また、それぞれ C+,C の放物線に対応することに注意。

 a1,a2 の方を固定して対応する座標平面上の点 (±x,y) における C+,C の接線の傾きを k1,k2 とする。

 k1=x2a1

 k2=x2a2

ここで

 y=x24a1+a1

 y=x24a2+a2

x2,y の連立方程式としてみて解くと

 x2=4a1a2

 y=a1+a2
 
 これを使うと


 k1k2=x24a1a2=1

 傾きの積が 1 となるから、C+C は直交する。


 また、xy 平面上の点のうち、x0 の点については、対応するC+,C の放物線がそれぞれ 1 つ定まります。

前提条件

 ここからは、平面全体に下向きに 1 の重力加速度がかかっている座標平面上で、座標 (x0,y0) から速度 v0=(vxvy) で発射した質点が、原点を焦点に持つ放物線上を弾性衝突し続けるときの軌道について考えます。
 場合分けを減らすため、x0<0,vx0 の場合のみ考えることにします。

 (x0,y0) に対応するC+,C の放物線をそれぞれ CP,CQ と呼ぶことにし、頂点をそれぞれ (0,P),(0,Q) します。(P>0,Q<0)

記号の説明 記号の説明

 x0,y0P,Q で表すとこうなります。

x0=2PQ

y0=P+Q


証明はここをクリック / タップ
 y0=x024P+P

 y0=x024Q+Q

を連立方程式とみて解く。まず y0 を消去する。

 x024P+P=x024Q+Q

 Qx02+4P2Q=Px02+4PQ2

 (PQ)x02=4PQ(PQ)

 x0=2PQx0<0,PQ

 元の式に代入して

 y0=P+Q
 

 次に、v0 を、CP(x0,y0) における接線右向きのベクトル VP と、CQ(x0,y0) における接線右向きのベクトル VQ に分解します。(図 7 参照)

 (x0,y0) における CP,CQ の傾きがそれぞれ QP,PQとなることから、適切に vp,vq をとることで次のように表すことができます。

VP=(vpPvpQ)

VQ=(vqQvqP)


証明はここをクリック / タップ
(x0,y0) における CP,CQ の接線の傾きを k1,k2 とする。

 k1=x02P=QP

 k2=x02Q=PQ
 
 となることから従う。


 ベクトルの成分を比較することで vx,vy,vp,vq の次の関係がわかります。

vx=vpP+vqQ

vy=vpQvqP

包絡線の検証

 原点を焦点に持つ放物線 y=x24A+A が放物線 y=y0+vy(xx0vx)12(xx0vx)2 と接するような A(0) を求めたい。y を消去して

 x24A+A=y0+vy(xx0vx)12(xx0vx)2

 両辺に 4Avx2 を乗じて

 x2vx2+4A2vx2=4Ay0vx2+4Avxvy(xx0)2A(xx0)2

 vx2x24A2vx2+4Ay0vx2+4Avxvy(xx0)2A(xx0)2=0

 (vx22A)x2+4A(vxvy+x0)x4A2vx2+4Ay0vx24Avxvyx02Ax02=0

 これを x2 次方程式と見れば、解は 2 つの放物線の交点の座標を表すのであるから、この方程式の判別式 D=0 となるような A を求めればよい。

 D4=4A2(vxvy+x0)2(vx22A)(4A2vx2+4Ay0vx24Avxvyx02Ax02)=2A(2A(vxvy+x0)2(vx22A)(2Avx2+2y0vx22vxvyx0x02))

 D=0,A0 より
 2A(vxvy+x0)2(vx22A)(2Avx2+2y0vx22vxvyx0x02)=0

 2Avx2vy2+4Ax0vxvy+2Ax02+2Avx42y0vx4+2vx3vyx0+vx2x024A2vx2+4Ay0vx24Ax0vxvy2Ax02=0

 vx2(4A22A(vy2+vx2+2y0)+2y0vx22vxvyx0x02)=0

 vx0 より

 4A22A(vy2+vx2+2y0)+2y0vx22vxvyx0x02=0

 A=vy2+vx2+2y0±(vy2+vx2+2y0)24(2y0vx22vxvyx0x02)4

 ここで 2y0vx22vxvyx0x020 のときは複号のどちらかで A=0 となり、不適となりますが、この初期条件のときは、yy0=vyxx0vx12(xx0vx)2(0,0) を通過する軌道、すなわち質点が原点を通る軌道となります。

 2y0vx22vxvyx0x020 でないときは複号のどちらでも A0 とはならず、質点は原点をとおりません。以下ではこのケースのみ考えることにします。

 
 つぎに、vx,vy,x0,y0vp,vq,P,Q で書き換える。(座標変換のアナロジー)

 vx2+vy2=(vpP+vqQ)2+(vpQvqP)2=vp2P+2vpvqPQvq2Qvp2Q2vpvqPQ+vq2P=(vp2+vq2)(PQ)

などを使って

 A=(vp2+vq2)(PQ)+2(P+Q)±((vp2+vq2)(PQ)+2(P+Q))24(2(P+Q)(vpP+vqQ)22(vpP+vqQ)(vpQvqP)(2PQ)+4PQ)4

 赤字の部分を抜き出して変形していく。

 2(P+Q)(vpP+vqQ)22(vpP+vqQ)(vpQvqP)(2PQ)+4PQ=2(P+Q)(vp2Pvq2Q+2vpvqPQ)+2(vp2PQvq2PQvpvqPvpvqQ)2PQ+4PQ=2P2vp22PQvq2+4PvpvqPQ+2PQvp22Q2vq2+4QvpvqPQ4PQvp2+4PQvq24PvpvqPQ4QvpvqPQ+4PQ=2P2vp22PQvq2+2PQvp22Q2vq24PQvp2+4PQvq2+4PQ=2P2vp22Q2vq22PQvp2+2PQvq2+4PQ=2(vp2P+vq2Q)(PQ)+4PQ

 元の式に戻すと

 A=(vp2+vq2)(PQ)+2(P+Q)±((vp2+vq2)(PQ)+2(P+Q))24(2(vp2P+vq2Q)(PQ)+4PQ)4=(vp2+vq2)(PQ)+2(P+Q)±(vp2+vq2)2(PQ)2+4(vp2+vq2)(P2Q2)+4(P+Q)216PQ8vp2P28vq2PQ+8vp2PQ+8vq2Q24=(vp2+vq2)(PQ)+2(P+Q)±(vp2+vq2)2(PQ)2+4(vp2+vq2)(PQ)2+4(PQ)2+8PQvp2+8PQvq28Q2vp28Q2vq28vp2P28vq2PQ+8vp2PQ+8vq2Q24=(vp2+vq2)(PQ)+2(P+Q)±(PQ)(vp2+vq2)2+4(vp2+vq2)+48vp24=(vp2+vq2)(PQ)+2(P+Q)±(PQ)(vp2+vq2)24(vp2vq2)+44

 A=(vp2+vq2)(PQ)+2(P+Q)±(PQ)(vp2+vq2)24(vp2vq2)+44

 ここで、vp,vq はいずれも 2 乗の形になっていることに注目してください。
 つまり、vpvp に置き換えたり、 vqvq に置き換えたりしても A2 解は不変、すなわち不変量ということになります。

符号を逆にする 符号を逆にする

 C+ 又は C で質点が反射した場合、vq 又は vp の符号が逆になることと合わせて、C+ 又は C で質点が反射しても、『原点を通り質点の軌跡に接する 2 つの放物線』は不変である」ことが言えたことになります。

焦点の原点からの距離や、準線の y 座標を求める

 質点の描く放物線の焦点を F0 と、準線の y 座標を Y0 として、F0,Y0x0,y0,vx,vy の式で表すことを考えます。

 まず頂点の座標を求めます。

 y=y0+vy(xx0vx)12(xx0vx)2=12vx2((x(x0+vxvy))2)+y0+vy22

と変形できるので頂点は (x0+vxvy,y0+vy22)

 x2 次の係数から焦点までの距離は vx22 とわかるので、

 Y0=y0+vx2+vy22

 F0(x0+vxvy,y0+vx2+vy22)

とわかります。

 つぎに、vx,vy,x0,y0vp,vq,P,Q で書き換えます。(座標変換のアナロジー)

Y0 について

 Y0=(P+Q)+(vp2+vq2)(PQ)2=(vp2+vq2)(PQ)+2(P+Q)2

 Y0=(vp2+vq2)(PQ)+2(P+Q)2

 vp,vq はいずれも 2 乗の形になっていることに注目してください。
 つまり、vpvp に置き換えたり、 vqvq に置き換えたりしても Y0 は不変、すなわち不変量ということになります。
 C+ 又は C で質点が反射した場合、vq 又は vp の符号が逆になることと合わせて、C+ 又は C で質点が反射しても、『質点が描く放物線の準線』は不変である」ことが言えたことになります。

 実はこのパート、変数を書き換えなくても不変量であることが証明できるのですが、記事の他の部分とのバランスを考えて変数を書き換えてから証明しています。

 また、よくみると Y0 は力学的エネルギーの和の式になっています。
 つまり、準線の位置は「運動エネルギーが全て位置エネルギーとなる高さ」と解釈することができます。

F0 について

 F0(2PQ+(vpP+vqQ)(vpQvqP),(P+Q)+(vpP+vqQ)2+(vpQvqP)22)

 x 座標について

 2PQ+(vpP+vqQ)(vpQvqP)=2PQ+vp2PQvq2PQvpvqPvpvqQ=(2vp2+vq2)PQvpvq(P+Q)

 y 座標について

 (P+Q)+(vpP+vqQ)2+(vpQvqP)22=(P+Q)+vp2P2vpvqPQ+vq2Qvp2Q2vpvqPQ+vq2P2=(2vp2+vq2)(P+Q)4vpvqPQ2

F0 の原点からの距離を r0 とすると

 r02=((2vp2+vq2)PQvpvq(P+Q))2+((2vp2+vq2)(P+Q)4vpvqPQ2)2=14(4(2vp2+vq2)2PQ+8(2vp2+vq2)PQvpvq(P+Q)+4vp2vq2(P+Q)2+(2vp2+vq2)2(P+Q)28(2vp2+vq2)(P+Q)vpvqPQ16vp2vq2PQ)=14(4(2vp2+vq2)2PQ+4vp2vq2(P+Q)2+(2vp2+vq2)2(P+Q)216vp2vq2PQ)=14(PQ)2(4vp2vq2+4+vp4+vq44vp2+4vq22vp2vq2)=14(PQ)2((vp2+vq2)24(vp2vq2)+4)

r0=(PQ)(vp2+vq2)24(vp2vq2)+42

r0=(PQ)(vp2+vq2)24(vp2vq2)+42

 vp,vq はいずれも 2 乗の形になっていることに注目してください。
 つまり、vpvp に置き換えたり、 vqvq に置き換えたりしても r0 は不変、すなわち不変量ということになります。
 C+ 又は C で質点が反射した場合、vq 又は vp の符号が逆になることと合わせて、C+ 又は C で質点が反射しても、『質点が描く放物線の焦点の原点からの距離』は不変である」ことが言えたことになります。

A,Y0,r0 の関係

 それでは、ここまでに得られた A,Y0,r0 の式を並べてみましょう。

 A=(vp2+vq2)(PQ)+2(P+Q)±(PQ)(vp2+vq2)24(vp2vq2)+44

 Y0=(vp2+vq2)(PQ)+2(P+Q)2

 r0=(PQ)(vp2+vq2)24(vp2vq2)+42

 明らかに次の関係があることがわかりますね!

 A=Y0±r02

 これで冒頭の性質を数式から確かめたことになります。

放物線上で反射する質点の描く放物線の焦点・準線・包絡線の関係(再掲)

 質点の描く放物線の「焦点の原点からの距離」及び「準線の y 座標」は反射前後で変化しない不変量となる。それぞれ r1,Y1 とおく。

 原点焦点で (0,Y1±r12) を頂点とする放物線は、どちらも質点の軌跡と必ず接する。反射を繰り返すと包絡線として現れる。

さいごに

 これで今回の記事はおしまいです。
 実は、先にご紹介したDesmosファイルですが、隠し機能がありまして

画像の名前 画像の名前

 「焦点・準線」フォルダの表示をオンにすることで、軌道の焦点・準線を確認することができますので、いろいろと試してみていただければと思います。

Desmosファイルへのリンク(再掲)

 当初、不思議な模様にみえていたものが、放物線の焦点と準線に注目してみると、背景にある実にシンプルで美しい関係が現れてきて、実に感動的でありました。
 また、この記事では数式で証明する流れでしたが、幾何的に証明することもできるのではないかと思います。
 例えば、衝突点を P 、衝突前後の軌道となる放物線の焦点をそれぞれ F1,F2 、ベースとなる放物線の焦点を Fb とおくと、「"PF1 の距離" と "PF2 の距離" は同じ」という性質がありますが、これは幾何的にも証明できそうです。また、この性質から、「F1F2 をつなぐ線分の垂直二等分線上に Fb があること」を導くことが可能ですし、そこから軌道の焦点が円を描くことも導くことができます。
 これらの現象を、力学的エネルギーとからめて物理的に解釈することも楽しそうです。
 もし何か面白い性質や新しい解釈を見つけたら教えていただければ幸いです!

投稿日:2024317
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

apu_yokai
apu_yokai
484
65353

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 包絡線のパターン
  3. ①下に凸な放物線上、 $Y_1>r_1$ のとき
  4. ②下に凸な放物線上、 $Y_1< r_1$ のとき
  5. ③上に凸な放物線上、 $Y_1>r_1$ のとき
  6. ④上に凸な放物線上、 $Y_1< r_1$ のとき
  7. 包絡線
  8. 放物線族
  9. 前提条件
  10. 包絡線の検証
  11. 焦点の原点からの距離や、準線の $y$ 座標を求める
  12. $A,Y_0,r_0$ の関係
  13. さいごに