CN105997036A - A heart rate calculating method based on the heart sound signal time domain features - Google Patents
A heart rate calculating method based on the heart sound signal time domain features Download PDFInfo
- Publication number
- CN105997036A CN105997036A CN201610300369.3A CN201610300369A CN105997036A CN 105997036 A CN105997036 A CN 105997036A CN 201610300369 A CN201610300369 A CN 201610300369A CN 105997036 A CN105997036 A CN 105997036A
- Authority
- CN
- China
- Prior art keywords
- heart rate
- time interval
- cardiechema signals
- sequence
- heart
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Cardiology (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention provides a heart rate calculating method based on the heart sound signal time domain features. The method comprises the steps of firstly obtaining heart sound signals; secondly, extracting local maximum points; thirdly, extracting points meeting the heart rate features; fourthly, calculating the heart rate. The method saves the processes of envelope extraction and combination of heart sound splitting, determines heart sound subsections based on the heart sound time domain features, reduces error by using the median method, and further reduces the time cost for heart rate calculation.
Description
Technical field
The present invention relates to electronic information medicine technology field, be specifically related to a kind of based on cardiechema signals temporal signatures
Rate calculation method.
Background technology
Heart sound refers to be clashed into what ventricle wall, main artery wall etc. caused by myocardial contraction, cardiac valve closedown and blood
The produced sound of vibration.It can be listened at thoracic wall certain position stethoscope, it is also possible to the instruments such as transducer
The mechanical vibration of record heart sound.The heart sound of one heart beat cycle can be divided into first heart sound (S1) second heart sound (S2)
(the most all can hear), third heart sound (S3 the most only can hear child and teenager), the 4th heart
Sound (S4 normal condition is rarely heard).
For comprising the cardiechema signals in multiple successive heartbeat cycle, then can be by positioning in each heart beat cycle
S1, S2 calculate heart rate value.The time interval such as utilizing adjacent two S1 (or S2) can calculate one
Individual instantaneous heart rate value, and the average of multiple instantaneous heart rate is exactly average heart rate value interior during this period of time.Therefore,
In the cardiechema signals of location, the position of S1, S2 becomes the key calculating heart rate based on cardiechema signals.
Existing rate calculation technology often rely on signal envelope extract, arrange Shannon threshold value, merge heart sound divide
Split and determine heart sound segmentation.
Summary of the invention
Instant invention overcomes the deficiencies in the prior art, it is provided that a kind of cardiotachometer based on cardiechema signals temporal signatures
Calculation method.
For achieving the above object, the present invention is by the following technical solutions:
A kind of rate calculation method based on cardiechema signals temporal signatures, it is characterised in that it includes following step
Rapid:
Step 1, obtains cardiechema signals；
Step 2, is divided into N section by cardiechema signals, extracts each section of interior maximum of points respectively, then from this N
M the point that amplitude is maximum is chosen in individual；
Step 3, sequentially in time by M some sequence, if that is first A of certain continuous three point, second
Point B, thirdly C meet the time interval time interval less than B to C of following condition: A to B, and
The heart rate value R drawn in the time interval of A to C in scope [30,240], then retains this R value, wherein
R=(d*60)/f, d represents AC time interval, and f represents cardiechema signals sample frequency；
Step 4, by all heart rate value R of obtaining in step 3 in magnitude order, takes after sequence in sequence
Median and two numbers adjacent with this median left and right, and draw median and two adjacent values around
Meansigma methods, this meansigma methods is heart rate.
Further technical scheme is, in described step 3, sequentially in time by M some sequence, if
That is first A of certain continuous three point, second point B, thirdly C are unsatisfactory for the time interval of following condition: A to B
Less than the time interval of B to C, and the heart rate value R drawn in the time interval of A to C is in scope [30,240]
In, then heart rate value R=(d*60)/f, d represent M put in time interval between adjacent 2, f
Represent cardiechema signals sample frequency, then be retained in the heart rate value R in scope [30,240].
Further technical scheme is, if after sequence, the median in sequence does not has and it in described step 4
Left neighbour or the number of right neighbour, then, when calculating meansigma methods, do not count the number of this left neighbour or right neighbour.
This method is based primarily upon following three temporal signatures, the cardiechema signals of first, heart beat cycle, leads to
Often the signal amplitude of its S1, S2 position is maximum, and second, S2 are less than with the time interval of previous S1
S2 and the time interval of later S1.3rd, heart rate value is usually located between 30 to 240.Then in conjunction with
Three above temporal signatures, then can preferably orient the position of multiple S1, S2 in one section of cardiechema signals, thus
Calculate heart rate value.
Compared with prior art, the invention has the beneficial effects as follows:
This method saves envelope extraction and merges splitting of heart sounds process, directly determines from heart sound temporal signatures
Heart sound segmentation, and utilize median method to reduce error.And then make the time cost calculating heart rate less.
Detailed description of the invention
A kind of rate calculation method based on cardiechema signals temporal signatures, it comprises the following steps:
Step 1, obtains cardiechema signals
Utilizing electronic stethoscope to gather cardiechema signals, acquisition time is 6 seconds；
Step 2, extracts local maximum point
6 seconds long cardiechema signals are divided into 48 sections, extract each section of interior maximum of points respectively, then from this 48
Choose 12 points of amplitude maximum individual o'clock；
Step 3, extracts the point meeting heart rate feature
Sequentially in time by the sequence of 12 points, if that is first A of certain continuous three point, second point B,
Thirdly C meets the time interval time interval less than B to C of following condition: A to B, and A to C
Time interval in the heart rate value R that draws in scope [30,240], then retain this R value, wherein R=(d
* 60)/f, d represent AC time interval, and f represents cardiechema signals sample frequency；
Step 4, calculates heart rate
By all heart rate value R of obtaining in step 3 in magnitude order, the median in sequence after sequence is taken
And two numbers adjacent with this median left and right, and draw the meansigma methods of median and two adjacent values around,
This meansigma methods is heart rate；
In step 3, sequentially in time by M some sequence, if that is first A of certain continuous three point, the
2 B, thirdly C are unsatisfactory for the time interval time interval less than B to C of following condition: A to B,
And the heart rate value R drawn in the time interval of A to C is in scope [30,240], then heart rate value R=(d
* 60)/f, d represent the time interval in M point between adjacent 2, and f represents cardiechema signals sampling frequency
Rate, then it is retained in the heart rate value R in scope [30,240]；
If not neighbour left with it or the number of right neighbour of the median in sequence after sequence, then calculate flat in step 4
During average, the number of this left neighbour or right neighbour is not counted.
The essence of the present invention is described in detail by above detailed description of the invention, but can not be to the guarantor of the present invention
The scope of protecting limits, it should be apparent that, under the enlightenment of the present invention, the art those of ordinary skill
Many improvement and modification can also be carried out, it should be noted that these improve and modify all to fall the present invention's
Within claims.
Claims (3)
1. a rate calculation method based on cardiechema signals temporal signatures, it is characterised in that it includes following step
Rapid:
Step 1, obtains cardiechema signals；
Step 2, is divided into N section by cardiechema signals, extracts each section of interior maximum of points respectively, then from this N
M the point that amplitude is maximum is chosen in individual；
Step 3, sequentially in time by M some sequence, if that is first A of certain continuous three point, second
Point B, thirdly C meet the time interval time interval less than B to C of following condition: A to B,
And the heart rate value R drawn in the time interval of A to C is in scope [30,240], then retain this R
Value, wherein R=(d*60)/f, d represents AC time interval, and f represents cardiechema signals sample frequency；
Step 4, by all heart rate value R of obtaining in step 3 in magnitude order, takes after sequence in sequence
Median and two numbers adjacent with this median left and right, and draw median and two around
The meansigma methods of adjacent value, this meansigma methods is heart rate.
Rate calculation method based on cardiechema signals temporal signatures the most according to claim 1, it is characterised in that
In described step 3, sequentially in time by M some sequence, if certain continuous three points that is first point
A, second point B, thirdly C are unsatisfactory for the time interval of following condition: A to B less than B to C's
The heart rate value R drawn in time interval, and the time interval of A to C is interior in scope [30,240],
Time interval between adjacent 2 during then heart rate value R=(d*60)/f, d represents M point,
F represents cardiechema signals sample frequency, then is retained in the heart rate value R in scope [30,240].
Rate calculation method based on cardiechema signals temporal signatures the most according to claim 1, it is characterised in that
If not neighbour left with it or the number of right neighbour of the median in sequence after sequence, then count in described step 4
When calculating meansigma methods, the number of this left neighbour or right neighbour is not counted.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610300369.3A CN105997036A (en) | 2016-05-09 | 2016-05-09 | A heart rate calculating method based on the heart sound signal time domain features |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610300369.3A CN105997036A (en) | 2016-05-09 | 2016-05-09 | A heart rate calculating method based on the heart sound signal time domain features |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105997036A true CN105997036A (en) | 2016-10-12 |
Family
ID=57099907
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610300369.3A Pending CN105997036A (en) | 2016-05-09 | 2016-05-09 | A heart rate calculating method based on the heart sound signal time domain features |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105997036A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109998514A (en) * | 2019-04-25 | 2019-07-12 | 四川长虹电器股份有限公司 | A kind of rate calculation method of irregular cardiechema signals |
| CN110123370A (en) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | The rate calculation method of first heart sound part omission cardiechema signals |
-
2016
- 2016-05-09 CN CN201610300369.3A patent/CN105997036A/en active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109998514A (en) * | 2019-04-25 | 2019-07-12 | 四川长虹电器股份有限公司 | A kind of rate calculation method of irregular cardiechema signals |
| CN110123370A (en) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | The rate calculation method of first heart sound part omission cardiechema signals |
| CN110123370B (en) * | 2019-05-30 | 2021-06-29 | 四川长虹电器股份有限公司 | Heart rate calculation method for missing heart sound signal of first heart sound part |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| De Carvalho et al. | Development of a Matlab software for analysis of heart rate variability | |
| CN109124610B (en) | Anti-interference method and device for non-invasive blood pressure measurement | |
| US20110257548A1 (en) | Method and system for processing heart sound signals | |
| CN103845079B (en) | A kind of detection method of the Doppler's fetal heart sound instantaneous heart rate based on blind separation | |
| WO2005099562A1 (en) | Non-invasive measurement of second heart sound components | |
| CN106073783B (en) | Method for extracting respiration rate from photoplethysmography wave | |
| CN104622440B (en) | The method and device of punctuate during a kind of extraction pulse wave | |
| US20230015562A1 (en) | Device and process for ecg measurements | |
| Xie et al. | Heart rate estimation from ballistocardiography based on hilbert transform and phase vocoder | |
| CN105997036A (en) | A heart rate calculating method based on the heart sound signal time domain features | |
| CN104382589B (en) | Fetal electrocardiogram separation extraction method based on partial resampling by segments | |
| JP2820721B2 (en) | QRS detector | |
| Yang et al. | Removal of pulse waveform baseline drift using cubic spline interpolation | |
| CN108814642B (en) | Heart sound positioning and heart rate calculating method of electronic stethoscope | |
| WO2018023698A1 (en) | Fetal-electrocardiogram separation method and device | |
| Ahlstrom et al. | A method for accurate localization of the first heart sound and possible applications | |
| CN110491504A (en) | A kind of acquisition methods of cardiechema signals medical guidelines data | |
| Hikmah et al. | A signal processing framework for multimodal cardiac analysis | |
| CN111528900A (en) | Heart sound segmentation method and device based on Butterworth filter and Shannon entropy method | |
| Varghees et al. | Two-channel heart sound segmentation framework using phonocardiogram and pulsatile signals | |
| Ari et al. | On a robust algorithm for heart sound segmentation | |
| Liu et al. | Algorithm of heart sound segmentation and feature extraction | |
| CN114366121A (en) | Double-dynamic threshold electrocardio QRS wave detection method based on digital morphology | |
| Wang et al. | Continuous blood pressure prediction using pulse features and Elman neural networks | |
| CN109998514A (en) | A kind of rate calculation method of irregular cardiechema signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161012 |