CN108113664A - Inhibit the method and system of notch filter ring in electrocardiosignal acquisition process - Google Patents
Inhibit the method and system of notch filter ring in electrocardiosignal acquisition process Download PDFInfo
- Publication number
- CN108113664A CN108113664A CN201711286389.0A CN201711286389A CN108113664A CN 108113664 A CN108113664 A CN 108113664A CN 201711286389 A CN201711286389 A CN 201711286389A CN 108113664 A CN108113664 A CN 108113664A
- Authority
- CN
- China
- Prior art keywords
- electrocardiosignal
- current sample
- sample time
- caching
- moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- 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/30—Input circuits therefor
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth 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)
- Cardiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
Inhibit the method and system of notch filter ring in a kind of electrocardiosignal acquisition process, electrocardiosignal is sampled in real time, obtain electrocardiosignal, the described method includes：Step is determined including caching step, the detection of jumbo jump signal and inhibition step, small size jump signal detecting step, trap step and caching moment.The method and system that the embodiment of the present invention proposes have the advantages that calculate efficient, realization simplicity, strong applicability.
Description
Technical field
The invention belongs to ECG's data compression technical fields, and in particular to inhibit trap in a kind of electrocardiosignal acquisition process
The method and system of wave filter ring.
Background technology
Heart can record heart in each cardiac cycle by electrocardiosignal (Electrocardiogram, ECG)
In succession excited by pacemaker, atrium, ventricle in each cardiac cycle, with bioelectric variation, bioelectric variation claims
For electrocardio, metamorphosis reflects human heart function, is the reference frame for reflecting heart state.
In the prior art, using electrocardioscanner gather human ecg signal, by filtering method remove baseline drift,
The interference signal that Hz noise, myoelectricity interference are adulterated when detections is obtained more clearly electrocardiosignal, is then carried using feature
The features such as phase, S ripples, T ripples, the shifting of ST sections, ST slope over 10 between method extraction electrocardiosignal QRS complex, RR are taken, to reflect cardiac function
Power.The signal amplitude scope of electrocardiosignal is about in 10 μ V-5mV, and frequency range is about in 0.05-100Hz.The typical heart
Electric signal waveform is as shown in Figure 1.
During cardiac pacing, original electro-cardiologic signals rising edge QR wave steepness is high and steep, there is mutation, as shown in Fig. 2, transverse axis represents in figure
Time, unit are millisecond ms, and the longitudinal axis is electrocardiosignal intensity, and unit is microvolt μ V.Due to the filter in ECG signal sampling instrument
The inherent characteristic of ripple device, jump signal near R ripple inflection points after notch filter by easily occurring ring, the generation collected
The electrocardiosignal of ring is as shown in figure 3, there is interference in it, it is impossible to normal reflection human heart function.
In order to inhibit to gather the ring in electrocardiosignal, generally use sef-adapting filter or related average method,
But these methods calculate complicated, computationally intensive, the real-time processing of unsuitable electrocardiosignal.
The content of the invention
In order to reduce operand, the adaptability for improving algorithm in electrocardiosignal acquisition process, the embodiment of the present invention proposes
A kind of electrocardiosignal inhibits notch filter ring method during using, includes the following steps：
Caching step：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1
The electrocardiosignal d [0] of trap output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the period
Length value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches the moment
ti-1Closest to current sample time ti；
Jumbo jump signal detection and inhibition step：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] is carried out
Jumbo jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], into small size jump signal detect walk
Suddenly；
Small size jump signal detecting step：It calculates in moment t1、...、ti-1、tiEach electrocardiosignal slope absolute value | k
[n]|：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval the heart
Change in electric rate threshold value L, and | k [i] | more than the electrocardiosignal change rate threshold value H of unit interval, then when judging present sample
Carve tiElectrocardiosignal d [i] there are small size mutation；
Trap step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first is used to it
Wave filter carries out trap output；
The caching moment determines step：If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then really
Determine current sample time tiBelong to the caching moment；
Return to step：Return cache step.
Further, in jumbo jump signal detects and inhibits step, current sample time t is obtainediOriginal electrocardiographicdigital letter
Number d0[i], previous caching moment ti-1Electrocardiosignal d [i-1], calculate current sample time electrocardiosignal slope k [i] be：
K [i]=(d0[i]-d[i-1])/(ti-ti-1),
If the absolute value of current sample time electrocardiosignal slope k [i] | k [i] | more than the heart of the unit interval of setting
Change in electric threshold value A, it is determined that the electrocardiosignal of current sample time be jumbo jump signal, current sample time tiThe heart
Electric signal d [i] is：
D [i]=d [i-1]+(k [i]/| k [i] |) * A；
If current sample time tiElectrocardiosignal slope k [i] absolute value be less than setting unit interval electrocardio
Signal intensity threshold value A, current sample time tiElectrocardiosignal d [i] be：D [i]=d0[i]。
Further, in jumbo jump signal detects and inhibits step, the electrocardiosignal change threshold A of every ms of setting exists
Value in the range of [800 μ V, 1200 μ V].
Further, if current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then in trap step
First filter is used to carry out trap output to it, wherein first filter is the wave filter without caching function.
Further, if current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then in trap step
Second filter is used to carry out trap output to it.
Further, after caching step, the method further includes time interval judgment step, obtains current sample time
tiWith previous caching moment ti-1Between time interval ti-ti-1If time interval ti-ti-1More than the time threshold of setting,
Then using the second notch filter to current sample time tiOriginal electro-cardiologic signals d0[i] carries out trap output, and determines current
Sampling instant tiBelong to the caching moment, return cache step carries out the ECG's data compression of subsequent time；If time interval ti-
ti-1No more than the time threshold of the setting, then detected into jumbo jump signal and inhibit step.
Further, in caching step, the time segment length of caching is 16ms.
Further, the time interval of electrocardiosignal sampling is constant duration 2ms.
Further, the electrocardiosignal change rate threshold value L of the unit interval takes in the range of [10 μ V/ms, 30 μ V/ms]
Value, the electrocardiosignal change rate threshold value H of the unit interval value in the range of [60 μ V/ms, 100 μ V/ms].
The embodiment of the present invention also provides the system for inhibiting notch filter ring in a kind of electrocardiosignal acquisition process, described
System includes caching component, the detection of jumbo jump signal and suppression component, small size jump signal detection part, trap component and eases up
It deposits and determines component constantly, wherein,
Cache component：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1
The electrocardiosignal d [0] of trap output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the period
Length value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches the moment
ti-1Closest to current sample time ti；
Jumbo jump signal detection and suppression component：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] is carried out
Jumbo jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], into small size jump signal detect walk
Suddenly；
Small size jump signal detection part：It calculates in moment t1、...、ti-1、tiEach electrocardiosignal slope absolute value | k
[n]|：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval the heart
Change in electric rate threshold value L, and | k [i] | more than the electrocardiosignal change rate threshold value H of unit interval, then when judging present sample
Carve tiElectrocardiosignal d [i] there are small size mutation；
Trap component：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first is used to it
Wave filter carries out trap output；
The caching moment determines component：If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then really
Determine current sample time tiBelong to the caching moment.
The method for inhibiting notch filter ring in electrocardiosignal acquisition process that the embodiment of the present invention proposes has meter
Calculate the advantages of efficient, realization is easy, strong applicability.
Description of the drawings
Fig. 1 is the exemplary waveforms of electrocardiosignal；
Fig. 2 is the waveform of the original electro-cardiologic signals before trap；
Fig. 3 is not use the inhibition notch filter alerting method of the embodiment of the present invention and the electrocardiosignal of ring occurs
Waveform；
Fig. 4 is the method stream for inhibiting notch filter ring in electrocardiosignal acquisition process of one embodiment of the invention
Cheng Tu；
Fig. 5 is the waveform for inhibiting the electrocardiosignal that notch filter alerting method obtains using the embodiment of the present invention；
Fig. 6 is the method stream for inhibiting notch filter ring in electrocardiosignal acquisition process of yet another embodiment of the invention
Cheng Tu.
Specific embodiment
Understand to make the object, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and reference
Attached drawing, the present invention is described in more detail.But those skilled in the art know, the invention is not limited in attached drawings and following reality
Apply example.
The embodiment of the present invention proposes a kind of method for inhibiting notch filter ring in electrocardiosignal acquisition process,
The electrocardiosignal of acquisition is handled before trap, the original electro-cardiologic signals before trap are for example as shown in Fig. 2, when transverse axis represents
Between, unit is millisecond ms, and the longitudinal axis is electrocardiosignal intensity, and unit is microvolt μ V.It will be appreciated by those skilled in the art that Fig. 2 is only
It is schematic diagram, human ecg signal obtains in real time, trap is carried out in real time using the method for the embodiment of the present invention, so as to inhibit
Notch filter ring.
As shown in figure 4, inhibit the side of notch filter ring in the electrocardiosignal acquisition process that the embodiment of the present invention proposes
Method in real time samples electrocardiosignal, obtains electrocardiosignal, includes the following steps：
Caching step：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1
The electrocardiosignal d [0] of trap output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the period
Length value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches the moment
ti-1Closest to current sample time ti。
Jumbo jump signal detection and inhibition step：Obtain current sample time tiOriginal electro-cardiologic signals d0It is [i], previous
Cache moment ti-1Electrocardiosignal d [i-1], the sampling time intervals of two electrocardiosignals is ti-ti-1, when calculating present sample
Carving electrocardiosignal slope k [i] is：
K [i]=(d0[i]-d[i-1])/(ti-ti-1),
If the absolute value of current sample time electrocardiosignal slope k [i] | k [i] | more than the heart of the unit interval of setting
Change in electric threshold value A, it is determined that the electrocardiosignal of current sample time be jumbo jump signal, current sample time tiThe heart
Electric signal d [i] is：
D [i]=d [i-1]+(k [i]/| k [i] |) * A,
So as to force that jumbo jump signal is inhibited, into small size jump signal detecting step.If present sample
Moment tiElectrocardiosignal slope k [i] absolute value be less than setting unit interval electrocardiosignal change threshold A, it is determined that
Current sample time tiElectrocardiosignal jumbo mutation, d [i]=d do not occur0[i], into small size jump signal detecting step.
The electrocardiosignal change threshold A of the unit interval of the setting preferably sets the electrocardiosignal change threshold A per ms
For 1000 μ V, can also in the range of [800 μ V, 1200 μ V] value.The electrocardiosignal variation threshold of the unit interval of the setting
The unit of value A is related to electrocardiosignal and chronomere.
It will be understood by those skilled in the art that the detection of jumbo jump signal is relatively easy, to the correction of jumbo jump signal/
Inhibit there is also a variety of methods, and be not limited to the aforesaid way of proposition of the embodiment of the present invention, such as can also be by comparing front and rear
The change rate of two sampled signals carries out jumbo jump signal detection, and then decides whether to be inhibited.
Small size jump signal detecting step：It calculates in moment t1、t2、...、ti-2、ti-1、tiEach electrocardiosignal slope
Absolute value | k [n] | (that is, | k [1] |, | k [2] | ..., | k [i-1] |, | k [i] |) be：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval the heart
Change in electric rate threshold value L=20uV/ms, and | k [i] | more than the electrocardiosignal change rate threshold value H=80uV/ of unit interval
Ms then judges current sample time tiElectrocardiosignal d [i] there are small size mutation.
Electrocardiosignal change rate threshold value L=20uV/ms, the electrocardiosignal change rate threshold of unit interval of above-mentioned unit interval
Value H=80uV/ms is the preferable example of the embodiment of the present invention, and the electrocardiosignal change rate threshold value L of the unit interval can be
Value in the range of [10 μ V/ms, 30 μ V/ms], the electrocardiosignal change rate threshold value H of the unit interval can be in [60 μ V/
Ms, 100 μ V/ms] in the range of value.
Trap step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first is used to it
After wave filter carries out trap output.If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then it is adopted
Trap output is carried out with second filter.
The caching moment determines step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, it is determined that
Current sample time tiIt is not belonging to the caching moment；If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then
Determine current sample time tiBelong to the caching moment；
Return to step：Return cache step carries out the processing of the electrocardiosignal of subsequent sampling.
That is, if current sample time tiElectrocardiosignal d [i] there are small size mutation, determine sampling instant tiIt is not belonging to delay
It deposits the moment, then to next sampling instant ti+1Electrocardiosignal be treated in journey, in its caching step, not to sampling when
Carve tiThe electrocardiosignal of trap output is cached.With the lasting progress of the method for the embodiment of the present invention, the electrocardio of caching is believed
Number it is that there is no the electrocardiosignals being slightly mutated.
Trap is carried out and using second to there is the electrocardiosignal being slightly mutated using first filter in the present embodiment
Wave filter is that in general, wave filter has caching in itself to there is the reason for electrocardiosignal being slightly mutated carries out trap
Or feedback, the previous output of wave filter, which can currently export wave filter, to be had an impact, and is handled respectively using two separate filters
Slightly mutation electrocardiosignal and normal electrocardiosignal, can be to avoid small size mutation electrocardiosignal to the trap shadow of normal electrocardiosignal
It rings, obtains better output effect.Certainly, if the first filter used without caching or can ignore above-mentioned shadow in itself
It rings, then first filter can also be used for normal electrocardiosignal and carry out trap.The wave filter used in the present embodiment can be with
For FIR filter or iir filter.
It should be noted that when starting to obtain electrocardiosignal, since the electrocardiosignal of caching is less, inhibit notch filter
The effect of device ring is limited, but the lasting progress of the method with the embodiment of the present invention, updates the data caching, can obtain preferably
Output image.Using the method for inhibiting notch filter ring in electrocardiosignal acquisition process of the embodiment of the present invention, energy
It is enough that the electrocardiosignal of acquisition is handled, it is effective to inhibit notch filter ring, obtain electrocardiosignal figure as shown in Figure 5
Picture, transverse axis represent the time, and unit is millisecond ms, and the longitudinal axis is electrocardiosignal intensity, and unit is microvolt μ V.
In yet another embodiment of the present invention, inhibition notch filter ring in a kind of electrocardiosignal acquisition process is provided
Method, electrocardiosignal is sampled in real time, obtain electrocardiosignal, include the following steps：
Caching step：The constant duration that the time interval of electrocardiosignal sampling is 2ms, caching current sample time tiIt
Multiple caching moment t in a preceding period 16ms0、t1、t2、...、t7Electrocardiosignal d [0], the d of trap output
[1]、...、d[7]。
Jumbo jump signal detection and inhibition step：To current sample time t8Original electro-cardiologic signals d0[8] carry out jumbo
Jump signal detection inhibits, and obtains current sample time t8Electrocardiosignal d [8], into small size jump signal detecting step.
Small size jump signal detecting step：It calculates in moment t1、t2、...、t7、t8Each electrocardiosignal slope absolute value
| k [n] | (that is, | k [1] |, | k [2] | ..., | k [7] |, | k [8] |) be：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1,2 ..., 8, if | k [1] |, | k [2] | ..., | k [7] | average value be less than the unit interval the heart
Change in electric rate threshold value L=20uV/ms, and | k [8] | more than the electrocardiosignal change rate threshold value H=80uV/ of unit interval
Ms then judges current sample time t8Electrocardiosignal d [i] there are small size mutation.
Trap step：Due to current sample time t8Electrocardiosignal d [8] there are small size mutation, therefore first is used to it
Wave filter carries out trap output.
The caching moment determines step：Current sample time t8Electrocardiosignal d [8] there are small size mutation, thus determine work as
Preceding sampling instant t8It is not belonging to the caching moment.
Return to step：Return cache step carries out subsequent sampling moment t9Electrocardiosignal d [9] processing.Carry out under
One moment t9Electrocardiosignal d [9] processing procedure in, in its caching step, due to sampling instant t8The caching moment is not belonging to,
The electrocardiosignal cached is seven caching moment t1、t2、...、t6、t7, trap output electrocardiosignal d [1], d [2] ...,
d[6]、d[7]。
As shown in fig. 6, in one more embodiment of the present invention, provide and inhibit trap filter in a kind of electrocardiosignal acquisition process
The method of ripple device ring, in real time samples electrocardiosignal, obtains electrocardiosignal, includes the following steps：
Caching step：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1
The electrocardiosignal d [0] of trap output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the period
Length value in the range of [10ms, 30ms].
Time interval judgment step：Obtain current sample time tiWith previous caching moment ti-1Between time interval ti-
ti-1If time interval ti-ti-1More than the time threshold 500ms of setting, then using the second notch filter to present sample
Moment tiElectrocardiosignal d [i] carry out trap output, and determine current sample time tiBelong to caching moment, return cache step
Suddenly, subsequent time t is carried outi+1ECG's data compression；If time interval ti-ti-1No more than 500ms, then into jumbo mutation
Signal detection and inhibition step.
Jumbo jump signal detection and inhibition step：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] is carried out
Jumbo jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], subsequently enter the inspection of small size jump signal
Survey step.
Small size jump signal detecting step：It calculates in moment t1、t2、...、ti-2、ti-1、tiEach electrocardiosignal slope
Absolute value | k [n] | (that is, | k [1] |, | k [2] | ..., | k [i-1] |, | k [i] |) be：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval the heart
Change in electric rate threshold value L=20uV/ms, and | k [i] | more than the electrocardiosignal change rate threshold value H=80uV/ of unit interval
Ms then judges current sample time tiElectrocardiosignal d [i] there are small size mutation.
Trap step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first is used to it
After wave filter carries out trap output.If current electrocardiosignal d [i] uses second filter there is no slightly mutation to it
Carry out trap output.
The caching moment determines step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, it is determined that
Current sample time tiIt is not belonging to the caching moment；If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then
Determine current sample time tiBelong to the caching moment；
Return to step：Return cache step carries out the processing of the electrocardiosignal of subsequent sampling.
The embodiment of the present invention also provides the system for inhibiting notch filter ring in a kind of electrocardiosignal acquisition process, including
It is true to cache component, the detection of jumbo jump signal and suppression component, small size jump signal detection part, trap component and caching moment
Determine component, wherein,
Cache component：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1
The electrocardiosignal d [0] of trap output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the period
Length value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches the moment
ti-1Closest to current sample time ti；
Jumbo jump signal detection and suppression component：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] is carried out
Jumbo jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], into small size jump signal detect walk
Suddenly；
Small size jump signal detection part：It calculates in moment t1、...、ti-1、tiEach electrocardiosignal slope absolute value | k
[n]|：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval the heart
Change in electric rate threshold value L, and | k [i] | more than the electrocardiosignal change rate threshold value H of unit interval, then when judging present sample
Carve tiElectrocardiosignal d [i] there are small size mutation；
Trap component：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first is used to it
Wave filter carries out trap output；
The caching moment determines component：If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, then really
Determine current sample time tiBelong to the caching moment.
The embodiment of the present invention also proposes a kind of computer readable storage medium, is stored thereon with computer program, the calculating
The step of above method is realized when machine program is executed by processor.
A kind of computer equipment can be run on a memory and on a processor including memory, processor and storage
The step of computer program, the processor realizes the above method when performing described program.
It will be understood by those skilled in the art that in flow charts represent or logic described otherwise above herein and/or
Step for example, being considered the order list for the executable instruction for being used to implement logic function, may be embodied in and appoint
In what computer-readable medium, for instruction execution system, device or equipment (such as computer based system, including processor
System or other can be from instruction execution system, device or equipment instruction fetch and the system executed instruction) use or with reference to this
A little instruction execution systems, device or equipment and use.For the purpose of this specification, " computer-readable medium " can be it is any can be with
Comprising, store, communicate, propagate or transmit program and perform system for instruction execution system, device or equipment or with reference to these instructions
System, device or equipment and the device used.
The more specific example (non-exhaustive list) of computer-readable medium includes following：It is connected up with one or more
Electrical connection section (electronic device), portable computer diskette box (magnetic device), random access memory (RAM), read-only memory
(ROM), erasable edit read-only storage (EPROM or flash memory), fiber device and portable optic disk is read-only deposits
Reservoir (CDROM).In addition, computer-readable medium can even is that the paper that can print described program on it or other are suitable
Medium, because can be for example by carrying out optical scanner to paper or other media, then into edlin, interpretation or if necessary with it
His suitable method is handled electronically to obtain described program, is then stored in computer storage.
It should be appreciated that each several part of the present invention can be realized with hardware, software, firmware or combination thereof.Above-mentioned
In embodiment, software that multiple steps or method can in memory and by suitable instruction execution system be performed with storage
Or firmware is realized.If for example, with hardware come realize in another embodiment, can be under well known in the art
Any one of row technology or combination thereof are realized：With for the logic gates to data-signal realization logic function
Discrete logic, have suitable combinational logic gate circuit application-specific integrated circuit, programmable gate array (PGA), scene
Programmable gate array (FPGA) etc..
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or the spy for combining the embodiment or example description
Point is contained at least one embodiment of the present invention or example.In the present specification, schematic expression of the above terms is not
Centainly refer to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be any
One or more embodiments or example in combine in an appropriate manner.
More than, embodiments of the present invention are illustrated.But the present invention is not limited to the above embodiments.It is all
Within the spirit and principles in the present invention, any modification, equivalent substitution, improvement and etc. done should be included in the guarantor of the present invention
Within the scope of shield.
Claims (10)
1. a kind of electrocardiosignal inhibits notch filter ring method during using, it is characterised in that：Include the following steps：
Caching step：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1Trap
The electrocardiosignal d [0] of output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the length of the period
The value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches moment ti-1Most
Close to current sample time ti；
Jumbo jump signal detection and inhibition step：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] carries out jumbo
Jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], into small size jump signal detecting step；
Small size jump signal detecting step：It calculates in moment t1、...、ti-1、tiEach electrocardiosignal slope absolute value | k [n]
|：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval electrocardio believe
Number change rate threshold value L, and | k [i] | more than the electrocardiosignal change rate threshold value H of unit interval, then judge current sample time ti
Electrocardiosignal d [i] there are small size mutation；
Trap step：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first filter is used to it
Carry out trap output；
The caching moment determines step：If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, it is determined that it is current
Sampling instant tiBelong to the caching moment；
Return to step：Return cache step.
2. the method as described in claim 1, it is characterised in that：In jumbo jump signal detects and inhibits step, obtain and work as
Preceding sampling instant tiOriginal electro-cardiologic signals d0[i], previous caching moment ti-1Electrocardiosignal d [i-1], calculate present sample when
Carving electrocardiosignal slope k [i] is：
K [i]=(d0[i]-d[i-1])/(ti-ti-1),
If the absolute value of current sample time electrocardiosignal slope k [i] | k [i] | more than the electrocardio letter of the unit interval of setting
Number change threshold A, it is determined that the electrocardiosignal of current sample time is jumbo jump signal, current sample time tiElectrocardio letter
Number d [i] is：
D [i]=d [i-1]+(k [i]/| k [i] |) * A；
If current sample time tiElectrocardiosignal slope k [i] absolute value be less than setting unit interval electrocardiosignal become
Change threshold value A, current sample time tiElectrocardiosignal d [i] be：D [i]=d0[i]。
3. method as claimed in claim 2, it is characterised in that：In jumbo jump signal detects and inhibits step, setting
The electrocardiosignal change threshold A of every ms values in the range of [800 μ V, 1200 μ V].
4. the method as described in claim 1, it is characterised in that：If current sample time tiElectrocardiosignal d [i] be not present
Slightly mutation, then it is used in trap step first filter carry out trap output, wherein first filter be without
The wave filter of caching function.
5. the method as described in claim 1, it is characterised in that：If current sample time tiElectrocardiosignal d [i] be not present
Slightly mutation then uses it in trap step second filter to carry out trap output.
6. method as claimed in claim 5, it is characterised in that：After caching step, the method further includes time interval
Judgment step obtains current sample time tiWith previous caching moment ti-1Between time interval ti-ti-1If time interval
ti-ti-1More than the time threshold of setting, then using the second notch filter to current sample time tiOriginal electro-cardiologic signals d0
[i] carries out trap output, and determines current sample time tiBelong to the caching moment, return cache step carries out subsequent time
ECG's data compression；If time interval ti-ti-1No more than the time threshold of the setting, then examined into jumbo jump signal
It surveys and inhibits step.
7. the method as described in claim 1, it is characterised in that：In caching step, the time segment length of caching is 16ms.
8. the method for claim 7, it is characterised in that：The time interval of electrocardiosignal sampling is constant duration 2ms.
9. method as described in one of claim 1-8, it is characterised in that：The electrocardiosignal change rate threshold value of the unit interval
L values in the range of [10 μ V/ms, 30 μ V/ms], the electrocardiosignal change rate threshold value H of the unit interval [60 μ V/ms,
100 μ V/ms] in the range of value.
10. inhibit the system of notch filter ring in a kind of electrocardiosignal acquisition process, it is characterised in that：The system comprises
It is true to cache component, the detection of jumbo jump signal and suppression component, small size jump signal detection part, trap component and caching moment
Determine component, wherein,
Cache component：Cache current sample time tiMultiple caching moment t in a period before0、t1、...、ti-1Trap
The electrocardiosignal d [0] of output, d [1] ..., d [i-1], wherein i be natural number more than or equal to 2, the length of the period
The value in the range of [10ms, 30ms], the caching moment t0、t1、...、ti-1It in chronological sequence arranges, caches moment ti-1Most
Close to current sample time ti；
Jumbo jump signal detection and suppression component：To current sample time tiThe original electro-cardiologic signals d of sampling0[i] carries out jumbo
Jump signal detection inhibits, and obtains current sample time tiElectrocardiosignal d [i], into small size jump signal detecting step；
Small size jump signal detection part：It calculates in moment t1、...、ti-1、tiEach electrocardiosignal slope absolute value | k [n]
|：
| k [n] |=| (d [n]-d [n-1])/(tn-tn-1) |,
Wherein n=1 ..., i, if | k [1] |, | k [2] | ..., | k [i-1] | average value be less than the unit interval electrocardio believe
Number change rate threshold value L, and | k [i] | more than the electrocardiosignal change rate threshold value H of unit interval, then judge current sample time ti
Electrocardiosignal d [i] there are small size mutation；
Trap component：If current sample time tiElectrocardiosignal d [i] there are small size mutation, then first filter is used to it
Carry out trap output；
The caching moment determines component：If current sample time tiElectrocardiosignal d [i] there is no slightly mutation, it is determined that it is current
Sampling instant tiBelong to the caching moment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711286389.0A CN108113664B (en) | 2017-12-07 | 2017-12-07 | Method and system for inhibiting ringing of notch filter in electrocardiosignal acquisition process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711286389.0A CN108113664B (en) | 2017-12-07 | 2017-12-07 | Method and system for inhibiting ringing of notch filter in electrocardiosignal acquisition process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108113664A true CN108113664A (en) | 2018-06-05 |
| CN108113664B CN108113664B (en) | 2020-09-18 |
Family
ID=62228971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711286389.0A Active CN108113664B (en) | 2017-12-07 | 2017-12-07 | Method and system for inhibiting ringing of notch filter in electrocardiosignal acquisition process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108113664B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109077714A (en) * | 2018-07-05 | 2018-12-25 | 广州视源电子科技股份有限公司 | Signal recognition method, device, equipment and storage medium |
| CN109567788A (en) * | 2018-11-29 | 2019-04-05 | 武汉中旗生物医疗电子有限公司 | A kind of Filtering of ECG Signal method removing ring |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080004904A1 (en) * | 2006-06-30 | 2008-01-03 | Tran Bao Q | Systems and methods for providing interoperability among healthcare devices |
| CN101953684A (en) * | 2009-07-17 | 2011-01-26 | 深圳市理邦精密仪器股份有限公司 | Power frequency interference filtering method for eliminating ringing |
| CN103813746A (en) * | 2011-04-22 | 2014-05-21 | 德尔格医疗系统有限公司 | Adaptive notch filter |
| CN104502925A (en) * | 2015-01-22 | 2015-04-08 | 西安航天华迅科技有限公司 | Narrowband interference resisting system and method based on adaptive signal processing |
| CN106725433A (en) * | 2016-12-30 | 2017-05-31 | 包磊 | A kind of cardioelectric monitor method and system |
-
2017
- 2017-12-07 CN CN201711286389.0A patent/CN108113664B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080004904A1 (en) * | 2006-06-30 | 2008-01-03 | Tran Bao Q | Systems and methods for providing interoperability among healthcare devices |
| CN101953684A (en) * | 2009-07-17 | 2011-01-26 | 深圳市理邦精密仪器股份有限公司 | Power frequency interference filtering method for eliminating ringing |
| CN103813746A (en) * | 2011-04-22 | 2014-05-21 | 德尔格医疗系统有限公司 | Adaptive notch filter |
| CN104502925A (en) * | 2015-01-22 | 2015-04-08 | 西安航天华迅科技有限公司 | Narrowband interference resisting system and method based on adaptive signal processing |
| CN106725433A (en) * | 2016-12-30 | 2017-05-31 | 包磊 | A kind of cardioelectric monitor method and system |
Non-Patent Citations (1)
| Title |
|---|
| 王吉鸣等: "一种滤除高采样心电工频干扰的改进算法", 《中国生物医学工程学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109077714A (en) * | 2018-07-05 | 2018-12-25 | 广州视源电子科技股份有限公司 | Signal recognition method, device, equipment and storage medium |
| CN109077714B (en) * | 2018-07-05 | 2021-03-23 | 广州视源电子科技股份有限公司 | Signal identification method, device, equipment and storage medium |
| CN109567788A (en) * | 2018-11-29 | 2019-04-05 | 武汉中旗生物医疗电子有限公司 | A kind of Filtering of ECG Signal method removing ring |
| CN109567788B (en) * | 2018-11-29 | 2021-08-20 | 武汉中旗生物医疗电子有限公司 | Electrocardiosignal filtering method for removing ringing |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108113664B (en) | 2020-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108158573B (en) | Electrocardiosignal noise reduction method based on adaptive threshold wavelet transformation | |
| CA1210075A (en) | Method for predicting ventricular tachycardia with a gain controlled high pass filter | |
| CA1210457A (en) | System and method for predicting ventricular tachycardia by adaptive high pass filter | |
| CN109907752B (en) | Electrocardiogram diagnosis and monitoring system for removing motion artifact interference and electrocardio characteristic detection | |
| US8112147B2 (en) | Method and apparatus for generating determination indexes for identifying ECG interfering signals | |
| CN1829553B (en) | Method and devices for performing cardiac waveform appraisal | |
| US20070197926A1 (en) | Method for analysing an intracardiac electrocardiogram and an electrophysiological system as well as a computer program product | |
| US5908393A (en) | Reducing noise in a biological signal | |
| CN108113665B (en) | Automatic noise reduction method for electrocardiosignal | |
| CN107361764B (en) | Method for rapidly extracting electrocardiosignal characteristic waveform R wave | |
| CN108720832B (en) | Electrocardiosignal processing method and device | |
| CN110367936B (en) | Electrocardiosignal detection method and device | |
| CN108113664A (en) | Inhibit the method and system of notch filter ring in electrocardiosignal acquisition process | |
| US20170281021A1 (en) | Heartbeat detection method and heartbeat detection device | |
| Singh et al. | ECG signal denoising based on empirical mode decomposition and moving average filter | |
| Li et al. | Application of an EMG interference filtering method to dynamic ECGs based on an adaptive wavelet-Wiener filter and adaptive moving average filter | |
| Pourazad et al. | Heart sound cancellation from lung sound recordings using adaptive threshold and 2D interpolation in time-frequency domain | |
| CN107280659A (en) | The processing method and system of a kind of electrocardiosignal | |
| CN106529425B (en) | A kind of the R wave extracting method and system of electrocardiosignal | |
| KR20120122856A (en) | R-peak detection method in ECG Signal and apparatus thereof | |
| Chen et al. | An ECG R-wave detection algorithm based on adaptive threshold | |
| CN113786202B (en) | Electrocardiogram characteristic starting point and ending point detection method, system, device and medium | |
| Naaz et al. | Feature extraction and analysis of ECG signal for cardiac abnormalities-A review | |
| CN115293219A (en) | Wavelet and kurtosis fused pulse signal denoising method | |
| Kim et al. | A biomedical signal segmentation algorithm for event detection based on slope tracing |
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 |