（その４１）をやり直し

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

θ=π/n、m=(n-1)/2として、掛谷定数の計算をリパラメトライズしておきたい。

ｎ芒星では三角おむすび状の図形がｎ個出現する。

原点からお結びの頂点までの距離をrとすると、おむすびの頂点は(x,y)=(rcosθ,rsinθ)

また、弧の半径をRとする。

(1,0),(cos2mθ,sin2mθ)で接線をひくとR=tan(mθ)=cot(θ/2)として,その中心は(1,R)となることがわかる

Rsinθ=2{cos(θ/2)}^2

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

(1,R)を焦点とする楕円・放物線・双曲線(x^2/a^2+y^2/b^2=1,y^2=4px,x^2/a^2-y^2/b^2=1)を求めたい

軸はy=tan(nθ)x

(1,0),(cos2mθ,sin2mθ)を通る。

(x,y)を標準的な座標(x^2/a^2+y^2/b^2=1,y^2=4px,x^2/a^2-y^2/b^2=1)とすると

X-1=[cosmθ,-sinmθ][x]

Y-R=[sinmθ, cosmθ][y]

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

極座標r=L/(1-ecos(t)),(r,t)のほうが使いやすいかもしれない

(1,0)を通るからt=-π/2-mθのとき R=L/(1-ecos(t)),L=R(1-ecos(t))

(cos2mθ,sin2mθ)を通るからt=π/2+mθのとき R=L/(1-ecos(t)),L=R(1-ecos(t))

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

角度の測り方を変えて

r=L/(1+ecos(t)),(r,t)のほうが使いやすい

t=[-(π/2-mθ),(π/2-mθ)]

S=1/2∫r^2dt

e=1のときしか計算は難しそうである。

r=L/(1+cost)=L/2{cos(t/2)}^2

r^2=L^2/4{cos(t/2)}^4・・・この不定積分は既知である

これで面積は測れるようになっらが、あとは交点が求められるかである・・・

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

(x-1)^2+(y-R)^2=R^2とy=tanθ・xの交点が(rcosθ,rsinθ)であるから

ｒはr^2cosθ^2-2rcosθ+1+r^2sinθ^2-2Rrsinθ+R^2=R^2の解

r^2-2r(cosθ+Rsinθ)+1=0

r=cosθ+Rsinθ-{(cosθ+Rsinθ)^2-1}^1/2

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

S1=1/2・(tanθ-θ)

S2=1/2・(1-rcosθ)・(R-tanθ)

S3=1/2・R^2arctan{(1-rcosθ)/(R-rsinθ)}

おむすびの面積は

S4=(S3-S1-S2)x2となる

S=π-ｎS4

L=1+r

S/L^2→?

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

S1→1/2・(θ^3/3+・・・)→0

n・S1→ｎ/2・(θ^3/3+・・・)→0

S2→1/2・(1-r)・(R-tanθ))

n・S2→n/2・(1-r)・R-1/2・(1-r)・(n・tanθ)→n/2・(1-r)・R-1/2・(1-r)・π

n・tanθ→π

