CN112702060A - 一种用于铯原子钟的微波幅度控制方法 - Google Patents

一种用于铯原子钟的微波幅度控制方法 Download PDF

Info

Publication number
CN112702060A
CN112702060A CN202011453080.8A CN202011453080A CN112702060A CN 112702060 A CN112702060 A CN 112702060A CN 202011453080 A CN202011453080 A CN 202011453080A CN 112702060 A CN112702060 A CN 112702060A
Authority
CN
China
Prior art keywords
amplitude
value
microwave
control method
transition probability
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.)
Pending
Application number
CN202011453080.8A
Other languages
English (en)
Inventor
马沛
董鹏玲
赵玉龙
汪东军
王骥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lanzhou Institute of Physics of Chinese Academy of Space Technology
Original Assignee
Lanzhou Institute of Physics of Chinese Academy of Space Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Lanzhou Institute of Physics of Chinese Academy of Space Technology filed Critical Lanzhou Institute of Physics of Chinese Academy of Space Technology
Priority to CN202011453080.8A priority Critical patent/CN112702060A/zh
Publication of CN112702060A publication Critical patent/CN112702060A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03LAUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
    • H03L7/00Automatic control of frequency or phase; Synchronisation
    • H03L7/26Automatic 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)
  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)

Abstract

本发明公开了一种用于铯原子钟的微波幅度控制方法,通过根据设定的微波信号幅度计算不同微波幅度下对应的最大原子跃迁几率,跃迁几率的最大值对应的幅度值即为控制的目标值,采用本发明提供的方法能够有效地提高铯原子钟的频率准确度和长期频率稳定度。

Description

