CN106026919A - 高精度晶体振荡器的守时补偿方法 - Google Patents
高精度晶体振荡器的守时补偿方法 Download PDFInfo
- Publication number
- CN106026919A CN106026919A CN201610324056.1A CN201610324056A CN106026919A CN 106026919 A CN106026919 A CN 106026919A CN 201610324056 A CN201610324056 A CN 201610324056A CN 106026919 A CN106026919 A CN 106026919A
- Authority
- CN
- China
- Prior art keywords
- error
- temperature
- crystal oscillator
- gps
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03B—GENERATION OF OSCILLATIONS, DIRECTLY OR BY FREQUENCY-CHANGING, BY CIRCUITS EMPLOYING ACTIVE ELEMENTS WHICH OPERATE IN A NON-SWITCHING MANNER; GENERATION OF NOISE BY SUCH CIRCUITS
- H03B5/00—Generation of oscillations using amplifier with regenerative feedback from output to input
- H03B5/02—Details
- H03B5/04—Modifications of generator to compensate for variations in physical values, e.g. power supply, load, temperature
Landscapes
- Oscillators With Electromechanical Resonators (AREA)
- Electric Clocks (AREA)
Abstract
本发明提供一种高精度晶体振荡器守时补偿方法,包括:通过FPGA采集板卡对晶体振荡器的频率进行处理,得到晶体振荡器与GPS秒脉冲之间的对时误差,测量晶体振荡器的工作温度,构造温度特性曲线,得到晶体振荡器温度特性模型,对温度补偿之后的对时误差数据进行卡尔曼滤波,得到构造的基准时钟,从而获得GPS误差特性曲线,在分离了温度误差及GPS误差之后,采用最小二乘法预测守时补偿值。本发明根据GPS接收机的误差特性、晶体振荡器的相位随机误差特性、温度频率偏移误差特性以及老化特性等,对晶振守时误差进行晶振温度‑频率变化特性模型的重构,并使用卡尔曼滤波的方法构建基准时钟,以此分离各误差源的特性,得到高精度的守时补偿值。
Description
技术领域
本发明涉及电力***继电保护及同步采样技术领域,具体而言涉及基于实时温度测量更新晶体振荡器温度模型参数,以此构造基准时钟模型并生成守时补偿值的方法。
背景技术
基于GPS的晶振驯服技术是利用GPS接收机输出的秒脉冲信号锁定本地时钟,当GPS模块失锁或者秒脉冲信号丢失之后,则主要依靠晶振自身的稳定度及预设的补偿量来保证晶振在守时阶段的准确度。因此,当晶振自身的稳定度一定时,对于守时阶段的补偿算法的研究具有重要意义。
目前对于温度补偿晶体振荡器的守时补偿算法主要是针对算法本身的改进,如采用改进kalman滤波或者滑动平均等算法来实现守时。由于kalman滤波对于已知误差模型的预测较为有效,滑动平均算法则主要应用于具有正态分布特性的误差中,而守时误差中包含各类误差特性,若将其统一处理,将会导致由于模型不准确而造成的预测偏差。另外还有利用历史数据采用拟合的方式预测守时补偿值。预测需基于准确的模型,若采用各误差特性相叠加的对时误差进行分析处理,会导致由于误差特性混叠从而加大预测算法的复杂度,同时降低预测准确度。因此迫切需要提出一种分离各个误差特性的守时方法。
发明内容
本发明的目的在于解决目前守时算法将误差统一建模处理而未考虑各误差源的分布特性导致守时精度降低的问题。
为解决上述技术问题,本发明所采用的技术方案如下:
一种高精度晶体振荡器守时补偿方法,包括以下步骤:
通过FPGA采集板卡对晶体振荡器的频率进行处理,得到晶体振荡器与GPS秒脉冲之间的对时误差,采用高精度温度变送器测量晶体振荡器的工作温度,构造温度特性曲线,得到晶体振荡器温度特性模型,对温度补偿之后的对时误差数据进行卡尔曼滤波,得到构造的基准时钟,从而获得GPS误差特性曲线,在分离了温度误差及GPS误差之后,采用最小二乘法预测守时补偿值。
在一些实施例中,具体步骤为:
第一步,以晶体振荡器的频率作为基准,对振荡器自身时钟进行计数,得到相对时标(程序运行开始累计的时间值)。同时,判断GPS接收机的时钟是否有效,若有效进一步判断GPS秒脉冲是否失准,得到以GPS秒脉冲为基准的晶体振荡器对时误差。
第二步,得到晶体振荡器的初始温度频率特性曲线。
选取测量的温度数据,将其中的最大值及最小值去除(避免温度传感器的突变量对测量造成较大影响),其余值取算术平均。同时对振荡器的频率处理方式也相同,将最大与最小温度值对应的频率去除,其余频率值取算术平均:
其中,Ni为第i秒的晶体振荡器频率计数值,n为计数次数。将频率与温度的采样值代入如下方程中进行求解:
式中,Δfi为对应温度Ti时测得的频率fi和f0的差值;ΔTi为实时温度Ti与T0的差值;ΔE频率修正量。通过引入频率修正量实现对三次函数的a0、b0、c0三个温度特性曲线的参数进行更新修正,得到晶体振荡器的初始温度频率特性曲线。
第三步,得到晶体振荡器温度频率特性曲线。
在滤除GPS接收机及振荡器频率偏差的综合误差曲线中的高频随机分量之后,采用滤波后的数据取负值作为基准时钟信号,再次采集温度及振荡器频率,进而计算出新的温度频率特性的各次参数。重复上述步骤,直到温度频率特性参数更新值小于设定范围,则停止更新。最终得到的温度频率特性曲线。
第四步,采用卡尔曼滤波得到基准时钟。
实际测得的误差量包含了GPS的随机抖动及晶体振荡器的各误差特性曲线,其表达式为
XGPS=ErrGPS+X(3)式中,XGPS表示总的误差曲线;ErrGPS为GPS接收机的随机抖动;X为晶体振荡器初始温度特性曲线造成的频率偏移。
采用卡尔曼滤波算法对XGPS进行预处理,。卡尔曼滤波的核心由5个基本公式组成,主要可以分为时间更新方程和测量更新方程。其中时间更新方程负责及时向前推算当前状态变量和误差协方差估计值,以便为下一个时间状态构造先验估计,测量更新方程将先验估计和新的测量变量结合以构造改进的后验估计。
时间更新方程:
状态更新方程:
式中,为k时刻***观测值,为k-1时刻对***
的控制量,A、B为***参数,Pk -为误差协方差,Q为观测随机误差的方差;Kk为卡尔曼增益,H为测量***参数,Zk为k时刻的测量值,为k时刻的预测值。
经过卡尔曼滤波处理可以滤除GPS接收机随机抖动中的高频分量,从而可以得到一个较为准确的时间基准。
第五步,构造GPS误差特性曲线。
对晶体振荡器计数模块进行两路并行频率计数,其中一路采用GPS为基准的频率计数得到Err1,另一路采用之前经过卡尔曼滤波处理之后的数据为基准进行频率计数得到Err2,其表达式为:
Err1=S-XGPS=S-ErrGPS-X (6)
Err2=S-X(7)式中S为第四步构建出的较为准确的基准时钟,Err1为GPS接收机有效时的守时误差,Err2表示GPS丢失之后,利用S作为基准时钟而产生的守时误差,XGPS表示GPS接收机与晶体振荡器的综合误差,ErrGPS表示GPS接收机的误差。
由于在守时***中采用两路并行频率计数的形式,Err2可以直接得到,Err1经过滤波及求差之后得到的。因此,由式(6)、式(7)可以得到GPS的误差ErrGPS为:
ErrGPS=Err2-Err1(8)
第六步,在重构了GPS误差特性曲线的基础上,根据式
X=XGPS-ErrGPS(9)
得到晶体振荡器总的误差特性曲线X,利用温度传感器及构建的晶体振荡器温度特性模型,对晶体振荡器温度特性进行补偿,得到温度补偿之后的晶体振荡器误差曲线为Xtemp。当经过温度补偿之后,该晶体振荡器误差模型就可以表示为:f0'=fB+f1(x0)+ΔT+f3(t)+δ(10)式中,ΔT为考虑到温度传感器的误差;f1(x)为晶体振荡器的随机相位偏差;f3(x)为晶体振荡器的老化偏差,δ为其他***误差;f0’为进行了温度补偿之后的晶体振荡器输出频率。
对式(10)中的f0’进行最小二乘预测。最小二乘法的表达式为:
式中,Xpre是采用最小二乘预测的守时数据。N为进行最小二乘算法采用的数据量。拟合得到的直线反映晶体振荡器频率的变化,同时也为生成补偿值提供依据。
应当理解,前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外,所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。
结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见,或通过根据本发明教导的具体实施方式的实践中得知。
附图说明
附图不意在按比例绘制。在附图中,在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见,在每个图中,并非每个组成部分均被标记。现在,将通过例子并参考附图来描述本发明的各个方面的实施例,其中:
图1是分频计数的流程图。
图2是基于GPS的晶体振荡器历史数据的示意图。
图3是最小二乘拟合结果的示意图。
图4是基于误差分离的改进守时算法仿真曲线的示意图。
图5是总仿真曲线的示意图。
图6是晶体振荡器守时结构图。
具体实施方式
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。
在本公开中参照附图来描述本发明的各方面,附图中示出了许多说明的实施例。本公开的实施例不必定意在包括本发明的所有方面。应当理解,上面介绍的多种构思和实施例,以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。
根据本发明的实施例,提出一种高精度温度补偿晶体振荡器的守时方法。该方法综合考虑了GPS接收机的误差、晶振振荡器的相位随机误差、温度频率偏移误差以及老化特性误差,利用对时误差的历史数据,通过误差反馈更新基准时钟源,得到精确的晶体振荡器温度特性模型,进一步以更新生成的时钟源为基准得到晶体振荡器误差特性,通过最小二乘法预测出守时补偿值。
在一些具体的实施例中,我们下面将结合附图1-附图6来更加详细地说明前述方法的实现。
第一步,以晶体振荡器的频率作为基准,对振荡器自身时钟进行计数,得到相对时标(程序运行开始累计的时间值)。同时,判断GPS接收机的时钟是否有效,若有效进一步判断GPS秒脉冲是否失准,得到以GPS秒脉冲为基准的晶体振荡器对时误差。程序流程图如图1所示。
第二步,得到晶体振荡器的初始温度频率特性曲线:
采用标称频率为25MHz的温度补偿晶体振荡器,其晶体为AT切型,温度频率偏移特性为三次函数:
式中,f0为参考温度点的频率,可以设置零度或者室温为参考温度点;Δf为实际温度与参考温度点之间的偏差;T为实测温度;T0是参考温度;A0、B0、C0分别为温度特性曲线的一次、二次、三次温度系数。因此,要得到AT切型晶体振荡器的温度特性曲线只需确定其三个温度系数即可。GPS接收机标准差为25ns。历史数据是基于GPS接收机的晶体振荡器输出频率偏差数据,其24小时内的历史数据如图2所示(1秒为间隔进行采样)。
选取测量的温度数据,采取连续取5s的数据,将其中的最大值及最小值去除(避免温度传感器的突变量对测量造成较大影响),其余值取算术平均。同时对振荡器的频率处理方式也相同,将最大与最小温度值对应的频率去除,其余频率值取算术平均:
其中,Ni为第i秒的晶体振荡器频率计数值,n为计数次数。将频率与温度的采样值代入如下方程中进行求解:
式中,Δfi为对应温度Ti时测得的频率fi和f0的差值;ΔTi为实时温度Ti与T0的差值;ΔE频率修正量。通过引入频率修正量实现对三次函数的a0、b0、c0三个温度特性曲线的参数进行更新修正,得到晶体振荡器的初始温度频率特性曲线。
第三步,得到晶体振荡器温度频率特性曲线:
采用卡尔曼滤波滤除GPS接收机及振荡器频率偏差的综合误差曲线中的高频随机分量,采用滤波后的数据取负值作为基准时钟信号,再次采集温度及振荡器频率,计算出新的温度频率特性的各次参数。重复上述步骤,直到温度频率特性参数更新值小于设定范围,则停止更新,最终得到温度频率特性曲线。为了避免得到病态方程,设置温度变化在0.1度以上时再进行数据更新。
第四步,采用卡尔曼滤波得到基准时钟:
实际测得的误差包含了GPS的随机抖动及晶体振荡器的各误差特性曲线,其表达式为:
XGPS=ErrGPS+X(4)式中,XGPS表示总的误差曲线;ErrGPS为GPS接收机的随机抖动;X为晶体振荡器使用初始温度特性曲线所造成的频率偏移。
采用卡尔曼滤波算法对XGPS进行预处理。卡尔曼滤波的核心由5个基本公式组成,主要可以分为时间更新方程和测量更新方程。其中时间更新方程负责及时向前推算当前状态变量和误差协方差估计值,以便为下一个时间状态构造先验估计,测量更新方程将先验估计和新的测量变量结合以构造改进的后验估计。
时间更新方程:
状态更新方程:
式中,为k时刻***观测值,为k-1时刻对***
的控制量,A、B为***参数,Pk -为误差协方差,Q为观测随机误差的方差;Kk为卡尔曼增益,H为测量***参数,Zk为k时刻的测量值,为k时刻的预测值。
经过卡尔曼滤波处理滤除GPS接收机随机抖动中的高频分量,从而得到一个较为准确的时间基准。
第五步,构造GPS误差特性曲线:
对晶体振荡器计数模块进行两路并行频率计数,其中一路采用GPS为基准的频率计数得到Err1,另一路采用之前经过卡尔曼滤波处理之后的数据为基准进行频率计数得到Err2,其表达式为:
Err1=S-XGPS=S-ErrGPS-X (7)
Err2=S-X(8)式中S为第四步构建出的较为准确的基准时钟,Err1为GPS接收机有效时的守时误差,Err2表示GPS丢失之后,利用S作为基准时钟而产生的守时误差,XGPS表示GPS接收机与晶体振荡器的综合误差,ErrGPS表示GPS接收机的误差。
由于在守时***中采用两路并行频率计数的形式,Err2可以直接得到,Err1经过滤波及求差之后得到的。因此,由式(7)、式(8)可以得到GPS的误差ErrGPS为:
ErrGPS=Err2-Err1(9)
第六步,在重构了GPS误差特性曲线的基础上,根据式
X=XGPS-ErrGPS(10)
得到晶体振荡器总的误差特性曲线X,利用温度传感器及构建的晶体振荡器温度特性模型,对晶体振荡器温度特性进行补偿。得到温度补偿之后的晶体振荡器误差曲线为Xtemp。当经过温度补偿之后,该晶体振荡器误差模型就可以表示为:f0'=fB+f1(x0)+ΔT+f3(t)+δ(11)式中,ΔT为考虑到温度传感器的误差;f1(x)为晶体振荡器的随机相位偏差;f3(x)为晶体振荡器的老化偏差,δ为其他***误差;f0’为进行了温度补偿之后的晶体振荡器输出频率。
对于直线拟合y=a*x+b
根据直线最小二乘拟合公式:
b=E(y)-a*E(x)
所以:
这里E(x)为求均值。
根据上述公式,对式(11)中的f0’进行最小二乘预测。最小二乘法的表达式为:
式中xi为ti时刻对应的晶体振荡器误差值,Xpre是采用最小二乘预测的守时数据,N为进行最小二乘算法采用的数据量。其中
拟合得到的直线反映晶体振荡器频率的变化,同时也为生成补偿值提供依据。通过多次实验反复调整式中的参数,以N=200为例对上述频率特性曲线进行最小二乘拟合,结果如图3所示。基于误差分离的守时算法的守时曲线如图4所示。
通过仿真来对比本发明提出的基于误差分离的改进守时算法与仅基于温度补偿的算法以及基于预测的守时算法,得到10小时守时曲线如图5所示。
由图5曲线(a)可见,由于传统的守时补偿方法没有考虑晶体振荡器的温度特性以及GPS的特性,其10小时守时误差为165微秒。在图5曲线(b)中是基于温度补偿之后的晶体振荡器守时误差曲线,由图可知,其10小时守时误差在45微秒以内。这是由于对晶体振荡器进行了建模,考虑了温度的影响,但是由于模型本身的误差,以及忽略了GPS的随机抖动,其守时精度仍然具有改进空间。图5曲线(c)采用本发明的守时算法,可见在基准信号GPS丢失10小时后,守时误差呈现一个抖动的下降趋势,基本保证10小时守时误差在2微秒以内。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (5)
1.一种晶体振荡器的守时补偿方法,其特征在于,包括以下步骤:
1)采用温度变送器测量晶体振荡器的工作温度,通过FPGA采集板卡对晶体振荡器的频率进行处理得到对时误差;
2)根据采集的工作温度数据,构造晶体振荡器的温度特性曲线,更新晶体振荡器的温度特性曲线参数,得到晶体振荡器温度特性模型并进行温度-频率误差补偿;
3)将温度-频率误差补偿之后的对时误差数据进行卡尔曼滤波,得到构造的基准时钟,从而获得GPS误差特性;
4)在分离了温度误差及GPS误差之后,采用最小二乘法预测守时补偿值。
2.根据权利要求1所述的晶体振荡器守时补偿方法,其特征在于,在所述步骤1)中,通过温度传感器获取晶体振荡器的工作温度,并在FPGA采集板卡设置分频电路,得到晶体振荡器的频率与GPS秒脉冲之间的对时误差。
3.根据权利要求2所述的晶体振荡器守时补偿方法,其特征在于,在根据前述步骤1)得到晶体振荡器的实时温度及其相应频率值后,在所述步骤2)中,按照下述方式处理进行温度-频率误差补偿:
先设定初始温度特性参数,然后对其进行修正以获取温度特性模型;其中在选取测量数据构造温度特性模型时,采用连续若干点的数据,将其中的最大值及最小值去除,然后将其余值取算术平均,同时对晶体振荡器的频率处理方式也相同,将最大与最小温度值对应的频率去除,然后将其余频率值取算术平均,具体处理如下:
其中,Ni为第i秒的晶体振荡器频率计数值,n为计数次数;
将频率与温度的采样值代入如下方程中进行求解:
式中,Δfi为对应温度Ti时测得的频率fi和f0的差值;f0为基准频率,ΔTi为实时温度Ti与T0的差值;ΔE频率修正量,通过引入该修正量实现对三次函数的a0、b0、c0三个温度特性曲线的参数进行更新修正,实现温度-频率误差补偿修正。
4.根据权利要求1所述的高精度晶体振荡器守时补偿方法,其特征在于,在所述步骤3)中,采用卡尔曼滤波得到基准时钟,具体如下:
测得的数据为基于GPS时钟的晶体振荡器频率计数值,即测得的数据为XGPS,其中包含了GPS的随机抖动及晶体振荡器的各误差特性曲线:
XGPS=ErrGPS+X (3)
式中,XGPS表示总的误差曲线,也即实际测得的误差量;ErrGPS为GPS接收机的随机抖动;X为晶体振荡器初始温度特性曲线造成的频率偏移;
采用卡尔曼滤波算法对XGPS进行预处理,卡尔曼滤波的核心由5个基本公式组成,包括时间更新方程和测量更新方程,其中时间更新方程负责及时向前推算当前状态变量和误差协方差估计值,以便为下一个时间状态构造先验估计,测量更新方程将先验估计和新的测量变量结合以构造改进的后验估计;
时间更新方程:
状态更新方程:
式中,为k时刻***观测值,为k-1时刻对***的控制量,A、B为***参数,为误差协方差,Q为观测随机误差的方差;Kk为卡尔曼增益,H为测量***参数,Zk为k时刻的测量值,为k时刻的预测值。
经过卡尔曼滤波处理之后的数据,滤除GPS接收机随机抖动中的高频分量,从而可以得到一个准确的时间基准。
5.根据权利要求4所述的高精度晶体振荡器的守时补偿方法,其特征在于:
在重构了GPS误差特性曲线的基础上,便可得到晶体振荡器总的误差特性曲线X:
X=XGPS-ErrGPS (6)
在得到晶体振荡器总的误差曲线之后,基于温度传感器及构建的晶体振荡器温度特性模型,对晶体振荡器温度特性进行补偿;
得到温度补偿之后的晶体振荡器误差曲线为Xtemp,当经过温度补偿之后,该晶体振荡器误差模型就可以表示为:
f0'=fB+f1(x0)+ΔT+f3(t)+δ (7)
式中,ΔT为考虑到温度传感器的误差,本文取温度变送器的分辨率对应的误差;f1(x0)为晶体振荡器的随机相位偏差;f3(t)为晶体振荡器的老化偏差,δ为其他***误差;f0’为进行了温度补偿之后的晶体振荡器输出频率;
然后,对式(7)中的f0’进行最小二乘预测,最小二乘法的表达式为:
式中,xi为ti时刻对应的晶体振荡器误差值,Xpre是采用最小二乘预测的守时数据,N为进行最小二乘算法采用的数据量,其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610324056.1A CN106026919B (zh) | 2016-05-16 | 2016-05-16 | 晶体振荡器的守时补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610324056.1A CN106026919B (zh) | 2016-05-16 | 2016-05-16 | 晶体振荡器的守时补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106026919A true CN106026919A (zh) | 2016-10-12 |
CN106026919B CN106026919B (zh) | 2019-05-07 |
Family
ID=57098421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610324056.1A Active CN106026919B (zh) | 2016-05-16 | 2016-05-16 | 晶体振荡器的守时补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106026919B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106706021A (zh) * | 2016-11-15 | 2017-05-24 | 中国电子科技集团公司第四十研究所 | 一种基于调制域分析仪的线性调频测量方法 |
CN107276582A (zh) * | 2017-05-17 | 2017-10-20 | 电子科技大学 | 一种基于模拟电路的温度补偿晶体振荡器 |
CN107621644A (zh) * | 2017-08-09 | 2018-01-23 | 上海移为通信技术股份有限公司 | 一种gps定位***的抗干扰控制方法 |
CN108900161A (zh) * | 2018-05-21 | 2018-11-27 | 阮晓萍 | 一种晶振器时钟温度误差补偿方法 |
CN109100930A (zh) * | 2017-06-21 | 2018-12-28 | 卡西欧计算机株式会社 | 电子设备、日期时间取得控制方法及记录介质 |
CN110262210A (zh) * | 2019-06-28 | 2019-09-20 | 北斗天汇(北京)科技有限公司 | 基于计数器的晶振守时方法 |
CN110412438A (zh) * | 2018-04-26 | 2019-11-05 | 三星电子株式会社 | 生成包括晶体振荡器的测试产品的晶体模型的方法和*** |
WO2019232610A1 (en) * | 2018-06-06 | 2019-12-12 | Microsemi Semiconductor Ulc | Non-linear oven-controlled crystal oscillator compensation circuit |
CN110740001A (zh) * | 2018-07-18 | 2020-01-31 | 北京松果电子有限公司 | 晶振频率校准方法、装置以及电子设备 |
CN111711446A (zh) * | 2020-06-10 | 2020-09-25 | 国网湖南省电力有限公司 | 一种利用gps信号驯服晶振频率的方法、***及介质 |
CN112213541A (zh) * | 2020-09-08 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | 一种应用于暂态录波型故障指示器的守时方法 |
CN112485520A (zh) * | 2020-12-03 | 2021-03-12 | 成都市精准时空科技有限公司 | 基于电压采样的绝对频差测量方法及***及装置及介质 |
CN113067578A (zh) * | 2021-03-17 | 2021-07-02 | 成都金诺信高科技有限公司 | 一种时钟源的精密补偿方法 |
CN113422593A (zh) * | 2021-07-05 | 2021-09-21 | 北京信息科技大学 | 滤波方法、滤波器、计算机可读存储介质、处理器和fpga |
CN114362140A (zh) * | 2021-12-01 | 2022-04-15 | 国网内蒙古东部电力有限公司通辽供电公司 | 一种适用于配电网量测装置的高精度守时方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079687A (zh) * | 2006-05-25 | 2007-11-28 | 上海欣泰通信技术有限公司 | 一种基于最小二乘法模型的时钟调整算法 |
CN102629103A (zh) * | 2012-03-30 | 2012-08-08 | 上海鸿晔电子科技有限公司 | 卫星授时***的网络时钟守时模块低阶调整方法 |
CN103117742A (zh) * | 2011-11-17 | 2013-05-22 | 沈阳工业大学 | Gps/北斗双模卫星时钟晶体振荡器驯服*** |
CN103427793A (zh) * | 2013-07-29 | 2013-12-04 | 国电南瑞科技股份有限公司 | 基于温度补偿的对时守时***及方法 |
CN104570717A (zh) * | 2013-10-25 | 2015-04-29 | 沈阳工业大学 | 基于gps/北斗和有限状态机的守时*** |
-
2016
- 2016-05-16 CN CN201610324056.1A patent/CN106026919B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079687A (zh) * | 2006-05-25 | 2007-11-28 | 上海欣泰通信技术有限公司 | 一种基于最小二乘法模型的时钟调整算法 |
CN103117742A (zh) * | 2011-11-17 | 2013-05-22 | 沈阳工业大学 | Gps/北斗双模卫星时钟晶体振荡器驯服*** |
CN102629103A (zh) * | 2012-03-30 | 2012-08-08 | 上海鸿晔电子科技有限公司 | 卫星授时***的网络时钟守时模块低阶调整方法 |
CN103427793A (zh) * | 2013-07-29 | 2013-12-04 | 国电南瑞科技股份有限公司 | 基于温度补偿的对时守时***及方法 |
CN104570717A (zh) * | 2013-10-25 | 2015-04-29 | 沈阳工业大学 | 基于gps/北斗和有限状态机的守时*** |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106706021B (zh) * | 2016-11-15 | 2019-10-22 | 中国电子科技集团公司第四十一研究所 | 一种基于调制域分析仪的线性调频测量方法 |
CN106706021A (zh) * | 2016-11-15 | 2017-05-24 | 中国电子科技集团公司第四十研究所 | 一种基于调制域分析仪的线性调频测量方法 |
CN107276582A (zh) * | 2017-05-17 | 2017-10-20 | 电子科技大学 | 一种基于模拟电路的温度补偿晶体振荡器 |
CN109100930A (zh) * | 2017-06-21 | 2018-12-28 | 卡西欧计算机株式会社 | 电子设备、日期时间取得控制方法及记录介质 |
CN107621644A (zh) * | 2017-08-09 | 2018-01-23 | 上海移为通信技术股份有限公司 | 一种gps定位***的抗干扰控制方法 |
CN110412438A (zh) * | 2018-04-26 | 2019-11-05 | 三星电子株式会社 | 生成包括晶体振荡器的测试产品的晶体模型的方法和*** |
CN108900161A (zh) * | 2018-05-21 | 2018-11-27 | 阮晓萍 | 一种晶振器时钟温度误差补偿方法 |
WO2019232610A1 (en) * | 2018-06-06 | 2019-12-12 | Microsemi Semiconductor Ulc | Non-linear oven-controlled crystal oscillator compensation circuit |
CN112272921B (zh) * | 2018-06-06 | 2024-04-23 | 美高森美半导体无限责任公司 | 非线性恒温晶体振荡器补偿电路 |
CN112272921A (zh) * | 2018-06-06 | 2021-01-26 | 美高森美半导体无限责任公司 | 非线性恒温晶体振荡器补偿电路 |
CN110740001A (zh) * | 2018-07-18 | 2020-01-31 | 北京松果电子有限公司 | 晶振频率校准方法、装置以及电子设备 |
CN110262210A (zh) * | 2019-06-28 | 2019-09-20 | 北斗天汇(北京)科技有限公司 | 基于计数器的晶振守时方法 |
CN110262210B (zh) * | 2019-06-28 | 2021-03-26 | 北斗天汇(北京)科技有限公司 | 基于计数器的晶振守时方法 |
CN111711446A (zh) * | 2020-06-10 | 2020-09-25 | 国网湖南省电力有限公司 | 一种利用gps信号驯服晶振频率的方法、***及介质 |
CN111711446B (zh) * | 2020-06-10 | 2023-08-15 | 国网湖南省电力有限公司 | 一种利用gps信号驯服晶振频率的方法、***及介质 |
CN112213541A (zh) * | 2020-09-08 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | 一种应用于暂态录波型故障指示器的守时方法 |
CN112485520A (zh) * | 2020-12-03 | 2021-03-12 | 成都市精准时空科技有限公司 | 基于电压采样的绝对频差测量方法及***及装置及介质 |
CN112485520B (zh) * | 2020-12-03 | 2024-03-22 | 成都市精准时空科技有限公司 | 基于电压采样的绝对频差测量方法及***及装置及介质 |
CN113067578A (zh) * | 2021-03-17 | 2021-07-02 | 成都金诺信高科技有限公司 | 一种时钟源的精密补偿方法 |
CN113422593A (zh) * | 2021-07-05 | 2021-09-21 | 北京信息科技大学 | 滤波方法、滤波器、计算机可读存储介质、处理器和fpga |
CN113422593B (zh) * | 2021-07-05 | 2024-04-26 | 北京信息科技大学 | 滤波方法、滤波器、计算机可读存储介质、处理器和fpga |
CN114362140A (zh) * | 2021-12-01 | 2022-04-15 | 国网内蒙古东部电力有限公司通辽供电公司 | 一种适用于配电网量测装置的高精度守时方法及装置 |
CN114362140B (zh) * | 2021-12-01 | 2023-11-21 | 国网内蒙古东部电力有限公司通辽供电公司 | 一种适用于配电网量测装置的高精度守时方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106026919B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106026919A (zh) | 高精度晶体振荡器的守时补偿方法 | |
CN1904642B (zh) | 用于补偿用作采样频率的本地时钟的漂移的装置和方法 | |
CN103970008B (zh) | 一种基于晶振误差补偿的守时方法 | |
CN103472259B (zh) | 一种硅微谐振式加速度计温度补偿方法 | |
CN103033186A (zh) | 一种用于水下滑翔器的高精度组合导航定位方法 | |
CN104702214B (zh) | 一种晶体振荡器频率补偿的方法 | |
CN103427793B (zh) | 基于温度补偿的对时守时***及方法 | |
CN106932125B (zh) | 一种硅谐振压力传感器的补偿方法 | |
CN103699001B (zh) | 利用恒温晶体振荡器实现的计时方法及*** | |
CN101562451A (zh) | 二级频标的精密驯服保持方法 | |
CN111650617A (zh) | 基于新息加权自适应不敏卡尔曼滤波的晶振频率驯服方法、***及介质 | |
CN105324647A (zh) | 用于标定或调整任意可振荡单元的方法 | |
US8224606B2 (en) | Measuring clock jitter | |
CN114584136A (zh) | 基于gnss的高稳晶振驯服与保持***和方法 | |
CN104142624B (zh) | 一种基于波形匹配的时间同步方法及*** | |
US20110226056A1 (en) | Method for simulating the operating behavior of a coriolis gyro | |
CN104199054A (zh) | 一种用于北斗卫星导航***共视数据的预处理方法 | |
CN115524725A (zh) | 基于差分观测与回归方程模型修正北斗和晶振秒时钟的方法 | |
CN112639630B (zh) | 用于估计用于对地震数据样本进行测年的时钟漂移的方法和*** | |
JP6770156B2 (ja) | 電子ウオッチに組み込まれた時間基準の平均周波数を調節するための方法 | |
Elster et al. | Uncertainty evaluation of dynamic measurements in line with GUM | |
CN103727947B (zh) | 基于ukf滤波的bds与gis深耦合定位方法和*** | |
CN102679984B (zh) | 基于极小化矢量距离准则的有限模型滤波方法 | |
CN104316047B (zh) | 一种利用gps自主提高敏感器数据时标精度的方法 | |
CN110618604B (zh) | 一种利用ntp辅助作源提高守时精度的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |