CN117572309B - 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 - Google Patents
基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 Download PDFInfo
- Publication number
- CN117572309B CN117572309B CN202311370124.4A CN202311370124A CN117572309B CN 117572309 B CN117572309 B CN 117572309B CN 202311370124 A CN202311370124 A CN 202311370124A CN 117572309 B CN117572309 B CN 117572309B
- Authority
- CN
- China
- Prior art keywords
- signal
- coil
- harmonic
- response
- phase
- 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.)
- Active
Links
- 239000006249 magnetic particle Substances 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000011002 quantification Methods 0.000 title abstract description 7
- 230000005284 excitation Effects 0.000 claims abstract description 99
- 238000012545 processing Methods 0.000 claims abstract description 66
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 abstract description 15
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000001514 detection method Methods 0.000 description 17
- 239000002245 particle Substances 0.000 description 12
- 230000000694 effects Effects 0.000 description 11
- 238000002659 cell therapy Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000005415 magnetization Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- WTFXARWRTYJXII-UHFFFAOYSA-N iron(2+);iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[O-2].[Fe+2].[Fe+3].[Fe+3] WTFXARWRTYJXII-UHFFFAOYSA-N 0.000 description 2
- 229940031182 nanoparticles iron oxide Drugs 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
- G01R33/1276—Measuring magnetic properties of articles or specimens of solids or fluids of magnetic particles, e.g. imaging of magnetic nanoparticles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/1031—Investigating individual particles by measuring electrical or magnetic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electrochemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Engineering & Computer Science (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Nanotechnology (AREA)
- Dispersion Chemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
本发明公开了基于谐波相位偏转的磁粒子谱仪，包括电磁线圈模块、激励电路模块、接收电路模块和信号处理设备；电磁线圈模块包括激励线圈、偏置线圈、接收线圈和补偿线圈，待测样本放置于接收线圈中；激励电路模块，用于控制激励线圈和偏置线圈的工作；接收线圈，用于接收待测样本在激励磁场下的第一响应谐波信号以及在叠加磁场下的第二响应谐波信号并通过接收电路模块将其发送给信号处理设备，以使得信号处理设备根据两个响应谐波信号计算得出待测样本的SPIONs浓度。本发明能够大大提高SPOINs的浓度测量的精确度，提高MPS的灵敏度。本发明还公开了基于谐波相位偏转的磁粒子谱仪的定量方法及存储介质。
Description
技术领域
本发明涉及生物传感领域，尤其涉及一种基于谐波相位偏转的磁粒子谱仪和定量方法。
背景技术
磁粒子谱仪(全称为：Magnetic particle spectroscopy,MPS)是一种非侵入式的、高灵敏的超顺磁性氧化铁纳米颗粒(全称为：Superparamagnetic iron oxidenanoparticles，SPIONs)探测技术。其中，MPS的系统使用高幅值的正弦交变磁场来激励SPIONs信号，在外界交变磁场的作用下，SPIONs的磁矩发生翻转，并产生非线性响应谐波信号，通过响应谐波信号的探测与分析，MPS可实现SPIONs的浓度以及其周围环境(如：粘度、温度、SPIONs结合状态等)的探测。目前MPS的一个主要的应用为：探测SPIONs所标记的细胞。其中，提高MPS的探测灵敏度可让其探测到数量更少的标记细胞，从而有助于更加准确地评估细胞疗法的疗效。
由于传统MPS的系统的灵敏度受到SPIONs磁化性质、系统馈通信号、探测线圈灵敏度、系统噪音等因素的限制，导致其灵敏度很难进一步提高，无法探测到数量更少的标记细胞，进而无法更加准确地评估细胞疗法的疗效。
发明内容
为了克服现有技术的不足，本发明的目的之一在于提供基于谐波相位偏转的磁粒子谱仪，其能够解决现有技术中磁粒子谱仪测量SPIONs浓度灵敏度不高导致细胞疗法评估疗效不准确等问题。
本发明的目的之二在于提供基于谐波相位偏转的磁粒子谱仪的定量方法，其能够解决现有技术中磁粒子谱仪测量SPIONs浓度灵敏度不高导致细胞疗法评估疗效不准确等问题。
本发明的目的之三在于提供一种计算机可读存储介质，其能够解决现有技术中磁粒子谱仪测量SPIONs浓度灵敏度不高导致细胞疗法评估疗效不准确等问题。
本发明的目的之一采用如下技术方案实现：
基于谐波相位偏转的磁粒子谱仪，包括激励电路模块、电磁线圈模块、接收电路模块和信号处理设备；
其中，所述电磁线圈模块包括激励线圈、偏置线圈、接收线圈和补偿线圈，并且所述偏置线圈、激励线圈、接收线圈和补偿线圈为共轴的同心螺线管；所述接收线圈、补偿线圈内嵌入激励线圈中，所述接收线圈、补偿线圈上下相离设置，线圈绕线方向相反；
所述激励电路模块的第一输出端与所述激励线圈电性连接，用于生成第一电信号并发送给所述激励线圈以使得所述激励线圈产生激励磁场；所述激励电路模块的第二输出端与所述偏置线圈电性连接，用于生成第二电信号并发送给所述偏置线圈以使得所述偏置线圈产生偏置磁场；
所述接收线圈与所述接收电路模块电性连接，进而与所述信号处理设备电性连接，用于感应并接收所述激励线圈激发的待测样本的响应谐波信号并将所述响应谐波信号通过所述接收电路模块发送给所述信号处理设备；所述补偿线圈，用于抵消所述接收线圈所感应到的激励馈通信号；所述待测样本放置于所述接收线圈中；
所述信号处理设备，用于根据所述响应谐波信号得出所述响应谐波信号的相位和幅值，并根据待测样本的SPIONs浓度与响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度；其中，所述响应谐波信号包括第一响应谐波信号和第二响应谐波信号，并且所述第一响应谐波信号为激励磁场下待测样本激发时所产生的响应谐波信号、第二响应谐波信号为激励磁场与偏置磁场的叠加磁场下待测样本激发时所产生的响应谐波信号。
进一步地，所述接收电路模块，还用于对所述响应谐波信号进行过滤以得出N次谐波信号，并对N次谐波信号放大后发送到所述信号处理设备；其中，N为奇数。
进一步地，所述接收电路模块，包括带通电路模块和前置放大器；其中，所述带通电路模块的一端与所述接收线圈电性连接、另一端与所述前置放大器电性连接，用于将接收到的第一响应谐波信号进行过滤以得出对应的N次谐波信号，以及将接收到的第二响应谐波信号进行过滤以得出对应的N次谐波信号；
所述前置放大器，还与所述信号处理设备电性连接，用于对所述第一响应谐波信号的N次谐波信号进行放大后发送给所述信号处理设备，以及对所述第二响应谐波信号的N次谐波信号进行放大后发送给所述信号处理设备。
进一步地，所述激励电路模块，包括信号发生器、功率放大器和直流电源；
其中，所述信号发生器的一端与所述功率放大器电性连接、另一端与所述激励线圈电性连接，用于生成第一电信号并将所述第一电信号通过所述功率放大器放大后发送给所述激励线圈以使得所述激励线圈产生激励磁场；所述第一电信号为正弦信号；
所述直流电源与所述偏置线圈电性连接，用于向所述偏置线圈通入第二电信号以使得所述偏置线圈产生偏置磁场；所述第二电信号为直流电信号。
进一步地，所述偏置线圈为亥姆霍兹线圈；所述接收线圈与所述补偿线圈的尺寸相同。
进一步地，所述信号处理设备，包括数据采集卡、虚拟信号发生器和终端设备；其中，虚拟信号发生器设于终端设备内部；
其中，所述数据采集卡的一端与所述接收电路模块电性连接、另一端与所述终端设备电性连接，用于获取所述响应谐波信号并将所述响应谐波信号发送给所述终端设备；
所述终端设备，用于根据所述第一响应谐波信号和所述第二响应谐波信号得出参考信号的参数信息，并根据所述参考信号的参数信息控制所述虚拟信号发生器生成所述参考信号，进而根据所述第一响应谐波信号、第二响应谐波信号、参考信号以及公式(1)计算得出待测样本的SPIONs浓度；所述参考信号的参数信息包括参考信号的频率、相位和幅值；
其中，
式中：φ'1为第一叠加信号的相位，第一叠加信号为参考信号与第一响应谐波信号的叠加信号；φ'2为第二叠加信号的相位，第二叠加信号为参考信号与第二响应谐波信号的叠加信号；Ar为参考信号的幅值；φN为系统随机噪声引起的相位噪声，为常数；a为通过激励磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；b为通过叠加磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数，C为SPIONs浓度。
进一步地，在待测样本放入所述接收线圈之前：
首先，通过所述激励电路模块控制所述激励线圈产生激励磁场时，所述接收线圈，用于感应并接收所述激励线圈激发的第一固有谐波信号并发送给所述接收电路模块，从而使得所述接收电路模块将所述第一固有谐波信号发送给所述信号处理设备；所述信号处理设备，用于根据所述第一固有谐波信号计算得出第一基准参考信号的参数信息并存储于系统中；所述第一基准参考信号的参数信息包括相位和幅值，并且第一基准参考信号的相位为第一固有谐波信号的相位+π，第一基准参考信号的幅值为第一固有谐波信号的幅值；
然后，通过所述激励电路模块控制所述激励线圈产生激励磁场、控制所述偏置线圈产生偏置磁场时，所述接收线圈，用于感应并接收所述激励线圈与偏置线圈共同激发的第二固有谐波信号并发送给所述接收电路模块，从而使得所述接收电路模块将所述第二固有谐波信号发送给所述信号处理设备；所述信号处理设备，用于根据所述第二固有谐波信号计算得出第二基准参考信号的参数信息并存储于系统中；所述第二基准参考信号的参数信息包括相位和幅值，并且第二基准参考信号的相位为第二固有谐波信号的相位+π，第二基准参考信号的幅值为第二固有谐波信号的幅值；
当所述信号处理设备接收到第一响应谐波信号时，根据第一基准参考信号的参数信息生成第一基准参考信号，进而根据第一基准参考信号对第一响应谐波信号中的第一固有谐波信号的相位影响进行抵消；
当所述信号处理设备接收到第二响应谐波信号时，根据第二基准参考信号的参数信息生成第二基准参考信号，进而根据第二基准参考信号对第二响应谐波信号中的第二固有谐波信号的相位影响进行抵消。
本发明的目的之二采用如下技术方案实现：
基于谐波相位偏转的磁粒子谱仪的定量方法，应用于如本发明的目的之一采用的基于谐波相位偏转的磁粒子谱仪，包括：
单磁场激励步骤：在激励磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第一响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场；
叠加场激励步骤：在叠加磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第二响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场、同时产生第二电信号以控制所述偏置线圈产生偏置磁场，进而使得所述激励磁场与偏置磁场叠加形成叠加磁场；
计算步骤：根据第一响应谐波信号、第二响应谐波信号得出每个响应谐波信号的相位和幅值，并根据待测样本的SPIONs浓度与响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度。
进一步地，所述计算步骤还包括：
参考信号生成步骤：根据第一响应谐波信号的相位、幅值和第二响应谐波信号的相位、幅值计算得出参考信号的相位和幅值，并根据参考信号的相位和幅值生成参考信号；
叠加步骤：将参考信号与第一响应谐波信号叠加生成第一叠加信号、将参考信号与第二响应谐波信号叠加生成第二叠加信号；
计算步骤：根据第一叠加信号、第二叠加信号和公式(1)计算得出待测样本的SPIONs浓度；
其中，
式中：φ'1为第一叠加信号的相位；φ'2为第二叠加信号的相位；Ar为参考信号的幅值；φN为系统随机噪声引起的相位噪声，为常数；a为通过激励磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；b为通过叠加磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；C为SPIONs浓度。
本发明的目的之三采用如下技术方案实现：
一种计算机可读存储介质，其上存储有定量程序，所述定量程序为计算机程序，所述定量程序被处理器执行时实现如本发明的目的之二采用的基于谐波相位偏转的磁粒子谱仪的定量方法的步骤。
相比现有技术，本发明的有益效果在于：
本发明通过利用SPIONs谐波信号的相位偏转，将浓度变化转换为信号的相位变化，实现对SPIONs浓度的检测，尤其适合于低浓度范围的检测，可大大提高MPS的检测灵敏度和检测精准度，进而提高细胞疗法评估疗效的准确度。
附图说明
图1为本发明提供的基于谐波相位偏转的磁粒子谱仪的模块示意图；
图2为图1中的激励电路模块、接收电路模块、信号处理设备的模块示意图，以及与电磁线圈模块的连接示意图；
图3为图1中的电磁线圈模块中激励线圈、接收线圈、补偿线圈和偏置线圈的结构剖面图；
图4为本发明提供的基于谐波相位偏转的磁粒子谱仪的定量方法流程图；
图5为图4中的步骤S3的流程图；
图6为本发明提供的磁粒子谱仪校准的流程图；
图7为本发明提供的各个信号矢量变化图。
图中：121、激励线圈；122、偏置线圈；123、接收线圈；124、补偿线圈。
具体实施方式
下面，结合附图以及具体实施方式，对本发明做进一步描述，需要说明的是，在不相冲突的前提下，以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。
实施例一
本发明利用SPIONs的响应谐波信号的相位偏转的原理，通过探测SPIONs在未施加偏置磁场时产生的响应谐波信号的相位和施加偏置磁场时产生的响应谐波信号的相位实现对待测样本的SPIONs浓度的定量。
具体地，如图1-2所示，本发明提供一种优选的实施例，基于谐波相位偏转的磁粒子谱仪，包括：激励电路模块、电磁线圈模块、接收电路模块和信号处理设备。
其中，激励电路模块与电磁线圈模块电性连接，用于产生电信号并将电信号通入到电磁线圈模块的电磁线圈中，以使得电磁线圈产生磁场，以便激发放置于线圈中的待测样本。
电磁线圈模块，与接收电路模块电性连接，用于在电信号的作用下产生激励磁场和/或偏置磁场，以获取待测样本激发时产生的响应谐波信号并将响应谐波信号发送给接收电路模块。
具体地，电磁线圈模块，包括激励线圈、偏置线圈、补偿线圈和接收线圈。更为具体地，如图3所示，激励线圈121、偏置线圈122、接收线圈123和补偿线圈124为共轴的同心螺线管，并且偏置线圈122设于激励线圈121的上下两端；接收线圈123、补偿线圈124内嵌入激励线圈121中，接收线圈123、补偿线圈124上下相离设置，规格相同，线圈绕线方向相反。
激励线圈与激励电路模块的第一输出端电性连接，用于接收激励电路模块发送的第一电信号以产生激励磁场。在实际的测量过程中，将待测样本置入接收线圈的内部，激励线圈在第一电信号的作用下产生激励磁场，待测样本的SPIONs会在激励磁场得作用下产生响应谐波信号，并被接收线圈感应到。
偏置线圈，与激励电路模块的第二输出端电性连接，用于接收激励电路模块的第二电信号并在第二电信号的作用下产生偏置磁场，以便与激励线圈产生的激励磁场进行叠加以形成叠加磁场，实现叠加磁场对待测样本的激发。
接收线圈，与接收电路模块电性连接，用于感应并接收激励线圈激发的响应谐波信号并将响应谐波信号发送给接收电路模块。补偿线圈，用于对接收线圈感应到的激励线圈激发的响应谐波信号中的激励溃通信号进行补偿。
其中，响应谐波信号包括第一响应谐波信号和第二响应谐波信号。第一响应谐波信号为激励线圈产生激励磁场时的接收线圈感应到的待测样本的SPIONs激发时产生的响应谐波信号，第二响应谐波信号为激励线圈产生激励磁场与偏置线圈产生偏置磁场叠加后形成叠加磁场时的接收线圈感应到的待测样本的SPIONs激发时产生的响应谐波信号。在待测样本放置于接收线圈中时，首先控制激励线圈产生激励磁场以对待测样本进行激发，进而获取在未施加偏置磁场时的接收线圈感应到的待测样本的SPIONs激发时产生的响应谐波信号，也即第一响应谐波信号；然后再控制激励线圈产生激励磁场、并同时控制偏置线圈产生偏置磁场形成叠加磁场以对待测样本进行激发，进而获取在施加偏置磁场时(也即：叠加磁场)的接收线圈感应到的待测样本的SPIONs激发时产生的响应谐波信号，也即第二响应谐波信号。由于磁粒子的相位在未施加偏置磁场时，和施加偏置磁场时所产生的响应谐波信号的相位会发生偏转，因此，本发明通过将待测样本中的SPIONs浓度的变化转换为信号的相位变化，以利用相位偏转实现对SPIONs浓度的检测。另外，由于SPIONs的浓度发生微小的变化，响应谐波信号的相位的变化量也会较大，产生的响应谐波信号的相位的变化量相对于幅值的变化量更为明显，因此，本发明利用相位的变化量实现SPIONs浓度的测量，并且尤其适合于测量和检测低浓度的SPIONs，可有效提升MPS对低浓度SPIONs的信号检测能力，也即提高了系统探测灵敏度，更容易探测到数量更少的标记细胞，有助于更加准确地评估细胞疗法的疗效。
其中，偏置线圈优选为亥姆霍兹线圈。
接收电路模块，与信号处理设备电性连接，用于将接收到的响应谐波信号发送给信号处理设备。具体地，在激励磁场下时，接收第一响应谐波信号并发送给信号处理设备；在叠加磁场下时，接收第二响应谐波信号并发送给信号处理设备。
更为具体地，接收电路模块，包括带通电路模块和前置放大器。带通电路模块的一端与接收线圈电性连接、另一端通过前置放大器与信号处理设备电性连接。
其中，带通电路模块，用于接收响应谐波信号并对响应谐波信号进行过滤以得出N次谐波信号并将N次谐波信号发送给前置放大器进行放大。
前置放大器，与信号处理设备电性连接，用于将放大后的N次谐波信号发送给信号处理设备。同理，每个响应谐波信号均会经过上次相同的处理，以得出对应的N次谐波信号并发送给信号处理设备。
优选地，N为奇数。
由于响应谐波信号中的偶次谐波信号较为不明显，因此，本次采用奇次谐波信号作计算处理。也即，通过带通电路模块来保留响应谐波信号中的某个奇次谐波信号，比如三次谐波信号、五次谐波信号，并抑制其他的谐波信号。
更为优选地，接收电路模块，还可采用锁相环放大器来实现相同的功能。
进一步地，激励电路模块，包括信号发生器、功率放大器和直流电源。
其中，功率放大器的一端与信号发生器电性连接、另一端与电磁线圈模块的激励线圈电性连接。信号发生器，用于产生正弦信号并发送给功率放大器放大，并将放大后的正弦信号发送到电磁线圈模块的激励线圈，以使得激励线圈产生激励磁场。更为具体地，功率放大器的一端通过电线与信号发生器电性连接、另一端通过电线与电磁线圈模块的激励线圈电性连接。
直流电源，与电磁线圈模块的偏置线圈电性连接，用于产生直流电以将直流电通入电磁线圈模块的偏置线圈以产生偏置磁场。本发明将对激励线圈和偏置线圈的通电控制分开，可分别控制激励磁场和偏置磁场的产生。
进一步地，信号处理设备，用于对接收到的响应谐波信号进行数据分析处理以实现对待测样本的SPIONs浓度的检测。具体地，信号处理设备，用于根据第一响应谐波信号、第二响应谐波信号生成参考信号，并将第一响应谐波信号和参考信号叠加生成第一叠加信号、将第二响应谐波信号和参考信号叠加生成第二叠加信号，并结合待测样本的SPIONs浓度和叠加信号的关联关系计算得出待测样本的SPIONs浓度。其中，更为优选地，信号处理设备所处理的响应谐波信号均为接收电路模块处理后的N次谐波信号。
更为具体地，信号处理设备，包括数据采集卡、虚拟信号发生器和终端设备。其中，数据采集卡，与接收电路模块、终端设备电性连接，用于接收经过接收电路模块处理后的响应谐波信号，并将其转换为对应的数字信号后存储到终端设备中。
虚拟信号发生器，用于生成参考信号。一般来说，虚拟信号发生器集成于终端设备中。其中，终端设备优选为PC端(电脑)。具体地，终端设备在接收到第一响应谐波信号和第二响应谐波信号后，首先根据响应谐波信号计算得出第一响应谐波信号的相位和幅值、第二响应谐波信号的相位和幅值，然后根据第一响应谐波信号的相位和幅值、第二响应谐波信号的相位和幅值计算得出参考信号的参数信息，并控制虚拟信号发生器生成参考信号。
其中，参考信号的参数信息包括参考信号的相位和幅值。具体地，参考信号的相位为第一响应谐波信号的相位与第二响应谐波信号的相位的和的一半；参考信号的幅值为第一响应谐波信号与第二响应谐波信号中的幅值的最小值。另外，本发明还规定参考信号的频率与响应谐波信号的频率相同。
设定，第一响应谐波信号的相位为φ1、幅值为A1，第二响应谐波信号的相位为φ2，幅值为A2：
则参考信号的相位为φr为：φr＝(φ1+φ2)/2。
参考信号的幅值为Ar为：Ar＝min(A1,A2)。
终端设备通过上述计算关系可计算得出参考信号的相位和幅值，并控制虚拟信号发生器生成参考信号。
终端设备，还用于将生成的参考信号和第一响应谐波信号进行叠加以得出第一叠加信号，以及参考信号和第二响应谐波信号进行叠加以得出第二叠加信号，然后结合待测样本的SPIONs浓度和响应谐波信号相位的关联关系计算得出待测样本的SPIONs浓度。
更为具体地，终端设备，还利用公式(1)实现对待测样本的SPIONs浓度的测量。其中，
式中：φ'1为第一叠加信号的相位；φ'2为第二叠加信号的相位；Ar为参考信号的幅值；φN为系统随机噪声引起的相位噪声，为常数；a为通过激励磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；b为通过叠加磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；C为SPIONs浓度。其中，a，b为常数，是系统固有参数，通过测量标准浓度的SPIONs信号获得。
另外，本领域技术人员可知，待测样本的SPIONs(粒子)浓度C与感应到的响应信号的电压幅值相关或与感应到的响应信号的电压谐波幅值相关，具体可由以下公式推导得出：
首先、粒子的磁化强度M与粒子浓度成正比，具体由公式(2)表示：
其中，M为粒子的磁化强度，c为粒子的浓度，为粒子的平均磁矩。
其次、感应电压u(t)与粒子的磁化强度M成正比，具体由公式(3)表示：
其中，t为时间，u0为与系统感应信号灵敏度相关的常数；
再者、对感应电压做傅里叶分解，可得到谐波幅值uk与感应电压u(t)成正比，具体由公式(4)表示：
其中，k为谐波次数，fE为谐波频率；由于uk与u(t)成正比，结合公式(3)可知，谐波幅值uk与待测样本的粒子浓度C成正比。
响应谐波信号的幅值与相位差相关，具体为如公式(5)和公式(6)：
u'1(t)＝U1 sin(2πkfEt+φ1)+Arsin(2πkfEt+φr) (5)，
u'2(t)＝U2 sin(2πkfEt+φ2)+Arsin(2πkfEt+φr) (6)。
其中，U1、φ1分别为激励磁场激发下粒子谐波信号的幅值和相位，U2、φ2分别为叠加磁场激发下粒子谐波信号的幅值和相位，Ar、φr为加入的参考信号幅值和相位，通常(不考虑系统噪音干扰的情况下)u'1(t)和u'2(t)分别为第一叠加信号和第二叠加信号。
信号变化的矢量图如图7所示。一般来说，参考信号的相位设置为：因此，根据三角形法则可推导出：
又因为U1、U2与SPIONs粒子的浓度成正比，因此可令
U1＝b×c
U2＝a×c
其中a、b均为常数，是系统固有参数，通过测量标准浓度的SPIONs信号获得。
因此，公式(11)可改写为相位差与SPIONs粒子浓度的关系式，也即公式(12)：
其中，U1∠φ1为施加激励磁场时第k次谐波的矢量，U2∠φ2为施加叠加磁场时第k次谐波的矢量；Ar∠φr为引入的参考信号的矢量；U′1∠φ′1为第一叠加信号的矢量；U′2∠φ′2为第二叠加信号的矢量。
在系统运行过程中会因系统随机噪声而引起相位噪声，因此，将系统随机噪声引起的相位噪声φN考虑进入后，可对公式(12)进行修正后可得到公式(1)：
优选地，由于系统本身可能存在固有谐波，为了将系统固有谐波对相位的影响消除，则在待测样本放入到接收线圈之前，需要对磁粒子谱仪进行校准。
在待测样本未放入接收线圈时：首先通过激励电路模块控制激励线圈产生激励磁场，此时，接收线圈，感应并接收激励线圈激发的第一固有谐波信号并发送给接收电路模块，并经过接收电路模块处理后发送到信号处理设备。
信号处理设备，根据第一固有谐波信号得出第一固有谐波信号的相位和幅值，并根据第一固有谐波信号的相位和幅值得出第一基准参考信号的相位和幅值，并将第一基准参考信号的相位和幅值存储到系统中。其中，第一基准参考信号的相位为第一固有谐波信号的相位+π，第一基准参考信号的幅值为第一固有谐波信号的幅值。
然后，通过激励电路模块控制激励线圈产生激励磁场、同时控制偏置线圈产生偏置磁场，此时，接收线圈，感应并接收激励线圈激发的第二固有谐波信号并发送给接收电路模块，并经过接收电路模块处理后发送到信号处理设备。
信号处理设备，根据第二固有谐波信号得出第二固有谐波信号的相位和幅值，并根据第二固有谐波信号的相位和幅值得出第二基准参考信号的相位和幅值，并将第二基准参考信号的相位和幅值存储到系统中。第二基准参考信号的相位为第二固有谐波信号的相位+π，第二基准参考信号的幅值为第二固有谐波信号的幅值。
在待测样本放入接收线圈中后，信号处理设备在获取第一响应谐波信号时，还根据系统中存储的第一基准参考信号的相位和幅值生成第一基准参考信号，并根据第一基准参考信号对第一响应谐波信号进行处理，以消除第一响应谐波信号中的第一固有谐波信号的相位影响。
通过校准计算得出第一基准参考信号，以便通过第一基准参考信号实现对第一响应谐波信号中所存在的第一固有谐波信号的相位影响的抵消，避免第一固有谐波信号对后续计算的干扰。具体地，数据采集卡接收到第一响应谐波信号发送给终端设备后，终端设备根据第一基准参考信号的相位和幅值控制虚拟信号发生器生成第一基准参考信号，这样，终端设备可通过第一基准参考信号实现对第一响应谐波信号中所存在的第一固有谐波信号的相位影响的抵消。
同理，信号处理设备在获取第二响应谐波信号时，还根据系统中存储的第二基准参考信号的相位和幅值生成第二基准参考信号，并根据第二基准参考信号以消除第二响应谐波信号中存在的第二固有谐波信号的相位影响。
通过校准计算得出第二基准参考信号，以便通过第二基准参考信号实现对第二响应谐波信号中所存在的第二固有谐波信号的相位影响的抵消，避免第二固有谐波信号对后续计算的干扰。具体地，数据采集卡接收到第二响应谐波信号发送给终端设备后，终端设备根据第二基准参考信号的相位和幅值控制虚拟信号发生器生成第二基准参考信号，这样，终端设备可通过第二基准参考信号实现对第二响应谐波信号中所存在的第二固有谐波信号的相位影响的抵消。
本发明提供的基于谐波相位偏转的磁粒子谱仪系统和定量方法，通过SPIONs奇次谐波在施加偏置磁场前后的相位偏转角度进行SPIONs浓度的定量计算。在没有施加偏置磁场时，SPIONs响应谐波信号与传统MPS相同；当施加偏置磁场且偏置磁场幅值大于临界值时，SPIONs响应谐波信号的谐波相位发生180度的偏转。然而，单纯利用该180度的偏转无法完成定量计算，因此本发明引入了虚拟参考信号，利用与虚拟参考信号的叠加信号、以及SPIONs浓度与叠加信号相位的数学关系实现SPIONs的定量测量。
该方法将磁粒子的浓度变化转变为信号的相位变化，同时利用了磁粒子相位未偏转和偏转时的SPIONs响应谐波信号对磁粒子的浓度进行检测；在低浓度范围内，相位变化具有比幅值变化更大的变化量；这意味着在低浓度范围内，即使SPIONs浓度发生微小的变化，相位的变化量也会较大，相比之下，幅值变化相对较小。这种效应对于测量和检测低浓度的SPIONs非常有用，因为它能够提供更高的信号灵敏度和检测精度。本发明通过在相位域中的测量，可以更容易地捕捉到微弱的信号变化，并且可以更准确地分析SPIONs浓度的变化情况，有效提升了MPS对低浓度SPIONs的信号检测能力，即系统的灵敏度。并且本发明提供的信号补偿方法显著减小了系统固有的相位噪音，从而使得MPS的探测灵敏度得到了显著提升。
实施例二
基于实施例一，本发明还提供基于谐波相位偏转的磁粒子谱仪的定量方法，应用于基于谐波相位偏转的磁粒子谱仪，如图4所示，包括：
步骤S1、在激励磁场存在的情况下，通过信号处理设备获取接收线圈感应到的待测样本的第一响应谐波信号。
其中，待测样本放置于接收线圈内；激励电路模块产生第一电信号以控制激励线圈产生激励磁场。
步骤S2、在叠加磁场存在的情况下，通过信号处理设备获取接收线圈感应到的待测样本的第二响应谐波信号。
其中，激励电路模块产生第一电信号以控制激励线圈产生激励磁场、同时产生第二电信号以控制偏置线圈产生偏置磁场，进而使得激励磁场与偏置磁场叠加形成叠加磁场。
步骤S3、根据第一响应谐波信号、第二响应谐波信号生成参考信号，并根据参考信号、第一响应谐波信号、第二响应谐波信号以及待测样本的SPIONs浓度与第一响应谐波信号、第二响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度。
具体地，如图5所示，步骤S3还包括：
步骤S31、根据第一响应谐波信号的相位、幅值和第二响应谐波信号的相位、幅值计算得出参考信号的相位和幅值。
具体地，参考信号的相位为第一响应谐波信号的相位与第二响应谐波信号的相位和的一半。参考信号的幅值为第一响应谐波信号的幅值与第二响应谐波信号的幅值中的最小值。
步骤S32、根据参考信号的相位和幅值生成参考信号。具体地，通过控制虚拟信号发生器生成参考信号。
步骤S33、将参考信号与第一响应谐波信号叠加生成第一叠加信号、将参考信号与第二响应谐波信号叠加生成第二叠加信号。
同时在生成第一叠加信号、第二叠加信号后，还分别计算得出各个叠加信号的相位和幅值。
步骤S34、根据第一叠加信号、第二叠加信号和公式(1)计算得出待测样本的SPIONs浓度。
其中，
式中：φ'1为第一叠加信号的相位；φ'2为第二叠加信号的相位；Ar为参考信号的幅值；φN为系统随机噪声引起的相位噪声，为常数；a为通过激励磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；b为通过叠加磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；C为SPIONs浓度。
本发明通过推导得出待测样本的SPIONs浓度与叠加信号的相位的关联关系，以实现对待测样本的SPIONs浓度的计算。
优选地，步骤S1之前，在待测样本未放入接收线圈中时，如图6所示，还包括：
步骤S11、在激励磁场的作用下，获取接收线圈感应到的第一固有谐波信号并根据第一固有谐波信号的相位和幅值计算得出第一基准参考信号的相位和幅值。
在待测样本未放入到接收线圈时，激励线圈产生的激励磁场会激发固有谐波信号，本发明通过获取固有谐波信号来生成基准参考信号，以抵消固有谐波信号的相位影响，避免对后续的计算造成干扰。
也即，设定第一固有谐波信号的相位为φ0，幅值为A0，则：
第一基准参考信号的相位为φ0+π，幅值为A0。
这样，当第一固有谐波信号与第一基准参考信号叠加时，补偿了固有谐波信号的相位干扰。也即，在后续测量的过程中，通过产生第一基准参考信号以实现对第一响应谐波信号中的第一固有谐波信号的相位影响的抵消。也即，步骤S1中获取接收线圈感应到的待测样本的第一响应谐波信号时，还根据系统中存储的第一基准参考信号的相位和幅值生成第一基准参考信号，以对第一响应谐波信号中的第一固有谐波信号的相位影响进行抵消。
步骤S12、在叠加磁场的作用下，获取接收线圈感应到的第二固有谐波信号并根据第二固有谐波信号的相位和幅值计算得出第二基准参考信号的相位和幅值。
同理，设定第二固有谐波信号的相位为φ00，幅值为A00，则第二基准参考信号的相位为：φ00+π，幅值为A00。也即，当第二固有谐波信号与第二基准参考信号叠加时，补偿激励线圈所导致的固有谐波信号的相位影响。也即，在后续测量的过程中，通过产生第二基准参考信号以实现对第二响应谐波信号中的第二固有谐波信号的相位影响的抵消。
也即，步骤S2中获取接收线圈感应到的待测样本的第二响应谐波信号时，还根据系统中存储的第二基准参考信号的相位和幅值生成第二基准参考信号，以对第二响应谐波信号中的第二固有谐波信号的相位影响进行抵消。
本发明在待测样本定量之前，还对磁粒子谱仪进行校准，以消除系统自身的固有谐波信号对相位的影响，保证后续计算的准确度。
本发明利用磁粒子的相位偏转，将磁粒子的浓度变化转换为信号的相位变化，实现对磁粒子的浓度的自动检测；同时，在低浓度范围内，响应谐波信号的相位变化量具有比幅值变化量更大的特性，这意味着在低浓度范围内，即使SPIONs浓度发生微小的变化，相位的变化量也会较大，相比之下，幅值变化相对较小。本发明可大大提高MPS的检测灵敏度和检测精准度。
实施例三
基于实施例一，本发明还提供一种计算机可读存储介质，其上存储有定量程序，所述定量程序为计算机程序，所述定量程序被处理器执行时实现以下步骤：
单磁场激励步骤：在激励磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第一响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场；
叠加场激励步骤：在叠加磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第二响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场、同时产生第二电信号以控制所述偏置线圈产生偏置磁场，进而使得所述激励磁场与偏置磁场叠加形成叠加磁场；
计算步骤：根据第一响应谐波信号、第二响应谐波信号得出每个响应谐波信号的相位和幅值，并根据待测样本的SPIONs浓度与响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度。
上述实施方式仅为本发明的优选实施方式，不能以此来限定本发明保护的范围，本领域的技术人员在本发明的基础上所做的任何非实质性的变化及替换均属于本发明所要求保护的范围。
Claims (8)
1.基于谐波相位偏转的磁粒子谱仪，其特征在于，包括激励电路模块、电磁线圈模块、接收电路模块和信号处理设备；
其中，所述电磁线圈模块包括激励线圈、偏置线圈、接收线圈和补偿线圈，并且所述偏置线圈、激励线圈、接收线圈和补偿线圈为共轴的同心螺线管；所述接收线圈、补偿线圈内嵌入激励线圈中，所述接收线圈、补偿线圈上下相离设置，线圈绕线方向相反；
所述激励电路模块的第一输出端与所述激励线圈电性连接，用于生成第一电信号并发送给所述激励线圈以使得所述激励线圈产生激励磁场；所述激励电路模块的第二输出端与所述偏置线圈电性连接，用于生成第二电信号并发送给所述偏置线圈以使得所述偏置线圈产生偏置磁场；
所述接收线圈与所述接收电路模块电性连接，进而与所述信号处理设备电性连接，用于感应并接收所述激励线圈激发的待测样本的响应谐波信号并将所述响应谐波信号通过所述接收电路模块发送给所述信号处理设备；所述补偿线圈，用于抵消所述接收线圈所感应到的激励馈通信号；所述待测样本放置于所述接收线圈中；
所述信号处理设备，用于根据所述响应谐波信号得出所述响应谐波信号的相位和幅值，并根据待测样本的SPIONs浓度与响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度；其中，所述响应谐波信号包括第一响应谐波信号和第二响应谐波信号，并且所述第一响应谐波信号为激励磁场下待测样本激发时所产生的响应谐波信号、第二响应谐波信号为激励磁场与偏置磁场的叠加磁场下待测样本激发时所产生的响应谐波信号；
所述信号处理设备，包括数据采集卡、虚拟信号发生器和终端设备；其中，虚拟信号发生器设于终端设备内部；
其中，所述数据采集卡的一端与所述接收电路模块电性连接、另一端与所述终端设备电性连接，用于获取所述响应谐波信号并将所述响应谐波信号发送给所述终端设备；
所述终端设备，用于根据所述第一响应谐波信号和所述第二响应谐波信号得出参考信号的参数信息，并根据所述参考信号的参数信息控制所述虚拟信号发生器生成所述参考信号，进而根据所述第一响应谐波信号、第二响应谐波信号、参考信号以及公式（1）计算得出待测样本的SPIONs浓度；所述参考信号的参数信息包括参考信号的频率、相位和幅值；其中，（1）；式中：/>为第一叠加信号的相位，第一叠加信号为参考信号与第一响应谐波信号的叠加信号；/>为第二叠加信号的相位，第二叠加信号为参考信号与第二响应谐波信号的叠加信号；/>为参考信号的幅值；/>为系统随机噪声引起的相位噪声,为常数；/>为通过激励磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数；/>为通过叠加磁场激发的单位浓度SPIONs的响应谐波信号的幅值，为常数，C为SPIONs浓度。
2.根据权利要求1所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，所述接收电路模块，还用于对所述响应谐波信号进行过滤以得出N次谐波信号，并对N次谐波信号放大后发送到所述信号处理设备；其中，N为奇数。
3.根据权利要求2所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，所述接收电路模块，包括带通电路模块和前置放大器；其中，所述带通电路模块的一端与所述接收线圈电性连接、另一端与所述前置放大器电性连接，用于将接收到的第一响应谐波信号进行过滤以得出对应的N次谐波信号，以及将接收到的第二响应谐波信号进行过滤以得出对应的N次谐波信号；所述前置放大器，还与所述信号处理设备电性连接，用于对所述第一响应谐波信号的N次谐波信号进行放大后发送给所述信号处理设备，以及对所述第二响应谐波信号的N次谐波信号进行放大后发送给所述信号处理设备。
4.根据权利要求1所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，所述激励电路模块，包括信号发生器、功率放大器和直流电源；其中，所述信号发生器的一端与所述功率放大器电性连接、另一端与所述激励线圈电性连接，用于生成第一电信号并将所述第一电信号通过所述功率放大器放大后发送给所述激励线圈以使得所述激励线圈产生激励磁场；所述第一电信号为正弦信号；所述直流电源与所述偏置线圈电性连接，用于向所述偏置线圈通入第二电信号以使得所述偏置线圈产生偏置磁场；所述第二电信号为直流电信号。
5.根据权利要求1所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，所述偏置线圈为亥姆霍兹线圈；所述接收线圈与所述补偿线圈的尺寸相同。
6.根据权利要求1所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，在待测样本放入所述接收线圈之前：
首先，通过所述激励电路模块控制所述激励线圈产生激励磁场时，所述接收线圈，用于感应并接收所述激励线圈激发的第一固有谐波信号并发送给所述接收电路模块，从而使得所述接收电路模块将所述第一固有谐波信号发送给所述信号处理设备；所述信号处理设备，用于根据所述第一固有谐波信号计算得出第一基准参考信号的参数信息并存储于系统中；所述第一基准参考信号的参数信息包括相位和幅值，并且第一基准参考信号的相位为第一固有谐波信号的相位+π，第一基准参考信号的幅值为第一固有谐波信号的幅值；然后，通过所述激励电路模块控制所述激励线圈产生激励磁场、控制所述偏置线圈产生偏置磁场时，所述接收线圈，用于感应并接收所述激励线圈与偏置线圈共同激发的第二固有谐波信号并发送给所述接收电路模块，从而使得所述接收电路模块将所述第二固有谐波信号发送给所述信号处理设备；所述信号处理设备，用于根据所述第二固有谐波信号计算得出第二基准参考信号的参数信息并存储于系统中；所述第二基准参考信号的参数信息包括相位和幅值，并且第二基准参考信号的相位为第二固有谐波信号的相位+π，第二基准参考信号的幅值为第二固有谐波信号的幅值；当所述信号处理设备接收到第一响应谐波信号时，根据第一基准参考信号的参数信息生成第一基准参考信号，进而根据第一基准参考信号对第一响应谐波信号中的第一固有谐波信号的相位影响进行抵消；当所述信号处理设备接收到第二响应谐波信号时，根据第二基准参考信号的参数信息生成第二基准参考信号，进而根据第二基准参考信号对第二响应谐波信号中的第二固有谐波信号的相位影响进行抵消。
7.基于谐波相位偏转的磁粒子谱仪的定量方法，应用于权利要求1-6中任意一项所述的基于谐波相位偏转的磁粒子谱仪，其特征在于，包括：单磁场激励步骤：在激励磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第一响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场；叠加场激励步骤：在叠加磁场存在的情况下，通过所述信号处理设备获取所述接收线圈感应到的待测样本的第二响应谐波信号；其中，所述待测样本放置于所述接收线圈内；所述激励电路模块产生第一电信号以控制所述激励线圈产生激励磁场、同时产生第二电信号以控制所述偏置线圈产生偏置磁场，进而使得所述激励磁场与偏置磁场叠加形成叠加磁场；计算步骤：根据第一响应谐波信号、第二响应谐波信号得出每个响应谐波信号的相位和幅值，并根据待测样本的SPIONs浓度与响应谐波信号的相位的关联关系计算得出待测样本的SPIONs浓度。
8.一种计算机可读存储介质，其上存储有定量程序，其特征在于，所述定量程序为计算机程序，所述定量程序被处理器执行时实现如权利要求7所述的基于谐波相位偏转的磁粒子谱仪的定量方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311370124.4A CN117572309B (zh) | 2023-10-20 | 2023-10-20 | 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311370124.4A CN117572309B (zh) | 2023-10-20 | 2023-10-20 | 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117572309A CN117572309A (zh) | 2024-02-20 |
| CN117572309B true CN117572309B (zh) | 2024-04-09 |
Family
ID=89888849
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311370124.4A Active CN117572309B (zh) | 2023-10-20 | 2023-10-20 | 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117572309B (zh) |
Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103477566A (zh) * | 2011-04-13 | 2013-12-25 | 高通股份有限公司 | 用于电动车辆的无线充电的天线对准和车辆导引 |
| CN106537113A (zh) * | 2014-04-15 | 2017-03-22 | 斯派超科学股份有限公司 | 粒子计数器和分类系统 |
| CN106740221A (zh) * | 2017-01-06 | 2017-05-31 | 中国计量大学 | 一种基于低频pwm整流器的v2g无线电能双向传输装置 |
| EP3525013A1 (en) * | 2018-01-23 | 2019-08-14 | KONSTRUKTA - TireTech a.s. | Detector and localizer of ferromagnetic inclusions in a continually passing medium |
| CN110705072A (zh) * | 2019-09-24 | 2020-01-17 | 郑州轻工业学院 | 基于Fokker-Planck的高精度磁纳米温度估计方法与系统 |
| CN111183364A (zh) * | 2017-08-16 | 2020-05-19 | 加利福尼亚大学董事会 | 脉冲磁粒子成像系统和方法 |
| CN112212996A (zh) * | 2020-10-10 | 2021-01-12 | 郑州轻工业大学 | 一种高频激励磁场中磁纳米粒子测温的谐波幅值-温度方法 |
| CN113260412A (zh) * | 2018-05-15 | 2021-08-13 | 霍洛比姆技术有限公司 | 利用具有时间相关驻波干涉和相干强度放大的全息能量传送（het）进行能量的精确传送 |
| DE102020118102A1 (de) * | 2020-07-09 | 2022-01-13 | Universität Zu Lübeck | Verfahren zur Erfassung der magnetischen Antwort einer Verteilung magnetisierbarer Partikel auf ein magnetisches Anregungsfeld |
| CN114246574A (zh) * | 2021-12-30 | 2022-03-29 | 北京航空航天大学 | 基于谐波正交投影的磁粒子成像方法、系统、设备 |
| KR20220055897A (ko) * | 2020-10-27 | 2022-05-04 | 광주과학기술원 | 자기 입자 분광기 |
| CN114696643A (zh) * | 2022-03-15 | 2022-07-01 | 华南理工大学 | 一种基于n次谐波和相位同步控制的负电阻 |
| CN115133665A (zh) * | 2022-06-01 | 2022-09-30 | 华北电力大学 | 一种轻量化、抗短路的大功率无线充电电路 |
| CN115607112A (zh) * | 2022-11-29 | 2023-01-17 | 暨南大学附属第一医院（广州华侨医院） | 一种基于光磁声的一体化智能成像系统及方法 |
| CN115656898A (zh) * | 2022-11-02 | 2023-01-31 | 北京航空航天大学 | 一种多频磁场激励下磁纳米粒子磁化响应及频谱测量装置与方法 |
| CN115868929A (zh) * | 2022-11-29 | 2023-03-31 | 暨南大学附属第一医院（广州华侨医院） | 一种面向肿瘤乏氧在体定量研究的成像装置及方法 |
| CN116195988A (zh) * | 2023-02-13 | 2023-06-02 | 重庆邮电大学 | 一种基于直流偏置和二次谐波检测的磁粒子成像系统 |
| CN116503507A (zh) * | 2023-06-26 | 2023-07-28 | 中国科学院自动化研究所 | 基于预训练模型的磁粒子图像重建方法 |
| CN116509366A (zh) * | 2023-03-21 | 2023-08-01 | 西安电子科技大学 | 基于有源补偿的二维磁粒子成像信号接收链路及获取方法 |
| CN116520214A (zh) * | 2023-03-21 | 2023-08-01 | 西安电子科技大学 | 基于垂直磁化的快速单谐波磁粒子成像设备及方法 |
| WO2023161520A1 (en) * | 2022-02-28 | 2023-08-31 | Julius-Maximilians-Universität Würzburg | A method and system for measuring a modified property of a sample comprising magnetic particles in liquid suspension |
| CN116763357A (zh) * | 2022-06-20 | 2023-09-19 | 暨南大学附属第一医院（广州华侨医院） | 一种血管内超声探头的在体导航方法 |
| CN116893374A (zh) * | 2023-08-21 | 2023-10-17 | 北京航空航天大学 | 一种高灵敏磁纳米粒子谱仪及其波形信号转换方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8884617B2 (en) * | 2008-06-23 | 2014-11-11 | The Regents Of The University Of California | Magnetic particle imaging devices and methods |
| FR3026193B1 (fr) * | 2014-09-19 | 2016-12-23 | Commissariat Energie Atomique | Magnetometre sans asservissement et a compensation des fluctuations de la pente de resonance en champ faible, reseau de magnetometres et procede de mesure |
| WO2017062821A1 (en) * | 2015-10-08 | 2017-04-13 | The University Of Florida Research Foundation, Inc. | Magnetic nanoparticle spectrometer |
| US20220252546A1 (en) * | 2021-02-11 | 2022-08-11 | Government Of The United States Of America, As Represented By The Secretary Of Commerce | System and method for measurment of magnetic nanoparticles |
-
2023
- 2023-10-20 CN CN202311370124.4A patent/CN117572309B/zh active Active
Patent Citations (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103477566A (zh) * | 2011-04-13 | 2013-12-25 | 高通股份有限公司 | 用于电动车辆的无线充电的天线对准和车辆导引 |
| CN106537113A (zh) * | 2014-04-15 | 2017-03-22 | 斯派超科学股份有限公司 | 粒子计数器和分类系统 |
| CN106740221A (zh) * | 2017-01-06 | 2017-05-31 | 中国计量大学 | 一种基于低频pwm整流器的v2g无线电能双向传输装置 |
| CN111183364A (zh) * | 2017-08-16 | 2020-05-19 | 加利福尼亚大学董事会 | 脉冲磁粒子成像系统和方法 |
| CN116879818A (zh) * | 2017-08-16 | 2023-10-13 | 加利福尼亚大学董事会 | 脉冲磁粒子成像系统和方法 |
| EP3525013A1 (en) * | 2018-01-23 | 2019-08-14 | KONSTRUKTA - TireTech a.s. | Detector and localizer of ferromagnetic inclusions in a continually passing medium |
| CN113260412A (zh) * | 2018-05-15 | 2021-08-13 | 霍洛比姆技术有限公司 | 利用具有时间相关驻波干涉和相干强度放大的全息能量传送（het）进行能量的精确传送 |
| CN110705072A (zh) * | 2019-09-24 | 2020-01-17 | 郑州轻工业学院 | 基于Fokker-Planck的高精度磁纳米温度估计方法与系统 |
| DE102020118102A1 (de) * | 2020-07-09 | 2022-01-13 | Universität Zu Lübeck | Verfahren zur Erfassung der magnetischen Antwort einer Verteilung magnetisierbarer Partikel auf ein magnetisches Anregungsfeld |
| CN112212996A (zh) * | 2020-10-10 | 2021-01-12 | 郑州轻工业大学 | 一种高频激励磁场中磁纳米粒子测温的谐波幅值-温度方法 |
| KR20220055897A (ko) * | 2020-10-27 | 2022-05-04 | 광주과학기술원 | 자기 입자 분광기 |
| CN114246574A (zh) * | 2021-12-30 | 2022-03-29 | 北京航空航天大学 | 基于谐波正交投影的磁粒子成像方法、系统、设备 |
| WO2023161520A1 (en) * | 2022-02-28 | 2023-08-31 | Julius-Maximilians-Universität Würzburg | A method and system for measuring a modified property of a sample comprising magnetic particles in liquid suspension |
| CN114696643A (zh) * | 2022-03-15 | 2022-07-01 | 华南理工大学 | 一种基于n次谐波和相位同步控制的负电阻 |
| CN115133665A (zh) * | 2022-06-01 | 2022-09-30 | 华北电力大学 | 一种轻量化、抗短路的大功率无线充电电路 |
| CN116763357A (zh) * | 2022-06-20 | 2023-09-19 | 暨南大学附属第一医院（广州华侨医院） | 一种血管内超声探头的在体导航方法 |
| CN115656898A (zh) * | 2022-11-02 | 2023-01-31 | 北京航空航天大学 | 一种多频磁场激励下磁纳米粒子磁化响应及频谱测量装置与方法 |
| CN115868929A (zh) * | 2022-11-29 | 2023-03-31 | 暨南大学附属第一医院（广州华侨医院） | 一种面向肿瘤乏氧在体定量研究的成像装置及方法 |
| CN115607112A (zh) * | 2022-11-29 | 2023-01-17 | 暨南大学附属第一医院（广州华侨医院） | 一种基于光磁声的一体化智能成像系统及方法 |
| CN116195988A (zh) * | 2023-02-13 | 2023-06-02 | 重庆邮电大学 | 一种基于直流偏置和二次谐波检测的磁粒子成像系统 |
| CN116509366A (zh) * | 2023-03-21 | 2023-08-01 | 西安电子科技大学 | 基于有源补偿的二维磁粒子成像信号接收链路及获取方法 |
| CN116520214A (zh) * | 2023-03-21 | 2023-08-01 | 西安电子科技大学 | 基于垂直磁化的快速单谐波磁粒子成像设备及方法 |
| CN116503507A (zh) * | 2023-06-26 | 2023-07-28 | 中国科学院自动化研究所 | 基于预训练模型的磁粒子图像重建方法 |
| CN116893374A (zh) * | 2023-08-21 | 2023-10-17 | 北京航空航天大学 | 一种高灵敏磁纳米粒子谱仪及其波形信号转换方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于SQUID的磁纳米粒子磁特性测量研究;张朴;成鹏;刘文中;;计量与测试技术;20160830(08);全文 * |
| 基于磁纳米测温的弱磁信号检测系统;李思;孙毅;张国良;;湖北民族学院学报(自然科学版);20171212(04);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117572309A (zh) | 2024-02-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101909768B1 (ko) | 표면 특성 검사 장치 및 표면 특성 검사 방법 | |
| CN107315150B (zh) | 一种正交基模磁通门传感器 | |
| KR100834846B1 (ko) | 자기전기 감수율 측정 시스템 | |
| US20230400534A1 (en) | Sensor using a field gradient in a given volume | |
| JP2009525481A (ja) | 基準ユニットを備える磁気センサデバイス | |
| EP3059596B1 (en) | Electric current measurement apparatus and electric current measurement method | |
| US6586930B1 (en) | Material thickness measurement using magnetic information | |
| Wang et al. | Multi-frequency imaging with non-linear calibration of magnetoresistance sensors for surface and buried defects inspection | |
| CN117572309B (zh) | 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 | |
| US7019518B2 (en) | Non-contacting apparatus and method for measurement of ferromagnetic metal thickness | |
| CN112513623B (zh) | 用于检测物质响应的方法和系统 | |
| CN111505556A (zh) | 测量磁通门传感器中的探头噪声的方法 | |
| CN112945427A (zh) | 利用巴克豪森效应测量焊缝处二维应力的方法及检测仪器 | |
| JP4917812B2 (ja) | 鉄系構造物の劣化診断方法 | |
| JP2015133397A (ja) | マグネトインピーダンスを用いた磁気センサー素子 | |
| CA2184770C (en) | Method and apparatus for correcting drift in the response of analog receiver components in induction well logging instruments | |
| KR101001616B1 (ko) | 가역투자율 측정 장치 | |
| US11035912B2 (en) | No-switching AC magnetic hall-effect measurement method | |
| KR101151536B1 (ko) | 와전류 신호 복조 모듈 및 이를 포함하는 와전류 신호 검사 장치 | |
| JPH07190991A (ja) | 変態率測定方法及び装置 | |
| JPH05281063A (ja) | 鋼材の張力測定装置 | |
| Wolframm et al. | PCB Coil Enables In Situ Calibration of Magnetoelectric Sensor Systems | |
| He et al. | High sensitive magnetic sensor with amorphous wire | |
| CN114720918B (zh) | 一种测量交变矢量磁场的方法及装置 | |
| RU2737030C1 (ru) | Петлескоп для исследования тонких магнитных пленок |
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 |