CN105342590A - Blood pressure enveloping estimation method based on oscillography - Google Patents
Blood pressure enveloping estimation method based on oscillography Download PDFInfo
- Publication number
- CN105342590A CN105342590A CN201510732538.6A CN201510732538A CN105342590A CN 105342590 A CN105342590 A CN 105342590A CN 201510732538 A CN201510732538 A CN 201510732538A CN 105342590 A CN105342590 A CN 105342590A
- Authority
- CN
- China
- Prior art keywords
- value
- threshold
- pressure
- blood pressure
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention provides a blood pressure enveloping estimation method based on oscillography. The blood pressure enveloping estimation method comprises the following steps: firstly selecting a dynamic pressure sequence within one range after the measurement, defining a fitting function and setting initial values of parameter vectors, seeking the parameter vectors enabling the errors to be the minimum by adopting the least square principle of a levenberg marquardt algorithm, and finally determining an enveloping curve. The invention provides the blood pressure enveloping estimation method based on the oscillography which is small in error and high in precision.
Description
Technical field
The present invention relates to human blood-pressure and measure medical domain, particularly relate to a kind of method of estimation of the blood pressure envelope based on oscillographic method.
Background technology
Blood pressure is one of physiological parameter of wanting of body weight for humans, and the method for Measure blood pressure is mainly divided into oscillographic method and Korotkoff's Sound method two class.Current most noinvasive electric sphygmomanometer all adopts oscillographic method.Electric sphygmomanometer is converted to electronic signal by the pressure transducer on cuff the pressure oscillation with pace, judges high hypotension by the amplitude characteristic of electronic signal.In the measuring process of oscillographic method, cuff slowly pressurizes.When cuff pressure is greater than diastolic pressure, pressure wave amplitude starts to increase.Along with cuff pressure rises to mean pressure, pressure wave amplitude is maximum.When cuff pressure is greater than mean pressure, pressure wave amplitude reduces gradually.When cuff pressure is greater than systolic pressure, pressure wave amplitude is reduced to minima.
In oscillographic method, cuff pressure is called static pressure, pressure wave is called dynamic pressure.In oscillographic method measuring process, sixty-four dollar question is the envelope curve going out pressure wave according to pressure wave sequential extraction procedures, and sets up the corresponding relation of static pressure and dynamic pressure envelope.After setting up static pressure and dynamic pressure envelope relation, finally can determine the high voltage/low voltage value of blood pressure according to the amplitude characteristic ratios method in oscillographic method technology or S method.
In actual application, dynamic pressure envelope has multiple acquiring method.Such as fitting of a polynomial, gaussian curve approximation etc.Publication number is CN102018507A application for a patent for invention, and disclosed " a kind of data analysis processing method based on oscillographic method Measure blood pressure ", just make use of the method for fitting of a polynomial, definition curve P (x)=a
0+ a
1x+...+a
nx
n, wherein a0, a1 ..., an is the coefficient needing to ask, n >=2.Dynamic pressure envelope curve is similar to bell shaped curve, and polynomial curve can not simulate bell dynamic pressure envelope when selecting lower order polynomial expressions in effective blood pressure span.Polynomial fitting curve characteristic can not be consistent with dynamic pressure envelope in effective interval.As shown in Figure 3, the curve that L1 representative polynomial fits to, can find out that fitting of a polynomial accurately can not estimate the maximum of dynamic pressure envelope.According to oscillographic method principle, the estimation of diastolic pressure and systolic pressure all needs the maximum using dynamic pressure envelope.Fitting of a polynomial envelope curve two ends trend and theory inconsistent.Dynamic pressure amplitude can approach a fixed value in envelope both sides in theory.
Although and Gauss curve fitting curve can ensure curve characteristic and dynamic pressure envelope curve basically identical.Ascending because the amplitude of cuff pulse wave produced in actual deflation course changes trend, descending again after arriving amplitude extremum, form asymmetric and nonlinear curve envelope trend.Due to the symmetry of Gaussian curve, can not be satisfied in actual measurement process.If use the actual dynamic pressure envelope of gaussian curve approximation, then the variance of actual measured results can be caused to become large.And the maximum to dynamic pressure envelope curve of Gauss curve fitting estimates to there is very important deviation with true maximum.As shown in Figure 3, L2 represents the curve that gaussian curve approximation becomes, and can see from figure, and Gauss curve fitting is more effective than fitting of a polynomial in high-pressure section curve fitting effect.But at the curve of low-pressure section Gauss curve fitting still higher than actual measured value.
Summary of the invention
Goal of the invention of the present invention is the limitation solving above-mentioned existing blood pressure envelope evaluation method, provides the method for estimation of the blood pressure envelope based on oscillographic method that a kind of error is little, precision is high.
The method of estimation of a kind of blood pressure envelope based on oscillographic method of the present invention, comprises the following steps:
S1, after measurement terminates, obtains the static pressure of dynamic pressure peak value and correspondence by computing, two threshold values of the scope of setting dynamic pressure amplitude are y
t1and y
t2, choose dynamic pressure sequences y
i(y
t1≤ y
i≤ y
t2), the static pressure sequence that dynamic pressure sequence pair is answered is x
i, corresponding functional relationship is y
i(x
i);
S2, if the function of matched curve is f (x
i, β), parameter vector β=(a
0, a
1..., a
i..., a
n, b, c), wherein a
0, a
1..., a
i..., a
n, b, c be parameter of curve variable, N is natural number, 0≤i≤N; With dynamic pressure sequences y selected in S1
ibased on, calculate iteration initial value a
0, a
1..., a
i..., a
n, b, c; Wherein, described a
0=y
m, y
mfor dynamic pressure sequences y
iin maximum; a
i=0,0 < i≤N; B=c
2; C=x
m, x
mfor y
mcorresponding static pressure; The formula of definition matched curve is
m be greater than 0 even number;
S3, perform levenbergmarquardt iterative algorithm and find optimum parameter vector, formula is
λ is each iteration stepping length, calculates β; Calculate
t1 is dynamic pressure sequences y
ithreshold value y
t1time i value, t2 is dynamic pressure sequences y
ithreshold value y
t2time i value;
S4, definition threshold value threshold and threshold_1, judge as | θ | < threshold or | E (β+θ)-E (β) | during < threshold_1, iteration terminates, proceed to next step, obtain optimum parameter vector β=(a simultaneously
0, a
1..., a
i..., a
n, b, c);
S5, preserves parameter of curve (a
0, a
1..., a
i..., a
n, b, c), extract envelope f (x
i, β).
Further, described y
t1value be y
i40% of sequence amplitude.
Further, the N value described in step S2 is 1.
Further, the M value described in step S5 is 4.
Further, the value of described threshold is 5-25.
Further, the value of described threshold is 10.
Further, the value of described threshold_1 is initial value (a
0, a
1..., a
i..., a
n, b, c) * 0.005.
After adopting above structure, the present invention compared with prior art, has the following advantages:
Contrast fitting of a polynomial of the prior art and Gauss curve fitting curve, curve fitting trend of the present invention is more consistent with notional result; Also the maximum of more approaching to reality measurement, can estimate systolic pressure and diastolic pressure more accurately according to the method for oscillographic method; Good fitting result can be obtained in low pressure and high-pressure section.Overcome cubic fitting model error greatly, accurate not problem, define the novel curve of a class, with in actual oscillographic method measuring process, extract dynamic pressure envelope, improve accuracy and stability that electric sphygmomanometer measures.
Accompanying drawing explanation
Fig. 1 is flow chart of steps of the present invention.
Fig. 2 is cuff pressure of the present invention and pressure wave curve.
Fig. 3 is fitting of a polynomial of the prior art, Gauss curve fitting curve and curve-fitting results figure of the present invention.
Detailed description of the invention
Below in conjunction with the drawings and specific embodiments, the invention will be further described.
As shown in Figure 1, the method for estimation of a kind of blood pressure envelope based on oscillographic method of the present invention, comprises the following steps:
S1, after measurement terminates, obtains the static pressure of dynamic pressure peak value and correspondence by computing, two threshold values of the scope of setting dynamic pressure amplitude are y
t1and y
t2, choose dynamic pressure sequences y
i(y
t1≤ y
i≤ y
t2), the static pressure sequence that dynamic pressure sequence pair is answered is x
i, corresponding functional relationship is y
i(x
i);
S2, if the function of matched curve is f (x
i, β), parameter vector β=(a
0, a
1..., a
i..., a
n, b, c), wherein a
0, a
1..., a
i..., a
n, b, c be parameter of curve variable, N is natural number, 0≤i≤N; With dynamic pressure sequences y selected in S1
ibased on, calculate iteration initial value a
0, a
1..., a
i..., a
n, b, c; Wherein, described a
0=y
m, y
mfor dynamic pressure sequences y
iin maximum; a
i=0,0 < i≤N; B=c
2; C=x
m, x
mfor y
mcorresponding static pressure; The formula of definition matched curve is
m be greater than 0 even number;
S3, perform levenbergmarquardt iterative algorithm and find optimum parameter vector, formula is
λ is each iteration stepping length, calculates β; Calculate
t1 is dynamic pressure sequences y
ithreshold value y
t1time i value, t2 is dynamic pressure sequences y
ithreshold value y
t2time i value;
S4, definition threshold value threshold and threshold_1, judge as | θ | < threshold or | E (β+θ)-E (β) | during < threshold_1, iteration terminates, proceed to next step, obtain optimum parameter vector β=(a simultaneously
0, a
1..., a
i..., a
n, b, c);
S5, preserves parameter of curve (a
0, a
1..., a
i..., a
n, b, c), extract envelope f (x
i, β).
Described y
t1value be y
i40% of sequence amplitude.
N value described in step S2 is 1.
M value described in step S5 is 4.
The value of described threshold is 5-25.
The value of described threshold is 10.
The value of described threshold_1 is initial value (a
0, a
1..., a
i..., a
n, b, c) * 0.005.
It should be noted that, sequence amplitude gets 40% is herein empirical value.Y of the present invention
t1with y
t2be empirical value, obtain by experiment.Y is set
t1with y
t2mainly to consider the coverage of envelope curve and the impact on matching, also need to consider oscillography ratio juris.The standard chosen is y
t1with y
t2comprise diastolic pressure and systolic pressure, try not to comprise the invalid data outside diastolic pressure and systolic pressure.
In addition, the value of N and M of the present invention, and threshold value threshold is the empirical value tested and obtain.
When the discrete experimental data of an acquired blood pressure measurement, fitting problems is exactly a kind of fundamental problems of application of continuous model that set up correspondence for so discrete measurement point, that be similar to.Because measure the discrete point existence obtained to be disturbed factor, the error between desired result can be reduced by the method being fitted to ideal curve.And in fitting problems, because each discrete point that the matched curve of object module can not be accurately excessively known, only require that matched curve substantially can describe the general outline of known discrete point distribution, and with the error of ideal curve try one's best little.。
The present invention in step s3, adopts levenbergmarquardt algorithm to find the parameter vector making error minimum.Levenbergmarquardt algorithm is a kind of attainable method of nonlinear least square method, and object is the theoretical value that program can be made within the regular hour to calculate method of least square calculate.The principle of least square is by the optimal parameter coupling of the quadratic sum searching data of minimum error.That is, the value of optimum parameter vector β is searched out.
Fig. 3 is fitting of a polynomial of the prior art, Gauss curve fitting curve carries out the curve-fitting results figure of the method for estimation of a kind of blood pressure envelope based on oscillographic method of testing with using the present invention.The vertical coordinate of Fig. 3 represents target blood pressure amplitude, and abscissa represents static pressure.In the experimental result shown in Fig. 3, curve L3 processes the envelope curve obtained.From in figure, envelope curve L3 of the present invention compared to polynomial fitting curve L1 and Gauss curve fitting curve L2, more close to theoretical value in trend; Further, curve of the present invention is more close to actual maximum; And no matter can obtain good fitting result in low pressure and high-pressure section.
Below only just most preferred embodiment of the present invention is described, but can not be interpreted as it is limitations on claims.The present invention is not limited only to above embodiment, and all various changes done in the protection domain of independent claims of the present invention are all in protection scope of the present invention.
Claims (7)
1., based on a method of estimation for the blood pressure envelope of oscillographic method, it is characterized in that: comprise the following steps:
S1, after measurement terminates, obtains the static pressure of dynamic pressure peak value and correspondence by computing, two threshold values of the scope of setting dynamic pressure amplitude are y
t1and y
t2, choose dynamic pressure sequences y
i(y
t1≤ y
i≤ y
t2), the static pressure sequence that dynamic pressure sequence pair is answered is x
i, corresponding functional relationship is y
i(x
i);
S2, if the function of matched curve is f (x
i, β), parameter vector β=(a
0, a
1..., a
i..., a
n, b, c), wherein a
0, a
1..., a
i..., a
n, b, c be parameter of curve variable, N is natural number, 0≤i≤N; With dynamic pressure sequences y selected in S1
ibased on, calculate iteration initial value a
0, a
1..., a
i..., a
n, b, c; Wherein, described a
0=y
m, y
mfor dynamic pressure sequences y
iin maximum; a
i=0,0 < i≤N; B=c
2; C=x
m, x
mfor y
mcorresponding static pressure; The formula of definition matched curve is
M be greater than 0 even number.；
S3, perform levenbergmarquardt iterative algorithm and find optimum parameter vector, formula is
λ is each iteration stepping length, calculates β; Calculate
T1 is dynamic pressure sequences y
ithreshold value y
t1time i value, t2 is dynamic pressure sequences y
ithreshold value y
t2time i value;
S4, definition threshold value threshold and threshold_1, judge as | θ | < threshold or | E (β+θ)-E (β) | during < threshold_1, iteration terminates, proceed to next step, obtain optimum parameter vector β=(a simultaneously
0, a
1..., a
i..., a
n, b, c);
S5, preserves parameter of curve (a
0, a
1..., a
i..., a
n, b, c), extract envelope.
2. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 1, is characterized in that: described y
t1value be y
i40% of sequence amplitude.
3. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 1, is characterized in that: the N value described in step S2 is 1.
4. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 1, is characterized in that: the M value described in step S5 is 4.
5. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 1, is characterized in that: the value of described threshold is 5-25.
6. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 5, is characterized in that: the value of described threshold is 10.
7. the method for estimation of a kind of blood pressure envelope based on oscillographic method according to claim 1, is characterized in that: the value of described threshold_1 is initial value (a
0, a
1..., a
i..., a
n, b, c) * 0.005.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510732538.6A CN105342590B (en) | 2015-11-02 | 2015-11-02 | A kind of method of estimation of the blood pressure envelope based on oscillographic method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510732538.6A CN105342590B (en) | 2015-11-02 | 2015-11-02 | A kind of method of estimation of the blood pressure envelope based on oscillographic method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105342590A true CN105342590A (en) | 2016-02-24 |
| CN105342590B CN105342590B (en) | 2018-06-08 |
Family
ID=55318772
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510732538.6A Expired - Fee Related CN105342590B (en) | 2015-11-02 | 2015-11-02 | A kind of method of estimation of the blood pressure envelope based on oscillographic method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105342590B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP4008245A1 (en) * | 2020-12-02 | 2022-06-08 | Koninklijke Philips N.V. | Apparatus for determining an indicator representative for a physiological parameter |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1077042A1 (en) * | 1999-08-18 | 2001-02-21 | Critikon Company, L.L.C. | Calculation of a quality index in blood pressure measurement |
| US6475154B1 (en) * | 2000-11-08 | 2002-11-05 | Taidoc Corp., Ltd. | Method and apparatus for non-invasive blood-pressure measurement |
| JP2005323853A (en) * | 2004-05-14 | 2005-11-24 | Osachi:Kk | Hemodynamics evaluating apparatus |
| JP3953073B2 (en) * | 2005-05-12 | 2007-08-01 | オムロンヘルスケア株式会社 | Oscillometric automatic blood pressure measuring device |
| CN101273886A (en) * | 2007-03-30 | 2008-10-01 | 通用电气公司 | Method of controlling inflation of a cuff in blood pressure determination |
| CN102018507A (en) * | 2010-11-25 | 2011-04-20 | 北京悦琦创通科技有限公司 | Data analysis processing method for measuring blood pressure by oscillography |
-
2015
- 2015-11-02 CN CN201510732538.6A patent/CN105342590B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1077042A1 (en) * | 1999-08-18 | 2001-02-21 | Critikon Company, L.L.C. | Calculation of a quality index in blood pressure measurement |
| US6475154B1 (en) * | 2000-11-08 | 2002-11-05 | Taidoc Corp., Ltd. | Method and apparatus for non-invasive blood-pressure measurement |
| JP2005323853A (en) * | 2004-05-14 | 2005-11-24 | Osachi:Kk | Hemodynamics evaluating apparatus |
| JP3953073B2 (en) * | 2005-05-12 | 2007-08-01 | オムロンヘルスケア株式会社 | Oscillometric automatic blood pressure measuring device |
| CN101273886A (en) * | 2007-03-30 | 2008-10-01 | 通用电气公司 | Method of controlling inflation of a cuff in blood pressure determination |
| CN102018507A (en) * | 2010-11-25 | 2011-04-20 | 北京悦琦创通科技有限公司 | Data analysis processing method for measuring blood pressure by oscillography |
Non-Patent Citations (1)
| Title |
|---|
| 李雪情: "基于示波法和高斯拟合的血压测量方法", 《传感技术学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP4008245A1 (en) * | 2020-12-02 | 2022-06-08 | Koninklijke Philips N.V. | Apparatus for determining an indicator representative for a physiological parameter |
| WO2022117471A1 (en) | 2020-12-02 | 2022-06-09 | Koninklijke Philips N.V. | Apparatus for determining an indicator representative for a fluid responsiveness parameter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105342590B (en) | 2018-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107822608B (en) | Pulse wave feature extraction method based on Gaussian mixture model | |
| US9119536B2 (en) | Pressure gauge, blood pressure gauge, method of determining pressure values, method of calibrating a pressure gauge, and computer program | |
| EP2393423B1 (en) | Calculating cardiovascular parameters | |
| Gao et al. | Estimation of pulse transit time as a function of blood pressure using a nonlinear arterial tube-load model | |
| US6841389B2 (en) | Method of determining concentration of glucose in blood | |
| EP3220811B1 (en) | Blood pressure and arterial compliance estimation from arterial segments | |
| CN106028917B (en) | Method for determining the pulse wave velocity in artery | |
| KR100820159B1 (en) | Blood pressure measuring method and apparatus | |
| EP3571984B1 (en) | Systems for creating a model for use in detecting a peripheral arterial pressure decoupling in a subject and for detecting peripheral arterial pressure decoupling using arterial pressure waveform data and said model | |
| Hahn et al. | Subject-specific estimation of central aortic blood pressure using an individualized transfer function: a preliminary feasibility study | |
| US20140066788A1 (en) | Methods and apparatus for determining pulse transit time as a function of blood pressure | |
| CN106974631A (en) | A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal | |
| CN102018507A (en) | Data analysis processing method for measuring blood pressure by oscillography | |
| CN107992452B (en) | Method, device, storage medium and equipment for calculating central hemodynamics index | |
| CN108903929A (en) | The modified method, apparatus of heart rate detection, storage medium and system | |
| CN111493855A (en) | System and method for non-invasive measurement of individualized cardiac output | |
| CN105342590A (en) | Blood pressure enveloping estimation method based on oscillography | |
| CN106137184A (en) | Electrocardiosignal QRS complex detection method based on wavelet transformation | |
| EP2416701B1 (en) | Monitoring peripheral decoupling | |
| WO2015180554A1 (en) | Method of determining stiffness index of arterial network and system thereof | |
| US20180199833A1 (en) | Device for calculating systolic blood pressure using pulse transit time and method therefor | |
| CN108937896B (en) | Central arterial blood pressure measuring method and equipment | |
| Poliński et al. | Can we rely on the pulse transit time-pressure relationship-models comparison | |
| CN109745020A (en) | Non-intrusion type blood vessel elasticity appraisal procedure based on blood flow magnetoelectric effect | |
| CN106580276B (en) | Pulse wave conduction time obtaining method based on correlation |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180608 Termination date: 20181102 |