一种用于铯原子钟的微波幅度控制方法
技术领域
本发明属于原子钟技术领域,具体涉及一种用于铯原子钟的微波幅度控制方法。
背景技术
铯原子钟是一级频标,有准确度高、长期稳定度好的特点,被广泛使用于授时、守时等技术领域。目前国产磁选态铯原子钟是针对这类有严格控温的实验室环境所研制。一旦失去精准控温的实验室环境,国产磁选态铯原子钟的长期稳定度将会出现恶化。现有的国产磁选态铯原子钟产品没有通过对微波幅度的控制补偿铯束管由于温度变化所引起的腔体失谐,从而降低磁选态铯原子钟对温度的敏感,提高其在非实验室环境下的长期频率稳定度。
发明内容
有鉴于此,本发明提供了一种用于铯原子钟的微波幅度控制方法,能够用于提高铯原子钟的频率准确度和长期频率稳定度,也可为其他类型原子钟的微波幅度控制提供参考。
本发明提供了一种用于铯原子钟的微波幅度控制方法,包括以下步骤:
设置所述铯原子钟的微波信号的幅度值,计算所述幅度值对应的原子跃迁几率,所述原子跃迁几率的最大值所对应的微波幅度值即为微波幅度控制的目标值。
进一步地,所述微波信号的幅度值的初值为微波信号源输出幅度最大值的1/50。
进一步地,所述计算所述幅度值对应的原子跃迁几率的过程包括以下步骤:
计算所述幅度值An对应的原子跃迁几率的最大值Imax n=F(An,vx),其中,An=A0+nAd,A0为所述幅度值的初始值,Ad为所述幅度值的固定步进,n为固定步进的系数,vx为微波信号源的频率值。
进一步地,所述固定步进Ad为微波信号源输出幅度最大值的1/20。
进一步地,所述微波幅度控制方法还包括以下步骤:
保持所述幅度值的初值不变,以微波信号源的标称频率v0为中心进行扫频;计算在所述幅度值下原子跃迁几率的最大值Imax n=F(An,vx);
记录最近的三个扫频时的最大值Imax n=F(An,vx1)、Imax(n+1)=F(An+1,vx2)、Imax(n+2)=F(An+2,vx3),直到满足下式时:
Figure BDA0002832186870000021
停止扫描,此时An+1即为所述目标值。
有益效果:
本发明通过根据设定的微波信号幅度计算不同微波幅度下对应的最大原子跃迁几率,跃迁几率的最大值对应的幅度值即为控制的目标值,采用本发明提供的方法能够有效地提高铯原子钟的频率准确度和长期频率稳定度。
附图说明
图1为本发明提供的一种用于铯原子钟的微波幅度控制方法的流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种用于铯原子钟的微波幅度控制方法,如图1所示,具体包括以下步骤:
步骤1、设置输入微波信号的幅度值。
输入的微波信号频率为9192631770Hz,其幅度初值由软件设定,例如幅度初值为A0,一般为微波信号源最大输出幅度的1/50。
设置微波信号幅度的大小,每次以固定步进增加幅度值Ad.
步骤2、计算不同微波幅度值下对应的最大原子跃迁几率。
铯原子的最大原子跃迁几率与铯束管的输出信号幅度成正比。因此,在计算不同微波幅度值下对应的最大原子跃迁几率时,需要通过电路采集铯束管输出信号幅度作为算法的判断依据。
在给定的微波幅度下,以铯原子跃迁频率9192631770Hz为中心频率,±50Hz内进行扫频,根据微波幅度与铯束管输出信号幅度的函数关系式进行计算,得到当前微波幅度下的铯束管输出信号幅度最大值。
再以固定步进,一般为微波信号源最大输出幅度的1/20。增加微波幅度值,再次计算当前微波幅度值下的铯束管输出信号幅度最大值。
以微波幅度的初值固定不变,微波信号源以其标称的频率v0为中心,在v0±vs范围内进行扫频,vs为扫频的范围。计算出在此微波幅度下原子跃迁几率的最大值,记为Imax 0=F(A0,vx),vx为与幅度最大值对应的频率值。
以固定步进Ad,每增加一次微波幅度值,记为An=A0+nAd
再以相同的扫频范围v0±vs扫频,计算在此微波功率下原子跃迁几率的最大值,记为Imax n=F(An,vx)。
计算原子跃迁几率最大时,对应的微波幅度值,该值即为软件控制的目标值。包含以下步骤:
记录最近的3个微波幅度点扫频时的最大值,分别为Imax n=F(An,vx1)、Imax(n+1)=F(An+1,vx2)、Imax(n+2)=F(An+2,vx3),直到第三个幅度点的原子跃迁几率最大值小于第二个幅度点的原子跃迁几率最大值,即:
Figure BDA0002832186870000041
停止扫描。
第二个幅度点An+1即为软件控制的微波幅度目标值。
由于环境温度的变化会造成铯束管微波腔体失谐,引起铯原子跃迁的中心频率发生迁移。通过计算铯原子跃迁几率最大时的微波幅度值,作为软件控制的目标值,能够补偿腔体失谐对频率的恶化,从而降低铯原子钟对环境温度的敏感性,提高铯原子钟的长期稳定度。本发明的使用能使铯原子钟的温度系数提升半个数量级。本发明通过控制微波幅度,补偿铯束管的腔体频率迁移,提高铯束管对温度的敏感程度,提高铯原子钟整机的温度系数,从而改善铯原子钟的环境适应性,提高铯原子钟整机在非温控环境下的长期稳定度。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种用于铯原子钟的微波幅度控制方法,其特征在于,包括以下步骤:
设置所述铯原子钟的微波信号的幅度值,计算所述幅度值对应的原子跃迁几率,所述原子跃迁几率的最大值所对应的微波幅度值即为微波幅度控制的目标值。
2.根据权利要求1所述的微波幅度控制方法,其特征在于,所述微波信号的幅度值的初值为微波信号源输出幅度最大值的1/50。
3.根据权利要求2所述的微波幅度控制方法,其特征在于,所述计算所述幅度值对应的原子跃迁几率的过程包括以下步骤:
计算所述幅度值An对应的原子跃迁几率的最大值Imaxn=F(An,vx),其中,An=A0+nAd,A0为所述幅度值的初始值,Ad为所述幅度值的固定步进,n为固定步进的系数,vx为微波信号源的频率值。
4.根据权利要求3所述的微波幅度控制方法,其特征在于,所述固定步进Ad为微波信号源输出幅度最大值的1/20。
5.根据权利要求1所述的微波幅度控制方法,其特征在于,所述微波幅度控制方法还包括以下步骤:
保持所述幅度值的初值不变,以微波信号源的标称频率v0为中心进行扫频;计算在所述幅度值下原子跃迁几率的最大值Imaxn=F(An,vx);
记录最近的三个扫频时的最大值Imaxn=F(An,vx1)、Imax(n+1)=F(An+1,vx2)、Imax(n+2)=F(An+2,vx3),直到满足下式时:
Figure FDA0002832186860000011
停止扫描,此时An+1即为所述目标值。
CN202011453080.8A 2020-12-11 2020-12-11 一种用于铯原子钟的微波幅度控制方法 Pending CN112702060A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011453080.8A CN112702060A (zh) 2020-12-11 2020-12-11 一种用于铯原子钟的微波幅度控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011453080.8A CN112702060A (zh) 2020-12-11 2020-12-11 一种用于铯原子钟的微波幅度控制方法

