"计算机测量与控制5+1+!5!-(!)"&'()
*+,-.,/0*-,(,1+2&'1+-'3
"#+30"#
收稿日期 +1++1+1>+"修回日期
+1++1+1.$基金项目 河北省科技重大专项(+++2131-Y
)$作者简介 孙同川(-332)
!男!河北石家庄人!硕士研究生!主要从事卫星导航时间频率方向的研究$引用格式 孙同川!王振岭!孙建设!等5基于]A 9=A &滤波的原子时算法研究,M -5计算机测量与控制!+1+!!!-(!)"+30+335
文章编号 ->.-0
,32 +1+! 1!1+301>""[R O  -15->,+> 658&;%5--a 0.>+ *#5+1+!51!510!""中图分类号 C -+.5-h +""文献标识码 U 基于C /3(/1滤波的原子时算法研究
孙同川- 王振岭- 孙建设+ 刘铁强-(-/中国电子科技集团第,0研究所!石家庄"1,11,-+
+/中国人民解放军>-.--部队!新疆喀什"200111
)摘要 守时系统的目标在于建立和保持一个稳定可靠的时间尺度!时间尺度算法正是基于此目标计算出一个频率稳定度'准确度'可靠性更高的时间尺度!时间尺度的算法本质就是综合守时系统内的原子钟!通过各原子钟与主钟的Ka -组观测钟差对K 台原子钟的权重和预测值进行估计+传统的加权平均算法会忽略发挥主要影响的噪声过程!更注重权重的合理分配来提高综合
原子时的稳定度!缺少对噪声的关注!针对守时系统实时性的需求!对原子钟噪声模型进行了研究!在频率预测过程中研究了]A 9=A &滤波和频率跳变检测的应用!并与传统加权平均算法进行了对比!仿真实验表明改进的算法提升了综合原子时的中长期
稳定度!其中-11天稳达到了,f -1a -0
数量级!既保留了U P -时间尺度连续'实时的良好特性!又避免了]A 9=A &算法发散性的问题!经实际测试可应用于小型守时实验室的守时系统构建$
关键词 守时系统+时间尺度算法+]A 9=A &滤波+预测值+频率跳变检测
6,0,/-8<'15+'(:8?:(,53;
'-:+<(@/0,9'1C /3(/1T :3+,-"^K P )&'8V D A &-!XU K ZY V F &9%&'-!"^KM %A &7V F +!H O ^Y %F c %A &'
-(-/,0*VT F 7F A $8V O &7*%*D *F )IN Q P N !"V %6%A <V D A &'
"1,11,-!N V %&A ++/N V %&AU 8A J F =:)IU F $)7#A 8F A &JU F $)J :
&A =%87!]A 7V %"200111!N V %&A )5B 0+-/8+"P V F ')A 9)I *%=F;F F #%&'7:7*F =%7*)F 7*A @9%7VA &J=A %&*A %&A 7*A @9F A &J $F 9%A @9F *%=F 78A 9F 5*V F *%=F 78A 9F A 9'
)$%*V =@A 7F J)&*V %7')A 9%7D 7F J *)8A 98D 9A *F *V F *%=F 78A 9F4%*VV %'V F $I $F c D F &8:7*A @%9%*:!A 88D $A 8:A &J $F 9%A @%9%*:5P V F F 77F &8F )I *V F *%=F 78A 9F A 9')$%*V =%7*)%&*F '$A *F *V F #$%=A $:89)8;7%&*V F *%=F ;F F #%&'7:7*F =!A &J F 7*%=A *F *V F4F %'V *A &J #$F J %8*F J E A 9D F )IKA *)=%889)8;7*V $)D 'V *V FKa -'$)D #7)I )@7F $E A *%)&89)8;J %I I F $F &8F 7@F *4F F &F A 8VA *)=%889)8;A &J=A %&89)8;5P V F *$A J %*%)&A 94F %'V *F J A E F $A 'F A 9')$%*V =%'&)$F 7*V F #$)8F 77)I #$)J D 8%&'*V F=A %&%&I 9D F &8F )&&)%7F !A &J #A :7=)$F A **F &*%)&*)*V F $F A 7)&A @9F J %7*$%@D *%)&)I4F %'V *7*)%=#$)E F *V F 7*A @%9%*:)I *V F 7:&*V F *%8A *)=%8*%=F 5O *9A 8;7A **F &*%)&*)&)%7F 5U %=%&'A **V F $F A 9B *%=F $F c
D %$F =F &*7)I *V F *%=F;F F #%&'7:7*F =!*V F A *)=%889)8;&)%7F=)J F 9%77*D J %F J !A &J*V FA ##9%8A *%)&)I]A 9=A &I %9*F $A &JI $F c D F &8:6D =#J F *F 8*%)&%77*D J %F J %&*V FI $F c D F &8:#
$F J %8*%)&#$)8F 77!A &J8)=#A $F J 4%*V*V F*$A J %*%)&A 94F %'V *F JA E F $A 'FA 9')$%*V =!*V F7%=D 9A *%)&$F 7D 9*77V )4*V A **V F %=#$)E F JA 9')$%*V =%=#$)E F 7*V F7*A @%9%*:A &J$F 9%A @%9%*:)I *V F7:&*V F *%8A *)=%8*%=F %&*V F=F J %D =A &J9)&'*
F $=!4V %8V $F A 8V F 7,f -1a -0)$J F $)I=A '&%*D J F %&-11J A :75O *&)*)&9:
$F *A %&7*V F '))J 8V A $A 8*F $%7*%87)I 8)&*%&D )D 7A &J $F A 9B *%=F %&U P -*%=F 78A 9F !@D *A 97)A E )%J 7*V F J %E F $'F &8F )I *V F]A 9=A &A 9')$%*V =5(:*V F #$A 8*%8A 9*F 7*%&'!%*8A &@F A ##9%F J *)*V F 7:
7*F =8)&7*$D 8B *%)&)I 7=A 99*%=F;F F #%&'
9A @)$A *)$%F 75C ,D
E '-90"*%=F;
F F #%&'7:7*F =+*%=F 78A 9F A 9')$%*V =+]A 9=A &I %9*F $%&'+F 7*%=A *F +I $F c D F &8:V )##%&'J F *F 8*%)&F "引言
时间对一个国家的经济'社会生活'军事等方面有着十分重要的作用!时间的稳定关系着国家和社会的安全稳定!时间频率体系建设作为国家的重大基础设施!具有非常重要的意义$卫星导航'通信互联'联合作战'反导拦截'金融'电力等对时间频率系统的准确度和可靠性的需求越来越高$只有一个高精度的时频体系!才能提供可靠
的时间频率服务$随着当前国产原子钟技术水平的突破,--!
为了提高自主可控能力!采用国产原子钟组进行原子时算法的研究具有重要意义$通常情况下!构建时间频率体系的关键是建立和保持系统时间!时间尺度的性能一方面受
时间频率系统内硬件(
比如原子钟)的性能水平的影响!一方面受所采用的时间尺度算法的影响,+-$
较为著名的时候尺度方法主要包括"目前国外普遍使
用的U H Z R "算法,!-'U P -算法,0-等!+1世纪21年代!针
苹果 范冰冰浴室戏
对传统加权平均算法采用U 99A &方差计算权重会忽视主要噪声过程之外的噪声过程的问题!美国学者(A $&F 7将]A 9=A &
滤波用于钟组噪声的处理提出了]A 9=A &算法,,>-!近年来!
国家授时中心对原子钟数据异常情况下如何计算时间尺度进行了研究!并采用改进的]A 9=A &滤波计算时间尺度取得
了比较好的结果,-0-!国家授时中心还对W )&J $A ;B N F #
F ;滤波进行研究!并研究分析这种方法在时间尺度计算中的应
"
投稿网址 444567689:
;<58)=Copyright ©博看网. All Rights Reserved.
第!期孙同川!等"基于]A9=A&滤波的原子时算法研究#+3,
"#用!这种方法计算出的综合时间尺度稳定度进一步提高,-2-$
U H Z R"算法和U P-算法的核心思想都是加权平均!具体
做法是根据某个计算周期内钟组内各个原子钟的稳定度来
确定每台原子钟在此计算周期内的权重!通关比较合理的
权重分配!使钟组内的噪声最小!以此来提高时间尺度的
稳定度$但是加权平均算法只考虑起主要影响的噪声过程!
而不论其他噪声过程如何!在计算时间尺度时并没有消除
或抑制钟组内的噪声!相反]A9=A&算法通过]A9=A&滤波
器对原子钟噪声建模消除或抑制钟组内的噪声$从文献
,.-可以知悉]A9=A&'U P-算法的短期稳定度较好!而
U H Z R"的长期稳定度较好!U P-算法和]A9=A&算法是根
据当前计算周期的钟差数据来预测下个计算周期的钟差!
具有实时性的特点!而U H Z R"算法则是采用过去一段时间
的钟差数据来计算时间尺度,2-!具有滞后性!造成滞后的
原因主要是因为U H R Z"算法的权重计算和速率预报需要采
用数据计算前一个月的钟差预测值$]A9=A&算法的稳定度
较好!但是存在发散性的问题,2-$
守时实验室的目标是产生一个稳定'准确'可靠的时
间尺度!在二级守时节点中!守时钟组类型单一(大多为
铯钟组)!针对此前我国大多数实验室采用的原子钟大多为
进口原子钟的问题!为了提高自主可控水平!经典的U H B
Z R"算法主要考虑原子时的长期稳定度!在二级守时节点
中综合原子时应同时注重综合原子时的短期稳定度和灵活
性!以便守时系统随时进行调整$
因此对加权平均算法进行了研究和改进!改进后的算
法将针对传统加权平均算法中的噪声问题!加入]A9=A&滤
波过程用于改善频率预测值并加入了频率跳变检测!可以
进一步抑制原子钟的噪声!提高综合原子时的稳定度$首
先研究分析加权平均算法的基本原理!然后对算法进行改
进!最后通过搭建的试验平台对改进的算法性能进行验证
分析$
G"原子钟噪声模型的研究
时间尺度算法目的就是研究原子钟间的噪声问题!通
过设定符合实际钟组情况的参数使得钟组内的噪声最小!
以此提高时间尺度的稳定度$而原子钟作为一种表征频率
的设备!其纸面读数与绝对时间之间必然存在一定的偏差!
同时原子钟的纸面读数还受到观测噪声的影响!因此为了
分析原子钟的噪声!首先要对其建立合适的噪声模型$
原子钟的观测方程可以表示为,+1-"
<(!7)%F-(!7)+!#1(7)(-)
式中!<(!
7
)为观测量+!#1(7)为观测噪声!1是标准
高斯白噪声!属于X C?!!表示噪声强度$
原子钟模型可用随机微分方程,23-描述"
J F-(!)%F+(!)#J!+!-#J W-(!)
J F+(!)%J#J!+!+#J W+(!
4)(+)
式中!F
-(!)!F
+
(!)表示原子钟的两个状态变量!F
-
示时差状态变量!F
+表示频率偏差的随机游走噪声!W
-
(!)!W
+
(!)是两个独立的维纳过程$为了方便分析仿真
实验时将此微分方程的初值设置为1$M/U/(A$&F7等人提
出原子钟噪声的幂律谱模型即原子钟噪声是五种噪声分量
线性叠加的结果!文献,---不仅描述了各个噪声特性还
给出了仿真方法!它的谱密度函数"
J O(3)%,+
-%=+
/-3-(!)
""不同的-取值对应不同类型的噪声!具体为,-.-"-e
a+是随机游走调频噪声(TX\?)!-ea-是闪变调频噪
声(\\?)!-e1是调频白噪声(X\?)!-e-是闪变调
相噪声(\C?)!-e+是调相白噪声(X C?)!不同原子时
之间的算法受到各种噪声不同程度的干扰!但这些噪声有
时不能避免!有时可以采用加入滤波器的方式进行消除$
U P-算法中一般都采用了指数滤波器!指数滤波器可以比
较有效的控制原子钟的调频白噪声和频率随机游走噪声$
本文算法采用的滤波器是卡尔曼形式的!卡尔曼滤波器是
最小均方意义下的最优估计!可以有效避免观测量中频差
随机游走噪声的影响$
加权平均算法在计算过程中用U99A&方差表征一台原子
钟的稳定度!U99A&方差主要与原子钟的噪声有关!每台钟
的权重大小主要与某个起主要作用的噪声过程有关!针对
加权平均算法中的权重和噪声过程问题!]A9=A&算法从估
值的角度出发!避免了权重的分配问题!对所选取的主钟
和标准时间的差做最小均方差意义下的最优估计$这种方
法通过对每台原子钟建立一个状态方程评估的噪声模型!
来构成]A9=A&滤波器$
假设原子钟的噪声模型符合"
N A(!)
O A(!
,-)%-,
,
-
1-
N A(!=-)
O A(!=-
,-)+1A(!=-)
.A(!=-
,-)(0)
*A%
C1
1[
,-A A%-!3!*(,)
式中!N
A
(!)'O
A
(!)为第!次测量时的时间和频率与理想
时间尺度的偏差+*
A
是系统噪声矩阵+.
A
为随机游走调频
噪声!其方差是/
A
+1
A
为白调频噪声!其方差是C
A
$
由]A9=A&滤波原理!得到钟组的动态模型为"
'(!)% '(!=-)+B(!)(>)
式中! 为转移矩阵!B(!)为系统噪声矩阵$
钟组的观测方程为"
.(!)%/F(!)+,(!)(.)
/
%
-=-131
-1=-31
88888
-113=
.
/
1
-(*=-)Q*
(2)
式中!/是测量矩阵!,(!)为测量噪声
N-(!+%)
O-(!+%)
N+(!+%)
O+(!+%)
8
N*(!+%)
O*(!+%
.
/
1
)
%
"
投稿网址 444567689:;<58)=
Copyright©博看网. All Rights Reserved.
""计算机测量与控制"第!-卷
#+3>"#-%-.+%
+
33331-%333311-3333333=
333
3333-%-.+%
+
3
3331-%
3
3
3
3
.
/011
1
-N -(!)O -(!)N +(!)O +(!)8N *(!)O *(!.
/01)+1
-
.-1
+.+
8
1*.
./01*(3
)""根据卡尔曼滤波模型和噪声矩阵!通过卡尔曼滤波迭
代计算得到每台钟的状态!具体迭代步骤如下"
0'(7X 7=-)%#(7)'(7=-X 7=-))(7X 7=-)%#(7))(7=->7=-)#,
(
7)+*(7)1%)(7>7=-)/,,/)(7>7=-)/,++-=
-0'(7X 7)%0'(7>7=-)+1,.(7)=/0'(7>7=--)(7.7)%(-=1/))(7.7=-
'(
)(-1)式中!)为'的误差方差矩阵+*为系统驱动噪声矩阵+1为卡尔曼增益!+为测量噪声矩阵$
H "综合原子时算法研究与改进
目前守时实验室常用的综合原子时算法主要有两大类"加权平均算法和各种滤波类算法
,-,-
$时间尺度产生最典型
的例子就是世界协调时(^P N )的产生!国际计量局((O B
C ?)采用U
H Z R "算法对全球的守时实验室数据加权平均!以此提高时间尺度的可靠性和稳定性$上文提到加权平均算法的基本原理都是通过给每台原子钟分配一个权重来调
整钟组内的噪声关系,
-
>-!以此来提高时间尺度的稳定性$但U H Z R "算法采用连续一个月的数据来预测前一个月的钟差值!具有滞后性参数调整设置极其不方便!需要的钟差数据时间太长!不适合在规模较小的守时实验室使用$
理想状态下!时间尺度算法的最主要步骤是根据钟组内*台原子钟!利用*a -组钟差对钟组内的原子钟进行权重调整!以及钟差预报等$而时间尺度算法就是通过计算每台钟的权重来调整钟组内的噪声关系使其对钟组产生的时间尺度影响最小化!由此生成的时间尺度比钟组中任何一个单独的原子钟具有更高的可靠性'稳定性和频率精度$钟组产生的时间尺度是采用数学方法计算得到的!并不是某台钟的时间也不是他们简单的集合!这就是综合原子时!它的物理实现是通过对钟组中任意一个原子钟进行适当的校正来实现的!如果没有测量噪声的影响!该值与所采用的原子钟无关$时间尺度算法的输入是每个原子钟与选取的主钟之间的钟差!原子钟作为一个物理系统产生一个频率!原子钟的时间是由频率经过人工推导得来的!频率才是真正的物理量!因此用来衡量原子钟性能的参数都包含了对频率的描述!原子时算法中还需要估计每个原子钟频率偏移参数$
类加权平均算法在计算时间尺度的过程上大致相同!U H R Z R "采用之前的钟差数据进行计算!U P -算法则是根
据当前的钟差值估计下一个时间的钟差值$在钟组内选定一台钟作为主钟!主钟的物理信号输入相位微跃器!微跃器可以调整主钟的物理信号相位!微跃器和主钟的组合称为组合钟!每台原子钟和组合钟的钟差值
都可以通过与主钟的钟差间接计算得出$物理主钟信号与相位微跃器调整后的主钟信号的钟差值可通过钟%与主钟的钟差间接计算得
出!每台钟的钟差估值为N 6A (!h %)!对钟差估值N 6A (!h %)加权平均计算出主钟与组合钟的钟差N 6(!h %)!对相位微跃器则参照这个值对主钟物理信号进行调整$
假设钟A 在时刻!的钟差和速率分别为N A (!)!O A (!)!测量时间间隔为%$则钟组内某台钟!h %时刻相对于组合钟的钟差可由下式计算"
N A (!+%)%N A (!)+O A (!+%)#%(--)""将计数器在时刻!h %测量得到的钟A 和主钟的钟差记为!A (
!h %)!则主钟与组合钟的钟差可由下式计算"N 6A (!+%)%N A (!+%)=!A (
!+%)(-+)""即"
N 6A (!+%)%N A (!)+O A (!)#%=!A (
!+%)(-!)""这样就通过钟A 的钟差和速率计算了组合钟相对主钟的钟差值$
假设钟组内有*台中!那么可以得到*a -个钟差值!在对其进行加权平均可以得到"
N A (
!+%)%,*=-
A %-
,
,A
#N 6A (!+%)-(-0
))式中!N 6A (!h %)为第A 台钟间接计算出的钟差估算值!,A
为它的权重!N A (!h %)是钟组内所有原子钟计算的得出的钟差值的加权平均!N A (
!h %)是相位微跃器调整主钟输出信号的相位的参照!由此便生成了综合原子时的物理信号$为了削弱噪声对时间尺度稳定性的影响!U P -算法需要引入一个指数滤波器!指数滤波器的计算方式为"
O A (
!+%)%3A (!+%)+7#O A (!)-+7
(-,
)式中!7是一个重要参数!这里取7e !1
!第A 台钟的频率估计值为3A !通过时间差分计算得到"
3A (
!+%)%N A (!+%)=N A (!)%(->)""由上述计算过程可以看出"U
P -算法钟每个原子钟有两个状态!原子钟相对主钟的时间偏移和频率偏移!频率偏移可以用来计算!但算法本身不能对频率偏移进行估计$改进的算法结合传统的U P -算法和]A 9=A &滤波来估计表示频率随机游走噪声和原子钟的频率偏移状态和方差$对于一个给定的原子钟!这种状态不是一种物理状态!而是在白噪声调制频率的情况下对频率偏移的数学估计!假设这种状态称之为%`&
!那么这个状态的方差就表示对偏移的数学估计的置信度$
将U P -算法中的钟A 的钟差写为"
N A (!+%)%N A (
!
)+O A (!)+9A %+
()
%(-.
)式中!9A 是频率漂移率$
"
投稿网址 444567689:
;<58)=Copyright ©博看网. All Rights Reserved.
第!期孙同川!等"基于]A 9=A &
滤波的原子时算法研究#+3."#
如果钟:在!h %时刻的刻度钟差测量值是原子钟N A :(!h %)!则钟:在!h %时刻的钟差预报偏差为"
N :(
!
+%)%,;
A %-
\A
,
N A
(!+%)=N A :(!+%)-(-2
)""原子钟的权重计算方式通过时间残差计算"
,A %++N (%)++A (
%)(-3)式中!++N (%)为钟组时间预测残差!++
A (%)为钟A 的时间预测残差$采用这种方式计算权重!是在最小方差意义下优化了钟组时间的稳定性$
频率偏移量可通过下式计算"
O
A (!+%)%F A (!+%)=F A (!)%
(+1)""把频率偏移量合并到钟A 当前平均频率偏移的指数滤波估计中"
O A (!+%)%-K A +
-,O A (!+%)+K A O A (!)-(+-)""K A 是钟A 的指数频率加权常数!K A 由钟A 的调频白噪
声和随机闪烁调频白噪声决定!用来衡量时间预测值的稳定性$
其中"
K A %-++=-+-!+0%+K A ;-.
+!%,
-,-
+
1
(++)""%K A ;是!O (%)到达最小值的时间!%1是用于计算!O (%)的最小时间$
接下来需要对状态G 进行频率偏移估计!文献,-+-和,-!-对原子钟的频率稳定度和钟差预测不确定度进行了估计!采用一个简单的卡尔曼形式对频率偏移量进行过
滤,+--!假设B (,)是测量噪声!对状态G 建模!状态G
包含随机游走噪声和固定漂移!那么G 不是一个由原子钟
直接产生的物理层面的频率!而是过滤了白调频噪声的之后的频率随机游走分量加上频率漂移$
将式(0)和(,
)代入卡尔曼方程,0-
!系统模型为"O A (!+%)%O A (!)+9A
%+.(%)(+!)""测量模型"
O A %O A (!)++A (
%)(+0
)式中!O A 是频率直接测量值!+A 是白噪声$
残差估计"
#(!+%)%#(!)+!+.#%(+,
)""其中"
#%++
A (%)##++A
(%)
+#(+>
)""!+.
是白噪声方差!++
A 是测量噪声!#(!)是G 的残差!#(!h %)是残差#(
!)的预测值$这相当于一个自适应卡尔曼滤波器!G 的残差可以从初始值进行变化!G 的指数滤波器参数也可随时间变化$
方程的稳态形式中"
#%.
+
+
-+0#++A (%).
+
=()
-(+.
)""此原子钟的指数频率加权常数为"
K %++(%)#+.
+
(+2
)""频率预测为"
O (!+%)%++A #
O (!)+##?O ++A (
%)+#(+3
)""由此可以看出卡尔曼形式也得到了类似U
P -中的指数滤波器形式,-3-!此方法继承了U P -算法对闪变频率的建模能力$
频率跳变的检测需要观测检查一段时间内的数据!考虑到实时性的操作!可以将某一区间内的平均频率偏
移估计值与该区间开始时的滤波估计值进行比较$将频率偏差的估计值作为离值进行测试!就可以确定频率跳变是否在最近发生$
通过测量得到的是时钟之间的相位差!由于受到频率闪烁噪声和随机游走噪声的影响!频率跳变会在某一时刻发生!有些情况下时钟本身小的频率波动会被当成随机噪声!为了区分这种现象!需要对频率跳变的大小进行限定$
本文的频率跳变检测方法是通过迭代每个时钟的测量范围来检测频率跳变!从当前测量时间到之后的某个时间!此时间间隔为%=%&!$=A L 是测量次数的最大值!$=A L 是由%=%&
确定的整数$在每次测量过程中!对每台时钟的频率跳变进行检测!定义一个时钟在测量时间间隔内的平均频率O A E '"
O A E '%
N =-=N =$
!=-=!=
$(!1
)""N a $是在!a $时刻的估计值!N a -是在!a -时刻的估计值$
通过将这个平均频率和当前测量时刻的估计值进行比较!
实际就是把固定时间区间上的实时频率估计值与该区间的平均频率估计值进行比较!如果两者频率差大于阈值!就可以认为时钟发生了跳变!具体判断方法为"
X O A E '=O
=$X %0
%=%&$
#=$+#=+++!+-(
)
N +!+
/
N #槡$(!-
)式中!!+-N 和!+
/
N 是白噪声和随机游走调频噪声$通过这种方式检测频率跳变的发生!当检测到某台时钟发生频率跳变时!此时算法会将发生频率跳变的时钟从这个计算周期中剔除!也就是把发生频率跳变的时钟的权
重变为1!此后的计算周期中对该时钟的频率继续观测!直到算法自适应学习了此时钟的新频率特性!然后将此时钟重新加入到原子时的计算当中$这种方法的好处在于原子钟发生频率跳变时!此台钟对综合原子时的影响会降低!有利于提高综合原子时的长期稳定度$
!"守时系统硬件设计及原理
二级守时节点可分为硬件和软件两个部分!守时系统的硬件部分主要包括"原子钟组'时频信号产生与分配单元'时差测量单元'溯源比对单元'数据存储分析处理和监控单元$硬件部分主要考虑物理信号产生问题'与上级
"
投稿网址 444567689:
;<58)=Copyright ©博看网. All Rights Reserved.
""计算机测量与控制"第!-
#+32"#节点或同级节点间的数据交互问题和可拓展性等+软件部分主要包括"溯源比对软件'综合原子时计算软件'数据采集分析软件和监控软件$稳定可靠的守时系统由硬件和软件共同支撑!缺一不可$守时系统模块图如图-所示
$
图-"硬件系统结构图本次试验的守时系统采用0台国产高性能铯原子钟组成的钟组!铯原子钟独立工作!不对原子钟的数据在频率和相位上进行改变$0台铯原子钟的原始信号直接输出到频率信号切换器和多通道计数器$0台铯钟输出的-C C "信号直接进入多通道计数器!多通道计数器的作用就是对钟差进行测量!测量结果通过^[C 协议端口传输至数据库$综合原子时软件部分每隔一定的周期会从数据库中提取需要的钟差数据进行综合原子时的计算!并将计算结果和每台钟相对于综合原子时的时差存入数据库$
在参与守时的钟组0台铯原子钟里面选用一个稳定度较好的原子钟作为主钟也就是参考钟$0台铯原子钟产生的-1?信号输入频率切换器!频率切换器输出主钟物理信号!随后进入相位微跃器!综合原子时软件对综合原子时进行计算!同时相位微跃器依据综合原子时与主钟的钟差对主钟物理信号的相位进行调整!同时对主钟产生的-1?S <
信号进行驾驭!完成主钟向综合原子时的溯源$#"仿真实验
采用本单位0台铯钟的数据对本文算法进行了仿真验证!铯钟数据得采样周期为>17!采样时间为-+1天数据点-.+211个!计算周期为!1天!软件算法使用N 语言完成!
运行环境是国产麒麟操作系统$钟差测量设备一般使用时
间间隔计数器("T >+1)!"T >+1的观测噪声方差为!+
e -
f -1a +1
$分别采用U P -算法和本文改进算法对综合原子时进
行了计算!并对两种方法的频率偏差值进行了分析+对频率跳变检测进行仿真+对本文算法计算得出的原子时稳定度验证分析$
首先对正常状态的铯钟频率预测值进行仿真分析!改进的方法对频率预测值有明显改善!频率偏差曲线更加平滑!波动更小如图+ !所示$
其中一台铯钟在第.>天铯钟频率发生跳变!图0是铯钟的频率跳变曲线!图,是由本文检测方法计算的平均频率估计值!通过平均频率计算的跳变值与阈值进行比较!可以判断出铯钟频率发生了跳变!因此将发生跳变的铯钟从当前原子时计算周期中剔除也即权重置1
$
图+"U P -
算法频率偏差
图!"
改进算法的频率偏差
图0"跳变铯钟的频率
最后根据0台铯钟组成的钟组系统!分别采用U P -算法和改进后的算法对时间尺度进行了计算!原子时稳定度曲线如图,所示$因为本文算法改善了频率预测值且在一定程度上抑制了频率跳变带来的影响!所以理论上来讲综合原子时的中长期稳定度比U P -算法更好!从图>的曲线可以看出理论与仿真基本符合$
"
投稿网址 444567689:
;<58)=Copyright ©博看网. All Rights Reserved.