CN114460506A - 基于变参数控制的磁共振信号相位闭环控制方法及*** - Google Patents

基于变参数控制的磁共振信号相位闭环控制方法及*** Download PDF

Info

Publication number
CN114460506A
CN114460506A CN202111585410.3A CN202111585410A CN114460506A CN 114460506 A CN114460506 A CN 114460506A CN 202111585410 A CN202111585410 A CN 202111585410A CN 114460506 A CN114460506 A CN 114460506A
Authority
CN
China
Prior art keywords
signal
magnetic resonance
digital
phase
resonance signal
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.)
Granted
Application number
CN202111585410.3A
Other languages
English (en)
Other versions
CN114460506B (zh
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.)
Beijing Automation Control Equipment Institute BACEI
Original Assignee
Beijing Automation Control Equipment Institute BACEI
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 Beijing Automation Control Equipment Institute BACEI filed Critical Beijing Automation Control Equipment Institute BACEI
Priority to CN202111585410.3A priority Critical patent/CN114460506B/zh
Publication of CN114460506A publication Critical patent/CN114460506A/zh
Application granted granted Critical
Publication of CN114460506B publication Critical patent/CN114460506B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/032Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种基于变参数控制的磁共振信号相位闭环控制方法及***,包括:基于DDS技术生成数字参考信号和磁场激励信号;获取原子气室的磁共振信号;对数字参考信号和磁共振信号分别进行滤波处理以及对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整;计算获取磁共振信号和数字参考信号之间的相位差;对磁共振信号和数字参考信号之间的相位差进行相位补偿;根据磁共振信号和磁场激励信号之间的相位差与期望相位差之间的差值调整数字参考信号和磁场激励信号的频率完成磁共振信号相位闭环控制。应用本发明的技术方案,以解决现有技术中相位解调方法的噪声干扰大且计算结果随环境温度变化具有较大的漂移的技术问题。

Description

基于变参数控制的磁共振信号相位闭环控制方法及***
技术领域
本发明涉及原子磁强计相位控制技术领域,尤其涉及一种基于变参数控制的磁共振信号相位闭环控制方法及***。
背景技术
碱金属原子磁强计利用电子自旋或核自旋在磁场中的进动测量磁场,具有精度高、体积小等优势,主要用于地球磁场测量、脑磁心磁等生物弱磁测量以及水下、水面和航空目标探测等方面,在国民经济建设和国防领域具有重要意思。碱金属原子磁强计测量地球磁场时,要求原子气室内的原子工作在共振状态,此时激励信号和磁共振信号频率相同,相位差为常值,原子磁强计具有最高的灵敏度和稳定性。因此,为提高磁场测量稳定性,首先要实现磁共振信号相位的高精度闭环控制。
传统的模拟相位解调控制方式是通过基于模拟乘法器的方法实现。由于模拟乘法器是全带宽范围内信号的处理,非有用信号频段的噪声会直接影响待测信号相位的解调计算,带来噪声干扰;而且模拟乘法器的计算结果随环境温度变化具有较大的漂移,难以满足原子磁强计工程环境应用需求。
发明内容
本发明提供了一种基于变参数控制的磁共振信号相位闭环控制方法及***,能够解决现有技术中相位解调方法的噪声干扰大且计算结果随环境温度变化具有较大的漂移的技术问题。
根据本发明的一方面,提供了一种基于变参数控制的磁共振信号相位闭环控制方法,基于变参数控制的磁共振信号相位闭环控制方法包括:基于DDS技术生成数字参考信号和磁场激励信号;根据磁场激励信号对原子气室进行磁场激励,获取原子气室的磁共振信号;对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,实现对数字参考信号和磁共振信号的变参数控制;根据滤波参考信号和滤波磁共振信号计算获取磁共振信号和数字参考信号之间的相位差;对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差;将磁共振信号和磁场激励信号之间的相位差与期望相位差相比较,当磁共振信号和磁场激励信号之间的相位差与期望相位差不相等时,调整数字参考信号和磁场激励信号的频率,直至磁共振信号和磁场激励信号之间的相位差等于期望相位差,完成磁共振信号相位闭环控制。
进一步地,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差具体包括:根据滤波参考信号和滤波磁共振信号进行数字互相关运算和数字自相关运算以获取第一自相关运算量、第二自相关运算量和互相关运算量,根据第一自相关运算量、第二自相关运算量和互相关运算量计算获取磁共振信号和参考信号的相位差。
进一步地,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差具体包括:计算获取***相位滞后和相位漂移,根据***相位滞后和相位漂移对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差。
进一步地,基于DDS技术生成数字参考信号和磁场激励信号具体包括:利用DDS内核单元生成第一数字信号和第二数字信号,将第一数字信号作为数字参考信号,对第二数字信号进行数模转换以获取磁场激励信号。
进一步地,数字参考信号和磁共振信号的滤波处理结果均可根据YN=a0XN-2+a1XN-1+a2XN+b1YN-1+b2YN-2来获取,其中,XN为当前时刻输入值,XN-1为前一时刻输入值,XN-2为前两时刻输入值,YN-1为前一时刻滤波结果,YN-2为前两时刻滤波结果,a0、a1、a2、b1和b2为滤波系数。
进一步地,第一自相关运算量XAA、第二自相关运算量XBB和互相关运算量XAB可根据
Figure BDA0003427702400000031
计算获取,其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。
进一步地,磁共振信号和数字参考信号之间的相位差Δθ可根据
Figure BDA0003427702400000032
计算获取。
进一步地,磁共振信号和磁场激励信号之间的相位差θ可根据
Figure BDA0003427702400000033
计算获取,其中,
Figure BDA0003427702400000034
为***相位滞后,
Figure BDA0003427702400000035
为相位漂移。
根据本发明的又一方面,提供了一种基于变参数控制的磁共振信号相位闭环控制***,磁共振信号相位闭环控制***使用如上所述的磁共振信号相位闭环控制方法进行磁共振信号相位闭环控制。
进一步地,磁共振信号相位闭环控制***包括DDS内核单元、变参数控制单元、模数转换器、数模转换器、数字相关运算单元、相位补偿单元和PID闭环控制单元,DDS内核单元用于生成第一数字信号和第二数字信号,第一数字信号为数字参考信号,数模转换器用于对第二数字信号进行数模转换以获取磁场激励信号,模数转换器用于将模拟量的磁共振信号转换为数字量的磁共振信号,滤波及变参数控制单元用于对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号以及根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,数字相关运算用于对滤波参考信号和滤波磁共振信号进行数字自相关运算和数字互相关运算以获取磁共振信号和数字参考信号之间的相位差,相位补偿单元用于对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差,PID闭环控制单元用于根据磁共振信号和磁场激励信号之间的相位差与期望相位差之间的大小调整数字参考信号和磁场激励信号的输出频率以完成磁共振信号相位闭环控制。
应用本发明的技术方案,提供了一种基于变参数控制的磁共振信号相位闭环控制方法,该方法基于DDS技术和数字相关运算,可实现磁共振相位的数字提取;变参数控制可有效控制有用信号的处理带宽,降低频带内其他噪声的干扰,提高磁场解算精度;而相位补偿可大幅降低***固有的相位滞后、相位漂移,提高磁场测量稳定性,保证***环境适应性。因此,本发明所提供的磁共振信号相位闭环控制方法与现有技术相比,实现了磁共振信号相位闭环控制,降低了频带内其他噪声的干扰,提高了磁场解算精度;且大幅降低***固有的相位滞后、相位漂移和温度漂移,提高了磁场测量稳定性、***的相位解算精度、稳定性及工程环境适应能力,保证了***的环境适应性。
附图说明
所包括的附图用来提供对本发明实施例的进一步的理解,其构成了说明书的一部分,用于例示本发明的实施例,并与文字描述一起来阐释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了根据本发明的具体实施例提供的基于变参数控制的磁共振信号相位闭环控制***的结构示意图。
其中,上述附图包括以下附图标记:
10、DDS内核单元;20、变参数控制单元;30、模数转换器;40、数模转换器;50、数字相关运算单元;60、相位补偿单元;70、PID闭环控制单元。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的,决不作为对本发明及其应用或使用的任何限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。同时,应当明白,为了便于描述,附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为授权说明书的一部分。在这里示出和讨论的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制。因此,示例性实施例的其它示例可以具有不同的值。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步讨论。
如图1所示,根据本发明的具体实施例提供了一种基于变参数控制的磁共振信号相位闭环控制方法,该基于变参数控制的磁共振信号相位闭环控制方法包括:基于DDS技术生成数字参考信号和磁场激励信号;根据磁场激励信号对原子气室进行磁场激励,获取原子气室的磁共振信号;对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,实现对数字参考信号和磁共振信号的变参数控制;根据滤波参考信号和滤波磁共振信号计算获取磁共振信号和数字参考信号之间的相位差;对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差;将磁共振信号和磁场激励信号之间的相位差与期望相位差相比较,当磁共振信号和磁场激励信号之间的相位差与期望相位差不相等时,调整数字参考信号和磁场激励信号的频率,直至磁共振信号和磁场激励信号之间的相位差等于期望相位差,完成磁共振信号相位闭环控制。
应用此种配置方式,提供了一种基于变参数控制的磁共振信号相位闭环控制方法,该方法基于DDS技术和数字相关运算,可实现磁共振相位的数字提取;变参数控制可有效控制有用信号的处理带宽,降低频带内其他噪声的干扰,提高磁场解算精度;而相位补偿可大幅降低***固有的相位滞后、相位漂移,提高磁场测量稳定性,保证***环境适应性。因此,本发明所提供的磁共振信号相位闭环控制方法与现有技术相比,实现了磁共振信号相位闭环控制,降低了频带内其他噪声的干扰,提高了磁场解算精度;且大幅降低***固有的相位滞后、相位漂移和温度漂移,提高了磁场测量稳定性、***的相位解算精度、稳定性及工程环境适应能力,保证了***的环境适应性。
在本发明中,为了实现磁共振信号相位闭环控制,首先需要基于DDS技术生成数字参考信号和磁场激励信号。在本发明中,基于DDS技术生成数字参考信号和磁场激励信号具体包括:利用DDS内核单元10生成第一数字信号和第二数字信号,将第一数字信号作为数字参考信号,对第二数字信号进行数模转换以获取磁场激励信号。其中,数字参考信号XA可根据
Figure BDA0003427702400000071
计算获取,其中,A为数字参考信号幅度,ω为数字参考信号频率,t为时间,θA为数字参考信号相位,
Figure BDA0003427702400000072
为数字参考信号噪声。
进一步地,在获取了数字参考信号和磁场激励信号之后,即可根据磁场激励信号对原子气室进行磁场激励,获取原子气室的磁共振信号。磁共振信号XB与数字参考信号XA同频,存在相位差,磁共振信号XB可根据
Figure BDA0003427702400000073
计算获取,其中,B为磁场激励信号幅度,ω为磁场激励信号频率,磁场激励信号频率与数字参考信号频率相同,t为时间,θB为磁场激励信号相位,
Figure BDA0003427702400000074
为磁场激励信号噪声。
在获取了原子气室的磁共振信号之后,为了降低数字参考信号XA和磁共振信号XB内的杂频噪声干扰,即可对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,实现对数字参考信号和磁共振信号的变参数控制。在本发明中,数字参考信号和磁共振信号的滤波处理结果均可根据YN=a0XN-2+a1XN-1+a2XN+b1YN-1+b2YN-2来获取,其中,XN为当前时刻输入值,XN-1为前一时刻输入值,XN-2为前两时刻输入值,YN-1为前一时刻滤波结果,YN-2为前两时刻滤波结果,a0、a1、a2、b1、b2为滤波系数。
在本发明中,在对数字参考信号和磁共振信号分别进行滤波处理后,滤波参考信号XAPF可表示为:XAPE=Asin(ωt+θAPE),滤波磁共振信号XBPF可表示为:XBPE=Bsin(ωt+θBPE),其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。
进一步地,在本发明中,在完成了数字参考信号和所述磁共振信号的变参数控制之后,即可根据滤波参考信号和滤波磁共振信号计算获取磁共振信号和数字参考信号之间的相位差。在本发明中,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差具体包括:根据滤波参考信号和滤波磁共振信号进行数字互相关运算和数字自相关运算以获取第一自相关运算量、第二自相关运算量和互相关运算量,根据第一自相关运算量、第二自相关运算量和互相关运算量计算获取磁共振信号和参考信号的相位差。
作为本发明的一个具体实施例,第一自相关运算量XAA、第二自相关运算量XBB和互相关运算量XAB可根据
Figure BDA0003427702400000091
计算获取,其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。磁共振信号和数字参考信号之间的相位差Δθ可根据
Figure BDA0003427702400000092
计算获取。
进一步地,在获取了磁共振信号和数字参考信号之间的相位差之后,即可对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和所述磁场激励信号之间的相位差。在本发明中,***需要求解的为磁共振信号与磁场激励信号之间的相位差,而数字参考信号与磁场激励信号之间存在相位滞后,且该滞后随温度漂移,因此需要对相位进行补偿。
在本发明中,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差具体包括:计算获取***相位滞后和相位漂移,根据***相位滞后和相位漂移对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差。具体地,磁共振信号和磁场激励信号之间的相位差θ可根据
Figure BDA0003427702400000093
计算获取,其中,
Figure BDA0003427702400000094
为***相位滞后,
Figure BDA0003427702400000095
为相位漂移。
进一步地,在获取了磁共振信号和磁场激励信号之间的相位差之后,即可将磁共振信号和磁场激励信号之间的相位差与期望相位差相比较,当磁共振信号和磁场激励信号之间的相位差与期望相位差不相等时,调整数字参考信号和磁场激励信号的频率,直至磁共振信号和磁场激励信号之间的相位差等于期望相位差,完成磁共振信号相位闭环控制。
根据本发明的另一方面,提供了一种基于变参数控制的磁共振信号相位闭环控制***,该磁共振信号相位闭环控制***使用如上所述的磁共振信号相位闭环控制方法进行磁共振信号相位闭环控制。
应用此种配置方式,提供了一种基于变参数控制的磁共振信号相位闭环控制***,该***基于DDS技术和数字相关运算,可实现磁共振相位的数字提取;变参数控制可有效控制有用信号的处理带宽,降低频带内其他噪声的干扰,提高磁场解算精度;而相位补偿可大幅降低***固有的相位滞后、相位漂移,提高磁场测量稳定性,保证***环境适应性。因此,本发明所提供的磁共振信号相位闭环控制***与现有技术相比,实现了磁共振信号相位闭环控制,降低了频带内其他噪声的干扰,提高了磁场解算精度;且大幅降低***固有的相位滞后、相位漂移和温度漂移,提高了磁场测量稳定性、***的相位解算精度、稳定性及工程环境适应能力,保证了***的环境适应性。
进一步地,在本发明中,磁共振信号相位闭环控制***包括DDS内核单元10、变参数控制单元20、模数转换器30、数模转换器40、数字相关运算单元50、相位补偿单元60和PID闭环控制单元70,DDS内核单元10用于生成第一数字信号和第二数字信号,第一数字信号为数字参考信号,数模转换器40用于对第二数字信号进行数模转换以获取磁场激励信号,模数转换器30用于将模拟量的磁共振信号转换为数字量的磁共振信号,滤波及变参数控制单元20用于对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号以及根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,数字相关运算用于对滤波参考信号和滤波磁共振信号进行数字自相关运算和数字互相关运算以获取磁共振信号和数字参考信号之间的相位差,相位补偿单元60用于对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差,PID闭环控制单元70用于根据磁共振信号和磁场激励信号之间的相位差与期望相位差之间的大小调整数字参考信号和磁场激励信号的输出频率以完成磁共振信号相位闭环控制。
为了对本发明有进一步地了解,下面结合图1对本发明所提供的基于变参数控制的磁共振信号相位闭环控制方法进行详细说明。
如图1所示,根据本发明的具体实施例提供了一种基于变参数控制的磁共振信号相位闭环控制方法,该方法首先基于DDS技术产生一路磁场激励信号,一方面作为模拟输出,激励原子气室内的碱金属原子,一方面作为数字参考信号,用于提取磁共振信号的相位信息;然后搭建高速数据采集***,采集磁共振信号,结合参考信号,利用数字相关算法、变参数控制等方式,提取磁共振信号的相位信息;而后通过相位补偿,抑制***内的相位滞后及漂移;最后通过PID控制,改变DDS电路的输出频率,实现相位的闭环控制。该方法具体包括如下步骤。
步骤一,基于DDS技术生成数字参考信号和磁场激励信号。利用DDS内核单元10生成第一数字信号和第二数字信号,将第一数字信号作为数字参考信号,对第二数字信号进行数模转换以获取磁场激励信号。磁场激励信号基于DDS技术,由FPGA和DAC产生。其中,其中,DDS技术相对成熟,本发明不做论述。数字参考信号XA可根据
Figure BDA0003427702400000121
计算获取,其中,A为数字参考信号幅度,ω为数字参考信号频率,t为时间,θA为数字参考信号相位,
Figure BDA0003427702400000122
为数字参考信号噪声。
步骤二,根据磁场激励信号对原子气室进行磁场激励,获取原子气室的磁共振信号。磁共振信号XB与数字参考信号XA同频,存在相位差,磁共振信号XB可根据
Figure BDA0003427702400000123
计算获取,其中,B为磁场激励信号幅度,ω为磁场激励信号频率,t为时间,θB为磁场激励信号相位,
Figure BDA0003427702400000124
为磁场激励信号噪声。
步骤三,对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,根据数字参考信号和磁共振信号的信号频率对每个运算周期内的数字参考信号和磁共振信号的滤波参数进行调整,实现对数字参考信号和磁共振信号的变参数控制。
在本实施例中,为了降低数字参考信号XA和磁共振信号XB内的杂频噪声干扰,即可对数字参考信号和磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,数字参考信号和磁共振信号的滤波处理结果均可根据YN=a0XN-2+a1XN-1+a2XN+b1YN-1+b2YN-2来获取,其中,XN为当前时刻输入值,XN-1为前一时刻输入值,XN-2为前两时刻输入值,YN-1为前一时刻滤波结果,YN-2为前两时刻滤波结果,a0、a1、a2、b1和b2为滤波系数。由于数字参考信号XA和磁共振信号XB的信号频率ω均为已知量,为了降低杂频噪声干扰,每个运算周期内,滤波系数a0、a1、a2、b1和b2可进行适应性调整,保持信号的窄带处理,带宽可设置为1kHz,实现变参数控制。
在对数字参考信号和磁共振信号分别进行滤波处理后,滤波参考信号XAPF可表示为:XAPE=Asin(ωt+θAPE),滤波磁共振信号XBPF可表示为:XBPE=Bsin(ωt+θBPE),其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。
由于滤波磁共振信号和滤波参考信号同频,因此滤波后,两信号相位差与滤波前相同,即θAPFBPF=θAB,不影响相位提取。
步骤四,根据滤波参考信号和滤波磁共振信号计算获取磁共振信号和数字参考信号之间的相位差。在本实施例中,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差具体包括:根据滤波参考信号和滤波磁共振信号进行数字互相关运算和数字自相关运算以获取第一自相关运算量、第二自相关运算量和互相关运算量,根据第一自相关运算量、第二自相关运算量和互相关运算量计算获取磁共振信号和参考信号的相位差。
其中,第一自相关运算量XAA、第二自相关运算量XBB和互相关运算量XAB可根据
Figure BDA0003427702400000131
计算获取,其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。基于第一自相关运算量XAA、第二自相关运算量XBB和互相关运算量XAB,磁共振信号和数字参考信号之间的相位差Δθ可根据
Figure BDA0003427702400000132
计算获取。
步骤五,对磁共振信号和数字参考信号之间的相位差进行相位补偿以获取磁共振信号和磁场激励信号之间的相位差。在本实施例中,***需要求解的为磁共振信号与模拟激励信号之间的相位差,而参考信号与模拟激励信号之间存在相位滞后,且该滞后随温度漂移,因此需要对相位进行补偿。
经分析,***相位滞后主要由固定时间延迟τ带来,因此相位滞后
Figure BDA0003427702400000141
可表示为:
Figure BDA0003427702400000142
其中,ω为磁共振信号与数字参考信号的频率,为已知量。
相位漂移
Figure BDA0003427702400000143
与温度T有关,相位漂移
Figure BDA0003427702400000144
可表示为:
Figure BDA0003427702400000145
其中,温度漂移可通过多个温度点试验条件下分别获取与多个温度点对应的相位漂移,通过多个温度点以及多个相位漂移拟合曲线,根据拟合曲线计算获取相位漂移。
在获取了***相位滞后和相位漂移后,磁共振信号和磁场激励信号之间的相位差θ可根据
Figure BDA0003427702400000146
计算获取,其中,
Figure BDA0003427702400000147
为***相位滞后,
Figure BDA0003427702400000148
为相位漂移。
步骤六,将磁共振信号和磁场激励信号之间的相位差与期望相位差相比较,当磁共振信号和磁场激励信号之间的相位差与期望相位差不相等时,调整数字参考信号和磁场激励信号的频率,直至磁共振信号和磁场激励信号之间的相位差等于期望相位差,完成磁共振信号相位闭环控制。在本实施例中,将磁共振信号和磁场激励信号之间的相位差θ与期望相位差相比较,差值带入PID闭环控制算法中,得出激励信号频率的变化量,改变数字参考信号和磁场激励信号的频率,直至位差θ等于期望相位差,由此实现***的相位闭环控制。
由上可知,本实施例利用FPGA内的DDS内核和高速DAC芯片实现磁场激励信号的产生;数据采集***由高速ADC芯片实现,其余变参数控制、数字相关算法、相位补偿及PID控制等均在FPGA内部实现。
综上所述,本发明提供了一种基于变参数控制的磁共振信号相位闭环控制方法,该方法适合于碱金属原子磁强计磁共振信号相位闭环控制,该方法采用DDS技术、数字相关运算、变参数控制以及相位补偿等方法,代替传统的模拟相位解调控制方法,降低了***噪声干扰、温度漂移,提高了***的相位解算精度、稳定性及工程环境适应能力。
为了便于描述,在这里可以使用空间相对术语,如“在……之上”、“在……上方”、“在……上表面”、“上面的”等,用来描述如在图中所示的一个器件或特征与其他器件或特征的空间位置关系。应当理解的是,空间相对术语旨在包含除了器件在图中所描述的方位之外的在使用或操作中的不同方位。例如,如果附图中的器件被倒置,则描述为“在其他器件或构造上方”或“在其他器件或构造之上”的器件之后将被定位为“在其他器件或构造下方”或“在其他器件或构造之下”。因而,示例性术语“在……上方”可以包括“在……上方”和“在……下方”两种方位。该器件也可以其他不同方式定位(旋转90度或处于其他方位),并且对这里所使用的空间相对描述作出相应解释。
此外,需要说明的是,使用“第一”、“第二”等词语来限定零部件,仅仅是为了便于对相应零部件进行区别,如没有另行声明,上述词语并没有特殊含义,因此不能理解为对本发明保护范围的限制。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,所述基于变参数控制的磁共振信号相位闭环控制方法包括:
基于DDS技术生成数字参考信号和磁场激励信号;
根据所述磁场激励信号对原子气室进行磁场激励,获取所述原子气室的磁共振信号;
对所述数字参考信号和所述磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号,根据所述数字参考信号和所述磁共振信号的信号频率对每个运算周期内的所述数字参考信号和所述磁共振信号的滤波参数进行调整,实现对所述数字参考信号和所述磁共振信号的变参数控制;
根据所述滤波参考信号和所述滤波磁共振信号计算获取所述磁共振信号和所述数字参考信号之间的相位差;
对所述磁共振信号和所述数字参考信号之间的相位差进行相位补偿以获取所述磁共振信号和所述磁场激励信号之间的相位差;
将所述磁共振信号和所述磁场激励信号之间的相位差与期望相位差相比较,当所述磁共振信号和所述磁场激励信号之间的相位差与期望相位差不相等时,调整所述数字参考信号和所述磁场激励信号的频率,直至所述磁共振信号和所述磁场激励信号之间的相位差等于所述期望相位差,完成磁共振信号相位闭环控制。
2.根据权利要求1所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,对所述磁共振信号和所述数字参考信号之间的相位差进行相位补偿以获取所述磁共振信号和所述磁场激励信号之间的相位差具体包括:根据所述滤波参考信号和所述滤波磁共振信号进行数字互相关运算和数字自相关运算以获取第一自相关运算量、第二自相关运算量和互相关运算量,根据所述第一自相关运算量、所述第二自相关运算量和所述互相关运算量计算获取所述磁共振信号和所述参考信号之间的相位差。
3.根据权利要求2所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,对所述磁共振信号和所述数字参考信号之间的相位差进行相位补偿以获取所述磁共振信号和所述磁场激励信号之间的相位差具体包括:计算获取***相位滞后和相位漂移,根据所述***相位滞后和所述相位漂移对所述磁共振信号和所述数字参考信号之间的相位差进行相位补偿以获取所述磁共振信号和所述磁场激励信号之间的相位差。
4.根据权利要求1至3中任一项所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,基于DDS技术生成数字参考信号和磁场激励信号具体包括:利用DDS内核单元(10)生成第一数字信号和第二数字信号,将所述第一数字信号作为数字参考信号,对所述第二数字信号进行数模转换以获取所述磁场激励信号。
5.根据权利要求4所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,所述数字参考信号和所述磁共振信号的滤波处理结果均可根据YN=a0XN-2+a1XN-1+a2XN+b1YN-1+b2YN-2来获取,其中,XN为当前时刻输入值,XN-1为前一时刻输入值,XN-2为前两时刻输入值,YN-1为前一时刻滤波结果,YN-2为前两时刻滤波结果,a0、a1、a2、b1和b2为滤波系数。
6.根据权利要求2所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,所述第一自相关运算量XAA、所述第二自相关运算量XBB和所述互相关运算量XAB可根据
Figure FDA0003427702390000031
计算获取,其中,XAPF为滤波参考信号,XBPF为滤波磁共振信号,A为数字参考信号幅度,B为磁共振信号幅度,θAPF为滤波后的数字参考信号相位,θBPF为滤波后的磁共振信号相位。
7.根据权利要求1所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,所述磁共振信号和所述数字参考信号之间的相位差Δθ可根据
Figure FDA0003427702390000032
计算获取。
8.根据权利要求7所述的基于变参数控制的磁共振信号相位闭环控制方法,其特征在于,所述磁共振信号和所述磁场激励信号之间的相位差θ可根据
Figure FDA0003427702390000033
计算获取,其中,
Figure FDA0003427702390000034
为***相位滞后,
Figure FDA0003427702390000035
为相位漂移。
9.一种基于变参数控制的磁共振信号相位闭环控制***,其特征在于,所述磁共振信号相位闭环控制***使用如权利要求1至8中任一项所述的磁共振信号相位闭环控制方法进行磁共振信号相位闭环控制。
10.根据权利要求9所述的基于变参数控制的磁共振信号相位闭环控制***,其特征在于,所述磁共振信号相位闭环控制***包括DDS内核单元(10)、变参数控制单元(20)、模数转换器(30)、数模转换器(40)、数字相关运算单元(50)、相位补偿单元(60)和PID闭环控制单元(70),所述DDS内核单元(10)用于生成第一数字信号和第二数字信号,所述第一数字信号为数字参考信号,所述数模转换器(40)用于对所述第二数字信号进行数模转换以获取所述磁场激励信号,所述模数转换器(30)用于将模拟量的磁共振信号转换为数字量的磁共振信号,所述滤波及变参数控制单元(20)用于对所述数字参考信号和所述磁共振信号分别进行滤波处理以获取滤波参考信号和滤波磁共振信号以及根据所述数字参考信号和所述磁共振信号的信号频率对每个运算周期内的所述数字参考信号和所述磁共振信号的滤波参数进行调整,所述数字相关运算用于对所述滤波参考信号和所述滤波磁共振信号进行数字自相关运算和数字互相关运算以获取所述磁共振信号和所述数字参考信号之间的相位差,所述相位补偿单元(60)用于对所述磁共振信号和所述数字参考信号之间的相位差进行相位补偿以获取所述磁共振信号和所述磁场激励信号之间的相位差,所述PID闭环控制单元(70)用于根据所述磁共振信号和所述磁场激励信号之间的相位差与期望相位差之间的大小调整所述数字参考信号和所述磁场激励信号的输出频率以完成磁共振信号相位闭环控制。
CN202111585410.3A 2021-12-22 2021-12-22 基于变参数控制的磁共振信号相位闭环控制方法及*** Active CN114460506B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111585410.3A CN114460506B (zh) 2021-12-22 2021-12-22 基于变参数控制的磁共振信号相位闭环控制方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111585410.3A CN114460506B (zh) 2021-12-22 2021-12-22 基于变参数控制的磁共振信号相位闭环控制方法及***

Publications (2)

Publication Number Publication Date
CN114460506A true CN114460506A (zh) 2022-05-10
CN114460506B CN114460506B (zh) 2023-08-15

Family

ID=81405032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111585410.3A Active CN114460506B (zh) 2021-12-22 2021-12-22 基于变参数控制的磁共振信号相位闭环控制方法及***

Country Status (1)

Country Link
CN (1) CN114460506B (zh)

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6198353B1 (en) * 1999-08-05 2001-03-06 Lucent Technologies, Inc. Phase locked loop having direct digital synthesizer dividers and improved phase detector
CN1343393A (zh) * 1999-03-12 2002-04-03 西门子公司 发射接收机
US6950957B1 (en) * 2000-09-11 2005-09-27 Adc Telecommunications, Inc. Phase comparator for a phase locked loop
CN1866747A (zh) * 2006-06-23 2006-11-22 北京航空航天大学 电阻拾振式硅微机械谐振传感器的数字锁相闭环
US20110012618A1 (en) * 2009-07-15 2011-01-20 Maxim Integrated Products, Inc. Method and apparatus for sensing capacitance value and converting it into digital format
JP2013170976A (ja) * 2012-02-22 2013-09-02 Anritsu Corp 移動体通信機器試験用信号発生装置およびその周波数制御方法
US20160084925A1 (en) * 2014-09-19 2016-03-24 Commissariat A L'energie Atomique Et Aux Energies Alternatives Magnetometer without slaving and with compensation for fluctuations in the resonance gradient in weak field, magnetometers network and measurement method
US9397675B1 (en) * 2015-07-31 2016-07-19 Shure Acquisition Holdings, Inc. Hybrid frequency synthesizer and method
JP2018007097A (ja) * 2016-07-05 2018-01-11 三菱電機株式会社 シンセサイザ回路
CN109521373A (zh) * 2017-09-19 2019-03-26 北京自动化控制设备研究所 一种基于dds的原子磁强计相位闭环控制方法
CN109709496A (zh) * 2017-10-26 2019-05-03 北京自动化控制设备研究所 一种量子传感器闭环控制***及相位误差补偿控制方法
CN109752671A (zh) * 2017-11-03 2019-05-14 北京自动化控制设备研究所 一种原子磁强计光频移的稳定控制***
CN110646751A (zh) * 2019-09-18 2020-01-03 北京自动化控制设备研究所 基于同相激励的标量原子磁强计闭环控制***及方法
CN110673670A (zh) * 2019-09-12 2020-01-10 北京自动化控制设备研究所 原子气室温度闭环控制方法及装置
CN110927635A (zh) * 2019-12-23 2020-03-27 北京航空航天大学 一种基于极弱磁测量的小型磁强计磁闭环控制***及设计方法
CN110988760A (zh) * 2019-12-06 2020-04-10 天津大学 一种Mx型铯光泵磁力仪的数字化信号检测***
CN111060860A (zh) * 2018-10-17 2020-04-24 北京自动化控制设备研究所 一种自旋系综磁共振相位高带宽高精度检测方法
CN111337868A (zh) * 2020-03-23 2020-06-26 北京自动化控制设备研究所 一种核自旋纵弛豫时间测量方法
CN212134928U (zh) * 2019-12-06 2020-12-11 天津大学 一种Mx型铯光泵磁力仪的数字化信号检测***

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1343393A (zh) * 1999-03-12 2002-04-03 西门子公司 发射接收机
US6198353B1 (en) * 1999-08-05 2001-03-06 Lucent Technologies, Inc. Phase locked loop having direct digital synthesizer dividers and improved phase detector
US6950957B1 (en) * 2000-09-11 2005-09-27 Adc Telecommunications, Inc. Phase comparator for a phase locked loop
CN1866747A (zh) * 2006-06-23 2006-11-22 北京航空航天大学 电阻拾振式硅微机械谐振传感器的数字锁相闭环
US20110012618A1 (en) * 2009-07-15 2011-01-20 Maxim Integrated Products, Inc. Method and apparatus for sensing capacitance value and converting it into digital format
JP2013170976A (ja) * 2012-02-22 2013-09-02 Anritsu Corp 移動体通信機器試験用信号発生装置およびその周波数制御方法
US20160084925A1 (en) * 2014-09-19 2016-03-24 Commissariat A L'energie Atomique Et Aux Energies Alternatives Magnetometer without slaving and with compensation for fluctuations in the resonance gradient in weak field, magnetometers network and measurement method
US9397675B1 (en) * 2015-07-31 2016-07-19 Shure Acquisition Holdings, Inc. Hybrid frequency synthesizer and method
JP2018007097A (ja) * 2016-07-05 2018-01-11 三菱電機株式会社 シンセサイザ回路
CN109521373A (zh) * 2017-09-19 2019-03-26 北京自动化控制设备研究所 一种基于dds的原子磁强计相位闭环控制方法
CN109709496A (zh) * 2017-10-26 2019-05-03 北京自动化控制设备研究所 一种量子传感器闭环控制***及相位误差补偿控制方法
CN109752671A (zh) * 2017-11-03 2019-05-14 北京自动化控制设备研究所 一种原子磁强计光频移的稳定控制***
CN111060860A (zh) * 2018-10-17 2020-04-24 北京自动化控制设备研究所 一种自旋系综磁共振相位高带宽高精度检测方法
CN110673670A (zh) * 2019-09-12 2020-01-10 北京自动化控制设备研究所 原子气室温度闭环控制方法及装置
CN110646751A (zh) * 2019-09-18 2020-01-03 北京自动化控制设备研究所 基于同相激励的标量原子磁强计闭环控制***及方法
CN110988760A (zh) * 2019-12-06 2020-04-10 天津大学 一种Mx型铯光泵磁力仪的数字化信号检测***
CN212134928U (zh) * 2019-12-06 2020-12-11 天津大学 一种Mx型铯光泵磁力仪的数字化信号检测***
CN110927635A (zh) * 2019-12-23 2020-03-27 北京航空航天大学 一种基于极弱磁测量的小型磁强计磁闭环控制***及设计方法
CN111337868A (zh) * 2020-03-23 2020-06-26 北京自动化控制设备研究所 一种核自旋纵弛豫时间测量方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
FANG JIANCHENG; WAN SHUANGAI: "Light intensity fluctuation insensitive faraday polarimeter for atomic magnetometer", 《IEEE 2011 10TH INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT & INSTRUMENTS》, pages 215 - 218 *
JIE QIN,JIANCHENG FANG,HAIFENG DONG: "Control scheme for chip-scale atomic magnetometer", 《2009 9TH INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT & INSTRUMENTS》, pages 1 - 4 *
房建成,万双爱,陈瑶: "Measurement of 129Xe frequency shift due to Cs–129Xe collisions", 《CHIN. PHYS. B》, pages 1 - 7 *
楚中毅等人。: "无自旋交换弛豫原子磁强计的主动磁补偿", 《光学精密工程》, pages 1 - 6 *
董海峰,房建成,周斌权,秦杰,万双爱: "MEMS原子自旋陀螺气室芯片加工设备与工艺研究", 《仪器仪表学报》, pages 1 - 5 *
薛帅,汤恩琼,万双爱,李新纯: "基于同步双口SRAM 的原子磁强计磁共振相位提取***设计", 《导航定位与授时》, pages 76 - 80 *
马锦贵,孙晓光,庞喜浪,万双爱: "基于DDS+PLL的光源调制与功率稳定控制方法", 《导航定位与授时》, pages 164 - 167 *

