読者です 読者をやめる 読者になる 読者になる

muttenz's blog

スイス星空だより

KO Gem 52年ぶりの極小観測

26日夜、mei/nekoさんがKO Gemの極小観測に成功なさいました!

Hoffmeisterは1965年1月5日にたった一回の極小観測をしただけです。実に52年ぶりの極小観測です。

23日にも試みてくださったのですが、観測に穴があいて極小をとらえられていませんでした。

ありがとうございました!

これが23日の。

f:id:muttenz:20170127061151p:plain

こちらが26日の

f:id:muttenz:20170127061223p:plain

素晴らしいです。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