CN114460504B - Online measurement and control system and method for line width of atomic magnetometer - Google Patents

Online measurement and control system and method for line width of atomic magnetometer Download PDF

Info

Publication number
CN114460504B
CN114460504B CN202111585407.1A CN202111585407A CN114460504B CN 114460504 B CN114460504 B CN 114460504B CN 202111585407 A CN202111585407 A CN 202111585407A CN 114460504 B CN114460504 B CN 114460504B
Authority
CN
China
Prior art keywords
signal
frequency
module
phase
measurement
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
Application number
CN202111585407.1A
Other languages
Chinese (zh)
Other versions
CN114460504A (en
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 CN202111585407.1A priority Critical patent/CN114460504B/en
Publication of CN114460504A publication Critical patent/CN114460504A/en
Application granted granted Critical
Publication of CN114460504B publication Critical patent/CN114460504B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • 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)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

The application provides an atomic magnetometer linewidth on-line measurement and control system and method, comprising a signal processing component and a measurement and control component, wherein an excitation control module is used for generating a sine voltage excitation signal, a signal preprocessing module is used for preprocessing an acquired nuclear magnetic resonance voltage signal to acquire the amplitude and phase of a resonance signal, and an instruction analysis module is used for analyzing an instruction output by the measurement and control component; the measurement and control assembly comprises a second data receiving and transmitting module, a data processing and storing module and a control module, wherein the data processing and storing module is used for calculating and obtaining the line width of the atomic magnetometer according to the amplitude of the processed nuclear magnetic resonance signal and the frequency of the current excitation signal sent by the first data receiving and transmitting module and storing the line width, and the control module is used for adjusting and controlling the amplitude and the sweep frequency range of the excitation signal in real time according to the amplitude and the phase of the resonance signal. By applying the technical scheme of the application, the technical problems of complex line width test system, lower test efficiency and low measurement accuracy in the prior art are solved.

Description

