ケプラーのデータを解析 きれいです！
これの続きです。
Period04で周波数を15個まで探してその周期で計算したライトカーブと観測値がかなり良く合うのでもう少し周波数を探させました。
それから変動の強い周波数のdetrendの位相図を作ることにしました。
ここでちょっと面白いのは、F1とF2の二つが周期にしてそれぞれ3.5日と7日強で、F1はF2のほぼ二倍ながら完全に2倍ではないのです。
まずF1(=P: 3.4851618428464 d.)detrendedの位相図です。
F2(=P: 7.0959151462743 d.)の位相図です。
ここでF1F2を同時にdetrendしてみたらどうなるかやってみました。いわばF1 F2以外の変動を取り除いた変動です。
ただし位相図を作るための周期を何にするか問題です。F1に相当する周期か、F2か？
まず思いつくのはF2の周期です。（F1とF2の両方の変動が入っているわけですからF2の長い周期にしなければ無理でしょう。）
それで位相図を作ると。
これはダメだ！
F1の周期の2倍にしてみよう。
これは綺麗です！この位相図は何を意味するのでしょうね？まるでEWみたいに見えますが。。。
同じことを僕があしかけ3年かけた観測データでやると。。。残念ながらこんなです。。。
比較する事自体無駄なことですが、自分の観測は一体なんだったのだろう。少々がっくり来ます。。。
周期不明の食変光星V563 Her の周期を探そう 続き
これの続きです。
一昨日の夜、久しぶりに晴れました。薄雲がある上に月も明るいので空全体が相当明るいのですが、23時半にあるはずのV563 Herの極小をとにかく観測してみました。ようやく極小時ごろに隣の家の屋根上に出てきたので観測開始。15等台と暗く、空が明るいのでスキャッターはひどいですが、極小後からの増光をはっきりと捉えられました。Hoffmeisterの2回めの極小観測が1964年だったので、おそらく53年ぶりの観測になるでしょうか？2時頃には完全に曇られてしまいました。
3月に見つけた周期がそれほど外れていないことが、今回の観測で検証できました。これから何度か第一、第二極小を観測して周期と元期をもっと正確に求めることが課題です。周期が6日あまりなので第一極小の次のチャンスは5月12日夜半すぎ。（11日から12日かけて）（第二極小は明日8日夜。天気は悪そう）
ASAS113804-0026.5のケプラーデータがあった！すごい！
セバスチャンからメールでこの星のケプラーデータがあるよと教えてくれました。
全体で80日弱の観測です。
このライトカーブをパッと見た時には、あ２晩ぐらいの観測だと思ったのですが、80日間の連続観測です。昼も夜もないし、空に雲はないし、こんな観測が出来るとは羨ましいです。
このデータをまずPDMで解析してみました。
周期0.01から0.1の範囲で周期0.02..が出ました。位相図。
追記 これを書いてから考えたら、ケプラーは0.02日ぐらいの間隔で観測しているのでその周期が出てしまったようです。
周期1日から10日の範囲でPDMで解析してみました。6.97..日という周期が出ました。
これはASASデータでも見つかっていた周期で、その位相図は、上とほとんど同じものができました。
ケプラーデータを今度はフーリエ解析してみました。
最初に出て来るスペクトルです。
ここからF15までやってみました。僕が見つけた周期とかなり合っています。
ところどころ違うのもありますが。
ここまでやってライトカーブに計算したカーブがどのくらいの合うか見てみました。
見事ですね！これの凹んだあたりを拡大してみると。
もう嫌になってしまう。すごい精度です。ただ連続測光という点ではこれはかなりまばらで、0.02日ぐらいの間隔です。
Period04を使ったフーリエ解析でPTS (periodic time shift) オプションを使う効果
昨日PTSを使った場合に見つかる周波数が変わったと書きましたが、グラフにしたほうがその効果、変化がわかりやすいかと思い作成しました。
周波数対amplitudeのグラフというのはライトカーブのスペクトルに当たるわけですね。
PTSなしの場合のスペクトル。
PTSありの場合のスペクトル。
10c/d以上の部分だけ拡大すると。
PTSなし。二つの周波数が同じ強さで突出しています。
PTSあり。一つだけが突出しています。
違いがはっきりわかるでしょうか？
ASAS113804-0026.5の「Periodic time shift」を使用したフーリエ解析
すでに先日ブログに載せた最初のフーリエ解析の結果です。
F9の16.3359がF4の16.27595とたった0.06弱しか違わないので、F9はperiodic time shiftのためでリアルではない可能性があるということで、expert modeでPTSの公式を使ったフーリエ解析をしました。
その結果を載せます。
F1からF8までの周波数は上の表とほとんど変わりませんが、amplitudeでは例えばほとんど同じ強さだったF2とF4が、F4は強くなり、F2は弱くなり差がはっきりしました。
F17以降が新しく計算しなおして見つかった周波数です。
周波数が10以上、周期にして0.1日以下で見つかった周波数とamplitudeをグラフにしました。強い順にF4 F2 F13 F19 F17 F18 F20 F22 です。
周波数10以下、周期にして0.1日以上のグラフです。圧倒的にF1が強いです。
PTSで見つかった周波数は上記の表中一番上PTSFとある0.0088 c/dで、周期にして113.5日です。この周期が物理的に何を意味するのかはわかりません。
下は表の全部の周波数とamplitudeをグラフにしたものです。
Period04のexpert mode
Oｓｍさんが早速エキスパートモードでPTSに切り替えるやり方を教えて下さいました。どうもありがとうございました！
その結果は次の機会に載せることにして、昨日のフーリエ解析の結果から作った位相図を載せます。２番めからの位相図、どれも同じように見えますが、周期と振幅が異なるのです。
F1のdetrended これが出来るのを長いこと待っていました！やっと３シーズン目でできた。観測した回数は延べ２９回です。
F2のdetrended
F4のdetrended
F7も十分のamplitudeがあります。
ASAS113804-0026.5の多重周期（フーリエ解析）
今年の観測データだけでフーリエ解析をやってみて、ある程度やり方に慣れたし、意味のありそうな多重周期が見つかるようなので、思い切って2015年、2016年のデータも付け加えてPeriod04でフーリエ解析をしました。
大変に嬉しい結果が出ました！
数字ばかりでちょっとおもしろくないかもしれませんが載せてみます。
まずF1です。なんと周期がASAS-3データからPDM解析で得られた周期3.483日に近い値3.485日がようやく出てきたのです。これでこの周期は確定でしょう。
そしてDSCTの多重周期がいくつか見つかりました。カッコ内は周期と振幅です。
F2(0.05804 d.0.0082mag)、F4(0.06144 d.0.0082mag)、F7(0.07546 d.0.0036mag)
F2とF4の振幅がほぼ同じなので時期によってはどちらかの周期が強く見えることがあるのでしょう。
問題はF9、周波数は16.33590668 /d.です。
ちなみにF2の方は 16.27595494 /d.です。
その差は僅かに0.06 /d弱！
Period04のマニュアルにこういう時は「注意すること」とあります（P60）。
このように近接した周期が見つかる場合は周期は基本的には変わらずほんの僅かに周期的にフェーズがずれていることがあるそうです。（例えばDSCTが連星で周期的に見かけ上フェーズがずれるような場合。）そういう時はPTS(periodic time shift)のオプションを使うと解決できるとあります。
ところがそのオプションを使うためにはexpert modeにしてそのオプションを選択するとあるのですが、expert modeにしてもそのオプションが出てこない！どなたかご存知ですか？とりあえず、今日はここまで。
