Archive for 7 月, 2009
Rで部分日食の食の入り/食の終わりの時間を計算してみる
火曜日, 7 月 21st, 2009ふと、部分日食の時間を計算で予測できないものかと考え、調べてみたところ山口大学 藤沢研究室に丁寧な解説がありました(PDF)。このドキュメントでは、日食の計算に必要なベッセル要素を数表から出しているのですが、NASAにベッセル要素の近似式パラメータがあったので、試しにそれを使ってみることにしました。作製した関数は以下の通り。
test <- function(phi0=35.658, lambda0=139.741, h=0) {
# NASAから取ってきたパラメータ (Besselian elements)
# http://eclipse.gsfc.nasa.gov/SEbeselm/SEbeselm2001/SE2009Jul22Tbeselm.html
tanf1 <- 0.0046014
tanf2 <- 0.0045784
deltaT <- 66.5
x <- function(t) {
0.239961 + 0.5563955 * t - 0.0000576 * t^2 - 0.0000094 * t^3
}
y <- function(t) {
-0.003276 - 0.1774579 * t - 0.0001344 * t^2 + 0.0000032 * t^3
}
d <- function(t) {
(20.26424 - 0.007874 [...]
ハクビシン(?) 目撃
金曜日, 7 月 10th, 20097/10 23:50ごろ、マンションの庭から”ケケケケケ、ケケケケケ”と奇妙な声がしました。最初は鳥かと思ったのですが、カミさんと覗いてみると、体長40cmほど、尻尾が20cm程度のタヌキのようなズングリした動物がゆっ [...]