Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

アステロイド曲線

アステロイド曲線

アステロイド曲線(asteroid)はギリシア語の「星のような形」に由来する図形で、星芒形あるいは星形ともよばれます。また4つの尖点を持つことから、四尖点形(tetracuspid)の名称でよばれることもあります。アステロイドは、図のように半径 a の大円の内側を半径 a/4 の小円が滑ることなく転がるときに、小円の円周上の点が描く軌跡として定義されます。

Excel アステロイド曲線の面積と弧の長さ

アステロイドの代数曲線
x2/3+y2/3=a2/3
によって与えられ、t を媒介変数として
x=acos3t,y=asin3t(a>0)
と表すこともできます。媒介変数 t
t:0π2π32π2π
を動くとき x,y
x:a0a0ay:0a0a0
を動いてちょうど一周します。アステロイドはスーパー楕円とよばれる曲線
|xa|k+|yb|k=1
において k=2/3,a=b とおいた形とみることもできます。

アステロイド曲線の面積

アステロイドの媒介変数表示
x=acos3t,y=asin3t(a>0)
を用いて曲線に囲まれた内側の面積を計算することができます。第 1 象限の面積を 4 倍すればいいので、
S=40aydx=4π20ydxdtdt
dx/dt を計算すると
dxdt=3acos2t(sint)
となるので、
S=12a20π2sin4tcos2tdt=12a20π2sin4t(1sin2t)dt=12a2(0π2sin4tdt0π2sin6t)
ここで an=0π2sinntdt とおくと、漸化式
a0=π2,a1=1,an=n1nan2
によって定積分を計算できます。a6=56a4 なので、
S=12a2(a456a4)=2a2a4=2a23412π2=38πa2
となります。

アステロイド曲線の長さ

x,y が媒介変数 t で表されたときの曲線の長さは
L=ab(dxdt)2+(dydt)2dt
で計算できます。この公式を使ってアステロイド
x=acos3t,y=asin3t(a>0)
の1周分の長さを計算してみます。
(dxdt)2+(dydt)2=9a2cos2tsin2t
となるので、
L=40π23acostsintdt=12a[12sin2t]0π2=6a
となります。

Excel でアステロイド曲線を描く

 アステロイド曲線は内トロコイド(hypotrochoid) とよばれる曲線群に含まれます。内トロコイドは
x=(rcrm)cosθ+rdcos(rcrmrmθ)y=(rcrm)sinθrdsin(rcrmrmθ)
で表され、rm=rd, rc=4rd とおいたときに、アステロイド曲線の方程式が得られます。

当サイトからトロコイドを描画するExcelワークシートをダウンロードできます:

≫ ファイルをダウンロードする

ダウンロードしたらファイルを開いて [内トロコイド] のシートを選択してください。
デフォルトでアステロイド曲線が描かれています。

Excel アステロイド (asteroid)

一般に、rm = rd, rc = 4rd に設定したときにアステロイド曲線となります。
たとえば、rc = 2, rm = rd = 0.5 としたときにもアステロイド曲線になります。
rc, rm, rd はスピンボタンで自由に変更できます:

Excel スピンボタン

興味のある人はスピンボタンを適当にいじって、アステロイド以外にも色々と試してみてください。たとえば、デフォルト設定から rc を 0.1 下げて 3.9 にするだけでも、アステロイドとは全く異なる曲線が描かれます:

Excel 内トロコイド (hypotrochoid)

天動説とアステロイド曲線

アステロイド曲線は大昔の天文学において「天動説」を説明するために発見された形です。地球を中心に惑星が回っていると仮定すると「惑星の逆行」という奇妙な動きを説明できません。そこで惑星が小さな円を描きながら地球の周りを動くという「周転円」という考え方で何とかこじつけます。周転円にもとづく天動説は数学的には矛盾がないので、「地動説」支持者にとっては手ごわい理論でした。数学の定理が間違った理論を支えていたこともあるのです。「周転円」の考え方がなければ、人類はもっと早くに「地動説」を受け入れていたかもしれません。

エクセルや数学に関するコメントをお寄せください