基于比时法的晶振频率测量建模与分析

摘 要:受器件老化、随机噪声等因素影响,晶振频率变化较复杂。以GPS秒脉冲作为测量标准,构建了晶振频率随时间变化的测量系统,通过对测量数据进行一元回归统计处理,分离出了晶振实际频率与其标称频率的相对偏差及晶振的各种随机误差,并分析了这两种误差对晶振准确度及稳定度的影响。该方法可为频率源误差测量分析提供借鉴作用。

关键词:GPS;晶振;频率;误差;回归分析

中图分类号:TM930 文献标识码:A

文章编号:1004-373X(2010)03-119-03

Modeling and Analysis of Crystal-oscillator Frequency Based on Measuring Time

LI Erpeng,WEN Kaizhang,WANG Weiwei

(Northwest Institute of Nuclear Technology,Xi′an,710024,China)

Abstract:The characteristics of oscillator frequency error are very complex owing to the aging and stochastic noise.The system measure model of crystal oscillator frequency error is designed with GPS as frequency standard.A method of regression analysis is used for disposing the oscillator frequency error data and separating the linearity excursion error from the stochastic error.Both two errors are analyzed in contribution of the crystal oscillator stabilization and precision.This mean is provided with reference in frequency criterion error of measure and analysis.

Keywords:GPS;crystal-oscillator;frequency;error;regression analysis

0 引 言

石英晶体振荡器受制造工艺、器件老化以及外部温度等因素影响,其实际频率值与标称频率值存在偏差。此外,晶体振荡器内部存在着各种噪声也会使频率值产生随机起伏,从而导致晶振频率的准确度和稳定度降低。当前,国内外在研究晶体振荡器的老化和随机噪声、分析晶振频率误差特性方面提出了很多方法,比较常见的有时间对数线性模型法、自适应滤波法及非线性时变预测法。前两种方法的缺点是参数较多,选择合适的参数较难[1,2],非线性时变预测法的模型描述能力强,但由于没有函数参数的显式表达式,不能将模型求解归结为参数求解问题,一般通过学习来逼近该函数,主要用于频率变化的预测[3]。

在上述方法的基础上,本文利用CPLD设计了┮恢值缏,该电路采用比时法[4,5]来测量晶振频率变化,并根据其频率随时间变化拟合曲线的特点,用线性回归法分析其频率误差特性。该方法模型简单,参数易于估计,可通过简洁的补偿方法消除晶振相对频偏,具有实际的应用价值。

1 晶振频率测量系统组成

测量系统由GPS接收机、晶振、时差测量模块、时钟产生模块、计算机数据采集处理组成。组成原理如┩1所示。

图1 测量系统组成

GPS接收机每秒输出1路TTL电平的标准秒脉冲(1 PPS),晶振是时差测量和时钟产生的频率源。时钟产生电路产生本地秒脉冲。时差测量电路测量GPS秒脉冲与本地秒脉冲的相位差值。

采用比时法测量晶振频率的系统工作原理如下:首先由晶振分频得到本地的1 Hz频率源,将GPS秒脉冲与本地晶振秒脉冲送入时差测量模块进行相位比较,得到两者的相位差信号,设计时间间隔计数器[6]对此相差闸门信号计数,计数值即为晶振频率相对于标准频率的偏差,反映了晶振频率的误差特征。时间间隔计数器每秒测量一次两者的偏差值,需要测量的频差数据量由计算机设定,测量结果传送到计算机进行数据统计,并对晶振频率误差特性进行分析。本地时钟产生、时差测量及数据采集电路模块等都由CPLD设计实现[7,8]。

2 模型建立

实验以GPS秒脉冲作为标准参考信号,采用比时法对晶振的输出频率进行测量。按照建立的测量系统,实际对某一晶振采集到30个计数值,如表1所示。

表1 实验数据

x/s12345678910

y/个79156235313392470548626704783

x/s11121314151617181920

y/个8619401 0181 0971 1751 2541 3321 4101 4881 566

x/s21222324252627282930

y/个1 6451 7221 8011 8791 9562 0342 1132 1932 2702 349

表1中x表示测量的时间(单位:s);y表示晶振脉冲计数个数。假设晶振在某秒计数值为M,它的计数周期为T,则MT为晶振秒脉冲与GPS秒脉冲的时间差值。例如在第6 s时测得计数值为470,则表示在┑6 s时晶振与GPS秒脉冲的时间差为470T。为了研究时间x与计数个数y之间的关系,用ORIGIN工具软件对数据进行拟合处理,得到的x,y关系曲线如图2所示。

图2 数据拟合曲线

从散点图可以看出,测量计数值和测量时间大致呈线性关系。据此假设这两个变量之间的内在关系是┮惶踔毕,这些点与直线的偏离是由于测量过程中其他一些随机因素的影响而引起的[9],这样可以假设这组测量数据有如下结构形式:

yt=β0+βxt+εt,t=1,2,…,N

(1)

式中:Е弄1,ε2,…,εN分别表示其他随机因素对变量y1,y2,…,yN影响的总和,一般假设它们是一组相互独立,并服从同一正态分布N(0,δ)的随机变量。变量x在实验中为自然数,表示具体的秒脉冲数值。这样,变量y表示实际所测得的晶振与标准频率的计数差值,它是服从N(β0+βxt,δ)У乃婊变量。

用最小二乘法[10]来估计参数Е陋0,β。设b0,b分别是参数β0,βУ淖钚《乘估计,于是得到一元线性回归的回归方程:

=b0+b

(2)

式中:b0,b是回归方程的回归系数,分别表示晶振相对于标准频率的初始误差和累积误差。应用最小二乘法可求得回归系数b,b0为:

b=N∑Nt=1xtyt-(∑Nt=1xt)(∑Nt=1yt)N∑Nt=1x2t-(∑Nt=1xt)2=lxylxx

(3)

b0=(∑Nt=1x2t)(∑Nt=1yt)-(∑Nt=1xt)(∑Nt=1xtyt)N∑Nt=1x2t-(∑Nt=1x┆t)2=-b