■ペリトロコイド曲線(その24)
黄色の線で示したペリトロコイド
x=(R−r)cosθ+Rcos((R−r)/Rθ)
y=(R−r)sinθ−Rsin((R−r)/Rθ)
では、長さ2Rの線分を1回転させることができる。(尖点が奇数の場合)
そこで、掛谷の定数のような面積の極限値を求めてみたい。
===================================
【1】2n-1尖点ペリトロコイドの面積
2n-1個の尖点をもつペリトロコイドは,パラメータθを用いて
x=(n−1)rcosθ+nrcos(1−1/n)θ
y=(n−1)rsinθ−nrsin(1−1/n)θ
θ(0、2nπ)
と記述されます.
θで微分すると
x’=−(n−1)rsinθ−(n−1)rsin(1−1/n)θ
y’=(n−1)rcosθ−(n−1)rcos(1−1/n)θ
ここで注意しなければならないことは,θは極座標(r,θ)のパラメータではないことです.そのため,
S=1/2∫r^2dθ r^2=x^2+y^2
として計算すると正しい値が得られません.
計算方法はいくつか考えられるのですが,
S=∫ydx=∫yx’dθ
S=∫xdy=∫xy’dθ
S=1/2∫(ydx-xdy)=1/2∫(yx’-xy’)dθ
その結果,ペリトロコイドの面積は
S=n(n−1)・πr^2
で表されることが計算されます.回転円の半径をR(=nr)とした場合は,
S=n(n−1)/n^2・πR^2
となります.
ここで求めたい値はR=1/2のときの面積ですから
S=(n−1)/4n・π
デルトイドではR=2rですから
S=π/8
以下、R=3rではS=π/6
R=4rではS=3π/16・・・これは半アステロイドの面積に等しい
R=5rではS=π/5
また,n→∞のとき
S→π/4
となります。
===================================
n=3,7のとき、ペリトロコイドは大円の直径を含みますが、n=5,9のときはそれよりも大きい線分を含むことができます。
y=0となるθを求めると
y=(n−1)rsinθ−nrsin(1−1/n)θ
(n−1)sinθ=nsin(1−1/n)θ
sinθ=n/(n−1)sin(n−1)/n)θ
sinθ=αsinθ/αを解く必要があるが、αは整数とならない。
そこで、t=θ/nと起きます。T=[0,2π]
x=(n−1)cosnt+ncos(n-1)t
y=(n−1)sinnt-nsin(n-1)t
n=2のとき,y=sin2t-2sint,x=cos2t+2cost
n=3のとき,y=2sin3t-3sin2t,x=2cos3t+3cos2t
n=4のとき,y=3sin4t-4sin3t,x=3cos4t+4cos3t
n=5のとき,y=4sin5t-5sin4t,x=4cos5t+5cos4t
===================================
間違いを訂正T=[0,2π]
[n=2]
2sintcost=2sint
cost=1
x=2(cost)^2-1+2cost
x=3・・・OKだがx=-1はどこへ行ったのだろうか?
===================================
[n=3]
6sintcost=2(-4(sint)^3+3sint)
3cost=-4(sint)^2+3
3cost=-4(1-(cost)^2)+3=4(cost)^2- 1
cost=1,-1/4
x=2cos3t+3cos2t=8(cost)^3-6(cost)+6(cost)^2-3
cost=1のときx=5・・・OK
cost=-1/4のときx=-1/8+3/2+3/8-3=-5/4・・・OK
x=1はどこへ行ったのだろうか?
===================================
[n=4]
-24(sint)^3cost+12sintcost=-16(sint)^3+12sint
-24(sint)^2cost+12cost=-16(sint)^2+12
-6(sint)^2cost+3cost=-4(sint)^2+3
-6{1-(cost)^2}cost+3cost=-4{1-(cost)^2}+3
-6cost+6(cost)^3-3+3cost=-1+4(cost)^2
6(cost)^3-3cost=-1+4(cost)^2
6(cost)^3-4(cost)^2-3cost+1=0
6z^3-4z^2-3z+1=0
(z-1)(6z^2+2z-1)=0
z=1,(-1+/-√7)/6
x=3cos4t+4cos3t=24(cost)^4-24(cost)^2+3+16(cost)^3-12cost
cost=1のときx=7・・・OK
cost=(-1+/-√7)/6のとき
(-1+√7)^2=8-2√7=2(4-√7)
(-1+√7)^3=-1+3√7-21+7√7=-22+10√7
(-1+√7)^4=4(23-8√7)
x=24・4(23-8√7)/6^4+16・(-22+10√7)/6^3-24・2(4-√7)/6^2-12(-1+√7)/6+3
x=2(23-8√7)/27+2・(-22+10√7)/27-36(4-√7)/27-2(-1+√7)+3
x=(-142+40√7)/27-2(-1+√7)+3=-0.26+1.37
(-1-√7)^2=8+2√7=2(4+√7)
(-1-√7)^3=-1-3√7-21-7√7=-22-10√7
(-1-√7)^4=4(23+8√7)
x=24・4(23+8√7)/6^4+16・(-22-10√7)/6^3-24・2(4+√7)/6^2-12(-1-√7)/6+3
x=2(23-8√7)/27+2・(-22+10√7)/27-36(4-√7)/27-2(-1+√7)+3
x=(-142-40√7)/27-2(-1-√7)+3=-0.26-1.37
(1,-1に近い解)が存在するが
x=-1はどこへ行ったのだろうか?
===================================
[n=4]
64(sint)^5-80(sint)^3+20sint=-40(sint)^3cost+20sintcost
16(sint)^5-20(sint)^3+5sint=-10(sint)^3cost+5sintcost
16(sint)^4-20(sint)^2+5=-10(sint)^2cost+5cost
16{1-(cost)^2}^2-20{1-(cost)^2}+5=-10{1-(cost)^2}cost+5cost
16-32(cost)^2+16(cost)^4-20+20(cost)^2+5=10cost-10(cost)^3+5cost
16(cost)^4-12(cost)^2+1=-10(cost)^3+15cost
16z^4+10z^3-12z^2-15z+1=0
(z-1)(16z^3+26z^2+14z-1)=0
z=1と(-1に近い解)が存在する
===================================