CN101869476B - Method for checking heartbeat - Google Patents
Method for checking heartbeat Download PDFInfo
- Publication number
- CN101869476B CN101869476B CN 200910134758 CN200910134758A CN101869476B CN 101869476 B CN101869476 B CN 101869476B CN 200910134758 CN200910134758 CN 200910134758 CN 200910134758 A CN200910134758 A CN 200910134758A CN 101869476 B CN101869476 B CN 101869476B
- Authority
- CN
- China
- Prior art keywords
- trough
- heart beating
- peak
- primary peak
- ecg signal
- 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
Images
Abstract
The invention relates to a method for checking a heartbeat. The method comprises the following steps: providing an electrocardiogram (ECG) signal; increasing the ECG signal to generate a processing signal; checking the first crest of the processing signal; defining the crest after the first crest to be a second crest, wherein the second crest is more than a first threshold value; and judging whether a heartbeat exists near the first crest according to the second crest.
Description
Technical field
The present invention relates to a kind of detection method, particularly a kind of detection heart beating method.
Background technology
Electrocardiogram (ECG) is an important indicator of assessment cardiovascular disease.Just can learn the time and interval of heart beating by the QRS ripple of ECG signal, and then learn whether heart beating is unusual.Fig. 1 is the sketch map of ECG signal.Heart beating each time can make ECG signal produce P ripple, QRS complex wave and T ripple in regular turn.The P ripple occurs in the section 110, and the QRS complex wave occurs in the section 120, and the T ripple occurs in the section 130.As shown in the figure, the QRS complex wave in the section 120 is the most obvious.
Yet in some cases, the QRS complex wave can't be obvious as Fig. 1.Fig. 2 a and Fig. 2 b show improper ECG signal.In Fig. 2 a, QRS complex wave and T ripple are reverse.In Fig. 2 b, the QRS complex wave is also different with normal condition.Therefore, when ECG signal is shown in Fig. 2 a or 2b, just can't correctly learn the time that heart beating takes place.
Summary of the invention
The invention provides a kind of method that detects heart beating, comprise the following steps, an ECG signal is provided; Strengthen this ECG signal, to produce a processing signals; Detect a primary peak of this processing signals; Crest after this primary peak is defined as a secondary peak, and this secondary peak is greater than a first threshold; And according to this secondary peak, learn near this primary peak, whether to have a heart beating.
For the features and advantages of the present invention can be become apparent, cited below particularlyly go out preferred embodiment, and conjunction with figs., be described in detail below:
Description of drawings
Fig. 1 is the sketch map of ECG signal.
Fig. 2 a and Fig. 2 b show improper ECG signal.
Fig. 3 a is that one of detection heart beating method of the present invention may embodiment.
Fig. 3 b is the relation between secondary peak and the heart beating.
The sketch map of Fig. 4 a and Fig. 4 b display process signal.
Fig. 5 a is another possibility embodiment of detection heart beating method of the present invention.
Fig. 5 b is the relation between secondary peak and second trough and the heart beating.
Fig. 6 is that one of enhancing step of the present invention may embodiment
[main element symbol description]
110,120,130: section;
S310～S350, S510～S570, S610～S630: step;
41～44: curve;
421～429,441,443,445,447: crest;
442,444,448: trough.
The specific embodiment
Fig. 3 a is that one of detection heart beating method of the present invention may embodiment.At first, provide an ECG signal (step S310).ECG signal may be shown in Fig. 2 a or 2b.
Then, strengthen (enhance) ECG signal, to produce a processing signals (step S320).The sketch map of Fig. 4 a and Fig. 4 b display process signal.Curve 41 and 43 is respectively the ECG signal shown in Fig. 2 a and Fig. 2 b, is the described ECG signal of step S310.After curve 41 and 43 is enhanced, just can obtain curve 42 and 44.Curve 42 and 44 is the described processing signals of step S320.
Detect a primary peak (step S330) of processing signals.In the present embodiment, primary peak needs greater than a first threshold.Therefore, with the crest of processing signals earlier with first threshold make comparisons (step S331).If the crest of processing signals not greater than first threshold, then continues to look for next crest (step S332) again.If the crest of processing signals greater than first threshold, then is defined as this crest primary peak (step S333).In other embodiments, primary peak can be less than first threshold.For example, so long as the crest of processing signals all can be defined by primary peak.
Then, the crest after the primary peak is defined as secondary peak (step S340).In the present embodiment, first and second crest is all greater than a first threshold, but is not in order to limit the present invention.In other embodiments, first and/or secondary peak less than first threshold.
According to secondary peak, learn after primary peak, whether to have a heart beating (step S350).For example, can learn after primary peak, whether to have heart beating according to quantity, amplitude or the position of secondary peak.
In a possibility embodiment, after defining primary peak, pick up counting, in order to judge the quantity of the secondary peak in cycle time.For example, in the 300ms after learning primary peak (person of outstanding talent's second), by calculating the quantity of secondary peak, just can learn near primary peak, whether to have the heart beating generation.
In the present embodiment, behind the execution of step S350, can get back to step S330, in order to redefine next primary peak.By repeating step S330～S350, just can learn the interval between heart beating.Relation between secondary peak and the heart beating below will be described.
Please refer to Fig. 3 b, Fig. 3 b is the relation between secondary peak and the heart beating.The shown curve of Fig. 3 b is identical with the curve 42 of Fig. 4 a, is the ECG signal (being processing signals) after strengthening.Wherein symbol N is the quantity of secondary peak.
Because the crest 421 of processing signals is greater than first threshold, so crest 421 is defined by primary peak Smax1.As shown in the figure, in a period of time behind crest 421, crest 422 and 423 is not greater than first threshold.Therefore, the secondary peak quantity N behind the primary peak Smax1 equals 0, and this has heart beating near being illustrated in primary peak Smax1.
Suppose that the crest of processing signals needs greater than first threshold, just can be defined by primary peak.Because the crest 422 and 423 of processing signals is not all greater than first threshold, so can't be defined by primary peak.
Because the crest 424 of processing signals is greater than first threshold, so crest 424 can be defined by primary peak Smax2.In a period of time behind crest 424, crest 425 is greater than first threshold, and crest 426 is not greater than first threshold, so have only crest 425 to be defined by secondary peak.Therefore, the secondary peak quantity N behind the primary peak Smax2 equals 1, and this near being illustrated in primary peak Smax2 heart beating takes place.
Similarly, the crest 427 of processing signals is greater than first threshold, so crest 427 is defined by primary peak Smax3.In a period of time behind crest 427, crest 428 is greater than first threshold, and crest 429 is not greater than first threshold, so have only crest 428 to be defined by secondary peak.Therefore, the secondary peak quantity N behind the primary peak Smax3 equals 1, and this near being illustrated in primary peak Smax3 heart beating takes place.
In the present embodiment, primary peak Smax1～Smax3 more greater than the 3rd threshold value (as 0.6), but is not in order to limit the present invention not only greater than first threshold (as 0.2).In other embodiments, primary peak Smax1～Smax3 can be less than first threshold.
Fig. 5 a is another possibility embodiment of detection heart beating method of the present invention.Fig. 5 a similar diagram 3a, difference has been more than Fig. 5 a step S540, S560.For simplicity, identical with Fig. 3 a step repeats no more.In Fig. 5 a, after defining primary peak (step S530), define first trough (step S540).In the present embodiment, first trough needs less than second threshold value.For reaching this purpose, the trough of processing signals and second threshold value can be made comparisons (step S541).If the trough of processing signals is not less than second threshold value, then look for the next trough (step S542) of processing signals again.When the trough of processing signals during less than second threshold value, then this trough is defined as first trough (step S543).
After having defined secondary peak, define second trough (step S560).In the present embodiment, the trough less than second threshold value just can be defined by second trough.At last, according to secondary peak and second trough, just can learn between primary peak and first trough whether have heart beating.In some may be implemented, can judge between this primary peak and first trough whether have heart beating by quantity, amplitude or the position of secondary peak and second trough.
Relation between secondary peak and second trough and the heart beating below will be described.Please refer to Fig. 5 b, Fig. 5 b shows the relation between secondary peak and second trough and the heart beating.The shown curve of Fig. 5 b is identical with the curve 44 of Fig. 4 b, is the ECG signal (being processing signals) after strengthening.Wherein symbol N is the quantity of secondary peak, and symbol n is the quantity of second trough.
Because the crest 441 of processing signals is greater than first threshold, so crest 441 is defined by primary peak Smax1.In addition, the trough 442 of processing signals is less than second threshold value, so trough 442 is defined by the first trough Smin1.In the present embodiment, behind the definition primary peak, just define first trough.
As shown in the figure, in a period of time behind primary peak Smax1, crest 445 and 447 is greater than first threshold, and crest 443 is not greater than first threshold.Therefore, the secondary peak quantity N behind the primary peak Smax1 equals 2, may not be the QRS ripple between this expression primary peak Smax1 and the first trough Smin1.In other embodiments, the secondary peak quantity N behind the primary peak Smax1 may be greater than 2.
In addition, in a period of time behind the first trough Smin1, trough 446 is less than second threshold value, and trough 444 and 448 is not less than second threshold value.Therefore, the second trough quantity n behind the first trough Smin1 equals 1, may not be the QRS ripple between this expression primary peak Smax1 and the first trough Smin1.In other situation, the second trough quantity n equals 1 or may be greater than more than 2, near this expression primary peak noise takes place, but not heart beating.Therefore, by the second trough quantity n, but between auxiliary judgment primary peak and first trough whether noise takes place.
Then, continue definition primary peak Smax1 primary peak afterwards.In the present embodiment, the crest of processing signals needs greater than first threshold, so crest 443 is not defined by primary peak.Because crest 445 is greater than greater than first threshold, so crest 445 definition primary peak Smax2.After primary peak Smax2, because trough 446 is less than second threshold value, so trough 446 is defined by the first trough Smin2.
After primary peak Smax2, because crest 447 is greater than first threshold, so can be defined by secondary peak.Therefore, the quantity N of the secondary peak after the primary peak Smax2 equals 1, and this is illustrated between primary peak Smax2 and the first trough Smin2 has heart beating.In addition, after the first trough Smin2, processing signals does not have the trough less than second threshold value, so the quantity n of second trough equals 0.Because the quantity N of the secondary peak after the primary peak Smax2 equals 1, therefore when the quantity n of second trough after the first trough Smin2 equals 0, determine between primary peak Smax2 and the first trough Smin2, to have heart beating more.
In addition, if secondary peak and second trough is equal in number at 0 o'clock, is illustrated between primary peak and first trough and does not have heart beating.But, if primary peak greater than the 3rd threshold value, even if secondary peak and second trough is equal in number in 0, has heart beating near still can being illustrated in primary peak.
Fig. 6 is that one of enhancing step of the present invention may embodiment.Utilize a filter process ECG signal (step S610).In a possibility embodiment, can utilize a low pass filter (low-pass filter) or a band filter (band-pass filter) to handle ECG signal.In addition, in the present embodiment, the frequency by wave filter is 0～5Hz, and the cut-off frequency of wave filter is 12～18Hz.
Utilize the ECG signal (step S620) after a shade (mask) is handled filter process, in order to the QRS complex wave of enhancing (enhance) ECG signal, and inhibition P involves the T ripple.This shade can be 3 rank or 3 rank above (for symmetry, all can get odd-order).If be 3 rank, can be [12-1] or [101].If be 5 rank, can be [1-14-1-1] or [2-1012].
Step S620 mainly is that emphasized point changes or slope variation.For the shade on p rank, p is odd number, can be expressed as p=2k+1.In a possibility embodiment, utilisation point strengthens (point enhance) form, the ECG signal after the processing filter process.In another possibility embodiment, can utilize gradient to strengthen (gradient enhance) form, the ECG signal after the processing filter process.
The shade that strengthens the some variation is shown below:
Wherein, M is shade, k relevant with shade exponent number p (p=2*k+1).
The shade that increases slope variation is shown below:
M(n)＝-k+n，n＝0～2k
ECG signal (step S630) after standardization (normalize) shade processing.In a possibility embodiment, after the process standardization, the maximum of signal will be 1, in order to highlight the QRS complex wave.For example, if original ECG signal shown in the curve 41 of Fig. 4 a the time then after through step S610～630, just can produce processing signals shown in curve 42, wherein in step S620, utilizes point to increase a shade.In another may embodiment, if original ECG signal shown in the curve 43 of Fig. 4 b the time then after through step S610～630, just can produce processing signals shown in curve 44, wherein in step S620, utilizes gradient to increase shade.
Because the QRS complex wave of ECG signal may be not obvious, so can be by strengthening ECG signal, the ECG signal that is enhanced (being above-mentioned processing signals).Then, detect the primary peak of processing signals, with the position of learning that the QRS complex wave is possible.By quantity, amplitude and the position of the secondary peak behind the primary peak, just can judge processing signals primary peak near, whether the QRS complex wave is arranged.
In other embodiments, also can detect first trough of processing signals.By quantity, amplitude and the position of second trough behind first trough, just can judge processing signals primary peak near, whether noise takes place.In addition, by quantity, amplitude and the position of second trough behind first trough, and also can judge between the primary peak of processing signals and first trough whether have the QRS complex wave.
Though the present invention with preferred embodiment openly as above; right its is not in order to limiting the present invention, those skilled in the art, without departing from the spirit and scope of the present invention; when can doing a little change and retouching, so protection scope of the present invention is as the criterion when looking the appended claims person of defining.
Claims (17)
1. one kind is detected the heart beating method, comprising:
One ECG signal is provided;
Strengthen this ECG signal, to produce a processing signals;
Detect a primary peak of this processing signals;
Crest after this primary peak is defined as secondary peak, and this secondary peak is greater than a first threshold; And
According to this secondary peak, learn near this primary peak, whether to have a heart beating;
This enhancing step comprises:
Utilize this ECG signal of a filter process;
Utilize this ECG signal after this filter process of a shade processing; And
This ECG signal after this shade processing of standardization,
Wherein utilize this ECG signal after this filter process of a shade processing, strengthening the QRS complex wave of this ECG signal, and suppress P and involve the T ripple.
2. detection heart beating method as claimed in claim 1 wherein according to quantity, amplitude or the position of this secondary peak, learns whether have this heart beating near this primary peak.
3. detection heart beating method as claimed in claim 1, wherein this primary peak is greater than this first threshold.
4. detection heart beating method as claimed in claim 1 wherein, according to the quantity of the secondary peak in cycle time behind this primary peak, learns whether have this heart beating near this primary peak.
5. detection heart beating method as claimed in claim 1 wherein when the quantity of this secondary peak equals 0 or 1, has this heart beating near being illustrated in this primary peak.
6. detection heart beating method as claimed in claim 1, also comprise: detect one first trough of this processing signals, this first trough occurs in after this primary peak.
7. detection heart beating method as claimed in claim 6, when the quantity of this secondary peak more than or equal to 2 the time, be illustrated between this primary peak and this first trough and do not have this heart beating.
8. detection heart beating method as claimed in claim 7 also comprises:
Trough after this first trough is defined as second trough, and this first and second trough is less than one second threshold value; And
According to this second trough, learn between this primary peak and first trough whether have this heart beating.
9. detection heart beating method as claimed in claim 8 wherein according to quantity, amplitude or the position of this second trough, learns between this primary peak and first trough whether have this heart beating.
10. detection heart beating method as claimed in claim 8 wherein, according to the quantity of second trough in cycle time behind this first trough, learns between this primary peak and first trough whether have this heart beating.
11. detection heart beating method as claimed in claim 8 wherein when the quantity of this second trough is not equal to 0, is illustrated between this primary peak and first trough and does not have this heart beating.
12. detection heart beating method as claimed in claim 11 wherein when the quantity of this second trough equals 1, is illustrated between this primary peak and first trough and does not have this heart beating.
13. detection heart beating method as claimed in claim 11, equal in number in 0 time when this secondary peak and second trough wherein is illustrated between this primary peak and first trough and do not have this heart beating.
14. detection heart beating method as claimed in claim 11, wherein work as this primary peak greater than one the 3rd threshold value, and this secondary peak and second trough is equal in number at 0 o'clock, has this heart beating near being illustrated in this primary peak, and the 3rd threshold value is greater than this first threshold.
15. detection heart beating method as claimed in claim 1, wherein this wave filter is low pass filter or band filter.
16. detection heart beating method as claimed in claim 1, wherein the frequency by this wave filter is 0～5Hz, and the cut-off frequency of this wave filter is 12～18Hz.
17. detection heart beating method as claimed in claim 1, wherein this shade is some enhanced forms or a gradient enhanced form.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200910134758 CN101869476B (en) | 2009-04-22 | 2009-04-22 | Method for checking heartbeat |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200910134758 CN101869476B (en) | 2009-04-22 | 2009-04-22 | Method for checking heartbeat |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101869476A CN101869476A (en) | 2010-10-27 |
| CN101869476B true CN101869476B (en) | 2013-07-17 |
Family
ID=42994697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200910134758 Active CN101869476B (en) | 2009-04-22 | 2009-04-22 | Method for checking heartbeat |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101869476B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109199355B (en) * | 2018-09-18 | 2021-09-28 | 深圳和而泰数据资源与云技术有限公司 | Heart rate information detection method and device and detection equipment |
| CN115225450A (en) * | 2022-09-20 | 2022-10-21 | 南京艾泰克物联网科技有限公司 | Multi-data machine room virtualization cluster management system based on edge computing |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5682900A (en) * | 1995-11-29 | 1997-11-04 | Hewlett-Packard Company | Method and apparatus for obtaining heartbeat measurements from a ECG waveform |
| TW333445B (en) * | 1996-01-25 | 1998-06-11 | Ind Tech Res Inst | Long term monitoring method for heart activities |
-
2009
- 2009-04-22 CN CN 200910134758 patent/CN101869476B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5682900A (en) * | 1995-11-29 | 1997-11-04 | Hewlett-Packard Company | Method and apparatus for obtaining heartbeat measurements from a ECG waveform |
| TW333445B (en) * | 1996-01-25 | 1998-06-11 | Ind Tech Res Inst | Long term monitoring method for heart activities |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101869476A (en) | 2010-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103750820B (en) | A kind of sleep quality monitoring method and device | |
| CN101946475B (en) | Radio sensor for detecting wireless microphone signals and method thereof | |
| JP5653791B2 (en) | FSK demodulation circuit | |
| CN109381181A (en) | The end-point detecting method of electrocardiosignal signature waveform | |
| CN201150533Y (en) | Blood pressure measuring device | |
| CN108261195A (en) | Real-time heart beat detection method and real-time heart beat detecting system | |
| CN100447855C (en) | Methods and devices for obtaining sampling clocks | |
| CN105468951A (en) | Method and device for identity recognition through electrocardiographic feature and wearable device | |
| CN112754441B (en) | Millimeter wave-based non-contact heartbeat detection method | |
| WO2019036908A1 (en) | Heart rate measuring method and apparatus, and electronic terminal | |
| CN101869476B (en) | Method for checking heartbeat | |
| Li et al. | Dynamical complexity detection in short-term physiological series using base-scale entropy | |
| CN110389678A (en) | Parallel detection touch device and its operation method | |
| CN103584854A (en) | Extraction method of electrocardiosignal R waves | |
| CN110974189A (en) | Method, device, equipment and system for detecting signal quality of pulse wave | |
| CN102438334B (en) | Method for carrying out modulation recognition through distributed cooperation of multiple sensor nodes | |
| JP2007206037A (en) | Signal measuring/analyzing apparatus | |
| CN104463057B (en) | A kind of magnetic card coding/decoding method and device | |
| JP4476183B2 (en) | OOB signal detection circuit | |
| CN105094413B (en) | The noise detecting method of touch device, touch control controller and touch panel | |
| Schmidt et al. | No evidence of nonlinear or chaotic behavior of cardiovascular murmurs | |
| JP5320837B2 (en) | Pulse measuring device, pulse measuring program, and pulse measuring method | |
| AU2012244118B2 (en) | Method and System for Identifying Events of Digital Signal | |
| WO2023130869A1 (en) | Respiratory signal acquisition method and apparatus, and computer device | |
| CN110151156B (en) | Fetal heart rate extraction method and system based on in-window automatic weighted average algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |