判定光变周期一种方法.pdfVIP

维普资讯 中 国 科 学 院 上 海 天 文 台 年 刊 2004年 第25期 ANNALSOFSHANGHAIOBSERVATORYACADEMIASINICA No．25，2004 判定光变周期的一种方法木 陶 隽 钱伯辰 梁兆刚 潘红鉴 (中国科学院上海天文台，上海200030) 提 要 Jurkevich方法通常用来发现光变曲线的周期，并且根据 极小值的大小来验证所发现的光 变周期可靠性。但是，这种验证方法又带有一定的不确定性，最佳的方法是把所发现的周期与观测 光变曲线拟合，得出每一个周期的振幅，以振幅大小来验证所发现的光变周期的可靠性。将此方法 用于0J287的历史光变曲线，得到其最主要的周期为11．9年，与其它结果相符合。 主题词：光变曲线一周期一Jurkevich方法一振幅检验 分类号：P157．7 1 引 言 发现和研究活动星系核的光变是一项极有意义的工作。较长时标的光变通常是由星系核 的进动或轨道运动产生的，因此可以利用所发现的周期来估计相应的物理参数。 Jurkevich方法 (Jurkevich1971)和Fourier分析都可以用来寻找时间序列的周期。 Jurkevich方法的一个优势就是不要求数据的时间间隔是等间距的，而一些周期分析的方法只 能处理等间隔的数据，这就往往需要通过插值来满足要求，由此带来一些不真实的信息。在实 际过程中大量数据因为各种条件限制不一定是等间隔的时间序列，所以Jurkevich方法在天文 上有较多的应用。Kidger等人 (1992)认为就发现光变曲线周期而言，Jurkevich方法要优于 Fourier分析。然而Jurkevich方法在理论上虽能正确地得出光变曲线的周期，但实际的光变曲 线中除了周期性的变化外还有非周期性的光变和噪声，因此我们在应用Jurkevich方法时发现 它仍有一些不足之处，本文将讨论这些不足之处，并提出一些改进的方法。 本文的讨论将用OJ287的历史光变曲线作为例子。OJ287是有名的活动星系核，它具有 长达一百多年的历史观测资料，其历史光变曲线呈现很强的周期性，并已有很多人对其进行过 研究。我们编辑了OJ287从 1891．04年到2001．98年的历史光变资料，共5085个观测。 Sillanp~ia 等人(1996)指出，由于近来二十多年的观测资料要远远多于以前，如果不对这些资 料作一些处理就直接进行统计的话，必然会使最近的资料的权重远远大于以前的，对统计结果 带来不良影响。一般按时间间隔对这些资料进行平均，得到平均的光变资料后再进行统计分 收稿日期：2004-03—29； 修回日期：2004434．26 上海市科学技术委员会项目(02QA14055)；中国科学院光学天文开放实验室资助课题。 维普资讯 98 中国科学 院上 海天文 台年刊 2004年 析。我们采用0．02年的时间间隔进行平均，得到634个V星等资料。图1是0J287的历史光 变曲线资料。 年份 图1OJ287的历史光变曲线资料(0．02年平均值) Fig．1 HistoriclightcurveofOJ287(theaveragevalueofeach0．02adata) 2 Jurkevich方法 Jurkevich方法是按照不同的试验周期P把数据折叠，然后把折叠好的资料按其相位从小 到大排列，然后分成m组。如第z组的离差平方和 有： ml = ∑ 一m (1) i=l 其中 是时间序列对应的强度值，。为第z组的平均值，