Publications (1)

Publication Number Publication Date
CN112702060A true CN112702060A (zh) 2021-04-23

Family

ID=75508808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011453080.8A Pending CN112702060A (zh) 2020-12-11 2020-12-11 一种用于铯原子钟的微波幅度控制方法

Country Status (1)

Country Link
CN (1) CN112702060A (zh)

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB875242A (en) * 1958-02-21 1961-08-16 Varian Associates Atomic stabilized frequency source
US5146184A (en) * 1991-08-01 1992-09-08 Hewlett-Packard Company Atomic clock system with improved servo system
US5852386A (en) * 1997-06-02 1998-12-22 Northrop Grumman Corporation Apparatus and method for microwave field strength stabilization in cell type atomic clocks
US20040233003A1 (en) * 2003-03-11 2004-11-25 William Happer Method and system for operating an atomic clock with reduced spin-exchange broadening of atomic clock resonances
CN101917178A (zh) * 2010-07-29 2010-12-15 北京交通大学 速度传感器输出信号模拟装置及方法
CN102704919A (zh) * 2012-05-24 2012-10-03 中国石油大学(北京) 一种用于旋转磁场井间测距的频变信号消噪方法及装置
CN102799102A (zh) * 2012-07-03 2012-11-28 华中科技大学 一种被动型cpt原子钟的温控参数自整定方法及装置
CN102944854A (zh) * 2012-11-16 2013-02-27 江汉大学 用于测量剩场量绝对值的测量系统及其测量方法
CN203151469U (zh) * 2013-01-31 2013-08-21 江汉大学 能够缩小原子钟的原子能级跃迁动态探测频率范围的系统
CN103338036A (zh) * 2013-07-12 2013-10-02 西安电子科技大学 基于相位群处理的原子钟频率信号链接控制方法
CN103823356A (zh) * 2014-03-07 2014-05-28 中国科学院武汉物理与数学研究所 基于pxi系统的被动型cpt原子钟实验装置及方法
CN104363017A (zh) * 2014-10-09 2015-02-18 兰州空间技术物理研究所 提高原子频标查询信号信噪比的方波矢量锁定放大方法
CN105242521A (zh) * 2015-11-13 2016-01-13 中国科学院武汉物理与数学研究所 一种实现微型cpt原子钟物理系统的装置及方法
JP2017092797A (ja) * 2015-11-13 2017-05-25 セイコーエプソン株式会社 原子発振器
CN106773610A (zh) * 2016-11-28 2017-05-31 北京工业大学 一种铯原子钟和氢钟频差预估方法
CN107561916A (zh) * 2017-10-10 2018-01-09 兰州空间技术物理研究所 一种基于法拉第激光抽运的铷原子微波钟
CN107783412A (zh) * 2017-11-03 2018-03-09 中国科学院国家授时中心 双光子跃迁铷原子钟
CN108614224A (zh) * 2018-04-03 2018-10-02 北京航天控制仪器研究所 一种用于cpt磁力仪的气室工作温度自动标定系统及方法
CN110018631A (zh) * 2019-04-10 2019-07-16 温州激光与光电子协同创新中心 应用法拉第激光抽运探测的铯原子微波频率标准及其实现方法
CN111404431A (zh) * 2020-04-10 2020-07-10 驰驱电气(嘉兴)有限公司 一种伺服电机控制参数自整定控制方法

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB875242A (en) * 1958-02-21 1961-08-16 Varian Associates Atomic stabilized frequency source
US5146184A (en) * 1991-08-01 1992-09-08 Hewlett-Packard Company Atomic clock system with improved servo system
US5852386A (en) * 1997-06-02 1998-12-22 Northrop Grumman Corporation Apparatus and method for microwave field strength stabilization in cell type atomic clocks
US20040233003A1 (en) * 2003-03-11 2004-11-25 William Happer Method and system for operating an atomic clock with reduced spin-exchange broadening of atomic clock resonances
CN101917178A (zh) * 2010-07-29 2010-12-15 北京交通大学 速度传感器输出信号模拟装置及方法
CN102704919A (zh) * 2012-05-24 2012-10-03 中国石油大学(北京) 一种用于旋转磁场井间测距的频变信号消噪方法及装置
CN102799102A (zh) * 2012-07-03 2012-11-28 华中科技大学 一种被动型cpt原子钟的温控参数自整定方法及装置
CN102944854A (zh) * 2012-11-16 2013-02-27 江汉大学 用于测量剩场量绝对值的测量系统及其测量方法
CN203151469U (zh) * 2013-01-31 2013-08-21 江汉大学 能够缩小原子钟的原子能级跃迁动态探测频率范围的系统
CN103338036A (zh) * 2013-07-12 2013-10-02 西安电子科技大学 基于相位群处理的原子钟频率信号链接控制方法
CN103823356A (zh) * 2014-03-07 2014-05-28 中国科学院武汉物理与数学研究所 基于pxi系统的被动型cpt原子钟实验装置及方法
CN104363017A (zh) * 2014-10-09 2015-02-18 兰州空间技术物理研究所 提高原子频标查询信号信噪比的方波矢量锁定放大方法
CN105242521A (zh) * 2015-11-13 2016-01-13 中国科学院武汉物理与数学研究所 一种实现微型cpt原子钟物理系统的装置及方法
JP2017092797A (ja) * 2015-11-13 2017-05-25 セイコーエプソン株式会社 原子発振器
CN106773610A (zh) * 2016-11-28 2017-05-31 北京工业大学 一种铯原子钟和氢钟频差预估方法
CN107561916A (zh) * 2017-10-10 2018-01-09 兰州空间技术物理研究所 一种基于法拉第激光抽运的铷原子微波钟
CN107783412A (zh) * 2017-11-03 2018-03-09 中国科学院国家授时中心 双光子跃迁铷原子钟
CN108614224A (zh) * 2018-04-03 2018-10-02 北京航天控制仪器研究所 一种用于cpt磁力仪的气室工作温度自动标定系统及方法
CN110018631A (zh) * 2019-04-10 2019-07-16 温州激光与光电子协同创新中心 应用法拉第激光抽运探测的铯原子微波频率标准及其实现方法
CN111404431A (zh) * 2020-04-10 2020-07-10 驰驱电气(嘉兴)有限公司 一种伺服电机控制参数自整定控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JIANXIANG WANG: "Design of a Cesium Atomic Clock 1PPS Signal Generation and Synchronization Module Based on FPGA", 《2020 INTERNATIONAL WORKSHOP ON ELECTRONIC COMMUNICATION AND ARTIFICIAL INTELLIGENCE (IWECAI)》 *
廉吉庆: "CPT铯原子钟微波信号源设计与分析", 《宇航计测技术》 *
张程洋: "数字化稀土超磁致伸缩换能器驱动系统的研制", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
马沛: "铯原子钟优化算法中微波幅度控制方法的研究", 《真空与低温》 *

