KO Gem 52年ぶりの極小観測
26日夜、mei/nekoさんがKO Gemの極小観測に成功なさいました!
Hoffmeisterは1965年1月5日にたった一回の極小観測をしただけです。実に52年ぶりの極小観測です。
23日にも試みてくださったのですが、観測に穴があいて極小をとらえられていませんでした。
ありがとうございました!
これが23日の。
こちらが26日の
素晴らしいです。Ngaさん、来年のBulletinに書いてくださいね!
昨日EASolverのスクリプトをブログに書きましたが、僕がおかしいと思っていたところのミスを、Katさんがご指摘くださって、直しました。ありがとうございました!カッコが1つ足りなかったのでした。その結果、その後の行で幾つか変更しなければならなかったので、一応、全スクリプトをもう一度載せておきます。
# mini: vector with m elements of minima(HJD) in ascending order
# Period to search between P1: minimum period & P2: maximum period
P1 <- 2
P2 <- 20
m <- length(mini)
d <- mini - mini[1]
d[m]
N1 <- floor(d[m] / P2)
N2 <- ceiling(d[m] / P1)
k <- N2 - N1
dev <- vector("numeric",k+1) # deviance
for (n in 1:(k+1)) { ここの( )が1組足りなかった
P <- d[m] / (N1 + n - 1)
N <- round(d / P)
res <- lm(d ~ N)
dev[n] <- deviance(res)
}
min(dev)
N0 <- N1 + which.min(dev) - 1 -1が入ります
N0 # N for the last minimum
P0 <- d[m] / N0 # found Period
N <- round(d /P0)
res <- lm(mini ~ N)
coef(res) # improved epoch & Period
residuals(res) # residuals for each minimum
ord <- order(dev)
P <- d[m] / (ord[1:3] + N1)
P # the best 3 periods