CN105147269B - A kind of noninvasive continuous BP measurement method - Google Patents

A kind of noninvasive continuous BP measurement method Download PDF

Info

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
conduction
Prior art date
Application number
CN201510335331.5A
Other languages
Chinese (zh)
Other versions
CN105147269A (en
Inventor
果淑颖
Original Assignee
江苏斯坦德利医疗科技有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 江苏斯坦德利医疗科技有限公司 filed Critical 江苏斯坦德利医疗科技有限公司
Priority to CN201510335331.5A priority Critical patent/CN105147269B/en
Publication of CN105147269A publication Critical patent/CN105147269A/en
Application granted granted Critical
Publication of CN105147269B publication Critical patent/CN105147269B/en

Links

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

A kind of noninvasive continuous BP measurement method

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.
CN201510335331.5A 2015-06-16 2015-06-16 A kind of noninvasive continuous BP measurement method CN105147269B (en)

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 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 (4)

* Cited by examiner, † Cited by third party
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

Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
EP3015067A4 (en) * 2013-06-28 2017-03-29 Murata Manufacturing Co., Ltd. Biological state-estimating device

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
基于个性化脉搏波传导参数的连续血压测量方法研究;董骁等;《中国医疗设备》;20141031;第29卷(第10期);24-27 *

Also Published As

Publication number Publication date
CN105147269A (en) 2015-12-16

Similar Documents

Publication Publication Date Title
WO2017097073A1 (en) Method and system for calculating blood vessel pressure difference and fractional flow reserve
Shao et al. Noncontact monitoring breathing pattern, exhalation flow rate and pulse transit time
US20140296727A1 (en) Detection of Parameters in Cardiac Output Related Waveforms
KR101068116B1 (en) Apparatus and method for sensing radial arterial pulses for noninvasive and continuous measurement of blood pressure
EP3305188A1 (en) Biosignal measurement device utilizing plurality of electrodes for measuring biosignal as touch sensors
Kurylyak et al. A Neural Network-based method for continuous blood pressure estimation from a PPG signal
CN102008296B (en) Device and method for measuring arterial blood pressures based on pulse wave signals and electrocardiosignals
US9028421B2 (en) Blood flow image diagnosing device
US8591428B2 (en) Method and apparatus for measuring blood volume
US8469895B2 (en) Deriving central aortic systolic pressure and analyzing arterial waveform data to derive central aortic systolic pressure values
US20120071767A1 (en) Pulmonary artery pressure estimator
US7309313B2 (en) Vascular disease examining system and bypass vascular diagnosing device
US20150141765A1 (en) Method and system to identify motion data associated with consistent electrical and mechanical behavior for a region of interest
JP2013517908A (en) Eliminating the effects of irregular cardiac cycles in determining cardiovascular parameters
JP5689116B2 (en) Method and apparatus for detecting and evaluating reactive hyperemia using segmental plethysmography
US5025795A (en) Non-invasive cardiac performance monitoring device and method
CN103549950B (en) Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring
US20090312648A1 (en) Adaptive Medical Image Acquisition System
CN102869303B (en) Visualization of myocardial infarct size in diagnostic ECG
US20030153837A1 (en) Non-invasive determination of left-ventricular pressure
US20090259266A1 (en) System for Heart monitoring, Characterization and Abnormality Detection
CN101357062B (en) Vital signs parameter evaluation device based on volume pulsation wave detection
US20130116580A1 (en) System for quality assessment of physiological signals and method thereof
Mase et al. Feasibility of cuff-free measurement of systolic and diastolic arterial blood pressure
CN103565429A (en) System and method for observing heart rate of passenger

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