搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于近似熵的斯隆数字化巡天中类星体光变复杂性分析

唐洁 刘晓琴

引用本文:
Citation:

基于近似熵的斯隆数字化巡天中类星体光变复杂性分析

唐洁, 刘晓琴

Analysis on complexity of optical variability based on approximate entropy in Sloan digital sky survey quasars

Tang Jie, Liu Xiao-Qin
PDF
HTML
导出引用
  • 光变是类星体的重要观测特征之一, 类星体在多个波段存在剧烈的光变现象. 光变非常复杂, 具有非线性特征. 以斯隆数字化巡天(Sloan digital sky survey, SDSS) stripe 82天区中的类星体为研究对象, 利用近似熵方法分析了类星体光变的复杂性. 首先应用模拟信号检验了近似熵方法对周期序列、白噪声序列、混沌序列和组合序列的区分能力, 验证了近似熵方法是一种识别不同类型时间序列的有效方法. 再计算了SDSS第7次释放数据中光谱证认过的类星体光变的近似熵, 并分析了它们的复杂性. 结果表明: SDSS类星体光变的近似熵值最大值为0.58, 类星体光变的复杂性介于周期序列和白噪声序列的复杂性之间, 近一半样本的复杂性与混沌序列基本一致.
    Variability is one of the most important observational features of quasars, and it is still not clear that the different quasars show different characteristic variability patterns. The optical variability of quasar is very complex, and optical variability has the non-linear characteristic of complex system. In this paper, the approximate entropy method is employed to analyze the complexities of variability in the Sloan digital sky survey (SDSS) stripe 82 quasars. Firstly, in order to show that the approximate entropy method has the effective ability to distinguish the different types of time sequences, the approximate entropy of periodic sequence, noise sequence, chaotic sequence and their mixed sequences are calculated by using the analog signals. The approximate entropy method proves to be an effective method to identify different types of time sequences. Then, we present the approximate entropy of optical variability of spectroscopically-confirmed quasars from the SDSS data release 7 spectroscopic catalog, and their complexities are analyzed. The results show that the maximum approximate entropy of quasars’ optical variability is only 0.58. The complexity of quasars’ optical variability is between the complexities of periodic sequence and white noise sequence. For nearly half of the samples, the complexities of their optical variability are basically consistent with the complexity of chaotic sequence. Quasars’ optical variability is neither completely periodic nor completely stochastic.
      通信作者: 唐洁, tj168@163.com
    • 基金项目: 国家自然科学基金(批准号: 11373008)、陕西省自然科学基金(批准号: 2016JM1028)和陕西省教育厅科研计划(批准号: 14JK1157)资助的课题.
      Corresponding author: Tang Jie, tj168@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11373008), the Program of the Natural Science Foundation of Shaanxi Province, China (Grant No. 2016JM1028), and the Scientific Research Program Funded by Shaanxi Provincial Education Department, China (Grant No. 14JK1157).
    [1]

    Fan J H, Zhang Y W, Qian B C, Tao J, Liu Y, Hua T X 2009 Astrophys. J. Suppl. Ser. 181 466Google Scholar

    [2]

    黄克谅 2005 类星体与活动星系核 (北京: 科学出版社) 第3页

    Huang K L 2005 Quasars and Active Galactic Nucleus (Beijing: Science Press) p3 (in Chinese)

    [3]

    Kawaguchi T, Mineshige S, Umemura M 2002 Astrophys. J. 504 671

    [4]

    MacLeod C L, ŽIvezić C S, Kochanek S, Kozłowski B, Kelly E, Bullock A, Kimball B, Sesar D, Westman K B 2010 Astrophys. J. 721 1014Google Scholar

    [5]

    Kasliwal P V, Vogeley S M, Richards T G 2015 Mon. Not. R. Astron. Soc. 451 4328Google Scholar

    [6]

    Guo H X, Wang J X, Cai Z Y, Sun M Y 2017 Astrophys. J. 847 132Google Scholar

    [7]

    Harko T, Leung C S, Mocanu G 2014 European Physical Journal C 74 2900Google Scholar

    [8]

    Charisi M, Bartos I, Haiman Z, Price-Whelan A M, Márka S 2015 Mon. Not. R. Astron. Soc. 45 4

    [9]

    D'Orazio D J, Haiman Z, Schiminovich D 2015 Nature 525 351Google Scholar

    [10]

    Bhatta G, Zola S, Stawarz Ł, Ostrowski M, Winiarski M, Ogłoza W, Dróżdż M, Siwak M, Liakos A, et al 2016 Astrophys. J. 832 47Google Scholar

    [11]

    Carrasco L, Dultzinhacyan D, Cruzgonzalez I 1985 Nature 314 146Google Scholar

    [12]

    Kidger M, Takalo L, Sillanpää A 1992 Astron. Asrtophys. 264 32

    [13]

    Valtonen M J, Lehto H J, Sillanpää A, Nilsson K, Mikkola S, Hudec R, Basta M, Teräsranta H, Haque S, Rampadarath H 2006 Astrophys. J. 646 34

    [14]

    Misra R, Harikrishnan K P, Ambika G, Kembhavia A K 2006 Advances in Space Research 643 1114

    [15]

    Harikrishnan K P, Misra R, Ambika G 2011 Research in Astron. Astrophys. 11 71Google Scholar

    [16]

    Emmanoulopoulos D 2007 Ph. D. Dissertation (Heidelber: Ruperto-Carola University of Heidelberg)

    [17]

    唐洁 2014 物理学报 63 049701Google Scholar

    Tang J 2014 Acta. Phys. Sin. 63 049701Google Scholar

    [18]

    郝柏林 2001 物理 30 466Google Scholar

    Hao B L 2001 Physics 30 466Google Scholar

    [19]

    Pincus S M, Huang W M 1992 Commun. Statist-Theory Meth. 21 3061Google Scholar

    [20]

    金宁德, 董芳, 赵舒 2007 物理学报 56 0720

    Jin N D, Dong F, Zhao S 2007 Acta. Phys. Sin. 56 0720

    [21]

    Lempel A, Ziv J 1976 IEEE Trans. on Inf. Theory 22 75Google Scholar

    [22]

    郭家梁, 钟宁, 马小萌, 张明辉, 周海燕 2016 物理学报 65 190501Google Scholar

    Guo J L, Zhong N, Ma X M, Zhang M H, Zhou H Y 2016 Acta. Phys. Sin. 65 190501Google Scholar

    [23]

    Pettitt A N, Stephens M A 1977 Technometrics 19 205Google Scholar

    [24]

    Hughes P A, Aller H D, Aller M F 1992 Astrophys. J. 396 469Google Scholar

    [25]

    Scargle J D 1982 Astrophys. J. 263 835Google Scholar

    [26]

    唐洁, 刘晓琴, 傅明星 2018 陕西理工大学学报(自然科学版) 33 87Google Scholar

    Tang J, Liu X Q, Fu M X 2018 Journal of Shaanxi University of Technology (Natural Science Edition) 33 87Google Scholar

  • 图 1  典型时间序列

    Fig. 1.  Typical time series

    图 2  典型时间序列的近似熵

    Fig. 2.  The approximate entropy of typical time series

    图 3  SDSS J012228.72-001332.0的光变曲线

    Fig. 3.  Light curve of SDSS J012228.72-001332.0

    图 4  SDSS J012228.72-001332.0的近似熵

    Fig. 4.  The approximate entropy of SDSS J012228.72-001332.0

    图 5  SDSS类星体的近似熵分布

    Fig. 5.  Approximate entropy distribution of SDSS quasars

    表 1  SDSS类星体的近似熵

    Table 1.  Approximate entropy of SDSS quasars

    近似熵值 0—0.1 0.1—0.2 0.2—0.3 0.3—0.4 0.4—0.5 0.5—0.6
    波段 总样本数 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/%
    i 6465 39 0.60 365 5.65 1848 28.58 3106 48.04 1077 16.66 30 0.47
    r 6547 40 0. 61 438 6.69 1898 28.99 3131 47.82 1003 15.32 37 0.57
    g 6439 34 0.53 374 5.81 1794 27.86 3153 48.97 1033 16.04 51 0.79
    z 6373 40 0.63 203 3.19 1438 22.56 3223 50.57 1400 21.97 69 1.08
    u 6092 27 0.44 190 3.12 1503 24.67 3117 51.17 1214 19.93 41 0.67
    下载: 导出CSV
  • [1]

    Fan J H, Zhang Y W, Qian B C, Tao J, Liu Y, Hua T X 2009 Astrophys. J. Suppl. Ser. 181 466Google Scholar

    [2]

    黄克谅 2005 类星体与活动星系核 (北京: 科学出版社) 第3页

    Huang K L 2005 Quasars and Active Galactic Nucleus (Beijing: Science Press) p3 (in Chinese)

    [3]

    Kawaguchi T, Mineshige S, Umemura M 2002 Astrophys. J. 504 671

    [4]

    MacLeod C L, ŽIvezić C S, Kochanek S, Kozłowski B, Kelly E, Bullock A, Kimball B, Sesar D, Westman K B 2010 Astrophys. J. 721 1014Google Scholar

    [5]

    Kasliwal P V, Vogeley S M, Richards T G 2015 Mon. Not. R. Astron. Soc. 451 4328Google Scholar

    [6]

    Guo H X, Wang J X, Cai Z Y, Sun M Y 2017 Astrophys. J. 847 132Google Scholar

    [7]

    Harko T, Leung C S, Mocanu G 2014 European Physical Journal C 74 2900Google Scholar

    [8]

    Charisi M, Bartos I, Haiman Z, Price-Whelan A M, Márka S 2015 Mon. Not. R. Astron. Soc. 45 4

    [9]

    D'Orazio D J, Haiman Z, Schiminovich D 2015 Nature 525 351Google Scholar

    [10]

    Bhatta G, Zola S, Stawarz Ł, Ostrowski M, Winiarski M, Ogłoza W, Dróżdż M, Siwak M, Liakos A, et al 2016 Astrophys. J. 832 47Google Scholar

    [11]

    Carrasco L, Dultzinhacyan D, Cruzgonzalez I 1985 Nature 314 146Google Scholar

    [12]

    Kidger M, Takalo L, Sillanpää A 1992 Astron. Asrtophys. 264 32

    [13]

    Valtonen M J, Lehto H J, Sillanpää A, Nilsson K, Mikkola S, Hudec R, Basta M, Teräsranta H, Haque S, Rampadarath H 2006 Astrophys. J. 646 34

    [14]

    Misra R, Harikrishnan K P, Ambika G, Kembhavia A K 2006 Advances in Space Research 643 1114

    [15]

    Harikrishnan K P, Misra R, Ambika G 2011 Research in Astron. Astrophys. 11 71Google Scholar

    [16]

    Emmanoulopoulos D 2007 Ph. D. Dissertation (Heidelber: Ruperto-Carola University of Heidelberg)

    [17]

    唐洁 2014 物理学报 63 049701Google Scholar

    Tang J 2014 Acta. Phys. Sin. 63 049701Google Scholar

    [18]

    郝柏林 2001 物理 30 466Google Scholar

    Hao B L 2001 Physics 30 466Google Scholar

    [19]

    Pincus S M, Huang W M 1992 Commun. Statist-Theory Meth. 21 3061Google Scholar

    [20]

    金宁德, 董芳, 赵舒 2007 物理学报 56 0720

    Jin N D, Dong F, Zhao S 2007 Acta. Phys. Sin. 56 0720

    [21]

    Lempel A, Ziv J 1976 IEEE Trans. on Inf. Theory 22 75Google Scholar

    [22]

    郭家梁, 钟宁, 马小萌, 张明辉, 周海燕 2016 物理学报 65 190501Google Scholar

    Guo J L, Zhong N, Ma X M, Zhang M H, Zhou H Y 2016 Acta. Phys. Sin. 65 190501Google Scholar

    [23]

    Pettitt A N, Stephens M A 1977 Technometrics 19 205Google Scholar

    [24]

    Hughes P A, Aller H D, Aller M F 1992 Astrophys. J. 396 469Google Scholar

    [25]

    Scargle J D 1982 Astrophys. J. 263 835Google Scholar

    [26]

    唐洁, 刘晓琴, 傅明星 2018 陕西理工大学学报(自然科学版) 33 87Google Scholar

    Tang J, Liu X Q, Fu M X 2018 Journal of Shaanxi University of Technology (Natural Science Edition) 33 87Google Scholar

  • [1] 苟竞, 刘俊勇, 魏震波, Gareth Taylor, 刘友波. 基于多尺度熵的电力能量流复杂性分析. 物理学报, 2014, 63(20): 208402. doi: 10.7498/aps.63.208402
    [2] 周杰, 王亚林, 菊池久和. 多天线信道空间衰落相关性近似算法及其复杂性研究. 物理学报, 2014, 63(23): 230205. doi: 10.7498/aps.63.230205
    [3] 向郑涛, 陈宇峰, 李昱瑾, 熊励. 基于多尺度熵的交通流复杂性分析. 物理学报, 2014, 63(3): 038903. doi: 10.7498/aps.63.038903
    [4] 欧建文, 郑永刚, 张雄. 振荡吸积盘的光变研究. 物理学报, 2014, 63(23): 239801. doi: 10.7498/aps.63.239801
    [5] 唐洁. 基于集合经验模态分解的类星体光变周期及其混沌特性分析. 物理学报, 2014, 63(4): 049701. doi: 10.7498/aps.63.049701
    [6] 孙克辉, 贺少波, 何毅, 尹林子. 混沌伪随机序列的谱熵复杂性分析. 物理学报, 2013, 62(1): 010501. doi: 10.7498/aps.62.010501
    [7] 张学清, 梁军. 基于EEMD-近似熵和储备池的风电功率混沌时间序列预测模型. 物理学报, 2013, 62(5): 050505. doi: 10.7498/aps.62.050505
    [8] 唐洁, 傅明星, 吴学兵. 基于结构函数方法的类星体证认. 物理学报, 2012, 61(21): 219501. doi: 10.7498/aps.61.219501
    [9] 何文平, 何涛, 成海英, 张文, 吴琼. 基于近似熵的突变检测新方法. 物理学报, 2011, 60(4): 049202. doi: 10.7498/aps.60.049202
    [10] 张伟超, 杨立军, 吕小青. 基于近似熵测度的铝合金P-MIG焊亚射流过渡自适应控制研究. 物理学报, 2011, 60(2): 020601. doi: 10.7498/aps.60.020601
    [11] 董富通, 张雄. 类星体3C345的光变周期特性. 物理学报, 2009, 58(11): 8116-8122. doi: 10.7498/aps.58.8116
    [12] 何文平, 王启光, 吴琼, 张文, 张勇. 滑动去趋势波动分析与近似熵在动力学结构突变检测中的性能比较. 物理学报, 2009, 58(4): 2862-2871. doi: 10.7498/aps.58.2862
    [13] 周磊, 龚志强, 支蓉, 封国林. 基于复杂网络研究中国温度变化的区域特征. 物理学报, 2009, 58(10): 7351-7358. doi: 10.7498/aps.58.7351
    [14] 王启光, 张增平. 近似熵检测气候突变的研究. 物理学报, 2008, 57(3): 1976-1983. doi: 10.7498/aps.57.1976
    [15] 庄建军, 宁新宝, 邹 鸣, 孙 飙, 杨 希. 两种熵测度在量化射击运动员短时心率变异性信号复杂度上的一致性. 物理学报, 2008, 57(5): 2805-2811. doi: 10.7498/aps.57.2805
    [16] 李 强, 王太勇, 冷永刚, 何改云, 何慧龙. 基于近似熵测度的自适应随机共振研究. 物理学报, 2007, 56(12): 6803-6808. doi: 10.7498/aps.56.6803
    [17] 金宁德, 董 芳, 赵 舒. 气液两相流电导波动信号复杂性测度分析及其流型表征. 物理学报, 2007, 56(2): 720-729. doi: 10.7498/aps.56.720
    [18] 曹 彪, 吕小青, 曾 敏, 王振民, 黄石生. 短路过渡电弧焊电流信号的近似熵分析. 物理学报, 2006, 55(4): 1696-1705. doi: 10.7498/aps.55.1696
    [19] 谢勇, 徐健学, 杨红军, 胡三觉. 皮层脑电时间序列的相空间重构及非线性特征量的提取. 物理学报, 2002, 51(2): 205-214. doi: 10.7498/aps.51.205
    [20] 尤峻汉, 程富华. 佘楞科夫线状发射和类星体光谱. 物理学报, 1980, 29(7): 927-936. doi: 10.7498/aps.29.927
计量
  • 文章访问数:  3040
  • PDF下载量:  26
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-22
  • 修回日期:  2019-05-13
  • 刊出日期:  2019-07-20

基于近似熵的斯隆数字化巡天中类星体光变复杂性分析

  • 1. 陕西理工大学物理与电信工程学院, 汉中 723001
  • 2. 汉中职业技术学院药学与医学技术系, 汉中 723002
  • 通信作者: 唐洁, tj168@163.com
    基金项目: 国家自然科学基金(批准号: 11373008)、陕西省自然科学基金(批准号: 2016JM1028)和陕西省教育厅科研计划(批准号: 14JK1157)资助的课题.

摘要: 光变是类星体的重要观测特征之一, 类星体在多个波段存在剧烈的光变现象. 光变非常复杂, 具有非线性特征. 以斯隆数字化巡天(Sloan digital sky survey, SDSS) stripe 82天区中的类星体为研究对象, 利用近似熵方法分析了类星体光变的复杂性. 首先应用模拟信号检验了近似熵方法对周期序列、白噪声序列、混沌序列和组合序列的区分能力, 验证了近似熵方法是一种识别不同类型时间序列的有效方法. 再计算了SDSS第7次释放数据中光谱证认过的类星体光变的近似熵, 并分析了它们的复杂性. 结果表明: SDSS类星体光变的近似熵值最大值为0.58, 类星体光变的复杂性介于周期序列和白噪声序列的复杂性之间, 近一半样本的复杂性与混沌序列基本一致.

