CN114440853A - Method for improving response speed of SERF (spin-relaxation free fiber) atomic spin gyroscope based on transient response calculation - Google Patents

Method for improving response speed of SERF (spin-relaxation free fiber) atomic spin gyroscope based on transient response calculation Download PDF

Info

Publication number
CN114440853A
CN114440853A CN202111580141.1A CN202111580141A CN114440853A CN 114440853 A CN114440853 A CN 114440853A CN 202111580141 A CN202111580141 A CN 202111580141A CN 114440853 A CN114440853 A CN 114440853A
Authority
CN
China
Prior art keywords
polarizability
electron
atomic spin
gyroscope
axis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111580141.1A
Other languages
Chinese (zh)
Other versions
CN114440853B (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN202111580141.1A priority Critical patent/CN114440853B/en
Publication of CN114440853A publication Critical patent/CN114440853A/en
Application granted granted Critical
Publication of CN114440853B publication Critical patent/CN114440853B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/58—Turn-sensitive devices without moving masses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Manufacturing & Machinery (AREA)
  • Gyroscopes (AREA)

Abstract

A method for improving the response speed of an SERF (spin-exchange reactor) atomic spin gyroscope based on transient response calculation comprises the steps of calibrating system parameters by utilizing the magnetic field response of the SERF atomic spin gyroscope, constructing a state space matrix of the SERF atomic spin gyroscope by utilizing the system parameters, collecting output signals when the SERF atomic spin gyroscope works, calculating the electron transverse polarization rate, calculating the nuclear transverse polarization rate by utilizing the electron transverse polarization rate, constructing a state vector by utilizing transient electrons and nuclear transverse polarization rate information, calculating the angular speed input, greatly improving the response speed of the SERF atomic spin gyroscope under the condition of rotation input, overcoming the problem of double-shaft coupling, and expanding the applicable range of the atomic gyroscope.

Description

Method for improving response speed of SERF (spin-relaxation free fiber) atomic spin gyroscope based on transient response calculation
Technical Field
The invention relates to an atomic gyroscope technology, in particular to a method for improving the response speed of an SERF (spin resonance filter) atomic spin gyroscope based on transient response calculation.
Background
Spin-Exchange Relaxation-Free (SERF) atomic Spin gyroscopes based on magnetic field, optical field and atomic interaction have ultra-high angular velocity measurement sensitivity21Theoretical sensitivity of atomic spin gyroscope with Ne atom as inert gas nucleus of 10-8°/s/Hz1/2Magnitude, 2.9 of maximum sensitivity over fiber optic gyroscopes-7°/s/Hz1/2. The SERF atomic spin gyroscope not only has ultrahigh angular velocity measurement sensitivity, but also has advantages in volume and cost compared with a gyroscope with the same precision, and is expected to be applied to a long-endurance and high-precision inertial navigation system, so that the SERF atomic spin gyroscope gradually becomes an important development direction of a novel inertial measurement technology. In addition, the SERF atomic spin inertial system with high sensitivity is widely applied to the leading edge research fields of charge space symmetry time symmetry breaking, abnormal interaction force detection, dark substance measurement and the like.
The development of the SERF atomic spin gyroscope has important significance in inertial navigation and front edge physics problem exploration. The atomic gyroscope indirectly measures the input magnitude of the angular velocity through measuring the steady state value of the transverse polarizability of electrons under the input of the angular velocity, and the response speed of the atomic gyroscope is determined by the Larmor precession frequency and the attenuation rate of inert gas nuclei. K-Rb + working at magnetic field compensation point21In the Ne inertial measurement unit, a first reference signal is generated,21the Larmor precession frequency of the Ne nucleus is about 0.366Hz in relation to the magnitude of the magnetic field it experiences, and the decay rate of its precession is about 1s-1. The angular velocity signal output needs to be subjected to a longer time (after the input of the fixed angular velocity is caused by the lower Larmor precession attenuation frequency)>5s) to reach a steady state valueTherefore, the output signal of the gyroscope is difficult to reach a steady state under the condition of complex or high-frequency rotation input, and the dynamic range of the input for measuring the angular speed is greatly limited. How to improve the dynamic performance and the response speed of the system on the premise of having measurement sensitivity becomes an urgent problem to be solved.
Disclosure of Invention
Aiming at the defects or shortcomings of the prior art, the invention provides a method for improving the response speed of an SERF (spin-nuclear reactor) atomic spin gyroscope based on transient response calculation, which comprises the steps of calibrating system parameters by utilizing the magnetic field response of the SERF atomic spin gyroscope, constructing a state space matrix of the SERF atomic spin gyroscope by utilizing the system parameters, collecting output signals when the SERF atomic spin gyroscope works, calculating the electron transverse polarizability, calculating the nuclear transverse polarizability by utilizing the electron transverse polarizability, constructing a state vector by utilizing transient electrons and nuclear transverse polarizability information, and calculating the angular speed input, so that the response speed under the rotation input of the SERF atomic spin gyroscope can be greatly improved, the problem of biaxial coupling is solved, and the applicable range of the atomic gyroscope is expanded.
The technical solution of the invention is as follows:
the method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation is characterized by comprising the following steps of:
step 1, calibrating a conversion coefficient;
step 2, calibrating the electronic relaxation rate and the optical frequency shift of the system;
step 3, testing and calibrating other parameters by the formants;
step 4, constructing a state space matrix;
step 5, resolving the electronic transverse polarizability in the working state;
step 6, resolving the transverse polarizability of the nuclei in the working state;
and 7, resolving system rotation input.
The step 1 comprises the following steps: heating an alkali metal gas chamber of an SERF atomic spin gyroscope to working temperature, compensating a magnetic field by adopting a magnetic field cross modulation compensation technology when atoms are polarized to a stable state by laser, and compensating the gyroscope by adopting the magnetic field cross modulation compensation technologyThe gyroscope works at a gyroscope compensation point; fixing a SERF atomic spin gyroscope in a static state of an inertial space, applying sine waves with amplitude of 0.3nT, frequency of 10mHz, 20mHz, 30mHz, 40mHz and 50mHz to an X axis and a Y axis respectively to obtain peak values of magnetic field response under different frequencies, drawing a frequency-peak value curve, fitting by adopting a straight line, wherein the slopes of the straight line are conversion coefficients K of X-axis and Y-axis signal output and electronic transverse polarizability respectivelyx、Ky。
The step 2 comprises the following steps: the SERF atomic spin gyroscope is fixed in a static state in an inertial space, a modulation square wave magnetic field with a peak-to-peak value of 0.3nT is input in the Y direction, and a magnetic field B in the Z direction is changedzThe magnitude of (2) is taken as an independent variable, and an output steady state solution difference value dS is measured to obtain an S-shaped curve, wherein the curve is as follows:
Figure BDA0003425804650000021
fitting to obtain the optical frequency shift L of the systemzTransverse electron relaxation rate
Figure BDA0003425804650000031
In which B iscBeing a point of self-compensation of the magnetic field, gammaeK is the gyromagnetic ratio of electrons and is a constant proportionality coefficient.
The step 3 comprises the following steps: fixing the static state of SERF atomic spin gyroscope in inertial space by changing different BzMeasuring the electron resonance peak according to the magnetic field size, and calculating the electron magnetic field BeNuclear magnetic field BnElectron polarizability
Figure BDA0003425804650000032
A slowdown factor Q, etc.
The step 4 comprises the following steps: and (3) after the basic parameters of the SERF atomic spin gyroscope are obtained through calibration in the steps 1 to 3, constructing a system matrix A and an input matrix B of the SERF atomic spin gyroscope.
The step 5 comprises the following steps: outputting X-axis signals S when SERF atomic spin gyroscope worksxY-axis signal output SyConversion coefficient K divided by X-axis and Y-axis signal output and transverse electric polarizability respectivelyx、KyObtaining the X-axis electron polarizability
Figure BDA0003425804650000033
And Y-axis electron polarizability
Figure BDA0003425804650000034
The step 6 comprises the following steps: the obtained X-axis electron polarizability
Figure BDA0003425804650000035
And Y-axis electron polarizability
Figure BDA0003425804650000036
Calculating the relative term of the nuclear spin transverse polarizability and the electron polarizability according to the system parameters
Figure BDA0003425804650000037
The step 7 comprises the following steps: the obtained X-axis electron nucleus polarizability
Figure BDA0003425804650000038
And polarizability of Y-axis electron nuclei
Figure BDA0003425804650000039
Calculating with the obtained system matrix A and input matrix B to obtain angular velocity input omegaxAnd omegay。
The invention has the following technical effects: the invention relates to a method for improving the response speed of an SERF (spin-relaxation free fiber) atomic spin gyroscope based on transient response calculation, which is characterized in that parameters such as the magnetic field amplitude-frequency response, the electron and nuclear relaxation rate, the electron and nuclear spin magnetic field and the electron spin progress slowing factor of a calibrated gyroscope are measured by utilizing the biaxial electron spin transverse output response of the SERF atomic gyroscope under the input of a magnetic field, and a system matrix and an input matrix are constructed. When the atomic gyroscope works, the biaxial electron spin transverse output transient response under the rotation input is collected, and the real-time rotation input can be obtained by resolving by using the calibration parameters. The method can greatly improve the response speed of the SERF atomic spin gyroscope under the condition of rotation input, overcomes the problem of biaxial coupling of the atomic gyroscope and enlarges the applicable range of the atomic gyroscope.
Drawings
FIG. 1 is a schematic flow chart of a method for calculating and increasing the response speed of a SERF atomic spin gyroscope based on transient response. Fig. 1 includes the following steps: step 1, calibrating a conversion coefficient, and measuring the conversion coefficient from the electronic transverse polarizability to signal output; step 2, testing an S curve (the S curve is an S-shaped curve), and calibrating the electronic relaxation rate and the optical frequency shift of the system; step 3, testing a resonance peak, and calibrating an electronic magnetic field, a nuclear magnetic field and a slow-down factor of the system; step 4, constructing a system matrix A and an input matrix B by using the system calibration parameters and the physical constants; step 5, collecting output signals when the gyroscope works, and calculating the transverse polarizability of electrons; step 6, resolving the nuclear transverse polarizability by utilizing the electronic transverse polarizability; and 7, constructing a state vector by using the transient electron and nuclear transverse polarizability information, and resolving the input of the angular velocity.
Fig. 2 is a schematic structural diagram of a system used for implementing the method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation. FIG. 2 includes a permalloy shielding cylinder, a ferrite, a magnetic field coil, an oven, and a gas chamber, a detection light path passing through the gas chamber along an X-axis and a pumping light path passing through the gas chamber along a Z-axis, a detection laser is disposed at a left end of the detection light path, a polarizer is disposed between the detection laser and the permalloy shielding cylinder, a Wollaston prism is disposed at a right end of the detection light path, 1/2 is disposed between the Wollaston prism and the permalloy shielding cylinder, the Wollaston prism is connected to a data collecting device through a differential amplifier, a polarization beam splitter is disposed on the pumping light path, 1/4 is disposed between the polarization beam splitter and a top of the permalloy shielding cylinder, a reflection side of the polarization beam splitter is connected to a control system through a photoelectric detector, and an incident side of the polarization beam splitter is connected to a liquid crystal light intensity stabilizer through a reflector, the liquid crystal light intensity stabilizer is respectively connected with the pumping laser and the control system.
Detailed Description
The invention is explained below with reference to the figures (fig. 1-2) and examples.
FIG. 1 is a schematic flow chart of a method for calculating and increasing the response speed of a SERF atomic spin gyroscope based on transient response. Fig. 2 is a schematic structural diagram of a system used for implementing the method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation. Referring to fig. 1 to 2, a method for calculating and improving the response speed of a SERF atomic spin gyroscope based on transient response is characterized by comprising the following steps: step 1, calibrating a conversion coefficient; step 2, calibrating the electronic relaxation rate and the optical frequency shift of the system; step 3, testing and calibrating other parameters by the formants; step 4, constructing a state space matrix; step 5, resolving the electronic transverse polarizability in the working state; step 6, resolving the nuclear transverse polarizability in the working state; and 7, resolving system rotation input.
The step 1 comprises the following steps: heating an alkali metal gas chamber of an SERF atomic spin gyroscope to a working temperature, compensating a magnetic field by adopting a magnetic field cross modulation compensation technology when the laser polarizes atoms to a stable state, and working the gyroscope at a gyroscope compensation point; fixing a SERF atomic spin gyroscope in a static state of an inertial space, applying sine waves with amplitude of 0.3nT, frequency of 10mHz, 20mHz, 30mHz, 40mHz and 50mHz to an X axis and a Y axis respectively to obtain peak values of magnetic field response under different frequencies, drawing a frequency-peak value curve, fitting by adopting a straight line, wherein the slopes of the straight line are conversion coefficients K of X-axis and Y-axis signal output and electronic transverse polarizability respectivelyx、Ky. The step 2 comprises the following steps: the SERF atomic spin gyroscope is fixed in a static state in an inertia space, a modulated square wave magnetic field with the peak value of 0.3nT is input in the Y direction, and the magnetic field B in the Z direction is changedzThe magnitude of (2) is taken as an independent variable, and an output steady state solution difference value dS is measured to obtain an S-shaped curve, wherein the curve is as follows:
Figure BDA0003425804650000051
fitting to obtain the optical frequency shift L of the systemzElectron relaxation rate of
Figure BDA0003425804650000052
In which B iscBeing a point of self-compensation of the magnetic field, gammaeK is the gyromagnetic ratio of electrons and is a constant proportionality coefficient.
The step 3 comprises the following steps: fixing the static state of SERF atomic spin gyroscope in inertial space by changing different BzMeasuring the electron resonance peak according to the magnetic field size, and calculating the electron magnetic field BeNuclear magnetic field BnElectron polarizability Pz eA slowdown factor Q, etc. The step 4 comprises the following steps: and (3) after the basic parameters of the SERF atomic spin gyroscope are obtained through calibration in the steps 1 to 3, constructing a system matrix A and an input matrix B of the SERF atomic spin gyroscope. The step 5 comprises the following steps: outputting X-axis signals S when SERF atomic spin gyroscope worksxY-axis signal output SyConversion coefficient K divided by X-axis and Y-axis signal output and transverse electric polarizability respectivelyx、KyObtaining the X-axis electron polarizability Px eAnd Y-axis electron polarizability Py e. The step 6 comprises the following steps: the obtained X-axis electron polarizability
Figure BDA0003425804650000053
And Y-axis electron polarizability
Figure BDA0003425804650000054
Calculating the relative term of the nuclear spin transverse polarizability and the electron polarizability according to the system parameters
Figure BDA0003425804650000055
. The step 7 comprises the following steps: the polarizability of the obtained X-axis electron nucleus
Figure BDA0003425804650000056
And polarizability of Y-axis electron nuclei
Figure BDA0003425804650000057
And calculating with the obtained system matrix A and input matrix B to obtain angular velocity input omegaxAnd omegay。
The invention provides a modeling and calibrating method of an SERF (spin-relaxation free) atomic spin gyroscope and a real-time transverse nuclear spin polarizability calculating method, overcomes the defect that the conventional SERF atomic spin gyroscope cannot nondestructively measure the nuclear polarizability in real time, measures the rotation input in real time by utilizing the transient response information of the atomic spin gyroscope, and greatly improves the response speed. In general, a SERF atomic spin gyroscope uses an electron transverse polarizability steady-state value as a rotation sensitive output, but the precession speed of an inert gas nucleus is slow, so that the electron transverse polarizability can reach a steady state after a long time. Such a measurement method has the following disadvantages: 1. the gyroscope of the measuring mode has slow response speed, and input cannot be accurately tracked for complex rotation input; 2. the measurement mode has the problem of sensitive shaft coupling, and the signal output change still exists for the rotation input change of a non-sensitive shaft. The above problems may cause that the rotation signal of the atomic spin gyroscope cannot be accurately measured in real time, so that the response speed of the gyroscope is slow and the bandwidth is narrowed. The present invention advantageously addresses the disadvantages of the above-described conventional approaches.
Fig. 2 illustrates the principle of a SERF atomic spin gyroscope and the xyz coordinate system setup.
Output signal S for an atomic spin gyroscopexAnd SyWe can convert the coefficient KxAnd KyConverting it into K-Rb atomic spin transverse polarizability
Figure BDA0003425804650000058
And
Figure BDA0003425804650000059
the kinetic equation for a SERF atomic spin gyroscope can be described approximately by the following Bloch equation:
Figure BDA0003425804650000061
where "x" represents the cross product of the vector, t is time,
Figure BDA0003425804650000062
are respectively K-Rb atom spin,21Spin polarizability of Ne nucleus, gammaeAnd gammanAre the spin-spin ratio of the electron spin and21ne nuclear spin gyromagnetic ratio, Q is an electron spin slowing factor. The pumping power and the photon polarization rate of the pumping light are respectively RpAnd Sp(ii) a The pumping rate and the photon polarization rate of the detection light are respectively RmAnd Sm. Electron spin and21the spin-exchange interaction between Ne nuclear spins polarizes21Ne atom having a spin exchange rate of
Figure BDA0003425804650000063
(ii) a Electron spin and21the spin-exchange interaction between Ne nuclear spins also induces electron spins with a spin-exchange rate of
Figure BDA0003425804650000064
. In the external magnetic field B ═ Bx,By,BzOptical frequency shift L ═ Lx,Ly,LzLambda M of inert gas equivalent magnetic fieldnPnAlkali metal electron equivalent magnetic field lambda MePeUnder the combined action of the rotation signal omega, the alkali metal electron spin and the inert gas nuclear spin generate precession. For describing relaxation process of longitudinal and transverse polarizability of atoms, longitudinal relaxation rate is introduced
Figure BDA0003425804650000065
And transverse relaxation rate
Figure BDA0003425804650000066
. When the angular velocity input and the magnetic field input are sufficiently small, the polarizability remains approximately constant, and there are:
Figure BDA0003425804650000067
the Bloch equation can be written linearized as:
Figure BDA0003425804650000068
the first two terms in the formula (3) are the change of the transverse electron spin polarizability along with time, the speed of the transverse electron spin polarizability to a steady state under the input of rotation or a magnetic field is less than 30ms, and the precession decay speed is about 1000 times faster than that of the inert gas nucleus Larmor. Relative to the slower precession of the inert gas nuclei Larmor, the precession of the electron spins can be approximately ignored and considered to be always in a steady state, i.e.:
Figure BDA0003425804650000069
combining vertical type (3) and formula (4) and neglecting pumping effect of nuclear spin on electron spin
Figure BDA0003425804650000071
The estimated value of the nuclear spin polarizability can be solved
Figure BDA0003425804650000072
And
Figure BDA0003425804650000073
Figure BDA0003425804650000074
in formula (5)
Figure BDA0003425804650000075
Which is a term related to the electron polarizability, can be calculated from the two lateral components of the electron polarizability,
Figure BDA0003425804650000076
the relevant items are entered for rotation. Spin electrons laterallyPolarizability of
Figure BDA0003425804650000077
Correlation term with electron polarizability in nuclear spin transverse polarizability
Figure BDA0003425804650000078
Defining state vectors as system state variables
Figure BDA0003425804650000079
Input vector u (t) ═ Ωx,Ωy,Bx,By]TThe state equation (3) of the linearized system may be organized in the form of a state space:
Figure BDA00034258046500000710
wherein the expression of the system matrix A is:
Figure BDA00034258046500000711
the expression of the input matrix B is:
Figure BDA00034258046500000712
and (3) carrying out matrix transformation on the formula (6), wherein the real-time change of the angular velocity input can be obtained by transient response calculation of the electron spin polarizability:
Figure BDA0003425804650000081
wherein B is-1Is the inverse of the input matrix B. Because of U (t) [. omega. ]x,Ωy,Bx,By]TTherefore, the first two items of U (t) can be taken to obtain the rotation information omegaxAnd omegay。
The method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation needs to calibrate the state of the SERF atomic spin gyroscope through 7 steps to obtain related parameters. And according to the obtained parameters, carrying out 7 steps of real-time calculation on angular speed input through data acquisition.
Step 1: conversion factor calibration
In the static state of an inertial space, the SERF atomic spin gyroscope is fixed, sine waves with the amplitude of 0.3nT and the frequencies of 10mHz, 20mHz, 30mHz, 40mHz and 50mHz are applied to the x axis and the y axis respectively, and peak values of magnetic field responses under different frequencies are obtained. Drawing a frequency-peak-to-peak value curve, fitting by adopting a straight line, wherein the slopes of the straight line are respectively the conversion coefficients K of the X-axis signal output, the Y-axis signal output and the electronic transverse polarizabilityx、Ky；
Step 2: electronic relaxation rate and optical frequency shift calibration of system
The SERF atomic spin gyroscope is fixed in a static state in an inertial space, a modulated square wave magnetic field with a peak-to-peak value of 0.3nT is input in the Y direction, and B is changedzThe magnitude of the magnetic field is used as an independent variable, and an S-shaped curve can be obtained by measuring the output steady state solution difference value dS. The following curves were used:
Figure BDA0003425804650000082
fitting to obtain the optical frequency shift L of the systemzElectron relaxation rate of
Figure BDA0003425804650000083
. Wherein B iscIs the magnetic field self-compensation point, and K is a constant proportionality coefficient.
And step 3: formant test calibration of other parameters
Fixing the static state of SERF atomic spin gyroscope in inertial space by changing different BzMeasuring the electron resonance peak according to the magnetic field size, and calculating the electron magnetic field BeNuclear magnetic field BnElectron polarizability Pz eA slowdown factor Q, etc.;
and 4, step 4: constructing a state space matrix
After basic parameters of the gyroscope are obtained through calibration in the three steps, a system matrix A and an input matrix B of the SERF atomic spin gyroscope are constructed according to the formulas (7) and (8);
and 5: resolution of operating state electron lateral polarizability
Outputting X-axis signals S when SERF atomic spin gyroscope worksxY-axis signal output SyConversion coefficient K divided by X-axis and Y-axis signal output and transverse electric polarizability respectivelyx、KyObtaining the X-axis electron polarizability
Figure BDA0003425804650000084
And Y-axis electron polarizability
Figure BDA0003425804650000091
Step 6: resolution of operating state nuclear transverse polarizability
The obtained X-axis electron polarizability
Figure BDA0003425804650000092
And Y-axis electron polarizability
Figure BDA0003425804650000093
Correlation term with system parameters to calculate electron polarizability in nuclear spin transverse polarizability in equation (5)
Figure BDA0003425804650000094
And 7: resolution of system rotational input
The obtained X-axis electron and nuclear polarizability
Figure BDA0003425804650000095
And Y-axis electron and nuclear polarizability
Figure BDA0003425804650000096
Substituting the obtained system matrix A and input matrix B into an equation (9) to calculate and obtain an angular velocity input omegaxAnd omegay。
Those skilled in the art will appreciate that the invention may be practiced without these specific details. It is pointed out here that the above description is helpful for the person skilled in the art to understand the invention, but does not limit the scope of protection of the invention. Any such equivalents, modifications and/or omissions as may be made without departing from the spirit and scope of the invention may be resorted to.

