作者:孙国佑
未经允许禁止转载!
本节介绍如何得到变星的各种参数。
一,计算周期。可以用多个软件做傅里叶计算,相互验证参考。(演示变星为我发现的XOSS-V J185558.93+720314.4)
1,P_search15f,下载地址:https://www.variablestarssouth.org/software-by-bob-nelson/,点击Period Search下载。 |
点击Open Data File导入测光数据,需要csv格式,一列时间,一列星等。 |
点击Fourir傅里叶计算(这个软件还提供了其他2种计算方法,自行尝试不在赘述,其他两个结果可以做相互验证),MIN MAX,周期区间可以调,方便选取。一般选取峰值即可,即可得出曲线,但多数情况还需要手动调整,这个软件最大的缺点就是非自动化,还需要手动选取,调整。有时候数据不太好,可能离峰值较远,才能选出像样的曲线来。 |
2,Period04,下载地址https://www.univie.ac.at/tops/period04/index.html,点击Period04-v1.2.0-Setup.exe下载。(个人建议用Period04结果作为周期计算基准,因为软件是自动的,然后用其他软件进行相互验证。) |
点击Import time string导入数据,数据为TXT格式,文件路径不能有中文。文件为一列时间,一列星等。 |
选择Fourier标签,点击Calculate,得出周期,振幅数据。F的倒数为周期,振幅A的2倍为变幅。所以这个数据最终结果是周期是1/F=1/14.7070125≈0.06799天,光变幅度为2*A=2*0.0269876717≈0.54星等。 |
3,vstar。新发现一个软件AAVSO开发的vstar(需要JAVA支持),极其小巧且强大,可以计算周期,画相位图。点击下载vstar。详细操作可以阅读VStar手册。下面我做简单介绍,详细使用需要自己琢磨。建议以vstar计算为基准,用其他周期计算软件做验证(最重要的是一个食双星跟变星的区分,食双星周期是变星的2倍。) |
点击New Star from file导入测光数据,CSV格式,一列时间,一列星等。 周期又多种方式计算,我习惯用蓝色框这个, 点ok 点Top Hits 选择周期 选中Power最大的一行周期,如图所示,period=0.209916天,点击New phase plot 画相位图。 这就是画出的相位图。实际操作过程中,要注意食双星跟其他变星的区别,假如这是一颗食双星,那么周期就应该乘以2,我们可以手动输入周期,进行画图。如图: 关于相位图,我只是做了演示,实际情况,食双星画2个周期,极小开始。变星画2个周期,极大开始。所以要搞清楚极大极小时间,见下面章节介绍。 |
二,变星分类
可根据光变曲线形态,周期,振幅,光谱型等参数进行分类,可参考本人撰写的变星的分类一文。 |
三,极大极小时间
软件Minima25d,下载地址:https://www.variablestarssouth.org/software-by-bob-nelson/,点击Minima下载。 |
运行Minima25c.exe,点File-Open,打开csv文件(注意:只需要第一次出现极大,或者极小的那一段测光数据即可,不需要全部测光数据)。极大最好用Fourier,极小值最好用Kwee and van Woerden。因为测周期时的那个演示变星,这个软件不好测,我用了其他的星,做演示。这个测的是极小值,选第4个Kwee and van Woerden,点go,会跳出一个曲线框框,右击鼠标选中你需要测量的区间,红框所示,然后点击Fit Kwee,测量极小时间。结果为2458233.30640,这个就是极小时间。 |
这个软件感觉极小值比较好用,极大值有时候不好使,而且软件经常出错,但目前没有其他软件可用,只能用这个先。我测极大的时候,Fourier用不了选不了测量区间,我手动选取了一段数据,然后用其他方法去算,这个方法大家可以借鉴,算不出来,或者出错,可以调整下数据只保留你需要的那一段数据,还有就是用其他的测量方法给你的测量进行参考佐证。具体的需要自行摸索。 |
重要提示:食双星相位图极小开始,画2个周期,Epoch对应为极小时间,所以食双星是计算极小时间。变星相位图是从极大开始,画2个周期,Epoch对应为极大时间,所以变星是计算极大时间。 |
四,相位图
画相图可由vstar获得! |
五,极小极大星等
也就是光变的星等区间,本文开头就提到测光使用CMunipack,不要用较差,直接用标准星星等,得出准确星等光变曲线。就是为了这个光变区间的确定。 |
因为没有相应的软件进行计算,所以我们可以进行估算。我建议估算的时候,用vstar抽出曲线,来估算星等区间,较为准确,同时要考虑各种软件得出振幅A值。 |
六,Rise/eclipse duration
如需了解VSX上报内容,联系站长邮箱dofur@hotmail.com。或加星明天文台qq群73444198,联系站长。 |