Online measurement and control system and method for line width of atomic magnetometer
Technical Field
The application relates to the technical field of atomic magnetometers, in particular to an online line width measurement and control system and method for an atomic magnetometer.
Background
The magnetic detection technology utilizes a plurality of magnetic sensors to detect the environment and the target magnetic field so as to obtain the abnormal magnetic field characteristics of the target and further identify the target. The atomic magnetometer detects a magnetic field to be detected by utilizing atomic spin precession, is suitable for geomagnetic field environments, has high detection precision, can measure the total field strength of the magnetic field, and is widely applied to the field of magnetic detection. The line width of the atomic magnetometer is the line width of the atomic nuclear magnetic resonance, and the sensitivity of the atomic magnetometer can be represented. However, existing test systems are complex and have low test efficiency.
Disclosure of Invention
The application provides an online line width measurement and control system and method for an atomic magnetometer, which can solve the technical problems of complex line width test system, lower test efficiency and low measurement precision in the prior art.
According to one aspect of the application, an atomic magnetometer linewidth online measurement and control system is provided, the atomic magnetometer linewidth online measurement and control system comprises a signal processing component and a measurement and control component, the signal processing component comprises an excitation control module, a signal preprocessing module, an instruction analysis module and a first data receiving and transmitting module, the excitation control module is used for generating sine voltage excitation signals with set amplitude and frequency, the signal preprocessing module is used for preprocessing acquired nuclear magnetic resonance voltage signals to acquire the amplitude and phase of resonance signals, the instruction analysis module is used for analyzing instructions output by the measurement and control component, and the first data receiving and transmitting module is used for realizing data transmission between the signal processing component and the measurement and control component; the measurement and control assembly comprises a second data receiving and transmitting module, a data processing and storing module and a control module, wherein the second data receiving and transmitting module is used for data transmission between the measurement and control assembly and the signal processing assembly, the data processing and storing module is used for calculating and obtaining the line width of the atomic magnetometer according to the amplitude of the processed nuclear magnetic resonance signal and the frequency of the current excitation signal sent by the first data receiving and transmitting module and storing, and the control module is used for adjusting and controlling the amplitude and the sweep frequency range of the excitation signal in real time according to the amplitude and the phase of the resonance signal.
Further, the measurement and control assembly further comprises an instruction module, wherein the instruction module is used for completing connection setting, magnetometer number selection, data protocol selection, light source parameter setting, excitation signal amplitude setting, magnetometer workflow and parameter setting and text file storage function setting.
Further, the measurement and control assembly further comprises a display module, wherein the display module is used for realizing real-time display of the data processing result of the data processing and storing module.
According to another aspect of the application, an atomic magnetometer linewidth online measurement and control method is provided, and the atomic magnetometer linewidth online measurement and control method uses the atomic magnetometer linewidth online measurement and control system to conduct linewidth online measurement and control.
Further, the online measurement and control method for the line width of the atomic magnetometer comprises the following steps: calculating rough frequency values of a first frequency point, a second frequency point and a third frequency point of the acquired reference signal; calculating and acquiring real parts and imaginary parts of the first frequency point, the second frequency point and the third frequency point based on rough frequency values of the first frequency point, the second frequency point and the third frequency point of the reference signal; calculating according to the real parts and the imaginary parts of the first frequency point, the second frequency point and the third frequency point to obtain the accurate frequency and the phase of the reference signal; the amplitude and phase of the resonance signal are obtained based on the exact frequency and phase calculation of the reference signal.
Further, the coarse frequency value of the first frequency point of the reference signal may be calculated according to fra_0_1=1000×n_zero×n/L, the coarse frequency value of the second frequency point of the reference signal may be calculated according to fra_0_2=1000×n (n_zero-1) ×n/L, the coarse frequency value of the third frequency point of the reference signal may be calculated according to fra_0_3 =1000×n_zero+1×n/L, where n_zero is a positive Zero crossing point, N is the number of sampling points of each window, and L is the signal length in each period.
Further, the signal length L in each period may be obtained by calculation from l= (l_int-l_int_lst) + (l_res_lst-l_res), where l_int=index_zero-index_zero_lst, l_res=a_pos/(a_pos-a_neg), l_int is an integer part of the signal length in any period, l_res is an integer part of the signal length in any period, l_int_lst is an integer part of the signal length in the last period of any period, l_res_lst is an integer part of the signal length in the last period of any period, index_zero is a sample point count value at the last Zero crossing point, index_zero_lst is a sample point count value at the last Zero crossing point of the last period, a_pos is a positive amplitude of the sample point signal at the last Zero crossing point, and a_neg is a negative amplitude of the sample point signal at the last Zero crossing point.
Further, the obtaining of the amplitude and phase of the resonance signal based on the accurate frequency and phase calculation of the reference signal specifically includes: generating standard sine and cosine signals according to the accurate frequency and the accurate phase of the reference signals; calculating and obtaining sine and cosine components based on standard sine and cosine signals; integrating the sine and cosine components and decomposing a difference formula to obtain a difference frequency component and a sum frequency component; performing low-pass filtering processing on the difference frequency component and the sum frequency component to reserve the difference frequency component; and calculating and acquiring the amplitude and the phase of the resonance signal according to the difference frequency component.
Further, the standard sine and cosine signal can be obtained by calculating sin_ref=sin (2×pi×freq_ref+phase_ref), cos_ref=cos (2×pi×freq_ref+phase_ref), wherein freq_ref is the precise frequency of the reference signal, and phase_ref is the Phase of the reference signal.
Further, the amplitude a_sig of the resonance signal may be obtained by calculating from a_sig=sqrt (r_sin+r_cos) and the phase_sig of the resonance signal may be obtained by calculating from phase_sig=arcsin (r_sin/a_sig), where r_sin and r_cos are difference frequency components.
By applying the technical scheme of the application, the online measurement and control system for the line width of the atomic magnetometer is provided, the functions of real-time resolving, displaying, sending control instructions and the like of the line width of the atomic magnetometer are completed through the measurement and control component, the functions of nuclear magnetic resonance signal acquisition pretreatment, sweep excitation signal generation output, communication with an upper computer measurement and control component and the like are completed through the signal processing circuit, and the line width test and control of the atomic magnetometer can be realized in real time under any environment.
Drawings
The accompanying drawings, which are included to provide a further understanding of embodiments of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the principles of the application. It is evident that the drawings in the following description are only some embodiments of the present application and that other drawings may be obtained from these drawings without inventive effort for a person of ordinary skill in the art.
FIG. 1 is a schematic diagram of an atomic magnetometer linewidth on-line measurement and control system provided in accordance with an embodiment of the present application;
FIG. 2 illustrates a schematic diagram of a measurement and control assembly provided in accordance with a specific embodiment of the present application;
fig. 3 shows a schematic diagram of a signal processing assembly provided in accordance with a specific embodiment of the present application.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. The following description of at least one exemplary embodiment is merely exemplary in nature and is in no way intended to limit the application, its application, or uses. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
The relative arrangement of the components and steps, numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present application unless it is specifically stated otherwise. Meanwhile, it should be understood that the sizes of the respective parts shown in the drawings are not drawn in actual scale for convenience of description. Techniques, methods, and apparatus known to one of ordinary skill in the relevant art may not be discussed in detail, but should be considered part of the specification where appropriate. In all examples shown and discussed herein, any specific values should be construed as merely illustrative, and not a limitation. Thus, other examples of the exemplary embodiments may have different values. It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further discussion thereof is necessary in subsequent figures.
As shown in fig. 1 to 3, according to a specific embodiment of the present application, an online measurement and control system for an atomic magnetometer linewidth is provided, where the online measurement and control system for an atomic magnetometer linewidth includes a signal processing module and a measurement and control module, the signal processing module includes an excitation control module, a signal preprocessing module, an instruction parsing module and a first data transceiver module, the excitation control module is used for generating sinusoidal voltage excitation signals with set amplitude and frequency, the signal preprocessing module is used for preprocessing acquired nuclear magnetic resonance voltage signals to obtain amplitude and phase of resonance signals, the instruction parsing module is used for parsing instructions output by the measurement and control module, and the first data transceiver module is used for implementing data transmission between the signal processing module and the measurement and control module; the measurement and control assembly comprises a second data receiving and transmitting module, a data processing and storing module and a control module, wherein the second data receiving and transmitting module is used for data transmission between the measurement and control assembly and the signal processing assembly, the data processing and storing module is used for calculating and obtaining the line width of the atomic magnetometer according to the amplitude of the processed nuclear magnetic resonance signal and the frequency of the current excitation signal sent by the first data receiving and transmitting module and storing, and the control module is used for adjusting and controlling the amplitude and the sweep frequency range of the excitation signal in real time according to the amplitude and the phase of the resonance signal.
By using the configuration mode, the system completes the functions of real-time calculation, display, control instruction sending and the like of the atomic magnetometer linewidth through the measurement and control component, completes the functions of nuclear magnetic resonance signal acquisition pretreatment, sweep excitation signal generation output, communication with an upper computer measurement and control component and the like through the signal processing circuit, and can realize real-time atomic magnetometer linewidth test and control under any environment.
Further, in the application, the measurement and control assembly further comprises an instruction module, wherein the instruction module is used for completing connection setting, magnetometer number selection, data protocol selection, light source parameter setting, excitation signal amplitude setting, magnetometer workflow and parameter setting and text file storage function setting. Specifically, the instruction module frames the front panel instruction buttons, data parameters and the like to form an instruction frame, and can complete connection setting, magnetometer number selection, data protocol selection, light source parameter setting, excitation signal amplitude setting, magnetometer workflow and parameter setting and text file storage function setting, wherein each function is realized through a parameter frame and buttons.
In addition, in the application, the measurement and control assembly also comprises a display module, wherein the display module is used for realizing the real-time display of the data processing result of the data processing and storage module, and comprises a numerical frame and graphic display. The control module can adjust and control the amplitude and the sweep frequency range of the excitation signal in real time according to the shape of the dispersion curve displayed in the display module, so that the measurement accuracy is maximized, and the proper sweep frequency range and the excitation signal amplitude are obtained. Specifically, when the dispersion curve is in a central symmetrical graph, the excitation signal amplitude and the sweep frequency range are reasonably set, and if the dispersion curve is asymmetrical, the excitation signal amplitude and the sweep frequency range are required to be adjusted; the amplitude of the excitation signal is adjusted to change the peak-to-valley difference value of the dispersion curve, and when the difference value is maximum, the line width measurement precision is highest, and the measurement result is more reliable.
For further understanding of the present application, the atomic magnetometer linewidth on-line measurement and control system provided by the present application is described in detail below with reference to fig. 1 to 3.
As shown in fig. 1 to 3, according to a specific embodiment of the present application, an online measurement and control system for line width of an atomic magnetometer is provided, where the online measurement and control system for line width of an atomic magnetometer includes a signal processing component and a measurement and control component, the measurement and control component includes a measurement and control terminal and an upper computer measurement and control software, and the measurement and control terminal runs the upper computer measurement and control software to complete functions of real-time calculation, display, control instruction sending, etc. of line width of the atomic magnetometer. The signal processing component comprises a signal processing circuit and signal processing software, and the signal processing circuit runs the signal processing software to complete the functions of nuclear magnetic resonance signal acquisition pretreatment, sweep frequency excitation signal generation and output, communication with the upper computer measurement and control software and the like.
The measurement and control terminal can be a desktop, a notebook or a tablet computer, an operating system is operated, and the upper computer measurement and control software can be realized by labview and comprises a flow control module, an instruction module, a data receiving and transmitting module, a data processing and storing module and a display module. The flow control module executes the modules in sequence, the data transceiver module interacts with the signal processing circuit, the instruction module frames a front panel instruction button, data parameters and the like to form an instruction frame, connection setting, magnetometer number selection, data protocol selection, light source parameter setting, excitation signal amplitude setting, magnetometer work flow and parameter setting and text file storage function setting can be completed, and each function is realized through a parameter frame and a button, wherein the connection setting comprises serial port selection and baud rate setting; the selection of the number of the magnetometers comprises the steps that a plurality of magnetometers can be connected, and the currently measured and controlled magnetometers are selected; data protocol selection includes selecting the frequency of the output transmission (20 Hz or 200 Hz), the content of the transmission (signal amplitude, measured magnetic field value, product temperature, product parameters, etc.); the light source parameter setting comprises temperature, current and PID control parameters of a laser of the atomic magnetometer can be changed; the excitation signal amplitude setting includes that the amplitude of the excitation signal can be altered, which can affect the measurement accuracy; the magnetometer workflow and parameter settings include: the magnetometer can be selected to work in an automatic mode or a manual mode, and whether to start functions such as magnetic field measurement, line width measurement and the like or not; the text file storage function setting includes selecting a storage folder, changing a file name, whether to restore, whether to stop storing, and the like.
The data processing and storing module is used for completing data processing and storing, completing a line width real-time calculating function, calculating line width, namely line width= |frequency value corresponding to the wave crest and wave trough and obtaining the difference between frequency points corresponding to the wave crest and wave trough of the nuclear magnetic resonance signal, simultaneously storing the data processing result in a text file in real time, transmitting the data processing result to the display module and displaying a dispersion curve in real time. The display module completes real-time display of the data processing result, including numerical frame and graphic display. The upper computer measurement and control software is a control module, and the measurement and control software can adjust and control the amplitude and the sweep frequency range of the excitation signal in real time according to the shape of the dispersion curve so as to maximize the measurement accuracy, and then the proper sweep frequency range and the excitation signal amplitude are obtained.
The signal processing circuit runs signal processing software, and the data receiving and transmitting module transmits a data preprocessing result to the upper computer at a fixed frequency (1 kHz) and can receive an instruction frame transmitted by the upper computer in real time. The signal processing circuit software is realized by c language and comprises a flow control module, a data receiving and transmitting module, an instruction analysis module, an excitation control module and a signal preprocessing module, wherein the signal processing circuit can transmit the amplitude value and the phase of a nuclear magnetic resonance signal and the frequency and the phase of a current excitation signal, and the upper computer draws a dispersion curve by taking the amplitude value as an ordinate and the frequency value as an abscissa. The flow control module executes the modules in sequence, and the signal processing circuit can complete corresponding functions according to the instructions; the data transceiver module interacts with the upper computer measurement and control software, the instruction analysis module finishes the analysis of the host control instruction, and the excitation control module generates sine voltage excitation signals with set amplitude and frequency and outputs the sine voltage excitation signals to the atomic magnetometer gauge head; the signal preprocessing module completes the processes of nuclear magnetic resonance signal phase locking calculation, filtering and the like. The signal preprocessing module can preprocess the acquired nuclear magnetic resonance voltage signals, and comprises front band-pass filtering, reference signal frequency phase calculation, nuclear magnetic resonance signal correlation calculation and rear low-pass filtering. The front band-pass filtering firstly filters the original resonance signals of the collected magnetometer and secondly filters the recovered excitation signals; only the signal in the sweep frequency range is reserved after filtering.
The measurement and control terminal and the signal processing circuit communicate by using an RS232 or RS422 serial port, wherein the baud rate of the RS232 serial port is 115200bps, and the baud rate of the RS422 serial port is 921600bps; the two carry out handshake according to a specified data transmission protocol; the measurement and control software can receive the data frame sent by the analysis signal processing circuit in real time, obtain the amplitude value of the nuclear magnetic resonance signal after pretreatment, and send a control instruction according to the test result.
According to another aspect of the application, an atomic magnetometer linewidth online measurement and control method is provided, and the atomic magnetometer linewidth online measurement and control method uses the atomic magnetometer linewidth online measurement and control system to conduct linewidth online measurement and control. The online line width measurement and control method of the atomic magnetometer comprises the following steps: calculating rough frequency values of a first frequency point, a second frequency point and a third frequency point of the acquired reference signal; calculating and acquiring real parts and imaginary parts of the first frequency point, the second frequency point and the third frequency point based on rough frequency values of the first frequency point, the second frequency point and the third frequency point of the reference signal; calculating according to the real parts and the imaginary parts of the first frequency point, the second frequency point and the third frequency point to obtain the accurate frequency and the phase of the reference signal; the amplitude and phase of the resonance signal are obtained based on the exact frequency and phase calculation of the reference signal.
By using the configuration mode, the method completes the functions of real-time calculation and display of the line width of the atomic magnetometer, sending a control instruction and the like through the measurement and control component, completes the functions of nuclear magnetic resonance signal acquisition pretreatment, sweep excitation signal generation and output, communication with an upper computer measurement and control component and the like through the signal processing circuit, and can realize real-time test and control of the line width of the atomic magnetometer under any environment. In addition, the line width calculation of the magnetometer is realized by the measurement and control method, and the accuracy and the calculation efficiency of the line width calculation can be effectively improved.
Specifically, the bandwidth of the front band-pass filter is the range of sweep-frequency excitation signals, which is set as [ f1, f2], the excitation control module generates sine voltage excitation signals with set amplitude and frequency, the excitation signals are used for enabling the atomic magnetometer to generate resonance, the reference signals are the recovered excitation signals, the real frequency and the Phase of the reference signals need to be calculated in real time, and the real frequency and the Phase of the reference signals are set as Freq_ref and phase_ref; the nmr signal correlation calculation is to calculate the amplitude and Phase of the original resonance signal based on the reference signal, and set the amplitude and Phase as a_sig and phase_sig, and set the original resonance signal as sin_sig=a_sig (2×pi×freq_ref+phase_sig).
And carrying out frequency phase calculation on the reference signal after band-pass filtering. As a specific embodiment of the present application, the calculation is performed with 1ms as a time window, and the sampling point of each window is set to N. Firstly, counting sampling data points, and setting the sampling data points as Index; the positive Zero crossing in the 1ms time window is calculated and set as N_zero, and the positive Zero crossing represents: the value of the original resonance signal at the previous moment is smaller than zero, and the value of the original resonance signal at the current moment is larger than or equal to zero. And recording the count value of the sampling point at the last Zero crossing point, setting the count value as index_zero, and setting the positive and negative amplitude values of the sampling point signals at the last Zero crossing point as A_Pos and A_Neg. The signal length within each cycle (1 ms per cycle in this embodiment) can be noted as:
L=(L_Int-L_Int_Lst)+(L_Res_Lst-L_Res)
L_Int=Index_Zero-Index_Zero_Lst
L_Res=A_Pos/(A_Pos-A_Neg)
wherein l_int is an integer part of the signal length in any period, l_res is a fractional part of the signal length in any period, l_int_lst is an integer part of the signal length in the last period of any period, l_res_lst is a fractional part of the signal length in the last period of any period, index_zero is a sampling point count value at the last Zero crossing point, index_zero_lst is a sampling point count value at the last Zero crossing point of the last period, a_pos is a positive amplitude of the sampling point signal at the last Zero crossing point, and a_neg is a negative amplitude of the sampling point signal at the last Zero crossing point.
The coarse frequency value of the first frequency point of the reference signal may be calculated and obtained according to fra_0_1=1000×n_zero×n/L, the coarse frequency value of the second frequency point of the reference signal may be calculated and obtained according to fra_0_2=1000×n/L (n_zero-1), the coarse frequency value of the third frequency point of the reference signal may be calculated and obtained according to fra_0_3 =1000×n_zero+1×n/L, where n_zero is a positive Zero crossing point, N is the number of sampling points of each window, and L is the signal length in each period.
According to the frequency values of the three frequency points, respectively calculating and obtaining the real part and the imaginary part of each frequency point, wherein the real part and the imaginary part are respectively:
dR1=sqrt(Re_1*Re_1+Im_1*Im_1);
dR2=sqrt(Re_2*Re_2+Im_2*Im_2);
dR3=sqrt(Re_3*Re_3+Im_3*Im_3);
Q=(dR2/dR3-2)/(dR2/dR3+1)+Fra_0_1。
wherein re_1 is the real part of the first frequency point, im_1 is the imaginary part of the first frequency point, re_2 is the real part of the second frequency point, im_2 is the imaginary part of the second frequency point, re_3 is the real part of the third frequency point, and im_3 is the imaginary part of the third frequency point.
Therefore, the accurate frequency for acquiring the reference signal can be calculated based on the frequency values of the three frequency points, and is as follows:
freq_ref=fs Q/N, where Fs is the sampling frequency and N is the sampling point in 1ms time.
The phases of the reference signals are: phase_ref=atan (im_1/re_1). The calculation is to obtain the accurate frequency and phase of the reference signal by performing FFT operation on three frequency points, wherein the operation on the three frequency points is used for estimation so as to improve the calculation accuracy of the frequency.
Further, after the accurate frequency and phase of the reference signal are obtained by calculation, the amplitude and phase of the resonance signal can be obtained based on the accurate frequency and phase of the reference signal. In the present application, the obtaining of the amplitude and phase of the resonance signal based on the accurate frequency and phase calculation of the reference signal specifically includes: generating standard sine and cosine signals according to the accurate frequency and the accurate phase of the reference signals; calculating and obtaining sine and cosine components based on standard sine and cosine signals; integrating the sine and cosine components and decomposing a difference formula to obtain a difference frequency component and a sum frequency component; performing low-pass filtering processing on the difference frequency component and the sum frequency component to reserve the difference frequency component; and calculating and acquiring the amplitude and the phase of the resonance signal according to the difference frequency component.
Specifically, a standard sine and cosine signal is generated according to the real frequency and the Phase freq_ref of the reference signal, wherein the standard sine and cosine signal can be obtained according to sin_ref=sin (2×pi×freq_ref+phase_ref), cos_ref=cos (2×pi×freq_ref+phase_ref), wherein freq_ref is the precise frequency of the reference signal, and phase_ref is the Phase of the reference signal.
Then, based on the standard sine and cosine signals, correlation operation is performed to obtain sine and cosine components, and the sine and cosine components are Rx_sin and Rx_cos.
Rx_sin=sin_ref*sin_sig=sin(2*pi*Freq_ref+Phase_ref)*A_sig*sin(2*pi*Freq_ref+Phase_sig)
Rx_cos=cos_ref*sin_sig=cos(2*pi*Freq_ref+Phase_ref)*A_sig*sin(2*pi*Freq_ref+Phase_sig)
After the integral sum-difference formula decomposition, the difference frequency component and the sum frequency component can be obtained, the difference frequency component is reserved by using low-pass filtering processing on the difference frequency component and the sum frequency component, and is recorded as R_sin and R_cos, and the amplitude and the phase of the original resonance signal are calculated based on the difference frequency component. The amplitude a_sig of the resonance signal may be obtained by calculating a_sig=sqrt (r_sin+r_cos×r_cos), and the Phase phase_sig of the resonance signal may be obtained by calculating phase_sig=arcsin (r_sin/a_sig), where r_sin and r_cos are difference frequency components.
In summary, the application provides an online measurement and control system for the line width of an atomic magnetometer, which completes the functions of real-time calculation, display, control instruction sending and the like of the line width of the atomic magnetometer through a measurement and control component, completes the functions of nuclear magnetic resonance signal acquisition pretreatment, sweep excitation signal generation output, communication with an upper computer measurement and control component and the like through a signal processing circuit, and can realize real-time measurement and control of the line width of the atomic magnetometer in any environment.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
In addition, the terms "first", "second", etc. are used to define the components, and are only for convenience of distinguishing the corresponding components, and the terms have no special meaning unless otherwise stated, and therefore should not be construed as limiting the scope of the present application.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims (8)

1. The online measurement and control method for the line width of the atomic magnetometer is characterized by using an online measurement and control system for the line width of the atomic magnetometer to conduct online measurement and control on the line width, and comprises the following steps:
calculating rough frequency values of a first frequency point, a second frequency point and a third frequency point of the acquired reference signal;
calculating and acquiring real parts and imaginary parts of a first frequency point, a second frequency point and a third frequency point of the reference signal based on rough frequency values of the first frequency point, the second frequency point and the third frequency point;
calculating and acquiring accurate frequency and phase of the reference signal according to real parts and imaginary parts of the first frequency point, the second frequency point and the third frequency point;
acquiring the amplitude and the phase of a resonance signal based on the accurate frequency and the phase of the reference signal;
the atomic magnetometer linewidth online measurement and control system comprises a signal processing component and a measurement and control component, wherein the signal processing component comprises an excitation control module, a signal preprocessing module, an instruction analysis module and a first data receiving and transmitting module, the excitation control module is used for generating sinusoidal voltage excitation signals with set amplitude and frequency, the signal preprocessing module is used for preprocessing acquired nuclear magnetic resonance voltage signals to acquire the amplitude and the phase of resonance signals, the instruction analysis module is used for analyzing instructions output by the measurement and control component, and the first data receiving and transmitting module is used for realizing data transmission between the signal processing component and the measurement and control component; the measurement and control assembly comprises a second data receiving and transmitting module, a data processing and storing module and a control module, wherein the second data receiving and transmitting module is used for data transmission between the measurement and control assembly and the signal processing assembly, the data processing and storing module is used for calculating and obtaining the line width of the atomic magnetometer according to the amplitude of the processed nuclear magnetic resonance signal and the frequency of the current excitation signal sent by the first data receiving and transmitting module and storing, and the control module is used for adjusting and controlling the amplitude and the frequency sweeping range of the excitation signal in real time according to the amplitude and the phase of the resonance signal.
2. The method for online measurement and control of line width of atomic magnetometer according to claim 1, wherein the coarse frequency value of the first frequency point of the reference signal is calculated and obtained according to fra_0_1=1000×n_zero×n/L, the coarse frequency value of the second frequency point of the reference signal is calculated and obtained according to fra_0_2=1000×n (n_zero-1) ×n/L, and the coarse frequency value of the third frequency point of the reference signal is calculated and obtained according to fra_0_3 =1000×n (n_zero+1) ×n/L, where n_zero is a positive Zero crossing point, N is the number of sampling points of each window, and L is the signal length in each period.
3. The method according to claim 2, wherein the signal length L in each period is calculated from l= (l_int-l_int_lst) + (l_res_lst-l_res), where l_int=index_zero-index_zero_lst, l_res=a_pos/(a_pos-a_neg), l_int is an integer part of the signal length in any period, l_res is an integer part of the signal length in any period, l_int_lst is an integer part of the signal length in the last period in any period, l_res_lst is a fractional part of the signal length in the last period in any period, index_zero is a count value of the sample point at the last Zero crossing point in the last period, a_pos is a count value of the sample point at the last Zero crossing point in the last period, and a_ps is a positive amplitude value of the sample point at the Zero crossing point in the last period.
4. The method for online measurement and control of line width of atomic magnetometer according to claim 3, wherein obtaining amplitude and phase of resonance signal based on accurate frequency and phase calculation of the reference signal specifically comprises:
generating standard sine and cosine signals according to the accurate frequency and the accurate phase of the reference signals;
calculating and obtaining sine and cosine components based on standard sine and cosine signals;
integrating the sine and cosine components and decomposing a difference formula to obtain a difference frequency component and a sum frequency component;
performing low-pass filtering processing on the difference frequency component and the sum frequency component to reserve the difference frequency component;
and calculating and acquiring the amplitude and the phase of the resonance signal according to the difference frequency component.
5. The method according to claim 4, wherein the standard sine-cosine signal is obtained by calculating sin_ref=sin (2×pi×freq_ref+phase_ref), cos_ref=cos (2×pi×freq_ref+phase_ref), wherein freq_ref is the precise frequency of the reference signal, and phase_ref is the Phase of the reference signal.
6. The method according to claim 5, wherein the amplitude a_sig of the resonance signal is obtained by calculating a_sig=sqrt (r_sin+r_cos) and the Phase phase_sig of the resonance signal is obtained by calculating phase_sig=arcsin (r_sin/a_sig), wherein r_sin and r_cos are difference frequency components.
7. The method of claim 1, wherein the measurement and control module further comprises an instruction module for implementing connection setup, magnetometer count selection, data protocol selection, light source parameter setup, excitation signal amplitude setup, magnetometer workflow and parameter setup, and text file storage function setup.
8. The online measurement and control method of the atomic magnetometer linewidth according to claim 7, wherein the measurement and control assembly further comprises a display module, and the display module is used for realizing real-time display of the data processing result of the data processing and storing module.
CN202111585407.1A 2021-12-22 2021-12-22 Online measurement and control system and method for line width of atomic magnetometer Active CN114460504B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111585407.1A CN114460504B (en) 2021-12-22 2021-12-22 Online measurement and control system and method for line width of atomic magnetometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111585407.1A CN114460504B (en) 2021-12-22 2021-12-22 Online measurement and control system and method for line width of atomic magnetometer

Publications (2)

Publication Number Publication Date
CN114460504A CN114460504A (en) 2022-05-10
CN114460504B true CN114460504B (en) 2023-08-15

Family

ID=81406014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111585407.1A Active CN114460504B (en) 2021-12-22 2021-12-22 Online measurement and control system and method for line width of atomic magnetometer

Country Status (1)

Country Link
CN (1) CN114460504B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104730484A (en) * 2015-03-23 2015-06-24 东南大学 Determination method for SERF of atomic spin magnetometer
CN105651649A (en) * 2016-01-27 2016-06-08 东南大学 Real-time online atomic density measuring method suitable for atom magnetometer
DE102016010324A1 (en) * 2015-08-20 2017-02-23 Leibniz-Institut für Photonische Technologien e. V. Arrangement for the high-resolution recording of magnetic field distributions or magnetic field changes
CN106597338A (en) * 2016-12-28 2017-04-26 北京航空航天大学 Method for measuring atomic transverse relaxation time based on electron resonance phase frequency analysis
US10901052B1 (en) * 2018-12-10 2021-01-26 AOSense, Inc. Atomic magnetometer

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104730484A (en) * 2015-03-23 2015-06-24 东南大学 Determination method for SERF of atomic spin magnetometer
DE102016010324A1 (en) * 2015-08-20 2017-02-23 Leibniz-Institut für Photonische Technologien e. V. Arrangement for the high-resolution recording of magnetic field distributions or magnetic field changes
CN105651649A (en) * 2016-01-27 2016-06-08 东南大学 Real-time online atomic density measuring method suitable for atom magnetometer
CN106597338A (en) * 2016-12-28 2017-04-26 北京航空航天大学 Method for measuring atomic transverse relaxation time based on electron resonance phase frequency analysis
US10901052B1 (en) * 2018-12-10 2021-01-26 AOSense, Inc. Atomic magnetometer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SERF原子磁强计最新进展及应用综述;王宇 等;光学仪器;全文 *

