ORBITLIFE V2.0.7 ORBIT DETERMINATION

練習ページ  ·  220P/McNAUGHT  ·  観測1,618件

はじめての軌道計算を、最後まで。

実在の彗星 220P/McNaught を題材に、観測データを取ってくるところから軌道が決まるまでを 1ステップずつたどります。途中で「これは壊れたのでは」と思う場面が一度だけ出てきますが、 そこが軌道計算のいちばん面白いところです。やめずに読み進めてください。

§00 MAP

全体の流れ。

軌道計算は、7つの手順を順に踏むだけです。難しいのは理屈ではなく、 途中で出てくる数字に驚かないことのほうです。

STEP 1—2
観測を取り込む[MPC から観測データを取得] で、天体名を入れるだけ。
STEP 3
3点を選ぶ離れた3点をクリックすると、仮の軌道がひとつ決まる。
STEP 4
驚く古い観測の残差が桁違いに大きい。ここで多くの人がやめてしまいます。
STEP 5
[手動] で強引に押していく1回ずつ押して、残差が降りてくるのを見る。
STEP 6
リジェクトする年で分けた2つのしきい値。なぜ2つなのかを説明します。
STEP 7
仕上げて保存する[手動] で収束させ、結果を残す。
§01 LOAD

まず観測を手元に取り込みます。

STEP 1

OrbitLife を起動する

特別な準備は要りません。惑星暦さえ受信してあれば、起動してすぐ計算に入れます。

STEP 2

MPC から観測データを取得する

[ファイル] → [MPC から観測データを取得] を選び、天体名として「220P」と入力します。 旧版ではあらかじめ用意した固定ファイルを読む形でしたが、いまは Minor Planet Center から その天体の観測をその場で取ってこられます。数秒で 1,618 件が一覧に並びました。

OrbitLife のメイン画面。左に軌道決定・軌道改良・リジェクトの設定、中央上に軌道要素の欄、中央下に1618件の観測一覧が並んでいる。残差の欄はまだ空欄
読み込み直後。残差の欄はまだ「-----.-」です。軌道がまだ無いので、 比べる相手がありません。

画面の見かた  ·  画面写真はクリック(タップ)で原寸に開きます

左の列が設定です。上から[軌道決定(3点選択)][軌道改良][改良][リジェクト]と 並んでいて、これが作業の順番そのものになっています。中央上が求まった軌道要素、 中央下が観測の一覧。それぞれの[ⓘ]を押すと、その場で解説が開きます。

§02 THREE POINTS

3点を選ぶと、仮の軌道が1つ決まります。

望遠鏡に写るのは天体の「方向」だけで、距離は写りません。それでも、ほどよく離れた3つの方向が あれば軌道はひとつに絞れます。200年前から変わっていない手順です。

STEP 3

一覧から3点をクリックする

クリックすると、その行の[印]欄に * が付きます。画面のいちばん下には、 いま選んだ観測の番号・光度・観測所名が出ます。一覧のいちばん最後(最終観測)から始めて、 古いほうへさかのぼりながら3点選びます。

観測一覧の末尾、2026年9月4日の最後の行を選択したところ。印の欄にアスタリスクが付いている
1点目。一覧のいちばん最後、2026年9月4日。ここが最終観測です。 まずここから選びます。
観測一覧の2026年8月14日の行を選択したところ。画面下部のステータスバーに観測番号1438、光度13.2G、観測所L65 Bredenkamp Observatory, Bremen と表示されている
2点目。そこから3週間さかのぼって 2026年8月14日。下端に 観測番号 1438 / 光度 13.2G / L65 Bredenkamp Observatory, Bremen と出ています。 どこの誰が撮った1点なのかが分かります。
3点目は、さらにさかのぼって 2026年7月11日を選びました。

POINT — 3点の選び方

