CN106772648A - Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio - Google Patents
Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio Download PDFInfo
- Publication number
- CN106772648A CN106772648A CN201611104938.3A CN201611104938A CN106772648A CN 106772648 A CN106772648 A CN 106772648A CN 201611104938 A CN201611104938 A CN 201611104938A CN 106772648 A CN106772648 A CN 106772648A
- Authority
- CN
- China
- Prior art keywords
- fid signal
- capacitance
- sensor
- tuning
- frequency
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
Method and system the present invention relates to improve rotary proton class sensor tuning precision and signal to noise ratio, including stimulus sensor exports the first FID signal；Preset Time is waited after the completion of excitation, the first FID signal is gathered, discrete data is generated；Space matrix is built according to discrete data, and uses the svd algorithm singular value decomposition is carried out to space matrix with cancelling noise, obtain reconstruct data；Reconstruct data are processed using fft algorithm, the corresponding first frequency value of maximum peak voltage in the first FID signal frequency spectrum is obtained；The inductance value of sensor and first frequency value are substituted into LC resonance equations and solves the first capacitance, and the first capacitance will be switched to by zero with the capacitance of the tuning capacitance of sensor parallel.The present invention realizes the tuning of rotary proton class sensor using the method that is combined by svd algorithm and fft algorithm, effectively overcomes that existing tuning algorithm tuned speed is slow, tuning precision is low, easily the defects such as detuning phenomena occur under interference environment.
Description
Technical field
The present invention relates to rotary proton class sensor technical field, rotary proton class sensor tuning essence is more particularly to improved
The method and system of degree and signal to noise ratio.
Background technology
Rotary proton class magnetometer is a kind of magnetic-field measurement for measuring slowly varying low-intensity magnetic field or System for Low DC Magnetic Field
Instrument, its sensor is that rotary proton class sensor is inductance element.Measuring principle is by inductance using certain shooting condition
Proton in the solution be active, proton can draw around external magnetic field i.e. magnetic field of the earth is stablized after removing shooting condition
More's precessional motion, produces FID (Free Induction Decay) signal, its precession frequency to be proportional to external magnetic field；Therefore, profit
Sense FID signal with inductance, amplified, shaping and measure frequency, you can obtain external magnetic field value.With other magnetic-field measurements
Technology is compared, rotary proton class magnetometer have high accuracy, high sensitivity the features such as, be widely used and space exploration, near-earth
The fields such as table detection, marine exploration, geomagnetic field measuring, military technology.Because the signal to noise ratio of FID signal is to weigh frequency-measurement accuracy
Key factor, therefore, to increase the signal to noise ratio that sensor exports FID signal, frequency-measurement accuracy is improved, need to be by variable capacitance and sensing
Device is in parallel, is tuned, and this variable capacitance is tuning capacitance.
At present, most rotary proton class magnetometers are used sensor tuning scheme for scanning method, preset capacitance method,
The automatic tracing of blindness and double measurement automatic follow-up tuning method.And the methodical core operating principle of institute is identical：Excite sensing
Device, switches tuning capacitance, the crest voltage of detection output FID signal, and corresponding tuning capacitance capacitance is just at maximum peak voltage
It is the values for tuning of sensor；Only difference is that the means of detection crest voltage.Chinese patent CN103995298A is announced
A kind of method of optimum choice proton magnetometer with humorous electric capacity, in that patent, it is first determined a fixed tuning capacitance
Capacitance, then gradually adjusts to determine final tuning capacitance capacitance.
Existing rotary proton class magnetometer still suffers from problems with the design of sensor tuning algorithm：1) speed is tuned
Degree is universal relatively slow, about needs several seconds, during the actual measurement in field, when change of magnetic field strength is larger, generally requires again to biography
Sensor is tuned, and very big inconvenience is brought to user.2) humorous FID signal pole is not matched somebody with somebody due to the output of rotary proton class magnetometer
Its is faint, it is easy to disturbed by noise signal, and signal to noise ratio is relatively low, once instrument is in the larger environment of interference, frequency spectrum easily occurs
Analytical error is larger to cause " imbalance ", causes the instrument cannot normal work.
The content of the invention
Carried it is an object of the present invention to provide a kind of method and system and one kind for improving rotary proton class sensor tuning precision
The method and system of high proton precession class sensor signal to noise ratio, solve the above-mentioned problems in the prior art.
The technical scheme that the present invention solves above-mentioned technical problem is as follows：
A kind of method for improving rotary proton class sensor tuning precision, comprises the following steps：
Step 1, stimulus sensor exports the first FID signal；
Step 2, waits Preset Time after the completion of excitation, gather first FID signal, generates discrete data；
Step 3, space matrix is built according to the discrete data, and using SVD (singular value decomposition) algorithms to the sky
Between matrix carry out singular value decomposition with cancelling noise, obtain reconstruct data；
Step 4, using data are reconstructed described in FFT (fast Fourier transform) algorithm process, obtains first FID signal frequently
The corresponding first frequency value of maximum peak voltage in spectrum；
Step 5, substitutes into the inductance value of the sensor and the first frequency value LC resonance equations and solves the first electric capacity
Value, and first capacitance will be switched to by zero with the capacitance of the tuning capacitance of the sensor parallel.
The beneficial effects of the invention are as follows：Realize that rotary proton class is passed using the method being combined by svd algorithm and fft algorithm
The tuning of sensor, effectively overcomes that existing tuning algorithm tuned speed is slow, tuning precision is low, easy imbalance occur existing under interference environment
As etc. defect；While tuning period is reduced, the tuning precision of rotary proton class sensor is improved, be capable of achieving Wide measuring range
Selection to tuning capacitance, improve later stage measurement signal signal to noise ratio, be applied to rotary proton class magnetometer, optical pumped magnetometer and
During nuclear magnetic resonance rotary proton FID signal imager etc. is by the instrument of rotary proton class sensor, instrument performance is effectively improved.
On the basis of above-mentioned technical proposal, the present invention can also do following improvement.
Further, the step 3 comprises the following steps：
Step 31, space matrix is built according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
Step 32, uses svd algorithm that singular value decomposition is carried out to the space matrix with cancelling noise, obtains reconstruct number
According to.
Further, also including step 6, stimulus sensor exports the second FID signal again, when waiting default after the completion of excitation
Between, obtain the second frequency value of second FID signal；The inductance value of the sensor and the second frequency value are substituted into LC
Resonance equation solves the second capacitance, and the capacitance of the tuning capacitance is switched into second electricity by first capacitance
Capacitance.
Beneficial effect using above-mentioned further scheme is：On the basis of above-mentioned once tuning, secondary tuning is carried out, entered
One step improves rotary proton class sensor tuning precision.
Further, the step 2 waits Preset Time after the completion of specifically including excitation, and first FID signal is carried out
Amplify, first FID signal after collection amplification generates the discrete data；
Obtained in the step 6 second FID signal second frequency value specifically include by second FID signal according to
It is secondary to be amplified and shaping, according to shaping after second FID signal obtain the second frequency value of second FID signal.
Beneficial effect using above-mentioned further scheme is：FID signal is amplified and/or shaping, further improved
Rotary proton class sensor tuning precision.
Another technical scheme of the invention is as follows：
A kind of method for improving rotary proton class sensor signal to noise ratio, using a kind of above-mentioned raising rotary proton class sensor
The method of tuning precision, also including step 7, stimulus sensor exports the 3rd FID signal again, when waiting default after the completion of excitation
Between, the 3rd FID signal is amplified, using the second frequency value as centre frequency to amplification after the described 3rd
FID signal is filtered, and shaping is carried out to filtered 3rd FID signal, obtains the 3rd FID letters after shaping
Number the 3rd frequency values, the frequency of the FID signal that the 3rd frequency values are detected as the sensor.
The beneficial effects of the invention are as follows：FID signal is entered based on the second frequency value obtained during above-mentioned secondary tuning
Row filtering, while tuning precision is improved, reduces circuit bandwidth, further improves rotary proton class sensor signal to noise ratio.
Another technical scheme of the invention is as follows：
A kind of system for improving rotary proton class sensor tuning precision, including controller, exciting circuit, collector and tune
Humorous circuit；
The controller, it is used to drive exciting circuit stimulus sensor to export the first FID signal；Waited after the completion of excitation
Preset Time, drives collector to gather first FID signal, generates discrete data；Space is built according to the discrete data
Matrix, and use the svd algorithm singular value decomposition is carried out to the space matrix with cancelling noise, obtain reconstruct data；Using FFT
Data are reconstructed described in algorithm process, the corresponding first frequency value of maximum peak voltage in the first FID signal frequency spectrum is obtained；Will
The inductance value of the sensor and the first frequency value substitute into LC resonance equations and solve the first capacitance, and drive tuning circuit
With the capacitance of the tuning capacitance of the sensor parallel first capacitance will be switched to by zero；
The exciting circuit, it is used for the stimulus sensor under the driving of the controller and exports the first FID signal；
The collector, it is used to gather first FID signal under the driving of the controller, generates dispersion number
According to；
The tuning circuit, its be used under the driving of the controller by with the tuning capacitance of the sensor parallel
Capacitance switches to first capacitance by zero.
On the basis of above-mentioned technical proposal, the present invention can also do following improvement.
Further, the controller includes space matrix generation module and SVD denoising modules；
The space matrix generation module, it is used to build space matrix according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
The SVD denoisings module, it is used to using svd algorithm to carry out singular value decomposition to the space matrix is made an uproar with rejecting
Sound, obtains reconstruct data.
Further, also including FPGA (field programmable gate array)；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the second FID signal again, and excitation is completed
After wait Preset Time, drive FPGA to measure the second frequency value of second FID signal；By the inductance value of the sensor and
The second frequency value substitutes into LC resonance equations and solves the second capacitance, and drives tuning circuit by the capacitance of the tuning capacitance
Second capacitance is switched to by first capacitance；
The exciting circuit, it is additionally operable to the stimulus sensor under the driving of the controller and exports the second FID signal；
The FPGA, its second frequency value for being used to be measured under the driving of the controller second FID signal；
The tuning circuit, it is additionally operable to the capacitance of the tuning capacitance under the driving of the controller by described
One capacitance switches to second capacitance.
Further, also including amplifying circuit and comparison circuit；
The amplifying circuit, it is used to wait Preset Time after the completion of encouraging, first FID signal is amplified,
For collector collection；It is additionally operable to wait Preset Time after the completion of excitation, and second FID signal is amplified,
For the comparison circuit shaping；
The comparison circuit, second FID signal after it is used for amplification carries out shaping, is surveyed for the FPGA
Amount；
The controller is STM32, and the collector is ADC collectors.
Beneficial effect using above-mentioned further scheme is：The high speed being combined using ADC collectors, STM32 and FPGA
Degree collection and treatment hardware structure, while rotary proton class sensor tuning precision is increased, are effectively ensured tuned speed.
Another technical scheme of the invention is as follows：
A kind of system for improving rotary proton class sensor signal to noise ratio, including a kind of above-mentioned raising rotary proton class sensor
The system of tuning precision, also including narrow-band filtering circuit；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the 3rd FID signal again；It is additionally operable to
The centre frequency of the narrow-band filtering circuit is set to the second frequency value；After it is additionally operable to drive FPGA measurement shapings
3rd frequency values of the 3rd FID signal, and the FID signal that the 3rd frequency values are detected as the sensor
Rate-adaptive pacemaker；
The amplifying circuit, it is additionally operable to wait Preset Time after the completion of excitation, the 3rd FID signal is put
Greatly, for the narrow-band filtering circuit filtering；
The narrow-band filtering circuit, its be used for using the second frequency value as centre frequency to amplification after the described 3rd
FID signal is filtered, for the comparison circuit shaping；
The comparison circuit, it is additionally operable to carry out shaping to filtered 3rd FID signal, is surveyed for the FPGA
Amount；
The FPGA, it is additionally operable to measure the of the 3rd FID signal after shaping under the driving of the controller
Three frequency values.
Brief description of the drawings
Fig. 1 is a kind of flow chart of the method for improving rotary proton class sensor tuning precision of the present invention；
Fig. 2 is a kind of theory diagram of the system for improving rotary proton class sensor signal to noise ratio of the present invention；
Fig. 3 is that the amplified circuit that ADC collectors are collected amplifies the untuned undressed spectrogram of FID signal；
Fig. 4 is after the amplified circuit untuned FID signal of amplification that ADC collectors are collected is processed through auto-correlation algorithm
Spectrogram；
Fig. 5 be the amplified circuit that collects of ADC collectors amplify untuned FID signal through svd algorithm of the present invention at
The spectrogram of gained reconstruct data composition after reason.
Specific embodiment
Principle of the invention and feature are described below in conjunction with accompanying drawing, example is served only for explaining the present invention, and
It is non-for limiting the scope of the present invention.
As shown in figure 1, described in the embodiment of the present invention 1 it is a kind of improve rotary proton class sensor tuning precision method,
Comprise the following steps：
Step 1, stimulus sensor exports the first FID signal；
Step 2, waits Preset Time after the completion of excitation, gather first FID signal, generates discrete data；
Step 3, builds space matrix, and the space matrix is carried out very using svd algorithm according to the discrete data
Different value is decomposed with cancelling noise, obtains reconstruct data；
Step 4, the reconstruct data are processed using fft algorithm, obtain peak-peak electricity in the first FID signal frequency spectrum
Press corresponding first frequency value；
Step 5, substitutes into the inductance value of the sensor and the first frequency value LC resonance equations and solves the first electric capacity
Value, and first capacitance will be switched to by zero with the capacitance of the tuning capacitance of the sensor parallel.
Wherein, the LC resonance equations are as follows：
Wherein, the f0, L and C be respectively frequency variable, inductance variable and electric capacity variable.
According to the difference of sensor type, using exciting circuit to solution input different frequency pulse signal where sensor
Carry out stimulus sensor output FID signal, the actuation duration is 400ms in specific implementation, after waiting 5ms after the completion of excitation, gather FID
Signal, beneficial to the interference excluded produced by circuit oscillation.
Svd algorithm is also singular value decomposition method, is the effective ways that a kind of data characteristics is extracted, based on nonlinear filtering side
Method, the value that it is decomposed reflects the inherent attribute of data, therefore can realize that the background in small target deteection suppresses and Denoising Problems.
Svd algorithm has successfully been incorporated into geophysical SNR estimation and compensation technology, is compared with other denoising methods, denoising effect
It is good, and useful signal distortion is small.The FID signal treatment of rotary proton class sensor output is under the jurisdiction of geophysical noise point
From the background fallen within small target deteection suppresses and Denoising Problems, therefore carries out denoising, energy to FID signal using svd algorithm
Enough obtain good denoising effect.
A kind of method for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 2, in the base of embodiment 1
On plinth, the step 3 comprises the following steps：
Step 31, space matrix is built according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
The DmIn each row vector correspondence space matrix in a point.
Step 32, uses svd algorithm that singular value decomposition is carried out to the space matrix with cancelling noise, obtains reconstruct number
According to.
A kind of method for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 3, in embodiment 1 or 2
On the basis of, also including step 6, stimulus sensor exports the second FID signal again, and Preset Time is waited after the completion of excitation, obtains
Take the second frequency value of second FID signal；The inductance value of the sensor and the second frequency value are substituted into LC resonance
The capacitance of equations second, and the capacitance of the tuning capacitance is switched into second electric capacity by first capacitance
Value.
A kind of method for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 4, in the base of embodiment 3
On plinth, the step 2 waits Preset Time after the completion of specifically including excitation, and first FID signal is amplified, and collection is put
First FID signal after big, generates the discrete data；
Obtained in the step 6 second FID signal second frequency value specifically include by second FID signal according to
It is secondary to be amplified and shaping, according to shaping after second FID signal obtain the second frequency value of second FID signal.
A kind of method for improving rotary proton class sensor signal to noise ratio described in the embodiment of the present invention 5, on the basis of embodiment 4
On, also including step 7, stimulus sensor exports the 3rd FID signal again, and Preset Time is waited after the completion of excitation, to described the
Three FID signals are amplified, using the second frequency value as centre frequency to amplification after the 3rd FID signal filter
Ripple, and shaping is carried out to filtered 3rd FID signal, obtain the 3rd frequency of the 3rd FID signal after shaping
Value, the frequency of the FID signal that the 3rd frequency values are detected as the sensor.
As shown in Fig. 2 described in the embodiment of the present invention 1 it is a kind of improve rotary proton class sensor tuning precision system,
Including controller, exciting circuit, collector and tuning circuit；
The controller, it is used to drive exciting circuit stimulus sensor to export the first FID signal；Waited after the completion of excitation
Preset Time, drives collector to gather first FID signal, generates discrete data；Space is built according to the discrete data
Matrix, and use the svd algorithm singular value decomposition is carried out to the space matrix with cancelling noise, obtain reconstruct data；Using FFT
Data are reconstructed described in algorithm process, the corresponding first frequency value of maximum peak voltage in the first FID signal frequency spectrum is obtained；Will
The inductance value of the sensor and the first frequency value substitute into LC resonance equations and solve the first capacitance, and drive tuning circuit
With the capacitance of the tuning capacitance of the sensor parallel first capacitance will be switched to by zero；
The exciting circuit, it is used for the stimulus sensor under the driving of the controller and exports the first FID signal；
The collector, it is used to gather first FID signal under the driving of the controller, generates dispersion number
According to；
The tuning circuit, its be used under the driving of the controller by with the tuning capacitance of the sensor parallel
Capacitance switches to first capacitance by zero.
A kind of system for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 2, in the base of embodiment 1
On plinth, the controller includes space matrix generation module and SVD denoising modules；
The space matrix generation module, it is used to build space matrix according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
The SVD denoisings module, it is used to using svd algorithm to carry out singular value decomposition to the space matrix is made an uproar with rejecting
Sound, obtains reconstruct data.
A kind of system for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 3, in embodiment 1 or 2
On the basis of, also including FPGA；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the second FID signal again, and excitation is completed
After wait Preset Time, drive FPGA to measure the second frequency value of second FID signal；By the inductance value of the sensor and
The second frequency value substitutes into LC resonance equations and solves the second capacitance, and drives tuning circuit by the capacitance of the tuning capacitance
Second capacitance is switched to by first capacitance；
The exciting circuit, it is additionally operable to the stimulus sensor under the driving of the controller and exports the second FID signal；
The FPGA, its second frequency value for being used to be measured under the driving of the controller second FID signal；
The tuning circuit, it is additionally operable to the capacitance of the tuning capacitance under the driving of the controller by described
One capacitance switches to second capacitance.
A kind of system for improving rotary proton class sensor tuning precision described in the embodiment of the present invention 4, in the base of embodiment 3
On plinth, also including amplifying circuit and comparison circuit；
The amplifying circuit, it is used to wait Preset Time after the completion of encouraging, first FID signal is amplified,
For collector collection；It is additionally operable to wait Preset Time after the completion of excitation, and second FID signal is amplified,
For the comparison circuit shaping；
The comparison circuit, second FID signal after it is used for amplification carries out shaping, is surveyed for the FPGA
Amount；
The controller is STM32, and the collector is ADC collectors.Because FID signal is exponentially decayed, therefore must
Suitable A/D sample rates and sampling number must be set, it is known that magnetic field of the earth scope is 20,000nT~100,000nT, according to magnetic
Revolve and can obtain the frequency range of FID signal for 850Hz~4 than formula, 300Hz, therefore, in specific embodiment, ADC collectors are set
Sample rate be 10kHz, sampling number be 2048 points, collection signal time be about 205ms, frequency resolution is 10kHz/
2048=4.88Hz.
A kind of system for improving rotary proton class sensor signal to noise ratio described in the embodiment of the present invention 5, on the basis of embodiment 4
On, also including narrow-band filtering circuit；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the 3rd FID signal again；It is additionally operable to
The centre frequency of the narrow-band filtering circuit is set to the second frequency value；After it is additionally operable to drive FPGA measurement shapings
3rd frequency values of the 3rd FID signal, and the FID signal that the 3rd frequency values are detected as the sensor
Rate-adaptive pacemaker；
The amplifying circuit, it is additionally operable to wait Preset Time after the completion of excitation, the 3rd FID signal is put
Greatly, for the narrow-band filtering circuit filtering；
The narrow-band filtering circuit, its be used for using the second frequency value as centre frequency to amplification after the described 3rd
FID signal is filtered, for the comparison circuit shaping；
The comparison circuit, it is additionally operable to carry out shaping to filtered 3rd FID signal, is surveyed for the FPGA
Amount；
The FPGA, it is additionally operable to measure the of the 3rd FID signal after shaping under the driving of the controller
Three frequency values.
Fig. 3, Fig. 4 and Fig. 5 are followed successively by tuning capacitance when being zero, and the amplified circuit that ADC collectors are collected amplifies not to be adjusted
After the undressed frequency spectrum of humorous FID signal, this FID signal are processed through existing auto-correlation algorithm (Auto Correlation)
Gained reconstructs the frequency spectrum that data are constituted after frequency spectrum and this FID signal are processed through svd algorithm of the present invention；Contrast can be found through the present invention
Frequency spectrum after svd algorithm treatment is distincter, and signal to noise ratio is higher.
Using same experiment porch, to conventional at present tuning method and this hair based on peak detection method and auto-correlation algorithm
The bright secondary tuning method based on svd algorithm is tuned the contrast experiment of precision and speed, with the timer of microcontroller STM32
Timing is carried out to the time used by three kinds of tuning algorithms.According to the design parameter of experiment porch, in same test magnetic field environment：Survey
Examination place magnetic field is about 49,323nT, i.e. 2100Hz；Sensor resonant -3dB frequency ranges：2072Hz~2128Hz, sensor
Inductance value 34mH, tuning capacitance capacitance scope is obtained with reference to LC resonance equations：168nF~173nF.Therefore, as long as three kinds of algorithms are obtained
The tuning capacitance capacitance for arriving is within this range, you can illustrate to tune successfully.In glitch-free environment, three kinds of methods are entered respectively
5 observation of row, its result is as shown in table 1.
The measurement result of the lower three kinds of methods of the noiseless environment of table 1
As can be seen from Table 1, from the point of view of the speed of tuning, because three kinds of algorithms are sequential programme structure, so program
The execution time is definite value, and the time used by peak detection method is about 5 times of auto-correlation algorithm and svd algorithm；Come from the precision of tuning
See, the tuning capacitance value obtained by three kinds of algorithms in the range of 168nF~173nF, more accurately.It can be said that it is bright,
Under noiseless environment, the performance of auto-correlation algorithm and svd algorithm is suitable.
In the environment of interference, 5 observation is carried out to above-mentioned three kinds of methods respectively, its result is as shown in table 2.
The measurement result of the lower three kinds of methods of the interference environment of table 2
As can be seen from Table 2, when external environment has interference, the tuning precision of peak detection method and auto-correlation algorithm is bright
Aobvious to decline, the tuning capacitance capacitance for obtaining can make rotary proton class sensor detuning phenomena occur, reduce the signal to noise ratio of signal, from
And cause that application has the instrument of rotary proton class sensor cannot normal work；And the test result of svd algorithm is not taken office
What influences, and tuning precision is high, speed is fast and with repeatability.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all it is of the invention spirit and
Within principle, any modification, equivalent substitution and improvements made etc. should be included within the scope of the present invention.
Claims (10)
1. it is a kind of improve rotary proton class sensor tuning precision method, it is characterised in that comprise the following steps：
Step 1, stimulus sensor exports the first FID signal；
Step 2, waits Preset Time after the completion of excitation, gather first FID signal, generates discrete data；
Step 3, builds space matrix, and carry out singular value to the space matrix using svd algorithm according to the discrete data
Decompose with cancelling noise, obtain reconstruct data；
Step 4, the reconstruct data are processed using fft algorithm, obtain maximum peak voltage pair in the first FID signal frequency spectrum
The first frequency value answered；
Step 5, substitutes into the inductance value of the sensor and the first frequency value LC resonance equations and solves the first capacitance, and
With the capacitance of the tuning capacitance of the sensor parallel first capacitance will be switched to by zero.
2. according to claim 1 it is a kind of improve rotary proton class sensor tuning precision method, it is characterised in that it is described
Step 3 comprises the following steps：
Step 31, space matrix is built according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
Step 32, uses svd algorithm that singular value decomposition is carried out to the space matrix with cancelling noise, obtains reconstruct data.
3. it is according to claim 1 or claim 2 it is a kind of improve rotary proton class sensor tuning precision method, it is characterised in that
Also include step 6, stimulus sensor exports the second FID signal again, Preset Time is waited after the completion of excitation, obtain described second
The second frequency value of FID signal；The inductance value of the sensor and the second frequency value are substituted into LC resonance equations and solves the
Two capacitances, and the capacitance of the tuning capacitance is switched into second capacitance by first capacitance.
4. according to claim 3 it is a kind of improve rotary proton class sensor tuning precision method, it is characterised in that it is described
Step 2 waits Preset Time after the completion of specifically including excitation, and first FID signal is amplified, the institute after collection amplification
The first FID signal is stated, the discrete data is generated；
The second frequency value that second FID signal is obtained in the step 6 specifically includes successively to enter second FID signal
Row amplifies and shaping, according to shaping after second FID signal obtain the second frequency value of second FID signal.
5. it is a kind of improve rotary proton class sensor signal to noise ratio method, it is characterised in that carried using one kind described in claim 4
The method of high proton precession class sensor tuning precision, also including step 7, stimulus sensor exports the 3rd FID signal again, swashs
Preset Time is waited after the completion of encouraging, the 3rd FID signal is amplified, using the second frequency value as centre frequency pair
The 3rd FID signal after amplification is filtered, and shaping is carried out to filtered 3rd FID signal, obtains shaping
3rd frequency values of the 3rd FID signal afterwards, the FID signal that the 3rd frequency values are detected as the sensor
Frequency.
6. it is a kind of improve rotary proton class sensor tuning precision system, it is characterised in that including controller, exciting circuit,
Collector and tuning circuit；
The controller, it is used to drive exciting circuit stimulus sensor to export the first FID signal；Waited after the completion of excitation default
Time, drive collector to gather first FID signal, generate discrete data；Space matrix is built according to the discrete data,
And use the svd algorithm singular value decomposition is carried out to the space matrix with cancelling noise, obtain reconstruct data；Using fft algorithm
The reconstruct data are processed, the corresponding first frequency value of maximum peak voltage in the first FID signal frequency spectrum is obtained；Will be described
The inductance value of sensor and the first frequency value substitute into LC resonance equations and solve the first capacitance, and drive the tuning circuit will be with
The capacitance of the tuning capacitance of the sensor parallel switches to first capacitance by zero；
The exciting circuit, it is used for the stimulus sensor under the driving of the controller and exports the first FID signal；
The collector, it is used to gather first FID signal under the driving of the controller, generates discrete data；
The tuning circuit, it is used for the capacitance with the tuning capacitance of the sensor parallel under the driving of the controller
First capacitance is switched to by zero.
7. according to claim 6 it is a kind of improve rotary proton class sensor tuning precision system, it is characterised in that it is described
Controller includes space matrix generation module and SVD denoising modules；
The space matrix generation module, it is used to build space matrix according to following first formula；
First formula is as follows：
Wherein, the DmIt is the space matrix, { x1,x2,x3…xn,xn+1…xn+m-1It is the discrete data；
The SVD denoisings module, it is used to using the svd algorithm singular value decomposition is carried out to the space matrix with cancelling noise,
Obtain reconstruct data.
8. according to claim 6 or 7 it is a kind of improve rotary proton class sensor tuning precision system, it is characterised in that
Also include FPGA；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the second FID signal again, after the completion of excitation etc.
Preset Time is treated, drives FPGA to measure the second frequency value of second FID signal；By the inductance value of the sensor and described
Second frequency value substitutes into LC resonance equations and solves the second capacitance, and drives tuning circuit by the capacitance of the tuning capacitance by institute
State the first capacitance and switch to second capacitance；
The exciting circuit, it is additionally operable to the stimulus sensor under the driving of the controller and exports the second FID signal；
The FPGA, its second frequency value for being used to be measured under the driving of the controller second FID signal；
The tuning circuit, it is additionally operable to the capacitance of the tuning capacitance under the driving of the controller by the described first electricity
Capacitance switches to second capacitance.
9. according to claim 8 it is a kind of improve rotary proton class sensor tuning precision system, it is characterised in that also wrap
Include amplifying circuit and comparison circuit；
The amplifying circuit, it is used to wait Preset Time after the completion of encouraging, first FID signal is amplified, is used for
The collector collection；It is additionally operable to wait Preset Time after the completion of excitation, and second FID signal is amplified, and is used for
The comparison circuit shaping；
The comparison circuit, second FID signal after it is used for amplification carries out shaping, for FPGA measurements；
The controller is STM32, and the collector is ADC collectors.
10. a kind of system for improving rotary proton class sensor signal to noise ratio, it is characterised in that including a kind of described in claim 9
The system for improving rotary proton class sensor tuning precision, also including narrow-band filtering circuit；
The controller, it is additionally operable to drive exciting circuit stimulus sensor to export the 3rd FID signal again；It is additionally operable to institute
The centre frequency for stating narrow-band filtering circuit is set to the second frequency value；It is additionally operable to drive FPGA measurement shaping after it is described
3rd frequency values of the 3rd FID signal, and the FID signal that the 3rd frequency values are detected as the sensor frequency
Output；
The amplifying circuit, it is additionally operable to wait Preset Time after the completion of excitation, the 3rd FID signal is amplified, uses
In the narrow-band filtering circuit filtering；
The narrow-band filtering circuit, its be used for using the second frequency value as centre frequency to amplification after the 3rd FID
Signal is filtered, for the comparison circuit shaping；
The comparison circuit, it is additionally operable to carry out shaping to filtered 3rd FID signal, for FPGA measurements；
The FPGA, its 3rd frequency for being additionally operable to be measured under the driving of the controller the 3rd FID signal after shaping
Rate value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611104938.3A CN106772648A (en) | 2016-12-05 | 2016-12-05 | Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611104938.3A CN106772648A (en) | 2016-12-05 | 2016-12-05 | Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106772648A true CN106772648A (en) | 2017-05-31 |
Family
ID=58874079
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611104938.3A Pending CN106772648A (en) | 2016-12-05 | 2016-12-05 | Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772648A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109597137A (en) * | 2018-12-28 | 2019-04-09 | 中国地质大学（武汉） | Overhauser magnetometer based on semiconductor magnetic sensor is quickly tracked with humorous method |
| CN110333546A (en) * | 2019-06-28 | 2019-10-15 | 中国地质大学（武汉） | Improve the system and method for proton magnetic precession Signal-to-Noise |
| CN110333545A (en) * | 2019-06-28 | 2019-10-15 | 中国地质大学（武汉） | Improve the system and method for rotary proton class sensor tuning precision and signal-to-noise ratio |
| CN110632535A (en) * | 2019-08-29 | 2019-12-31 | 中国地质大学（武汉） | Noise suppression system and method suitable for proton magnetic precession signal |
| CN113655534A (en) * | 2021-07-14 | 2021-11-16 | 中国地质大学（武汉） | Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1841083A (en) * | 2005-04-01 | 2006-10-04 | 西门子公司 | Method for avoiding linear phase error in magnetic resonance spectroscopy |
| CN101839968A (en) * | 2010-04-09 | 2010-09-22 | 中国地质大学（武汉） | OVERHAUSER magnetometer |
| CN101963656A (en) * | 2009-07-24 | 2011-02-02 | 西门子迈迪特(深圳)磁共振有限公司 | Complementary field detection method and device |
| US20120191417A1 (en) * | 2010-09-24 | 2012-07-26 | University Of Southern California | Noise reduction for spectroscopic signal processing |
| CN103620437A (en) * | 2011-06-27 | 2014-03-05 | 皇家飞利浦有限公司 | Bone mri using an ultrashort echo time pulse sequence with fid and multiple gradient echo acquisition and water-fat separation processing |
| CN104808251A (en) * | 2015-05-18 | 2015-07-29 | 中国地质大学（武汉） | Method for improving frequency measuring precision of larmor signal of Overhauser magnetometer and circuit thereof |
-
2016
- 2016-12-05 CN CN201611104938.3A patent/CN106772648A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1841083A (en) * | 2005-04-01 | 2006-10-04 | 西门子公司 | Method for avoiding linear phase error in magnetic resonance spectroscopy |
| CN101963656A (en) * | 2009-07-24 | 2011-02-02 | 西门子迈迪特(深圳)磁共振有限公司 | Complementary field detection method and device |
| CN101839968A (en) * | 2010-04-09 | 2010-09-22 | 中国地质大学（武汉） | OVERHAUSER magnetometer |
| US20120191417A1 (en) * | 2010-09-24 | 2012-07-26 | University Of Southern California | Noise reduction for spectroscopic signal processing |
| CN103620437A (en) * | 2011-06-27 | 2014-03-05 | 皇家飞利浦有限公司 | Bone mri using an ultrashort echo time pulse sequence with fid and multiple gradient echo acquisition and water-fat separation processing |
| CN104808251A (en) * | 2015-05-18 | 2015-07-29 | 中国地质大学（武汉） | Method for improving frequency measuring precision of larmor signal of Overhauser magnetometer and circuit thereof |
Non-Patent Citations (1)
| Title |
|---|
| HUAN LIU等: "Research on a secondary tuning algorithm based on SVD ＆ STFT for FID signal", 《MEASUREMENT SCIENCE AND TECHNOLOGY》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109597137A (en) * | 2018-12-28 | 2019-04-09 | 中国地质大学（武汉） | Overhauser magnetometer based on semiconductor magnetic sensor is quickly tracked with humorous method |
| CN109597137B (en) * | 2018-12-28 | 2020-09-04 | 中国地质大学（武汉） | Overhauser magnetometer tracking and tuning method based on semiconductor magnetic sensor |
| CN110333546A (en) * | 2019-06-28 | 2019-10-15 | 中国地质大学（武汉） | Improve the system and method for proton magnetic precession Signal-to-Noise |
| CN110333545A (en) * | 2019-06-28 | 2019-10-15 | 中国地质大学（武汉） | Improve the system and method for rotary proton class sensor tuning precision and signal-to-noise ratio |
| CN110632535A (en) * | 2019-08-29 | 2019-12-31 | 中国地质大学（武汉） | Noise suppression system and method suitable for proton magnetic precession signal |
| CN110632535B (en) * | 2019-08-29 | 2021-07-13 | 中国地质大学（武汉） | Noise suppression system and method suitable for proton magnetic precession signal |
| CN113655534A (en) * | 2021-07-14 | 2021-11-16 | 中国地质大学（武汉） | Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition |
| CN113655534B (en) * | 2021-07-14 | 2022-05-17 | 中国地质大学（武汉） | Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772648A (en) | Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio | |
| EP3455936B1 (en) | Systems, methods and programs for denoising signals using wavelets | |
| JP4114989B2 (en) | Method for compensating the magnetic field of a magnetic resonance system | |
| US4431968A (en) | Method of three-dimensional NMR imaging using selective excitation | |
| JP3857752B2 (en) | Method and apparatus for collecting MR data | |
| US6420872B1 (en) | Probe for detecting a transient magnetic resonance signal, wherein the ratio of the Q of the probe to the Q of the resonance signal is relatively large | |
| CN110333545A (en) | Improve the system and method for rotary proton class sensor tuning precision and signal-to-noise ratio | |
| CN104237820B (en) | A kind of single sweep obtains the method that magnetic resonance two dimension J decomposes spectrum | |
| CN104142484B (en) | With two HF motivational techniques of resonant frequency, the method and apparatus of collection MR data | |
| US10222439B2 (en) | System and method for spiral multislab magnetic resonance imaging | |
| Liu et al. | Research on a secondary tuning algorithm based on SVD & STFT for FID signal | |
| JPH0113054B2 (en) | ||
| US4626784A (en) | NMR imaging device | |
| CN110632535B (en) | Noise suppression system and method suitable for proton magnetic precession signal | |
| CN111044960B (en) | Solid nuclear magnetic resonance quantitative detection method and device based on successive cross polarization | |
| JP6023084B2 (en) | Nuclear magnetic resonance imaging apparatus and nuclear magnetic resonance imaging method | |
| CN102028469A (en) | Digital radio-frequency receiver for magnetic resonance imaging system and signal processing method thereof | |
| US20110288802A1 (en) | Spectrometer, Measuring Apparatus, and Method of Data Processing | |
| CN113655534B (en) | Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition | |
| US4424488A (en) | Nuclear magnetic resonance apparatus | |
| CN110333546A (en) | Improve the system and method for proton magnetic precession Signal-to-Noise | |
| US7064544B1 (en) | Method and system of scaling MR spectroscopic data acquired with phased-array coils | |
| Liu et al. | An Improved Tuning Control Algorithm Based on SVD for FID Signal | |
| CN105807243B (en) | A kind of gradient shimming method based on phase code | |
| CN208420163U (en) | A kind of nmr probe performance testing device |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |