阅读说明：本技术 面向高频交变磁场的金刚石nv色心磁力仪频率测量方法 (Diamond NV color center magnetometer frequency measurement method facing high-frequency alternating magnetic field ) 是由 郑培祥 罗钐 周涛 廖海清 周文俊 于 2021-06-04 设计创作，主要内容包括：本发明属于电子测量技术领域,具体涉及面向高频交变磁场的金刚石NV色心磁力仪频率测量方法。本发明针对脉冲操控速度的限制,采用将硬件上的限制转化为软件实现的方式,利用在相位积累为零时的时刻为交变磁场周期整数倍这一原理,将问题转化为一个求解有整数约束的优化模型；然后运用搜索算法对其中一个未知量进行赋值,有效解决未知量个数比方程个数多的问题,最终通过待测磁场的最小正周期得出待测交变磁场的频率。相比现有技术,本发明具有高普适性,可用于任意频率的交变磁场测量；耗时短,在几十毫秒的量级；可用于地磁场的精确测量,并在空间磁测、资源勘探、军事探潜以及地磁匹配导航都具有极重要的意义。(The invention belongs to the technical field of electronic measurement, and particularly relates to a frequency measurement method of a diamond NV color center magnetometer for a high-frequency alternating magnetic field. Aiming at the limitation of pulse control speed, the invention adopts a mode of converting the limitation on hardware into software realization and converts the problem into an optimization model with integral constraint for solving by utilizing the principle that the time when the phase accumulation is zero is integral multiple of the period of the alternating magnetic field; and then, assigning a value to one unknown quantity by using a search algorithm, effectively solving the problem of more unknown quantity numbers than the equation number, and finally obtaining the frequency of the alternating magnetic field to be detected through the minimum positive period of the magnetic field to be detected. Compared with the prior art, the method has high universality and can be used for measuring alternating magnetic fields with any frequency; the time consumption is short, and is in the order of tens of milliseconds; the geomagnetic field sensor can be used for accurately measuring the geomagnetic field and has extremely important significance in space magnetic measurement, resource exploration, military exploration and geomagnetic matching navigation.)

1. The method for measuring the frequency of the diamond NV color center magnetometer facing the high-frequency alternating magnetic field is characterized by comprising the following steps of:

step 1, measuring each time point t with zero phase superposition under the condition of not applying pi pulse₁，t₂...t_MThat is, the time from the initial time to the time point is integral multiple of the period of the alternating magnetic field to be measured, and each time point which is zero corresponds to the minimum positive period multiple of the magnetic field to be measured and is n₁，n₂...n_MSetting the period of the alternating magnetic field to be measured as T, wherein the three parameters comprise the following equation sets:

t₁＝n₁T，

t₂＝n₂T，

......

t_m＝n_mT，

......

t_M＝n_MT，

wherein M is the total quantity M of electricity at the moment when the measured phase accumulation is zero and is more than or equal to 7, and M is any one of 1-M;

step 2, let n_m1, bySolving the minimum positive period T of the alternating magnetic field;

step 3, solving n through the equation group in the step 1₁，n₂，...，n_M；

Step 4, solving n in step 3₁，n₂，...，n_MIf one is not an integer, let n_m＝n_m+1, return to step 2 until n₁～n_MEnding when the values are integers, wherein the corresponding T is the final minimum positive period T;

step 5, according to the T finally obtained in the step 4 and a formula f_mCalculating and outputting f as 1/T_mI.e. the frequency f of the alternating magnetic field to be measured finally_mAnd (6) obtaining the result.

Technical Field

The invention belongs to the technical field of electronic measurement, and particularly relates to a frequency measurement method of a diamond NV color center magnetometer for a high-frequency alternating magnetic field.

Background

With the progress of science and technology, the field of application of the alternating magnetic field is more and more extensive, and the development of various subject fields needs the measurement technology of the alternating magnetic field: such as biotechnology, industrial technology, communication technology, civil engineering, etc. The traditional alternating magnetic field measuring method comprises an electromagnetic induction method, a Hall effect method, a magnetoresistance effect method, a magnetic saturation method, a magnetic resonance method, a magneto-optical method and the like. However, the traditional magnetism measuring methods have some problems: the resolution of a magnetic measurement system is not high, most systems have strict requirements on the external environment, and the existence of noise causes great interference to the traditional method. In order to overcome the defects of the traditional magnetism measuring method, the diamond nitrogen-vacancy (NV) color center magnetism measuring method is gradually discovered by the world and rapidly developed, and becomes one of the popular magnetism measuring methods in recent years.

The existing NV color center magnetism measuring method also has some defects. The defect is embodied in the measurement of an alternating magnetic field, in the process of measuring the magnetic field by using a spin quantum interferometer with an NV color center, which is firstly proposed by Taylor in 2008, the measurement of the alternating magnetic field is realized by using a pulse sequence, because the accumulated phase of the alternating magnetic field in a period of time is small, an extra pi pulse is applied in the evolution process, so that the phase is overturned in a half period, the phase of the alternating field is continuously accumulated in the process, and when the accumulated phase reaches a certain value, the pulse speed is matched with the frequency of the alternating field, so that the frequency of the alternating magnetic field to be measured can be obtained. However, when the frequency of the alternating magnetic field to be measured increases, for example, to GHz, the method suffers from a bottleneck due to the limitation of the pulse control speed. Therefore, a technology which can control a high-frequency-band pulse sequence and well overcome the defects of the existing alternating magnetic field measurement technology is not available at present.

Disclosure of Invention

Aiming at the problems or the defects, the invention provides a frequency measurement method of a diamond NV color center magnetometer for a high-frequency alternating magnetic field, aiming at solving the problems of the existing diamond NV color center magnetism measurement method, and aiming at solving the defects by an integer constrained search magnetism measurement method.

A method for measuring the frequency of a diamond NV color center magnetometer facing a high-frequency alternating magnetic field is shown in figure 1 and comprises the following steps:

step 1, measuring each time point t with zero phase superposition under the condition of not applying pi pulse₁，t₂…t_M(namely time points of integral multiple of the period of the alternating magnetic field), the time from the initial time to the time points is integral multiple of the period of the alternating magnetic field to be measured, each time point which is zero corresponds to the minimum positive period multiple of the magnetic field to be measured and is n₁,n₂…n_M. Setting the period of the alternating magnetic field to be measured as T, wherein the three parameters comprise the following equation sets:

t₁＝n₁T,

t₂＝n₂T,

......

t_m＝n_mT,

......

t_M＝n_MT,

wherein M is the total quantity M of electricity at the moment when the measured phase accumulation is zero and is more than or equal to 7, and M is any one of 1-M.

Step 2, let n_m1, byThe minimum positive period T of the alternating magnetic field is determined.

Step 3, solving n through the equation group in the step 1₁,n₂,…,n_M。

Step 4, solving n in step 3₁,n₂,…,n_MIf one is not an integer, let n_m＝n_m+1, return to step 2 until n₁～n_MEnding when the values are all integers, wherein the corresponding T is the finally obtained minimum positive period T.

Step 5, according to the T finally obtained in the step 4 and a formula f_mCalculating and outputting f as 1/T_mI.e. is the mostFrequency f of alternating magnetic field to be measured_mAnd (6) obtaining the result.

Aiming at the limitation of the pulse control speed, the invention adopts a mode of converting the limitation on hardware into software realization and converts the problem into an optimization model with integral constraint for solving by utilizing the principle that the time when the phase accumulation is zero is integral multiple of the period of the alternating magnetic field. And then, assigning a value to one unknown quantity by using a search algorithm, effectively solving the problem of more unknown quantity numbers than the equation number, and finally obtaining the frequency of the alternating magnetic field to be detected through the minimum positive period of the magnetic field to be detected. Compared with the prior art, the method solves the problem that the method fails due to the limitation of the pulse control speed in the process of increasing the frequency of the alternating magnetic field to be detected, breaks through the bottleneck of the prior method, and has obviously better performance in all aspects than the prior art.

In the invention, step 1 is to list the equation set to be solved, and steps 2-4 are the cyclic part of the integer constraint problem by giving n_mAssign a value of n_mTake all positive integers from 1, using n_mSolving for T, solving for n_mA value of (M ═ 1.., M). When n is_mWhen M is an integer, i.e., e (0 is an error in the following model) is exactly 0 (minimum), the requirement is satisfied.

The specific optimization model of the invention is as follows:

min:e

n_m1,2

WhereinThe rounding-down operator. From the above model, it can be seen that the error e can be satisfied as small as possible (preferably zero) according to practical requirementsAll n_mAre all approximated to integers. Finding the minimum positive period T of the alternating magnetic field under the condition, and finally obtaining the minimum positive period T according to a formula f_mDetermining the frequency f of the alternating magnetic field to be measured 1/T_m。