なるべく新しい観測から選びます。順序は、1. 最終観測(一覧のいちばん最後)  ·  2. それより 2週間〜1か月前  ·  3. さらにその 2週間〜1か月前。ここで使ったのは 2026年9月4日・8月14日・7月11日の3点です。
間隔を空けるのは、時間的に近すぎる3点は空の上でほぼ一直線に並んでしまい、そこから 決まる軌道が少しの誤差で大きくぶれるからです。三角形は、頂点同士が離れているほど 形が崩れにくいのと同じ理屈です。
では、なぜ古いほうへ広げず、新しいほうに寄せるのか。このあとの摂動計算は、 最終観測より後の「40 で割り切れる日」を元期とし、その日の軌道要素として いま求めた要素を置くところから始まります。3点が最終観測に近いほど、元期までの 距離が短くて済みます。100年前の3点で決めてしまうと、現在まで100年ぶんの元期を 移動させる積分計算が余計に必要になり、そのぶん待つことになります。
観測を広く使うのは、このあとの改良の段階です。最初の3点は、新しいほうへ。

3点そろうと、その場で暫定軌道が計算されて中央上に出ます。

3点から暫定軌道が求まった画面。軌道要素とともに軌道分類が木星族彗星 TJ=3.00 と表示され、2026年7月付近の観測の残差は1秒角前後
q = 1.5591987 AU、e = 0.5003751、i = 8°.1258、P = 5.51年。 軌道分類: 木星族彗星 (TJ = 3.00) まで自動で付きます。 選んだ3点の近く(2026年7月)の残差は、もう1秒角前後です。

行に付いている淡い色分けは、観測所ごとの傾向を見るためのものです。ここでは使いませんので、 気にせず進めてください。

§03 THE SHOCK

ここで、たいていの人が「壊れた」と思います。

STEP 4

一覧の先頭までスクロールしてみる

いまの解は 2026年7月〜9月の、たった2か月足らずの3点から作ったものです。 では、いちばん古い 2004年の観測はどう合っているのか。一覧をいちばん上まで戻してみます。

観測一覧の先頭、2004年5月の観測。残差RAが36380秒角、残差Decが8400秒角前後という桁違いの値が並び、赤枠で強調されている
残差 RA が 36380″、残差 Dec が 8400″。 秒角ではピンと来ませんが、36380″ は約10度。満月20個ぶんです。

これは失敗ではありません

2か月の弧から引いた軌道を、22年前まで外挿しているのですから、10度ずれるのはむしろ 当たり前です。壊れているのは軌道ではなく、まだ古い観測を1点も使っていないという ことだけです。ここで「うまくいかなかった」と思ってやめてしまうのが、 いちばんもったいない。

そしてもうひとつ大事なこと。この段階でリジェクトしてはいけません。 いま2秒角で切れば、2004年から2025年までの観測が丸ごと消えます。残るのは選んだ3点の まわりだけ。それは「よく合う解」ではなく、「合わない観測を全部捨てた解」です。

§04 FORCE IT

[手動] を、強引に押していきます。

STEP 5

1回ずつ、残差が降りてくるのを見る

[改良]の枠にある [手動] を押します。1回押すと1回だけ改良します。何も選び直さず、 何も捨てず、ただ押します。10度ずれていた古い観測が、押すたびに一気に引き寄せられてきます。 数回で、全期間が1秒角の世界に入ります。

POINT — 軌道改良は [手動] が通常です

軌道改良は、ふだんはこの [手動] だけで進めます。1回押すごとに残差が少しずつ 下がり、多くの場合は2〜3回で収束します。収束したら残差の傾向を見て、相対的に 大きいものをリジェクトし、また [手動] を2〜3回。この繰り返しで計算は終わります。
1回ずつ止まるので、平均残差がちゃんと下がっているか、要素があらぬ方向へ飛んで いないかを確かめながら進められます。
[自動] は「楽をするための [手動]」ではありません。観測数が極めて少ないときと、 ミニムーンのとき——この2つの場面のための道具です。§07 で説明します。

手動改良を数回行ったあとの画面。2004年の観測の残差が0.3から2.2秒角に収まり、平均残差1.083931秒角と表示されている。右側に各惑星への最接近の一覧が出ている
[手動] を数回押したあと。2004年の残差が 0″.3〜2″.2 まで降りて、 Mean residual 1.083931 arc seconds。22年ぶんの観測が、1本の軌道に乗りました。
MEAN
軌道要素タブの Mean 0.00414 arc seconds は、この改良で要素がどれだけ動いたかを表します。これが十分小さくなれば収束です。
MEAN RESIDUAL
こちらが平均残差。いま使っている観測が、この軌道と平均でどれだけ食い違っているかです。
最接近
右の欄に、各惑星と月への最接近が自動で並びます。木星に 3.46 AU まで近づいていることが分かります。木星族彗星らしい姿です。

押しても数字が変わらないときは、[計算ログ]タブを見てください。要素変化量を自動で 調整した、といった経過が書いてあります。もう一度 [手動] を押せば続きから進みます。

§05 REJECT

リジェクトの欄が2つある理由。

左下の[リジェクト]の枠には、数字を入れるところが3つあります。はじめて見ると 「なぜこんな指定のしかたなのか」と思うはずです。ここには理由があります。

STEP 6

年で分けて、しきい値を2本立てる

リジェクトの設定欄の拡大。「1800 年以前 999 秒角」「以降 2.0 秒角」と入力され、下に[リジェクト]ボタンがある
境目の年と、その前後それぞれのしきい値。

読み方はこうです。

1800 年以前
1800年までの観測は、残差がこの値(″)を超えたら外す。
以降
1800年より後の観測は、残差がこの値(″)を超えたら外す。

なぜ年で分けるのか — ハレー彗星を思い浮かべてください

1P/Halley の観測は、紀元前まで遡ることができます。近いところだけ数えても、 1607年のケプラー、1682年のハレー自身、1759年、1835年、1910年、そして1986年。 この間に観測の道具は、肉眼と星図 → 眼視で目盛りを読む → 写真乾板 → CCD と変わりました。精度も、分の単位から秒角の1/10まで、けた違いに変わっています。

ここに一律「2秒角を超えたら外す」を当てたらどうなるか。 古い観測は1点残らず消えます。正しく観測されていても、当時の道具ではその精度が 出ないからです。そして古い観測が消えるということは、軌道を長い時間で押さえる支点を 自分から捨てるということです。

だからしきい値は1本ではなく2本。 「古い観測には古い観測なりの物差しを当てる」ための欄です。 既定値は境目が1820年、古い側が6″、新しい側が2″。1820年にしてあるのは、 写真観測が広まる前の眼視観測をひとまとめに扱えるからで、 天体ごとに計算者が決め直して使うものです。

今回の 220P/McNaught は 2004年以降の観測しかありません。ですから境目の年はどこでも 同じことで、古い側のしきい値も効きません。図では境目を1800年、古い側を大きな値にして 「古い側は使わない」という意味にしてあります。効いているのは 以降 2.0″ だけです。

[リジェクト] ボタンを押すと、しきい値を超えた観測に印が付き、次の改良から外れます。

リジェクト実行後の画面。外された観測の残差が括弧付きで表示され、印の欄にRが付いている。平均残差は0.742197秒角に下がり赤枠で強調されている
外れた観測は残差が括弧付きになり、[印]欄に R が付きます。 平均残差は 0.742197 arc seconds まで下がりました。
印 *
軌道決定に使った3点
印 R
しきい値を超えて自動的に外れた観測
印 F
一覧をダブルクリックして手で外した観測
印 X
観測所コードが見つからない観測
印 Wn
ウエイトを変えた観測(重くした行は紺地に黄文字で塗られます)

外すことは、消すことではありません

リジェクトは印を付けるだけです。しきい値を緩めて [リジェクト] を押し直せば、 いつでも戻ってきます。思い切って厳しくしてみて、行きすぎたら戻す、で構いません。

ただし、切りすぎに注意

厳しくしすぎて観測が極端に減ると、軌道要素がわずかな偏った情報で決まってしまい、 かえって不安定になります。平均残差の小ささだけを追わないでください。 何個外して何個残ったかを、そのつど見てください。

進め方の定石は「緩めに始めて、収束してきたら少しずつ厳しくする」です。 いきなり目標のしきい値を当てず、6″ あたりから 4″、3″、2″ と降ろし、 各段で改良を挟みます。CCD の高精度測定だけなら 1〜2″、眼視や古い写真乾板を含むなら 3〜5″ あたりが目安です。

§06 FINISH

仕上げて、残します。

STEP 7

収束させて、保存する

リジェクトのあとも、[手動] を2〜3回押して収束させます。残差が下がりきり、 外れる観測がもう増えなければ、そこが終わりです。220P/McNaught のように観測が 1600件あれば、[自動] の出番はありません。

[ファイル] → [軌道改良結果を保存] で結果を書き出します。同じ解をあとから そっくり復元したいときは、17桁の再現用ファイルとして保存しておいてください。 どの設定で何回改良し、どこでリジェクトしたかは[計算ログ]と操作ログに残ります。

SOLUTION220P/McNaught  ·  Epoch 2026 Oct. 7.0 TT
T            2026 June 14.11548 TT
q            1.5593046 AU
e            0.5003058
Peri.        180.56086     (2000.0)
Node         150.08417
Incl.          8.12082
P            5.51 年
分類         木星族彗星 (TJ = 3.00)

Mean residual  0.742197"

軌道要素の右上に名前が入っているのは[計算者署名]です。誰が計算した解なのかを 結果に残せます。

§07 AUTO

では、[自動] はいつ使うのか。

[自動] は、[手動] が効かない相手のための道具

ここまで [自動] をほとんど使いませんでした。観測が十分にある天体なら、それで足ります。 [自動] が要るのは、[手動] では進めない2つの場面です。この使い分けは、 軌道改良を行うすべての画面に共通します。

その1 — 観測数が極めて少ないとき

観測が数えるほどしかない天体では、[手動] を押しても残差が大きくなったり小さくなったりを ランダムに繰り返します。下がっていきません。
こういうときは [自動停止]にチェックを入れて [自動] を押します。 それまでで最小の残差になったところで止まります。もっと良い解を探したければ、 もう一度 [自動] を押します。前回止まった残差より小さいものが見つかれば、そこでまた止まります。 これを繰り返して、望む残差の軌道に行き着きます。
計算するというより、引き当てるに近い作業です。運よく非常に小さな残差の軌道が 出た場合でも、計算者が求めていたのはその軌道です。そこで止めて、軌道要素を示してください。
繰り返すほど、次に止まるまでの間隔は延びていきます。数十分待っても最小を更新できなく なったら、そのあたりがその観測から引き出せる限界です。

その2 — ミニムーンの軌道を求めるとき

ミニムーン(惑星に一時的にとらえられた天体)では、1回の改良で残差がわずかしか下がりません。 十分小さくするのに 数百回 [手動] を押すことになります。
こういうときは [自動停止]のチェックを外して [自動] を押します。時間はかかりますが、 残差がわずかずつ下がっていくのが画面で見えます。下がらなくなったら [中止] を押します。

言い換えると、[自動] は「[手動] を楽にするボタン」ではありません。 残差が素直に下がる相手は [手動]、下がり方が素直でない相手は [自動]。 そう覚えておいてください。

§08 NEXT

この先へ。

STEP 8

軌道が求まったところは、ゴールではありません

光度曲線で増光の様子を追ったり、軌道図で太陽系の中の位置を 確かめたり、地球儀の上で最接近点を眺めたり。そこから先の機能は 機能のページにまとめてあります。

もう一歩踏み込みたい方にはミニムーンの地心軌道図を。 2024 PT5 が地球へ巻きついて去っていく軌跡を、観測の読み込みから自分で出すまでの手順です。 強制リジェクトと元期の移動の、いちばん分かりやすい使いどころでもあります。

そして、求まった数値をどう扱うかについては最後の判断のページを 読んでおいてください。とくに非重力効果を解いたときに A₁ が負になった場合の話は、 (1566) Icarus の記録に実例があります。