用WWZ加权小波分析方法计算S5 0716+714的周期与研究.docVIP

用WWZ加权小波分析方法计算S5 0716+714的周期研究 摘要:本文采用的研究方法是加权小波Z变换Weighted wavelet Z-transform ，WWZ)，它是一种基小波分析原理等间隔数据的方法收集BL Lac天体S5 0716+714B、V、R、I波段较为完备的观测数据，了其一定时期内光变曲线。WWZ对S5 0716+714以上四个波段光数据进行分析发现5 0716+714存在的光变周期为 关键词: 小波分析；WWZ变换BL Lac 天体 引言 星系核是现代天体物理学中的一个热门研究领域星系核简称AGNs(Active Galactic Nucleus)星系核的中有关天体物理的基础性问题如能量的、辐射机制等是活动星系核的一个重要中发现具有许多的极端特性BL Lac天体又是类星体中的一个重要代表其观测上表现出偏振光变以及热辐射等特征1]。S5 0716+714是BL Lac天体之一，备受各界科学研究者的关注。对于活动星系而言光变周期的研究是了解其重要特性和本质的一个非常重要的方法。通过对光变周期等特性的研究2]。周期光变是类星体的一个普遍的特征在几年来也有不少研究组提出了多种研究类星体长周期光变的方法如解释长周期光变的理论模型双黑洞喷流模型进动吸积盘模型等 用来寻找计算天体周期光变的方法有很多经常用到的有时间序列的功率分析方法3]，Jurkevich方法4]，相位分析法自相关函数分析法等5]。虽然以上这些方法在研究天体周期光变中比较广泛但是也存在许多的条件限制和局限性在天体光周期的研究中将会有所困难。因此，本文中将使用一种较为新颖的方法加权小波Z变换对BL Lac天体S5 0716+714的光变周期进行计算和研究。 WWZ是一种基于小波分析原理产生的能够用于非规则数据处理的新方法[6]。析方法是一种非平稳信号分析方法它和功率谱分析方法都为傅立叶变换的周期分析方法。但是小波分析方法是将信号表示为时间尺度-频率尺度的域解析。因此在整个频率域中，具有多频率分辨的特点。由于小波分析的主要核心部分是多分辨率分析因此较为适合分析非平稳信号所以用WWZ来分析非平稳信号是一种较好的方法。 本文将采用WWZ变换对S5 0716+714的个光学波段周期数据进行计算研究最后得出S5 0716+714的光变周期的部分变化特性。 加权小波Z变换WWZ) 2.1 小波分析原理 传统的信号理论，是建立在傅立叶分析基础上的，而傅立叶变换作为一种全局性的变化，其有一定的局限性。小波变换与傅立叶变换相比，是一个时间和频域的局域变换因而能有效地从信号中提取信息，通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析（Multiscale Analysis），解决了傅立叶变换不能解决的许多困难问题。非常适合用于探测夹带瞬间反常现象的信号。 小波变换的定义7]为 设，则它的连续小波变换为： （1） 其中a为伸缩尺度，b为平移参数。 将其离散化形式为： （2） 称为小波函数，它是由小波母函数经过平移b和伸缩a得到。是一种长度有限，平均值为0的波形。它具有快速衰减的特性。通常所用的小波母函数为Morlet小波，它的形式为8] ： （3） 是衰减因子，当取较大值时上式第二项可以省略，因此有简化的Morlet小波： （4） Torrence et al 指出对于Morlet小波伸缩尺度a和傅立叶变换的频率之间有关系[9]： （5） 2.2 WWZ原理 分析方法能够较好的处理非平稳信号，但是在处理不规则的非等间距数据，有一的缺陷。因为（）和（）近似，在实际的应用中一定会出现一定的偏差，尤其是在处理非等间距数据的时候。，小波变换的结果受到边界效应的影响。又对于天文观测信号，会受到观测季节、天气、月相等因素的影响，观测到的数据往往是非等间距的，间距还会很大，间距还可能周期性的出现。这些，不论是对于傅变换还是小波变换，都会给分析研究带来相当大的困难。处理间距数据问题中，在使用傅叶或者小波变换，我们通常采用的方是插值方法，而对于插值法所到的数据是否真实可靠又是一个很大的影响。 处理非等间距数据，使用傅立叶处理时变换中会出现周期。不少的研究者在研究过程中提出了的方法处理其中的周期问题。 G提出了CLEANEST10]，能较好的处理研究过程中的多周期信号问题，有效的消除其中的伪周期。 小波变换处理等间距的过程中，Foster G又提出了向量投影的思想，把小波变换是向量的投影，那么研究的结果将会在很大程度上得到改善不但可以更准