用加权小波分析方法计算的周期研究模板.docVIP

资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。

用WWZ加权小波分析方法计算S50716+714的周期研究

摘要:本文采用的研究方法是加权小波Z变换(WeightedwaveletZ-transform,WWZ),它是一种基于小波分析原理处理非等间隔数据的方法。经过收集BLLac天体S50716+714光学B、V、R、I四个波段较为完备的观测数据,得到了其一定时期内的光变曲线。利用WWZ对S50716+714以上四个波段的光变数据进行计算分析研究,发现S50716+714存在的光变周期为3.4±0.1年。

关键词:小波分析;WWZ变换;BLLac天体

1引言

活动星系核是现代天体物理学中的一个热门研究领域,活动星系核简称AGNs(ActiveGalacticNucleus)。在有关活动星系核的研究中,涉及到很多有关天体物理的基础性问题,如能量的产生、辐射机制等。类星体是活动星系核的一个重要子类,在研究中发现,类星体具有许多的极端特性。而BLLac天体又是类星体中的一个重要代表,因为其在观测上表现出高光度,高偏振,快速光变以及非热辐射等特征[1]。S50716+714是BLLac天体之一,备受各界科学研究者的关注。对于活动星系而言,光变周期的研究是了解其重要特性和本质的一个非常重要的方法。经过对光变周期等特性的研究,我们能够了解天体辐射区域内部物理条件和结构的变化,辐射介质以及辐射区域外物理环境的变化和辐射过程的改变等[2]。周期光变是类星体的一个普遍的特征,在近几年来也有不少研究组提出了多种研究类星体长周期光变的方法,如解释长周期光变的理论模型,双黑洞喷流模型,进动吸积盘模型等。

用来寻找计算天体周期光变的方法有很多,经常见到的有时间序列的功率分析方法[3],Jurkevich方法[4],相位分析法,自相关函数分析法等[5]。虽然以上这些方法在研究天体周期光变中比较广泛,可是也存在许多的条件限制和局限性,在天体光变周期的研究中将会有所困难。因此,本文中将使用一种较为新颖的方法加权小波Z变换对BLLac天体S50716+714的光变周期进行计算和研究。

WWZ是一种基于小波分析原理产生的能够用于非规则数据处理的新方法[6]。小波分析方法是一种非平稳信号分析方法,它和功率谱分析方法都为傅立叶变换的周期分析方法。可是小波分析方法是将信号表示为时间尺度-频率尺度的域解析。因此在整个频率域中,具有多频率分辨的特点。由于小波分析的主要核心部分是多分辨率分析,因此较为适合分析非平稳信号,因此用WWZ来分析非平稳信号是一种较好的方法。

本文将采用WWZ变换,对S50716+714的四个光学波段周期数据进行计算研究,最后得出S50716+714的光变周期的部分变化特性。

2加权小波Z变换(WWZ)

2.1小波分析原理

传统的信号理论,是建立在傅立叶分析基础上的,而傅立叶变换作为一种全局性的变化,其有一定的局限性。小波变换与傅立叶变换相比,是一个时间和频域的局域变换,因而能有效地从信号中提取信息,经过伸缩和平移等运算功能对函数或信号进行多尺度细化分析(MultiscaleAnalysis),解决了傅立叶变换不能解决的许多困难问题。小波变换是一种时间与尺度分析的方法,它具有时频局部化的功能,在低频部分具有较高的频率分辨率和较低的时间分辨率,在高频部分具有较高的时间分辨率和较低的频率分辨率,非常适合用于探测夹带瞬间反常现象的信号。

小波变换的定义[7]为:

设,则它的连续小波变换为:

(1)

其中a为伸缩尺度,b为平移参数。

将其离散化形式为:

(2)

称为小波函数,它是由小波母函数经过平移b和伸缩a得到。是一种长度有限,平均值为0的波形。它具有快速衰减的特性。一般所用的小波母函数为Morlet小波,它的形式为[8]:

(3)

是衰减因子,当取较大值时上式第二项能够省略,因此有简化的Morlet小波:

(4)

Torrenceetal指出对于Morlet小波伸缩尺度a和傅立叶变换的频率之间有关系[9]:

(5)

2.2WWZ原理

小波分析方法能够较好的处理非平稳信号,可是在处理不规则的非等间距数据时,有一定的缺陷。因为在(2)式和(1)式较为近似,在实际的应用中