杨海燕，李孝攀，党兴菊

（昭通学院，物理与信息工程学院，云南 昭通 657000）

0 引言

为分析部分活动星系核可能存在的周期性光变，天文学家们提出了一些可行的周期分析算法。这些算法都充分考虑了天文观测数据非均匀、大误差的特点。例如：Jurkevich于1971年提出的一种专门用于分析非等间距数据的期望均方误差统计方法(Jurkevich方法)[1]；Scargle于1981年在Lomb算法[2]基础上提出的著名的Lomb-Scargle(LS)周期图法[3]。Stellingwerf于1978年提出的状态弥散最小化算法(Phase Dispersion Minimization,PDM)算法[4]；，Ferraz-Mello于1981年提出的时间补偿离散傅里叶变换(Date-compensated Discrete Fourier Transform, DCDFT)算 法[5]；Simonetti于1985年提出的结构函数（Structure Functions, SF）法[6]；Edelson 和 Krolik于1988年提出的离散相关函 数 (Discrete Correlation Function, DCF)算 法[7]；Alexander于1997年对DCF算法进行了改进和发展，得到的Z变换离散相关函数(Z-transformed Discrete Correlation Function, ZDCF)算法[8]。Foster于1996年提出的对非平稳信号分析有独特的优势的加权小波Z变换(Weighted Wavelet Z-transform,WWZ)算法[9]。以上方法被广泛应用于活动星系核周期性光变的研究中，对天体辐射机制的研究、辐射区的限定、天体质量的估算等起到了重要的促进作用。

天文观测的局限性会导致观测数据不连续性。Jurkevich方法在不等间距数据处理方面有着独特的优势，被广泛用于活动星系核长周期光变的研究。1998年，樊等[10]使用Jurkevich方法研究了NGC4151自有数据记载以来的光学资料，发现NGC4151的光学波段可能有（14.04±0.8）a的周期性光变；2000年，苏[11]也应用Jurkevich方法对PKS2251+158近100年光变数据进行研究，得到该目标在光学B波段的光变周期约为12.39 a；2001年，樊等[12]再次使用Jurkevich方法研究了3C273自1887年至1997年以来的光学数据，发现光变曲线具有明显的2.0 a光变周期。近年，Li等[13]使用LS方法和Jurkevich方法研究了3C 454.3的多波段光变曲线，发现该目标在红外、光学、软X射线和γ射线波段的光变均存在(210.8±12)d的周期；……