CN115128936B - 一种铷原子钟数字伺服的相检和积分装置及方法 - Google Patents
一种铷原子钟数字伺服的相检和积分装置及方法 Download PDFInfo
- Publication number
- CN115128936B CN115128936B CN202210617840.7A CN202210617840A CN115128936B CN 115128936 B CN115128936 B CN 115128936B CN 202210617840 A CN202210617840 A CN 202210617840A CN 115128936 B CN115128936 B CN 115128936B
- Authority
- CN
- China
- Prior art keywords
- phase detection
- signal
- bits
- accumulator
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 48
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 title claims abstract description 38
- 229910052701 rubidium Inorganic materials 0.000 title claims abstract description 37
- 230000010354 integration Effects 0.000 title claims abstract description 25
- 238000013139 quantization Methods 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims abstract description 21
- 230000001629 suppression Effects 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 15
- 238000009825 accumulation Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 abstract description 14
- 230000008569 process Effects 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 238000012935 Averaging Methods 0.000 abstract description 2
- 239000013078 crystal Substances 0.000 description 21
- 238000005516 engineering process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method 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/26—Automatic control of frequency or phase; Synchronisation using energy levels of molecules, atoms, or subatomic particles as a frequency reference
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
一种铷原子钟数字伺服的相检和积分装置及方法,包括一个“带符号控制的快速累加器”:实现铷原子钟伺服的相检和积分功能;实现对模数和数模转换过程的量化噪声的抑制。“带符号控制”指的是将光检交流误差信号Sn(k)与参考电压P相乘算法简化为符号控制单元实现相检功能;“快速”指的是采用远高于参考电压P的累加频率fclk;“累加器”指的是用长度为N位的累加器实现积分功能,在低n位对Sn'(k)和DN(k)进行累加并溢出至高位,截取累加值DN(k+1)的高m位作为直流纠偏信号Vm(k+1)。本方法仅使用少量处理器资源,在寄存器发生单粒子翻转效应时,可通过锁频环路实时纠正错误数据,同时通过长周期平均的方法降低错误数据的影响,大幅增强抗空间单粒子翻转能力。
Description
技术领域
本发明涉及一种铷原子钟数字伺服的相检和积分装置及方法,属于电子信息技术领域。
背景技术
在铷原子钟内,数字伺服不仅有调试便利和简化硬件资源的优点,还可以克服模拟伺服中存在的开关通断不足、器件温度敏感、积分参数漂移等难题。但数字伺服需要完成信号在模数及数模之间转化:①模数转换使用AD器件,将物理***输出的一路光检电压转为数字量的交流误差信号;②数模转换使用DA器件,将数字量的纠偏信号转为一路电压控制晶振输出频率。在模拟电压与数字信号相互转化过程中引入的量化噪声大小与AD、DA器件的有效位数直接相关:位数越少、量化噪声越大。过大的量化噪声引起信噪比下降,导致频率稳定度性能恶化。因此,对于较高频率稳定度性能要求的铷钟产品,需要采用足够多有效位数的AD、DA器件,不论从成本上还是器件获取渠道上存在较大限制。同时,让AD、DA器件更多的有效位数在电路中发挥作用,也需要在硬件电源、布局和优化算法等相关硬件软件资源上引入繁琐的处理措施。在以往的软件认知中,高频率稳定度需求的铷钟产品不适用采用数字伺服。随着数字技术不断发展,AD的过采样技术和DA的过采样调制技术均能大幅抑制量化噪声,可通过软件的办法来降低对硬件AD和DA位数的需求。特别是在航天领域,数字处理器使用了大量寄存器,存在空间单粒子翻转效应风险,软件中寄存器的减少对于星载铷原子钟特别重要。
综上,现有铷原子钟的数字伺服检相和积分的软件功能配套的硬件AD、DA器件的位数需求较高;检相和积分方法较复杂;在空间环境应用中存在寄存器存在单粒子翻转效应的风险。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种铷原子钟数字伺服的相检和积分装置及方法,解决了模拟伺服中存在的开关通断不足、器件温度敏感、积分参数漂移等问题。
本发明的技术解决方案是:一种铷原子钟数字伺服的相检和积分装置,包括符号控制单元和累加器;
所述符号控制单元接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号并作为累加器输入;
所述累加器接收相检信号以及时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号。
进一步地,所述符号控制的方法包括全相检方法和半相检方法;
所述全相检方法包括:
利用交流误差信号全时段信号,在参考电压的高低电平控制下对交流误差信号赋予正值或者负值,得到相检信号Sn'(k)。
进一步地,所述半相检方法包括:
利用交流误差信号一半时段信号,选择在参考电压的某一电平下,将交流误差信号简化为0,得到为0的相检信号Sn'(k),保持参考电压另一个电平下的相检信号Sn'(k)不变。
进一步地,N>n,N>m。
进一步地,所述对相检信号进行N位累加时,不进行清零,不对溢出数据处理,通过铷原子钟锁频环路形成直流纠偏信号的动态平衡。
进一步地,对N位累加器进行累加,低n位溢出数据进位到高m位,形成对直流纠偏信号低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制。
进一步地,所述截取累加器高位m位作为直流纠偏信号,对检相信号进行高于Nyquist频率10倍及以上采样累加,并进行2N-m倍的平滑处理。
根据所述的一种铷原子钟数字伺服的相检和积分装置实现的相检和积分方法,包括:
接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号;
接收时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号。
进一步地,所述符号控制的方法全相检方法和半相检方法;
所述全相检方法包括:
利用交流误差信号全时段信号,在参考电压的高低电平控制下对交流误差信号赋予正值或者负值,得到相检信号Sn'(k);
所述半相检方法包括:
利用交流误差信号一半时段信号,选择在参考电压的某一电平下,将交流误差信号简化为0,得到为0的相检信号Sn'(k),保持参考电压另一个电平下的相检信号Sn'(k)不变;
N>n,N>m;
所述对相检信号进行N位累加时,不进行清零,不对溢出数据处理,通过铷原子钟锁频环路形成直流纠偏信号的动态平衡;
对N位累加器进行累加,低n位溢出数据进位到高m位,形成对直流纠偏信号低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制;
所述截取累加器高位m位作为直流纠偏信号,对检相信号进行高于Nyquist频率10倍及以上采样累加,并进行2N-m倍的平滑处理。
一种计算机可读存储介质,所述的计算机可读存储介质存储有计算机程序,所述的计算机程序被处理器执行时实现一种铷原子钟数字伺服的相检和积分方法的步骤。
本发明与现有技术相比的优点在于:数字伺服的相检和积分软件功能与AD过采样技术、DA过采样调制技术融合到一个带符号控制的快速累加器中;在锁频环路的闭环情况下,带符号控制的快速累加器抑制了量化噪声,降低了AD、DA器件的位数需求;一个简单的特殊构成累加器,简化了清零、进位、多级累加、滤波等功能;本方法的累加器不存在计数清零、累加器的快速刷新和铷原子钟的锁频环路自我纠偏的特性相结合,能大幅降低单粒子风险。具体为:
(1)本发明通过累加器快速累加,降低了对配套的AD、DA硬件的位数需求;
(2)本发明通过一级累加器,***稳定性好,减少了处理器资源的使用;
(3)本发明通过寄存器不断实时刷新,实现锁频环路实现空间环境下单粒子效应的自我纠偏,大幅增强抗空间单粒子翻转能力。
附图说明
图1数字伺服软件与整机的关系示意图;
图2带符号控制的快速累加器原理框图;
图3 N位累加器运行原理;
图4快速采样时钟改变噪声功率谱分布的原理示意图;
图5、6应用本发明的铷钟整机在常压下的频率稳定度实测图。
具体实施方式
为了更好的理解上述技术方案,下面通过附图以及具体实施例对本申请技术方案做详细的说明,应当理解本申请实施例以及实施例中的具体特征是对本申请技术方案的详细的说明,而不是对本申请技术方案的限定,在不冲突的情况下,本申请实施例以及实施例中的技术特征可以相互组合。
以下结合说明书附图对本申请实施例所提供的一种铷原子钟数字伺服的相检和积分装置及方法做进一步详细的说明,具体实现方式可以包括(如图1~6所示):
数字伺服软件在铷原子钟锁频环路内的关系示意如图1所示。频综输出的6.8…GHz微波信号在物理***中激发铷原子能级跃迁,影响到物理***中的光强变化。光强变化经光电转换为光检电压,电压经隔直交流放大和AD采样,转换为k时刻的n位数字交流误差信号Sn(k)。调制微波的信号经过相位处理后生成参考信号P。Sn(k)与参考电压P经过本软件处理后生成直流纠偏信号Vm(k),Vm(k)经过DA转换为一路控制电压控制晶振的输出频率调整微波信号的频率,实现铷钟的锁频闭环。
将伺服的基本功能相检积分与AD过采样技术、DA过采样调制技术融合到一个带符号控制的快速累加器中。
基于此,本发明提出一种铷原子钟数字伺服的相检和积分装置,包括符号控制单元和累加器;
所述符号控制单元接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号并作为累加器输入;
所述累加器接收相检信号以及时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号。
进一步,在一种可能实现的方式中,所述符号控制的方法包括全相检方法和半相检方法;
所述全相检方法包括:
利用交流误差信号全时段信号,在参考电压的高低电平控制下对交流误差信号赋予正值或者负值,得到相检信号Sn'(k)。
在一种可能实现的方式中,所述半相检方法包括:
利用交流误差信号一半时段信号,选择在参考电压的某一电平下,将交流误差信号简化为0,得到为0的相检信号Sn'(k),保持参考电压另一个电平下的相检信号Sn'(k)不变。
进一步,N>n,N>m。
在一种可能实现的方式中,所述对相检信号进行N位累加时,不进行清零,不对溢出数据处理,通过铷原子钟锁频环路形成直流纠偏信号的动态平衡。
进一步,在一种可能实现的方式中,对N位累加器进行累加,低n位溢出数据进位到高m位,形成对直流纠偏信号低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制。
可选的,在一种可能实现的方式中,所述截取累加器高位m位作为直流纠偏信号,对检相信号进行高于Nyquist频率10倍及以上采样累加,并进行2N-m倍的平滑处理。
在图1中的软件部分可以简化为一个“带符号控制的快速累加器”,原理如图2所示。基于与图2相同的发明构思,本发明还提供一种铷原子钟数字伺服的相检和积分方法,包括:
接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号;
接收时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号。
具体的,该方法包括以下步骤:
S1,“带符号控制”是将交流误差信号Sn(k)与参考电压P相乘算法简化为符号控制单元。符号控制方法可分为全相检和半相检两种方法:全相检利用了Sn(k)全时段信号,即在参考电压P的高低电平控制下对Sn(k)赋予正值或者负值,得到相检信号Sn'(k)。半相检与全相检的不同点是利用了Sn(k)一半时段信号,即选择在参考电压P的某一电平下,强制将Sn(k)简化为0,得到相检信号Sn'(k)也为0,保持参考电压P另一个电平下的Sn'(k)不变。具体的控制取值方法须确保锁频环路硬件与软件符合规律:当微波信号中心频率f低于铷原子跃迁谱线中心频率f0时,输出的直流纠偏信号控制基于晶振的微波信号中心频率f变大;当微波信号中心频率f高于铷原子跃迁谱线中心频率f0时,输出的直流纠偏信号控制基于晶振的微波信号中心频率f变小。“快速”是采用远高于参考电压P的累加频率fclk,可由晶振频率分频获得。远高于指高于10倍及以上。
S2,“累加器”为N位累加器实现积分功能,运行过程如图3所示。在低n位对Sn'(k)和DN(k)进行累加,累加值溢出到高位;截取累加值DN(k+1)的高m位作为直流纠偏信号Vm(k+1)。在锁频环路正常工作状态下,m位直流纠偏信号Vm(k+1)处于纠偏值附近小幅波动,解决了数据累加溢出的问题。
由图2可知,***函数H(Z)=DN(Z)/Sn'(Z)=1-Z-1,存在一个极点0,并在单位圆内,是一个稳定***。
进一步地,所述模数过程量化噪声的抑制是将采集的n位交流误差信号Sn(k)经相检后生成的Sn'(k)信号快速长周期累加,取高m位数据等同于取平均值,形成对信号的过采样低通滤波,实现带内的量化噪声抑制。这也要求硬件AD器件的采样速率高于或等于累加频率fclk。
进一步地,所述数模过程量化噪声的抑制是将通过N位的累加器的快速累加,低n位溢出数据快速进位到高m位,形成对直流纠偏信号Vm(k)低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制。这也要求硬件DA器件的转换速率高于或等于累加频率fclk。
以下结合说明书附图对本申请实施例做进一步详细的说明,具体实现方式可以包括:由图2可知,本发明构成源自于参数的繁琐计算,需要编程人员对铷钟锁频环路具备初步的知识储备,根据锁频环路硬件设计状态确定符号控制单元取值的方法与5个基本参数:***时钟fclk、Sn(k)位数n、直流纠偏信号Vm(k)位数m、累加器位数N、DN(k)的初始值DN(0)。其约束条件及获取方法如下所示,可根据实际锁频环路设计方法调整相应顺序或采用其他方法:
1)、相检的基本算法是Sn'(k)=P×Sn(k),可分为半相检和全相检两种。具体赋值方法如表1所示:
表1不同相检方法对应的Sn'(k)取值方法
i.全相检。当P电平为‘0’状态时,设定P=-1,此时乘法器Sn'(k)=P×Sn(k)结果为-Sn(k);当P电平为‘1’状态时,设定P=1,此时乘法器Sn'(k)=P×Sn(k)结果为Sn(k)。因此乘法器简化为控制器:当P为‘0’时,Sn(k)取负值;当P为‘1’时,Sn(k)取正值。
ii.半相检。当P电平为‘0’状态时,设定P=0,此时乘法器Sn'(k)=P×Sn(k)结果应为0;当P电平为‘1’状态时,设定P=1,此时乘法器Sn'(k)=P×Sn(k)结果为Sn(k)。因此乘法器简化为控制器:当P为‘0’时,Sn(k)取负值;当P为‘1’时,Sn(k)取正值。
iii.值得注意的是,Sn(k)本应是交流信号,但有些AD器件采样后会带有直流分量。在全相检方法中,通过后级累加器能去除直流分量;但在半相检方法中,需要提前去除直流分量,若不去除直流分量,需在符号控制单元中将取值’0’改变为一个固定值用于抵消直流分量的影响。
2)、获得n位的交流误差信号Sn(k)。在***时钟fclk触发下,k时刻光检电压通过AD转为n位的交流误差信号,n的取值受限于AD最大位数。例如采用12位AD,则n≤12。一般情况下,n越大,Sn(k)的量化噪声越小。
3)、获得m位的直流纠偏信号Vm(k)。在k时刻截取累加器的高m位作为直流纠偏信号Vm(k),m取值需要根据铷钟的性能要求,满足以下两个条件:
i.过小的m无法保证晶振的锁定,即微波频率中心f距离物理***中心频率f0的差值应远小于谱线带宽δf:f-f0<<δf。这要求晶振在压控电压控制下的频率步进足够小,即m足够大;
例如:当δf为500Hz,若取f-f0≤5Hz,晶振在压控电压控制下的频率步进当f晶为100MHz,Vc=10V,K晶=10Hz,可算得m≥5。
ii.m取值决定了数字伺服相检积分的量化误差,m越高相检积分的量化误差越小,经过采样调制后可大幅抑制量化误差对锁频环路频移的影响,其极限抑制能力可按公式1进行计算获得:
式中:
为相对频移;
f晶为晶振频率;
KV为晶振压控斜率;
Vc为晶振压控电压范围;
m为截取累加高位结果送给DA的Vm(k)位数;
A0为锁频环路增益。
例如:Vc=10V,m=10,f晶=100MHz,Kv=10Hz/V,A0=5×105;经公式1,可算出由m带来的量化误差经过采样调制后,造成的相对频移抑制能力极限值为1.95E-15。
4)、设置累加器长度N。累加器长度N由三部分组成:高m位、中m'位、低n位,即N=m+m'+n。高m位为输出压控信号Vm(k)的位数,低n位为交流误差信号Sn(k)的位数。m'值由锁频环路的时间常数T决定,其关系如公式2所示:
式中:
T为锁频环路时间常数;
fclk为***时钟频率;
m为输出压控信号Vm(k)的位数。
例如:锁频环路参数n=10、m=10、fclk=200KHz、T=10s;代入公式2,可算出m'=10.95,取整后m'=11。此时累加器长度N=31位。
5)、获得累加器初始值DN(0)。在启动时刻,DN(0)的高m位Vm(0)经DA转换后的电压控制晶振频率在锁频环路的捕捉带内即可。由此晶振压控输出频率使得锁频环路的正确建立。
6)、设置***时钟fclk。在数模和模数转化技术中,AD过采样和DA过采样调制方法均是建立在快速时钟的基础上。如图4所示,利用功率密度谱的矩形面积来表示奈奎斯特采样与过采样的量化噪声的能量。A1(f)表示奈奎斯特采样的功率密度谱,A2(f)表示过采样的功率密度谱。两种采样方法的噪声功率密度谱面积一样,仅分布上有所区别,即总的量化噪声功率不变。噪声功率谱从带宽为2fP的A1(f)变为带宽为fclk/2的A2(f)。当fclk>>fP,噪声幅值大幅降低,同时锁频环路的带宽非常窄,一般在秒量级上,可滤除锁频环路带宽外的噪声。在满足AD、DA及处理器等硬件响应速率条件下,时钟越高,对带内量化噪声抑制能力越强。通过快速时钟的触发,将光检电压转变为交流误差信号Sn(k),在Sn(k)经相检转换为Sn'(k)后,再用快速长位数累加器滤除Sn'(k)的带外噪声,实现带内信噪比的提升。但在后续计算累加器位数N时,在固定的锁频环路时间常数T下,fclk越高累加器位数N越长,对处理器的资源占用越多。因此,需要根据实际情况选择适当的参数。就***时钟fclk给出了公式3与公式4的近似估计方法,但不限于此方法:
i.根据AD过采样方法增强信号的信噪比,优于物理***输出交流误差信号的信噪比。可由公式3近似获得fclk频率最小值判决条件f'clk:
式中:
fp为参考信号频率;
σs锁频环路的频率稳定度;
τ为频率采样时间;
Q为原子跃迁谱线频差与线宽的比值;
n为数字交流误差信号Sn(k)位数。
例如:锁频环路参数Q=14000000、n=10、fp=80Hz;经公式3,可算出f′clk≥9KHz。
ii.根据DA一阶过采样调制的噪声幅值与晶振参数、频移的关系。可由公式4近似获得fclk频率最小值判决条件f"clk。
式中:
为相对频移误差;
Vc为晶振压控电压范围;
KV为晶振压控斜率;
τ为频率采样时间;
f晶为晶振频率;
m为截取累加高位结果送给DA的Vm(k)位数。
例如:τ=1s,Vc=10V,m=10,f晶=100MHz,Kv=10Hz/V,经公式4,可算出f"clk≥6.2KHz。
iii.取f'clk、和f"clk的大值作为fclk下限,即fclk≥f'clk且fclk≥f"clk。实际应用中,还需考虑n、N、m三个数值之间关系对fclk的附带要求:若Sn'(k)为大幅值起伏信号,则N需要足够大才能避免Sn'(k)幅值变化直接进入Vm(k),由此需要更高频率的时钟,建议实际应用中fclk≥50KHz。另,存在一组特殊关系:fclk=2M×fp,M为正整数。此特殊关系下,累加器在任何k时刻输出的数据为参考电压P的整周期倍数的相检信号的累加值。
用VHDL语句对图2逻辑进行描述,全相检的例子如下:
注:
En 启动信号(建议在锁频环路硬件完成预热后启动)
fclk ***时钟
P 参考信号(要求已完成相位调整)
Sn n位的交流误差信号
Vm m位的直流纠偏信号
DN N位的累加器寄存器,其初始值需要保证晶振输出频率在锁频环路的捕捉带内。
通过本发明,可大幅降低与配套的AD与DA硬件的位数需求。实际测试结果见图5,在本发明的应用下,12位的AD、DA器件能满足铷钟整机优于的频率稳定度性能需求,满足铷钟整机优于1E-12的频率准确度调整精度需求。在极限情况下取AD位数为1位(即不使用AD器件),通过处理器直接将交流误差信号转换为1位Sn(k)在保证铷钟正常锁定的同时,对频率稳定度的影响小于1E-11/1s。
综上所述,本发明的核心是一个带符号控制的快速累加器,使用极少的处理器的资源特别是少量的寄存器,对时钟速率要求不高,可用于各类处理器实现软件功能。软件中寄存器的不断自我刷新与锁频环路对输入端寄存器的纠偏作用相结合,能大幅增强算法的抗单粒子能力,特别适用于星载铷原子钟。
本申请技术方案中,过采样和过采样调制技术降低量化噪声,降低软件对配套硬件AD、DA的位数需求。将多阶***简化为一阶***,增强软件的稳定性。极简的软件构成,降低软件占用资源特别是寄存器的数量。通过锁频环路实时纠正软件中寄存器错误数据;快速累加并取高位数据的方法相当于在长周期内取平均值,可大幅降低寄存器错误数据瞬间对频率的影响;增强了数字伺服软件的抗单粒子翻转能力。
综上,本发明可以解决铷原子钟数字伺服电路对硬件AD、DA的位数需求较高的问题;现有数字伺服检相和积分方案较为复杂性占用资源多的问题;空间环境应用中检相和积分中寄存器存在单粒子翻转的问题。
本申请提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机指令,当所述计算机指令在计算机上运行时,使得计算机执行图1所述的方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (5)
1.一种铷原子钟数字伺服的相检和积分装置,其特征在于:包括符号控制单元和累加器;
所述符号控制单元接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号并作为累加器输入;
所述累加器接收相检信号以及时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号;
所述符号控制的方法包括全相检方法和半相检方法;
所述全相检方法包括:
利用交流误差信号全时段信号,在参考电压的高低电平控制下对交流误差信号赋予正值或者负值,得到相检信号Sn'(k);
所述半相检方法包括:
利用交流误差信号一半时段信号,选择在参考电压的某一电平下,将交流误差信号简化为0,得到为0的相检信号Sn'(k),保持参考电压另一个电平下的相检信号Sn'(k)不变;
N>n,N>m;
所述对相检信号进行N位累加时,不进行清零,不对溢出数据处理,通过铷原子钟锁频环路形成直流纠偏信号的动态平衡;
对N位累加器进行累加,低n位溢出数据进位到高m位,形成对直流纠偏信号低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制。
2.根据权利要求1所述的一种铷原子钟数字伺服的相检和积分装置,其特征在于,所述截取累加器高位m位作为直流纠偏信号,对检相信号进行高于Nyquist频率10倍及以上采样累加,并进行2N-m倍的平滑处理。
3.根据权利要求1所述的一种铷原子钟数字伺服的相检和积分装置实现的相检和积分方法,其特征在于,包括:
接收n位数字交流误差信号和参考电压信号,进行交流误差信号与参考电压相乘的符号控制,生成相检信号;
接收时钟信号,以时钟信号的频率对相检信号进行N位累加后,截取累加器高位m位作为直流纠偏信号。
4.根据权利要求3所述的方法,其特征在于,所述符号控制的方法包括全相检方法和半相检方法;
所述全相检方法包括:
利用交流误差信号全时段信号,在参考电压的高低电平控制下对交流误差信号赋予正值或者负值,得到相检信号Sn'(k);
所述半相检方法包括:
利用交流误差信号一半时段信号,选择在参考电压的某一电平下,将交流误差信号简化为0,得到为0的相检信号Sn'(k),保持参考电压另一个电平下的相检信号Sn'(k)不变;
N>n,N>m;
所述对相检信号进行N位累加时,不进行清零,不对溢出数据处理,通过铷原子钟锁频环路形成直流纠偏信号的动态平衡;
对N位累加器进行累加,低n位溢出数据进位到高m位,形成对直流纠偏信号低位的过采样调制,通过锁频环路的窄带特性实现带内的量化噪声抑制;
所述截取累加器高位m位作为直流纠偏信号,对检相信号进行高于Nyquist频率10倍及以上采样累加,并进行2N-m倍的平滑处理。
5.一种计算机可读存储介质,所述的计算机可读存储介质存储有计算机程序,其特征在于,所述的计算机程序被处理器执行时实现如权利要求3或4所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210617840.7A CN115128936B (zh) | 2022-06-01 | 2022-06-01 | 一种铷原子钟数字伺服的相检和积分装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210617840.7A CN115128936B (zh) | 2022-06-01 | 2022-06-01 | 一种铷原子钟数字伺服的相检和积分装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115128936A CN115128936A (zh) | 2022-09-30 |
CN115128936B true CN115128936B (zh) | 2024-05-31 |
Family
ID=83377834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210617840.7A Active CN115128936B (zh) | 2022-06-01 | 2022-06-01 | 一种铷原子钟数字伺服的相检和积分装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115128936B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491451A (en) * | 1994-11-17 | 1996-02-13 | Ball Corporation | Method and apparatus for reduction of atomic frequency standard phase noise |
US5663733A (en) * | 1995-08-28 | 1997-09-02 | Trimble Navigation Limited | Digital bandwidth compression for optimum tracking in satellite positioning system receiver |
CN101018056A (zh) * | 2007-02-16 | 2007-08-15 | 中国科学院武汉物理与数学研究所 | 铷原子频标数字伺服装置 |
CN201008146Y (zh) * | 2007-02-16 | 2008-01-16 | 中国科学院武汉物理与数学研究所 | 铷原子频标数字伺服装置 |
CN101951261A (zh) * | 2010-10-08 | 2011-01-19 | 中国科学院武汉物理与数学研究所 | 一种被动型相干布居数囚禁原子频标的全数字伺服装置 |
CN202231700U (zh) * | 2011-08-24 | 2012-05-23 | 江汉大学 | 用于铷原子频标的伺服电路及铷原子频标 |
CN103326717A (zh) * | 2013-05-10 | 2013-09-25 | 西安空间无线电技术研究所 | 一种铷钟扫描捕获辅助锁定方法 |
CN111373886B (zh) * | 2013-05-10 | 2015-06-10 | 西安空间无线电技术研究所 | 一种星载铷钟fsk电路 |
CN105811973A (zh) * | 2016-03-17 | 2016-07-27 | 江汉大学 | 一种被动型铷原子频标 |
CN107483050A (zh) * | 2017-07-12 | 2017-12-15 | 西安空间无线电技术研究所 | 一种基于实时跟踪技术的原子钟平稳切换*** |
CN113949486A (zh) * | 2021-09-23 | 2022-01-18 | 武汉正维电子技术有限公司 | 基于符号累加和相关运算的ads_b信号解析方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6946915B2 (en) * | 2003-03-17 | 2005-09-20 | Xiaopin Zhang | Maximally digitized fractional-N frequency synthesizer and modulator with maximal fractional spurs removing |
-
2022
- 2022-06-01 CN CN202210617840.7A patent/CN115128936B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491451A (en) * | 1994-11-17 | 1996-02-13 | Ball Corporation | Method and apparatus for reduction of atomic frequency standard phase noise |
US5663733A (en) * | 1995-08-28 | 1997-09-02 | Trimble Navigation Limited | Digital bandwidth compression for optimum tracking in satellite positioning system receiver |
CN101018056A (zh) * | 2007-02-16 | 2007-08-15 | 中国科学院武汉物理与数学研究所 | 铷原子频标数字伺服装置 |
CN201008146Y (zh) * | 2007-02-16 | 2008-01-16 | 中国科学院武汉物理与数学研究所 | 铷原子频标数字伺服装置 |
CN101951261A (zh) * | 2010-10-08 | 2011-01-19 | 中国科学院武汉物理与数学研究所 | 一种被动型相干布居数囚禁原子频标的全数字伺服装置 |
CN202231700U (zh) * | 2011-08-24 | 2012-05-23 | 江汉大学 | 用于铷原子频标的伺服电路及铷原子频标 |
CN103326717A (zh) * | 2013-05-10 | 2013-09-25 | 西安空间无线电技术研究所 | 一种铷钟扫描捕获辅助锁定方法 |
CN111373886B (zh) * | 2013-05-10 | 2015-06-10 | 西安空间无线电技术研究所 | 一种星载铷钟fsk电路 |
CN105811973A (zh) * | 2016-03-17 | 2016-07-27 | 江汉大学 | 一种被动型铷原子频标 |
CN107483050A (zh) * | 2017-07-12 | 2017-12-15 | 西安空间无线电技术研究所 | 一种基于实时跟踪技术的原子钟平稳切换*** |
CN113949486A (zh) * | 2021-09-23 | 2022-01-18 | 武汉正维电子技术有限公司 | 基于符号累加和相关运算的ads_b信号解析方法及*** |
Non-Patent Citations (5)
Title |
---|
Digital servo system based on FPGA for optically pumped magnetometer;Sheng Zhou 等;《2015 Joint Conference of the IEEE International Frequency Control Symposium & the European Frequency and Time Forum》;1-3 * |
基于FPGA的星载铷钟数字伺服电路设计;何冬;刘昶;屈勇晟;贺玉玲;;空间电子技术(05);4-7 * |
基于单片机的锁相程控频率合成技术应用于铷频标的研究;王玉泰, 邢宝玲, 王前虹, 赵建顺;济南大学学报(自然科学版);20011025(03);19-21 * |
数字技术在星载铷原子钟的应用;刘昶 等;《第十二届中国卫星导航年会论文集——S05 空间基准与精密定位》;1-9 * |
铷原子频标数字伺服***;曹远洪;张贤谊;康松柏;钟达;梅刚华;;波谱学杂志(04);125-131 * |
Also Published As
Publication number | Publication date |
---|---|
CN115128936A (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI384760B (zh) | 全數位鎖相迴路 | |
CN108667458B (zh) | 能够消除来自σ-δ调制器的量化噪声的分数n数字pll | |
US10911054B2 (en) | Digital-to-time converter (DTC) assisted all digital phase locked loop (ADPLL) circuit | |
JP5347534B2 (ja) | 位相比較器、pll回路、及び位相比較器の制御方法 | |
TW201412023A (zh) | 具有增益校正之鎖相迴路、用於鎖相迴路之增益量測方法、校正方法及抖動量測方法 | |
CN109995360B (zh) | 抑制扰动的锁相环 | |
US7482880B2 (en) | Frequency modulated output clock from a digital frequency/phase locked loop | |
US10979059B1 (en) | Successive approximation register analog to digital converter based phase-locked loop with programmable range | |
CN102210103A (zh) | 利用σ-δ调制器控制的分频器进行的频率变换 | |
JP2004507916A (ja) | ノイズシェーピング・デジタル周波数合成 | |
CN112134558A (zh) | 具有锁频环的全数字锁相环(adpll) | |
CN110324038B (zh) | 使用锁相环来快速建立斜坡生成 | |
JP2004260791A (ja) | フラクショナルn周波数シンセサイザ及びシンセサイズ方法 | |
CN115128936B (zh) | 一种铷原子钟数字伺服的相检和积分装置及方法 | |
JPH0795055A (ja) | ディジタル位相同期装置 | |
JPH04212522A (ja) | 周波数合成装置 | |
Brandonisio et al. | Noise-shaping all-digital phase-locked loops | |
US11088698B2 (en) | Phase-locked loop circuit | |
CN109391267B (zh) | 使用带adcs和dac的数字plls时的抖动减少技术 | |
CN113852370A (zh) | 一种相位抖动补偿方法、模块和数字锁相环 | |
CN101667832B (zh) | 分数-n锁相装置与方法 | |
US11632118B2 (en) | Closed-loop oscillator based sensor interface circuit | |
US20080129565A1 (en) | Apparatus for and method of sigma-delta modulation | |
US8494105B1 (en) | Apparatus and method for providing digital representation of time difference between clocks | |
US11876523B1 (en) | DPLL timing normalization |
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 |