CN106814249A - 一种用于光泵磁力仪频率计的频率测量方法 - Google Patents
一种用于光泵磁力仪频率计的频率测量方法 Download PDFInfo
- Publication number
- CN106814249A CN106814249A CN201710113556.5A CN201710113556A CN106814249A CN 106814249 A CN106814249 A CN 106814249A CN 201710113556 A CN201710113556 A CN 201710113556A CN 106814249 A CN106814249 A CN 106814249A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- clock
- meter
- reference 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
- G01R23/10—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into a train of pulses, which are then counted, i.e. converting the signal into a square wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Power Engineering (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明提供了一种用于光泵磁力仪频率计的频率测量方法,所述的频率测量方法在被测时间段的开始与结束点,均由GPS授时信号进行校正,使用了GPS授时信号对频率计的参考频率源进行修正,将频率源误差大幅降低,进而得出由GPS授时信号校正后的准确输出时间,然后再通过计数信息计算精确的频率值,由于GPS授时信号的精度大大超过普通晶振的精度,所以利用本方法实现的频率计,其精度大大超越一般的频率计精度。
Description
技术领域
本发明涉及频率测量技术领域,具体涉及一种用于光泵磁力仪频率计的频率测量方法。
背景技术
高精度磁场检测在空间探测、地面勘探、海洋勘探、人体生物磁场测量等方面发挥着重要的作用,所以一直是各个国家研究的热点领域。由于中国在磁场检测领域起步较晚,仍与发达国家有较大差距。
光泵磁力仪是利用某种元素原子在外磁场作用下的塞曼***、能级跃迁,以及磁共振原理制成的精密磁场测量仪器,已被广泛应用于地球物理勘探、地震监测与预测、航空磁测绘、海事打捞及军事用途等领域。
所述的光泵磁力仪一般可以分为两种,即跟踪式光泵磁力仪和自激式光泵磁力仪。其工作原理为:由光泵作用排列好的原子磁矩,在特定频率的交变电磁场的作用下,又将产生共振吸收作用,打乱原子的排列情况。发生共振吸收现象的电磁场的频率与样品所在点的外磁场强度成一比例关系,故测定这一频率就可以测出外磁场的值。
目前,市面上所使用的测量频率的频率计,其参考频率源都由晶振来提供,时间计数也是由晶振来提供,但是在实际情况中,晶振的精度想当有限,而且受环境温度及供电质量的影响,晶振的精度会有下降,误差较为随机,因此会导致频率计的输出精度不高。
发明内容
本发明的目的在于,为克服光泵磁力仪频率计因采用晶振作为参考频率源,导致频率测量误差较大的技术问题,本发明提供一种用于光泵磁力仪频率计的频率测量方法,利用本发明的方法实现的频率计,其精度大大超越一般的频率计精度。
为实现上述目的,本发明提供的一种用于光泵磁力仪频率计的频率测量方法,包括:
步骤1)选择GPS授时信号作为频率测量的基准信号,并以频率计的晶振提供参考频率信号;
步骤2)在GPS授时信号的整个输出周期内累计获得参考频率信号和被测信号的时钟计数,利用GPS授时信号校正参考频率信号的时钟计数,进一步获得修正后的实际输出时间;
步骤3)将被测信号的时钟计数与实际输出时间之比,获得被测信号的频率值。
作为上述技术方案的进一步改进,所述步骤2)具体包括:
步骤101)在基准信号上升沿来临后,开始监测被测信号的上升沿T1,在上升沿T1来临后,对参考频率信号以及被测信号进行时钟计数,当基准信号下降沿来临后,重新开始监测被测信号的上升沿T2,在上升沿T2来临后,停止对参考频率信号以及被测信号进行时钟计数,且上升沿T2所在的信号不计入计数内;
步骤102)利用下述公式计算获得修正后的实际输出为高电平时参考频率信号的时钟计数T:
其中,R为未修正前实际输出为高电平时参考频率信号的时钟计数,P为GPS基准信号周期内参考频率信号时钟计数的实际输出值,W为GPS基准信号周期内参考频率信号时钟计数的理论输出值;
步骤103)利用步骤102)中的时钟计数T,计算获得修正后的实际输出时间表示为T*N,其中,N表示参考频率信号的时钟周期。
本发明的一种用于光泵磁力仪频率计的频率测量方法优点在于:
本发明提供的频率测量方法,在被测时间段的开始与结束点,均由GPS授时信号进行校正,得出由GPS授时信号校正后的实际时间,然后再通过计数信息计算精确的频率值,由于GPS授时信号的精度大大超过普通晶振的精度,所以利用本方法实现的频率计,其精度大大超越一般的频率计精度。
附图说明
图1为本发明提供的用于光泵磁力仪频率计的频率测量方法流程图。
图2为利用本发明的频率测量方法对各信号进行时钟计数的操作图。
具体实施方式
下面结合附图和实施例对本发明所述的一种用于光泵磁力仪频率计的频率测量方法进行详细说明。
所述的光泵磁力仪一般可以分为两种,即跟踪式光泵磁力仪和自激式光泵磁力仪。以跟踪式激光铯光泵磁力仪为例,铯光泵抽运以铯原子的能级在磁场中产生塞曼效应为基础,利用射频激励铯光谱灯,产生的光谱波长为894nm的D1线和852nm的D2线,当波长为894nm的光持续照射时,D1光作用到铯吸收泡里的铯原子上,将其抽运到F=4能级,完成原子态制备过程。
此时,若有一个外加恒定磁场,F=4的能级上将产生原子磁矩的空间量子化取向,塞曼能级***为9条子能级,这些子能级间的跃迁则为塞曼跃迁,其跃迁频率则为原子磁矩绕外加恒定磁场进动的拉莫尔频率,在忽略其高阶项时,该频率的大小与外加恒定磁场的大小呈线性关系,即
上式中:f0为拉莫尔跃迁频率;B为被测外磁场强度;γ为磁旋比。
对于铯原子,γ/(2π)为3.49857Hz/nT,则被测外磁场强度B和拉莫尔跃迁频率f0有如下关系:
B=f0/3.49857 (2)
故可以通过测量频率实现对磁场的测量。
为了实现拉莫尔跃迁频率的高精度测量,本发明提供了一种用于光泵磁力仪频率计的频率测量方法。该方法使用了GPS授时信号对频率计的参考频率源进行修正,将频率源误差大幅降低。
对于大多数的频率计使用的方法,都是将被测信号与参考信号进行对比得出被测信号的周期时间,具体方法为:在被测信号间隔约20毫秒(此时间可根据实际情况自由设定)的两个上升沿之间的时间段内,通过对参考频率信号进行计数,参考频率由频率计的晶振提供,并对被测信号周期数进行计数,得出被测信号周期所消耗的时间,从而计算出被测信号的频率。
而在本发明提供的频率测量方法中,在被测时间段的开始与结束点,均由GPS授时信号接收芯片接收到的授时信号进行校正,得出由GPS授时信号校正后的实际时间,然后再通过计数信息计算精确的频率值,由于GPS授时信号的精度大大超过普通晶振的精度,所以利用本方法的频率计精度大大超越一般的频率计精度。
如图1所示,本发明提供的一种用于光泵磁力仪频率计的频率测量方法,该方法具体包括:
步骤1)选择GPS授时信号作为频率测量的基准信号,并以频率计的晶振提供参考频率信号;
步骤2)在GPS授时信号的整个输出周期内累计获得参考频率信号和被测信号的时钟计数,利用GPS授时信号校正参考频率信号的时钟计数,进一步获得修正后的实际输出时间;
步骤3)将被测信号的时钟计数与实际输出时间之比,获得被测信号的频率值。
基于上述频率测量方法,所述步骤2)具体包括:
步骤101)在基准信号上升沿来临后,开始监测被测信号的上升沿T1,在上升沿T1来临后,对参考频率信号以及被测信号进行时钟计数,当基准信号下降沿来临后,重新开始监测被测信号的上升沿T2,在上升沿T2来临后,停止对参考频率信号以及被测信号进行时钟计数,且上升沿T2所在的信号不计入计数内;
步骤102)利用下述公式计算获得修正后的实际输出为高电平时参考频率信号的时钟计数T:
其中,R为未修正前实际输出为高电平时参考频率信号的时钟计数,P为GPS基准信号周期内参考频率信号时钟计数的实际输出值,W为GPS基准信号周期内参考频率信号时钟计数的理论输出值;
步骤103)利用步骤102)中的时钟计数T,计算获得修正后的实际输出时间表示为T*N,其中,N表示参考频率信号的时钟周期。
实施例一
在本实施例中,以激光铯光泵磁力仪的输出频率为例,利用本发明提供的方法对Larmor进动频率实现高精度转换。将铯光泵输出频率经过光耦隔离后送入FPGA中,对输入信号进行高精度的计时和计数,并通过GPS接收芯片接收到的GPS时钟信号对参考时钟进行修正,以减少由于环境变化造成的时钟频率误差。具体的操作说明如下:
首先使用高精度晶振给FPGA提供高质量频率源。如图2所示,参考频率使用20M晶振来提供频率源,被测频率信号为30KHz–300KHz之间,基准信号为20Hz信号,其信号源为GPS信号。
在基准信号上升沿来临后,开始监测被测信号上升沿T1,被测信号上升沿来临后,对参考频率信号以及被测信号开始进行时钟计数,当基准信号下降沿来临后,重新开始监测被测信号上升沿T2,T2来临后,停止对参考频率信号以及被测信号的时钟计数,最后一个上升沿T2所在的信号不计入计数内。
由于20Hz信号源为GPS信号,在每个基准信号周期内,20M参考频率的信号数理论上为1M个,但由于环境变化和晶振自身精度限制等原因,晶振频率及长期漂移性能与GPS时钟精度还是有一定差距。所以使用GPS时间信号对参考频率信号的时钟计数进行修正,修正方法为式(3)。
其中,R为未修正前实际输出为1(高电平)时参考频率信号的时钟计数,P为GPS基准信号周期内参考频率时钟计数的实际输出值,T为修正后的实际输出为1时参考频率的时钟数。
然后计算T*50ns即为修正后的实际输出时间。所以实际得到的被测频率为:
其中,C为被测信号的时钟计数。
经过计算可得出,若GPS授时信号准确,晶振输出频率为线性,误差最大值为20Hz/20M,即1ppm。
本发明经过实际验证,使用FPGA对铯光泵输出频率进行高精度检测,并由GPS时间对***进行校时,可得到高精度的磁场测量数据。经过多次验证,本发明的频率测量方法能够达到预期性能,可以为铯光泵Larmor进动频率进行稳定高精度的转换,且使用FPGA芯片为A3P060,能够使得硬件电路结构更为简单,占用面积小,可靠性高,可以方便的集成到其他***中。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种用于光泵磁力仪频率计的频率测量方法,其特征在于,包括:
步骤1)选择GPS授时信号作为频率测量的基准信号,并以频率计的晶振提供参考频率信号;
步骤2)在GPS授时信号的整个输出周期内累计获得参考频率信号和被测信号的时钟计数,利用GPS授时信号校正参考频率信号的时钟计数,进一步获得修正后的实际输出时间;
步骤3)将被测信号的时钟计数与实际输出时间之比,获得被测信号的频率值。
2.根据权利要求1所述的用于光泵磁力仪频率计的频率测量方法,其特征在于,所述步骤2)具体包括:
步骤101)在基准信号上升沿来临后,开始监测被测信号的上升沿T1,在上升沿T1来临后,对参考频率信号以及被测信号进行时钟计数,当基准信号下降沿来临后,重新开始监测被测信号的上升沿T2,在上升沿T2来临后,停止对参考频率信号以及被测信号进行时钟计数,且上升沿T2所在的信号不计入计数内;
步骤102)利用下述公式计算获得修正后的实际输出为高电平时参考频率信号的时钟计数T:
其中,R为未修正前实际输出为高电平时参考频率信号的时钟计数,P为GPS基准信号周期内参考频率信号时钟计数的实际输出值,W为GPS基准信号周期内参考频率信号时钟计数的理论输出值;
步骤103)利用步骤102)中的时钟计数T,计算获得修正后的实际输出时间表示为T*N,其中,N表示参考频率信号的时钟周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710113556.5A CN106814249A (zh) | 2017-02-28 | 2017-02-28 | 一种用于光泵磁力仪频率计的频率测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710113556.5A CN106814249A (zh) | 2017-02-28 | 2017-02-28 | 一种用于光泵磁力仪频率计的频率测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106814249A true CN106814249A (zh) | 2017-06-09 |
Family
ID=59112218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710113556.5A Pending CN106814249A (zh) | 2017-02-28 | 2017-02-28 | 一种用于光泵磁力仪频率计的频率测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106814249A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107870262A (zh) * | 2017-11-01 | 2018-04-03 | 中国科学院上海微***与信息技术研究所 | 一种基于gps授时的高精度光泵磁力仪计频装置及方法 |
CN110297199A (zh) * | 2019-07-05 | 2019-10-01 | 中国地质大学(武汉) | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1862263A (zh) * | 2005-09-19 | 2006-11-15 | 华为技术有限公司 | 一种快速频率测量***及方法 |
CN102998970A (zh) * | 2012-11-30 | 2013-03-27 | 重庆市电力公司江津供电局 | 一种基于crio平台的高精度对时同步方法 |
CN103901271A (zh) * | 2012-12-25 | 2014-07-02 | 东莞市泰斗微电子科技有限公司 | 一种频率测试方法及*** |
CN103913987A (zh) * | 2014-04-26 | 2014-07-09 | 广西电网公司电力科学研究院 | Gps授时***及其获得精确时间基准的方法 |
CN203858351U (zh) * | 2013-11-20 | 2014-10-01 | 武汉理工大学 | 一种铯光泵磁共振信号的高精度检测装置 |
CN104407510A (zh) * | 2014-11-17 | 2015-03-11 | 江汉大学 | 一种授时的方法和装置 |
CN104485956A (zh) * | 2014-12-24 | 2015-04-01 | 天津七六四通信导航技术有限公司 | 高稳定度晶体振荡器的调试方法 |
-
2017
- 2017-02-28 CN CN201710113556.5A patent/CN106814249A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1862263A (zh) * | 2005-09-19 | 2006-11-15 | 华为技术有限公司 | 一种快速频率测量***及方法 |
CN102998970A (zh) * | 2012-11-30 | 2013-03-27 | 重庆市电力公司江津供电局 | 一种基于crio平台的高精度对时同步方法 |
CN103901271A (zh) * | 2012-12-25 | 2014-07-02 | 东莞市泰斗微电子科技有限公司 | 一种频率测试方法及*** |
CN203858351U (zh) * | 2013-11-20 | 2014-10-01 | 武汉理工大学 | 一种铯光泵磁共振信号的高精度检测装置 |
CN103913987A (zh) * | 2014-04-26 | 2014-07-09 | 广西电网公司电力科学研究院 | Gps授时***及其获得精确时间基准的方法 |
CN104407510A (zh) * | 2014-11-17 | 2015-03-11 | 江汉大学 | 一种授时的方法和装置 |
CN104485956A (zh) * | 2014-12-24 | 2015-04-01 | 天津七六四通信导航技术有限公司 | 高稳定度晶体振荡器的调试方法 |
Non-Patent Citations (4)
Title |
---|
张军 等: "基于GPS的时间频率实时校准***", 《现代科学仪器》 * |
张谨 等: "基于FPGA的铯光泵磁力仪频率计设计", 《海洋测绘》 * |
潘宇婷 等: "基于GPS授时的电网频率测量技术的研究", 《电源技术》 * |
罗必露 等: "基于FPGA的全场景试验***主时钟终端频率校准方法", 《电测与仪表》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107870262A (zh) * | 2017-11-01 | 2018-04-03 | 中国科学院上海微***与信息技术研究所 | 一种基于gps授时的高精度光泵磁力仪计频装置及方法 |
CN107870262B (zh) * | 2017-11-01 | 2019-10-15 | 中国科学院上海微***与信息技术研究所 | 一种基于gps授时的高精度光泵磁力仪计频装置及方法 |
CN110297199A (zh) * | 2019-07-05 | 2019-10-01 | 中国地质大学(武汉) | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
CN110297199B (zh) * | 2019-07-05 | 2021-08-27 | 中国地质大学(武汉) | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101692163B (zh) | 一种频率标准远程校准方法及其*** | |
CN103176400B (zh) | 智能电表时钟校准方法 | |
CN107870262B (zh) | 一种基于gps授时的高精度光泵磁力仪计频装置及方法 | |
CN201466800U (zh) | 改进的时间同步测试仪 | |
CN106885998B (zh) | 一种提高铯光泵磁共振信号测频精度和速度的方法和电路 | |
CN105182273A (zh) | 一种检定流水线时间量值的校准*** | |
CN103699001B (zh) | 利用恒温晶体振荡器实现的计时方法及*** | |
CN107272821A (zh) | 实时时钟校准方法及装置、存储介质、电子设备 | |
CN106814249A (zh) | 一种用于光泵磁力仪频率计的频率测量方法 | |
CN105978555A (zh) | 具有温度补偿的实时时钟计时精度修正电路及方法 | |
CN101930211A (zh) | 一种基于gps秒脉冲的时钟源装置及其控制方法 | |
CN105785402A (zh) | 一种用于gnss信号模拟器时延校准的***及方法 | |
CN105137751A (zh) | 一种计量生产调度平台时间量值的校准***及其校准方法 | |
US20160246324A1 (en) | Method and device for clock calibration and corresponding apparatus | |
CN107679260B (zh) | 一种静止轨道遥感卫星星上时统精度确定方法 | |
US20110205107A1 (en) | Positioning device, positioning method and storage medium storing program | |
US9470726B2 (en) | Temperature compensated real-time clock | |
CN103901271B (zh) | 一种频率测试方法及*** | |
US10422809B2 (en) | Electronic apparatus and program | |
CN101833079B (zh) | 全球定位***原始测量值实时转为精确均匀采样的方法 | |
CN104111481B (zh) | 同步时钟相位差测量***和方法 | |
CN103034116A (zh) | 一种提高石英计时器时间计时准确度的方法 | |
CN104460313A (zh) | 用于提供高精度大步长时间同步信号的gps授时装置 | |
CN115616541A (zh) | 一种声纳浮标参考时钟校准方法 | |
CN104614981B (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: 20170609 |