English Abstract

    • 作为活动性最强活动星系核, 类星体除了强发射线、较大红移、高光度等观测特征外, 还具有快速、不规则、剧烈的光变现象[1,2], 从长波的射电波段到高能的伽玛波段都观测到较强的光变信号. 对光变的研究可以为我们提供许多重要信息, 有利于更深入地了解光变产生的物理机制. 获得了很多观测现象支持的吸积盘模型认为[3], 吸积盘结构的改变, 吸积盘自身的不稳定性和变化的吸积率都有可能导致类星体光变现象的发生. 类星体的光变资料主要来源于长时间的多波段测光监测, 目前已经积累了大量的观测资料, 为进一步研究类星体光变提供了便利条件.

      类星体被发现后不久, 其辐射流量的变化就被观测到了. 已有研究对类星体光变现象从不同角度做了大量的分析. MacLeod等[4]认为类星体光变可以利用随机游走模型来拟合, 但Kasliwal等[5]却认为该模型并非适合所有的类星体, Guo等[6]认为随机游走模型不能极好的拟合部分类星体光线曲线. 解释光变现象发生常用的吸积盘不稳定模型也是将光变视为随机事件, 盘通过黏滞耗散和随机因素同外部介质产生作用[3,7]. 随着观测技术水平的提高和一些大型望远镜的使用, 在对一些长期监测的类星体光变资料分析时发现光变呈显出准周期现象[8,9], 如OJ 287长期光变曲线存在约12年的长周期[10,11], 运用该周期推测OJ 287下一次可能爆发的时间, 其后观测到的爆发时间与该预测时间非常吻合[12,13]. 随着非线性动力学理论的发展, 非线性动力学分析方法被引入到类星体光变的研究中来[1416]. 如果光变具有混沌特性, 非线性机制就可以被用来解释类星体光变.

      类星体光变资料里隐含有线性和非线性的成分, 是极其复杂的[17]. 将不规则的光变现象视为复杂的, 仅是一种假设, 简单的定性分析缺乏说服力. 类星体光变是不是属于复杂系统? 若其存在复杂性, 复杂到何种程度? 这些问题值得进一步深入研究. 复杂性研究处于有序与混沌的边缘, 现在主要指的是对非线性现象的研究[18]. 混沌动力学理论可以进行复杂性分析, 但混沌分析方法一般要求数据量非常大, 少则上千, 多则上万, 对观测资料的质量要求也较高, 要想获得较完整的类星体光变的资料需要占用望远镜较长的观测时间. 尽管与过去相比, 目前类星体的监测项目增多, 但监测样本数还是较少, 观测时间的跨度也较短.

      斯隆数字化巡天(Sloan digital sky survey, SDSS)项目提供了大量类星体测光监测数据, SDSS第7次释放数据(data release 7, DR7)释放有光谱证认过的类星体5个波段10年观测的数据, 但单个类星体样本光变资料数据量少. 常用的、对时间序列长度要求相对较低的复杂性分析方法有近似熵[19,20]、Lempel-Ziv算法[21]、样本熵[22]和Kolmogrov熵[23]等方法. 超过100个数据量时, 利用近似熵方法就能获得比较可靠的估计值, 并且没有粗粒化处理要求, 这样可以保留原始数据的有效信息[19], 近似熵方法优于后两种方法. 本文利用近似熵方法计算和度量类星体光变的复杂性, 分析类星体光变的复杂度有助于我们更好地了解光变的复杂特征, 揭示复杂光变的非线性动力学特性.

    • 混沌动力学理论常用Kolmogrov熵来度量动力学系统的复杂性, 要获得Kolmogrov熵值, 理论上要求被分析的对象数据量足够多, 且该方法抗噪能力差, 对于信噪比低的时间序列较难获得可靠的熵值, 因此最好别叠加有噪声, 但实测数据很难达到这些要求. Pincus和Huang[19]提出了对数据量长度要求相对较低, 有较好的抗干扰、抗噪能力的近似熵(approximate entropy, ApEn). 不同于Kolmogrov熵, 近似熵不需要利用大量数据重构奇异吸引子, 主要是从统计角度来度量时间序列产生新模式的概率大小, 概率越大复杂性越大, 近似熵值也越大, 随机成分增多; 近似熵值越小, 时间序列周期成分增多.

      若存在一个给定长度的一维时间序列$X{\rm{ = }}$$\left\{ {x\left( n \right), n = 1, 2, \cdots , N} \right\}$, 如模式维数m和相似容限r已经确定, 构建一个m维的向量空间:

      $\begin{split} &{{X}}\left( n \right) = \left[ {x\left( n \right), x\left( {n + 1} \right), \cdots , x\left( {n + m - 1} \right)} \right],\\ & n = 1, 2, \cdots , N - m + 1. \end{split}$

      将任意向量${{X}}\left( i \right)$与其他向量${{X}}\left( j \right)$的对应元素差值最大的记为

      $\begin{split} &{D_{ij}} = \max \left| {x\left( {i + k} \right) - x\left( {j + k} \right)} \right|, \;\;\\ &\text{且}\;i \ne j, k = 0, 1, 2, \cdots , m - 1. \end{split}$

      对每个向量${{X}}\left( i \right)$计算${D_{ij}}$小于相似容限r的数目, 将该数目与向量总数$N - m$的比记为$A_i^m\left( r \right)$. 将模式维数增加1, 按上面的方法计算维数$m + 1$$A_i^{m + 1}\left( r \right)$, 时间序列的近似熵Sapen可表示为[19]

      $ {S_{{\rm{apen}}}}\left( {m, r} \right) = - \frac{1}{{N - m}}\sum\limits_{i = 1}^{N - m} {\ln \frac{{A_i^{m + 1}\left( r \right)}}{{A_i^m\left( r \right)}}} . $

    • 为了检验近似熵对于不同时间序列的区分能力, 选取周期序列、混沌序列、白噪声序列和组合序列, 分析5种典型时间序列的近似熵, 通过分析典型不同类型时间序列在复杂度的取值范围, 可以作为类星体光变属于哪种类型时间序列的参照.

      5种典型时间序列如图1所示, 其中周期序列(用Sine表示)为$y = \sin \left( {2{\text{π}} t/10} \right)$; 混沌序列(用Logistic表示)为$x\left( {n + 1} \right) = \omega x\left( n \right)\left[ {1 - x\left( n \right)} \right]$, 其中初始值$x\left( 0 \right) = 0.7, \;\;\omega = 0.75$, 由该方程产生混沌序列; 白噪声序列 (用White noise表示)由均值为0、方差为1的随机序列组成; 组合序列包括混沌序列叠加周期序列(用Logistic+Sine表示)和混沌序列叠加白噪声序列(用Logistic+White noise表示).

      图  1  典型时间序列

      Figure 1.  Typical time series

      对于时间序列长度的选取, Pincus和Huang[19]通过他们的大量研究, 建议数据点数最好取100以上, 一般也不要超过5000. 我们将时间序列长度都取900, 模式维数m = 2, 相似容限r取0.05—0.50, 5种典型时间序列的近似熵分析结果见图2. 从近似熵值分布可以看出: 周期序列的近似熵最小, 白噪声序列的近似熵最大, 混沌序列的近似熵值介于周期序列和白噪声序列之间, 混沌序列叠加周期序列比混沌序列的近似熵值略大, 相似容限小于0.35时, 混沌序列叠加白噪声序列比白噪声序列的近似熵值小; 相似容限大于0.35时, 混沌序列叠加白噪声序列比白噪声序列的大. 对于不同类型的时间序列, 近似熵值按周期序列、混沌序列、混沌序列叠加周期序列、混沌序列叠加白噪声序列到白噪声序列的顺序逐渐增大, 近似熵值越大复杂性也越高, 周期序列、混沌序列、白噪声序列的近似熵值差别很大, 各自分别在不同的区间, 白噪声序列是混沌序列的近似熵值近5倍, 白噪声序列和混沌序列叠加白噪声序列复杂度差别不大, 白噪声序列近似熵值接近组合序列, 近似熵值都超过1.

      图  2  典型时间序列的近似熵

      Figure 2.  The approximate entropy of typical time series

      这些结果说明近似熵能有效将混沌序列、周期序列、白噪声序列和不同成分组合序列区分开, 复杂性能够区分不同随机程度的时间序列, 验证了近似熵可以描述时间序列随机程度, 对混沌有一定的识别能力, 较好地表征了非线性结构的复杂性, 且复杂度随着随机成分增加而增大.

    • SDSS项目使用2.5 m专用望远镜, 装备了先进的大视场扫描成像装置, 观测底板的640条光纤可以同时拍摄640个目标, 从2000年项目开始巡天观测后, 已经积累了许多高质量的u, g, r, i, z波段测光数据.

      为了展示近似熵方法在度量SDSS类星体光变复杂性的能力, 我们先以类星体SDSS J012228.72-001332.0为例来分析该方法的可行性. 由于SDSS有很多观测任务, 望远镜观测时间有限, 不可能长时间观测某个类星体, 这些类星体大多是间隔90—270 d才有测光数据. 为了便于计算近似熵, 缺失值用插值方法求得, 图3给出了类星体SDSS J012228.72-001332.0原始和插值处理后的5个波段光变曲线.

      图  3  SDSS J012228.72-001332.0的光变曲线

      Figure 3.  Light curve of SDSS J012228.72-001332.0

      在计算近似熵的值前, 需要先选取模式维数m、相似容限r和时间序列的长度N. 对于N的选取, Pincus和Huang[19]通过他们的大量研究, 建议数据点数最好取100以上, 一般也不要超过5000, 本文对原观测数据每隔10 d插值一个数据, 数据量完全满足条件. 对于模式维数m, 一般选取m = 2, 相对于m = 1, 可以在时间序列的联合概率重构过程中, 获得更多更详细的信息; 若取m大于2, 相似容限r要求比较大, 但这样会导致丢失许多信息. 对于参数相似容限r, Pincus和Huang[19]认为相似容限r取0.10—0.25倍原始数据的标准差, 可以估计出比较理想的统计特性.

      图4给出了类星体SDSS J012228.72-001332.0光变资料的相似容限-近似熵图, 在相似容限为0.2前, 近似熵值随相似容限的增加而增大, 相似容限r在增为0.2倍原始数据的标准差后, 5个波段光变资料的近似熵值基本保持一个比较稳定的值. g波段近似熵最小, 在0.3左右振荡; u波段近似熵最大, 最大值接近0.35, 相似容限在区间0.4—0.6时, 近似熵值随相似容限增加有减小的趋势. 因此, 在本文后面的大样本近似熵值分析时, 将相似容限r取为0.2倍原始数据的标准差. 从单个样本近似熵分析可以看出, 5个波段光变资料的近似熵值相差不多, 同一类星体不同波段光变曲线复杂性基本相同, 这可能是因为它们的光变曲线波形非常相似.

      图  4  SDSS J012228.72-001332.0的近似熵

      Figure 4.  The approximate entropy of SDSS J012228.72-001332.0

      本文选取的类星体样本来源于SDSS DR7中的stripe 82天区光谱证认过的类星体, Macleod等[4]搜集整理提供了9258个类星体测光数据, 为了获得可靠的分析结果, 要求每隔360 d需要至少一次的观测数据, 由于2000—2003年观测数据太少, 起始时间选择从2003年开始, 并除去测光误差大于0.1个星等的数据点, 并且要求这些数据大致均匀分布, 按这些条件再选择观测数据点数超过20的作为样本. 最后, 我们得到i, r, g, z, u波段的样本数分别是6465, 6547, 6439, 6373和6092.

      将模式维数和相似容限分别取m = 1, r = 0.2对选取的样本进行近似熵分析, 获得的结果如图5表1所示. 可以看出, 类星体光变的近似熵值都处于0—0.6范围内, 最大值分别是0.552, 0.553, 0.575, 0.571和0.580, 最小值都为0.008, 接近于0. 白噪声序列和混有白噪声序列的混沌序列都比光变的近似熵最大值超过2倍, 因此类星体样本光变不可能完全是白噪声序列, 这是因为类星体观测总会有噪声干扰, 但其占的比例不大. 近一半的近似熵值分布在0.3—0.4的范围内, 远大于周期序列, 又远小于白噪声序列, 和混沌序列在初始值取为$x\left( 0 \right) = 0.7, \;\;\omega = 0.75$时的近似熵值非常接近, 说明近一半的类星体光变既不是完全周期性的也不是完全白噪声的系统, 而是一个与混沌序列近似的、带有非线性结构的、具有混沌特性的动力学系统.

      近似熵值 0—0.1 0.1—0.2 0.2—0.3 0.3—0.4 0.4—0.5 0.5—0.6
      波段 总样本数 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/% 样本数 百分数/%
      i 6465 39 0.60 365 5.65 1848 28.58 3106 48.04 1077 16.66 30 0.47
      r 6547 40 0. 61 438 6.69 1898 28.99 3131 47.82 1003 15.32 37 0.57
      g 6439 34 0.53 374 5.81 1794 27.86 3153 48.97 1033 16.04 51 0.79
      z 6373 40 0.63 203 3.19 1438 22.56 3223 50.57 1400 21.97 69 1.08
      u 6092 27 0.44 190 3.12 1503 24.67 3117 51.17 1214 19.93 41 0.67

      表 1  SDSS类星体的近似熵

      Table 1.  Approximate entropy of SDSS quasars

      图  5  SDSS类星体的近似熵分布

      Figure 5.  Approximate entropy distribution of SDSS quasars

      近似熵值的大小和被分析的时间序列复杂程度存在成正比例的关系, 时间序列越复杂, 其近似熵值也就越大, 时间序列越趋近于随机性, 反之, 近似熵值越小, 越趋于周期性. 类星体光变的近似熵值低于0.2的样本数很少, 最多的r波段才占到7.3%, 最少的u波段只有3.56%, 说明周期性光变占主导的类星体比例很少. 近似熵值高于0.5的类星体样本数也少, 除z波段刚超过1%, 其余波段1%都不到, 说明所有类星体光变的近似熵值都偏小, 类星体光变趋于随机性概率较少. 同一类星体样本不同波段, 近似熵值都比较接近, 说明同一类星体光变的复杂性基本一致. 不同类星体近似熵值不尽相同, 说明类星体光变的非线性结构是有区别的, 可能光变产生的原因不一样引起的.

    • 自从类星体发现以来, 剧烈的光变现象引起了天文学家的广泛关注. 尽管现在观测技术有很大进步, 但现有的观测设备仍然无法分辨类星体的中心区域, 对类星体光变的深入研究有助于了解其中心区域辐射的物理机制. 常用分析类星体光变资料的方法如结构函数[24]、功率谱[25]、经验模态分解[26]等, 应用这些方法如要获得可信度高的结论, 对观测数据量要求较高, 这就对数据分析方法提出了要求. 本文采用的近似熵只需要100个以上的数据点就能获得比较可靠的复杂度, 完全满足数据分析的需要.

      本文利用近似熵方法分析了混沌序列、周期序列、白噪声序列和不同成分组合序列的近似熵, 发现近似熵值能有效区分不同随机程度的时间序列, 验证了近似熵对不同复杂性的时间序列有较好的识别能力, 金宁德等[20]也用模拟信号获得了跟我们一样的结论. 本文进一步基于SDSS stripe 82天区的类星体测光数据分析了类星体光变的复杂性, 从单个源的分析结果看来, 5个波段光变曲线波形非常相似, 它们的近似熵值相差较小, 同一类星体不同波段光变曲线复杂性基本相同. 对6000多个大样本的类星体分析发现, 所有的近似熵值都低于0.6, 最大值为0.58, 说明类星体光变中白噪声占的比例都不高. 近一半的类星体光变和Logistic方程产生的混沌序列的复杂性是一样的, 说明一部分类星体光变以混沌成分占主导. 小部分近似熵值接近于0, 这说明周期性成分占主导的类星体比例非常少. 类星体光变不可能是完全周期性的, 也不可能是白噪声, 因为若是完全周期性或白噪声, 近似熵值是0或超过1, 也不可能是完全混沌的, 因为观测中总夹杂着噪声的影响. 因此类星体光变可能是周期性成分, 混沌成分和噪声叠加在一起, 并以3种成分中的一种占主导.

      用结构函数、功率谱等周期分析方法分析类星体长期光变资料, 发现光变具有周期性[1,813], 这和我们的结论并不矛盾, 这些文献只提取了类星体光变中隐含的部分周期信息, 忽略了其他部分的有效信息. MacLeod等[4]用随机游走模型来描述光变, 发现该模型能较好地拟合类星体光变曲线, 但Kasliwal等[5]和Guo等[6]认为随机游走模型在描述类星体光变中还是存在部分缺陷. Misra等[14]运用混沌动力学理论分析了GRS 1915+105 的光变曲线, 获得的结果是GRS 1915+105光变也存在混沌性叠加随机性的现象. 唐洁[17]应用集合经验模态分解方法将类星体光变资料分解成周期项、混沌项和趋势项, 认为类星体光变是由周期成分、混沌成分和趋势成分叠加而成. 这些已有的研究成果与我们分析的结论是比较一致的.

      感谢美国SDSS项目提供的类星体光变数据.

参考文献 (26)

目录

    /

    返回文章
    返回