Also Published As

Publication number Publication date
CN114460506B (zh) 2023-08-15

Similar Documents

Publication Publication Date Title
Berger et al. 200 and more NMR experiments: a practical course
CN109814163B (zh) 一种基于扩展补偿模型的航磁张量数据抑噪方法及***
CA2108036A1 (en) Method for correcting spectral and imaging data and for using such corrected data in magnet shimming
Liu et al. Research on a secondary tuning algorithm based on SVD & STFT for FID signal
CN109856690B (zh) 基于混合范数拟合的航磁梯度张量数据抑噪方法及***
Liu et al. Noise characterization for the FID signal from proton precession magnetometer
CN112924910A (zh) 一种基于原位磁强计的屏蔽桶内剩磁测量方法
Ge et al. Short-time and high-precision measurement method for Larmor frequency of marine overhauser sensor
CN110554333A (zh) 一种钾光泵磁力仪的光磁共振信号数字化检测***
CN110568384B (zh) 一种用于超灵敏原子磁力计的主动式磁补偿方法
KR101878570B1 (ko) 병렬 MR 데이터 캡처를 위한 수정된 TrueFISP 시퀀스
Tan et al. Nonlinear compensation method based on data for increasing absolute measurement precision of FID signal
Yan et al. MEMS accelerometer calibration denoising method for Hopkinson bar system based on LMD-SE-TFPF
CN109597137B (zh) 基于半导体磁传感器的Overhauser磁力仪跟踪配谐方法
Liu et al. Design and implementation of a tuning-matching framework for a high-sensitivity broad band proton precession magnetometer sensing coil
CN114460506A (zh) 基于变参数控制的磁共振信号相位闭环控制方法及***
Liu et al. High-precision sensor tuning of proton precession magnetometer by combining principal component analysis and singular value decomposition
CN112698255B (zh) 一种星载磁通门磁强计
CN112526415B (zh) 一种基于线性调频信号的磁屏蔽系数快速测量方法及装置
CN113341359B (zh) 一种Overhauser磁力仪磁测数据置信水平评价方法
EP0148362B1 (en) Method of obtaining pseudofiltering effect in process of accumulation and nuclear magnetic resonance spectrometry utilizing same
JPH0217477A (ja) スペクトルパラメータ決定方法及び装置
Tan et al. A frequency measurement method using rising-falling edge of square wave for increasing proton magnetometer precision
Wang et al. A fluxgate magnetometer for navigation and sensing: Noise character and digital filtering
Liu et al. Frequency estimation enhancement for industrial free induction decay signals under low SNR via Hankelization and modified covariance

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