In summary, the integer constraint-based search magnetism measurement method provided by the invention solves the problem that a pulse generator cannot control high-precision time on hardware, and can obtain an accurate measurement result through multiple measurements (namely, the measurement number M is greater than a certain value). Meanwhile, the invention has high universality and can be used for measuring alternating magnetic fields with any frequency. In addition, the method has short time consumption, is in the magnitude of tens of milliseconds, belongs to a simple and quick solution method, can be used for accurately measuring the geomagnetic field, and has extremely important significance in space magnetic measurement, resource exploration, military exploration and geomagnetic matching navigation.

Drawings

FIG. 1 is a flow chart of the present invention;

FIG. 2 is a graph of measurement accuracy as a function of M for the examples;

FIG. 3 is a graph of error e as a function of M for an example;

FIG. 4 is a graph of mean convergence time (in seconds) as a function of M in a simulation experiment of an embodiment;

FIG. 5 is a graph of the measurement accuracy of the example as a function of the frequency of the alternating magnetic field.

Detailed Description

The invention is described in further detail below with reference to the figures and examples.

Aiming at an alternating magnetic field, a common high-frequency magnetic sensor is an induction type magnetometer, the frequency response is above several hertz and below dozens or hundreds of kHz based on the Faraday electromagnetic induction principle, the alternating magnetic field above MHz is still difficult to respond, and the requirement of high-precision measurement of the geomagnetic field is difficult to meet. The method can improve the problems of the existing NV color center magnetism measuring method, essentially improves the measuring precision because the method only needs to use integer constraint for analysis, and can also effectively reduce the time of measuring the alternating magnetic field, reduce the complexity and solve some defects of the existing method.

The true bookIn the examples: the simulation adopts Monte Carlo experiment, firstly generates 10 randomly⁴T (i.e. the period of the alternating magnetic field to be measured) and then, on the basis of each T, randomly measuring M time points at which the phase is accumulated to 0. The measurement of T is carried out by utilizing the steps 1-5, and the measurement accuracy of the invention can be checked by judging whether the value of T randomly generated each time is equal to the value of T finally calculated.

Fig. 2 shows the change of the measurement accuracy of the present embodiment on T in the process of increasing M. It can be seen that the calculation accuracy is only about 60% when M is 2. As M increases, the measurement accuracy increases rapidly and reaches 99% when M is 7 and 100% after M > 7. Meanwhile, in the process of increasing M, the change situation of the size of the error e in the optimization model is also observed, as shown in FIG. 3. It can be seen from the figure that at M <5, the trend of e monotonically decreasing with M is very clear. After M >6 the error e approaches 0. The variation trend of the graph is basically consistent with that of FIG. 2, and the variation trend is monotonous when M <7, and approaches to the optimal value when M is larger than or equal to 7.

Fig. 4 shows the change in the mean convergence time during the increase of M. It can be seen that when M <7, as M increases, the average convergence time also increases. When M is more than or equal to 7 and fluctuates between 0.04 and 0.05 seconds, the average convergence time is considered to be stable. Through the simulation experiments, the measurement accuracy and the error e reach the optimal values when M is more than or equal to 7, and the average convergence time also reaches a stable value at the moment. Therefore, to ensure the accuracy of the algorithm, the value of M should be 7 or more, i.e. at least 7 points in time at which the phase accumulation is zero should be measured. At this time, the algorithm takes about 0.04 to 0.05 seconds, the time consumption is small, the algorithm cost is not large, and the method can be considered as rapid measurement.

Finally, we tested the situation under different T values, i.e. the situation where the frequency fm of the magnetic field to be measured is different, as shown in fig. 5. Several different frequency values were tested in simulation experiments, varying from 121.3MHz to 4123.9 MHz. As can be seen from the graph, when M is constant, the measurement accuracy rate is basically not changed along with fm, and a stable value can be maintained. When fm is constant but M is varied, the measurement accuracy is monotonically increasing with M, consistent with fig. 2. Therefore, the measurement accuracy of the method is related to M but is basically unrelated to fm, which shows that the method is insensitive to the frequency fm of the magnetic field to be measured, can be applied to the alternating magnetic field with any frequency, is not limited by the frequency, and has universality.