CN107219750A - 改善原子钟稳定度的延时线实时噪声修正方法 - Google Patents
改善原子钟稳定度的延时线实时噪声修正方法 Download PDFInfo
- Publication number
- CN107219750A CN107219750A CN201710358764.1A CN201710358764A CN107219750A CN 107219750 A CN107219750 A CN 107219750A CN 201710358764 A CN201710358764 A CN 201710358764A CN 107219750 A CN107219750 A CN 107219750A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- delay line
- atomic clock
- stability
- 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
- 238000002715 modification method Methods 0.000 title claims abstract description 10
- 230000003111 delayed effect Effects 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 10
- 238000012805 post-processing Methods 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04F—TIME-INTERVAL MEASURING
- G04F5/00—Apparatus for producing preselected time intervals for use as timing standards
- G04F5/14—Apparatus for producing preselected time intervals for use as timing standards using atomic clocks
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/08—Details of the phase-locked loop
- H03L7/099—Details of the phase-locked loop concerning mainly the controlled oscillator of the loop
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
一种改善原子钟稳定度的延时线实时噪声修正方法,包括原有的闭环锁定系统的频率信号输出端即功分器的另一输出端设置延时线,该延时线之后连接一个调制器);输出微波信号;计算出误差信号转化成对应的电压;延时一个鉴相周期Tc输出;处理得出误差频率信号输出修正后的相对频率起伏信号y′i‑1。本发明可以改善喷泉钟的稳定度,还获得喷泉钟占空比对稳定度的影响,占空比越大,鉴相能力越强,稳定度越好。
Description
技术领域
本发明涉及原子钟,特别是一种改善原子钟稳定度的延时线实时噪声修正方法,并且可以清晰地观察到喷泉钟占空比对喷泉钟的影响。
背景技术
原子钟应用于时间的计量,可以看作是本机振荡器和量子系统的组合。在锁定环路中,晶振作为参考源为微波综合器提供信号,微波综合器输出的微波信号经由原子钟的量子系统鉴频,原子的跃迁几率差反映了鉴频输出的频率误差信号,再将误差信号转换成电压信号反馈给晶振,实现整个环路的锁定。
在闭环锁定的情况下,原子钟的长期稳定度受限于各种物理效应的影响。对锁定之后的信号进行噪声修正的后处理对于改善原子钟稳定度具有重要意义。
传统的原子钟把探测到的频率误差信号直接反馈给本机振荡器,降低频移的方法只能靠提高物理系统的性能和对原子及其所处环境的更精准的控制,采用后处理的方式对最终输出的信号的噪声的降低,但这些后处理的方法存在时延性,没有一种很好的改善稳定度的实时噪声补偿方法。
发明内容
本发明的目的在于解决上述现有技术的局限,提出一种改善原子钟稳定度的延时线实时噪声修正方法。该方法可以改善喷泉钟的稳定度,还获得喷泉钟占空比对稳定度的影响,占空比越大,鉴相能力越强,稳定度越好。
本发明的技术解决方案如下:
一种改善原子钟稳定度的延时线实时噪声修正方法,该方法包括以下步骤:
①在原有的闭环锁定系统的频率信号输出端即功分器的另一输出端设置延时线,该延时线之后连接一调制器;
②开启原子钟,使原子钟连续运行,本机振荡器输出端输出的微波频率信号经所述的微波综合器、功分器后分为两路输出微波信号;
③一路微波信号经原子钟鉴频后的误差信号输入所述的计算机,由计算机计算出误差信号转化成对应的电压,并通过反馈控制电路给本机振荡器的电压控端,将本机振荡器的输出频率锁定在原子钟的钟频率上,实现闭环锁定,此时所述的微波综合器输出的频率稳定度反映了原子钟的稳定度;
④所述的功分器的另一路输出相对频率起伏信号yi的平均值经所述的延时线延时一个鉴相周期Tc甚至更多输出;其中,是相对频率起伏信号,f0为原子钟的原子跃迁频率,fi为微波综合器输出频率;
⑤经延时线延时后上一周期的相对频率起伏信号平均值和由计算机处理得出误差频率信号同时输入所述的调制器中,所述的调制器按下列公式进行噪声后处理修正,输出修正后的相对频率起伏信号
其中:d是占空比,Tr是原子自由演化时间,TC是原子钟运行周期。
本发明的技术效果如下:
1、本发明由于运用了延时线实时噪声修正方法,使得原子钟的稳定度得到了明显改善。
2、本发明由于运用了延时线实时噪声修正方法,发现原子钟的频率稳定度与占空比的关系,占空比越大,稳定度越好。
附图说明
图1是本发明延时线实时噪声修正法系统框图;
图2是噪声离散化原理图。
图3是在模拟延迟线噪声修正得到的在不同占空比下与未加延时线噪声修正的频率稳定度结果,其中,图3(a)(b)(c)(d)分别为占空比d=0.1,d=0.5,d=0.7,d=1情形下利用延迟线噪声修正法得到的频率稳定度。
图4是延时线实时噪声修正法下稳定度与占空比d=0.3,d=0.6,d=0.9,d=1关系图。其中,(a)为在延时线实时噪声修正法下占空比d=1频率稳定度曲线,(b)为在延时线实时噪声修正法下占空比d=0.6频率稳定度曲线,(c)为在延时线实时噪声修正法下占空比d=0.9频率稳定度曲线,(d)为在延时线实时噪声修正法下占空比d=0.3频率稳定度曲线。
具体实施方式
下面结合附图和实施例,对本发明作进一步说明。
我们进行了数字模拟验证。一般情况下,最终输出信号的稳定度主要受本机振荡器和量子系统的噪声的影响。而本机振荡器所引入的噪声为白噪声,量子系统引入的噪声中量子投影噪声起主要作用,故在进行数字模拟时选取的噪声为白噪声,附加了量子投影噪声。量子投影噪声比白噪声小一个量级。图2是噪声离散化原理图。
在模拟实验当中,用stable32产生数组的大小都是106的两组噪声数组<yi>、<yi-Q>,考虑最简单的情况,<yi>、<yi-Q>的统计特性不随d变化,满足:
σy-MV(τ)=2E-12τ-1/2σQPN(τ)=4E-13τ-1/2
其中,σy-MV(τ)是整个系统中微波(本机振荡器、微波综合器)系统引入的噪声,σQPN(τ)是量子投影噪声,该噪声在整个量子系统引入的噪声中贡献最大。
对整个系统运行时间内的噪声数据做离散化处理,假设每个数据采样时间是1,原子钟的锁定周期为10,则一个锁定周期内共有10数据。对于占空比为d,喷泉钟每个周期鉴频的有效采样个数为10d。
本发明实施延时线噪声修正方法的具体步骤如下:
①该方法在已有的原子钟闭环锁定运行的基础上,对最终输出的频率信号进行修正,请参见图1,图1是本发明延时线实时噪声修正法系统框图,原子钟闭环锁定运行系统包括原子钟1、计算机2、反馈控制电路3、本机振荡器4、微波综合器5、功分器6,且计算机2能够可以实现多路输出,在原有的闭环锁定系统的频率信号输出端即功分器6的另一输出端设置延时线7,该延时线7之后连接一调制器8,开启原子钟1,使原子钟1连续运行;
②本机振荡器4输出端输出的微波频率信号经所述的微波综合器5、功分器6后分为两路;
③由功分器6输出的微波信号一路经原子钟1鉴频后的误差信号输入所述的计算机2,由计算机2计算出误差信号转化成对应的电压,并通过反馈控制电路3输入给本机振荡器4的电压控端,将本机振荡器4的输出频率锁定在原子钟1的钟频率上,实现闭环锁定,此时所述的微波综合器4输出的频率稳定度反映了原子钟1的稳定度;
④所述的功分器6的另一路输出相对频率起伏信号yi的平均值经所述的延时线(7)延时一个鉴相周期Tc甚至更多,输出上一周期的信号平均值其中,是相对频率起伏信号,f0为原子钟的原子跃迁频率,fi为微波综合器输出频率;
⑤该延时后的信号和由计算机(2)处理得出误差信号同时输入到所述的调制器(8)中,按下列公式进行噪声后处理修正输出:
其中:d=Tr/Tc,d是占空比,Tr是原子自由演化时间,Tc是原子钟运行周期;
⑥整个系统持续运行,当原子钟(1)运行一段时间后,关闭整个系统;
⑦对系统运行时间内的经所述的调制器(8)修正后的噪声信号进行阿兰方差统计处理,获得稳定度。
在实际实验中可采用光纤作为延时线:
假设每个周期为1s,则理论上需要延时线的距离L:
L=v×t=2×108×1=2×108m
将输出的信号噪声数据离散化处理,将每个锁定周期Tc均分为N等分,如图2所示,在模拟实验当中,用stable32产生数组的大小都是106的两组噪声数组<yi>、<yi-Q>,考虑最简单的情况,<yi>、<yi-Q>的统计特性不随d变化,满足:
σy-My(τ)=2E-12τ-1/2σQPw(τ)=4E-13τ-1/2
其中,σy-MV(τ)是整个系统中微波(本机振荡器、微波综合器)系统引入的噪声,σQPN(τ)是量子投影噪声,该噪声在整个量子系统引入的噪声中贡献最大。
对整个系统运行时间内的噪声数据做离散化处理,假设每个数据采样时间是1,原子钟的锁定周期为10,则一个锁定周期内共有10数据。对于占空比为d,喷泉钟每个周期鉴频的有效采样个数为10d。
鉴相时间Tr对应其中的M等分,因此占空比d在这里可以表示为d=M/N。
其中:i是代表第几个周期,由于我们是要对所有的在系统运行时间的周期都要作上述的离散化处理,任何周期内公式都一样,为方便起见,我们用i来表示;
⑤该延迟后的信号及由计算机(2)处理得出误差信号输入到调制器(8)中,按下列方式进行噪声后处理修正:
其中:d=Tr/Tc,d是占空比,Tr是原子自由演化时间,Tc是原子钟运行周期,表示在噪声模拟时模拟其他操作引入的附加噪声,如探测时引入的噪声;
⑥整个系统持续运行,当原子钟(1)运行一段时间后,关闭整个系统;
⑦对系统运行时间内的经调制器(8)修正后的噪声信号进行阿兰方差统计处理,获得稳定度。
其中:是噪声修正后的稳定度,是噪声未修正时系统的稳定度,是其他操作引入的噪声造成的稳定度影响。
模拟结果表明延时线实时噪声修正方法确实可以显著改善喷泉钟的频率稳定度。具体的结果如图3、图4所示,并且可以清晰地观察到随着占空比d的增大,稳定度逐渐改善,当d=1时,改善最佳,长稳可以达到量子投影噪声极限。
Claims (1)
1.一种改善原子钟稳定度的延时线实时噪声修正方法,该方法包括以下步骤:
①在原有的闭环锁定系统的频率信号输出端即功分器(6)的另一输出端设置延时线(7),该延时线(7)之后连接一个调制器(8);
②开启原子钟(1),使原子钟(1)连续运行,本机振荡器(4)输出端输出的微波频率信号经所述的微波综合器(5)、功分器(6)后分为两路输出微波信号;
③一路微波信号经原子钟(1)鉴频后的误差信号输入所述的计算机(2),由计算机(2)计算出误差信号转化成对应的电压,并通过反馈控制电路(3)给本机振荡器(4)的电压控端,将本机振荡器(4)的输出频率锁定在原子钟(1)的钟频率上,实现闭环锁定,此时所述的微波综合器(4)输出的频率稳定度反映了原子钟(1)的稳定度;
④所述的功分器(6)的另一路输出相对频率起伏信号yi的平均值经所述的延时线(7)延时一个鉴相周期Tc甚至更多后输出,其中,是相对频率起伏信号,f0为原子钟的原子跃迁频率,fi为微波综合器输出频率;
⑤经延时线(7)延时后输出的上一周期的相对频率起伏信号平均值和由计算机(2)处理得出误差频率信号同时输入所述的调制器(8)中,所述的调制器(8)按下列公式对上个周期的相对频率起伏信号进行噪声后处理修正,输出修正后的相对频率起伏信号y′i-1:
<mrow>
<msubsup>
<mi>y</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>dy</mi>
<mi>i</mi>
<mi>e</mi>
</msubsup>
</mrow>
其中:d=Tr/TC,d是占空比,Tr是原子自由演化时间,TC是原子钟运行周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358764.1A CN107219750B (zh) | 2017-05-19 | 2017-05-19 | 改善原子钟稳定度的延时线实时噪声修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358764.1A CN107219750B (zh) | 2017-05-19 | 2017-05-19 | 改善原子钟稳定度的延时线实时噪声修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107219750A true CN107219750A (zh) | 2017-09-29 |
CN107219750B CN107219750B (zh) | 2019-09-20 |
Family
ID=59945226
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710358764.1A Active CN107219750B (zh) | 2017-05-19 | 2017-05-19 | 改善原子钟稳定度的延时线实时噪声修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219750B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783412A (zh) * | 2017-11-03 | 2018-03-09 | 中国科学院国家授时中心 | 双光子跃迁铷原子钟 |
CN110928173A (zh) * | 2018-09-20 | 2020-03-27 | 中国计量科学研究院 | 减小喷泉钟频率偏移不确定度的调节装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799102A (zh) * | 2012-07-03 | 2012-11-28 | 华中科技大学 | 一种被动型cpt原子钟的温控参数自整定方法及装置 |
CN203164615U (zh) * | 2013-02-16 | 2013-08-28 | 江汉大学 | 脉冲式被动型铷原子钟 |
CN105680857A (zh) * | 2016-01-06 | 2016-06-15 | 中国科学院上海光学精密机械研究所 | 原子钟的实时误差补偿方法 |
EP2897001A3 (en) * | 2014-01-15 | 2016-06-29 | The Boeing Company | Multi-level/multi-threshold/multi-persistency GPS/GNSS atomic clock monitoring |
CN106444343A (zh) * | 2016-11-28 | 2017-02-22 | 中国科学院武汉物理与数学研究所 | 一种微型原子钟数字控制电路装置及方法 |
-
2017
- 2017-05-19 CN CN201710358764.1A patent/CN107219750B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799102A (zh) * | 2012-07-03 | 2012-11-28 | 华中科技大学 | 一种被动型cpt原子钟的温控参数自整定方法及装置 |
CN203164615U (zh) * | 2013-02-16 | 2013-08-28 | 江汉大学 | 脉冲式被动型铷原子钟 |
EP2897001A3 (en) * | 2014-01-15 | 2016-06-29 | The Boeing Company | Multi-level/multi-threshold/multi-persistency GPS/GNSS atomic clock monitoring |
CN105680857A (zh) * | 2016-01-06 | 2016-06-15 | 中国科学院上海光学精密机械研究所 | 原子钟的实时误差补偿方法 |
CN106444343A (zh) * | 2016-11-28 | 2017-02-22 | 中国科学院武汉物理与数学研究所 | 一种微型原子钟数字控制电路装置及方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783412A (zh) * | 2017-11-03 | 2018-03-09 | 中国科学院国家授时中心 | 双光子跃迁铷原子钟 |
CN107783412B (zh) * | 2017-11-03 | 2019-11-19 | 中国科学院国家授时中心 | 双光子跃迁铷原子钟 |
CN110928173A (zh) * | 2018-09-20 | 2020-03-27 | 中国计量科学研究院 | 减小喷泉钟频率偏移不确定度的调节装置 |
CN110928173B (zh) * | 2018-09-20 | 2021-10-22 | 中国计量科学研究院 | 减小喷泉钟频率偏移不确定度的调节装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107219750B (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9348868B2 (en) | Event processing with XML query based on reusable XML query template | |
CN103219946B (zh) | 极坐标发射器、调频路径及方法、参考相位产生器及方法 | |
JP3817550B2 (ja) | 外部early/late入力端子を有するクロック・データ・リカバリ・システム | |
CN104836580A (zh) | 任意相位轨迹频率合成器 | |
CN107219750A (zh) | 改善原子钟稳定度的延时线实时噪声修正方法 | |
CN101420294B (zh) | 一种时钟锁相环控制方法及装置 | |
CN102394593A (zh) | 一种整型化的最小均方lms自适应滤波器和方法 | |
WO2017107901A1 (zh) | 扩频时钟产生装置和生成扩频时钟信号的方法 | |
CN105391447A (zh) | 一种相位插值扩频时钟产生方法 | |
CN110504972A (zh) | 一种基于fpga的采样率转换方法、装置及数模转换方法、装置 | |
CN105830340A (zh) | 用于传送信号的电路和方法 | |
CN113010502B (zh) | 数据质量稽核方法、装置、设备和存储介质 | |
CN105354008A (zh) | 一种随机数生成器的输出电路及输出方法 | |
CN101310435B (zh) | 多路径角度调制器增益的连续校准方法 | |
CN102394640A (zh) | 延时锁定环电路及快速锁定算法 | |
US9429979B2 (en) | Circuit and method for producing USB host working clock | |
CN102325111B (zh) | 一种产生gfsk基带信号的方法 | |
CN102130685A (zh) | 一种原子频标仿真方法和仿真器 | |
CN103888105A (zh) | 一种宽频率连续可调的脉宽波数字产生方法及系统 | |
CN201359721Y (zh) | 一种高速除法器 | |
US11303287B1 (en) | Phase coherent frequency synthesis | |
US20070200642A1 (en) | Clock translator and parallel to serial converter | |
US9269345B2 (en) | Signal processing | |
CN105136167A (zh) | 一种光纤陀螺动态分频方法 | |
CN103312343A (zh) | 发射器及其降低频率偏移的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |