CN104545913A - Adaptive compensation method and system for environmental magnetic field suppression - Google Patents

Adaptive compensation method and system for environmental magnetic field suppression Download PDF

Info

Publication number
CN104545913A
CN104545913A CN201310492564.7A CN201310492564A CN104545913A CN 104545913 A CN104545913 A CN 104545913A CN 201310492564 A CN201310492564 A CN 201310492564A CN 104545913 A CN104545913 A CN 104545913A
Authority
CN
China
Prior art keywords
data
sigma
centerdot
magnetic field
rmj
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310492564.7A
Other languages
Chinese (zh)
Other versions
CN104545913B (en
Inventor
张树林
王永良
孔祥燕
谢晓明
李华
曾佳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Man Di medical instruments (Shanghai) Co., Ltd.
Original Assignee
Shanghai Institute of Microsystem and Information Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Microsystem and Information Technology of CAS filed Critical Shanghai Institute of Microsystem and Information Technology of CAS
Priority to CN201310492564.7A priority Critical patent/CN104545913B/en
Publication of CN104545913A publication Critical patent/CN104545913A/en
Application granted granted Critical
Publication of CN104545913B publication Critical patent/CN104545913B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7296—Specific aspects of physiological measurement analysis for compensation of signal variation due to stress unintentionally induced in the patient, e.g. due to the stress of the medical environment or examination

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Psychiatry (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Developmental Disabilities (AREA)
  • Signal Processing (AREA)
  • Child & Adolescent Psychology (AREA)
  • Psychology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The invention provides an adaptive compensation method and system for environmental magnetic field suppression. The method includes acquiring signal channel signals and reference channel signals output by a detection system to acquire the signal channel data and reference channel data; extracting the reference channel data from an initial time point through a window with the fixed length, increasing gradually according to the initial time point from 0 by adopting 1 as the gradient, and acquiring data of a plurality of compensation channels sequentially; establishing an error function containing a compensation coefficient, and acquiring the compensation coefficient by solving the error function according to the minimum variance criterion; performing environmental noise suppression on the signal channel data according to the compensation coefficient. Once a parameter of the length of the window is set, an initial value of a subtraction coefficient is omitted, the convergence process is avoided, the optimal subtraction coefficient can be acquired stably in an adaptive manner, and the gradient output with fine performance can be implemented.

Description

A kind of self-adapting compensation method and system being applied to environmental magnetic field suppression
Technical field
The invention belongs to signal detection technique field, relate to a kind of environmental magnetic field suppressing method, particularly relate to a kind of self-adapting compensation method and the system that are applied to environmental magnetic field suppression.
Background technology
Vital movement is richly stored with electromagnetic information behind.Biological magnetic, derive from the space magnetic field that bioelectric produces, its intensity is very faint, and as become the Typical strengths of popular feeling magnetic to only have 100pT, brain magnetic is then more faint, only has hundred fT magnitudes.Relative to biological magnetic, environmental magnetic field is very strong, and the Typical strengths as magnetic field of the earth is 30-50 μ T, and the change in urban environment magnetic field also reaches hundreds of nT.In order to detect faint biological magnetic signal in so strong environmental magnetic field, except the high-sensitive Magnetic Sensor of needs, also need high performance environmental magnetic field suppressing method.
At present, the main suppression means of environmental magnetic field have magnetic shield room and gradiometer technology.The environmental magnetic field that magnetic shield room adopts metal (as aluminium and the copper) shielding of the mu-metal (as permalloy) of high magnetic permeability and high conductivity extraneous, although the most effectively, involve great expense.Comparatively speaking, gradiometer technology is most economical and is most widely used.Gradiometer, is similar to a kind of spatial filter, is to produce response to magnetic field space gradient, according to different response exponent numbers, has single order, second order, more high-order etc.Ideally, n gradiometer is to uniform field with until the gradient field response on n-1 rank is zero.But desirable gradiometer is non-existent, due to the unbalanced error of self, actual hardware gradiometer inevitably produces certain response to the gradient magnetic of uniform magnetic field and low order, causes the reduction of gradiometer squelch performance.
In order to improve the detection performance of gradiometer, the most frequently used is composite noise suppression at present.Consider the signalling channel of the biological magnetic of detection and multiple references in the certain space distance basic detection background magnetic field that is separated by, composite noise suppresses to can be understood as utilization and subtracts each other coefficient, the output of deducting signalling channel is gone, to construct the method for better squelch capability gradient by reference to the linear combination exported.Thus, based on certain signal and reference channel, subtract each other the squelch performance that coefficient directly determines composition gradient.At present, the implementation method of subtracting each other coefficient is that electronics subtracts each other, software subtracts each other and self-adaptation variable coefficient.Electronics subtracts each other and to subtract each other with software is a kind of mode of fixed coefficient, implement the simplest, but when environmental magnetic field feature changes in time, the performance of determining coefficient will reduce greatly.Compared to fixed coefficient, self-adaptation variable coefficient can environmentally changes of magnetic field automatically adjust and subtracts each other coefficient, thus ensures certain squelch performance.But current self-adaptation variable coefficient adopts the method for setting initial value and convergence automatically, causes subtracting each other coefficient and occur very large error within initial a period of time, and there is unstable characteristic.
Above scaling method respectively has relative merits, seeks high stable and adaptive compensation method, detects have great importance for neighbourhood noise suppression and biological magnetic signal.
Summary of the invention
The shortcoming of prior art in view of the above, the object of the present invention is to provide a kind of self-adapting compensation method and the system that are applied to environmental magnetic field suppression, the method of setting initial value and convergence is automatically adopted for solving existing self-adaptation variable coefficient, cause subtracting each other coefficient and within initial a period of time, occur very large error, the problem of self-adaptation instability.
For achieving the above object and other relevant objects, the invention provides and be a kind ofly applied to the self-adapting compensation method and system that environmental magnetic field suppresses, wherein:
Be applied to the self-adapting compensation method that environmental magnetic field suppresses, comprise: gather signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports; The window utilizing length to be n is from initial time point n 0start to extract described reference channel data, along with initial time point n 0from 0, be increased to N-n+1 gradually with 1 for gradient, obtain m compensation channels data successively, wherein N is the data total length of single passage; Instrument error function wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; Be that criterion solves error function acquisition penalty coefficient k with minimum variance i, 1≤i≤m; Utilize penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1,0≤n 0≤ m-1, for n-th of the signalling channel data after compensation 0+ a data point.
Preferably, the length n of described window immobilizes, and n is lowest mean square data length used.
Preferably, be that criterion solves error function acquisition penalty coefficient k with minimum variance idetailed process comprise: when error function f reaches minimum, noise suppression effect is best, then have can obtain thus about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ,
Solve described penalty coefficient matrix, obtain with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i.
Preferably, described detection system comprises biological magnetic detection system, and described biological magnetic detection system is provided with signalling channel and reference channel.
Preferably, when the product of m and n equals the total length of all reference channel data, be the situation that self-adaptation determines coefficient.
Be applied to the adaptive equalization system that environmental magnetic field suppresses, comprise: data acquisition module, gather signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports; Slip scan module, is connected with described data acquisition module, and the window utilizing length to be n is from initial time point n 0start to carry out slip scan extraction to described reference channel data, each distance of sliding is 1 data point, and obtain m compensation channels data successively, the product of m and n is less than or equal to the number of data points of described reference channel data; Error function constructing module, is connected with slip scan module respectively with described data acquisition module, and the error function constructed is wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; 1≤i≤m; Penalty coefficient acquisition module, is connected with described error function constructing module, is that criterion solves described error function acquisition penalty coefficient k with minimum variance i; Signal compensation module, is connected with described penalty coefficient acquisition module, utilizes penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1, for n-th of the signalling channel data after compensation 0+ a data point.
Preferably, the length n of described window immobilizes, initial time point n 0increase gradually for gradient with 1 from 0, n is lowest mean square data length used.
Preferably, described penalty coefficient acquisition module comprises: condition setting module, Matrix Solving module; Condition setting module is connected with described error function constructing module, will according to minimum variance principle be set as 0, obtain about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ;
Matrix Solving module is connected with described condition setting module, solves described penalty coefficient matrix, obtains with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i, wherein N is the data total length of single passage.
Preferably, described detection system comprises biological magnetic detection system, and described biological magnetic detection system is provided with signalling channel and reference channel.
As mentioned above, self-adapting compensation method and the system being applied to environmental magnetic field suppression of the present invention, has following beneficial effect:
The present invention only need arrange a parameter n, subtracts each other coefficient k without the need to setting iinitial value, there is not the process of convergence, namely Absorbable organic halogens obtains adaptively and best subtracts each other coefficient k i, realize high performance gradient and export.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet being applied to the self-adapting compensation method that environmental magnetic field suppresses of the present invention.
Fig. 2 is the structural representation being applied to the adaptive equalization system that environmental magnetic field suppresses of the present invention.
Element numbers explanation
210 data acquisition modules
220 slip scan modules
230 error function constructing modules
240 penalty coefficient acquisition modules
241 condition setting modules
242 Matrix Solving modules
250 signal compensation modules
Embodiment
Below by way of specific instantiation, embodiments of the present invention are described, those skilled in the art the content disclosed by this instructions can understand other advantages of the present invention and effect easily.The present invention can also be implemented or be applied by embodiments different in addition, and the every details in this instructions also can based on different viewpoints and application, carries out various modification or change not deviating under spirit of the present invention.
Refer to accompanying drawing.It should be noted that, the diagram provided in the present embodiment only illustrates basic conception of the present invention in a schematic way, then only the assembly relevant with the present invention is shown in graphic but not component count, shape and size when implementing according to reality is drawn, it is actual when implementing, and the kenel of each assembly, quantity and ratio can be a kind of change arbitrarily, and its assembly layout kenel also may be more complicated.
Below in conjunction with embodiment and accompanying drawing, the present invention is described in detail.
Embodiment
The present embodiment provides a kind of self-adapting compensation method being applied to environmental magnetic field and suppressing, and as shown in Figure 1, the method comprises:
Gather signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports.Described detection system includes but not limited to biological magnetic detection system, and described biological magnetic detection system is provided with signalling channel and reference channel.The signalling channel signal that detection system exports and reference channel signal normally simulating signal, sampled by A/D, the signalling channel signal of simulation and reference channel signal can be converted into signalling channel data and the reference channel data of numeral.
The window utilizing length to be n is from initial time point n 0start to extract described reference channel data, along with initial time point n 0from 0, be increased to N-n+1 gradually with 1 for gradient, obtain m compensation channels data successively, wherein N is the data total length of single passage; The length n of described window immobilizes, and n is lowest mean square data length used, and the product of m and n is less than or equal to the number of data points of described reference channel data.
Instrument error function wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; Be that criterion solves error function acquisition penalty coefficient k with minimum variance i, 1≤i≤m.Further, be that criterion solves error function acquisition penalty coefficient k with minimum variance idetailed process comprise: when error function f reaches minimum, noise suppression effect is best, then have can obtain thus about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ,
Solve described penalty coefficient matrix, obtain with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i, wherein N is the data total length of single passage.
Utilize penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1,0≤n 0≤ m-1, for n-th of the signalling channel data after compensation 0+ a data point.
The present embodiment also provides a kind of adaptive equalization system being applied to environmental magnetic field and suppressing, this system can realize the self-adapting compensation method being applied to environmental magnetic field and suppressing of the present invention, but this implement device being applied to the self-adapting compensation method that environmental magnetic field suppresses includes but not limited to bucking-out system provided by the invention.
As shown in Figure 2, the adaptive equalization system being applied to environmental magnetic field suppression described in comprises: data acquisition module 210, slip scan module 220, error function constructing module 230, penalty coefficient acquisition module 240, signal compensation module 250.
Described data acquisition module 210 gathers the signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports.Described detection system includes but not limited to biological magnetic detection system, and described biological magnetic detection system is provided with signalling channel and reference channel.Described data acquisition module 210 includes but not limited to A/D analog input card.
Described slip scan module 220 is connected with described data acquisition module 210, and the window utilizing length to be n is from initial time point n 0start to carry out slip scan extraction to described reference channel data, each distance of sliding is 1 data point, and obtain m compensation channels data successively, the product of m and n is less than or equal to the number of data points of described reference channel data.The length n of described window immobilizes, initial time point n 0increase gradually for gradient with 1 from 0.When the product of m and n equals the total length of all reference channel data, be the situation that self-adaptation determines coefficient.That is, self-adapting compensation method of the present invention is expanded, length of window is set to total data length, then can obtain single penalty coefficient, be self-adaptation and determine coefficient method.
Described error function constructing module 230 is connected with slip scan module 220 respectively with described data acquisition module 210, and the error function constructed is wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; 1≤i≤m;
Described penalty coefficient acquisition module 240 is connected with described error function constructing module 230, is that criterion solves described error function acquisition penalty coefficient k with minimum variance i.
Further, described penalty coefficient acquisition module 240 comprises: condition setting module 241, Matrix Solving module 242.Condition setting module 241 is connected with described error function constructing module 230, will according to minimum variance principle be set as 0, obtain about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ;
Matrix Solving module 242 is connected with described condition setting module 241, solves described penalty coefficient matrix, obtains with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i.
Described signal compensation module 250 is connected with described penalty coefficient acquisition module 240, utilizes penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1, for n-th of the signalling channel data after compensation 0+ a data point.Work as n 0during+n> total amount of data, namely stop compensating.
The present invention is based on signalling channel and the reference channel data of detection system, construct penalty function, set certain window length and slide in time, calculating the penalty coefficient in each window length, utilize penalty coefficient to realize high performance gradient and export.Specifically from initial point n 0start the data of getting stationary window length n, be that coefficient (i.e. penalty coefficient) is subtracted each other in criterion calculating with minimum variance, pointwise slides backward window and (namely fixes n, adjustment n 0), calculate and subtract each other coefficient k in each window i, can obtain and change and the penalty coefficient k of adaptive change in time i, thus acquisition changes the output of high performance gradient in time, realizes the squelch performance optimized.Feature of the present invention is the mode calculation compensation coefficient utilizing window sliding, and its advantage only has unique optimum configurations, implements simple, and can obtain the penalty coefficient of adaptive stabilizing.
The present invention compared with prior art, only has a parameter n needing to arrange, subtracts each other coefficient k without the need to setting iinitial value, there is not the process of convergence, namely Absorbable organic halogens obtains adaptively and best subtracts each other coefficient k i, realize high performance gradient and export.
In sum, the present invention effectively overcomes various shortcoming of the prior art and tool high industrial utilization.
Above-described embodiment is illustrative principle of the present invention and effect thereof only, but not for limiting the present invention.Any person skilled in the art scholar all without prejudice under spirit of the present invention and category, can modify above-described embodiment or changes.Therefore, such as have in art usually know the knowledgeable do not depart from complete under disclosed spirit and technological thought all equivalence modify or change, must be contained by claim of the present invention.

Claims (9)

1. be applied to environmental magnetic field suppress a self-adapting compensation method, it is characterized in that, described in be applied to environmental magnetic field suppress self-adapting compensation method comprise:
Gather signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports;
The window utilizing length to be n is from initial time point n 0start to extract described reference channel data, along with initial time point n 0from 0, be increased to N-n+1 gradually with 1 for gradient, obtain m compensation channels data successively, wherein N is the data total length of single passage;
Instrument error function wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; Be that criterion solves error function acquisition penalty coefficient k with minimum variance i, 1≤i≤m;
Utilize penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1,0≤n 0≤ m-1, for n-th of the signalling channel data after compensation 0+ a data point.
2. the self-adapting compensation method being applied to environmental magnetic field and suppressing according to claim 1, is characterized in that: the length n of described window immobilizes, and n is lowest mean square data length used.
3. the self-adapting compensation method being applied to environmental magnetic field and suppressing according to claim 2, is characterized in that, is that criterion solves error function acquisition penalty coefficient k with minimum variance idetailed process comprise:
When error function f reaches minimum, noise suppression effect is best, then have
∂ f ∂ k i = 0 ,
Can obtain thus about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ,
Solve described penalty coefficient matrix, obtain with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i.
4. the self-adapting compensation method being applied to environmental magnetic field and suppressing according to claim 1, it is characterized in that: described detection system comprises biological magnetic detection system, described biological magnetic detection system is provided with signalling channel and reference channel.
5. the self-adapting compensation method being applied to environmental magnetic field and suppressing according to claim 1, is characterized in that: when the product of m and n equals the total length of all reference channel data, is the situation that self-adaptation determines coefficient.
6. be applied to environmental magnetic field suppress an adaptive equalization system, it is characterized in that, described in be applied to environmental magnetic field suppress adaptive equalization system comprise:
Data acquisition module, gathers signalling channel signal and reference channel signal acquisition signalling channel data and reference channel data that a detection system exports;
Slip scan module, is connected with described data acquisition module, and the window utilizing length to be n is from initial time point n 0start to carry out slip scan extraction to described reference channel data, each distance of sliding is 1 data point, and obtain m compensation channels data successively, the product of m and n is less than or equal to the number of data points of described reference channel data;
Error function constructing module, is connected with slip scan module respectively with described data acquisition module, and the error function constructed is wherein, V sjfor a jth data point of signalling channel data, V rijbe a jth data point of i-th compensation channels data, k ibe penalty coefficients corresponding to i-th compensation channels data; 1≤i≤m;
Penalty coefficient acquisition module, is connected with described error function constructing module, is that criterion solves described error function acquisition penalty coefficient k with minimum variance i;
Signal compensation module, is connected with described penalty coefficient acquisition module, utilizes penalty coefficient k ito described signalling channel data V scarry out neighbourhood noise suppression, obtain the signalling channel data after compensating 0≤a≤n-1, for n-th of the signalling channel data after compensation 0+ a data point.
7. the adaptive equalization system being applied to environmental magnetic field and suppressing according to claim 6, is characterized in that: the length n of described window immobilizes, initial time point n 0increase gradually for gradient with 1 from 0, n is lowest mean square data length used.
8. the adaptive equalization system being applied to environmental magnetic field and suppressing according to claim 7, it is characterized in that, described penalty coefficient acquisition module comprises:
Condition setting module, is connected with described error function constructing module, will according to minimum variance principle be set as 0, obtain about penalty coefficient k imatrix
Σ j = n 0 n + n 0 - 1 V r 1 j 2 Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 1 j · V r 2 j Σ j = n 0 n + n 0 - 1 V r 2 j 2 . . . Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . . . . . . . Σ j = n 0 n + n 0 - 1 V r 1 j · V rmj Σ j = n 0 n + n 0 - 1 V r 2 j · V rmj . . . Σ j = n 0 n + n 0 - 1 V rmj 2 k 1 k 2 . . . k m = Σ j = n 0 n + n 0 - 1 V r 1 j · V sj Σ j = n 0 n + n 0 - 1 V r 2 j · V sj . . . Σ j = n 0 n + n 0 - 1 V rmj · V sj ;
Matrix Solving module, is connected with described condition setting module, solves described penalty coefficient matrix, obtains with n 0for starting point, data length is the penalty coefficient k of n i; Data length n fixes, at n 0from 0 with 1 for gradient be increased to N-n+1 gradually time, time dependent adaptive equalization coefficient k can be obtained i, wherein N is the data total length of single passage.
9. the adaptive equalization system being applied to environmental magnetic field and suppressing according to claim 6, it is characterized in that: described detection system comprises biological magnetic detection system, described biological magnetic detection system is provided with signalling channel and reference channel.
CN201310492564.7A 2013-10-18 2013-10-18 Adaptive compensation method and system for environmental magnetic field suppression Active CN104545913B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310492564.7A CN104545913B (en) 2013-10-18 2013-10-18 Adaptive compensation method and system for environmental magnetic field suppression

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310492564.7A CN104545913B (en) 2013-10-18 2013-10-18 Adaptive compensation method and system for environmental magnetic field suppression

Publications (2)

Publication Number Publication Date
CN104545913A true CN104545913A (en) 2015-04-29
CN104545913B CN104545913B (en) 2017-05-03

Family

ID=53063769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310492564.7A Active CN104545913B (en) 2013-10-18 2013-10-18 Adaptive compensation method and system for environmental magnetic field suppression

Country Status (1)

Country Link
CN (1) CN104545913B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105640538A (en) * 2016-04-01 2016-06-08 中国科学院上海微系统与信息技术研究所 Full-tensor magnetocardiography probe and producing method thereof
CN106343999A (en) * 2016-10-10 2017-01-25 中国科学院上海微系统与信息技术研究所 Magnetocardiography, compensation and optimization method based on same, system and server
CN111413651A (en) * 2020-03-30 2020-07-14 中国科学院上海微系统与信息技术研究所 Compensation method, device and system for total magnetic field and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4837514A (en) * 1988-07-14 1989-06-06 Atlantic Richfield Company Method of reducing noise in electromagnetic geophysical exploration data
CN101882964A (en) * 2010-06-12 2010-11-10 桂林电子科技大学 De-noising method of transient electromagnetic detecting echo signal
CN102988038A (en) * 2012-12-14 2013-03-27 中国科学院上海微系统与信息技术研究所 First-order gradient compensation module and method for unmask magnetocardiography

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4837514A (en) * 1988-07-14 1989-06-06 Atlantic Richfield Company Method of reducing noise in electromagnetic geophysical exploration data
CN101882964A (en) * 2010-06-12 2010-11-10 桂林电子科技大学 De-noising method of transient electromagnetic detecting echo signal
CN102988038A (en) * 2012-12-14 2013-03-27 中国科学院上海微系统与信息技术研究所 First-order gradient compensation module and method for unmask magnetocardiography

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
陈亮等: "超导梯度计及其对环境噪声的抑制能力", 《功能材料与器件学报》 *
黄玉等: "基于FLANN和最小二乘的磁梯度计误差校正", 《仪器仪表学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105640538A (en) * 2016-04-01 2016-06-08 中国科学院上海微系统与信息技术研究所 Full-tensor magnetocardiography probe and producing method thereof
CN106343999A (en) * 2016-10-10 2017-01-25 中国科学院上海微系统与信息技术研究所 Magnetocardiography, compensation and optimization method based on same, system and server
CN106343999B (en) * 2016-10-10 2019-04-19 中国科学院上海微系统与信息技术研究所 Magnetocardiograph, compensation optimizing method, system and server based on it
CN111413651A (en) * 2020-03-30 2020-07-14 中国科学院上海微系统与信息技术研究所 Compensation method, device and system for total magnetic field and storage medium
CN111413651B (en) * 2020-03-30 2021-04-13 中国科学院上海微系统与信息技术研究所 Compensation method, device and system for total magnetic field and storage medium

Also Published As

Publication number Publication date
CN104545913B (en) 2017-05-03

Similar Documents

Publication Publication Date Title
CN103190898B (en) Cardiac magnetic signal noise adaptive filtering and eliminating design method
Fang et al. Noise reduction in lidar signal based on discrete wavelet transform
CN104793253B (en) Aviation electromagnetic data de-noising method based on mathematical morphology
CN103837884B (en) Based on the digital core pulse signal trapezoidal shaping algorithm of time-domain analysis
CN103454621B (en) Method for denoising broadband radar target echoes based on matching pursuit
CN104545913A (en) Adaptive compensation method and system for environmental magnetic field suppression
CN106845334A (en) A kind of innovative noise extracting method based on mathematical morphology
CN103308829B (en) A kind of GIS single Partial discharge signal extracts and trigger instants method of adjustment
CN105807241A (en) Exponential signal denoising method achieved by means of prior information
CN104133199A (en) Weak-breathing-signal enhancement method used for life detection radar
CN104133248B (en) A kind of high fidelity sound wave interference drawing method
Kimiaefar et al. Random noise attenuation by Wiener-ANFIS filtering
CN106548780A (en) A kind of compressed sensing reconstructing method of voice signal
CN103064119B (en) A kind of fixing high sampling rate method for designing of time domain transient electromagnetic detection equipment
CN111697952B (en) Method and system for adjusting pulse width based on digital PZC system
CN105375992B (en) Based on gradient operator and the morphologic frequency spectrum cavity-pocket detection method of mathematics
CN111650655B (en) Non-negative matrix factorization supervised transient electromagnetic signal noise reduction method
CN102928872B (en) A kind of method improving dam, beach sand seismic reservoir recognition performance and describe precision
CN104132884B (en) The immediate processing method of a kind of signal base line in signal processing system and device
Shang et al. Chaotic SVD method for minimizing the effect of exponential trends in detrended fluctuation analysis
CN104007477B (en) Ground nuclear magnetic resonance three-dimensional inversion method
CN102624660A (en) Four-weighted fractional Fourier transform-based narrowband interference suppression method
Zhou et al. A data-dependent Fourier filter based on image segmentation for random seismic noise attenuation
CN104462679A (en) Design method of digital filter shaping circuit of radiation detection front-end readout system
CN102799757A (en) Weak signal extraction method for removing interferences of strong trend term and transient-state pulse

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180912

Address after: 201108 221, room 1, 3688 Kim do road, Minhang District, Shanghai.

Patentee after: Man Di medical instruments (Shanghai) Co., Ltd.

Address before: 200050 No. 865, Changning Road, Shanghai, Changning District

Patentee before: Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences