判定光变周期的一种方法-上海天文台.PDFVIP

　　中　国　科　学　院　上　海　天　文　台　年　刊 ２００４年　第２５期　　 ＡＮＮＡＬＳＯＦＳＨＡＮＧＨＡＩＯＢＳＥＲＶＡＴＯＲＹ ＡＣＡＤＥＭＩＡ ＳＩＮＩＣＡ　　Ｎｏ．２５，　２００４ 判定光变周期的一种方法 陶　隽　钱伯辰　梁兆刚　潘红鉴 （中国科学院上海天文台，上海２０００３０） 提　　　要 ２ Ｊｕｒｋｅｖｉｃｈ方法通常用来发现光变曲线的周期，并且根据 极小值的大小来验证所发现的光 Ｖ ｍ 变周期可靠性。但是，这种验证方法又带有一定的不确定性，最佳的方法是把所发现的周期与观测 光变曲线拟合，得出每一个周期的振幅，以振幅大小来验证所发现的光变周期的可靠性。将此方法 用于ＯＪ２８７的历史光变曲线，得到其最主要的周期为１１．９年，与其它结果相符合。 主题词：光变曲线 — 周期 — Ｊｕｒｋｅｖｉｃｈ方法 — 振幅检验 分类号：Ｐ１５７．７ １　引　　言 发现和研究活动星系核的光变是一项极有意义的工作。较长时标的光变通常是由星系核 的进动或轨道运动产生的，因此可以利用所发现的周期来估计相应的物理参数。 ［１］ Ｊｕｒｋｅｖｉｃｈ方法 （Ｊｕｒｋｅｖｉｃｈ１９７１）和 Ｆｏｕｒｉｅｒ分析都可以用来寻找时间序列的周期。 Ｊｕｒｋｅｖｉｃｈ方法的一个优势就是不要求数据的时间间隔是等间距的，而一些周期分析的方法只 能处理等间隔的数据，这就往往需要通过插值来满足要求，由此带来一些不真实的信息。在实 际过程中大量数据因为各种条件限制不一定是等间隔的时间序列，所以Ｊｕｒｋｅｖｉｃｈ方法在天文 ［２］ 上有较多的应用。Ｋｉｄｇｅｒ等人 （１９９２）认为就发现光变曲线周期而言，Ｊｕｒｋｅｖｉｃｈ方法要优于 Ｆｏｕｒｉｅｒ分析。然而Ｊｕｒｋｅｖｉｃｈ方法在理论上虽能正确地得出光变曲线的周期，但实际的光变曲 线中除了周期性的变化外还有非周期性的光变和噪声，因此我们在应用Ｊｕｒｋｅｖｉｃｈ方法时发现 它仍有一些不足之处，本文将讨论这些不足之处，并提出一些改进的方法。 本文的讨论将用ＯＪ２８７的历史光变曲线作为例子。ＯＪ２８７是有名的活动星系核，它具有 长达一百多年的历史观测资料，其历史光变曲线呈现很强的周期性，并已有很多人对其进行过 研究。我们编辑了ＯＪ２８７从 １８９１．０４年到２００１．９８年的历史光变资料，共 ５０８５个观测。 ［３］ Ｓｉｌｌａｎｐ 等人（１９９６）指出，由于近来二十多年的观测资料要远远多于以前，如果不对这些资 料作一些处理就直接进行统计的话，必然会使最近的资料的权重远远大于以前的，对统计结果 得到平均的光变资料后再进行统计分 带来不良影响。一般按时间间隔对这些资料进行平均， 收稿日期：２００４０３２９；　修回日期：２００４０４２６ 上海市科学技术委员会项目（０２ＱＡ１４０５５）；中国科学院光学天文开放实验室资助课题。 ９８ 中 国科 学 院上 海 天 文 台年 刊　　　　　　　　　　　　　２００４年 析。我们采用０．０２年的时间间隔进行平均，得到６３４个Ｖ星等资料。图１是ＯＪ２８７的历史光 变曲线资料。 图１　ＯＪ２８７的历史光变曲线资料（０．０２年平均值） Ｆｉｇ．１　ＨｉｓｔｏｒｉｃｌｉｇｈｔｃｕｒｖｅｏｆＯＪ２８７（ｔｈｅａｖｅｒａｇｅｖａｌｕｅｏｆｅａｃｈ０．０２ａｄａｔａ） ２　Ｊｕｒｋｅｖｉｃｈ方法 Ｊｕｒｋｅｖｉｃｈ方法是按照不同的试验周期Ｐ把数据折叠，然后把折叠好的资料按其相位从小 ２ 到大排列，然后分成ｍ组。如第ｌ组的离差平方和Ｖ