Similar Documents

Publication Publication Date Title
US8237482B2 (en) Circuit and method for generating a clock signal
US20130066578A1 (en) Method and device for determining zero-rate offset of a gyroscope, and system comprising the device
CN105387859A (zh) Mems传感器组合温度漂移误差补偿方法
EP1935089A2 (en) Systems, methods and devices for dual closed loop modulation controller for nonlinear rf amplifier
CN104898550A (zh) 动态伺服系统基于滑模扩张状态观测器的复合控制方法
EP3197050A1 (en) Method for controlling operation of moving average filter
US9582018B1 (en) Automatic gain compression detection and gain control for analog front-end with nonlinear distortion
CN112702060A (zh) 一种用于铯原子钟的微波幅度控制方法
CN109781341B (zh) 气压传感器的数据校准方法及装置
US10707047B2 (en) Measuring device and measuring method
US9127966B2 (en) Physical quantity detection apparatus and inspection method of physical quantity detection apparatus
EP0499174A1 (en) Automatic skew calibration for multi-channel signal sources
JP2019215226A (ja) 情報処理装置、閉磁路演算方法、および閉磁路演算システム
CN107860973B (zh) 一种应用于频谱分析仪的频响补偿方法及系统
CN113359191A (zh) 一种恒温晶振的实时校正方法和电磁接收机
JP2009038571A (ja) 線形増幅器
Kamenský et al. Correction of ADC errors by additive iterative method with dithering
EP4130661A1 (en) Vibration-type angular velocity sensor
Hawam et al. The resolution uncertainty associated with digital indications revisited: the inclusion of the quantization effect and the impact of noise presence in the estimation process
Lopez et al. Towards decision making under interval uncertainty
CN114136537B (zh) 压力传感器
CN113670336B (zh) 石英挠性加速度计标度因数温度系数补偿特性确定方法
KR20040081272A (ko) 입력신호의 히스토그램 정보에 기초한 주파수 오차검출장치 및 검출방법
US20200021242A1 (en) Oven controlled crystal oscillator device, and frequency compensation method therefor
CN114879114B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210423