Also Published As

Publication number Publication date
CN114460504A (en) 2022-05-10

Similar Documents

Publication Publication Date Title
CN100462882C (en) Emulation test method of aerospace optical remote sensing machine focusing control circuit
CN106556811B (en) A kind of high-precision magnetic survey method for testing
CN101609136A (en) The method of voltage transformer secondary voltage drop tester and measurement phase angle difference thereof
CN104808165B (en) A kind of table automatic Verification method of electricity three
CN201935964U (en) Automatic intelligent detection device for quartz oscillator instrument
CN114460504B (en) Online measurement and control system and method for line width of atomic magnetometer
CN107367925A (en) A kind of multichannel rubidium clock automatic checkout system and method
CN216593446U (en) Closed loop fiber optic gyroscope bandwidth test system
CN101393240A (en) Electric network harmonic on-line supervisory circuit based on virtual instrument and working method thereof
CN205607423U (en) Automatic identification equipment of automobile instrument pointer indicated value
CN107089299B (en) A kind of measuring system and its measurement method of ship inclination angle
CN102721878A (en) Intelligentized automatic checkout device of quartz crystal oscillator instrument
CN115656910B (en) Remote calibration system, method and equipment for mutual inductor calibration instrument
CN109143389B (en) A kind of three-dimensional industrial frequency interference source for nuclear-magnetism quantifies orienting device and measurement method
CN100498229C (en) Method for processing periodic error in inertial components
CN109470887A (en) A kind of earthquake instrument system accelerometer Special testing device and method
CN110987000A (en) Method for accurately measuring zero-bias magnetic field sensitivity of laser gyroscope
CN114640825A (en) Intelligent terminal safety monitoring system
CN113001546B (en) Method and system for improving motion speed safety of industrial robot
CN210294400U (en) Transponder transmission module test equipment
CN102829941B (en) A kind of automatic mode of vibrational system calibrating
CN107991631B (en) Measuring method of magnetic signal measuring device independent of phase
CN202502024U (en) Virtual metering calibration instrument of fatigue testing machine
CN113188605B (en) Physical experiment method based on smart phone
CN102539266B (en) Metering and calibrating virtual instrument of fatigue testing machine and metering and calibrating method for fatigue test

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