MCMC法でEAの要素を求める際の、パラメータやその入れ方を工夫してみた。

Katさんのご意見では、MCMC法でやる場合にパラメータの数が増えるとそれだけ正しそうな解になかなか到達しないので、EAの場合、副極小はフェーズ0.5と固定してもよいのではということです。今までのスクリプトでは食外から極小への傾きがパラメータとして入…

2015-01-26

EAの副極小のフェーズもパラメータにしてみた

下のスクリプトのようにやってみた。p[8]は副極小が0.5からどれだけずれているかというパラメータ。 fphの代わりになんとかph <- {(x-p[2])/p[1]}%%1だけで出来ないかと abs{ph-round(ph)}とやってみた。しかしこれで計算させたらどうしても数値がおかしくな…

2015-01-25

MCMC法でKaiV20の要素を改良つづき

KatさんがMCMC法でEAのタイプの周期を求めることをすすめて下さり、この間ブログに書いたようにやってみました。 MCMC法でKaiV20の要素を改良 - muttenz's blog そのスクリプトの fph <- 0.5 - abs((((x-p[2])/p[1]) %% 1) - 0.5) yxep <- ifelse(fph < p[3]…

2015-01-22

動物たちの知能

天文とは全く関係ないのだが、昨日3SATテレビで見た番組が面白かったので書いてみたくなった。この番組、これから6日間はここをクリックすると見られるはず。「追記」 Katさんから、開始から３０秒位経過したところでストップしてしまうとのご指摘を頂きま…

2015-01-21

ぱなりさんの新変光星　続き

昨日ブログを書いて寝て、夜中に起きた時にそうか、パナリさんの観測時刻はUTではなくJSTだったのだと気が付きました。なので、この星の減光はUTの21時半頃ではなく、JSTの21時半頃だったのです。うお座がそんな明け方に見えるはずがないのに、時差などから…

2015-01-20

ぱなりさんの新変光星？　続き

前原さんのKWSでデータを見てみました。 http://kws.cetus-net.org/~maehara/VSdata.py?object=001321%2B0545.1&resolver=&p_band=V&plot=1&submit=Send+query ところどころ減光しているように見えます。目立って暗い12.5等のデータはその前後を見るとどうも…

2015-01-18

ぱなりさんの新変光星？　続き

先程、vsolj-obsのデータベースでPSCDVの観測を調べたらかなりの量が観測されています。Njhさん、Iohさん、Ngaさん、Sizさん、Kisさんと皆さん頑張っておられます。撮影されたフレームからぱなりさんの新変光星の変光が検出できないでしょうか？かなりお手…

2015-01-18

ぱなりさんの新変光星？

ぱなりラボラトリー通信で12月31日に「DV Pscのそばに新変光星？」という記事が出ました。ぱなりラボラトリー通信 - Yahoo!ブログぱなりさんが変光に気が付かれてからちょっと日が経ってしまったとのことですが、この変光の具合を見ると周期がほとんど１日…

2015-01-17

AK Eri の要素をKWSのデータも加えてMCMC法で改良

KWSのデータだけでMCMC法を試みたのですが、データ数が少なくてなかなか値が収束しない。そこで、得られた食外光度11.9566とASASの食外光度11.9999の差をKWSのデータの光度に加え、ASASとKWSのデータを一緒にしてからMCMC法で周期などを求めてみた。結果の…

2015-01-16

MCMC法でAK Eriの要素を改良

このブログでエリダヌス座であまり観測されていない食変光星の一つとしてAK Eriを紹介しました。食変光星AK Eri - muttenz's blog 食変光星AK Eri - muttenz's blog その時にGCVS（VSXも同じ）の要素などで作ったphase plotがあまり綺麗にならないので、少…