CN106772134B - A kind of apparatus and method of automatic field compensation - Google Patents
A kind of apparatus and method of automatic field compensation Download PDFInfo
- Publication number
- CN106772134B CN106772134B CN201710118620.9A CN201710118620A CN106772134B CN 106772134 B CN106772134 B CN 106772134B CN 201710118620 A CN201710118620 A CN 201710118620A CN 106772134 B CN106772134 B CN 106772134B
- Authority
- CN
- China
- Prior art keywords
- axis
- coil
- gradient coil
- magnetic field
- order gradient
- 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/0017—Means for compensating offset magnetic fields or the magnetic flux to be measured; Means for generating calibration magnetic fields
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a kind of devices for automatic field compensation, the first magnetoresistive chip, the second magnetoresistive chip, third magnetoresistive chip and the 4th magnetoresistive chip including being used to measure magnetic field value in external environment, it further include fixed cell and support unit, further include the three-dimensional Helmholtz coil, three-dimensional First-order Gradient coil and second order gradient coil unit for shimming, further includes acquisition unit, feedback unit and rotary unit.Also disclose a kind of method in automatic compensation magnetic field.Field compensation range of the present invention is big, and structure is simple, easy to operate.It is applicable in application based on the physical field research by the detection of downfield environmental signal and atomic magnetic force, is of great significance.
Description
Technical field
The invention belongs to magnetic field shimming fields, are more particularly to a kind of device of automatic field compensation, further relate to one kind certainly
The method of moving field compensation, the related science suitable for downfield environment are studied.
Background technique
Magnetic field is objective reality, for example, the magnetic field of earth surface is about 0.5 Gauss, geographical location and environment are not
Together, its size can change with direction.
Now many scientific research needs carry out in downfield or nearly zero magnetic field environment, such as SERF (spin-
Exchange-relaxation-free) atomic magnetic force meter could effectively work when its probe is under 10nT environment, because
This, it needs to shield the magnetic field in external environment using materials such as multilayer permalloys.
The shielding cylinder of permalloy production limits the application range of atomic magnetic force meter in volume, weight, structure.Example
For the detection of underwater different magnetic, the detection in earth's magnetic field, the detection of deep airborne magnetic, this just needs one for such as satellite or aircraft
The extraneous useless magnetic field of kind active shield, the device for obtaining extraneous useful magnetic signal.
The method that some automatic compensation magnetic fields had been developed in recent years and measure magnetic field.For example, S.J.Seltzer et al.
In " Unshielded three-axis vector operation of a spin-exchange-relaxation-free
It proposes that magnetic field modulation can be passed through in atomic magnetometer " [Appl.Phys.Lett.85 (20), 4804 (2004)]
Mode compensate magnetic field so that central area reaches nearly null field.Jiancheng Fang et al. is in " In situ triaxial
magnetic field compensation for the spin-exchange-relaxation-free atomic
It proposes and can finally be detected by field sweep in magnetometer " [Rev.Sci.Instrum.83,103104 (2012)]
The mode of light intensity compensates magnetic field.
Although current field compensation mode is able to solve some problems of application aspect.But they still have it is certain
Limitation.For example, compensating magnetic field by way of modulation, this requires the magnetic field of magnetometer probe working environment to be less than 10nT.
Therefore, develop new technology and method to be still a problem to be solved, the present invention will provide automatic compensation magnetic in response to this problem
The apparatus and method of field.
Summary of the invention
It is an object of the invention in view of the problems of the existing technology, provide a kind of device of automatic field compensation,
A kind of method that also proposed automatic field compensation.Solve the problems, such as that field compensation range is small, also solution atomic magnetic force meter is small
The problem of type.
To achieve the goals above, apparatus of the present invention use following technical scheme:
A kind of device of automatic field compensation, including turntable further include the three-dimensional Helmholtz of setting on a spinstand
Coil and three-dimensional First-order Gradient coil are horizontally arranged with the lucite tube of tubular on turntable, are arranged on lucite tube
Z2 gradient coil, lucite tube is interior to be fixed with supporting table by support rod, and four positioned at rectangular surfaces are provided in supporting table
The first magnetoresistive chip, the second magnetoresistive chip, third magnetoresistive chip and the 4th magnetoresistive chip of apex angle, the folder of rectangular surfaces and horizontal plane
Angle is 45 degree, and three-dimensional Helmholtz coil includes X-axis Helmholtz coil, Y-axis Helmholtz coil and Z axis Helmholtz's line
Circle, three-dimensional First-order Gradient coil includes X-axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil, Z2 gradient
The central axis of coil, lucite tube central axis, the central axis of Z axis Helmholtz coil, Z axis First-order Gradient coil
Central axis is overlapped, Z2 gradient coil central point, three-dimensional Helmholtz coil central point, three-dimensional First-order Gradient hub of a spool point,
Lucite tube central point, rectangular surfaces central point are overlapped.
It further include for being exported to the first magnetoresistive chip, the second magnetoresistive chip, third magnetoresistive chip and the 4th magnetoresistive chip
Signal carry out the differential amplifier circuit of differential amplification, differential amplifier circuit is connected by the first multi-way contral card and computer
It connects, computer is connect with the second multi-way contral card, the second multi-way contral card and X-axis First-order Gradient coil, Y-axis First-order Gradient
Coil, Z axis First-order Gradient coil, Z2 gradient coil connection, the second multi-way contral card also pass through power amplifier respectively with the X-axis last of the twelve Earthly Branches
Mu Huozi coil, Y-axis Helmholtz coil are connected with Z axis Helmholtz coil.
It is provided with rectangle on turntable as described above and places block, rectangle is placed and offers bar shaped placing groove, bar shaped on block
The cross section of placing groove is up-side down triangle, and lucite tube is placed in bar shaped placing groove.
X-axis Helmholtz coil, Y-axis Helmholtz coil, Z axis Helmholtz coil include concentrically as described above
A pair of of coil of axis, the central point coincidence of X-axis Helmholtz coil, Y-axis Helmholtz coil, Z axis Helmholtz coil, X
Axis Helmholtz coil, Y-axis Helmholtz coil, Z axis Helmholtz coil central axis be mutually perpendicular to, X-axis Helmholtz
Coil, Y-axis Helmholtz coil, the outer diameter of Z axis Helmholtz coil are equal, X-axis Helmholtz coil, Y-axis Helmholtz's line
It is fixed to each other between circle, Z axis Helmholtz coil by acrylic plastic stent and plastic screw.
X-axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil include concentrically as described above
A pair of of coil of axis, the central point coincidence of X-axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil, X
Axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil central axis be mutually perpendicular to, X-axis First-order Gradient
Coil, Y-axis First-order Gradient coil, the outer diameter of Z axis First-order Gradient coil are equal, X-axis First-order Gradient coil, Y-axis First-order Gradient line
It encloses, the distance D of each pair of coil of Z axis First-order Gradient coil and outside diameter d ratio are 0.886, X-axis First-order Gradient coil, Y-axis one
It is fixed between rank gradient coil, Z axis First-order Gradient coil by acrylic plastic stent and plastic screw.
Z2 gradient coil as described above include first annular coil to and the second loop coil pair, first annular coil pair
Two loop coils distance and outer diameter ratio be 1.2, the distance and outer diameter of two loop coils of the second loop coil pair
Ratio is 0.3, first annular coil to and the central axis of the second loop coil pair, lucite tube central axis be overlapped.
First annular coil as described above to and the second loop coil to the wire casing for being arranged at organic glass pipe outer wall
In.
A kind of method of automatic field compensation, comprising the following steps:
Step 1, for X-axis Helmholtz coil, Y-axis Helmholtz coil, Z axis Helmholtz coil, successively load is set respectively
Fixed each given voltage tests X-axis Helmholtz coil, Y-axis Helmholtz coil, Z axis using the first magnetoresistive chip respectively
Helmholtz coil is worth corresponding magnetic field value size in each given voltage, is fitted X-axis Helmholtz coil, Y-axis last of the twelve Earthly Branches nurse respectively
Suddenly hereby coil, the given voltage of Z axis Helmholtz coil and magnetic field ratio relation；
The each of setting is successively loaded respectively for X-axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil
It is a to give constant current, X-axis First-order Gradient coil, Y-axis First-order Gradient coil, one ladder of Z axis are tested respectively using the first magnetoresistive chip
Spend coil in the corresponding magnetic field value size of each given current value, respectively linear fit X-axis First-order Gradient coil, one ladder of Y-axis
Spend the ratio relation to constant current and magnetic field of coil, Z axis First-order Gradient coil；
It successively loads for Z2 gradient coil to constant current, using the test of the first magnetoresistive chip under each given current condition
Corresponding magnetic field value size, the ratio relation to constant current and magnetic field of linear fit Z2 gradient coil；
Step 2, the rotation speed for setting turntable, X-axis First-order Gradient coil, Y-axis First-order Gradient coil, one ladder of Z axis
Spend coil, X-axis Helmholtz coil, Y-axis Helmholtz coil, Z axis Helmholtz coil, Z2 gradient coil, the first magnetic resistance core
Piece, the second magnetoresistive chip, third magnetoresistive chip, the 4th magnetoresistive chip are as turntable rotates in the horizontal plane；
Step 3, initialized target magnetic field value B0Matrix, magnetic field of the goal value B0The matrix that matrix is one 1 × 12；
Step 4 measures magnetic field using the first magnetoresistive chip, the second magnetoresistive chip, third magnetoresistive chip, the 4th magnetoresistive chip
Value Bn, the number of n representative measurement, acquisition measurement magnetic field value BnMatrix measures magnetic field value BnMatrix is 1 × 12 matrix, measurement
Magnetic field value BnElement in matrix respectively correspond from top to bottom the first magnetoresistive chip measure X axis, Y-axis, Z axis to magnetic field
Value, the second magnetoresistive chip measure X axis, Y-axis, Z axis to magnetic field value, third magnetoresistive chip measures X axis, Y-axis, Z axis
To magnetic field value and the 4th magnetoresistive chip measure X axis, Y-axis, Z axis to magnetic field value；
Step 5 judges magnetic field of the goal value B0Matrix subtracts measurement magnetic field value BnMatrix of differences e after matrixnModulus value whether
In error range, if in error range, step 8 is gone to；If not going to step 6 in error range；
Step 6, by matrix of differences enIt carries out PID arithmetic and obtains output matrix Pn, output matrix PnFor 1 × 12 matrix；
Step 7, by output matrix Pn, the given voltage of X-axis Helmholtz coil and ratio relation, the Y-axis last of the twelve Earthly Branches nurse in magnetic field
Suddenly hereby the ratio relation in the given voltage of coil and magnetic field, the given voltage of Z axis Helmholtz coil and magnetic field ratio relation,
X-axis First-order Gradient coil to the ratio relation in constant current and magnetic field, the ratio to constant current and magnetic field of Y-axis First-order Gradient coil
Value relationship, Z axis First-order Gradient coil to the ratio relation in constant current and magnetic field, Z2 gradient coil to constant current and magnetic field
Ratio relation is input in the linprog function in Matlab, uses the L1 norm programming evaluation min ∑ based on linear programming
Vi×wi+Ij×NjOptimization problem is minimized, optimization matrix is obtained, first three element optimized in matrix corresponds respectively to the X-axis last of the twelve Earthly Branches
Mu Huozi coil, Y-axis Helmholtz coil, the Optimized loading voltage in Z axis Helmholtz coil, optimize latter four in matrix
Element correspond respectively to X-axis First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil, in Z2 gradient coil
Optimized loading electric current, goes to step 4, wherein i ∈ { 1,2,3 }, j ∈ { 1,2,3,4 }, w1、w2、w3Respectively indicate X-axis Hai Muhuo
The hereby voltage weight of coil, Y-axis Helmholtz coil, Z axis Helmholtz coil；V1、V2、V3Respectively indicate X-axis Helmholtz's line
The on-load voltage of circle, Y-axis Helmholtz coil, Z axis Helmholtz coil；N1、N2、N3、N4Respectively indicate X-axis First-order Gradient line
The current weights of circle, Y-axis First-order Gradient coil, Z axis First-order Gradient coil, Z2 gradient coil, I1、I2、I3、I4Respectively indicate X-axis
First-order Gradient coil, Y-axis First-order Gradient coil, Z axis First-order Gradient coil, Z2 gradient coil loading current, linprog function
Constraint condition are as follows: 10-3V≤|Vi|≤10V；10-6A≤|Ij|≤2×10-2A；Magnetic field error △ B=| Bn-B0|≤10nT；
Uniformity of magnetic field is less than or equal to 10-3。
The present invention has the advantages that compared with the existing technology
Field compensation range of the present invention is big, and structure is simple, easy to operate.It is applicable in the object based on the detection of downfield environmental signal
Area research and atomic magnetic force meter application are managed, is of great significance.
Detailed description of the invention
Fig. 1 is general structure schematic diagram of the invention；
The program flow diagram of Fig. 2 field compensation of the present invention；
Fig. 3 is the main view that rectangle places block and lucite tube coupling part；
Fig. 4 is the left view that rectangle places block and lucite tube coupling part；
Fig. 5 is the right view that rectangle places block and lucite tube coupling part；
In figure: 1-X axis First-order Gradient coil；2-Z axis First-order Gradient coil；3-Y axis First-order Gradient coil；4-Z axis last of the twelve Earthly Branches nurse
Hereby coil suddenly；5-X axis Helmholtz coil；6-Y axis Helmholtz coil；The first magnetoresistive chip of 7-；8-Z2 gradient coil；9-
Two magnetoresistive chips；10- supporting table；11- rectangle places block；The first fixed link of 12-；The second fixed link of 13-；14- third fixed link；
15- turntable；16- computer；The first data line of 17-；18- the first multi-way contral card；The second data line of 19-；20- difference is put
Big circuit；21- third data line；The 4th data line of 22-；The 5th data line of 23-；24- power amplifier；The 6th data line of 25-；
26- the second multi-way contral card；The 7th data line of 27-；The first pcb board of 28-；The second pcb board of 29-；30- third magnetoresistive chip；
The 4th magnetoresistive chip of 31-；32- lucite tube；33- steady arm；The 4th fixed link of 34-；The 5th fixed link of 35-；36- the 6th is solid
Fixed pole；The 7th fixed link of 37-；The 8th fixed link of 38-；39- terminal box.
Specific embodiment
Present invention is further described in detail with reference to the accompanying drawing:
Embodiment 1:
A kind of device of automatic field compensation, including turntable 15 further include the three-dimensional last of the twelve Earthly Branches nurse being arranged on turntable 15
Hereby coil and three-dimensional First-order Gradient coil suddenly are horizontally arranged with the lucite tube 32 of tubular, lucite tube on turntable 15
Z2 gradient coil 8 is arranged on 32, lucite tube 32 is interior to be fixed with supporting table 10 by support rod, is provided in supporting table 10
Positioned at the first magnetoresistive chip 7 of four apex angles of rectangular surfaces, the second magnetoresistive chip 9, third magnetoresistive chip 30 and the 4th magnetic resistance core
The angle of piece 31, rectangular surfaces and horizontal plane is 45 degree, and three-dimensional Helmholtz coil includes X-axis Helmholtz coil 5, Y-axis last of the twelve Earthly Branches nurse
Hereby coil 6 and Z axis Helmholtz coil 4, three-dimensional First-order Gradient coil include X-axis First-order Gradient coil 1, Y-axis First-order Gradient suddenly
Coil 3, Z axis First-order Gradient coil 2, central axis, 32 central axis of lucite tube, Z axis Helmholtz of Z2 gradient coil 8
The central axis of coil 4, the central axis of Z axis First-order Gradient coil 2 are overlapped, 8 central point of Z2 gradient coil, three-dimensional Helmholtz
Hub of a spool point, three-dimensional First-order Gradient hub of a spool point, 32 central point of lucite tube, rectangular surfaces central point are overlapped.
A kind of device of automatic field compensation further includes for the first magnetoresistive chip 7, the second magnetoresistive chip 9, third magnetic
Hinder the differential amplifier circuit 20 of the signal progress differential amplification of chip 30 and the output of the 4th magnetoresistive chip 31, differential amplifier circuit 20
It is connect by the first multi-way contral card 18 with computer 16, computer 16 is connect with the second multi-way contral card 26, more than second
Channel control card 26 and X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, Z2 gradient coil 8
Connection, the second multi-way contral card 26 also by power amplifier 24 respectively with X-axis Helmholtz coil 5, Y-axis Helmholtz coil
6 and Z axis Helmholtz coil 4 connect.X-axis, Y-axis, Z axis are mutually perpendicular to, and plane where X-axis, Z axis is horizontal plane.
It is provided with rectangle on turntable 15 and places block 11, rectangle is placed and offers bar shaped placing groove on block 11, and bar shaped is placed
The cross section of slot is up-side down triangle, and lucite tube 32 is placed in bar shaped placing groove.
X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4 are including concentrically axis
A pair of of coil, the central point coincidence of X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4, the X-axis last of the twelve Earthly Branches
Mu Huozi coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4 central axis be mutually perpendicular to, X-axis Helmholtz's line
Circle 5, Y-axis Helmholtz coil 6, the outer diameter of Z axis Helmholtz coil 4 are equal, X-axis Helmholtz coil 5, Y-axis Helmholtz
It is fixed to each other between coil 6, Z axis Helmholtz coil 4 by acrylic plastic stent and plastic screw.
X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2 are including concentrically axis
A pair of of coil, the central point coincidence of X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, X-axis one
Rank gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2 central axis be mutually perpendicular to, X-axis First-order Gradient line
Circle 1, Y-axis First-order Gradient coil 3, the outer diameter of Z axis First-order Gradient coil 2 are equal, X-axis First-order Gradient coil 1, Y-axis First-order Gradient
Coil 3, Z axis First-order Gradient coil 2 each pair of coil distance D and outside diameter d ratio be 0.886, X-axis First-order Gradient coil 1,
It is fixed between Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2 by acrylic plastic stent and plastic screw.
Z2 gradient coil 8 include first annular coil to and the second loop coil pair, two rings of first annular coil pair
The distance and outer diameter ratio of shape coil are 1.2, and the distance and outer diameter ratio of two loop coils of the second loop coil pair are
0.3, first annular coil to and the central axis of the second loop coil pair, 32 central axis of lucite tube be overlapped.
First annular coil to and the second loop coil to being arranged in the wire casing of 32 outer wall of lucite tube.
Embodiment 2:
A kind of device for automatic field compensation, including fixed cell, support unit, three-dimensional Helmholtz coil, three
Tie up First-order Gradient coil, second order gradient coil unit, acquisition unit, feedback unit, rotary unit.
Fixed cell includes that rectangle places block 11, the first fixed link 12, the second fixed link 13, third fixed link the 14, the 4th
Fixed link 34, the 5th fixed link 35, the 6th fixed link 36, the 7th fixed link 37, the 8th fixed link 38, it is whole that rectangle places block 11
It for rectangle, is fixed on support rod 33, the top surface that rectangle places block 11 is provided with the placing groove of strip, the cross section of placing groove
For up-side down triangle, lucite tube 32 is placed on the placing groove that rectangle is placed on block 11, and diameter is 100mm on turntable 15
The first fixed link 12, the second fixed link 13, the 4th fixed link 34 and the 5th fixed link 35, Y-axis Hai Muhuo are symmetric on circle
Hereby coil 6 is supported on fixed turntable 15 by the second fixed link 13 and the 4th fixed link 34, and lucite tube 32 passes through the
One fixed link 12 and the 5th fixed link 35 are supported on fixed turntable 15, and diameter is on 300mm circle in phase on turntable 15
Third fixed link 14, the 6th fixed link 36, the 7th fixed link 37, the 8th fixed link 38 with symmetric mode distribution, third are fixed
Bar 14, the 6th fixed link 36, the 7th fixed link 37, the 8th fixed link 38 bottom end connect with turntable 15, third fixed link 14,
6th fixed link 36, the 7th fixed link 37, the 8th fixed link 38 top connect with Y-axis First-order Gradient coil 3.
Support unit includes the supporting table 10 and cylindrical support for placing the first pcb board 28, the second pcb board 29
Bar 33, the length × width × height of supporting table 10 are 40 × 40 × 10mm3Rectangular block, the central point of 10 upper surface of supporting table is located at three
The place 10mm immediately below Helmholtz coil central point is tieed up, the first pcb board 28 and the second pcb board 29 are consolidated by using plastic screw
10 both ends of supporting table are scheduled on, the both ends of cylindrical support bar 33 are connected through a screw thread supporting table 10 and turntable 15 respectively.
Signal acquisition unit includes for surveying magnetic field the first magnetoresistive chip 7 in environment, the second magnetoresistive chip 9, third magnetic resistance
Chip 30, the 4th magnetoresistive chip 31 are used for transmission the first magnetoresistive chip 7, the second magnetoresistive chip 9, third magnetoresistive chip the 30, the 4th
The first pcb board 28, the second pcb board 29 for the magnetic signal that magnetoresistive chip 31 exports, the first magnetoresistive chip 7 and the second magnetoresistive chip 9
It is welded on the first pcb board 28, third magnetoresistive chip 30 and the 4th magnetoresistive chip 31 are welded on the second pcb board 29, the first magnetic
Chip 7, the second magnetoresistive chip 9, third magnetoresistive chip 30 and the 4th magnetoresistive chip 31 are hindered in distributed rectangular and place plane and water
Plane be in 45 degree of angles, the first pcb board 28, the second pcb board 29 length × it is wide be 40 × 40mm, the central point of the first pcb board 28
It is overlapped with the line of 29 central point of the second pcb board with the axis of Z axis Helmholtz coil 4, the central point of line and three-dimensional last of the twelve Earthly Branches nurse
Hereby the central point of coil is overlapped suddenly.
Signal acquisition unit further includes differential amplifier circuit 20 and the first multi-way contral truck 18, the first magnetoresistive chip 7,
Second magnetoresistive chip 9, third magnetoresistive chip 30, the 4th magnetoresistive chip 31 acquisition magnetic signal pass through respectively third data line 21 with
Differential amplifier circuit 20 connects, and differential amplifier circuit 20 is connect by the second data line 19 with the first multi-way contral truck 18,
First multi-way contral truck 18 is connect by the first data line 17 with computer 16.
Signal feedback unit further includes the second multi-way contral card 26, power amplifier 24, and power amplifier 24 is by the second multi-pass
The output voltage of road control card 26 acts on X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4,
X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, Z2 gradient coil 8 pass through the 4th number respectively
It is connect according to line 22 with the second multi-way contral card 26, X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz
Coil 4 is connect by the 5th data line 23 with power amplifier 24 respectively, and power amplifier 24 passes through the 6th data line more than 25 and second
Channel control card 26 connects, and the second multi-way contral card 26 is connect by the 7th data line 27 with computer 16.
X-axis, Y-axis, Z axis are mutually perpendicular to, and plane where X-axis, Z axis is horizontal plane.
Three-dimensional Helmholtz coil includes X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil
4, X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4 include a pair of of line of concentrically axis
Circle, the central point coincidence of X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4, X-axis Helmholtz
Coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4 central axis be mutually perpendicular to, X-axis Helmholtz coil 5, Y
Axis Helmholtz coil 6, the outer diameter of Z axis Helmholtz coil 4 are equal, X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6,
It is fixed by acrylic plastic stent with plastic screw between Z axis Helmholtz coil 4.
Three-dimensional First-order Gradient coil includes X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil
2, X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2 include a pair of of line of concentrically axis
Circle, the central point coincidence of X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, X-axis First-order Gradient
Coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2 central axis be mutually perpendicular to, X-axis First-order Gradient coil 1, Y
Axis First-order Gradient coil 3, the outer diameter of Z axis First-order Gradient coil 2 are equal, and the outside diameter d of distance D and coil between each pair of coil meet
D=0.886d relationship, X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3 pass through sub- gram between Z axis First-order Gradient coil 2
Power plastic stent is fixed with plastic screw.The center of the central point of three-dimensional First-order Gradient coil and three-dimensional Helmholtz coil
Point is overlapped.
Second order gradient coil unit includes the Z2 gradient coil 8 for generating second order gradient magnetic, and Z2 gradient coil includes
First annular coil to and the second loop coil pair, the distance D1 and loop coil of two loop coils of first annular coil pair
Outside diameter d 1 meet D1=1.2d1 relationship, the distance D2 of the second loop coil pair and the outside diameter d 2 of loop coil meet D2=
0.3d2 relationship, first annular coil to and the second loop coil pair central axis, 32 central axis of lucite tube, the Z axis last of the twelve Earthly Branches
The central axis of Mu Huozi coil 4, Z axis First-order Gradient coil 2 central axis be overlapped, it is the central point of Z2 gradient coil 8, organic
The central point coincidence of 32 central point of glass tube, the central point of three-dimensional Helmholtz coil, three-dimensional First-order Gradient coil, it is first annular
Coil to and the second loop coil in the coil fixing groove being each attached on lucite tube 32.
Rotary unit includes the turntable 15 without magnetic for one-dimensional rotation, and for connecting 21 signal of third data line
Output end, connection 22 signal input part of the 4th data line, the terminal box 39 for connecting 23 signal input part of the 5th data line, turntable
It is connected between 15 and terminal box 39 by bearing.
First magnetoresistive chip 7, the second magnetoresistive chip 9, third magnetoresistive chip 30, the 4th magnetoresistive chip 31 measure external environment
In magnetic field, magnetic signal is changed into voltage signal；Differential amplifier circuit 20 amplifies voltage signal, the voltage letter of amplification
It number is obtained by the first multi-way contral card 18.Computer determines voltage and current value in every group of coil；Voltage and electric current pass through the
The output of two multi-way contral cards 26；The voltage value of output makees Z axis Helmholtz coil 4, X-axis Helmholtz by power amplifier 24
Coil 5, Y-axis Helmholtz coil 6, the current value of output directly act on X-axis First-order Gradient coil 1, Z axis First-order Gradient coil
2, Y-axis First-order Gradient coil 3,8 coil of Z2 gradient coil.
Component is described as follows in apparatus of the present invention:
Z axis Helmholtz coil 4.It is made of Z axis Helmholtz coil 4 two copper wire wound coils, when work provides Z axis
Magnetic field.The enameled wire that model AIW/QZY+XY-2/220 outer diameter is 0.67mm is constituted.
X-axis Helmholtz coil 5.Structure identical as Z axis Helmholtz coil 4, identical size, identical enameled wire coiling,
X-axis magnetic field is provided when work.
Y-axis Helmholtz coil 6.Structure identical as Z axis Helmholtz coil 4, identical size, identical enameled wire coiling,
Y-axis magnetic field is provided when work.Three-dimensional Helmholtz's line is collectively formed with Z axis Helmholtz coil 4 and X-axis Helmholtz coil 5
Circle, apparatus of the present invention typically use: the central point of each pair of coil is overlapped, and the center line of three pairs of coils is mutually perpendicular to, each pair of line
The outer diameter of circle is equal.
Z axis First-order Gradient coil 2.It is made of Z axis First-order Gradient coil 2 two copper wire wound coils, when work provides Z axis
First-order Gradient magnetic field.The enameled wire that model AIW/QZY+XY-2/220 outer diameter is 0.67mm is constituted.
X-axis First-order Gradient coil 1.Structure identical as Z axis First-order Gradient coil 2, identical size, identical enameled wire coiling,
X-axis First-order Gradient magnetic field is provided when work.
Y-axis First-order Gradient coil 3.Structure identical as Z axis First-order Gradient coil 2, identical size, identical enameled wire coiling,
Y-axis First-order Gradient magnetic field is provided when work.Constitute three-dimensional one together with Z axis First-order Gradient coil 2 and X-axis First-order Gradient coil 1
Rank gradient coil, apparatus of the present invention typically use: the central point of each pair of coil is overlapped, and the center line of three pairs of coils mutually hangs down
Directly, the outer diameter of each pair of coil is equal.
Z2 gradient coil.The enameled wire that model AIW/QZY+XY-2/220 outer diameter is 0.33mm is constituted, for compensating ring
Z2 gradient magnetic in border.
First magnetoresistive chip.Model HMC1053, for measuring the magnetic field in environment.
Second magnetoresistive chip.It is identical as the first magnetoresistive chip model, effect.
Third magnetoresistive chip.It is identical as the first magnetoresistive chip model, effect.
4th magnetoresistive chip.It is identical as the first magnetoresistive chip, the second magnetoresistive chip model, effect.With the first magnetoresistive chip,
Second magnetoresistive chip, third magnetoresistive chip together, are placed in quadrangle.
Supporting table.Material is aluminium material, rectangle, for installing, placing the first magnetoresistive chip, the second magnetoresistive chip, third
Magnetoresistive chip, the 4th magnetoresistive chip.
Support rod.PEEK material, cylindrical, there is screw thread at both ends.
Rectangle places block.PEEK material, for installing, fixing the Z2 gradient on branch lucite tube and lucite tube
Coil.
First fixed link.PEEK material, for fixing Y-axis Helmholtz coil, fixed lucite tube and organic glass
The Z2 gradient coil installed, placed above glass pipe.
Second fixed link.It is identical as the first fixed link material, and the 4th fixed link collective effect, for fixing Y-axis last of the twelve Earthly Branches nurse
Hereby coil suddenly.
Third fixed link.It is identical as the first fixed link material, it is total with the 6th fixed link, the 7th fixed link, the 8th fixed link
Same-action, for fixing Y-axis First-order Gradient coil.
4th fixed link.It is identical as the first fixed link material, and the second fixed link collective effect, for fixing Y-axis last of the twelve Earthly Branches nurse
Hereby coil suddenly.
5th fixed link.It is identical as the first fixed link material, and the first fixed link collective effect, for fixing Y-axis last of the twelve Earthly Branches nurse
Hereby coil suddenly, installation, the Z2 gradient coil placed above fixed lucite tube and lucite tube.
6th fixed link.It is identical as the first fixed link material, it is total with third fixed link, the 7th fixed link, the 8th fixed link
Same-action, for fixing Y-axis single order coil.
7th fixed link.It is identical as the first fixed link material, it is total with third fixed link, the 6th fixed link, the 8th fixed link
Same-action, for fixing Y-axis single order coil.
8th fixed link.It is identical as the first fixed link material, it is total with third fixed link, the 6th fixed link, the 7th fixed link
Same-action, for fixing Y-axis First-order Gradient coil.
Turntable without magnetic.Material is aluminium and stainless steel, nonmagnetic, one-dimensional precise electric controlled rotating.
Terminal box.It is connected for line.
Computer.Model apocalypse M7150.
First data line.Computer-internal slot is used for transmission acquisition data.
First multi-way contral card.Model PCIe6251 is used for data acquisition control.
Second data line.Model SHC68-68-EPM Cable, is used for transmission acquisition data.
Differential amplifier circuit.It is made of the difference amplifier for amplifying voltage signal.
Third data line.Model BVR-1.5, transmission acquisition data.
4th data line.Model is identical as third data line.
Power amplifier.Power amplification circuit, function are that voltage remains unchanged Current amplifier.
5th data line.Model is identical as third data line.
6th data line.Model SH68-68-D1Cable is used for transmission feedback voltage data.
Second multi-way contral card.Model PCI6704 is exported for feedback voltage and electric current.
7th data line.Model is identical as the first data line.
Embodiment 3:
A method of it is carried out using embodiment 1 or embodiment 2 for automatic field compensation, comprising the following steps:
Step 1,
Successively loaded for X-axis Helmholtz coil 5 given voltage -2V, -1.5V, -1V, -0.5V, 0V, 0.5V, 1V,
1.5V, 2V are worth corresponding magnetic field value size, linear fit X-axis last of the twelve Earthly Branches nurse in each given voltage using the test of the first magnetoresistive chip 7
The hereby ratio relation of the given voltage of coil 5 and magnetic field suddenly；
Successively loaded for Y-axis Helmholtz coil 6 given voltage -2V, -1.5V, -1V, -0.5V, 0V, 0.5V, 1V,
1.5V, 2V are worth corresponding magnetic field value size, linear fit Y-axis last of the twelve Earthly Branches nurse in each given voltage using the test of the first magnetoresistive chip 7
The hereby ratio relation of the given voltage of coil 6 and magnetic field suddenly；
Successively loaded for Z axis Helmholtz coil 4 given voltage -2V, -1.5V, -1V, -0.5V, 0V, 0.5V, 1V,
1.5V, 2V are worth corresponding magnetic field value size, linear fit Z axis last of the twelve Earthly Branches nurse in each given voltage using the test of the first magnetoresistive chip 7
The hereby ratio relation of the given voltage of coil 4 and magnetic field suddenly；
For X-axis First-order Gradient coil 1 successively load to constant current -0.05A, -0.04A, -0.03A, 0A, 0.01A,
0.02A, 0.03A, 0.04A, 0.05A are big in the corresponding magnetic field value of each given current value using the test of the first magnetoresistive chip 7
It is small, the ratio relation to constant current and magnetic field of linear fit X-axis First-order Gradient coil 1；
For Y-axis First-order Gradient coil 3 successively load to constant current -0.05A, -0.04A, -0.03A, 0A, 0.01A,
0.02A, 0.03A, 0.04A, 0.05A are big in the corresponding magnetic field value of each given current value using the test of the first magnetoresistive chip 7
It is small, the ratio relation to constant current and magnetic field of linear fit Y-axis First-order Gradient coil 3；
For Z axis First-order Gradient coil 2 successively load to constant current -0.05A, -0.04A, -0.03A, 0A, 0.01A,
0.02A, 0.03A, 0.04A, 0.05A are big in the corresponding magnetic field value of each given current value using the test of the first magnetoresistive chip 7
It is small, the ratio relation to constant current and magnetic field of linear fit Z axis First-order Gradient coil 2；
For Z2 gradient coil 8 successively load to constant current -0.005A, -0.004A, -0.003A, 0A, 0.001A,
0.002A, 0.003A, 0.004A, 0.005A use the first magnetoresistive chip 7 test corresponding magnetic under each given current condition
Field size, the ratio relation to constant current and magnetic field of linear fit Z2 gradient coil 8.
The electric current of first annular coil pair is the setting times in the electric current of the second loop coil centering in Z2 gradient coil 8
Number, set multiple as 3.3 times, is that 8 loading current I'(of Z2 gradient coil shows the second loop coil to constant current I') to adding
Electric current I'(is carried to constant current I'), first annular coil 3.3 times of electric current I' given to the electric current I'(for loading 3.3 times).First
The distance and outer diameter ratio of two loop coils of loop coil pair are 1.2, two loop coils of the second loop coil pair
Distance and outer diameter ratio are 0.3.
Step 2, rotation speed of the setting without magnetic rotation turntable 15, X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z
Axis First-order Gradient coil 2, X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4, Z2 gradient coil
8, the first magnetoresistive chip 7, the second magnetoresistive chip 9, third magnetoresistive chip 30, the 4th magnetoresistive chip 31 are with no magnetic rotation turntable 15
It rotates in the horizontal plane.
Step 3, initialized target magnetic field value B0Matrix, magnetic field of the goal value B0The matrix that matrix is one 1 × 12, target magnetic
Field value B0Each element value is 10nT in matrix.
Step 4 is surveyed using the first magnetoresistive chip 7, the second magnetoresistive chip 9, third magnetoresistive chip 30, the 4th magnetoresistive chip 31
Obtain magnetic field value Bn(magnetic field value BnThe magnetic field value obtained for n-th measurement), obtain measurement magnetic field value BnMatrix measures magnetic field value Bn
Matrix is 1 × 12 matrix, measures magnetic field value BnElement in matrix respectively corresponds the first magnetoresistive chip 7 from top to bottom and measures
X axis, Y-axis, Z axis to magnetic field value, the second magnetoresistive chip 9 measure X axis, Y-axis, Z axis to magnetic field value, third magnetic resistance
Chip 30 measure X axis, Y-axis, Z axis to magnetic field value, the 4th magnetoresistive chip 31 measures X axis, Y-axis, Z axis to magnetic field
Value.
Step 5 judges magnetic field of the goal value B0Matrix subtracts measurement magnetic field value BnThe modulus value e of the matrix of differences of matrixnWhether
In error range, if going to step 8 in error range；If not going to step 6 in error range.
Step 6, by magnetic field of the goal value B0Matrix subtracts magnetic field value BnThe matrix of differences e of matrixnIt is defeated to carry out PID arithmetic acquisition
Matrix P outn, can be by matrix of differences enThe pid algorithm being input in computer 16, pid algorithm can be by calling in LabVIEW
PID module realizes that pid algorithm realizes data processing,Wherein en=
Bn-B0, en-1=Bn-1-B0(wherein Bn-1Obtained magnetic field value matrix is sampled for (n-1)th time), pass through the first magnetoresistive chip 7, second
Magnetoresistive chip 9, third magnetoresistive chip 30, the 4th magnetoresistive chip 31 measure magnetic field value to adjust Proportional coefficient KP, integral coefficient KI、
Differential coefficient KD, scaling up COEFFICIENT KPIt can accelerate the response of system, but excessive KP, biggish overshoot can be generated, and generate
Concussion, makes stability degenerate, and increases integral coefficient KIOvershoot is advantageously reduced, keeps system more stable, while can also be eliminated
Static error increases differential coefficient KDThe response for being conducive to quickening system, increases stability, but system have to disturbance it is more sensitive
Response；Finally obtain output matrix Pn, output matrix PnFor 1 × 12 matrix.
Step 7, by output matrix Pn, the given voltage of X-axis Helmholtz coil 5 and ratio relation, the Y-axis last of the twelve Earthly Branches nurse in magnetic field
Hereby the ratio of the ratio relation in the given voltage of coil 6 and magnetic field, the given voltage of Z axis Helmholtz coil 4 and magnetic field closes suddenly
System, X-axis First-order Gradient coil 1 to the ratio relation in constant current and magnetic field, Y-axis First-order Gradient coil 3 give constant current and magnetic
Ratio relation, Z axis First-order Gradient coil 2 to the ratio relation in constant current and magnetic field, Z2 gradient coil 8 give constant current
It is input in the linprog function in Matlab with the ratio relation in magnetic field, using based on linear programming (linear
Programming, LP) L1 norm programming evaluation min ∑ Vi×wi+Ij×NjThis minimizes optimization problem, obtains optimization square
Battle array compensates magnetic field by the voltage or current value in optimization coil, first three element optimized in matrix corresponds respectively to X
Axis Helmholtz coil 5, Y-axis Helmholtz coil 6, the Optimized loading voltage in Z axis Helmholtz coil 4, optimize in matrix
Rear four elements correspond respectively to X-axis First-order Gradient coil 1, Y-axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, Z2 ladder
The Optimized loading electric current in coil 8 is spent, step 4 is gone to.Wherein, { 1,2,3 } i ∈, j ∈ { 1,2,3,4 }, w1、w2、w3Table respectively
Show the voltage weight of X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4, is taken in the present embodiment
1；V1、V2、V3Respectively indicate the load electricity of X-axis Helmholtz coil 5, Y-axis Helmholtz coil 6, Z axis Helmholtz coil 4
Pressure；N1、N2、N3、N4Respectively indicate X-axis First-order Gradient coil 1Y axis First-order Gradient coil 3, Z axis First-order Gradient coil 2, Z2 gradient
The current weights of coil 8 take 1 in the present embodiment.I1、I2、I3、I4Respectively indicate X-axis First-order Gradient coil 1, one ladder of Y-axis
Spend the loading current of coil 3, Z axis First-order Gradient coil 2, Z2 gradient coil 8.By calling in Matlab in experiment
Linprog function, the constraint condition of linprog function are as follows: 1,10-3V≤|Vi|≤10V, 2,10-6A≤|Ij|≤2×10-2A,
3, magnetic field error △ B=| Bn-B0|≤10nT, 4, uniformity of magnetic field be less than or equal to 10-3。
Step 8, whether halt system, if so, halt system, if it is not, then going to step 4.
Claims (1)
1. a kind of method of automatic field compensation, using the device of automatic field compensation, the device of automatic field compensation includes rotation
Turntable (15) further includes the three-dimensional Helmholtz coil being arranged on turntable (15) and three-dimensional First-order Gradient coil, turntable
(15) it is horizontally arranged with the lucite tube (32) of tubular on, is arranged with Z2 gradient coil (8) on lucite tube (32), it is organic
It is fixed with supporting table (10) in glass tube (32) by support rod, four apex angles positioned at rectangular surfaces is provided in supporting table (10)
The first magnetoresistive chip (7), the second magnetoresistive chip (9), third magnetoresistive chip (30) and the 4th magnetoresistive chip (31), rectangular surfaces with
The angle of horizontal plane is 45 degree, and three-dimensional Helmholtz coil includes X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6)
With Z axis Helmholtz coil (4), three-dimensional First-order Gradient coil includes X-axis First-order Gradient coil (1), Y-axis First-order Gradient coil
(3), Z axis First-order Gradient coil (2), central axis, lucite tube (32) central axis, Z axis last of the twelve Earthly Branches nurse of Z2 gradient coil (8)
The hereby central axis coincidence of the central axis of coil (4), Z axis First-order Gradient coil (2) suddenly, Z2 gradient coil (8) central point, three
Tie up Helmholtz coil central point, three-dimensional First-order Gradient hub of a spool point, lucite tube (32) central point, rectangular surfaces central point
It is overlapped, which comprises the following steps:
Step 1 is distinguished successively for X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6), Z axis Helmholtz coil (4)
The each given voltage for loading setting, tests X-axis Helmholtz coil (5), Y-axis last of the twelve Earthly Branches nurse using the first magnetoresistive chip (7) respectively
Hereby coil (6), Z axis Helmholtz coil (4) are fitted the X-axis last of the twelve Earthly Branches in the corresponding magnetic field value size of each given voltage value respectively suddenly
Mu Huozi coil (5), Y-axis Helmholtz coil (6), the given voltage of Z axis Helmholtz coil (4) and the ratio in magnetic field close
System；
For X-axis First-order Gradient coil (1), Y-axis First-order Gradient coil (3), Z axis First-order Gradient coil (2), successively load is set respectively
Fixed is each to constant current, tests X-axis First-order Gradient coil (1), Y-axis First-order Gradient line respectively using the first magnetoresistive chip (7)
(3), Z axis First-order Gradient coil (2) are enclosed in the corresponding magnetic field value size of each given current value, respectively linear fit X-axis single order
Gradient coil (1), Y-axis First-order Gradient coil (3), Z axis First-order Gradient coil (2) the ratio relation to constant current and magnetic field；
It successively loads for Z2 gradient coil (8) to constant current, is tested using the first magnetoresistive chip (7) in each given current condition
Under corresponding magnetic field value size, the ratio relation to constant current and magnetic field of linear fit Z2 gradient coil (8)；
Step 2, the rotation speed for setting turntable (15), X-axis First-order Gradient coil (1), Y-axis First-order Gradient coil (3), Z axis
First-order Gradient coil (2), X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6), Z axis Helmholtz coil (4), Z2 ladder
Spend coil (8), the first magnetoresistive chip (7), the second magnetoresistive chip (9), third magnetoresistive chip (30), the 4th magnetoresistive chip (31)
As turntable (15) rotates in the horizontal plane；
Step 3, initialized target magnetic field value B0Matrix, magnetic field of the goal value B0The matrix that matrix is one 1 × 12；
Step 4 uses the first magnetoresistive chip (7), the second magnetoresistive chip (9), third magnetoresistive chip (30), the 4th magnetoresistive chip
(31) magnetic field value B is measuredn, the number of n representative measurement, acquisition measurement magnetic field value BnMatrix measures magnetic field value BnMatrix is one 1
× 12 matrixes measure magnetic field value BnElement in matrix respectively corresponds the first magnetoresistive chip (7) from top to bottom and measures X axis, Y-axis
To, Z axis to magnetic field value, the second magnetoresistive chip (9) measure X axis, Y-axis, Z axis to magnetic field value, third magnetoresistive chip
(30) measure X axis, Y-axis, Z axis to magnetic field value and the 4th magnetoresistive chip (31) measure X axis, Y-axis, Z axis to magnetic
Field value；
Step 5 judges magnetic field of the goal value B0Matrix subtracts measurement magnetic field value BnMatrix of differences e after matrixnModulus value whether accidentally
In poor range, if going to step 8 in error range；If not going to step 6 in error range；
Step 6, by matrix of differences enIt carries out PID arithmetic and obtains output matrix Pn, output matrix PnFor 1 × 12 matrix；
Step 7, by output matrix Pn, the given voltage of X-axis Helmholtz coil (5) and ratio relation, the Y-axis Hai Muhuo in magnetic field
Hereby the ratio relation in the given voltage of coil (6) and magnetic field, the given voltage of Z axis Helmholtz coil (4) and magnetic field ratio
Relationship, X-axis First-order Gradient coil (1) to the ratio relation in constant current and magnetic field, the given electricity of Y-axis First-order Gradient coil (3)
Stream and the ratio relation in magnetic field, ratio relation, the Z2 gradient coil (8) to constant current and magnetic field of Z axis First-order Gradient coil (2)
The ratio relation to constant current and magnetic field, be input in the linprog function in Matlab, use the L1 based on linear programming
Norm programming evaluation min ∑ Vi×wi+Ij×NjOptimization problem is minimized, optimization matrix is obtained, optimizes first three member in matrix
Element corresponds respectively to X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6), the optimization in Z axis Helmholtz coil (4)
On-load voltage, rear four elements optimized in matrix correspond respectively to X-axis First-order Gradient coil (1), Y-axis First-order Gradient coil
(3), Z axis First-order Gradient coil (2), the Optimized loading electric current in Z2 gradient coil (8), go to step 4, wherein i ∈ 1,2,
3 }, { 1,2,3,4 } j ∈, w1、w2、w3Respectively indicate X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6), Z axis last of the twelve Earthly Branches nurse
The hereby voltage weight of coil (4) suddenly；V1、V2、V3Respectively indicate X-axis Helmholtz coil (5), Y-axis Helmholtz coil (6), Z
The on-load voltage of axis Helmholtz coil (4)；N1、N2、N3、N4Respectively indicate X-axis First-order Gradient coil (1), Y-axis First-order Gradient
The current weights of coil (3), Z axis First-order Gradient coil (2), Z2 gradient coil (8), I1、I2、I3、I4Respectively indicate X-axis single order
Gradient coil (1), Y-axis First-order Gradient coil (3), Z axis First-order Gradient coil (2), Z2 gradient coil (8) loading current,
The constraint condition of linprog function are as follows: 10-3V≤|Vi|≤10V；10-6A≤|Ij|≤2×10-2A；Magnetic field error △ B=| Bn-
B0|≤10nT；Uniformity of magnetic field is less than or equal to 10-3。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710118620.9A CN106772134B (en) | 2017-03-01 | 2017-03-01 | A kind of apparatus and method of automatic field compensation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710118620.9A CN106772134B (en) | 2017-03-01 | 2017-03-01 | A kind of apparatus and method of automatic field compensation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106772134A CN106772134A (en) | 2017-05-31 |
| CN106772134B true CN106772134B (en) | 2019-05-21 |
Family
ID=58960594
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710118620.9A Active CN106772134B (en) | 2017-03-01 | 2017-03-01 | A kind of apparatus and method of automatic field compensation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772134B (en) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950516B (en) * | 2017-03-24 | 2019-06-14 | 中国科学院上海微***与信息技术研究所 | A kind of weak swirl magnetic field measuring device and method |
| CN107677976B (en) * | 2017-09-26 | 2019-07-16 | 中国科学院武汉物理与数学研究所 | A kind of self-adaptive kernel magnetic resonance gradient preemphasis Waveform generating apparatus and method |
| CN108896945B (en) * | 2018-06-19 | 2020-11-13 | 山东航天电子技术研究所 | Sensitivity index calibration method of high-sensitivity atomic magnetometer |
| CN109358302B (en) * | 2018-09-25 | 2020-02-18 | 中国科学院武汉物理与数学研究所 | Passive magnetic shielding-free atomic magnetometer device and magnetism measuring method |
| CN111060858B (en) * | 2018-10-17 | 2021-12-10 | 北京自动化控制设备研究所 | Method for generating high-uniformity magnetic field and gradient composite magnetic field in magnetic shielding barrel |
| CN110568384B (en) * | 2019-08-27 | 2020-08-18 | 中国科学院武汉物理与数学研究所 | Active magnetic compensation method for ultra-sensitive atomic magnetometer |
| EP4130893A4 (en) * | 2020-03-31 | 2024-04-24 | JEOL Ltd. | Tri-axial magnetic field correction coil, physical package, physical package for optical lattice clock, physical package for atomic clock, physical package for atom interferometer, physical package for quantum information processing device, and physical package system |
| CN111596244B (en) * | 2020-05-18 | 2022-04-12 | 武汉中科牛津波谱技术有限公司 | Multichannel separation matrix type shimming coil of nuclear magnetic resonance spectrometer and design method |
| CN112345806B (en) * | 2020-10-27 | 2024-03-15 | 中国电力科学研究院有限公司 | Magneto-resistive element test fixture and test system |
| CN113189528B (en) * | 2021-06-30 | 2021-09-21 | 季华实验室 | Reading method, device, equipment and medium applied to miniaturized atomic magnetometer |
| CN114264984B (en) * | 2021-11-22 | 2024-04-12 | 上海科技大学 | Method and system for measuring weak magnetic field of zero-field optical pump atomic magnetometer |
| CN114770410B (en) * | 2022-05-24 | 2023-04-25 | 中国科学院地质与地球物理研究所 | Fluxgate magnetometer helmholtz coil installation frock |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101251584A (en) * | 2008-04-09 | 2008-08-27 | 武汉大学 | Three-axial magnetometer correcting method and three-axial magnetic gradient correcting method |
| CN101285690A (en) * | 2008-05-26 | 2008-10-15 | 浙江大学 | Optical fibre gyroscope magnetic field-temperature sensitivity test method |
| CN101341418A (en) * | 2005-12-20 | 2009-01-07 | 皇家飞利浦电子股份有限公司 | Magnetic resonance scanner with a longitudinal magnetic field gradient system |
| CN101422365A (en) * | 2008-12-16 | 2009-05-06 | 中国科学院物理研究所 | Equalizing pulse external magnetic field adjustment device and method in high-temperature SQUID application |
| CN101762792A (en) * | 2008-12-24 | 2010-06-30 | 北京瑞达思科技开发公司 | Integrated device of ultra-low magnetic space and magnetic field reproduction |
| CN101893693A (en) * | 2010-07-16 | 2010-11-24 | 中国科学院上海微***与信息技术研究所 | Magnetic-field dynamic compensation system and methods based on spatial correlation |
| CN102290188A (en) * | 2011-04-19 | 2011-12-21 | 中国科学院电工研究所 | Winding device for shimming superconducting coil and and winding method thereof |
| CN102579048A (en) * | 2012-02-21 | 2012-07-18 | 大连理工大学 | Method for controlling direction of rotation axis and rotation direction of space universal superposition rotating magnetic field |
| CN102653719A (en) * | 2012-04-28 | 2012-09-05 | 吉林大学 | Cell culture device capable of generating multiple magnetic fields and culture method |
| CN103076581A (en) * | 2011-10-25 | 2013-05-01 | 西门子公司 | Control apparatus with differential control in a magnetically coupled coil system for an amplifier |
| JP2013148406A (en) * | 2012-01-18 | 2013-08-01 | Alps Electric Co Ltd | Magnetic sensor |
| CN103245928A (en) * | 2013-05-23 | 2013-08-14 | 中国科学院上海微***与信息技术研究所 | Method and device for uniform magnetic field and one-order gradient magnetic field with adjustable directions |
| CN104134511A (en) * | 2014-07-22 | 2014-11-05 | 中国科学院电工研究所 | Zero magnetic space system device and adjustment method thereof |
| CN104237819A (en) * | 2013-06-09 | 2014-12-24 | 泰州市中益新型纺织科技有限公司 | Magnet adjusting and uniformity compensating device |
| CN104545875A (en) * | 2013-10-18 | 2015-04-29 | 中国科学院上海微***与信息技术研究所 | Full-tensor space gradient compensation method and system used for biomagnetic detection |
| CN204679618U (en) * | 2015-06-05 | 2015-09-30 | 中国地震局地球物理研究所 | A kind of three-dimensional orthogonal coil device of field homogeneity distribution |
| CN104983386A (en) * | 2015-05-21 | 2015-10-21 | 大连理工大学 | Linear polarization phase correction method for space universal rotating field azimuth errors |
| CN105044636A (en) * | 2014-04-16 | 2015-11-11 | 西门子公司 | Method, system and magnetic resonance apparatus for compensating for inhomogeneities in the magnetic field |
| CN105203088A (en) * | 2015-09-16 | 2015-12-30 | 中国电子科技集团公司第四十九研究所 | Thee-dimensional magnetic-induction magnetic compass |
| CN103926627B (en) * | 2014-04-23 | 2016-08-17 | 吉林大学 | Underwater carrier ground magnetic tricomponent measuring method |
| CN105929458A (en) * | 2016-03-21 | 2016-09-07 | 吉林大学 | Aeromagnetic field vector detecting device and detecting method |
| CN106341975A (en) * | 2016-09-14 | 2017-01-18 | 北京航空航天大学 | Hybrid magnetic shielding device based on high temperature superconducting coil |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2528123A (en) * | 2014-07-11 | 2016-01-13 | Univ Cape Town | Correcting for main magnetic field drift in MRI scanners |
-
2017
- 2017-03-01 CN CN201710118620.9A patent/CN106772134B/en active Active
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101341418A (en) * | 2005-12-20 | 2009-01-07 | 皇家飞利浦电子股份有限公司 | Magnetic resonance scanner with a longitudinal magnetic field gradient system |
| CN101251584A (en) * | 2008-04-09 | 2008-08-27 | 武汉大学 | Three-axial magnetometer correcting method and three-axial magnetic gradient correcting method |
| CN101285690A (en) * | 2008-05-26 | 2008-10-15 | 浙江大学 | Optical fibre gyroscope magnetic field-temperature sensitivity test method |
| CN101422365A (en) * | 2008-12-16 | 2009-05-06 | 中国科学院物理研究所 | Equalizing pulse external magnetic field adjustment device and method in high-temperature SQUID application |
| CN101762792A (en) * | 2008-12-24 | 2010-06-30 | 北京瑞达思科技开发公司 | Integrated device of ultra-low magnetic space and magnetic field reproduction |
| CN101893693A (en) * | 2010-07-16 | 2010-11-24 | 中国科学院上海微***与信息技术研究所 | Magnetic-field dynamic compensation system and methods based on spatial correlation |
| CN102290188A (en) * | 2011-04-19 | 2011-12-21 | 中国科学院电工研究所 | Winding device for shimming superconducting coil and and winding method thereof |
| CN103076581A (en) * | 2011-10-25 | 2013-05-01 | 西门子公司 | Control apparatus with differential control in a magnetically coupled coil system for an amplifier |
| JP2013148406A (en) * | 2012-01-18 | 2013-08-01 | Alps Electric Co Ltd | Magnetic sensor |
| CN102579048A (en) * | 2012-02-21 | 2012-07-18 | 大连理工大学 | Method for controlling direction of rotation axis and rotation direction of space universal superposition rotating magnetic field |
| CN102653719A (en) * | 2012-04-28 | 2012-09-05 | 吉林大学 | Cell culture device capable of generating multiple magnetic fields and culture method |
| CN103245928A (en) * | 2013-05-23 | 2013-08-14 | 中国科学院上海微***与信息技术研究所 | Method and device for uniform magnetic field and one-order gradient magnetic field with adjustable directions |
| CN104237819A (en) * | 2013-06-09 | 2014-12-24 | 泰州市中益新型纺织科技有限公司 | Magnet adjusting and uniformity compensating device |
| CN104545875A (en) * | 2013-10-18 | 2015-04-29 | 中国科学院上海微***与信息技术研究所 | Full-tensor space gradient compensation method and system used for biomagnetic detection |
| CN105044636A (en) * | 2014-04-16 | 2015-11-11 | 西门子公司 | Method, system and magnetic resonance apparatus for compensating for inhomogeneities in the magnetic field |
| CN103926627B (en) * | 2014-04-23 | 2016-08-17 | 吉林大学 | Underwater carrier ground magnetic tricomponent measuring method |
| CN104134511A (en) * | 2014-07-22 | 2014-11-05 | 中国科学院电工研究所 | Zero magnetic space system device and adjustment method thereof |
| CN104983386A (en) * | 2015-05-21 | 2015-10-21 | 大连理工大学 | Linear polarization phase correction method for space universal rotating field azimuth errors |
| CN204679618U (en) * | 2015-06-05 | 2015-09-30 | 中国地震局地球物理研究所 | A kind of three-dimensional orthogonal coil device of field homogeneity distribution |
| CN105203088A (en) * | 2015-09-16 | 2015-12-30 | 中国电子科技集团公司第四十九研究所 | Thee-dimensional magnetic-induction magnetic compass |
| CN105929458A (en) * | 2016-03-21 | 2016-09-07 | 吉林大学 | Aeromagnetic field vector detecting device and detecting method |
| CN106341975A (en) * | 2016-09-14 | 2017-01-18 | 北京航空航天大学 | Hybrid magnetic shielding device based on high temperature superconducting coil |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106772134A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772134B (en) | A kind of apparatus and method of automatic field compensation | |
| Ripka | Electric current sensors: a review | |
| CN111077595B (en) | Method for correcting angle error of superconducting magnetic measurement system and storage medium | |
| CN105203978B (en) | A kind of losing lock of SQUID Magnetic Sensors resets compensation device and method | |
| CN101865982B (en) | Device and method for measuring ambipolar magnetic moment of spatial magnet | |
| US20150369887A1 (en) | Methods and apparatuses related to instrumentation for magnetic relaxometry measurements | |
| CN210222235U (en) | Magnetic sensor's test system | |
| CN102053282A (en) | Borehole three-component magnetic gradient measuring device | |
| CN109633490B (en) | Calibration method of full-tensor magnetic gradient measurement assembly | |
| CN110118948A (en) | A kind of the resultant field measurement method and device of based superconductive quantum inteferometer | |
| CN109358302B (en) | Passive magnetic shielding-free atomic magnetometer device and magnetism measuring method | |
| Sapunov et al. | Ground overhauser DNP geophysical devices | |
| Serson | Instrumentation for induction studies on land | |
| JP2020148760A (en) | measuring device | |
| CN107966670B (en) | Superconducting full-tensor detection device and superconducting full-tensor detection method | |
| CN109283476A (en) | The low frequency intrinsic noise test macro and test method of Magnetic Sensor | |
| Zhi et al. | Digital fluxgate magnetometer for detection of microvibration | |
| CN101923152A (en) | Room temperature calibration method for equivalent error area of gradiometer | |
| CN113155004A (en) | Device and method for calibrating magnetic parameters of reinforcement cage length magnetic detector | |
| Aittoniemi et al. | On balancing superconducting gradiometric magnetometers | |
| Ge et al. | A Novel Coil-Based Overhauser Vector Magnetometer for the Automatic Measurement of Absolute Geomagnetic Total Field and Directions | |
| Post et al. | Design and construction of a magnetic field simulator for cubesat attitude control testing | |
| Qian et al. | Measurement and modeling of magnetic declination error associated with horizontal drift | |
| Deng et al. | Design and development of distributed magnetic measurement system based on hall array | |
| Jiang et al. | Experimental research on the directivity of overhauser magnetometer |
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 |