星空研究Memo

ここは某天文屋の外部記憶装置である。

Julian Day

2010-05-04 23:31:30 | R言語
Rの練習。
JDへの計算。

-----------
Y <- 2010                     # year
M <- 5                          # month
D <- 4.8474                 # day

if (M <= 2)

     {
      Y <- Y - 1
      M <- M + 12
     }

a <- trunc(Y / 100)

b <- trunc(2 - a + (a / 4))

ja <- trunc(365.25 * (Y + 4716))

jb <- trunc(30.6001 * (M + 1))
jc <- (D + b - 1524.5)

jd <- ja + jb + jc

     options(digits=14)     # 14桁表示
jd                                      # 結果

-----------

2010年2月1日で計算すると
2455228.5
てな感じの結果が出ました。合っていそうです。
時間・分・秒を日に直して計算してみると(ex. 2010年5月4.8474日UT)
2455321.3474
と出でた。
概ね良さそうかな。

脈動型の分類

2010-05-04 06:55:24 | 独り言
脈動型の分類記号を見ても、いつもパッとしなく
「あれ?これどんな星や?」
となります(単に勉強不足と怒られそうですが;)。

wikiにそれなりにまとまっていて、大変助かります☆
これはやはり、変光星研究会の何方かがまとめられたのだろうか?!