一种超极化惰性气体FID信号共振频率的精确测量方法

本发明旨在提出适用于原子磁强计领域的一种超极化惰性气体fid信号共振频率的精确测量方法，属于原子分子及光物理、信号检测与分析领域。

背景技术：

1、随着量子技术和磁场传感的迅速发展，磁场测量的精度和灵敏度不断提高，极具潜力的原子磁强计得到了迅猛的发展。原子磁强计不同于传统磁强计，其依赖于原子的自旋，通过光学手段操纵原子，进而根据原子的自旋状态探测环境磁场。其中，fsp磁强计由光泵磁强计发展而来，通过光学探测原子自旋极化的自由自旋进动获得磁场信息，而其中典型的一种信号即自由感应衰减(fid)信号。该类磁强计避免了反馈驱动类磁强计的系统读出误差，保留了通过测量频率进而测量磁场机制的优点，是一种极具发展潜力的原子磁强计。

2、自由感应衰减(fid，free induction decay)信号是核磁共振或电子顺磁共振领域的由围绕磁场方向的非平衡核自旋磁化产生的可观察共振信号，假设原子系综的磁化矢量被极化至z轴，在外磁场的作用下绕z轴做拉莫尔进动，当对系综施加一个频率接近原子拉莫尔进动频率的射频脉冲时，即可产生自由感应衰减。通过检测样品周围检测线圈中的电压，进而数字化，即可得到实验中的fid信号。fid信号常被用来测量原子磁强计中样品的弛豫或环境磁场、磁场梯度等，j.murday通过分析fid信号测量了磁场梯度，并分析了系统误差主要来源于共振条件、有限射频脉冲的时间延迟、梯度和圆柱轴的非正交性等；e.breschi等人利用fid信号校准了三维磁线圈，并测量了残余磁场的三轴分量。

3、fid信号是一种振荡衰减的信号，其频率反映了原子的拉莫尔进动频率，进而反映了环境磁场的大小，由于此处的环境磁场通常为原子磁强计中屏蔽后的磁场，幅值较小，常在nt量级，因此对频率分辨率的要求较高。将信号直接傅里叶变换得到的频率分辨率较低，由于存在栅栏效应，频谱中两个相邻点之间的频率无法测得，直接影响了磁场测量的精度，因此对fid信号的频谱进行细化估计是必要的。

4、为了解决fid信号的频谱分辨率不足，本发明提出一种适用于屏蔽环境下的超极化惰性气体fid信号共振频率的精确测量方法。本方法针对有限时长的超极化惰性气体fid时域信号引起的频谱分辨率较低的问题，提出了局部dft(discrete fourier transform，离散傅里叶变换)和抛物线插值的细化估计方法，该方法减小了栅栏效应带来的误差，能够进一步提升有限时长超极化惰性气体fid信号的频率分辨率，对提升磁场测量的精度具有重要意义。

技术实现思路

1、本发明目的：为了设计一种适合于有限时常、固定采样率的系统的频谱细化方法，本发明聚焦于超极化惰性气体的fid信号频率提取，提出了一种超极化惰性气体fid信号共振频率的精确测量方法，属于一种固定采样点、高精度、快速的频率细化估计方法，本方法主要包括两个部分：频谱分辨率细化、频率估计，该方法有助于实现精确获取fid信号的频率，为磁场真值的准确获取提供保障。

2、本发明的技术解决方案如下：

3、一种超极化惰性气体fid信号共振频率的精确测量方法，其特征在于，包括以下步骤：

4、步骤1，从超极化惰性气体fid信号采集实验装置得到自由感应衰减fid信号，对fid信号进行fft快速傅里叶变换得到fft频谱，从fft频谱中确定目标频段；

5、步骤2，对目标频段进行局部放大，找到幅值最大的两个点，分别为k1和k2；

6、步骤3，利用dft离散傅里叶变换在k1和k2之间进行频率细化；

7、步骤4，找到细化后幅值最大的三个点，进行抛物线拟合；

8、步骤5，算出抛物线导数为零的横坐标点，即为超极化惰性气体fid信号共振频率的精确估计解。

9、所述步骤4中包括以下抛物线公式：

10、y＝ax2+bx+c

11、其中y表示纵坐标幅值，x表示横坐标频率，a、b、c表示抛物线表达式系数，通过细化后幅值最大的三个点p1(x1,y1),p2(x2,y2),p3(x3,y3)建立以下联立方程求解a、b和c：

12、y1＝ax12+bx1+c

13、y2＝ax22+bx2+c

14、y3＝ax32+bx3+c

15、由此得到拟合后一条开口向下的抛物线函数y＝ax2+bx+c。

16、所述步骤5中包括算式：y′＝2ax0+b＝0

17、y’表示抛物线函数y的导数，x0表示抛物线导数为零的横坐标点，x0即为超极化惰性气体fid信号共振频率的精确估计解。

18、定义频谱提高的分辨率为横坐标x0与细化前幅值最大点k1对应的横坐标f(k1)之间的差，即δf′＝|x0-f(k1)|，用于衡量频谱细化和频率估计的效果，当该值较小时，证明频率真值在细化前位于谱线附近，当该值较大时，证明频率真值在细化前距离谱线较远。

19、所述步骤1中的目标频段在0～10hz之间。

20、所述步骤1中的目标频段为fft频谱的频谱峰值区域。

21、所述步骤2中的k1和k2为相邻峰值的频率点，且k1＜k2。

22、所述步骤3中包括以下dft的公式：

23、

24、其中x(k)表示频域信号，n为采样点数，n为采样点序号，x(n)表示时域信号，e为自然常数，j为虚数单位，k为频率点序号；

25、采用连续代替离散思想，在x(k1)、x(k2)之间进行插值细化，具体内插的值为，

26、

27、其中d为细化倍数，由于k1、k2为相邻点，k1在左侧，因此k1+1＝k2，细化后的点对应的标号为k1＝k10,k11,k12,k13,...,k1d＝k2，各点之间的频率差为δf/d，其中δf为细化前的频率分辨率，公式为，

28、δf＝fs/n

29、其中fs为采样率，n为采样点数。

30、所述超极化惰性气体fid信号采集实验装置包括沿z轴正向依次串联的泵浦光发射器、第一半波片、第一扩束镜、第二半波片、第一楔形块、第二楔形块、第一格兰-汤普逊棱镜、第一1/4波片和原子气室，所述原子气室外围向外依次设置有线圈和磁屏蔽桶。

31、所述超极化惰性气体fid信号采集实验装置包括沿x轴负向依次串联的探测光发射器、第二扩束镜、第三半波片、第三楔形块、第四楔形块、第二格兰-汤普逊棱镜、原子气室和第四半波片，所述第四半波片通过第一分束器的透射侧连接第二光电二极管进行信号输出，第一分束器的反射侧连接第一光电二极管进行信号输出。

32、本发明的技术效果如下：本发明一种超极化惰性气体fid信号共振频率的精确测量方法，属于原子分子及光物理、信号检测与分析领域，对于超极化惰性气体的fid信号，由于信号本身随时间振荡衰减，获取其频谱需要在有限时间内进行采样，而采样率无法无限提高，因此得到的频谱分辨率受限，进而导致磁场测量的精度受影响。本发明通过先后进行细化频谱和频率估计，显著提高了频率分辨率，对提高磁场测量精度具有重要意义。

33、本发明的特点如下：(1)提供了一种适用于低频信号，高分辨率频谱的估计方法；(2)本方法有助于利用fid进行磁场的高精度测量；(3)精确的fid信号频率和磁场测量为原子磁强计的残余磁场测量提供保障。