一种基于混频迭代的磁共振信号频率测量方法与流程

1.本发明属于信号测量领域，具体涉及一种基于混频迭代的磁共振信号频率测量方法。


背景技术：

2.基于核磁共振原理的质子磁力仪、overhauser磁力仪被广泛运用于地磁勘探、考古发掘、地磁异常检测等领域。这类磁力仪在进行弱磁测量过程中，通过传感器输出的磁共振信号的频率来反映被测磁场的强度，因此对磁共振信号频率的精确测量，对微弱磁场的精确测量至关重要。
3.由于磁力仪传感器直接输出磁共振信号极其微弱，处于微伏及以下量级，因此从磁共振信号接收到传输的整个过程中会混入各种各样的噪声，而噪声则是影响磁力仪测量的信号频率精度的关键因素之一。
4.面对上述问题，为了提高弱磁测量的精确度，相关领域研究人员提出了诸多频率测频算法，其中常用的方法包括基于快速傅里叶变换（fft）的插值测频法、多通道等精度频率测量算法、多通道插值方法等，这些频率测量算法不同程度的降低了噪声对磁共振信号频率测量精度的影响，获得了相对较高的精度，但在抗干扰能力、频率测量精度等方面仍存在进一步提升的空间。


技术实现要素：

5.本发明的目的在于针对现有技术存在的上述问题，提供一种基于混频迭代的磁共振信号（fid信号）的频率测量方法。
6.本发明的上述目的通过以下技术手段实现：一种基于混频迭代的磁共振信号频率测量方法，包括以下步骤：步骤s1、将磁共振信号检测装置输出的待测磁共振信号分为两路，一路进行模数转换模块采样获得离散fid信号，进一步获得第一路粗测频率f
fft
；另一路转化为方波后经过多通道等精度频率测量方法进行测量，获得第二路粗测频率f
equ
，选取其中一路的粗测频率作为最终粗测频率f
coa
；步骤s2、对离散fid信号进行相位调整，得到初始相位为零的fidz信号，并生成初始频率为f
coa
初始相位为零的余弦信号；步骤s3、将余弦信号与fidz信号进行数字混频、低通滤波、以及算术平均得到混频滤波输出；步骤s4、以预设频率步进f
step
调整余弦信号频率并实施步骤s3直至相邻两次混频滤波输出的符号相异；步骤s5、结合二分法对余弦信号频率进行进一步调整并实施步骤s3直至混频滤波小于设定参数val，计算离散fid信号的最终输出频率f
out
。
7.如上所述步骤s1中第一路粗测频率通过对离散fid信号采用快速傅里叶变换求解
获得。
8.如上所述步骤s1中最终粗测频率f
coa
通过以下步骤获得：将第一路粗测频率f
fft
和第二路粗测频率f
equ
相减并取绝对值，当该绝对值大于傅里叶变换的频谱分辨率时，选取第一路粗测频率f
fft
作为最终粗测频率f
coa
，当该绝对值小于等于傅里叶变换的频谱分辨率时，选取第二路粗测频率f
equ
作为最终粗测频率f
coa
。
9.如上所述频谱分辨率为，为离散fid信号的采样率，为采样点数。
10.如上所述步骤s3中混频滤波输出基于以下公式表示：，其中，时间参数a=-1/t2，meai为混频滤波输出，t为时间，k1为fidz信号的初始幅度，t2为磁共振信号横向弛豫时间，w为余弦信号频率与fidz信号频率的角频率差，w=2π(f
fid-f
std
)，f
fid
为fidz信号频率，f
std
为余弦信号频率。
11.如上所述步骤s3中的低通滤波运算的上限截止频率小于最终粗测频率f
coa
。
12.如上所述步骤s4包括以下步骤：步骤s4.1、预设f
step
为频率步进，将步骤s3中的余弦信号的频率增加f
step
，重复步骤s3获得混频滤波输出；步骤s4.2、记当前实施步骤s3获得混频滤波输出为mea2，对应的余弦信号为当前余弦信号；记前次实施步骤s3获得混频滤波输出为mea1；步骤s4.3：如果混频滤波输出mea1与混频滤波输出mea2符号相同，则将当前余弦信号的频率加f
step
后重复实施步骤s3后返回步骤s4.2；如果混频滤波输出mea1与混频滤波输出mea2符号相异，则将混频滤波输出mea1和混频滤波输出mea2中大于零的记为混频滤波输出mea
l
，将混频滤波输出mea1和混频滤波输出mea2中小于零的记为混频滤波输出mear，混频滤波输出mea
l
和混频滤波输出mear对应的余弦信号频率分别为f
l
与fr，进入步骤s5。
13.如上所述步骤s5包括以下步骤：步骤s5.1：选取余弦信号频率f
l
与余弦信号频率fr的中间值频率fc，将步骤s3中的余弦信号的频率替换为中间值频率fc，实施步骤s3获得的混频滤波输出记为meac并进入步骤s5.2；步骤s5.2：如果混频滤波输出meac与混频滤波输出mea
l
符号相异则将混频滤波输出meac的值赋予混频滤波输出mear，中间值频率fc的值赋予余弦信号频率fr；如果混频滤波输出meac与混频滤波输出mear符号相异则将混频滤波输出meac的值赋予混频滤波输出mea
l
，中间值频率fc的值赋予余弦信号频率f
l
；若混频滤波输出mea
l
和混频滤波输出mear的绝对值都小于设定参数val，则进入步骤s5.3；否则，返回步骤s5.1；步骤s5.3：根据最后一次实施步骤s5.2获得的混频滤波输出mea
l
、混频滤波输出mear、余弦信号频率f
l
、以及余弦信号频率fr，利用线性拟合计算可得离散fid信号的最终输出频率f
out
。
14.如上所述步骤s5.3中最终输出频率f
out
基于以下公式获得：
。
15.本发明相对于现有技术具有以下有益效果：1、在现有频率测量的基础上进一步提高了磁共振信号频率测量的精度。
16.2、基于本发明的混频滤波特性，可有效提高对低信噪比磁共振信号的频率测量的稳定性。
17.3、基于提高磁共振信号频率测量的精度，可以使质子磁力仪、overhauser磁力仪的磁场测量绝对精度与灵敏度得到较大提升，从而更好地运用于地质勘探、地质灾害监测等领域。
18.4、为其他涉及单一信号频率测量的相关领域提供参考思路。
附图说明
19.图1为本发明的流程框图。
20.图2为角频率差w与混频滤波输出mea的函数关系图。
21.图3为离散fid信号示意图。
22.图4为余弦信号std示例图。
具体实施方式
23.为了便于本领域普通技术人员理解和实施本发明。下面结合实施实例对本发明做进一步详细说明，应该理解此处描述的实施示例仅用于说明和解释本发明，并不限于限定本发明。
24.一种基于混频迭代的磁共振信号频率测量方法，具体包括以下步骤：步骤s1、将磁共振信号检测装置输出的待测磁共振信号分为两路，一路待测磁共振信号利用模数转换模块adc对信号进行数字化，采样率f
acq
为100khz，采样点数p
nm
为100000个点，得到离散fid信号，如图3所示，同时采用快速傅里叶变换求解待测磁共振信号频率，记为第一路粗测频率f
fft
；另一路待测磁共振信号经过比较器处理转化为方波，采用多通道等精度频率测量方法进行测量，计算得出第二路粗测频率f
equ
。将第一路粗测频率f
fft
和第二路粗测频率f
equ
相减并取绝对值，当该绝对值大于傅里叶变换的频谱分辨率（）时，选取第一路粗测频率f
fft
作为最终粗测频率f
coa
，当该绝对值小于等于傅里叶变换的频谱分辨率（）时，选取第二路粗测频率f
equ
作为最终粗测频率f
coa
。
25.当（1）当（2）多通道等精度频率测量方法基于以下文章：质子磁力仪fid信号的高精度频率测量算法（ahigh-precisionfrequencymeasurementalgorithmforfidsignalofprotonmagnetometer），ieee仪器与测量学报（ieeetransactionsoninstrumentationandmeasurement），第65卷，第4期，2016年4月。
26.本实施例中测得第一路粗测频率f
fft
为2101hz，第二路粗测频率f
equ
为2174.92hz，|f
fft-f
equ
|=73.92hz；因此令最终粗测频率f
coa
等于第一路粗测频率f
fft
,即为2101hz。
27.步骤s2、将步骤s1中得到的离散fid信号进行相位调整，使其成为初始相位为零的正弦震荡衰减信号，得到调整后的fidz信号，生成初始频率为f
coa
、初始相位为零的余弦信号。
28.记余弦信号为std：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）其中，t为时间，f
std
为余弦信号频率。
29.fidz信号基于以下公式：
ꢀꢀꢀꢀꢀꢀꢀꢀ
（4）其中，t为时间，k1为fidz信号的初始幅度，t2为磁共振信号横向弛豫时间，f
fid
为fidz信号频率，noise为fidz信号中的噪声成分。
30.步骤s3、将余弦信号与fidz信号进行数字混频、低通滤波、以及算术平均得到混频滤波输出mea，低通滤波运算的上限截止频率应当小于最终粗测频率f
coa
，滤除其中包含的频率大于最终粗测频率f
coa
的信号成分。
31.上述混频滤波输出mea基于以下公式表示：
ꢀꢀ
（5）其中，时间参数a=-1/t2，w为余弦信号频率与fidz信号频率的角频率差，且w=2π(f
fid-f
std
)，依据上式可知，有且只有在w等于0时，有mea等于0，即表明只有余弦信号频率f
std
与待测的fidz信号频率f
fid
相等时，混频滤波运算的输出mea为0。
32.因此，通过对混频滤波运算的连续迭代，不断调整余弦信号的频率，找到令混频滤波输出mea等于0或者无限接近0的频率，即为所被测的离散fid信号的真实频率。
33.步骤s4、以预设参数f
step
调整余弦信号频率并实施步骤s3直至相邻两次混频滤波输出的符号相异，具体包括以下步骤：步骤s4.1、预设f
step
为频率步进，本实施例中，f
step
等于0.5hz，将步骤s3中的余弦信号的频率增加f
step
，即余弦信号的频率为2101.5hz，重复步骤s3获得混频滤波输出。
34.步骤s4.2、记当前实施步骤s3获得混频滤波输出为mea2，对应的余弦信号为当前余弦信号；记前次实施步骤s3获得混频滤波输出为mea1，mea1的大小为0.106922，mea2的大小为-0.326397，步骤s4.3、判断mea1与mea2是否符号相异（一正一负则为相异）；如果符号相同，则将当前余弦信号的频率加f
step
后重复实施步骤s3后返回步骤s4.2；如果符号相异，则表明待测频率值存在于两个符号相异的混频滤波运算输出所对应的输入频率之间，则将混频滤波输出mea1和混频滤波输出mea2中大于零的记为混频滤波输出mea
l
，将混频滤波输出mea1和混频滤波输出mea2中小于零的记为混频滤波输出mear，混频滤波输出mea
l
和混频滤波输出mear对应的余弦信号频率分别为f
l
与fr，进入步骤s5。
35.在本实施例中，混频滤波输出mea1与混频滤波输出mea2异号，因此，令混频滤波输出mea
l
等于混频滤波输出mea1，令混频滤波输出mear等于混频滤波输出mea2，余弦信号频率f
l
与余弦信号频率fr对应的频率分别为2101hz与2101.5hz。
36.步骤s5、结合二分法对余弦信号频率进行进一步调整并实施步骤s3直至混频滤波小于设定参数val，计算离散fid信号的最终输出频率f
out
，具体包括以下步骤：步骤s5.1、选取余弦信号频率f
l
与余弦信号频率fr的中间值频率fc，将步骤s3中的余弦信号的频率替换为中间值频率fc，实施步骤s3获得的混频滤波输出记为meac并进入步骤s5.2。
37.步骤s5.2、判断混频滤波输出meac与混频滤波输出mea
l
、混频滤波输出mear的符号关系，如果混频滤波输出meac与混频滤波输出mea
l
符号相异则将混频滤波输出meac的值赋予混频滤波输出mear，中间值频率fc的值赋予余弦信号频率fr；如果混频滤波输出meac与混频滤波输出mear符号相异则将混频滤波输出meac的值赋予混频滤波输出mea
l
，中间值频率fc的值赋予余弦信号频率f
l
；若混频滤波输出mea
l
和混频滤波输出mear的绝对值都小于设定参数val，本实施例中，val等于0.001，则进入步骤s5.3；否则，返回步骤s5.1；在步骤s5.2中，通过反复迭代每次将零点所在的区间减小至上一次迭代的一半，以达到逼近过零点的目的。其中设定参数val的大小，与最终测量精度和计算量相关，其值的选取应根据实际情况而定。
38.步骤s5.3、根据最后一次实施步骤s5.2获得的混频滤波输出mea
l
、混频滤波输出mear、余弦信号频率f
l
、以及余弦信号频率fr，利用线性拟合计算可得离散fid信号的最终输出频率f
out
，如下式所示：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（6）对于同一个离散fid信号，t、k1、以及a都为常数，即表明当w=0时，d(mea)/dw等于一个常数，
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）由此可知在过零点附近为线性函数，因此采用线性拟合的方式求得过零点，线性拟合表达式如式7所示。
39.本实施例中，第10次获得的mea
l
与mear的绝对值都小于0.001，计算过程中个参数值如下表1所示：表1迭代过程参数表
40.如表1所示，完成第9次迭代运算后mea
l
与mear的绝对值满足要求，此时对应频率f
l
与fr分别为2101.105469hz与2101.106445hz。
41.将f
l
、fr、mea
l
、mear的值代入上式，可得最终计算频率2101.105641hz。
42.需要指出的是，本发明中所描述的具体实施例仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例作各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或超越所附权利要求书所定义的范围。

