CN110108302A - A method of atomic group is improved to throwing precision - Google Patents
A method of atomic group is improved to throwing precision Download PDFInfo
- Publication number
- CN110108302A CN110108302A CN201910422528.0A CN201910422528A CN110108302A CN 110108302 A CN110108302 A CN 110108302A CN 201910422528 A CN201910422528 A CN 201910422528A CN 110108302 A CN110108302 A CN 110108302A
- Authority
- CN
- China
- Prior art keywords
- image
- cold atom
- magneto
- coordinate
- optical trap
- 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
Links
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
This application discloses a kind of raising atomic groups to the method for throwing precision, it include Magneto-Optical Trap and photosensitive camera in cold atom interferometer, comprising: step 1 suitable for cold atom interferometer, when the cooling light of Magneto-Optical Trap is opened and gradient magnetic is closed, the background image of Magneto-Optical Trap is obtained using photosensitive camera；Step 2, using photosensitive camera, gradually after cooling light and gradient magnetic are opened and fall into prisoner's cold atom cloud, multiple target images of Magneto-Optical Trap are obtained；Step 3, according to background image and multiple target images, using difference algorithm, the imprison position of cold atom cloud for every target image is calculated；Step 4, the coordinate jitter error for imprisoning position in multiple target images is calculated, and according to coordinate jitter error, calculates the first adjustment coefficient of Magneto-Optical Trap.By the technical solution in the application, be conducive to the accuracy for improving determining cold atom cloud imprison position, optimize the interference effect of cold atom interferometer and the accuracy of parameter adjustment.
Description
Technical field
This application involves the technical fields of cold atom interference gyro instrument, in particular to a kind of raising atomic group to throwing
The method of precision.
Background technique
Interference is formed using the fluctuation of atom, and further constitutes the interferometric loop for surrounding certain area, utilizes Sa
Lattice Neck effect measures the atomic gyroscope of velocity of rotation and acceleration, and resolution ratio and measurement accuracy can than optical gyroscope
To greatly improve.Currently, there are mainly two types of such atomic gyroscopes: hot atomic beam gyroscope and cold atom gyroscope.Cold atom
Gyroscope needs to reduce atomic group cooling temperature, to improve the stability of cold atom cloud imprison position, and adjusts cold atom cloud pair
It throws pointing accuracy, atomic group and the detection performance of cold atom gyroscope is improved to throwing speed.
In practical work process, cold atom cloud initial position locating before casting every time can be because of the number of Trapping of Atoms
Mesh, imprison light polarization, imprison light frequency, imprison luminous intensity, the slight change of factors of imprison magnetic field strength etc. and it is different, and
The direction that cold atom cloud is cast every time also will receive the influence of above-mentioned factor.Initial position and the deviation for casting direction will lead to original
The repeatability of sub- interference effect is deteriorated, and the detection error of gyroscope is caused to increase.
And in the prior art, in order to monitor the initial position of cold atom cloud, high resolution camera can be used in prisoner
The cold atom cloud of taboo state is taken pictures, and then according to the related algorithm of image procossing, analyzes the variation of cold atom cloud position.?
During this, on the one hand, when determining cold atom cloud mass center, only treat cold atom cloud as common point source hot spot, do not have
There is the decaying for considering cold atom cloud brightness, leading to its center coordination, there are deviations.On the other hand, judging that cold atom cloud flight is inclined
It is indirectly counter to push away cold atom cloud usually according to the fluorescence signal change information of detection phase after interference fringe and interference when poor
Projectile deviation, this deviation judgment method lacks the monitoring to cold atom cloud intermediate flight process, it is difficult to judge to fly quasi-ly
The deviation of row process and the reason of cause the deviation, is unfavorable for the reduction of deviation.
Summary of the invention
The purpose of the application is: the accuracy for determining cold atom cloud imprison position is improved, optimization cold atom interferometer
The accuracy of interference effect and parameter adjustment.
The technical solution of the application is: provide a kind of raising atomic group to the method for throwing precision, it is dry suitable for cold atom
Interferometer includes Magneto-Optical Trap and photosensitive camera in cold atom interferometer, and the imaging area of photosensitive camera is right against Magneto-Optical Trap, this method
Include: step 1, when the cooling light of Magneto-Optical Trap is opened and gradient magnetic is closed, the background of Magneto-Optical Trap is obtained using photosensitive camera
Image；Step 2, using photosensitive camera, gradually after cooling light and gradient magnetic are opened and fall into prisoner's cold atom cloud, magneto-optic is obtained
Multiple target images of trap；Step 3, it is calculated using difference algorithm for every mesh according to background image and multiple target images
Logo image, cold atom cloud imprison position, wherein imprison position is the position that Magneto-Optical Trap imprisons cold atom cloud；Step 4, it counts
The coordinate jitter error that position is imprisoned in multiple target images is calculated, and according to coordinate jitter error, calculate Magneto-Optical Trap first is adjusted
Integral coefficient, the first adjustment coefficient include the first optical path coefficient and the first magnetic field coefficient, wherein the meter of the first adjustment coefficient f ()
Calculate formula are as follows:
σpos(u)=f (P (u), A (u), F (u), C (u), G (u)), u=1,2,3 ... U
In formula, σposIt (u) is the u corresponding coordinate jitter errors of target image, P (u) is that u target images are corresponding
Cooling light light intensity, A (u) is the projectile angle of corresponding cold atom cloud, and F (u) is corresponding cooling light frequency, and C (u) is phase
The field supply intensity answered, G (u) are corresponding magnetic field gradient.
It in any of the above-described technical solution, further, specifically includes in step 3: step 31, according to background image, pressing
According to the sequence of target image, a target image is successively chosen, using difference algorithm, calculates the pixel of pixel in target image
Revision value generates revision image, wherein the calculation formula of pixel revision value are as follows:
In formula, PtFor pixel revision value, P0For the pixel value of pixel in background image, PpFor picture in pth target image
The pixel value of vegetarian refreshments, p=1,2 ..., N；
Step 32, the maximum pixel of pixel value in revision image is successively chosen, is denoted as central point, and with the center of selection
Centered on point, the image-region ROI of revision image is determined, wherein image-region includes whole hot spot；Step 33, according to figure
As the pixel value in the ROI of region, the facula mass center of image-region ROI is calculated, and facula mass center is denoted as to the imprison of target image
Position.
It is further, specific in step 33 further include: step 331, with facula mass center in any of the above-described technical solution
(x0,y0) it is coordinate origin, equally spacedly construct the lateral bisector and longitudinal bisector of image-region ROI, wherein laterally etc.
Separated time and the item number K of longitudinal bisector are odd number, and (K-1)/2 article lateral bisector and (K-1)/2 article longitudinal bisector
By facula mass center；Step 332, it is horizontal to be calculated using Gauss curve fitting model according to the pixel value in image-region ROI for any bar
Transverse center coordinate in bisector, and calculate longitudinal center's coordinate in the bisector of any bar longitudinal direction, wherein Gauss is quasi-
The calculation formula of molding type are as follows:
In formula, B is pixel value, and w is the coordinate value of lateral bisector or longitudinal bisector, and a is amplitude, and b is transverse center
Coordinate or longitudinal center's coordinate, c are fit radius；Step 333, according to the average value of transverse center coordinate and longitudinal center's coordinate
Average value, generate spot center, by spot center be denoted as imprison position.
In any of the above-described technical solution, further, this method further include: step 5, after falling into prisoner's cold atom cloud, gradually
The projectile image of cold atom cloud is obtained, position coordinates of the cold atom cloud in projecting process are calculated, according to position coordinates, fitting life
At trajectory of projection of the cold atom cloud in projecting process；Step 6, the track registration of trajectory of projection is calculated, and according to track weight
The right second adjustment coefficient for calculating Magneto-Optical Trap, wherein second adjustment coefficient includes the second optical path coefficient and the second magnetic field coefficient,
The calculation formula of second adjustment coefficient g () are as follows:
G (Δ (u '), T1 (u '), T2 (u '))=σdiff(u '), u '=1,2,3 ... U '
In formula, σdiff(u ') is the track registration of the u ' articles trajectory of projection, and Δ (u ') is that corresponding cooling light frequency is lost
Harmonic quantity, T1 (u ') are corresponding cooling light light application time, and T2 (u ') is the corresponding magnetic fields time.
The beneficial effect of the application is:
Due to improving the positioning accuracy to atomic group initial position, and can be according to atomic group in flight course in difference
The position at moment calculates heading, flying speed and the diffusion velocity of atomic group, thus the technical solution of the application can
Atomic group imprison positional stability, atomic group are accurately measured to throwing pointing accuracy, atomic group to throwing rate uniformity, Jin Erti
The performance of high gyroscope.And by taking pictures the moment for control camera, above-mentioned measurement can be obtained within an airborne period
As a result, greatly improving conventional efficient.In addition, the technical solution of the application can also be according to measurement result in real time to optical path, magnetic
Field system carries out feedback control, the projecting process of adjusting and optimizing atomic group.
Detailed description of the invention
The advantages of above-mentioned and/or additional aspect of the application, will become bright in combining description of the following accompanying drawings to embodiment
It shows and is readily appreciated that, in which:
Fig. 1 is the schematic flow diagram according to the raising atomic group of one embodiment of the application to the method for throwing precision；
Fig. 2 is the curve analogous diagram of the target image and revision image according to one embodiment of the application；
Fig. 3 is the Gaussian curve analogous diagram according to one embodiment of the application；
Fig. 4 is the schematic diagram according to the bisector of one embodiment of the application；
Fig. 5 is the bisector pixel value simulation curve figure according to one embodiment of the application；
Fig. 6 is the function curve between light intensity and jitter error function according to one embodiment of the application.
Specific embodiment
It is with reference to the accompanying drawing and specific real in order to be more clearly understood that the above objects, features, and advantages of the application
Mode is applied the application is further described in detail.It should be noted that in the absence of conflict, the implementation of the application
Feature in example and embodiment can be combined with each other.
In the following description, many details are elaborated in order to fully understand the application, still, the application may be used also
To be implemented using other than the one described here other modes, therefore, the protection scope of the application is not by described below
Specific embodiment limitation.
As shown in Figure 1, present embodiments provide it is a kind of raising atomic group to throw precision method, be suitable for cold atom interfere
Instrument includes Magneto-Optical Trap, photosensitive camera in cold atom interferometer, wherein the imaging area of photosensitive camera is right against Magneto-Optical Trap, cold original
The both ends of sub- interferometer are provided with Magneto-Optical Trap, and Magneto-Optical Trap is generally made of three-dimensional space stationary field and gradient magnetic, wherein three-dimensional
Space stationary field is mutually perpendicular to (six beams) by three two-by-two, with specific circular polarization state bear it is detuning formed to penetrating laser beam, and
This is denoted as cooling light (wherein having coupled time pump light) to laser beam is penetrated, gradient magnetic is produced by a pair of reversed Helmholtz coil
It is raw.When the cooling light of six beams and gradient magnetic are opened, the area Xian Fu can be formed in Magneto-Optical Trap, center can fall into prisoner's atom, be formed
Cold atom cloud, wherein the position that cold atom cloud is formed is the initial position of cold atom cloud.It will fall into prisoner's in the Magneto-Optical Trap of both ends again
Cold atom cloud is carried out to throwing, to realize the interference of cold atom cloud.The initial position of cold atom cloud is directly shot by photosensitive camera
When, it will receive the interference of cooling light in Magneto-Optical Trap, therefore, this method comprises:
Step 1, it is opened in the cooling light of Magneto-Optical Trap and gradient magnetic is closed, i.e., Magneto-Optical Trap is fallen into before prisoner's cold atom cloud, benefit
The background image that Magneto-Optical Trap is obtained with photosensitive camera, at this point, not will form cold atom cloud in Magneto-Optical Trap, but is cooled light
It influences, still cooled illumination is bright for the partial region in Magneto-Optical Trap, and this partial dot bright area, can be to the initial of cold atom cloud
The determination of position (imprison position) interferes.
Step 2, it using photosensitive camera, is gradually opened in the cooling light and gradient magnetic of Magneto-Optical Trap and falls into prisoner's cold atom cloud
Afterwards, multiple target images of Magneto-Optical Trap are obtained；
Specifically, in the step 2, in the case where the system parameter of cold atom interferometer is constant, it is contemplated that cold atom
There are waves for the parameters such as initial position, shape, size, the cold atom cloud brightness of cold atom cloud that group is formed in Magneto-Optical Trap
It is dynamic, therefore, cooling light and gradient magnetic are opened, atomic group is carried out in Magneto-Optical Trap to fall into prisoner, after falling into prisoner's cold atom cloud, utilizes sense
Light camera obtains a target image before cold atom cloud is cast in Magneto-Optical Trap, after cold atom cloud projectile, repeats again above-mentioned
Process obtains the target image of next Magneto-Optical Trap, i.e., throws after multiple formation cold atom clouds of acquisition Magneto-Optical Trap and in cold atom cloud
Target image before penetrating.
Step 3, according to background image and multiple target images, using difference algorithm, calculate it is for every target image,
The imprison position of cold atom cloud, wherein imprison position is the position that Magneto-Optical Trap imprisons cold atom cloud；
Further, it in the step 3, specifically includes:
Step 31, according to background image, according to the sequence of target image, a target image is successively chosen, using difference
Algorithm calculates the pixel revision value of pixel in target image, generates revision image, wherein the calculation formula of pixel revision value
Are as follows:
In formula, PtFor pixel revision value, P0For the pixel value of pixel in background image, PpFor picture in pth target image
The pixel value of vegetarian refreshments, p=1,2 ..., N.
Specifically, it sets shown in the pth target image such as Fig. 2 (a) got, due to the interference of cooling light, target figure
There is the higher noise spots 201 of a large amount of, brightness as in, are revised by above-mentioned pixel value, available pth target image
Revised revision image, as shown in Fig. 2 (b).
Step 32, the maximum pixel of pixel value in revision image is successively chosen, is denoted as central point, and with the center of selection
Centered on point, the image-region ROI of revision image is determined, wherein image-region includes whole hot spot, image-region ROI's
Size is m × n.
In a kind of implementation in the present embodiment, determine image-region ROI's method particularly includes:
As unit of any revision image, the maximum pixel of pixel value in the revision image is chosen, central point is denoted as,
The central point is denoted as to the center of hot spot marquee, according to the size of pixel value in the revision image, adjusts hot spot marquee
Size generates hot spot marquee, so that choosing the hot spot of generation comprising hot spot whole in revision image in hot spot marquee
Frame is denoted as the image-region ROI of revision image, wherein hot spot marquee can be in rectangle, square, ellipse and circle
One kind, in the present embodiment select rectangle as hot spot marquee.
In another implementation in the present embodiment, determine image-region ROI's method particularly includes:
As unit of any revision image, and set pixel value threshold value.It is maximum to choose pixel value in the revision image
Pixel is denoted as central point, which is denoted as to the center of hot spot marquee, according to the size of pixel value in the revision image
With the pixel value threshold value of setting, selected pixels value is greater than or equal to the pixel of pixel value threshold value, and the pixel selected is remembered
Make hot spot, adjust the size of hot spot marquee, so that generating hot spot comprising hot spot whole in revision image in hot spot marquee
The hot spot marquee of generation is denoted as the image-region ROI of revision image by marquee.Wherein, hot spot marquee can for rectangle,
One of square, ellipse and circle, select rectangle as hot spot marquee in the present embodiment.
Step 33, according to the pixel value in image-region ROI, the facula mass center of image-region ROI is calculated, and by hot spot matter
The heart is denoted as the imprison position of target image；
Specifically, in the step 33, facula mass center (x0,y0) calculation formula are as follows:
In formula, F (i, j) is the pixel value of pixel (i, j) in image-region ROI, and m, n are respectively image-region ROI
Width and height.
It is analyzed by the pixel value to the pixel of identical abscissa or identical ordinate in image-region ROI, into
It crosses a large amount of data statistics to obtain, for identical abscissa (ordinate), as being gradually increased for its ordinate (abscissa) (subtracts
It is small), the variation of pixel value meets Gaussian Profile, as shown in Fig. 3 (a), therefore, attempts to use Gauss curve fitting model, in hot spot
On the basis of mass center, the imprison position of image-region ROI is revised, then cold atom interferometer is adjusted by revised imprison position
Parameter so that the detection error of cold atom interferometer further decreases, wherein the brightness of revised image-region ROI is bent
Shown in line such as Fig. 3 (b).
Therefore, in a preferred embodiment of the present embodiment, in the step 33 further include:
Step 331, with facula mass center (x0,y0) it is coordinate origin, equally spacedly construct the lateral equal part of image-region ROI
Line and longitudinal bisector, wherein the item number K of lateral bisector and longitudinal bisector is odd number, and (K-1)/2 article lateral equal part
Line and (K-1)/2 article longitudinal bisector pass through facula mass center；
Step 332, any bar transverse direction equal part is calculated using Gauss curve fitting model according to the pixel value in image-region ROI
Transverse center coordinate on line, and
Calculate longitudinal center's coordinate in the bisector of any bar longitudinal direction, wherein the calculation formula of Gauss curve fitting model are as follows:
In formula, B is pixel value, and w is the coordinate value of lateral bisector or longitudinal bisector, and a is amplitude, and b is transverse center
Coordinate or longitudinal center's coordinate, c are fit radius；
Step 333, according to the average value of the average value of transverse center coordinate and longitudinal center's coordinate, spot center is generated,
Spot center is denoted as imprison position.
Specifically, in the present embodiment, the quantity for setting lateral bisector and longitudinal bisector is respectively 7, composition
Bisector is as shown in figure 4, take by coordinate origin (x0,y0) horizontal line of (facula mass center), between three etc. above coordinate origin
Horizontal line, total 7 lateral bisectors take by coordinate origin (hot spot matter three below horizontal line and coordinate origin at equal intervals
The heart) vertical line, three of coordinate origin right three of vertical line and coordinate origin left vertical lines at equal intervals at equal intervals, amount to 7
Longitudinal bisector.
By taking lateral bisector as an example, in image-region ROI, the value range of the abscissa i of lateral bisector is 1,
2 ..., m, for each transverse direction bisector, any one abscissa i, pixel value B fixed corresponding to one for determining, therefore,
By Gauss curve fitting model can corresponding brightness curve can find out the transverse direction of 7 lateral bisectors as shown in figure 5, therefore
Its average value is denoted as the abscissa x of imprison position by centre coordinate b1, b2 ..., b7center.Similarly, imprison position can be calculated
The ordinate y setcenter.To sum up, the coordinate for imprisoning position is (xcenter,ycenter)。
Step 4, the coordinate jitter error that position is imprisoned in multiple target images is calculated, and according to coordinate jitter error, meter
The first adjustment coefficient of Magneto-Optical Trap is calculated, the first adjustment coefficient includes the first optical path coefficient and the first magnetic field coefficient, wherein first adjusts
The calculation formula of integral coefficient f () are as follows:
σpos(u)=f (P (u), A (u), F (u), C (u), G (u)), u=1,2,3 ... U
In formula, σposIt (u) is the u corresponding coordinate jitter errors of target image, P (u) is that u target images are corresponding
Cooling light light intensity, A (u) is the projectile angle of corresponding cold atom cloud, and F (u) is corresponding cooling light frequency, and C (u) is phase
The field supply intensity answered, G (u) be corresponding magnetic field gradient, wherein the first optical path coefficient include cooling light light intensity P (u) and
Cooling light frequency F (u), the first magnetic field coefficient include the projectile angle A (u) and field supply intensity C (u) of cold atom cloud.
Specifically, in the step 4, coordinate jitter error σposCalculation formula are as follows:
In formula, σx=std (xcenter), it is abscissa xcenterStandard deviation in the direction x, σy=std (ycenter), it is vertical seat
Mark ycenterStandard deviation in the direction y.
In the present embodiment, using the method gradually scanned, coordinate jitter error σ is determinedposThe smallest parameter combination, will
This parameter combination is denoted as the first adjustment coefficient.By taking light intensity P (u) as an example, it is illustrated to the first adjustment coefficient is calculated, by first
Remaining coefficient in regulation coefficient immobilizes, and successively increases light intensity P (u), after obtaining U target images, constructs light intensity P (u)
With coordinate jitter error σposBetween function, as shown in Figure 6.Pass through the analysis to above-mentioned function, it can be deduced that min coordinates are trembled
Dynamic error σposUnder, corresponding light intensity P.And then by being repeatedly scanned with, obtain making coordinate jitter error σposThe smallest coefficient combination,
Determine the first adjustment coefficient.
It should be noted that some coefficients (such as cooling light frequency F) are can be by means of computer in actual mechanical process
Program realizes automatic scanning, and depends on if some coefficients (as cast angle A) and manually change scan values.
Further, method further include:
Step 5, after falling into prisoner's cold atom cloud, the projectile image of cold atom cloud is gradually obtained, cold atom cloud is calculated and was casting
Position coordinates in journey, according to position coordinates, fitting generates trajectory of projection of the cold atom cloud in projecting process；
Specifically, in the step 5, by image obtain method, set cold atom cloud in cold atom cloud interferometer in
Duration flight is 20ms, is taken pictures at interval of 1ms to cold atom cloud, using the method in above-mentioned steps 1 to 8, is calculated
Position coordinates in cold atom projecting process, then multiple position coordinates of the cold atom cloud got in 20ms are fitted,
Generate trajectory of projection.
Step 6, the track registration of trajectory of projection is calculated, and calculates the second adjustment system of Magneto-Optical Trap according to track registration
Number, wherein second adjustment coefficient includes the second optical path coefficient and the second magnetic field coefficient, and the calculating of second adjustment coefficient g () is public
Formula are as follows:
G (Δ (u '), T1 (u '), T2 (u '))=σdiff(u '), u '=1,2,3 ... U '
In formula, σdiff(u ') is the track registration of the u ' articles trajectory of projection, and Δ (u ') is that corresponding cooling light frequency is lost
Harmonic quantity, T1 (u ') are corresponding cooling light light application time, and T2 (u ') is the corresponding magnetic fields time, wherein the second optical path system
Number includes cooling light frequency mismatching angle Δ (u ') and cooling light light application time T1 (u '), when the second magnetic field coefficient includes magnetic fields
Between T2 (u ').
Specifically, it in the step 6, is fitted using trajectory of projection of the parabola to cold atom cloud.Theoretically, it throws
The vertex of object line is exactly the center of interference region, that is, corresponds to the effect position of the second beam Raman light in cold atom interferometer
It sets.Due to the projectile mechanism of cold atom in cold atom interferometer, i.e., there are two cold atom clouds to cast in opposite directions respectively, therefore, by two
The distance between the parabolical vertex postion (A) of item and postion (B) are defined as track spacing, i.e. track spacing diff=
postion(A)-postion(B).In cold atom cloud projecting process, need to carry out repeatedly to cast test, and generate multipair
There is shake between trajectory of projection, therefore, the shake standard deviation of spacing diff in track between multipair trajectory of projection is denoted as track
Registration σdiff, i.e., for any bar trajectory of projection, its available track registration.
Using the method gradually scanned, track registration σ is determineddiffIn the case that (u ') is optimal, corresponding parameter combination,
By this parameter combination, it is denoted as second adjustment coefficient.
The technical solution for having been described in detail above with reference to the accompanying drawings the application, present applicant proposes a kind of raising atomic groups to throwing
The method of precision is suitable for cold atom interferometer, includes Magneto-Optical Trap and photosensitive camera in cold atom interferometer, comprising: step 1,
When the cooling light of Magneto-Optical Trap is opened and gradient magnetic is closed, the background image of Magneto-Optical Trap is obtained using photosensitive camera；Step 2,
Using photosensitive camera, gradually after cooling light and gradient magnetic are opened and fall into prisoner's cold atom cloud, multiple targets of Magneto-Optical Trap are obtained
Image；Step 3, it according to background image and multiple target images, using difference algorithm, calculates for every target image, cold
The imprison position of atomic group, imprison position are the position that Magneto-Optical Trap imprisons cold atom cloud；Step 4, it calculates in multiple target images
The coordinate jitter error of position is imprisoned, and according to coordinate jitter error, calculates the first adjustment coefficient of Magneto-Optical Trap.Pass through the application
In technical solution, be conducive to improve the accuracy of determining cold atom cloud imprison position, improve the interference effect of cold atom interferometer
The accuracy of fruit and parameter adjustment.
Step in the application can be sequentially adjusted, combined, and deleted according to actual needs.
Unit in the application device can be combined, divided and deleted according to actual needs.
Although disclosing the application in detail with reference to attached drawing, it will be appreciated that, these descriptions are only exemplary, not
For limiting the application of the application.The protection scope of the application may include not departing from this Shen by appended claims
It please be in the case where protection scope and spirit for various modifications, remodeling and equivalent scheme made by inventing.
Claims (4)
1. a kind of atomic group that improves is to the method for throwing precision, which is characterized in that be suitable for cold atom interferometer, the cold atom is dry
It include Magneto-Optical Trap and photosensitive camera in interferometer, the imaging area of the photosensitive camera is right against the Magneto-Optical Trap, this method comprises:
Step 1, when the cooling light of the Magneto-Optical Trap is opened and gradient magnetic is closed, the magnetic is obtained using the photosensitive camera
The background image of ligh trap；
Step 2, it using the photosensitive camera, is gradually opened in the cooling light and the gradient magnetic and falls into prisoner's cold atom cloud
Afterwards, multiple target images of the Magneto-Optical Trap are obtained；
Step 3, it is calculated using difference algorithm for mesh described in every according to the background image and multiple described target images
Logo image, cold atom cloud imprison position, wherein the imprison position is the position that the Magneto-Optical Trap imprisons the cold atom cloud
It sets；
Step 4, the coordinate jitter error for imprisoning position described in multiple described target images is calculated, and is shaken according to the coordinate
Error, calculates the first adjustment coefficient of the Magneto-Optical Trap, and the first adjustment coefficient includes the first optical path coefficient and the first magnetic field
Coefficient, wherein the calculation formula of the first adjustment coefficient f () are as follows:
σpos(u)=f (P (u), A (u), F (u), C (u), G (u)), u=1,2,3 ... U
In formula, σpos(u) the corresponding coordinate jitter error of the target image is opened for u, P (u) is the u targets
The light intensity of the corresponding cooling light of image, A (u) are the projectile angle of the corresponding cold atom cloud, and F (u) is corresponding cold
But light frequency, C (u) are corresponding field supply intensity, and G (u) is corresponding magnetic field gradient.
2. improving atomic group as described in claim 1 to the method for throwing precision, which is characterized in that specifically included in step 3:
Step 31, according to the background image, according to the sequence of the target image, a target image is successively chosen,
Using the difference algorithm, the pixel revision value of pixel in the target image is calculated, generates revision image, wherein described
The calculation formula of pixel revision value are as follows:
In formula, PtFor the pixel revision value, P0For the pixel value of pixel in the background image, PpFor pth Zhang Suoshu target
The pixel value of pixel in image, p=1,2 ..., N；
Step 32, the maximum pixel of pixel value in the revision image is successively chosen, is denoted as central point, and with selected institute
Centered on stating central point, the image-region ROI of the revision image is determined, wherein described image region includes whole hot spot；
Step 33, according to the pixel value in the ROI of described image region, the facula mass center of described image region ROI is calculated, and by institute
State the imprison position that facula mass center is denoted as the target image.
3. improving atomic group as claimed in claim 2 to the method for throwing precision, which is characterized in that specifically also wrapped in step 33
It includes:
Step 331, with the facula mass center (x0, y0) it is coordinate origin, equally spacedly construct the transverse direction of described image region ROI
Bisector and longitudinal bisector, wherein the item number K of the transverse direction bisector and longitudinal bisector is odd number, and (K-
1)/2 article lateral bisector and (K-1)/2 article longitudinal bisector pass through the facula mass center；
Step 332, it according to the pixel value in the ROI of described image region, using Gauss curve fitting model, calculates described in any bar laterally
Transverse center coordinate in bisector, and
Calculate longitudinal center's coordinate in longitudinal direction bisector described in any bar, wherein the calculation formula of the Gauss curve fitting model
Are as follows:
In formula, B is the pixel value, and w is the coordinate value of the lateral bisector or longitudinal bisector, and a is amplitude, and b is
Transverse center coordinate or longitudinal center's coordinate, c are fit radius；
Step 333, it according to the average value of the average value of the transverse center coordinate and longitudinal center's coordinate, generates in hot spot
The spot center is denoted as the imprison position by the heart.
4. improving atomic group as described in claim 1 to the method for throwing precision, which is characterized in that this method further include:
Step 5, after falling into prisoner's cold atom cloud, the projectile image of the cold atom cloud is gradually obtained, the cold atom cloud is calculated
Position coordinates in projecting process, according to the position coordinates, fitting generates throwing of the cold atom cloud in projecting process
Penetrate track；
Step 6, the track registration of the trajectory of projection is calculated, and calculates the of the Magneto-Optical Trap according to the track registration
Two regulation coefficients, wherein the second adjustment coefficient includes the second optical path coefficient and the second magnetic field coefficient, the second adjustment system
The calculation formula of number g () are as follows:
G (Δ (u '), T1 (u '), T2 (u '))=σdiff(u '), u '=1,2,3 ... U '
In formula, σdiff(u ') is the track registration of trajectory of projection described in the u ' articles, and Δ (u ') is that corresponding cooling light frequency is detuning
Amount, T1 (u ') are corresponding cooling light light application time, and T2 (u ') is the corresponding magnetic fields time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910422528.0A CN110108302B (en) | 2019-05-21 | 2019-05-21 | Method for improving atom group polishing precision |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910422528.0A CN110108302B (en) | 2019-05-21 | 2019-05-21 | Method for improving atom group polishing precision |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110108302A true CN110108302A (en) | 2019-08-09 |
| CN110108302B CN110108302B (en) | 2023-01-17 |
Family
ID=67491345
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910422528.0A Active CN110108302B (en) | 2019-05-21 | 2019-05-21 | Method for improving atom group polishing precision |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110108302B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111323742A (en) * | 2020-03-21 | 2020-06-23 | 哈尔滨工程大学 | Phase interferometer based on curve array and direction finding method thereof |
| CN112729298A (en) * | 2020-12-24 | 2021-04-30 | 华中光电技术研究所（中国船舶重工集团公司第七一七研究所） | Multi-gyroscope atomic group position synchronous control method |
| CN116738590A (en) * | 2023-04-03 | 2023-09-12 | 中国计量科学研究院 | Optimization method and device for cold atom magneto-optical trap system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8237105B1 (en) * | 2011-02-03 | 2012-08-07 | Northrop Grumman Guidance & Electronics Company, Inc. | Magneto-optical trap for cold atom beam source |
| CN103134949A (en) * | 2013-01-08 | 2013-06-05 | 上海卫星工程研究所 | Satellite-borne cold atom imprisoning-type acceleration measurement method |
| CN103472494A (en) * | 2013-09-23 | 2013-12-25 | 中国科学院武汉物理与数学研究所 | Sensor and method for measuring gravitational potential three-order differential quotient based on atom interference effect |
| US8853613B1 (en) * | 2013-09-24 | 2014-10-07 | Honeywell International Inc. | Magnetic field coils for magneto-optical trap |
| CN203881958U (en) * | 2014-06-12 | 2014-10-15 | 中国科学院武汉物理与数学研究所 | Level gravity gradient measurement sensor based on cold atomic beam interferometer |
| CN104537656A (en) * | 2014-12-22 | 2015-04-22 | 中国科学院国家授时中心 | Detection method for decentered Gaussian beams emitted by fiber beam expanding collimating lens barrel |
| CN104880614A (en) * | 2015-06-09 | 2015-09-02 | 华南师范大学 | Microwave electric field intensity meter based on cold Rydberg atom interferometer and measuring method thereof |
| CN107462234A (en) * | 2017-09-01 | 2017-12-12 | 中国科学院武汉物理与数学研究所 | It is a kind of that northern measuring apparatus and measuring method are sought based on cold atom interference technique |
| CN109708674A (en) * | 2018-12-12 | 2019-05-03 | 中国船舶重工集团公司第七一七研究所 | A kind of cold atom cloud control method that quickly sunken prisoner and high frequency are cast in opposite directions |
-
2019
- 2019-05-21 CN CN201910422528.0A patent/CN110108302B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8237105B1 (en) * | 2011-02-03 | 2012-08-07 | Northrop Grumman Guidance & Electronics Company, Inc. | Magneto-optical trap for cold atom beam source |
| CN103134949A (en) * | 2013-01-08 | 2013-06-05 | 上海卫星工程研究所 | Satellite-borne cold atom imprisoning-type acceleration measurement method |
| CN103472494A (en) * | 2013-09-23 | 2013-12-25 | 中国科学院武汉物理与数学研究所 | Sensor and method for measuring gravitational potential three-order differential quotient based on atom interference effect |
| US8853613B1 (en) * | 2013-09-24 | 2014-10-07 | Honeywell International Inc. | Magnetic field coils for magneto-optical trap |
| CN203881958U (en) * | 2014-06-12 | 2014-10-15 | 中国科学院武汉物理与数学研究所 | Level gravity gradient measurement sensor based on cold atomic beam interferometer |
| CN104537656A (en) * | 2014-12-22 | 2015-04-22 | 中国科学院国家授时中心 | Detection method for decentered Gaussian beams emitted by fiber beam expanding collimating lens barrel |
| CN104880614A (en) * | 2015-06-09 | 2015-09-02 | 华南师范大学 | Microwave electric field intensity meter based on cold Rydberg atom interferometer and measuring method thereof |
| CN107462234A (en) * | 2017-09-01 | 2017-12-12 | 中国科学院武汉物理与数学研究所 | It is a kind of that northern measuring apparatus and measuring method are sought based on cold atom interference technique |
| CN109708674A (en) * | 2018-12-12 | 2019-05-03 | 中国船舶重工集团公司第七一七研究所 | A kind of cold atom cloud control method that quickly sunken prisoner and high frequency are cast in opposite directions |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111323742A (en) * | 2020-03-21 | 2020-06-23 | 哈尔滨工程大学 | Phase interferometer based on curve array and direction finding method thereof |
| CN112729298A (en) * | 2020-12-24 | 2021-04-30 | 华中光电技术研究所（中国船舶重工集团公司第七一七研究所） | Multi-gyroscope atomic group position synchronous control method |
| CN116738590A (en) * | 2023-04-03 | 2023-09-12 | 中国计量科学研究院 | Optimization method and device for cold atom magneto-optical trap system |
| CN116738590B (en) * | 2023-04-03 | 2024-07-05 | 中国计量科学研究院 | Optimization method and device for cold atom magneto-optical trap system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110108302B (en) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110108302A (en) | A method of atomic group is improved to throwing precision | |
| US10901090B2 (en) | TOF camera system and a method for measuring a distance with the system | |
| Shatsky et al. | The mass ratio distribution of B-type visual binaries in the Sco OB2 association | |
| CN106767529B (en) | The automatic focusing method and system of laser facula identification and laser-projector | |
| KR102255017B1 (en) | Method for calibrating an image capture sensor comprising at least one sensor camera using a time coded pattern target | |
| CN108370438A (en) | The depth camera component of range gating | |
| WO2018196303A1 (en) | Projector calibration method and apparatus based on multi-directional projection | |
| CN106546216A (en) | Distance measuring method and device, camera and mobile terminal | |
| KR20070052756A (en) | Enhancement of aimpoint in simulated training systems | |
| GB2508958A (en) | Size measurement using synthesised range image | |
| CN109373897A (en) | A kind of measurement method based on laser virtual ruler | |
| CN108291962A (en) | Laser radar apparatus | |
| CN111025032B (en) | Aerial beam measuring system and method based on lift-off platform | |
| CN108426700B (en) | A kind of gravity is directed toward the detection method of influence on camera lens optical axis | |
| CN106290256A (en) | Quantitative background schlieren method based on video measuring | |
| CN111508011A (en) | Depth data calibration method of flight time camera | |
| CN109631751A (en) | A kind of output of high-frequency without dead zone cold atom interferometer | |
| CN107850423A (en) | For measurement apparatus, system and the manufacture method of the shape for measuring destination object | |
| CN1327240C (en) | Pulse echo treating method and device for laser range finder | |
| CN109708674A (en) | A kind of cold atom cloud control method that quickly sunken prisoner and high frequency are cast in opposite directions | |
| CN106017864B (en) | Testing device and testing method for characteristic parameters of swing mirror | |
| CN110926373A (en) | Structured light plane calibration method and system under railway foreign matter detection scene | |
| CN105066983B (en) | The detection method and device of atomic group chilling temperature, flying speed and movement locus | |
| CN108243337A (en) | A kind of imaging sensor dynamic transfer function test phase alignment means | |
| CN105445570B (en) | Wide-angle lens camera signal-to-noise ratio testing system and measuring method |
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 |