Claims (8)

1. The method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation is characterized by comprising the following steps of:
step 1, calibrating a conversion coefficient;
step 2, calibrating the electronic relaxation rate and the optical frequency shift of the system;
step 3, testing and calibrating other parameters by the formants;
step 4, constructing a state space matrix;
step 5, resolving the electronic transverse polarizability in the working state;
step 6, resolving the nuclear transverse polarizability in the working state;
and 7, resolving system rotation input.
2. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as claimed in claim 1, wherein the step 1 comprises: heating an alkali metal gas chamber of an SERF atomic spin gyroscope to a working temperature, compensating a magnetic field by adopting a magnetic field cross modulation compensation technology when the laser polarizes atoms to a stable state, and working the gyroscope at a gyroscope compensation point; fixing a SERF atomic spin gyroscope in a static state of an inertial space, applying sine waves with amplitude of 0.3nT, frequency of 10mHz, 20mHz, 30mHz, 40mHz and 50mHz to an X axis and a Y axis respectively to obtain peak values of magnetic field response under different frequencies, drawing a frequency-peak value curve, fitting by adopting a straight line, wherein the slopes of the straight line are conversion coefficients K of X-axis and Y-axis signal output and electronic transverse polarizability respectivelyx、Ky。
3. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as claimed in claim 1, wherein the step 2 comprises: the SERF atomic spin gyroscope is fixed in a static state in an inertial space, a modulation square wave magnetic field with a peak-to-peak value of 0.3nT is input in the Y direction, and a magnetic field B in the Z direction is changedzThe magnitude of (2) is taken as an independent variable, and an output steady state solution difference value dS is measured to obtain an S-shaped curve, wherein the curve is as follows:
Figure FDA0003425804640000011
fitting to obtain the optical frequency shift L of the systemzElectron relaxation rate of
Figure FDA0003425804640000012
Wherein B iscBeing a point of self-compensation of the magnetic field, gammaeIs the electron gyromagnetic ratio, and K is a constant proportionality coefficient.
4. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as claimed in claim 1, wherein the step 3 comprises: fixing the static state of the SERF atomic spin gyroscope in the inertial space by changing different BzMeasuring the electron resonance peak according to the magnetic field size, and calculating the electron magnetic field BeNuclear magnetic field BnElectron polarizability Pz eA slowdown factor Q, etc.
5. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as claimed in claim 1, wherein the step 4 comprises: and (3) after the basic parameters of the SERF atomic spin gyroscope are obtained through calibration in the steps 1 to 3, constructing a system matrix A and an input matrix B of the SERF atomic spin gyroscope.
6. The method of claim 1The method for improving the response speed of the SERF atomic spin gyroscope based on transient response calculation is characterized in that the step 5 comprises the following steps: outputting X-axis signals S when SERF atomic spin gyroscope worksxY-axis signal output SyConversion coefficient K divided by X-axis and Y-axis signal output and transverse electric polarizability respectivelyx、KyObtaining the X-axis electron polarizability Px eAnd Y-axis electron polarizability
Figure FDA0003425804640000021
7. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as recited in claim 1, wherein the step 6 comprises: the obtained X-axis electron polarizability
Figure FDA0003425804640000022
And Y-axis electron polarizability
Figure FDA0003425804640000023
Calculating the relative term of the nuclear spin transverse polarizability and the electron polarizability according to the system parameters
Figure FDA0003425804640000024
8. The method for improving the response speed of a SERF atomic spin gyroscope based on transient response calculation as recited in claim 1, wherein the step 7 comprises: the obtained X-axis electron nucleus polarizability
Figure FDA0003425804640000025
And polarizability of Y-axis electron nuclei
Figure FDA0003425804640000026
Calculating with the obtained system matrix A and input matrix B to obtain angular velocity input omegaxAnd omegay。
CN202111580141.1A 2021-12-22 2021-12-22 Method for improving response speed of SERF atomic spin gyroscope based on transient response calculation Active CN114440853B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111580141.1A CN114440853B (en) 2021-12-22 2021-12-22 Method for improving response speed of SERF atomic spin gyroscope based on transient response calculation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111580141.1A CN114440853B (en) 2021-12-22 2021-12-22 Method for improving response speed of SERF atomic spin gyroscope based on transient response calculation

Publications (2)

Publication Number Publication Date
CN114440853A true CN114440853A (en) 2022-05-06
CN114440853B CN114440853B (en) 2023-08-25

Family

ID=81363395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111580141.1A Active CN114440853B (en) 2021-12-22 2021-12-22 Method for improving response speed of SERF atomic spin gyroscope based on transient response calculation

Country Status (1)

Country Link
CN (1) CN114440853B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115265512A (en) * 2022-07-28 2022-11-01 北京航空航天大学 Rotation modulation method for SERF (spin exchange fiber) atomic spin gyroscope

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901939A (en) * 2012-10-16 2013-01-30 北京航空航天大学 Precise control method of atom spin SERF (Self-Exchange Relaxation-Free) state for stabilizing atom spin device
CN107515000A (en) * 2017-08-21 2017-12-26 北京航空航天大学 A kind of modulation system twin shaft atomic spin gyroscope
CN108445428A (en) * 2018-04-11 2018-08-24 北京航空航天大学 A kind of SERF atom magnetometers electronic polarizability measurement method
CN109373989A (en) * 2018-10-12 2019-02-22 北京航空航天大学 A kind of SERF atomic spin gyro nuclear spin self compensation point closed loop control method
CN112556678A (en) * 2020-11-24 2021-03-26 北京航空航天大学 Method for measuring nuclear polarizability of atomic spin gyroscope based on adiabatic fast channel

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901939A (en) * 2012-10-16 2013-01-30 北京航空航天大学 Precise control method of atom spin SERF (Self-Exchange Relaxation-Free) state for stabilizing atom spin device
CN107515000A (en) * 2017-08-21 2017-12-26 北京航空航天大学 A kind of modulation system twin shaft atomic spin gyroscope
CN108445428A (en) * 2018-04-11 2018-08-24 北京航空航天大学 A kind of SERF atom magnetometers electronic polarizability measurement method
CN109373989A (en) * 2018-10-12 2019-02-22 北京航空航天大学 A kind of SERF atomic spin gyro nuclear spin self compensation point closed loop control method
CN112556678A (en) * 2020-11-24 2021-03-26 北京航空航天大学 Method for measuring nuclear polarizability of atomic spin gyroscope based on adiabatic fast channel

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
裴宏宇 等: "基于数值微分计算的SERF原子自旋惯性测量动态仿真", 《导航定位与授时》 *
陈东营 等: "无自旋交换弛豫原子自旋陀螺非线性特性实验研究", 《光学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115265512A (en) * 2022-07-28 2022-11-01 北京航空航天大学 Rotation modulation method for SERF (spin exchange fiber) atomic spin gyroscope
CN115265512B (en) * 2022-07-28 2024-04-09 北京航空航天大学 Rotation modulation method for SERF atomic spin gyroscope

Also Published As

Publication number Publication date
CN114440853B (en) 2023-08-25

Similar Documents

Publication Publication Date Title
CN103438877B (en) A kind of inertia based on SERF atomic spin effect and magnetic field integral measurement method
CN109738837B (en) Residual magnetic field in-situ compensation method for single-beam SERF atomic magnetometer
JP5854735B2 (en) Nuclear magnetic resonance imaging apparatus and nuclear magnetic resonance imaging method
JP6222974B2 (en) Optical pumping magnetometer and magnetic sensing method
CN107192633A (en) Under a kind of SERF states in on-line measurement atom magnetometer air chamber alkali metal density method
CN108717168B (en) Scalar magnetic field gradient measuring device and method based on light field amplitude modulation
CN112924910B (en) In-situ magnetometer-based method for measuring residual magnetism in shielding barrel
CN113075594B (en) Electronic polarizability double-shaft in-situ measurement system and method for SERF atomic magnetometer
CN106017451A (en) High-precision measurement method for magnetic field compensation inertial angular rate based on of SERF (spin-exchange relaxation free) atomic device
CN103969604A (en) Radio-frequency atom magnetometer and method for measuring nuclear magnetic resonance (NMR) signal by same
CN102721827B (en) Optical fiber accelerometer
CN107525946A (en) Acceleration measurement method and device based on atomic interference in optical waveguide
CN109631959B (en) Atomic spin precession detection device and method based on optical fiber Sagnac interference
CN114440853B (en) Method for improving response speed of SERF atomic spin gyroscope based on transient response calculation
CN107490775B (en) Triaxial coil constant and non-orthogonal angle integrated measurement method
Jiang et al. Techniques for measuring transverse relaxation time of xenon atoms in nuclear-magnetic-resonance gyroscopes and pump-light influence mechanism
CN111060853B (en) Three-dimensional magnetic field in-situ measurement method based on electron paramagnetic resonance-nuclear magnetic resonance
CN111964658B (en) Nuclear magnetic resonance gyroscope closed-loop magnetic resonance method driven by rotating field
CN108534770B (en) A kind of129Xe-Rb spin exchange rate rapid measurement method
CN111060747A (en) High-sensitivity nuclear spin precession detection method based on electron spin
CN116626558A (en) Magnetic field gradient error analysis and suppression method based on zero bias magnetic field sensitivity
CN115856727A (en) SERF system and method for simultaneously measuring rotation and magnetic field by using diatomic air chamber
Zou et al. A novel calibration method research of the scale factor for the all-optical atomic spin inertial measurement device
CN112683995B (en) Alkali metal electron polarizability measuring method based on mixed alkali metal
CN112683994B (en) Inert gas nuclear polarizability measuring method based on alkali metal inert gas mixing

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