技术特征：
1.一种基于混频迭代的磁共振信号频率测量方法，包括以下步骤：步骤s1、将磁共振信号检测装置输出的待测磁共振信号分为两路，一路进行模数转换模块采样获得离散fid信号，进一步获得第一路粗测频率f
fft
；另一路转化为方波后经过多通道等精度频率测量方法进行测量，获得第二路粗测频率f
equ
，选取其中一路的粗测频率作为最终粗测频率f
coa
；步骤s2、对离散fid信号进行相位调整，得到初始相位为零的fidz信号，并生成初始频率为f
coa
初始相位为零的余弦信号；步骤s3、将余弦信号与fidz信号进行数字混频、低通滤波、以及算术平均得到混频滤波输出；步骤s4、以预设频率步进f
step
调整余弦信号频率并实施步骤s3直至相邻两次混频滤波输出的符号相异；步骤s5、结合二分法对余弦信号频率进行进一步调整并实施步骤s3直至混频滤波小于设定参数val，计算离散fid信号的最终输出频率f
out
。2.根据权利要求1所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s1中第一路粗测频率通过对离散fid信号采用快速傅里叶变换求解获得。3.根据权利要求1所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s1中最终粗测频率f
coa
通过以下步骤获得：将第一路粗测频率f
fft
和第二路粗测频率f
equ
相减并取绝对值，当该绝对值大于傅里叶变换的频谱分辨率时，选取第一路粗测频率f
fft
作为最终粗测频率f
coa
，当该绝对值小于等于傅里叶变换的频谱分辨率时，选取第二路粗测频率f
equ
作为最终粗测频率f
coa
。4.根据权利要求3所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述频谱分辨率为，为离散fid信号的采样率，为采样点数。5.根据权利要求4所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s3中混频滤波输出基于以下公式表示：，其中，时间参数a=-1/t2，meai为混频滤波输出，t为时间，k1为fidz信号的初始幅度，t2为磁共振信号横向弛豫时间，w为余弦信号频率与fidz信号频率的角频率差，w=2π(f
fid-f
std
)，f
fid
为fidz信号频率，f
std
为余弦信号频率。6.根据权利要求1所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s3中的低通滤波运算的上限截止频率小于最终粗测频率f
coa
。7.根据权利要求1所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s4包括以下步骤：步骤s4.1、预设f
step
为频率步进，将步骤s3中的余弦信号的频率增加f
step
，重复步骤s3获得混频滤波输出；步骤s4.2、记当前实施步骤s3获得混频滤波输出为mea2，对应的余弦信号为当前余弦信号；记前次实施步骤s3获得混频滤波输出为mea1；步骤s4.3：如果混频滤波输出mea1与混频滤波输出mea2符号相同，则将当前余弦信号的
频率加f
step
后重复实施步骤s3后返回步骤s4.2；如果混频滤波输出mea1与混频滤波输出mea2符号相异，则将混频滤波输出mea1和混频滤波输出mea2中大于零的记为混频滤波输出mea
l
，将混频滤波输出mea1和混频滤波输出mea2中小于零的记为混频滤波输出mea
r
，混频滤波输出mea
l
和混频滤波输出mea
r
对应的余弦信号频率分别为f
l
与f
r
，进入步骤s5。8.根据权利要求7所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s5包括以下步骤：步骤s5.1：选取余弦信号频率f
l
与余弦信号频率f
r
的中间值频率f
c
，将步骤s3中的余弦信号的频率替换为中间值频率f
c
，实施步骤s3获得的混频滤波输出记为mea
c
并进入步骤s5.2；步骤s5.2：如果混频滤波输出mea
c
与混频滤波输出mea
l
符号相异则将混频滤波输出mea
c
的值赋予混频滤波输出mea
r
，中间值频率f
c
的值赋予余弦信号频率f
r
；如果混频滤波输出mea
c
与混频滤波输出mea
r
符号相异则将混频滤波输出mea
c
的值赋予混频滤波输出mea
l
，中间值频率f
c
的值赋予余弦信号频率f
l
；若混频滤波输出mea
l
和混频滤波输出mea
r
的绝对值都小于设定参数val，则进入步骤s5.3；否则，返回步骤s5.1；步骤s5.3：根据最后一次实施步骤s5.2获得的混频滤波输出mea
l
、混频滤波输出mea
r
、余弦信号频率f
l
、以及余弦信号频率f
r
，利用线性拟合计算可得离散fid信号的最终输出频率f
out
。9.根据权利要求8所述一种基于混频迭代的磁共振信号频率测量方法，其特征在于，所述步骤s5.3中最终输出频率f
out
基于以下公式获得：。

技术总结
本发明公开了一种基于混频迭代的磁共振信号频率测量方法，利用等精度频率测量方法与快速傅里叶变换相结合粗测信号，将余弦信号与FIDz信号进行数字混频、低通滤波、以及算术平均，调整余弦信号频率，迭代寻找混频滤波的过零点所在区间，采用二分法不断缩小过零点区间逼近过零点，最终采用线性拟合的方式求得最终输出频率。本发明基于混频滤波特性，可有效提高对低信噪比磁共振信号的频率测量的稳定性。高对低信噪比磁共振信号的频率测量的稳定性。高对低信噪比磁共振信号的频率测量的稳定性。


技术研发人员：刘朝阳 匡澳元 陈俊飞 李文帅 冯继文
受保护的技术使用者：湖北光谷实验室
技术研发日：2023.07.19
技术公布日：2023/8/24