■

CRTSからのデータはHJDになおしてからペランソPeranso（Period analysis software)で周期解析をやっています。どんな周期か全くわからないときはANOVAで0.1ｄから0.5ｄを4000分割して計算させてこれはと思える周期のあたりをもっと詳しく計算させてますが、これがあまり良く当たらない！いっぺんに40000分割ぐらいにして計算させると正しそうな周期が出やすくなるがものすごく時間がかかる。

今朝は思いついて、KatさんのR言語による周期解析、「pergrm」を使ってみた。これだと計算が速い！0.1-0.5ｄを４０００００（四十万！）分割して計算させたら30秒くらい計算していたけど、一発で正しい周期を出した。下の例はASAS111707+0402.9。データはCRTSのCSS。

pergrmの場合、EWのように山が二つあるようなライトカーブの周期解析はひとつ下の周期を求めてから２倍したほうが良さそう。上の図から0.170193ｄが一番強そうなのでその２倍の周期でPhase Plotを作ると下のようになる。あまりスッキリしないライトカーブだが。。。周期はASAS-3の周期0.340386ｄと完全に合っていた！

他のピークでももう少し探してみようか？