雑記
« Previous EntriesRで部分日食の食の入り/食の終わりの時間を計算してみる
火曜日, 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程度のタヌキのようなズングリした動物がゆっ [...]
汎用ロジックICでD級アンプ (3)
火曜日, 5 月 19th, 2009この辺(その1、その2)の続き。
74AC04が出す100MHz程度のリンギングが、スパイクノイズの原因であることが分かりました。RCスナバ(10Ω+1000pF)やフェライトビーズを入れて対処しようとしたものの、なか [...]
ダイオード2つで三角波から正弦波を近似
火曜日, 5 月 12th, 2009XR-2206の後段にあるsine wave shaperに興味を持ち、似たようなものを作れないかと考えて調べていたら、ダブリン工科大学の講義資料を発見。XR-2206に近いものはスライド26にある差動増幅器ですが、ス [...]
汎用ロジックICでD級アンプ (2)
土曜日, 5 月 9th, 2009この辺の続き。
現在、こんな感じです。
こいつに8Ωのダミーロードを付け、WaveGeneで作った1kHz (正確には65536ポイントFFT用の1000.488kHz)を突っ込んだときの波形です。100mVくらい [...]
「電子工作の広場SNS」に入会しました
水曜日, 4 月 22nd, 2009以前から気になっていた電子工作の広場SNSに入会してみました。ハンドルネームはここと同じ”kmgs”です。秋月やらPIC、自作オーディオ等、面白そうなコミュニティがあったので、早速参加してみました。
SNSは初めてなの [...]