CN106725415B - Method and device for processing electrophysiological signals - Google Patents
Method and device for processing electrophysiological signals Download PDFInfo
- Publication number
- CN106725415B CN106725415B CN201611029763.4A CN201611029763A CN106725415B CN 106725415 B CN106725415 B CN 106725415B CN 201611029763 A CN201611029763 A CN 201611029763A CN 106725415 B CN106725415 B CN 106725415B
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- electricity physiological
- physiological signal
- electricity
- 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
- 238000000034 method Methods 0.000 title abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 51
- 238000003672 processing method Methods 0.000 claims abstract description 12
- 230000005611 electricity Effects 0.000 claims description 106
- 230000003595 spectral effect Effects 0.000 claims description 56
- 238000001228 spectrum Methods 0.000 abstract description 8
- 230000007613 environmental effect Effects 0.000 abstract 3
- 239000008280 blood Substances 0.000 description 27
- 210000004369 blood Anatomy 0.000 description 27
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical group [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 23
- 229910052760 oxygen Inorganic materials 0.000 description 23
- 239000001301 oxygen Substances 0.000 description 23
- 238000001514 detection method Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 108010054147 Hemoglobins Proteins 0.000 description 5
- 102000001554 Hemoglobins Human genes 0.000 description 5
- 108010064719 Oxyhemoglobins Proteins 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 125000002915 carbonyl group Chemical group [*:2]C([*:1])=O 0.000 description 3
- 102000010445 Lactoferrin Human genes 0.000 description 2
- 108010063045 Lactoferrin Proteins 0.000 description 2
- 108010061951 Methemoglobin Proteins 0.000 description 2
- 230000007831 electrophysiology Effects 0.000 description 2
- 238000002001 electrophysiology Methods 0.000 description 2
- CSSYQJWUGATIHM-IKGCZBKSSA-N l-phenylalanyl-l-lysyl-l-cysteinyl-l-arginyl-l-arginyl-l-tryptophyl-l-glutaminyl-l-tryptophyl-l-arginyl-l-methionyl-l-lysyl-l-lysyl-l-leucylglycyl-l-alanyl-l-prolyl-l-seryl-l-isoleucyl-l-threonyl-l-cysteinyl-l-valyl-l-arginyl-l-arginyl-l-alanyl-l-phenylal Chemical compound C([C@H](N)C(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CS)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](CC=1C2=CC=CC=C2NC=1)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CC=1C2=CC=CC=C2NC=1)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CC(C)C)C(=O)NCC(=O)N[C@@H](C)C(=O)N1CCC[C@H]1C(=O)N[C@@H](CO)C(=O)N[C@@H]([C@@H](C)CC)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CS)C(=O)N[C@@H](C(C)C)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](C)C(=O)N[C@@H](CC=1C=CC=CC=1)C(O)=O)C1=CC=CC=C1 CSSYQJWUGATIHM-IKGCZBKSSA-N 0.000 description 2
- 229940078795 lactoferrin Drugs 0.000 description 2
- 235000021242 lactoferrin Nutrition 0.000 description 2
- 208000004003 siderosis Diseases 0.000 description 2
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- CSJDCSCTVDEHRN-UHFFFAOYSA-N methane;molecular oxygen Chemical compound C.O=O CSJDCSCTVDEHRN-UHFFFAOYSA-N 0.000 description 1
- -1 oxyhemoglobin Chemical compound 0.000 description 1
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/316—Modalities, i.e. specific diagnostic methods
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The embodiment of the invention discloses a method and a device for processing electrophysiological signals, wherein the processing method comprises the following steps: acquiring an electrophysiological signal and an environmental electrical signal; determining the frequency range of the frequency of a power frequency interference signal generated by the environmental electric signal to the electrophysiological signal; sampling the electrophysiological signal according to a sampling rate that lies within the frequency range; determining the frequency of the spectrum aliasing signal corresponding to the electrophysiological signal and the environmental electric signal according to the sampling rate in the frequency range and the frequency of the power frequency interference signal; removing the spectrally aliased signal from the sampled electrophysiological signal in dependence on the frequency of the spectrally aliased signal. The removal of the power frequency interference in the electrophysiological signal is realized under the condition of not increasing hardware resources.
Description
Technical field
The present embodiments relate to signal processing technology field, a kind of processing method more particularly to electricity physiological signal and
Device.
Background technique
Hz noise refers to that frequency is the interference that the line voltage of 50Hz radiates formation in the form of an electromagnetic wave.Due to user
Acquisition Circuit and power supply circuit etc. cannot cause power supply to form effective pressure difference to Acquisition Circuit, therefore well altogether in equipment
Electric power signal enters Acquisition Circuit and produces interference.Hz noise can interfere electrical equipment and electronic equipment, cause
Equipment is operating abnormally.
Due to the problem of being frequently encountered removal Hz noise in the measurement of electricity physiological signal.Hz noise will affect letter
Number detection, cause electricity physiological signal calculate inaccuracy.The minimizing technology of usual Hz noise is mostly from the angle of Stochastic signal processing
Degree sets out, although effect is pretty good, often there is algorithm complexity height, it is difficult to the disadvantages of understanding.In addition, there are also Hz noises
Minimizing technology be to improve 2 times or more of the signal sampling rate to power frequency interference signals frequency, application sample theorem passes through design work
Frequency trapper is removed, but can occupy a large amount of hardware resource.
Summary of the invention
In view of this, the present invention proposes a kind for the treatment of method and apparatus of electricity physiological signal, realize in electricity physiological signal
The removal of middle Hz noise, resource occupation is few, and design is simple.
In a first aspect, the treating method comprises the embodiment of the invention provides a kind of processing method of electricity physiological signal:
Obtain electricity physiological signal and environment electric signal；Determine the environment electric signal to Hz noise caused by the electricity physiological signal
Frequency range where the frequency of signal；The electricity physiological signal is adopted according to the sample rate being located in the frequency range
Sample；According to the frequency of sample rate and the power frequency interference signals in the frequency range determine the electricity physiological signal with it is described
The frequency of the corresponding spectral aliasing signal of environment electric signal；According to the frequency of the spectral aliasing signal from the electro physiology after sampling
The spectral aliasing signal is removed in signal.
Further, determine the environment electric signal to the frequency of power frequency interference signals caused by the electricity physiological signal
The frequency range at place includes: to calculate power frequency interference signals according to the frequency of power frequency component and the frequency of the environment electric signal
Frequency range where frequency.
Further, the frequency spectrum is removed from the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal
Aliasing signal includes: to carry out trap to the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal.
Further, the frequency of the electricity physiological signal and the signal frequency after the spectral aliasing are inconsistent.
Second aspect, the embodiment of the invention provides a kind of processing unit of electricity physiological signal, described device includes: to obtain
Unit, for obtaining electricity physiological signal and environment electric signal；Determination unit is connected with the acquiring unit, determines the environment
Electric signal is to the frequency range where the frequency of power frequency interference signals caused by the electricity physiological signal；Sampling unit, with institute
It states determination unit to be connected, for sampling the electricity physiological signal according to the sample rate being located in the frequency range；Meter
Unit is calculated, is connected with the sampling unit, for according to the sample rate and the power frequency interference signals in the frequency range
Frequency determines the frequency of electricity physiological signal spectral aliasing signal corresponding with the environment electric signal；Processing unit, with institute
It states computing unit to be connected, for removing the frequency from the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal
Compose aliasing signal.
Further, the determination unit is specifically used for: according to the frequency of the frequency of power frequency component and the environment electric signal
Rate calculates the frequency range where the frequency of power frequency interference signals.
Further, the processing unit is specifically used for: according to the frequency of the spectral aliasing signal to the electricity after sampling
Physiological signal carries out trap.
Further, the frequency of the electricity physiological signal and the signal frequency after the spectral aliasing are inconsistent.
In the embodiment of the present invention, by calculating the frequency of power frequency interference signals in electricity physiological signal, and Hz noise is calculated
Frequency range where the frequency of signal is sampled according to sample rate within the scope of this, in the electricity physiological signal after removal sampling
Spectral aliasing signal, the spectral aliasing signal be sample after power frequency interference signals.It realizes to work in electricity physiological signal
The removal of frequency interference signal, design is simple, and resources occupation is few.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, of the invention other
Feature, objects and advantages will become more apparent upon:
Fig. 1 is the flow chart of one of embodiment of the present invention one electricity physiological signal processing method；
Fig. 2 a is the flow chart of one of embodiment of the present invention two electricity physiological signal processing method；
Original signal figure when Fig. 2 b is the 80Hz sample rate in the embodiment of the present invention two；
Signal spectrum figure when Fig. 2 c is the 80Hz sample rate in the embodiment of the present invention two；
Original signal figure when Fig. 2 d is the 250Hz sample rate in the embodiment of the present invention two；
Signal spectrum figure when Fig. 2 e is the 250Hz sample rate in the embodiment of the present invention two；
Fig. 3 is the structure chart of one of embodiment of the present invention three electricity physiological signal processing unit.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It should also be noted that, for the ease of retouching
It states, only some but not all contents related to the present invention are shown in the drawings.It also should be noted that for the ease of saying
It is bright, example related to the present invention is shown in following embodiment, these examples are only as the principle for illustrating the embodiment of the present invention
It is used, it is not intended as the restriction to the embodiment of the present invention, meanwhile, these exemplary specific values can be according to different application environments
It is different and different with the parameter of device or component.
The treating method and apparatus of the electricity physiological signal of the embodiment of the present invention, which can be run on, is equipped with Windows (Microsoft
Company exploitation operating system platform), Android (Google exploitation the operation for Portable movable smart machine
System platform), the iOS operating system platform for Portable movable smart machine of exploitation (Apple Inc.), Windows
The terminal of the operating systems such as Phone (operating system platform for Portable movable smart machine of Microsoft's exploitation)
In, which can be desktop computer, laptop, mobile phone, palm PC, tablet computer, digital camera, digital vedio recording
Any one in machine etc..
Embodiment one
Fig. 1 is the flow chart of the processing method of one of the embodiment of the present invention one electricity physiological signal, and the present embodiment can fit
For electricity physiological signal situation to be treated, this method can be executed by electricity physiological signal processing unit, which can lead to
The mode for crossing hardware and/or software realizes that the device is typically configured in the equipment for providing image processing function.Specifically, this
The processing method of electricity physiological signal includes: step S110, step S120, step S130, step S140 and step in embodiment
S150。
Step S110 obtains electricity physiological signal and environment electric signal.
Specifically, the electricity physiological signal that acquisition is to be processed, measures its frequency, frequency is denoted as H1.Electricity physiological signal is reflection
A kind of important information of the various physiological status of human body, electricity physiological signal include but is not limited to electrocardio, myoelectricity and electro-ocular signal.Electricity is raw
Reason signal is the signal of interest of underlying disease information, the complexity with time domain.Environment electric signal is obtained, its frequency is measured, is remembered
For H2.Environment electric signal refers to, the electric signal of the other equipment in environment locating for electrophysiologicalsignal signal detection equipment, the electric signal
Electricity physiological signal is interfered.Optionally, electricity physiological signal includes but is not limited to blood oxygen signal, electrophysiologicalsignal signal detection equipment
Including but not limited to blood oxygen measuring device, environment electric signal include but is not limited to the incandescent modulating signal in environment.
Step S120 determines the environment electric signal to the frequency of power frequency interference signals caused by the electricity physiological signal
The frequency range at place.
Specifically, power frequency refers generally to the frequency of alternating current, it is denoted as H3, it is 50Hz in China, also there are 60Hz in other countries.Environment
Electric signal causes environment electric signal to produce power frequency interference signals under the influence of power frequency component by mains-supplied, and the power frequency is dry
The frequency for disturbing signal is denoted as H4, which can be to the measurement generation interference in electricity physiological signal in environment.According to
The frequency H of environment electric signal2, calculate the frequency H of the power frequency interference signals of generation4, and determine the frequency where power frequency interference signals
Range H5~H6。
Step S130 samples the electricity physiological signal according to the sample rate being located in the frequency range.
Specifically, according to the frequency range H of Hz noise5~H6Determine sample rate H0, choose sample rate H0It is dry positioned at power frequency
The frequency range H disturbed5~H6It is interior.According to sample rate H0Electricity physiological signal is sampled.Sampling is the sample for acquiring analog signal
This, i.e., will the time is upper, all continuous analog signal in amplitude, in the effect of sampling pulse, it is upper discrete (on the time to be converted into the time
Have fixed intervals) but amplitude on still continuous discrete analog signal.So sampling is also known as the discretization process of waveform.It is per second
The sample number of clock is called sample rate, and sample rate is higher, and the fidelity of signal is higher, but quantify after information content amount of storage
It is bigger.According to sample rate H0Electricity physiological signal is sampled, i.e., by electricity physiological signal from digital signal be converted on the time from
It dissipates, continuous discrete analog signal in amplitude.
Step S140 determines the electricity according to the frequency of sample rate and the power frequency interference signals in the frequency range
The frequency of physiological signal spectral aliasing signal corresponding with the environment electric signal.
Specifically, according to sample rate H0With the frequency H of power frequency interference signals4, calculate electricity physiological signal and environment electric signal pair
Frequency H after the spectral aliasing answered7.Frequency spectrum refers to representation of the signal of a ken under frequency domain, and frequency spectrum can indicate
One signal is made of the string wave of which frequency, it is also seen that the information such as size and phase of each frequency string wave.When adopting
When sample rate is lower than twice of Hz noise, spectral aliasing occurs.The case where spectral aliasing refers to, adjacent frequency spectrum stacks,
Step S150 removes the frequency spectrum according to the frequency of the spectral aliasing signal from the electricity physiological signal after sampling
Aliasing signal.
Specifically, sampling after electricity physiological signal in spectral aliasing signal be power frequency interference signals after over-sampling
Signal, the spectral aliasing signal in electricity physiological signal according to the frequency of the spectral aliasing signal, after removal sampling.
In the embodiment of the present invention, by calculating the frequency of power frequency interference signals in electricity physiological signal, and Hz noise is calculated
Frequency range where the frequency of signal is sampled according to sample rate within the scope of this, in the electricity physiological signal after removal sampling
Spectral aliasing signal, the spectral aliasing signal be sample after power frequency interference signals.It realizes to work in electricity physiological signal
The removal of frequency interference signal, design is simple, and resources occupation is few.
Embodiment two
Fig. 2 a is a kind of flow chart of the processing method of electricity physiological signal in the embodiment of the present invention two.The present embodiment is in reality
On the basis of applying example one, step S120 and step S150 are further illustrated.In step S120, the environment electric signal pair is determined
Frequency range where the frequency of power frequency interference signals caused by the electricity physiological signal includes: the frequency according to power frequency component
The frequency range where the frequency of power frequency interference signals is calculated with the frequency of the environment electric signal.In step S150, according to institute
It includes: according to the frequency that the frequency for stating spectral aliasing signal removes the spectral aliasing signal from the electricity physiological signal after sampling
The frequency for composing aliasing signal carries out trap to the electricity physiological signal after sampling.Specifically, in the present embodiment electricity physiological signal place
Reason method includes: step S210, step S220, step S230, step S240 and step S250.
Step S210 obtains electricity physiological signal and environment electric signal.
Step S220 calculates the frequency of power frequency interference signals according to the frequency of power frequency component and the frequency of the environment electric signal
Frequency range where rate.
Specifically, according to the frequency H of power frequency component3, the frequency H of environment electric signal2, consider a certain specific environment telecommunications
Number frequency, determine how the frequency of power frequency component influences environment electric signal, according to the influence relationship calculate power frequency interference signals
Frequency H4, and determine the frequency range H where the frequency of the power frequency interference signals5~H6。
Step S230 samples the electricity physiological signal according to the sample rate being located in the frequency range.
Step S240 determines the electricity according to the frequency of sample rate and the power frequency interference signals in the frequency range
The frequency of physiological signal spectral aliasing signal corresponding with the environment electric signal.
Step S250 carries out trap to the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal, from adopting
The spectral aliasing signal is removed in electricity physiological signal after sample.
Specifically, design trapper carries out trap to the electricity physiological signal after sampling.Trapper is also known as bandstop filter,
It is a kind of resonance circuit, or perhaps a kind of inductor of automatic switch, it can be according to signal using it on antenna works
Frequency automatically prolongs or shortens the length of antenna.Dedicated for eliminating certain garbage signals to reduce pair in radio receiver
The filter of the interference of useful signal is usually used in removing the interference signal of a certain specific frequency, apply anti-interference in telecommunication signal
In.The trapper of spectral aliasing signal frequency point is designed, realizes the removal to the power frequency interference signals in electricity physiological signal.
Optionally, in the equipment of detection electricity physiological signal, original equipment has default want during use, to sample rate
It asks.According to the requirement in former design to sample rate, liter sampling or down-sampled is carried out to the signal by trapper trap after
Meet requirement of the former design to sample rate.The sample rate calculated when if removing Hz noise is higher than intrinsic sample rate, right
The later signal of trap carries out down-sampled；The sample rate calculated when if removing Hz noise is right lower than intrinsic sample rate
The later signal of trap carries out a liter sampling.
Optionally, the frequency of the electricity physiological signal and the signal frequency after the spectral aliasing are inconsistent.
Specifically, the frequency of electricity physiological signal is typically small, with the signal frequency after spectral aliasing not in a frequency range, because
This will not have an impact original signal.Optionally, when electricity physiological signal is blood oxygen signal, frequency is 0.3HZ~5HZ,
Under the interference of environment light source incandescent lamp, the signal frequency after spectral aliasing is 30Hz/20Hz/20Hz/30Hz.After spectral aliasing
Signal frequency will not have an impact blood oxygen signal in the frequency of removal power frequency interference signals far from electricity physiological signal.
The embodiment of the present invention calculates power frequency interference signals according to the frequency of power frequency component and the frequency of the environment electric signal
Frequency, and determine the frequency where frequency range.Trapper is designed, to the frequency of the Frequency point of the power frequency interference signals after sampling
Rate carries out trap, the power frequency interference signals after removing spectral aliasing.By the design of trapper, electricity physiological signal intermediate frequency is realized
The removal of power frequency interference signals after composing aliasing.
In order to be more readily understood the statement of scheme, based on embodiment one and embodiment two, optionally, electro physiology is believed
It number is that blood oxygen signal, scheme when being incandescent lamp of blood oxygen measuring device and environment electric signal of detection device of electricity physiological signal carry out
It illustrates.It is worth noting that, following parameter value is a kind of exemplary values, and there may also be other values, specific number
Value is different, and the detection accuracy of blood oxygen signal is different, can according to circumstances carry out related adjustment.
Blood usually contains a plurality of types of hemoglobins, including oxyhemoglobin, reduced hemoglobin, carbon oxygen are blood red
Albumen and ferrihemoglobin.Wherein, comparision contents are low (1%) under normal circumstances for carbonyl haemoglobin and ferrihemoglobin.
Blood oxygen saturation is often referred to the percentage that oxyhemoglobin in blood accounts for hemoglobin, due to carbonyl haemoglobin and siderosis
Comparision contents are low under normal circumstances for Lactoferrin, so often ignoring carbonyl haemoglobin and siderosis when calculating blood oxygen saturation
The content of Lactoferrin.
The method of measurement blood oxygen saturation is to irradiate finger by feux rouges (660nm) and infrared light (940nm) at present, by
It is smaller than reduced hemoglobin absorption coefficient under feux rouges in oxyhemoglobin, system is absorbed than reduced hemoglobin under infrared light
Number is big, and when oxyhemoglobin changes of contents in blood of human body, feux rouges and infrared signal amplitude can change, and utilize this
A variation can calculate detecting sphygmus and blood oxygen saturation, be a kind of not damaged blood oxygen measuring method.Ring as where measuring device
The artificial light sources in border will affect signal detection there are reasons such as Hz noises, and oximeter is caused not calculated accurately really, influence examining for the state of an illness
Disconnected and treatment.Hz noise is removed in the case of measuring blood oxygen saturation using blood oxygen measuring device.
Specifically, the frequency H of blood oxygen signal1For 0.3Hz~5Hz, power frequency H3For 50Hz/60Hz, locating for blood oxygen measurement device
Incandescent lamp in environment is interference, and incandescent lamp has light and shade variation twice within a power frequency component period, so the work generated
Frequency interference H4For 100Hz/120Hz.Using spectral aliasing effect, the method for reducing sample rate is taken, rather than improves sample rate,
Frequency of the Hz noise frequency far from blood oxygen signal after making aliasing.In order to allow after 50Hz/60Hz and 100Hz/120Hz aliasing
Frequency is simultaneously far from blood oxygen signal frequency 0.3Hz~5Hz, near the tetra- frequency midpoints 50Hz/60Hz and 100Hz/120Hz
Frequency H5~H6(75Hz~90Hz) is used as sample rate.It is desirable can be by the sampling of frequency 80Hz that 1 second eliminates as the condition that meets
Rate, at this point, the frequency after 50Hz/60Hz and 100Hz/120Hz aliasing is respectively 30Hz/20Hz/20Hz/40Hz.Design trap
Device, separately designs 30Hz, and the digital trap of tri- Frequency points of 20Hz, 40Hz is removed the Hz noise after aliasing, together
When since three Frequency points are all far from the frequency of blood oxygen signal, so will not be had an impact to blood oxygen signal.Then it is set according to original
The requirement of meter carries out liter sampling or the down-sampled requirement with the former design of satisfaction to sample rate to the signal after trap.If
The sample rate of former design requirement is greater than 80Hz, then carries out a liter sampling, if the sample rate of former design requirement is less than 80Hz, carries out
It is down-sampled.By the removal to application blood oxygen measurement device to the power frequency interference signals in blood oxygen signal measurement, realize not
Under conditions of increasing hardware resource, it is removed using the power frequency interference signals that spectral aliasing effect is gone in blood oxygen signal.
Optionally, original signal figure when Fig. 2 b is 80Hz sample rate, Fig. 2 c are 80Hz sampling rate signal spectrograms, and Fig. 2 d is
Original signal figure when 250Hz sample rate, Fig. 2 e are 250Hz sampling rate signal spectrograms.Curve 260 in Fig. 2 b is sample rate
When 80Hz, original signal amplitude curve；When curve 270 in Fig. 2 c is sample rate 80Hz, sampling rate signal spectrogram.100Hz
Light interference the spectral aliasing after over-sampling be 20Hz, from Fig. 2 c as can be seen that original signal in interference signal quilt
It efficiently separates, to design the trapper removal Hz noise of respective frequencies point.Curve 280 in Fig. 2 d is sample rate 250Hz
When, original signal amplitude curve, when the curve 290 in Fig. 2 e is sample rate 250Hz, sampling rate signal spectrogram.Due to sampling
Rate 250Hz is greater than the 2 times or more of power frequency interference signals, and spectral aliasing does not occur.It can be seen that by Fig. 2 e when sample rate is
When 250Hz, it is also able to achieve the removal to Hz noise, but high sampling rate can occupy a large amount of hardware resource, design is complicated.
Embodiment three
Fig. 3 is the structure chart of the processing unit of one of the embodiment of the present invention three electricity physiological signal.The device is suitable for
The processing method of the electricity physiological signal provided in the embodiment of the present invention one and the embodiment of the present invention two is executed, which specifically wraps
It includes: acquiring unit 310, determination unit 320, sampling unit 330, computing unit 340 and processing unit 350.
Acquiring unit 310, for obtaining electricity physiological signal and environment electric signal.
Determination unit 320 is connected with acquiring unit 310, determines that the environment electric signal produces the electricity physiological signal
Frequency range where the frequency of raw power frequency interference signals.
Sampling unit 330 is connected with determination unit 320, is located at the sample rate in the frequency range to institute for basis
Electricity physiological signal is stated to be sampled.
Computing unit 340 is connected with sampling unit 330, for according to the sample rate and the work in the frequency range
The frequency of frequency interference signal determines the frequency of electricity physiological signal spectral aliasing signal corresponding with the environment electric signal.
Processing unit 350 is connected with computing unit 340, for according to the frequency of the spectral aliasing signal after sampling
Electricity physiological signal in remove the spectral aliasing signal.
Further, determination unit 320 is specifically used for: according to the frequency of the frequency of power frequency component and the environment electric signal
Calculate the frequency range where the frequency of power frequency interference signals.
Further, processing unit 350 is specifically used for: raw to the electricity after sampling according to the frequency of the spectral aliasing signal
It manages signal and carries out trap.
Further, the frequency of the electricity physiological signal and the signal frequency after the spectral aliasing are inconsistent.
In the embodiment of the present invention, by calculating the frequency of power frequency interference signals in electricity physiological signal, and Hz noise is calculated
Frequency range where the frequency of signal is sampled according to sample rate within the scope of this, in the electricity physiological signal after removal sampling
Spectral aliasing signal, the spectral aliasing signal be sample after power frequency interference signals.It realizes to work in electricity physiological signal
The removal of frequency interference signal, design is simple, and resources occupation is few.
Obviously, it will be understood by those skilled in the art that side provided by any embodiment of the invention can be performed in the said goods
Method has the corresponding functional module of execution method and beneficial effect.
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that
The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation,
It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention
It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also
It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.
Claims (6)
1. a kind of processing method of electricity physiological signal characterized by comprising
Obtain electricity physiological signal and environment electric signal；
Determine the environment electric signal to the frequency model where the frequency of power frequency interference signals caused by the electricity physiological signal
It encloses；
The electricity physiological signal is sampled according to the sample rate being located in the frequency range；
The electricity physiological signal and institute are determined according to the frequency of sample rate and the power frequency interference signals in the frequency range
State the frequency of the corresponding spectral aliasing signal of environment electric signal, the letter after the frequency of the electricity physiological signal and the spectral aliasing
Number frequency is not in the same frequency range；
The spectral aliasing signal is removed from the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal.
2. the processing method of electricity physiological signal according to claim 1, which is characterized in that determine the environment electric signal pair
Frequency range where the frequency of power frequency interference signals caused by the electricity physiological signal includes:
The frequency where the frequency of power frequency interference signals is calculated according to the frequency of power frequency component and the frequency of the environment electric signal
Range.
3. the processing method of electricity physiological signal according to claim 1, which is characterized in that according to the spectral aliasing signal
Frequency remove the spectral aliasing signal from the electricity physiological signal after sampling and include:
Trap is carried out to the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal.
4. a kind of processing unit of electricity physiological signal characterized by comprising
Acquiring unit, for obtaining electricity physiological signal and environment electric signal；
Determination unit is connected with the acquiring unit, determines the environment electric signal to work caused by the electricity physiological signal
Frequency range where the frequency of frequency interference signal；
Sampling unit is connected with the determination unit, for raw to the electricity according to the sample rate being located in the frequency range
Reason signal is sampled；
Computing unit is connected with the sampling unit, for according in the frequency range sample rate and the Hz noise
The frequency of signal determines the frequency of electricity physiological signal spectral aliasing signal corresponding with the environment electric signal, and the electricity is raw
The frequency of signal and the signal frequency after the spectral aliasing are managed not in the same frequency range；
Processing unit is connected with the computing unit, raw from the electricity after sampling for the frequency according to the spectral aliasing signal
The spectral aliasing signal is removed in reason signal.
5. the processing unit of electricity physiological signal according to claim 4, which is characterized in that the determination unit is specifically used
In:
The frequency where the frequency of power frequency interference signals is calculated according to the frequency of power frequency component and the frequency of the environment electric signal
Range.
6. the processing unit of electricity physiological signal according to claim 4, which is characterized in that the processing unit is specifically used
In:
Trap is carried out to the electricity physiological signal after sampling according to the frequency of the spectral aliasing signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611029763.4A CN106725415B (en) | 2016-11-15 | 2016-11-15 | Method and device for processing electrophysiological signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611029763.4A CN106725415B (en) | 2016-11-15 | 2016-11-15 | Method and device for processing electrophysiological signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106725415A CN106725415A (en) | 2017-05-31 |
| CN106725415B true CN106725415B (en) | 2019-10-18 |
Family
ID=58970379
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611029763.4A Active CN106725415B (en) | 2016-11-15 | 2016-11-15 | Method and device for processing electrophysiological signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106725415B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019144406A1 (en) * | 2018-01-29 | 2019-08-01 | 深圳迈瑞生物医疗电子股份有限公司 | Blood oxygen measurement system and method |
| CN109124610B (en) * | 2018-09-26 | 2021-07-30 | 深圳星脉医疗仪器有限公司 | Anti-interference method and device for non-invasive blood pressure measurement |
| CN110930973A (en) * | 2019-12-19 | 2020-03-27 | 上海宏桐实业有限公司 | Electrophysiological signal noise control method |
| CN116602643B (en) * | 2023-07-20 | 2023-11-17 | 成都晨电智能科技有限公司 | Sampling method and circuit for power frequency interference resistant signal and double-electrode heart rate sensor |
| CN118707172A (en) * | 2024-08-28 | 2024-09-27 | 宁德时代新能源科技股份有限公司 | Signal testing system, signal testing method, signal testing device and oscilloscope |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1930588A (en) * | 2004-03-12 | 2007-03-14 | 皇家飞利浦电子股份有限公司 | Adaptive sampling along edges for surface rendering |
| CN101548885B (en) * | 2009-04-17 | 2010-12-08 | 南京大学 | Method for eliminating power frequency interfering signals in electrophysiological signals |
| US20110301439A1 (en) * | 2010-06-08 | 2011-12-08 | AliveUSA LLC | Wireless, ultrasonic personal health monitoring system |
| CN102940489B (en) * | 2012-11-29 | 2014-12-17 | 天津开发区奥金高新技术有限公司 | Filter designing method of weak electrophysiology signal and filtering method |
| US9521979B2 (en) * | 2013-03-15 | 2016-12-20 | Medtronic, Inc. | Control of spectral agressors in a physiological signal monitoring device |
| US20140355001A1 (en) * | 2013-05-28 | 2014-12-04 | Stratus Devices, Inc. | Measuring Deflection in an Optical Fiber Sensor by Comparing Current and Baseline Frames of Speckle Interference Patterns |
| CN103598886B (en) * | 2013-11-28 | 2015-10-28 | 中山大学 | Based on the method for dynamic urine volume in the long-pending urine process of model compensation method monitoring |
| CN106031633A (en) * | 2015-03-07 | 2016-10-19 | 上海贝瑞电子科技有限公司 | An electrocardiogram monitoring method and system |
| CN104783780B (en) * | 2015-04-13 | 2017-11-24 | 深圳市飞马与星月科技研究有限公司 | ECG De method and device |
| CN105266800B (en) * | 2015-12-02 | 2018-05-25 | 广东工业大学 | One kind is based on fetus electrocardio blind separation under Low SNR |
| CN106096579A (en) * | 2016-06-22 | 2016-11-09 | 天津理工大学 | A kind of method of ECG signal processing |
-
2016
- 2016-11-15 CN CN201611029763.4A patent/CN106725415B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106725415A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106725415B (en) | Method and device for processing electrophysiological signals | |
| Smital et al. | Real-time quality assessment of long-term ECG signals recorded by wearables in free-living conditions | |
| CN201104882Y (en) | Blood oxygen saturation measurement mechanism | |
| US9504401B2 (en) | Atrial fibrillation analyzer and program | |
| CN104602594B (en) | Pulse wave detection method and pulse wave detection device | |
| JP6052027B2 (en) | Pulse wave detection device, pulse wave detection program, and pulse wave detection method | |
| Zhang et al. | VLSI friendly ECG QRS complex detector for body sensor networks | |
| US9662023B2 (en) | Robust heart rate estimation | |
| JP2004514493A (en) | Method and apparatus for estimating physiological parameters from physiological signals | |
| US8295901B2 (en) | Spectral analysis for a more reliable determination of physiological parameters | |
| CN105997043B (en) | A kind of pulse frequency extracting method based on wrist wearable device | |
| Lai et al. | Lightweight wrist photoplethysmography for heavy exercise: motion robust heart rate monitoring algorithm | |
| US20160029968A1 (en) | Tracking slow varying frequency in a noisy environment and applications in healthcare | |
| WO2013040399A1 (en) | Method and apparatus for denoising of physiological signals | |
| Li et al. | Comparison and noise suppression of the transmitted and reflected photoplethysmography signals | |
| CN112494001B (en) | PPG signal quality evaluation method and device, wearable device and storage medium | |
| Kumar et al. | Reference signal less Fourier analysis based motion artifact removal algorithm for wearable photoplethysmography devices to estimate heart rate during physical exercises | |
| CN111820870A (en) | Biological image processing method and physiological information detection device | |
| KR102451623B1 (en) | Method and Apparatus for Comparing Features of ECG Signal with Difference Sampling Frequency and Filter Methods for Real-Time Measurement | |
| KR101786014B1 (en) | SPO2 measurement system based on wrist-type photoplethysmography and method thereof | |
| TWI504378B (en) | Denoising method and apparatus of pulse wave signal and pulse oximetry | |
| CN113397523B (en) | Respiration rate detection method and device and medical equipment | |
| US11534075B2 (en) | Heart activity monitoring during physical exercise | |
| CN105796051B (en) | Three-dimensional physiology-detecting system and its operating method | |
| KR102688647B1 (en) | Signal processing method, signal filtering apparatus and signal processing unit |
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 |