振荡器频率稳定度的小波反演的研究.pdfVIP

振荡器频率稳定度的小波反演研究 柯熙政 中国科学院陕西天文台西安710600电话：029—3890449 cn E-mail：kexzh@mS．sxso．ac 摘要：本文根据小渡方差和Allan方差的联系，提出了由振荡嚣频率稳定度反演振荡器 频率的新方法．最后对数值模拟的误差进行了讨论． 关键词：振荡嚣频率稳定度小波变换 1 引言 原子钟是振荡器的一种。实际的振荡器不但表现出一个在频率标称值附 近的随机变化，不失一般性，我们可以将振荡器的额率表示为¨ 撕)=儿(f)+盯f+“ 儿是频率源的标称值，a是频率源的老化系数，儿O)是各种不确定的随机因 素的影响，是一个真正的随机过程振荡器频率稳定度的表征首先遇到的问题 是频率不可能瞬时测量，一般我们讨论持续时间为r的N次测量，其测量间 隔为T=k．一小测量间隙为r—r，测量序列可以表示为 n“，f)=抗，元，…，雨} (2) N次取样方差为 毋∽，T。r)=击耋魄一◇t)7 c，) 司伍tf)={(◇¨-歹ty) (4) 取样方差就与N．T．r有关，N=2，r=r时的Allan方差已经被推荐为频率稳定 度的表征量．振荡器的噪声可以用分形布郎运动函数来描述f“，虽然分形布 郎运动函数是非平稳的，但其增量函数确是平稳的。Allan方差实质是与过 程的增量有关的． 2小波方差 ．对随机过程y=几)E尺，其小波方差定义为 z∞=三卜，o叫‰ (5) 19 式中w，(口，6)=三，∞去y(譬卜是函数．厂o)的小波变换，n是小波尺度，6是 平移因子，y(．)是小波基函数。对离散小波变换，当小波基函数为Haar基 时，其小渡方差就是舢1an方差D】． 3振荡器频率稳定度反演的数学基础 平稳随机过程由其一阶和二阶矩确定，但非平稳过程的二阶矩是时变 的。一般的根据过程的均值和方差的模拟方法，仅能保证所,-f,j1．勘x的lili．,t凡过程 具有所要求的均值和方差，不能保证达到振荡器的频率稳定度．原子钟的噪 声一般呈幂律谱分布。它可由其核函数积分而得到14j． 3．1I-Iaar基小波方差和AJlan方差 对随机过程y=航废，…，霸)，Haar基小波变换系数是一时间序列 抚一只，只一只，…，霜一歹。，其方差可组成一个时间序列_乃，，…，侈^，一歹_Ⅳ一声， 其均值就是小波方差： 一：={(眵2一箩1n，‰节¨，) (6) 小波方差的标准差“是方差序列的均值，在Haar基下，小波方差与Allan 方差是等价的．式中且表示在测量持续期间r内的平均值．对同一个原子钟， 当测量持续时间不同时，可以得到不同的频率稳定度袁征值． 3．2频率稳定度的反演 欲模拟Allan方差为，：的振荡器噪声，假定取样方差为(0，2棚间正态分 布的随机数，这可很容易由同余法等计算机模拟方法得到： 眵2-Y17，…，侈’Ⅳ一罗^r一声，与之相对应，具有频率稳定度蠢的振荡器的的频率 采样应为：妒2一，1)，…，伊Ⅳ一歹’_Ⅳ一1)对其进行逆小波变换就可得到振荡器测量 值的一个实现p-．，，·：，，·，}，_n_-}i-关系 铲，，『，：，…，于，)=∥rt眵2—，1)，…，舻_Ⅳ一歹’Ⅳ一1) (7) 式申胛一表示逆小波变换． 3．3频率采样的内插与外推 对测量持续周期r。，我们可以根据Allan方差∥(2，丁，，r。)，由(6)式可以得 到模拟测量序列：靓臼。lF，：∽l