晶振频率测量系统建模与数据分析-

论文导读：统组成测量系统由GPS接收机、晶振、时差测量模块、时钟产生模块、计算机数据采集处理组成。组成原理如图1所示。图1测量系统组成GPS接收机每秒输出1路TTL电平的标准秒脉冲(1PPS),晶振是时差测量和时钟产生的频率源。时钟产生电路产生本地秒脉冲。时差测量电路测量GPS秒脉冲与本地秒脉冲的相位差
【摘 要】以GPS秒脉冲作为测量标准,构建了晶振频率随时间变化的测量系统,通过对测量数据进行一元回归统计处理,分离出了晶振实际频率与其标称频率的相对偏差及晶振的各种随机误差,并分析了这两种误差对晶振准确度及稳定度的影响。该方法可为频率源误差测量分析提供借鉴作用。
【关键词】GPS;晶振;频率;误差;回归分析

1、晶振频率测量系统组成

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


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

2、模型建立

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

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

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

式中: 分别表示其他随机因素对变量y1,y2,…,,yN影响的总和,一般假设它们是一组相互独立,并服从同一正态分布N(0, )的随机变量。变量x在实验中为自然数,表示具体的秒脉冲数值。这样,变量y表示实际所测得的晶振与标准频率的计数差值,它是服从 的随机变量。
用最小二乘法来估计参数 。设b0,b分别是参数 的最小二乘估计,于是得到一元线性回归的回归方程:

式中:b0,b是回归方程的回归系数,分别表示晶振相对于摘自：毕业论文评语www.7ctime.com
标准频率的初始误差和累积误差。应用最小二乘法可求得回归系数b,b0为:

3、数据分析与处理

3.1 回归系数估计

为了定量分析数据,从而确定晶振频率误差的组成,首先对上述测得的数据进行归一化处理。实际测量中得到的是晶振脉冲的计数个数,设测量系统所用晶振频率为10MHz,可将计数数据转化为晶振相对于标准时间每秒的时间之差。
用Matlab对归一化数据进行处理,依照最小二乘原理,得到计数时间x与时间差值y的均值,以及x的自相关、x和y的互相关、y的自相关及回归方程如表2所示。

从表2得到回归方程为:

从回归方程可以看出,给定一个时间x值,就可得到相对应的时间差值。例如, 表示在第15s时晶振相对于标准时间的误差为t= 从回归方程可预测此晶振在24h相对于标准时间的时间误差值为:
以上求得了回归方程,但是该方程是否基本上符合y与x之间的客观规律,是否符合晶振频率误差变化的实际特点,还需要对回归方程做进一步的分析。在回归分析法中,通常采用方差分析法对回归方程的显著性进行检验,其实质是将N个测量值的影响从数量上区分开,然后用F检验法对所求回归方程进行显著性检验。

3.2晶振频率误差数据方差分析及显著性检验

测量值y1,y2,…,,yN之间的差异(称为变差)是由两方面的原因引起的。一是自变量取值的不同,二是其他因素(包括试验误差、随机误差等)的影响。为了对回归方程进行检验,把两者所引起的变差从y的总变差中分解出来。根据上述数据可得:

式中:U称为回归平方和,它反映了在y的总变差中由于x和y的线性关系而引起的y变化的部分;Q称为残余平方和,即所有测量点距回归直线的残余误差平方和。
若总的平方和由N项组成,其自由度就为N-1,总的离差平方和的自由度可分为回归平方和的自由度vU和残余平方和的自由度vQ之和,即:
vS=vU+vQ(6)
在一元线性回归中,vU=1,vS=N-1,则Q的自由度vQ=N-2。由回归平方和与残余平方和的意义可知,一个回归方程是否显著,也就是y与x的线性关系是否密切,取决于U及Q的大小,U愈大Q愈小,说明y与x的线性关系愈密切。通常可以采用F检验法来对方程进行显著性检验。对于一元线性回归,将U及Q的值代入上式得到统计量F:
由F分布表可以查出, 可认为回归是高度显著的。
残余平方和Q除以它的自由度vQ所得商:

称为残余方差,它可以看作排除了x对y的线性影响后,衡量y随机波动大小的一个估计量。残余方差的平方根:
称为残余标准差,它可用来衡量所有随机因素对y一次性测量平均变差的大小,愈小,则回归直线的精度愈高。
把平方和及自由度进行分解的方差分析数据结果归纳在一个表格中,如表3所示。

从表3可以看出,在30s时间内,晶振实际频率与其标称频率的相对偏差引起的误差平方和为,其他各种随机因素引起的误差的平方和为 。晶振相对频偏引起的误差远远大于其他因素引起的误差,表明了晶振频率误差贡献主要来源于实际频率与标称频率之间的频偏。通过频率修正即可得到一种高稳定度的频标源。
4结语
以上分析结果表明,该晶振频率误差主要来源于实际频率与标称频率的相对偏差,随机误差对晶振整体误差的贡献很小。通过对晶振标称频率值进行在线补偿可以消除其相对频偏,从而获得一种具有较高稳定度的频率源,可为需要时间显示的场所提供高精度的时间服务。
参考文献
曾祥君,尹项根,林干,等.晶振信号同步GPS信号产生高精度时钟的方法及实现[J].电力系统自动化,2003,27(8):49-51.
寇艳红,张其善.GPS接收机中晶振误差模拟方法[J].电子与信息学报,2004,26(8):1320-1324.
注：文章内所有公式及图表请用PDF形式查看。