S3→1/2・R^2arctan{(1-r)/(R)・{1+(rsinθ)/R+(rsinθ)^2/R^2+・・・}

=1/2・R^2・{(1-r)/(R)・{{1+(rsinθ)/R+(rsinθ)^2/R^2+・・・}-1/3{(1-r)/(R)}^3{1+(rsinθ)/R+(rsinθ)^2/R^2+・・・}^3+・・・}

=1/2・(1-r)・R+1/2・(1-r)(rsinθ)+1/2・(1-r)(rsinθ)^2/R・・・-1/6・(1-r)^3/R+・・・}

n・S3→n/2・(1-r)・R+n/2・(1-r)(rsinθ)-n/6(1-r)^3/R+・・・→n/2・(1-r)・R+1/2・(1-r)rπ-1/12・(1-r)^3π+・・・

n・sinθ→π

n/R→π/2

r→3-√8

n・(S3-S1-S2)→1/2・(1-r)(1+r)π-1/12・(1-r)^3π =1/12・(1-r){6(1+r)-(1-r)^2}=1/12・(1-r){5+8r-r^2}π

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

2nS4=1/6・(1-r){5+8r-r^2}π→1/6・(-2+√8)(12-2√8)π

S=π-2nS4=1/6・(6-(-40+16√8)=1/6・(46-16√8)

S/L^2=S/(4-√8)^2=S/(24-8√8)=1/24・(5-2√2)

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

弧の中点における接線の長さを求めておきたい

弧の中点までの距離は

L=(1+R^2)^1/2-R=1/cosmθ-tanmθ=(1-sinmθ)/cosmθ

弧の中点は(Lcosmθ,Lsinmθ)

接線の方程式は

y=-1/tanmθ(x-Lcosmθ)+Lsinmθ=-1/tanmθ・x+Lcosmθ/tanmθ+Lsinmθ

=-1/tanmθ・x+L/sinmθ

=-1/R・x+L/sinmθではなく・・・・

(rcosθ、rsinθ)と(rcos(π-2θ),rsin(π-2θ))＝(-rcos2θ,rsin2θ)を結ぶ直線

y=(rsin2θ-rsinθ/(-rcos2θ-rcosθ)・(x-rcosθ)+rsinθ

y=-sinθ(2cosθ-1)/(cosθ+1)(2cosθ-1)・(x-rcosθ)+rsinθ

y=-sinθ/(cosθ+1)・(x-rcosθ)+rsinθ

y=-tan(θ/2)・(x-rcosθ)+rsinθ

y=-tan(θ/2)・x+rcosθtan(θ/2)+rsinθ

(x-1)^2+(y+R)^2=R^2

x^2-2x+1+y^2+2Ry=0との交点は

x^2-2x+1+(-tan(θ/2)・x+rcosθtan(θ/2)+rsinθ)^2+2R(-tan(θ/2)・x+rcosθtan(θ/2)+rsinθ)=0

x^2{1+(tan(θ/2)^2}+

x{-2-2tan(θ/2)(rcosθtan(θ/2)+rsinθ)-2Rtan(θ/2)}

+1+(rcosθtan(θ/2)+rsinθ)^2+2R(rcosθtan(θ/2)+rsinθ)

の小さいほうで計算される

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

極限では

R=tan(mθ)=cot(θ/2)→∞

x^2{1+(tan(θ/2)^2}+

x{-2-2tan(θ/2)(rcosθtan(θ/2)+rsinθ)-2Rtan(θ/2)}

+1+(rcosθtan(θ/2)+rsinθ)^2+2R(rcosθtan(θ/2)+rsinθ)

→x^2+x{-2-2Rtan(θ/2)}+1+2R(rcosθtan(θ/2)+rsinθ)=0

→x^2+x{-4}+1+2r+4r=0

→x=2-(4-6r-1)^1/2

弧の長さは2((x-(rcosθ-2rcos2θ)/2)^2+(y-(rsinθ+rsin2θ)/2)^2)^1/2

→2{(x-r/2)}^2+y^2}^1/2

yは-R+{R^2-(x-1)^2}^1/2

=-R+R{1-(x-1)^2/R^2}^1/2

=-R+R{1-1/2・(x-1)^2/R^2-1/8(x-1)^4/R^4-・・・}→0

弧の長さは2|x-r/2|となる

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

比較する弦の位置が少しわかったことにより、計算結果もわずかに変化が見られた

n=3: 0.39678

n=5: 0.311747

n=7: 0.292665

n=9: 0.286286

n=11: 0.285617

n=13: 0.285228

n=15: 0.284992

n=17: 0.284826

n=19: 0.284699

n=21: 0.284622

n=23: 0.284552

n=25: 0.284502

n=27: 0.284506

n=29: 0.284461

n=31: 0.284404

n=33: 0.28446

n=101: 0.284583

n=201: 0.28335

n=301: 0.284443

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

R=tan(mθ)=cot(θ/2)

sinmθ=cosθ,cosmθ=sinθ

L=(1-cosθ)/sinθ=tan(θ/2)

L/R→tan(θ/2)/cot(θ/2)→0

LR→tan(θ/2)/tan(θ/2)→1

A=1/(sinmθ)^2→1

-B=-{-2-L/(Rsinmθ)}→2

B^2-AC{-2-L/(Rsinmθ)}^2-1/(sinmθ)^2(1+(L/sinmθ)^2+2R(L/sinmθ)}→4-2

x=2-√2

y=0

n→∞のとき

弧の長さの2乗は4(2-√2)^2=4(6-4√2)に収束する

これは(1+r)^2に等しい→S/L^2=1/24・(5-2√2）

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