时间和时钟3.docVIP

关于时间服务应用的简述 频率源应用简述： 频率源按照应用的特征来分，大体可分为两种类型：一种做时间服务，一种做时钟服务。两种应用的特性作如下比较： 应用 功能 工作频段 表现形式 典型器件 频率稳定性 特征 表征的指标 时间服务 同步 基带 时间戳+时间信息单位 原子钟 漂移特性 单调 长稳，中稳 时钟服务 调制解调 射频 时钟信号、载波 低相噪OCXO 抖动特性 高斯 短稳，相噪 作为时间应用，通常由于原子钟价格比较昂贵，所以比较实际的应用是采用GPS/北斗+OCXO来实现空口同步。在此应用方案中关键指标是： 1 GPS/北斗接受模块本省输出的PPS的精度和稳定度，通常用输出PPS的时长值正态分布的方差α来表示，一般在15ns-100ns之间。 2 需要关注 OCXO 的变化特性，亦即在守时（丢失频标）过程中，频率源可能发生的频率变化。 通常守时能力可以表达为：⊿T=T0+∫0T ⊿fdt,其中T0为初⊿f相差（时差），(0,T)为积分时长（即丢失频标时长）, ⊿f为OCXO 在丢失频标这段时间里频率变化特性函数。这个函数通常包含频率相对于时间，温度，加速度，辐射场，噪声等变化特性，如下图所示： 在实际用过程中，主要考虑时间很温度的影响，对于时间服务器设计的一系列算法也是围绕这些特性设计的。 系统结构 程序过程： 1 开机上电 系统自上电起，有预热过程，此时OCXO 的恒温槽处于加温状态，频率在逐渐上升，一般2-3分钟后OVEN 工作进入稳定状态（电流下降到正常工作水平），在这一过程中，OCXO 频率变化有如下特性 如上图可见，前2小时里OCXO 的频率变化仍然剧烈，此阶段如果丢失频标，建议守时的频率标签采用前一次稳定工作时获得的频率数据。 2 计数 10mhz的OCXO进入FPGA之后建议先倍频到100mhz。 FPGA中的频率计数时长建议在（1024s，4096s）之间取决于系统设计指标。注意COUNTER的计数务必是连续的，不是分段加总。如1024s计数到100mhz的脉冲个数应该是：124000000000， 如果用二进制数表示大概要37位。 3 频率标签计算 对于高精度的要求，比如要守时到24小时1.5us以内，或者系统计算的频率精度到2.5*10-11以上时，需要考虑对COUNTER直接获得的频率标签进行修正。如下示意图： F=f2+(f2-f1)/3 对于低精度守时要求比如1小时守时1us或者系统的频率精度在1*10-10以下的，可以不考虑校正COUNTER 中的读数。 4 PPS 输出（跟踪，守时） 输出端的设计主要是要考虑输出PPS的精确度和稳定度，分为两个阶段考虑，在跟踪阶段，此时频标仍然存在，但是输出的PPS 如果直接参照(或者依据判别)从GPS模块获得的PPS的沿的话，那么GPS 模块输出抖动将会被全部引入到输出端，这样系统输出的也会出现100NS的高斯分布的抖动。建议沿的判别依据OCXO晶振周期的沿。通常OCXO的沿抖动特性会非常好（一般可以到几个ps），但是由OCXO合成的PPS 可能和标准的PPS的脉宽有一定差值。这个差值可以由COUNTER中的读书获得，如 COUNNTER 读到1024s的晶振周期总数为 1024X100000000+4（100mhz的计数频率）.这个时候如果再用100000000个晶振周期代表1pps的话，那么每1024s和UTC会相差4*10ns=40ns。若要校准这个误差，需要在1024s里平均插入40ns，即没256s插入一个晶振周期。 5 补偿及高阶算法 高阶算法的核心思想，是根据晶体的变化特性，推演出在丢失频标这段时间里由于频率变化引起的累积时差，然后做一个补偿。目前这个补偿基本上是针对频率随时间的变化特性（老化率），对于温度的变化特性基本上还没有成熟的算法可以有效实现。 频率对于时间的变化特性可以由如下公式体现： F=Aln(Bt+1)+c (其中F 表示频率，t表示时间， A,B,C是系数) 如下图所示 其中A,B,C 系数随每一个晶振都会不同，所以不能用一个固定的系数在所有的晶振上。 要获得每一个晶振具体的A,B,C的系数，需要正对晶振的特性做回归分析。由于直接做对数回归比较困难 建议将对数回归转化为线性回归，可以做如下转化： 若BT1 则 T=ln（Bt）=lnB+lnt ? F=AlnB+Alnt+C=Alnt+AlnB+C ? 令T=lnt 则 F=AT+AlnB+C 此时做线性回归即可得出 A, B,C的值。 做线性回归可以采用最小二乘法。（不过一般的高级语言的开发工具里应该带有回归的库函数） 算法注意点： 以下注意点1-7都是基本优化没有对OCXO频率的变化做补偿。首先应根据设计需求做一个理论上的量化推演，以变电站主钟为例，假定设计指