CN106019179B - A kind of quick locking method of CPT magnetometers - Google Patents
A kind of quick locking method of CPT magnetometers Download PDFInfo
- Publication number
- CN106019179B CN106019179B CN201610339095.9A CN201610339095A CN106019179B CN 106019179 B CN106019179 B CN 106019179B CN 201610339095 A CN201610339095 A CN 201610339095A CN 106019179 B CN106019179 B CN 106019179B
- Authority
- CN
- China
- Prior art keywords
- point
- frequency
- searching
- cpt
- peak
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of quick locking methods of CPT magnetometers, Searching point is estimated using slope, the signal curve slope searched near starting point is calculated by starting step size, then the corresponding frequency values of Searching point are inferred according to the amplitude characteristic of curve, reduces operand, shortens search time, due to actual variable searching step, so the fine search suitable for signal peak, can improve search speed, improve the whole work efficiency of system.This method applies also for a variety of optical sensor signals search and locking.
Description
Technical field
The present invention relates to a kind of quick locking method of CPT magnetometers more particularly to magnetometer CPT peak-to-peak signal peak values are fixed
Position.
Background technology
With the development of high-resolution laser spectroscopy and quantum optical research, one based on atom and laser interaction
Serial nonlinear optical phenomena be used to prepare the atomic resonance signal of narrow linewidth, to realize the height to magnetic field using the signal
Precision measure.Wherein Coherent Population Trapping imprison (coherent population trapping, CPT) magnetometer is the one of the principle
Kind typical case.
In atom magnetometer system, atom optics sensor can detect alkali metal there are external magnetic field
The CPT signal curves of atomic medium, control system are handled the signal by the algorithm of setting, accurate to determine CPT letters
Number corresponding frequency values in peak of curve place, and by the frequency-splitting between calculation side peak and main peak come calculating magnetic field intensity.
Since the precision and speed of peak-seeking algorithm directly affect the resolution ratio and response speed that magnetometer detects magnetic field, the peaks CPT
The design and optimum choice of peak-seeking algorithm routine just seem extremely important.And as geomagnetic survey system wants precision and responsiveness
Continuous promotion is asked, higher requirements are also raised for the quick stationkeeping ability of signal-data processing and Frequency point to control system.
In general, the CPT signal frequencies search of CPT magnetometers is realized by Slope Method.Slope Method peak-seeking is basis
The slope variation rule of the peaks CPT model curve finds peak position.Selecting frequency deviation in CPT models is shown in Fig. 1
A series of slope values that dw is calculated when being -100~100.As can be seen that the small range frequency-splitting near peak value is calculated
The corresponding slope value gone out is gradually increased, and when frequency-splitting is just zero, that is, at peaking when obtains
Slope value is equal to 0, and the frequency values symbol of peak value left and right side is opposite.It can be found according to this changing rule of slope
The iteration speed situation of peak value, slope algorithm can be embodied by model experiment arranged below:
(1) starting point and fixed step size are set.Starting point in typical case may be configured as dw0=-100, fixed step size
For d_w=3.
(2) investigation is determined.Slope value D (H0) and D (H1) size for comparing at front and back 2 points, if D (H0)<D(H1),
Then direction is just.
(3) it iterates to calculate.By fixed step size Continuous plus frequency-splitting go forward side by side line direction judge, if what is currently calculated is oblique
Rate value and initial slope value jack per line, then continue positive searching, otherwise reversed to find.
The peak-seeking iteration result of the algorithm is as shown in Figure 2:
When operation iterates to 33 step, dw=0 has arrived at peak value at this time.Since the iterative algorithm remains
It is updated the data according to fixed step size, so entering frequency locking process at peaking with rear, the state of Frequency Locking is attached in dw=0
The nearly oscillation for keeping slight amplitude.Amplitude is equal to the step value of setting.The advantage of slope peak-seeking method is can be according to practical feelings
Condition is arranged suitable step-length and quickly positions peak value.But if continuing to keep scanning mode after reaching target peak,
Appear in the fixed value oscillation at peak value.The the speed of peak-seeking the fast, and oscillation when being maintained at desired value is more violent.
Invention content
The technology of the present invention solves the problems, such as:Overcome the shortcomings of existing algorithm, a kind of quick frequency locking side of CPT magnetometers is provided
Method, can be according to initiating searches point come next Searching point in quick positioning searching section, more to the hunting speed of peak of curve
Soon, stability to the search precision higher of peak point, and after locking also higher.
The technical scheme is that:A kind of quick locking method of CPT magnetometers, steps are as follows:
Step 101:Establish alkali metal atom CPT peak-to-peak signal curve mathematic modelsWherein, H is light work(
Rate, Δ f are frequency-splitting, Δ f=f-f0, wherein f is the current frequency actually measured, frequency centered on f0, H0For power background amplitude, h
For signal amplitude, v is line width;IfTake the first approximation model of the mathematical model:
Step 102:Observed quantity is set as H ',
Step 103:Given initial value Δ f0, with the point for the first Searching point, corresponding observed quantity is
Step 104:Give positive increment δ f0, there is Δ f at this time1=Δ f0+δf0, which is the second Searching point, corresponding observation
Amount is
Step 105:Judge direction;If H ' (Δ f1)>H′(Δf0), the direction of search and ideal orientation are on the contrary, be anti-at this time
To Δ f need to be chosen again1=Δ f0-δf0, and seek againIf H ' (Δ f1)<H′(Δf0), table
Bright front direction of working as is forward direction, then is not necessarily to carry out operation bidirectional;
Step 106:Straight line is constructed, i.e., by two Searching points (Δ f0, H ' (Δ f0)) and (Δ f1, H ' (Δ f1)) construction straight line
Equation is
Step 107:Seek next Searching point on the constructed straight line of step 106;Work as H0For known quantity when, be located at next
A Frequency point, have in linear equation H ' (Δ f)=0, then the expression formula of the point be:This point
Corresponding observed quantity isWork as H0For unknown quantity when, it is assumed that next Frequency point has H' in linear equation
(Δ f)=Hp=h, then the expression formula of the point beΔ f is judged at this timemWith Δ f1Symbol
Number, if the two jack per line, this puts corresponding observed quantity and is set asOtherwise, to HpIt carries out
Update, makesThen update Δ fmWith H'(Δs fm);
Step 108:Data update;Take Δ f1=Δ fm, H ' (Δ f1)=H ' (Δ fm) and Δ f0=Δ f1,H′(Δf0)=
H′(Δf1), step 106 is re-executed to step 108;When search point at peaking with after, into frequency-locked state.
The advantages of the present invention over the prior art are that:The method of the present invention on the one hand can be according to initiating searches point come fast
Next Searching point in fast positioning searching section, faster to the hunting speed of peak of curve, to the search precision higher of peak point,
And the stability also higher after locking.On the other hand, which remains to change with less in the case where initial step length is smaller
Peak value is locked for calculation times.The algorithm is the variant of Slope Method peak-seeking, the feelings known to the background value size of peak value model
Using particularly convenient under condition.
Description of the drawings
Fig. 1 is the slope value that frequency-splitting is corresponding to -100~100Hz ranges in the peaks CPT model;
Fig. 2 is the peak-seeking iteration situation (error amount 100, step-length 3) of Slope Method;
Fig. 3 is the algorithm performs experimental result (error amount 100, step-length 3) of the present invention;
Fig. 4 is the iteration situation experimental result (error amount 100, step-length 3) of the present invention.
Specific implementation mode
The present invention is a kind of quick locking method for CPT magnetometers, and steps are as follows:
Step 101:Establish alkali metal atom CPT peak-to-peak signal curve mathematic modelsWherein, H is
Luminous power, Δ f are frequency-splitting, Δ f=f-f0.Wherein f is the current frequency actually measured, f0Centered on frequency, H0For power
Background magnitude, h are signal amplitude, and v is line width, can use constant 500Hz, ifTake the first approximation of the mathematical model
Model:
Step 102:Observed quantity is set as H ',
Step 103:Given initial value Δ f0, with the point for the first Searching point, corresponding observed quantity is
Step 104:Give positive increment δ f0, there is Δ f at this time1=Δ f0+δf0, which is the second Searching point, corresponding observation
Amount is
Step 105:Judge direction.If H ' (Δ f1)>H′(Δf0), the direction of search and ideal orientation are on the contrary, be anti-at this time
To Δ f need to be chosen again1=Δ f0-δf0, and seek againIf H ' (Δ f1)<H′(Δf0), table
Bright front direction of working as is forward direction;
Step 106:Straight line is constructed, i.e., by two Searching points (Δ f0, H ' (Δ f0)) and (Δ f1, H ' (Δ f1)) constitute it is straight
Line equation is
Step 107:Seek next Searching point on the constructed straight line of step 106.This step is divided into two kinds of situations:(1) when
H0For known quantity when, be located at next Frequency point, have in linear equation H ' (Δ f)=0, then the expression formula of the point be:This puts corresponding observed quantity(2) work as H0For unknown quantity when, it is assumed that it is next
A Frequency point has H'(Δs f)=Hp=h in linear equation, then the expression formula of the point isAt this time
Judge Δ fmWith Δ f1Symbol, if the two jack per line, this puts corresponding observed quantity and is set as
Otherwise, to HpIt is updated, makesThen update Δ fmWith H'(Δs fm)。
Step 108:Data update;Take Δ f1=Δ fm, H ' (Δ f1)=H ' (Δ fm) and Δ f0=Δ f1,H′(Δf0)=
H′(Δf1), step 106 is re-executed to step 108;When Searching point at peaking with after, iterative process is lasting, into lock
Frequency state.
In a particular embodiment, it is assumed that in the ideal case, H0For known quantity when:
(1) the first approximation model of the peaks CPT model is:H=50, v=are set at this time
500;
(2) observed quantity is then set as H ',
(3) setting initial value Δ f0=-100, with the point for first Searching point;Corresponding to observed quantity is
(4) positive increment δ f are given0=3, there is Δ f at this time1=Δ f0+δf0, i.e. Δ f1=-97 points are second Searching point;
Corresponding observed quantity is
(5) due to H ' (Δ f1)<H′(Δf0), show that when front direction be forward direction;
(6) straight line is constructed.By two Searching points (Δ f0, H ' (Δ f0)) and (Δ f1, H ' (Δ f1)) constitute linear equation be
(7) next Searching point is sought.It is located at next Frequency point, there is H ' (Δ f)=0, the then point in linear equation
Expression formula isThis puts corresponding observed quantity
(8) data update.Take Δ f1=Δ fm, H ' (Δ f1)=H ' (Δ fm) and Δ f0=Δ f1, H ' (Δ f0)=H ' (Δs
f1), step (6) is re-executed to step (8).
Experimental result is as shown in Figure 3 and Figure 4.From experimental result as can be seen that Fig. 3, which is shown, is being added power background
The search process shown after 16950, Searching point quickly approaches peak point, subsequent slow-down when starting, until reaching peak
Value.Fig. 4 is shown when iterations are 11, has found peak value, has hereafter continued iteration, has been locked in peak position.
For fast algorithm when iteration chooses next close keyed end, closer to keyed end, then step-length is smaller, compares for Slope Method not
It can vibrate by step-length when reaching keyed end, relatively stablize.
Claims (1)
1. a kind of quick locking method of CPT magnetometers, it is characterised in that steps are as follows:
Step 101:Establish alkali metal atom CPT peak-to-peak signal curve mathematic modelsWherein, H is luminous power,
Δ f is frequency-splitting, Δ f=f-f0, wherein f is the current frequency actually measured, f0Centered on frequency, H0For power background amplitude, h
For signal amplitude, v is line width;IfTake the first approximation model of the mathematical model:
Step 102:Observed quantity is set as H ',
Step 103:Given initial value Δ f0, with the point for the first Searching point, corresponding observed quantity is
Step 104:Give positive increment δ f0, there is Δ f at this time1=Δ f0+δf0, which is the second Searching point, and corresponding observed quantity is
Step 105:Judge direction;If H ' (Δ f1)>H′(Δf0), at this time the direction of search and ideal orientation on the contrary, be it is reversed,
Δ f need to be chosen again1=Δ f0-δf0, and seek againIf H ' (Δ f1)<H′(Δf0), show to work as
Front direction is forward direction, then is not necessarily to carry out operation bidirectional;
Step 106:Straight line is constructed, i.e., by two Searching points (Δ f0, H ' (Δ f0)) and (Δ f1, H ' (Δ f1)) construction linear equation be
Step 107:Seek next Searching point on the constructed straight line of step 106;Work as H0For known quantity when, be located at next frequency
Rate point, have in linear equation H ' (Δ f)=0, then the expression formula of the point be:This point corresponds to
Observed quantity beWork as H0For unknown quantity when, it is assumed that next Frequency point has H'(Δs in linear equation
F)=Hp=h, then the expression formula of the point beΔ f is judged at this timemWith Δ f1Symbol,
If the two jack per line, this puts corresponding observed quantity and is set asOtherwise, to HpIt carries out more
Newly, makeThen update Δ fmWith H'(Δs fm);
Step 108:Data update;Take Δ f1=Δ fm, H ' (Δ f1)=H ' (Δ fm) and Δ f0=Δ f1,H′(Δf0)=H ' (Δs
f1), step 106 is re-executed to step 108;When search point at peaking with after, into frequency-locked state.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610339095.9A CN106019179B (en) | 2016-05-19 | 2016-05-19 | A kind of quick locking method of CPT magnetometers |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610339095.9A CN106019179B (en) | 2016-05-19 | 2016-05-19 | A kind of quick locking method of CPT magnetometers |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019179A CN106019179A (en) | 2016-10-12 |
CN106019179B true CN106019179B (en) | 2018-09-18 |
Family
ID=57095168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610339095.9A Active CN106019179B (en) | 2016-05-19 | 2016-05-19 | A kind of quick locking method of CPT magnetometers |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019179B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111044943B (en) * | 2019-12-24 | 2022-04-19 | 北京航天控制仪器研究所 | Multi-spectrum closed-loop locking method and system for CPT magnetometer |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101847994A (en) * | 2010-05-05 | 2010-09-29 | 中国科学院武汉物理与数学研究所 | Method and device for implementing Ramsey-CPT atomic frequency standard by microwave periodic On-Off modulation VCSEL |
WO2011079664A1 (en) * | 2009-12-31 | 2011-07-07 | 上海舜宇海逸光电技术有限公司 | System and method for detecting magneto-optic with optical fiber |
CN103744034A (en) * | 2013-12-30 | 2014-04-23 | 浙江大学 | Difference method for improving sensitivity and absolute precision of CPT (Coherent Population Trapping) atom magnetometer |
CN104317185A (en) * | 2014-11-21 | 2015-01-28 | 中国科学院武汉物理与数学研究所 | Method and device for obtaining Ramsey-CPT spectrum |
CN104410416A (en) * | 2014-11-27 | 2015-03-11 | 江汉大学 | Coherent population trapping atomic frequency standard and method for improving short-term stability thereof |
CN104502867A (en) * | 2014-12-11 | 2015-04-08 | 北京航天控制仪器研究所 | VCSEL (vertical cavity surface emitting laser) laser tube parameter automatic adjustment method for CPT (coherent population trapping) magnetometer system |
CN105049040A (en) * | 2015-06-24 | 2015-11-11 | 中国科学院武汉物理与数学研究所 | Method for correcting output frequency of CPT (Coherent Population Trapping) atomic clock through GNSS(Global Navigation Satellite System) |
-
2016
- 2016-05-19 CN CN201610339095.9A patent/CN106019179B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011079664A1 (en) * | 2009-12-31 | 2011-07-07 | 上海舜宇海逸光电技术有限公司 | System and method for detecting magneto-optic with optical fiber |
CN101847994A (en) * | 2010-05-05 | 2010-09-29 | 中国科学院武汉物理与数学研究所 | Method and device for implementing Ramsey-CPT atomic frequency standard by microwave periodic On-Off modulation VCSEL |
CN103744034A (en) * | 2013-12-30 | 2014-04-23 | 浙江大学 | Difference method for improving sensitivity and absolute precision of CPT (Coherent Population Trapping) atom magnetometer |
CN104317185A (en) * | 2014-11-21 | 2015-01-28 | 中国科学院武汉物理与数学研究所 | Method and device for obtaining Ramsey-CPT spectrum |
CN104410416A (en) * | 2014-11-27 | 2015-03-11 | 江汉大学 | Coherent population trapping atomic frequency standard and method for improving short-term stability thereof |
CN104502867A (en) * | 2014-12-11 | 2015-04-08 | 北京航天控制仪器研究所 | VCSEL (vertical cavity surface emitting laser) laser tube parameter automatic adjustment method for CPT (coherent population trapping) magnetometer system |
CN105049040A (en) * | 2015-06-24 | 2015-11-11 | 中国科学院武汉物理与数学研究所 | Method for correcting output frequency of CPT (Coherent Population Trapping) atomic clock through GNSS(Global Navigation Satellite System) |
Non-Patent Citations (2)
Title |
---|
基于CPT原理原子钟的频率控制***;张廷凯,等;《宇航计测技术》;20090831;第29卷(第4期);第36-40页 * |
小型CPT原子钟的数字电路优化研究;张旭,等;《宇航计测技术》;20150630;第35卷(第3期);第36-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106019179A (en) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11977956B2 (en) | Performing a calibration process in a quantum computing system | |
Peirce et al. | Optimal control of quantum-mechanical systems: Existence, numerical approximation, and applications | |
US7979143B2 (en) | Apparatus and method for proportional-integral-derivative control | |
EP2502317B1 (en) | Advanced laser wavelength control | |
CN106019179B (en) | A kind of quick locking method of CPT magnetometers | |
WO2022121809A1 (en) | Method and device for calibrating frequency of superconducting qubit, and readable storage medium | |
Magdziarz et al. | Limit theorems and governing equations for Lévy walks | |
EP2215637A1 (en) | Iterative feedback tuning in a scanning probe microscope | |
US20150248995A1 (en) | Radiofrequency adjustment for instability management in semiconductor processing | |
JP6304461B1 (en) | Motor control device | |
CN104865889A (en) | Numerical Control Device Having Function Of Calculating Frequency Characteristic Of Control Loop | |
JP2013213721A (en) | Object movement controlling device and fourier transform infrared spectrophotometer | |
RU2619815C1 (en) | Method and system for compensating drift of solid wave gyro | |
WO2023216663A1 (en) | Method and apparatus for determining rod position of control rod of pressurized water reactor, and rod position measurement system | |
TWI793371B (en) | Devices and methods for coating control using forward parameter correction and enhanced reverse engineering and related computer-readable medium | |
CN107065537B (en) | A kind of horizontal low speed that pushes away of AUV is without the motion control method in the case of output | |
US20220261679A1 (en) | Probe of qubit-resonator dispersive shift using ac stark shift | |
CN106685605A (en) | Method for detecting weak signal frequency value in background of strong noise by chaotic synchronization | |
JP2014135851A (en) | Resonant frequency search device, resonant frequency search method and program | |
CN113432608B (en) | Generalized high-order CKF method based on maximum correlation entropy for INS/CNS integrated navigation system | |
CN112600070B (en) | Laser accurate frequency stabilization system and method for SERF magnetometer | |
CN108225373B (en) | Large misalignment angle alignment method based on improved 5-order cubature Kalman | |
CN112068423B (en) | Design method of finite power attraction repetitive controller | |
CN109708669B (en) | Accelerometer signal processing system and method | |
Doncescu et al. | Displacement estimation with an optical feedback interferometer using an evolutionary algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |