■ダンデリンの2球問題(その10)
円柱を斜めに切ると切り口は楕円になる.円錐曲線の学び初めに,円錐を斜めに切ると切り口は楕円になることを学んだが,切り口は少しずつ先すぼまりになって卵形になるのではと思われた方は少なくないのではなかろうか?
切断に用いた平面と円錐の両方に接する頂点側と底面側の2つの球(ダンデリンの球)を考えると,球は楕円の2つの焦点と接することから,切り口は楕円になることを証明することができるのであるが,それを直観するのは難しい.そこで,・・・
====================================
円錐面をx^2+y^2=R^2zとしてもよいのであるが,母線の傾きを考えて
x^2+y^2=R^2z^2
とすると,母線はz=±y/Rとなる.
平面(z=ay+b)の傾きが母線よりと小さい場合(a<1/R),
x^2+y^2/(1+a^2)=R^2(ay/√(1+a^2)+b)^2
したがって,
x^2+(1−a^2R^2)/(1+a^2){y−abR√(1+a^2)/(1−a^2R^2)}^2=b^2R^2/(1+a^2)
x^2/α^2+(y−γ)^2/β^2=1
α=bR/√(1−a^2R^2),β=bR{√(1+a^2)}/(1−a^2R^2),γ=aRβ
γ=aRβ=aRbR√(1+a^2)/(1−a^2R^2)
となって,切り口は楕円である.卵形ではなかった.
α^2-β^2=(bR)^2/(1-a^2R^2)-(bR)^2(1+a^2)/(1-a^2R^2)^2=(bR)^2 /(1-a^2R^2)^2・{(1-a^2R^2)-(1+a^2)}
-α^2+β^2=(bR)^2/(1-a^2R^2)-(bR)^2(1+a^2)/(1-a^2R^2)^2=(abR)^2 (1+R^2)/(1-a^2R^2)^2
y軸上に焦点があることに注意
e=(abR)√(1+R^2)/(1-a^2R^2)
βe=bR{√(1+a^2)}/(1−a^2R^2)(abR)・√{(1+R^2)} /(1-a^2R^2)
βe=(ab^2R^2)・√(1+R^2) √(1+a^2) /(1-a^2R^2)^2
γ=aRβ=aRbR√(1+a^2)/(1−a^2R^2)
====================================
(x,y,z)→(X,Y,Z)の座標軸の回転を考えると,tanθ=a,cosθ=1/√(1+a^2),sinθ=a/√(1+a^2)であるから
x=X
y=1/√(1+a^2)・Y-a/√(1+a^2)・Z
z=a/√(1+a^2)・Y+1/√(1+a^2)・Z+b
平面(z=ay+b)に代入すると
a/√(1+a^2)・Y+1/√(1+a^2)・Z+b=a/√(1+a^2)・Y-a^2/√(1+a^2)・Z+b
これを解いてZ=0を得る
x=X
y=1/√(1+a^2)・Y-a/√(1+a^2)・Z
z=a/√(1+a^2)・Y+1/√(1+a^2)・Z+b
において、Z=0とした式
x=X
y=1/√(1+a^2)・Y
z=a/√(1+a^2)・Y+bを
x^2+y^2=R^2z^2
に代入すると
X^2+Y^2/(1+a^2)=R^2(aY/√(1+a^2)+b)^2
X^2/α^2+(Y−γ)^2/β^2=1
====================================
z=±y/R、y>0
z=ay+bを考える
円の方程式をy^2+(z-c)^2=r^2とすると,
小さいほうの円ではc/√(1+1/R^2)=r=(b-c)/√(1+a^2)→c(1/√(1+1/R^2)+1/√(1+a^2))=b/√(1+a^2))→c、rが求まる
c(R/√(1+R^2)+1/√(1+a^2))=b/√(1+a^2))
c={b/√(1+a^2)}/(R/√(1+R^2)+1/√(1+a^2))
c={b√(1+R^2)}/(R√(1+a^2)+√(1+R^2))
c={b√(1+R^2)}(R√(1+a^2)-√(1+R^2))/(R^2(1+a^2)-(1+R^2))
c={bR√(1+a^2)√(1+R^2)-b(1+R^2)}/(R^2a^2-1))
c-b={bR√(1+R^2)√(1+a^2)-b(1+R^2)}/(R^2a^2-1)-b
c-b=(bR√(1+R^2)√(1+a^2)-b(1+R^2+R^2a^2-1)}/(R^2a^2-1)
c-b={bR√(1+R^2)√(1+a^2)-bR^2(a^2+1)}/(R^2a^2-1)
c-b=bR{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)
r=c/√(1+1/R^2)=cR/√(1+R^2)=R{bR√(1+a^2)√(1+R^2)-b(1+R^2)}/(R^2a^2-1))(1+R^2)
接点を(y0,z0)とするとz=ay+b,z-c=-1/R・y→焦点
ay+b=-1/R・y+c
y=(c-b)/(a+1/R)=bR^2{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)(aR+1)
=(aR-1)bR^2{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)^2
z=ay+b・・・Z=0
x=X
y=1/√(1+a^2)・Y
z=a/√(1+a^2)・Y+bを代入すると
y=1/√(1+a^2)・Y
Y==(aR-1)bR^2{√(1+R^2)(1+a^2)-R(a^2+1)√(1+a^2)}/(R^2a^2-1)^2
Y=(aR-1)bR^2(1+a^2){√(1+R^2)-R√(1+a^2)}/(R^2a^2-1)^2
Y=-(aR-1)bR^3(1+a^2)√(1+a^2)}/(R^2a^2-1)^2+(aR-1)bR^2(1+a^2)√(1+R^2)/(R^2a^2-1)^2
Y=-(aR-1)bR^3(1+a^2)√(1+a^2)}/(R^2a^2-1)^2+(aR-1)bR^2(1+a^2)√(1+R^2)/(aR+1)(aR-1)(R^2a^2-1)
Y=-(aR-1)bR^3(1+a^2)√(1+a^2)}/(R^2a^2-1)^2+bR^2(1+a^2)√(1+R^2)/(aR+1)(R^2a^2-1)
これが-βe+γになっているようには見えない
βe=(ab^2R^2)・√(1+R^2) √(1+a^2) /(1-a^2R^2)^2
γ=aRβ=aRbR√(1+a^2)/(1−a^2R^2)
====================================
X=x
Y=1/√(1+a^2)・y+a/√(1+a^2)・(z-b)
Z=-a/√(1+a^2)・y+1/√(1+a^2)・(z-b)
に代入してみる
z=ay+bを代入するとZ=0
Y=(a^2y+1)/√(1+a^2)
Y√(1+a^2)=(a^2y+1)
y=(c-b)/(a+1/R)=bR^2{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)(aR+1)
y=(c-b)/(a+1/R)=bR^2{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)(aR+1)
=(aR-1)bR^2{√(1+R^2)√(1+a^2)-R(a^2+1)}/(R^2a^2-1)^2
(a^2y)={(aR-1)a^2bR^2√(1+R^2)√(1+a^2)-(aR-1)a^2bR^2R(a^2+1)}/(R^2a^2-1)^2
(a^2y+1)={(aR-1)a^2bR^2√(1+R^2)√(1+a^2)-(aR-1)a^2bR^2R(a^2+1)+(R^2a^2-1)^2}/(R^2a^2-1)^2
簡単にはなりそうにない
====================================
確認のため、z=ay+bと母線z=-y/Rの交点を求めてみたい。
y=b/(a+1/R)=bR/(aR+1)
z=-bR/(aR+1)
y=1/√(1+a^2)・Y
z=a/√(1+a^2)・Y+b
に代入するとY=b√(1+a^2)/(a+1/R)=bR√(1+a^2)/ (aR+1)≠β=bR{√(1+a^2)}/(1−a^2R^2)
Y=1/√(1+a^2)・y+a/√(1+a^2)・(z-b)
Z=-a/√(1+a^2)・y+1/√(1+a^2)・(z-b)
に代入すると
Y=1/√(1+a^2)・y+a/√(1+a^2)・(z-b)=1/√(1+a^2)・{bR/(aR+1)+a(-Rb/(aR+1)-b)}
=1/√(1+a^2)・{bR/(aR+1)+a(-Rb/(aR+1)-b)}=1/√(1+a^2)・(bR-abR-b(aR+1)}/(aR+1)≠β
そもそもの計算が誤りだったことになる
====================================
接点を(y0,z0)とするとz=-y/R,z-c=a・y→準線
y=1/√(1+a^2)・Y-a/√(1+a^2)・Z
z=a/√(1+a^2)・Y+1/√(1+a^2)・Z+b
Y=1/√(1+a^2)・y+a/√(1+a^2)・z
Z-b=-a/√(1+a^2)・y+1/√(1+a^2)・z
大きいほうの円ではc/√(1+1/R^2)=r=(c-b)/√(1+a^2)→c(1/√(1+1/R^2)-1/√(1+a^2))=-b/√(1+a^2))
接点を(y0,z0)とするとz=ay+b,z-c=-1/R・y→焦点
接点を(y0,z0)とするとz=y/R,z-c=a・y→準線
====================================
x^2/a^2+y^2/b^2=1、e=√(a^2-b^2)
焦点(±ae,0)
準線(x=±a/e)
X^2/α^2+(Y−γ)^2/β^2=1、e=√(α^2-β^2)
焦点(±ae,0)+γ→aRβ-a√(α^2-β^2)
準線(x=±a/e)+γ
α=bR/√(1−a^2R^2),β=bR√(1+a^2)/(1−a^2R^2),γ=aRβ
====================================