CN105147269B - A kind of noninvasive continuous BP measurement method - Google Patents
A kind of noninvasive continuous BP measurement method Download PDFInfo
- Publication number
- CN105147269B CN105147269B CN201510335331.5A CN201510335331A CN105147269B CN 105147269 B CN105147269 B CN 105147269B CN 201510335331 A CN201510335331 A CN 201510335331A CN 105147269 B CN105147269 B CN 105147269B
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- point
- blood pressure
- electrocardiosignal
- value
- 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
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The present invention relates to a kind of noninvasive continuous BP measurement method, comprise the following steps：Gather the electrocardio of subject and the synchronizing signal of pulse wave；The synchronizing signal of electrocardio and pulse wave to subject is smoothed, and detection R wave of electrocardiosignal peak value point and pulse wave characteristic point carry out position mark, calculates the conduction ratios of the conduction time and corresponding pulse wave that obtain pulse wave；The conduction time and conduction ratios of gained pulse wave are brought into personalized blood pressure parameter model, the continuous systolic pressure for obtaining subject is calculated.The present invention introduces new parameter conduction ratios in existing non-invasive blood pressure model, establishes the personalizing parameters model for each subject, can obtain more accurate blood pressure result.This method is simple to operate, and accuracy is high, is that clinical monitoring of blood pressure, domestic monitoring and early warning and the monitoring of depressor provide important evidence.
Description
Technical field
The present invention relates to a kind of noninvasive continuous blood pressure assay method, and in particular to one kind is based on pulse wave translation time and biography
Lead the continuous BP measurement method of ratio.
Background technology
Angiosthenia is the important physiological parameter for reflecting human circulation function, while being also clinically to diagnose the illness, observe and control
Therapeutic effect, the important evidence for carrying out Index for diagnosis.The blood pressure that noninvasive continuous BP measurement can be monitored in each cardiac cycle becomes
Change, therefore noninvasive continuous BP measurement has ten in clinical monitoring of blood pressure, domestic monitoring and early warning and the monitoring of depressor
Divide important meaning.
In recent years, the noninvasive continuous blood of blood pressure was calculated using the conduction time (PWTT) or conduction of velocity (PWV) of pulse wave
Pressure measuring method is increasingly becoming each scholar's primary study direction.The conduction time of pulse wave is that arterial pulse wave is passed from aorta petal
The time required for tip is multicast to, when calculating generally using ecg-r wave crest as starting point, synchronous pulse wave characteristic point is eventually
Point.Conduction of velocity is the spread speed of pulse wave in the artery.Obtaining the electrocardio of subject and pulse wave synchrodata
On the basis of, it is possible to achieve for the reckoning of blood pressure.
Existing result of study shows that systolic pressure is linear relationship with conduction time PWTT, and sets up model SBP=b1+
B2*PWTT can realize noninvasive continuous blood pressure measuring.But the model is the blood pressure parameter unified model set up by different subjects, its
Accuracy is vulnerable to the influence of individual difference.
The content of the invention
In order to overcome above-mentioned technological deficiency, the present invention provides a kind of new noninvasive continuous BP measurement method, i.e., existing
New parameter-conduction ratios (dc) are introduced on the basis of blood pressure parameter model, while individually setting up personalization for each subject
Blood pressure parameter model, so that the accuracy of continuous BP measurement is significantly improved, and this method is simple and efficient, it is easy to operate, be conducive to
It is widely popularized.
To achieve these goals, the present invention is adopted the following technical scheme that：
A kind of continuous BP measurement method, comprises the following steps：
(1) synchronizing signal of the electrocardio of collection subject and pulse wave；
(2) electrocardio and the synchronizing signal of pulse wave to subject is smoothed, and detects R wave of electrocardiosignal peak value point
Position mark is carried out with pulse wave characteristic point, the conduction ratio of the conduction time and corresponding pulse wave that obtain pulse wave is calculated
Rate；
(3) conduction time and conduction ratios of gained pulse wave are brought into personalized blood pressure parameter model, calculating is obtained
The continuous systolic pressure of subject.
The present invention introduces pulse transit ratio parameter in existing continuous BP measurement method, and set up for it is each by
The personalized blood pressure parameter model of examination person, the accuracy of measurement can be improved by methods described.
In measuring method of the present invention, the detection method of the ecg-r wave wave crest point is as follows：In collection electrocardiosignal
First two seconds, the slope minimum value of the segment signal is obtained using algorithm, take the numerical value 0.6 times is threshold value thr, finds what is collected
Less than threshold value, itself slope, more than the point of threshold value, is starting point by the point to previous slope in electrocardiosignal, forward 20 points
Interior maximizing point, as ecg-r wave wave crest point.
In measuring method of the present invention, the pulse wave characteristic point is main crest rising edge greatest gradient point.The arteries and veins
The detection method of wave characteristic of fighting point is as follows：With the first ecg-r wave ripple detected in above-mentioned ecg-r wave crest point detecting method
Peak dot is starting point, and latter ecg-r wave wave crest point is terminal, and the greatest gradient for finding out pulse wave in this period clicks through rower
Note, as pulse wave characteristic point.
In measuring method of the present invention, the pulse wave translation time (PWTT) is pulse wave characteristic point and ecg-r wave
The alternate position spike of wave crest point.
In measuring method of the present invention, the pulse transit ratio (DC) is pulse wave translation time and cardiac electrical cycle
Ratio, i.e. DC=PWTT/T；Wherein cardiac electrical cycle (T) is time of first R ripples crest to secondary R wave crest location.
In measuring method of the present invention, the conduction ratios of the pulse wave account for the cardiac electrical cycle for pulse wave translation time
Ratio.As shown in Figure 1.Cardiac electrical cycle (T) is time of first R ripples crest to secondary R wave crest location.
In measuring method of the present invention, the personalized blood pressure parameter model is set up by following methods：
(1) synchronizing signal of the electrocardio of collection subject and pulse wave, is smoothed, and detects R wave crest of electrocardiosignal
Value point and pulse wave characteristic point carry out position mark, calculate the conduction of the conduction time and corresponding pulse wave that obtain pulse wave
Ratio；
(2) the interruption blood pressure of subject is measured, by interruption pressure value, the conduction time of gained pulse wave and conduction ratios generation
Enter in basic model SBP=b1+b2 × PWTT+b3 × DC, obtain parameter b1, b2, b3, so as to set up personalized blood pressure parameter mould
Type.
In measuring method of the present invention, the smoothing processing is the conventional signal processing method smothing filtering in this area
Method.
The present invention introduces new parameter-conduction ratios in existing non-invasive blood pressure model, establishes for each tested
The personalizing parameters model of person, can obtain more accurate blood pressure result.This method is simple to operate, and accuracy is high, is clinical
Monitoring of blood pressure, domestic monitoring and early warning and the monitoring of depressor provide important evidence.
Brief description of the drawings
Fig. 1 is the calculating schematic diagram of pulse wave translation time (PWTT).
Embodiment
Following examples are used to illustrate the present invention, but are not limited to the scope of the present invention.
A kind of noninvasive continuous BP measurement method of embodiment 1
Step is as follows：
(1) synchronizing signal of the electrocardio of collection subject and pulse wave；
(2) electrocardio and the synchronizing signal of pulse wave to subject is smoothed, and detects R wave of electrocardiosignal peak value point
Gone forward side by side line position tagging with pulse wave characteristic point, the conduction time that obtains pulse wave is calculated by the two alternate position spike and corresponding
The conduction ratios of pulse wave；
(3) conduction time and conduction ratios of gained pulse wave are brought into personalized blood pressure parameter model, calculating is obtained
The continuous systolic pressure of subject.
Wherein, the personalized blood pressure parameter model is set up by following methods：
(1) blood pressure of subject is measured, and gathers the electrocardio of subject and the synchronizing signal of pulse wave, is smoothly located
Reason, detection R wave of electrocardiosignal peak value point and pulse wave characteristic point carry out position mark, calculating obtain the conduction time of pulse wave with
And the conduction ratios of corresponding pulse wave；Repeat experiment at least three times；
(2) by gained pressure value, the conduction time of pulse wave and conduction ratios substitute into basic model SBP=b1+b2 ×
In PWTT+b3 × DC, parameter b1, b2, b3 are obtained, personalized blood pressure parameter model is set up.
The detection method of the ecg-r wave wave crest point is as follows：At first two seconds of collection electrocardiosignal, it is somebody's turn to do using algorithm
The slope minimum value of segment signal, take the numerical value 0.6 times is threshold value thr, finds previous point in the electrocardiosignal collected oblique
Less than threshold value, itself slope, more than the point of threshold value, is starting point by the point to rate, forward maximizing point in 20 points, the as heart
Electric R ripples wave crest point.
The pulse wave characteristic point is main crest rising edge greatest gradient point.The detection method of the pulse wave characteristic point is such as
Under：The first ecg-r wave wave crest point detected using in above-mentioned ecg-r wave crest point detecting method is starting point, latter electrocardio R
Ripple wave crest point is terminal, and the greatest gradient point for finding out pulse wave in this period is marked, as pulse wave characteristic point.
The pulse transit ratio (DC) is pulse wave translation time and the ratio of cardiac electrical cycle, i.e. DC=PWTT/T；
Wherein cardiac electrical cycle (T) is time of first R ripples crest to secondary R wave crest location.
Method validation is tested
Pass through exercise testing measuring method more of the present invention and the accuracy of traditional noninvasive continuous BP measurement method.
The basic model of the noninvasive continuous BP measurement method of tradition is：SBP=b1+b2 × PWTT；The present invention is noninvasive continuous
The basic model of blood pressure measuring method is：SBP=b1+b2 × PWTT+b3 × DC.Exercise testing is divided into be carried out twice, is once
" model experiment " of model is set up, is " examining experiment " for the second time.Quality to two kinds of measuring methods is contrasted.As a result such as table
1, shown in table 2.
Wherein, error calculation formula is：
Table 1 " model experiment " mean error
As shown in Table 1, for each subject, the contraction obtained in " model experiment " using the inventive method is held up
The systolic pressure error that difference is obtained all less than existing measuring method.Wherein, the mean error obtained using the inventive method is averaged
It is worth for 4.7542mmHg, the AME 6.0120mmHg obtained less than existing measuring method.
Table 2 " examining experiment " mean error
As shown in table 2, in " examining experiment ", the systolic pressure error obtained using the method for the invention is most of all small
The systolic pressure error obtained in existing measuring method.Wherein, the mean error average value that the method for the invention is obtained is
5.8646mmHg, the mean error average value 7.0066mmHg obtained less than existing measuring method.
The T assays of table 3
| Group | H | P | Conventional method average value | The inventive method average value |
| " model experiment " | 1 | 0.0076 | 6.0120 | 4.7542 |
| " examining experiment " | 1 | 0.0289 | 7.0066 | 5.8646 |
As shown in table 3, T inspections are carried out using matlab programmings to the result of " model experiment ", as a result (received for H=1
Assuming that), P=0.0076<α=0.05, illustrates that two grouping error discriminations are very high, with significant difference.Pass through " model experiment "
The average value of mean error compares, and the result that T is examined can determine systolic pressure of the method for the invention than existing method
Error is small.Similarly, same conclusion is also obtained for " examine and test ".
Comprehensive " model experiment " and the comparative result of the mean error average value of " examine and test ", and T inspections twice
As a result, we may safely draw the conclusion：Measuring method of the present invention compares the measuring method of existing single parameter in accuracy side
Face obtains very big raising.
Although above the present invention is described in detail with a general description of the specific embodiments,
On the basis of the present invention, it can be made some modifications or improvements, this will be apparent to those skilled in the art.Cause
This, these modifications or improvements, belong to the scope of protection of present invention without departing from theon the basis of the spirit of the present invention.
Claims (9)
1. a kind of personalized blood pressure parameter method for establishing model, it is characterised in that comprise the following steps：
(1) synchronizing signal of the electrocardio of collection subject and pulse wave, is smoothed, and detects R wave of electrocardiosignal peak value point
Position mark is carried out with pulse wave characteristic point, the conduction ratio of the conduction time and corresponding pulse wave that obtain pulse wave is calculated
Rate；
(2) blood pressure of subject is measured, pressure value, the conduction time of gained pulse wave and conduction ratios are substituted into model SBP=
In b1+b2 × PWTT+b3 × DC, parameter b1, b2, b3 are obtained, it is determined that personalized blood pressure parameter model；Wherein, SBP represents blood pressure
Value；PWTT represents the conduction time of pulse wave；DC represents the conduction ratios of pulse wave；
The conduction ratios of the pulse wave are pulse wave translation time and the ratio of cardiac electrical cycle.
2. personalized blood pressure parameter method for establishing model according to claim 1, it is characterised in that described in step (1)
The detection method of R wave of electrocardiosignal peak value point is as follows：In two seconds after collection electrocardiosignal starts, the section is obtained using algorithm
The slope minimum value of signal, take the numerical value 0.6 times is threshold value thr, finds previous slope in the electrocardiosignal collected
Less than threshold value, itself slope is more than the point of threshold value, is starting point by the point, forward maximizing point, as electrocardio in 20 points
Signal R crest value points.
3. personalized blood pressure parameter method for establishing model according to claim 1 or 2, it is characterised in that in step (1),
The detection method of the pulse wave characteristic point is as follows：It is latter using the R ripples crest value of first electrocardiosignal that detects as starting point
The R ripples crest value of individual electrocardiosignal is terminal, and the greatest gradient point for finding out pulse wave in this period is marked, as arteries and veins
Wave characteristic of fighting point.
4. personalized blood pressure parameter method for establishing model according to claim 1 or 2, it is characterised in that in step (1),
The conduction time of the pulse wave is pulse wave characteristic point and the alternate position spike of R wave of electrocardiosignal peak value point.
5. personalized blood pressure parameter method for establishing model according to claim 3, it is characterised in that described in step (1)
The conduction time of pulse wave is pulse wave characteristic point and the alternate position spike of R wave of electrocardiosignal peak value point.
6. the personalized blood pressure parameter method for establishing model according to claim 1,2 or 5, it is characterised in that the electrocardio
Cycle is time interval of the n-th R ripples crest to the N+1 R ripple crest location.
7. personalized blood pressure parameter method for establishing model according to claim 3, it is characterised in that the cardiac electrical cycle is
Time interval of the n-th R ripples crest to the N+1 R ripple crest location.
8. personalized blood pressure parameter method for establishing model according to claim 4, it is characterised in that the cardiac electrical cycle is
Time interval of the n-th R ripples crest to the N+1 R ripple crest location.
9. personalized blood pressure parameter method for establishing model according to claim 1, it is characterised in that wherein in step (1),
The detection method of the R wave of electrocardiosignal peak value point is as follows：At first two seconds of collection electrocardiosignal, this section letter is obtained using algorithm
Number slope minimum value, it is threshold value thr to take 0.6 times of the numerical value, finds previous slope in the electrocardiosignal collected small
It is more than the point of threshold value in threshold value and itself slope, is starting point by the point, forward maximizing point in 20 points, as electrocardio is believed
Number R crest values point；
The detection method of the pulse wave characteristic point is as follows：Using the R ripples crest value of first electrocardiosignal that detects as starting point,
The R ripples crest value of latter electrocardiosignal is terminal, and the greatest gradient point for finding out pulse wave in this period is marked, i.e.,
For pulse wave characteristic point；
The conduction time of the pulse wave is pulse wave characteristic point and the alternate position spike of R wave of electrocardiosignal peak value point；
The conduction ratios of the pulse wave are pulse wave translation time and the ratio of cardiac electrical cycle；The cardiac electrical cycle is n-th R
Time interval of the ripple crest to the N+1 R ripple crest location.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510335331.5A CN105147269B (en) | 2015-06-16 | 2015-06-16 | A kind of noninvasive continuous BP measurement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510335331.5A CN105147269B (en) | 2015-06-16 | 2015-06-16 | A kind of noninvasive continuous BP measurement method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105147269A CN105147269A (en) | 2015-12-16 |
| CN105147269B true CN105147269B (en) | 2017-10-13 |
Family
ID=54788582
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510335331.5A Active CN105147269B (en) | 2015-06-16 | 2015-06-16 | A kind of noninvasive continuous BP measurement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105147269B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107788965A (en) * | 2016-09-05 | 2018-03-13 | 京东方科技集团股份有限公司 | A kind of determination method and device of blood pressure |
| CN106725400A (en) * | 2016-11-24 | 2017-05-31 | 南昌大学 | A kind of Novel blood-pressure meter for merging electrocardiosignal and impulse wave form qualitative assessment |
| CN106974631B (en) * | 2017-03-17 | 2020-01-07 | 深圳大学 | Blood pressure measuring method and device based on pulse wave waveform and electrocardiosignal |
| CN107126201A (en) * | 2017-03-31 | 2017-09-05 | 悦享趋势科技（北京）有限责任公司 | Continuous blood pressure detection method, equipment and the device of non-invasive |
| WO2019019119A1 (en) * | 2017-07-27 | 2019-01-31 | Vita-Course Technologies (Hainan) Co., Ltd. | Systems and methods for determining blood pressure of subject |
| CN107854123B (en) * | 2017-12-13 | 2021-02-02 | 清华大学深圳研究生院 | Sleeveless continuous blood pressure monitoring method and device |
| CN110652318B (en) * | 2019-07-19 | 2022-09-13 | 飞依诺科技股份有限公司 | Measurement method and system for obtaining arteriosclerosis index based on ultrasonic equipment |
| CN112870543B (en) * | 2021-02-26 | 2022-06-21 | 中山大学 | Drug release control device, control method thereof, and computer-readable storage medium |
| CN113288092B (en) * | 2021-05-06 | 2022-08-19 | 福建师范大学 | Blood pressure detection method and terminal for extracting pulse wave based on video |
| CN113456043B (en) * | 2021-07-08 | 2023-05-26 | 军事科学院系统工程研究院卫勤保障技术研究所 | Continuous blood pressure detection method and device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6616613B1 (en) * | 2000-04-27 | 2003-09-09 | Vitalsines International, Inc. | Physiological signal monitoring system |
| CN1849998A (en) * | 2006-05-26 | 2006-10-25 | 中国人民解放军空军航空医学研究所 | Method and apparatus for continuously measuring blood pressure |
| CN101288587A (en) * | 2008-04-17 | 2008-10-22 | 天津同阳科技发展有限公司 | Watch strap human blood pressure non-invasive continuous detection device |
| CN101327121A (en) * | 2007-06-22 | 2008-12-24 | 香港中文大学 | Physiological parameter measurement mechanism |
| CN103393415A (en) * | 2013-08-20 | 2013-11-20 | 王卫东 | Method for measuring continuous changing blood pressure |
| CN103637788A (en) * | 2013-12-02 | 2014-03-19 | 清华大学 | Real-time blood pressure measuring device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7544168B2 (en) * | 2004-09-30 | 2009-06-09 | Jerusalem College Of Technology | Measuring systolic blood pressure by photoplethysmography |
| CN105377137A (en) * | 2013-06-28 | 2016-03-02 | 株式会社村田制作所 | Biological state-estimating device |
-
2015
- 2015-06-16 CN CN201510335331.5A patent/CN105147269B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6616613B1 (en) * | 2000-04-27 | 2003-09-09 | Vitalsines International, Inc. | Physiological signal monitoring system |
| CN1849998A (en) * | 2006-05-26 | 2006-10-25 | 中国人民解放军空军航空医学研究所 | Method and apparatus for continuously measuring blood pressure |
| CN101327121A (en) * | 2007-06-22 | 2008-12-24 | 香港中文大学 | Physiological parameter measurement mechanism |
| CN101288587A (en) * | 2008-04-17 | 2008-10-22 | 天津同阳科技发展有限公司 | Watch strap human blood pressure non-invasive continuous detection device |
| CN103393415A (en) * | 2013-08-20 | 2013-11-20 | 王卫东 | Method for measuring continuous changing blood pressure |
| CN103637788A (en) * | 2013-12-02 | 2014-03-19 | 清华大学 | Real-time blood pressure measuring device |
Non-Patent Citations (1)
| Title |
|---|
| 基于个性化脉搏波传导参数的连续血压测量方法研究;董骁等;《中国医疗设备》;20141031;第29卷(第10期);24-27 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105147269A (en) | 2015-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105147269B (en) | A kind of noninvasive continuous BP measurement method | |
| Miao et al. | Multi-sensor fusion approach for cuff-less blood pressure measurement | |
| CN102270264B (en) | Physiological signal quality evaluation system and method | |
| CN104382571B (en) | A kind of measurement blood pressure method and device based on radial artery pulse wave conduction time | |
| WO2019161609A1 (en) | Method for analyzing multi-parameter monitoring data and multi-parameter monitor | |
| US20100280396A1 (en) | System for Cardiac Pathology Detection and Characterization | |
| CN103565427B (en) | Paracycle physiological signal characteristic point detection | |
| CN103479383B (en) | Device for analyzing heart sound signals, and intelligent heart stethoscope provided with device for analyzing heart sound signals | |
| CN107854123B (en) | Sleeveless continuous blood pressure monitoring method and device | |
| CN106725396A (en) | A kind of noninvasive blood pressure measuring device by shooting and measuring method based on double pulse waves | |
| CN108354597A (en) | A kind of rapid blood pressure computational methods based on the extraction of optimal wave | |
| CN104523262A (en) | Sleep quality detection method based on electrocardiosignals | |
| CN107692988A (en) | The method of monitoring of blood pressure method and graphic software platform based on intelligent platform | |
| CN107145697A (en) | The treating method and apparatus of blood oxygen saturation | |
| CN107049289A (en) | Based on radial artery ripple without cuff sphygmomanometer and measuring method | |
| CN110897631A (en) | Real-time pregnancy monitoring device and method | |
| CN109938719A (en) | A kind of Driver Fatigue Detection based on physiological parameter | |
| CN105326477B (en) | A kind of breathing and heartbeat synchronization method for extracting signal and device for plurality of medical signal acquisition | |
| Foroozan et al. | Robust beat-to-beat detection algorithm for pulse rate variability analysis from wrist photoplethysmography signals | |
| CN104545883A (en) | Electrocardiosignal-based sleep quality detection equipment and service thereof | |
| CN108685605B (en) | Catheter tip positioning method and system | |
| CN108652604A (en) | It is a kind of based on electrocardiosignal without air bag blood pressure detecting method and system | |
| CN111000544B (en) | Method and system for constructing hybrid continuous blood pressure measurement model based on PPG waveform | |
| CN103767694B (en) | Method for accurately extracting cuff pressure shockwave | |
| Hu et al. | Feature extraction and choice in PCG based on Hilbert Transfer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A non-invasive continuous blood pressure measurement method Effective date of registration: 20230717 Granted publication date: 20171013 Pledgee: Jiangsu Jiangyin Rural Commercial Bank Co.,Ltd. high tech Zone sub branch Pledgor: JIANGSU STANDARD MEDICAL TECHNOLOGY Co.,Ltd. Registration number: Y2023980048651 |