MdyさんのASASSN-14acのデータを周期解析してみた

vsolj-obs 43788でMdyさんがASASSN-14acの観測報告をなさっていて、しかも彼のブログでそのデータのきれいなライトカーブが見られたので、練習のためにと周期解析を試させていただいた。

観測報告をRでjd_magのデータに変換した。それをプロットするとこうなります。

これを見ると、全体に弧を描いているので、Katさんのスクリプトにあるlfitnを使って２次曲線近似で平らにしてみた。

ところどころ、外れたようなデータがあるので（特に終わりの方で）、＋－0.08以内のデータだけを残した。（幅狭く切り取りすぎたでしょうか？）

このデータをPDMで周期解析。そのグラフは

二つピークがあります。それぞれ、

0.0467 d. と0.0940 d.でした。

これらの周期でphase plot を描くと

0.0467 d.

0.0940 d.

vsolj-alert 16830でKatさんが書いておられたところでは、シングルピークの場合は周期0.0472あるいはダブルピークならその2倍ということでした。

練習にやってみて、意外に良い値が得られて喜んでいます。