CN111685749A - Construction method of pulse pressure wave mathematical model - Google Patents

Construction method of pulse pressure wave mathematical model Download PDF

Info

Publication number
CN111685749A
CN111685749A CN202010557719.0A CN202010557719A CN111685749A CN 111685749 A CN111685749 A CN 111685749A CN 202010557719 A CN202010557719 A CN 202010557719A CN 111685749 A CN111685749 A CN 111685749A
Authority
CN
China
Prior art keywords
value
pressure wave
pulse pressure
pulse
pressure
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
Application number
CN202010557719.0A
Other languages
Chinese (zh)
Other versions
CN111685749B (en
Inventor
郑昕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN202010557719.0A priority Critical patent/CN111685749B/en
Publication of CN111685749A publication Critical patent/CN111685749A/en
Application granted granted Critical
Publication of CN111685749B publication Critical patent/CN111685749B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/02028Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02116Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave amplitude
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Vascular 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 discloses a construction method of a pulse pressure wave mathematical model, which uses a pressure sensor to collect a plurality of pulse pressure wave curves with maximum amplitude and synchronous/quasi-synchronous non-invasive blood pressure data; the pulse pressure wave curve data is refined and least square nonlinear fitting is carried out according to constraint conditions, so that a nonlinear system capable of fully and accurately reflecting the hemodynamic characteristics of a data source is obtained; and finally, constructing a pulse pressure wave mathematical model for extracting the hemodynamic parameters through the quantitative conversion of the waveform and the blood pressure data.

Description

Construction method of pulse pressure wave mathematical model
Technical Field
The invention belongs to the technology of noninvasive sign data measurement, in-vitro diagnosis and health monitoring, and particularly relates to a construction method of a pulse pressure wave mathematical model.
Background
The hemodynamic parameters can definitely reflect the individual physiological and pathological characteristics of a human body, can comprehensively and accurately provide reliable basis for cardiovascular system disease prevention and control, disease diagnosis, etiological analysis, trend speculation, diagnosis and treatment optimization and early high-risk symptom identification, and the acquisition modes of the hemodynamic parameters are invasive and noninvasive. In an invasive mode, a doctor needs to carry out surgical catheterization to connect an invasive monitor with a cardiovascular system in a professional environment of a medical institution, and a series of hemodynamic parameters are directly obtained in real time; non-invasive methods are used to indirectly obtain different hemodynamic parameters using various specialized medical diagnostic devices such as CT, MRI, PDE, sphygmomanometer, and electrocardiograph.
A series of hemodynamic parameters are acquired by using a conventional sphygmomanometer and small in-vitro pulse wave acquisition equipment, so that the practical application limitations of the two methods can be effectively overcome, and the convenient application under wide scenes such as daily monitoring can be realized. The invention is one of the representatives of Chinese patent with publication number CN1121798A and name of 'a method and device for detecting and analyzing dynamic parameters of cardiovascular function'. However, the methods have the defects of asynchronous data of different sensors in a mixed mode, manual intervention for data acquisition, lack of precise control on basic data and the like, so that the accuracy and the stability of the measurement result of the hemodynamic parameters are not ideal.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a method for constructing a pulse pressure wave mathematical model.
In order to solve the technical problems, the invention provides the following technical scheme:
the invention provides a method for constructing a pulse pressure wave mathematical model, which comprises the following steps:
s10, collecting and screening pulse pressure wave curves, which comprises the following steps:
s11, using pressure sensor and clock, when pulse pressure wave amplitude H is maximum, HWMaxAt least 5 continuous and complete pulse pressure wave curves are collected and placed in a pulse pressure wave curve coordinate system taking the amplitude H of the pulse pressure wave as a vertical axis and the time T as a horizontal axis;
s12, extracting the waveform height value H of each pulse pressure wave curveiSum waveform width value Ti(ii) a Wherein, Ti=ti2-ti1,ti2And ti1Respectively time values of an end point and a starting point of the ith pulse pressure wave curve; respectively calculating the average value H-of the height of the waveform and the average value T-of the width of the waveform by using all the data excluding the maximum value and the minimum value;
s13, reserving and satisfying | H less than or equal to 0i-H-/H-is less than or equal to 10% and 0. ltoreq. Ti-all pulse pressure wave curves for the condition T-/T-less than or equal to 10%;
s20, constructing a pulse pressure wave nonlinear system, specifically comprising the following steps:
s21, marking D on each pulse pressure wave curve according to the known definition1~D77 known feature points in total; the inflection point F of the curve with the same position and the common characteristic on each pulse pressure wave curve1~FnLabeling as individual feature points; according to the set contour definition and constraint conditions, selecting a plurality of curve control points Z on both sides of all inflection points and between two inflection points1~ZnSo as to effectively inhibit the overproof deformation of the curve;
s22, using the known characteristic point D on each pulse pressure wave curve1、D2、D7Constructing triangles for the vertices and finding the geometric centers K of all triangles1~Kn(ii) a Translating all pulse pressure wave curves after the first one to the first one until K2~KnAnd K1Completely overlapping;
s23, on the overlapped curve group, using the smallest circle to surround the same name point of each curve to form D1~D7、F1~Fn、Z1~ZnA plurality of homonymous point discrete areas; finding the weighted median center D0 for each homonym discrete region1~D07、F01~F0n、Z01~Z0n
S24, according to the basic rule that the abscissa value is from small to large, each group is more than or equal to 3 points, and at least 1 coincident point is arranged on two adjacent sections of curves, carrying out least square nonlinear fitting on the points to obtain a continuous pulse pressure wave nonlinear system;
s30, measuring blood pressure to obtain pressure value F when pulse pressure wave appearsA MinPressure value F when pulse pressure wave disappearsB MinAnd a pressure value F at maximum amplitudeW MaxCorresponding DP diastolic, SP systolic and MAP mean arterial pressures;
s40, constructing a pulse pressure wave mathematical model, specifically comprising the following steps:
s41, converting the known characteristic point D0 of the pulse pressure wave nonlinear system1Coincident with the ordinate axis, the feature point D0 was replaced on the ordinate by the DP and SP values in S381And D02Amplitude value H of1And H2Carrying out quantitative conversion from the amplitude H to the blood pressure value BP on the ordinate;
s42, establishing the blood pressure value BP as the vertical axisThe pulse pressure nonlinear system coordinate system with the horizontal axis as the time T redefines and names the known characteristic point D001~D007And the individual character point F001~F00nAnd finishing the construction of the pulse pressure wave digital model.
As a preferred embodiment of the present invention, the step S30 of measuring blood pressure may be any one of the known methods of synchronous/quasi-synchronous blood pressure measurement, or the following non-invasive blood pressure measurement method synchronized with pulse pressure wave acquisition, wherein the method of acquiring pulse pressure wave with a pressure sensor and a clock and synchronously performing blood pressure measurement specifically includes the following steps:
s31, slowly and vertically pushing the pressure sensor to the radial artery, carrying out the 1 st measurement when the pulse pressure wave is observed to appear, and reading the pressure value F1And the pulse pressure wave amplitude value H1
S32, defining f as a step pressure value and pressurizing step by step for subsequent measurement; 2 nd and 3 rd measurements were made at the 1/2f step value;
s33, performing subsequent measurement according to the step value f after the 3 rd measurement, and when the pulse pressure wave amplitude value observed in the q-th measurement is smaller than that observed in the previous measurement, namely Hq<Hq-1After the completion of the measurement, the qth was performed by decompressing at 1/2f+1Secondary measurement;
s34, according to q-1、q、q+1Three groups of measurement data are used for calculating the maximum amplitude value H of the pulse pressure waveW MaxAnd corresponding pressure value FW Max(ii) a With a sampling rate of 200Hz or more and applying FW MaxAt least 5 continuous and complete pulse pressure wave curves are acquired and provided to the step S10 in real time locally or in the cloud; then continuing to carry out subsequent measurement according to the step value f;
s35, when the q is+nThe secondary measurement observes that the pulse pressure wave amplitude is less than 85% of the maximum amplitude, namely Hq+n≤HWAt 85%, performing subsequent measurements at 1/2f step until the pulse pressure wave disappears;
s36, using the three groups of measurement data of 1 st, 2 nd and 3 rd times to carry out zero point convergence on the longitudinal axis, and calculating the occurrence H of pulse pressure waveA MinPressure value F at → 0A Min(ii) a Using the last three groups of measurement data to carry out zero point convergence of the longitudinal axis and calculate the disappearance H of the pulse pressure waveB MinPressure value F at → 0B Min
S37, completing pulse wave pressure/blood pressure conversion calculation in real time at local or cloud:
P=F/e
p is a blood pressure value, F is a pressure value and is a pulse wave pressure/blood pressure conversion factor;
wherein: n + [ (T-T)/a + (D-D)/b ] } c
N is the effective area value when the sensor measures the pulse wave pressure, T is the measured height value T is the system reference height value, D is the measured wrist circumference value D is the system reference wrist circumference value, a, b, c are the system constant value;
by FA Min、FB Min、FW MaxRespectively calculating the diastolic pressure DP, the systolic pressure SP and the mean arterial pressure MAP corresponding to the above steps:
DP=FA Min/
SP=FB Min/
MAP=FW Max/
s38, since the diastolic blood pressure DP, the systolic blood pressure SP and the mean arterial pressure MAP are independent results calculated from independent measurements, they can be used with each other as known:
and carrying out cross check on formulas such as MAP (SP +2 xDP)/3 or MAP (DP + I/3) (DP-SP) and the like, and weighting adjustment by using a least square method to obtain the most probable values of the diastolic pressure DP, the systolic pressure SP and the average arterial pressure MAP.
Compared with the prior art, the invention has the following beneficial effects:
1. synchronously acquiring blood pressure and pulse pressure wave data by using the same pressure sensor in the same time dimension, strictly keeping the original correlation among the data and eliminating system errors;
2. the automatic measurement/discrimination mode thoroughly eliminates accidental errors caused by human intervention;
3. optimizing, recursion and adjustment are carried out on the large redundant observation data, and the high precision and high reliability of the pulse pressure wave digital model are ensured; therefore, the accuracy of calculating a series of hemodynamic parameters by directly reading or extracting basic data subsequently is ensured;
4. the noninvasive acquisition method of the hemodynamic parameters with the result precision equal to or very similar to that of the invasive method can be realized by carrying out model refinement on the benchmarks which are synchronous with invasive data;
5. a series of hemodynamic parameters are obtained conveniently, accurately and immediately, and independent real-time local and cloud wide-area application of a user is supported.
Drawings
FIG. 1 is a schematic diagram of a set of pulse pressure wave curves with maximum amplitude and continuity in a pulse pressure wave curve processing coordinate system according to the present invention.
FIG. 2 is a graphical representation of the pulse pressure wave curves retained after screening in accordance with the present invention.
FIG. 3 is a schematic diagram of labeling feature points and finding geometric centers according to the present invention.
FIG. 4 is a schematic diagram of the same-name discrete areas after coincidence of the pulse pressure wave curves of the present invention.
FIG. 5 is a schematic diagram of finding the center point of the weighted median from the homonym discrete region according to the present invention.
FIG. 6 is a non-linear plot of the pulse pressure waves generated in the present invention.
FIG. 7 is a schematic diagram of the coordinate system and the measurement process of blood pressure measurement and the mark points in the present invention.
Fig. 8 is a schematic diagram of a mathematical model in the hemodynamic coordinate system of the present invention.
Detailed Description
The preferred embodiments of the present invention will be described in conjunction with the accompanying drawings, and it will be understood that they are described herein for the purpose of illustration and explanation and not limitation.
In order to achieve the object of the present invention, in one embodiment of the present invention, a method for constructing a pulse pressure wave mathematical model is provided, which includes the following steps:
s10, collecting and screening pulse pressure wave curves, which comprises the following steps:
s11, using pressure sensor and clock, when pulse pressure wave amplitude H is maximum, HW MaxAt least 5 continuous and complete pulse pressure wave curves are acquired and placed in a pulse pressure wave curve coordinate system taking the pulse pressure wave amplitude H as a vertical axis and the time T as a horizontal axis as shown in figure 1;
s12, extracting the waveform height value H of each pulse pressure wave curveiSum waveform width value Ti(ii) a Wherein, Ti=ti2-ti1,ti2And ti1Respectively time values of an end point and a starting point of the ith pulse pressure wave curve; respectively calculating the average value H-of the height of the waveform and the average value T-of the width of the waveform by using all the data excluding the maximum value and the minimum value;
s13, reserving and satisfying | H less than or equal to 0i-H-/H-is less than or equal to 10% and 0. ltoreq. TiAll pulse pressure wave curves for the-T-/T-10% condition are shown in FIG. 2;
s20, constructing a pulse pressure wave nonlinear system, specifically comprising the following steps:
s21, marking D on each pulse pressure wave curve according to the known definition1~D77 known feature points in total; the inflection point F of the curve with the same position and the common characteristic on each pulse pressure wave curve1~FnLabeling as individual feature points; according to the set contour definition and constraint conditions, selecting a plurality of curve control points Z on both sides of all inflection points and between two inflection points1~ZnThe superscalar deformation with the effective suppression curve is shown in FIG. 3;
s22, using the known characteristic point D on each pulse pressure wave curve1、D2、D7Constructing triangles for the vertices and finding the geometric centers K of all triangles1~Kn(ii) a Translating all pulse pressure wave curves after the first one to the first one until K2~KnAnd K1Full coincidence is shown in fig. 4;
s23, on the overlapped curve group, using the smallest circle to surround the same name point of each curve to form D1~D7、F1~Fn、Z1~ZnA plurality of homonymous point discrete areas; finding the weighted median center D0 for each homonym discrete region1~D07、F01~F0n、Z01~Z0nAs shown in fig. 5;
s24, according to the basic rule that the abscissa value is from small to large, each group is more than or equal to 3 points, and at least 1 coincident point is arranged on two adjacent sections of curves, performing least square nonlinear fitting on the points to obtain a continuous pulse pressure wave nonlinear system as shown in figure 6;
s30, measuring blood pressure to obtain pressure value F when pulse pressure wave appearsA MinPressure value F when pulse pressure wave disappearsB MinAnd a pressure value F at maximum amplitudeW MaxThe corresponding DP diastolic, SP systolic and MAP mean arterial pressures are shown in figure 7;
s40, constructing a pulse pressure wave mathematical model, specifically comprising the following steps:
s41, converting the known characteristic point D0 of the pulse pressure wave nonlinear system1Coincident with the ordinate axis, the feature point D0 was replaced on the ordinate by the DP and SP values in S381And D02Amplitude value H of1And H2Carrying out quantitative conversion from the amplitude H to the blood pressure value BP on the ordinate;
s42, establishing a pulse pressure nonlinear system coordinate system with the blood pressure value BP as the vertical axis and the time T as the horizontal axis, redefining and naming the known characteristic points D001~D007And the individual character point F001~F00nAnd finishing the construction of the pulse pressure wave digital model.
On the pulse pressure wave digital model constructed by the invention, a series of hemodynamic parameters are calculated by directly reading or extracting basic data through the nonlinear relation between two basic elements of blood pressure and time or the linear relation of line segments between characteristic points, the pressure integral of an ejection period, a diastole and a stroke period and the like.
According to the actual situation, step S30 is to measure the blood pressure, and any one of the known methods can be selected to perform synchronous/quasi-synchronous blood pressure measurement, or the following non-invasive blood pressure measurement method synchronized with the pulse pressure wave acquisition is adopted, and the method performs the pulse pressure wave acquisition with the pressure sensor and the clock and synchronously performs the blood pressure measurement, and specifically includes the following steps:
s31, slowly and vertically pushing the pressure sensor to the radial artery, carrying out the 1 st measurement when the pulse pressure wave is observed to appear, and reading the pressure value F1And the pulse pressure wave amplitude value H1
S32, defining f as 5-10 mmHg after the step pressure value is converted, and pressurizing step by step for subsequent measurement; 2 nd and 3 rd measurements were made at the 1/2f step value;
s33, performing subsequent measurement according to the step value f after the 3 rd measurement, and when the pulse pressure wave amplitude value observed in the q-th measurement is smaller than that observed in the previous measurement, namely Hq<Hq-1After the completion of the measurement, the qth was performed by decompressing at 1/2f+1Secondary measurement;
s34, according to q-1、q、q+1Three groups of measurement data are used for calculating the maximum amplitude value H of the pulse pressure waveW MaxAnd corresponding pressure value FW Max(ii) a With a sampling rate of 200Hz or more and applying FW MaxAt least 5 continuous and complete pulse pressure wave curves are acquired and provided to the step S10 in real time locally or in the cloud; then continuing to carry out subsequent measurement according to the step value f;
s35, when the q is+nThe secondary measurement observes that the pulse pressure wave amplitude is less than 85% of the maximum amplitude, namely Hq+n≤HWAt 85%, performing subsequent measurements at 1/2f step until the pulse pressure wave disappears;
s36, using the three groups of measurement data of 1 st, 2 nd and 3 rd times to carry out zero point convergence on the longitudinal axis, and calculating the occurrence H of pulse pressure waveA MinPressure value F at → 0A Min(ii) a Using the last three groups of measurement data to carry out zero point convergence of the longitudinal axis and calculate the disappearance H of the pulse pressure waveB MinPressure value F at → 0B Min
S37, completing pulse wave pressure/blood pressure conversion calculation in real time at local or cloud:
P=F/
p is a blood pressure value, F is a pressure value and is a pulse wave pressure/blood pressure conversion factor;
wherein: n + [ (T-T)/a + (D-D)/b ] } c
N is the effective area value when the sensor measures the pulse wave pressure, T is the measured height value T is the system reference height value, D is the measured wrist circumference value D is the system reference wrist circumference value, a, b, c are the system constant value;
by FA Min、FB Min、FW MaxRespectively calculating the diastolic pressure DP, the systolic pressure SP and the mean arterial pressure MAP corresponding to the above steps:
DP=FA Min/
SP=FB Min/
MAP=FW Max/
s38, since the diastolic blood pressure DP, the systolic blood pressure SP and the mean arterial pressure MAP are independent results calculated from independent measurements, they can be used with each other as known:
and carrying out cross check on formulas such as MAP (SP +2 xDP)/3 or MAP (DP + 1/3) (DP-SP), and weighting adjustment by using a least square method to obtain the most probable values of diastolic pressure DP, systolic pressure SP and average arterial pressure MAP.
In order to further optimize the implementation effect of the present invention, in another embodiment of the present invention, in order to enhance and strengthen the implementation effect of the present invention, a known method, such as time domain feature resolution, is further applied and not limited, so as to achieve the purpose of obtaining the mapping point of the physiological features, which is not obvious in local features of the waveform and cannot be accurately located or even ignored.
In order to further optimize the implementation effect of the invention, in another embodiment of the invention, calibration and refinement of the noninvasive pulse pressure wave digital model in the big data background can be realized by continuously synchronizing the benchmarks with a certain number of invasive hemodynamic parameters.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that changes may be made in the embodiments and/or equivalents thereof without departing from the spirit and scope of the invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (2)

1. A method for constructing a pulse pressure wave mathematical model is characterized by comprising the following steps:
s10, collecting and screening pulse pressure wave curves, which comprises the following steps:
s11, using pressure sensor and clock, when pulse pressure wave amplitude H is maximum, HW MaxAt least 5 continuous and complete pulse pressure wave curves are collected and placed in a pulse pressure wave curve coordinate system taking the amplitude H of the pulse pressure wave as a vertical axis and the time T as a horizontal axis;
s12, extracting the waveform height value H of each pulse pressure wave curveiSum waveform width value Ti(ii) a Wherein, Ti=ti2-ti1,ti2And ti1Respectively time values of an end point and a starting point of the ith pulse pressure wave curve; respectively calculating the average value H-of the height of the waveform and the average value T-of the width of the waveform by using all the data excluding the maximum value and the minimum value;
s13, reserving and satisfying | H less than or equal to 0i-H-/H-is less than or equal to 10% and 0. ltoreq. Ti-all pulse pressure wave curves for the condition T-/T-less than or equal to 10%;
s20, constructing a pulse pressure wave nonlinear system, specifically comprising the following steps:
s21, marking D on each pulse pressure wave curve according to the known definition1~D77 known feature points in total; the inflection point F of the curve with the same position and the common characteristic on each pulse pressure wave curve1~FnLabeling as individual feature points; according to the set contour definition and constraint conditions, selecting a plurality of curve control points Z on both sides of all inflection points and between two inflection points1~ZnSo as to effectively inhibit the overproof deformation of the curve;
s22, using the known characteristic point D on each pulse pressure wave curve1、D2、D7Constructing triangles for the vertices and finding the geometric centers K of all triangles1~Kn(ii) a Translating all pulse pressure wave curves after the first one to the first one until K2~KnAnd K1Completely overlapping;
s23, on the overlapped curve group, using the smallest circle to surround the same name point of each curve to form D1~D7、F1~Fn、Z1~ZnA plurality of homonymous point discrete areas; finding the weighted median center D0 for each homonym discrete region1~D07、F01~F0n、Z01~Z0n
S24, according to the basic rule that the abscissa value is from small to large, each group is more than or equal to 3 points, and at least 1 coincident point is arranged on two adjacent sections of curves, carrying out least square nonlinear fitting on the points to obtain a continuous pulse pressure wave nonlinear system;
s30, measuring blood pressure to obtain pressure value F when pulse pressure wave appearsA MinPressure value F when pulse pressure wave disappearsB MinAnd a pressure value F at maximum amplitudeW MaxCorresponding DP diastolic, SP systolic and MAP mean arterial pressures;
s40, constructing a pulse pressure wave mathematical model, specifically comprising the following steps:
s41, converting the known characteristic point D0 of the pulse pressure wave nonlinear system1Coincident with the ordinate axis, the feature point D0 was replaced on the ordinate by the DP and SP values in S381And D02Amplitude value H of1And H2Carrying out quantitative conversion from the amplitude H to the blood pressure value BP on the ordinate;
s42, establishing a pulse pressure nonlinear system coordinate system with the blood pressure value BP as the vertical axis and the time T as the horizontal axis, redefining and naming the known characteristic points D001~D007And the individual character point F001~F00nAnd finishing the construction of the pulse pressure wave digital model.
2. The method for constructing a mathematical model of pulse pressure wave according to claim 1, wherein the step S30 is performed to measure blood pressure, optionally by any one of the known methods, or by the following non-invasive blood pressure measurement methods synchronized with the pulse pressure wave acquisition, wherein the method uses a pressure sensor and a clock to perform the pulse pressure wave acquisition and perform the blood pressure measurement synchronously, and comprises the following steps:
s31, slowly and vertically pushing the pressure sensor to the radial artery, carrying out the 1 st measurement when the pulse pressure wave is observed to appear, and reading the pressure value F1And the pulse pressure wave amplitude value H1
S32, defining f as a step pressure value and pressurizing step by step for subsequent measurement; 2 nd and 3 rd measurements were made at the 1/2f step value;
s33, performing subsequent measurement according to the step value f after the 3 rd measurement, and when the pulse pressure wave amplitude value observed in the q-th measurement is smaller than that observed in the previous measurement, namely Hq<Hq-1After the completion of the measurement, the qth was performed by decompressing at 1/2f+1Secondary measurement;
s34, according to q-1、q、q+1Three groups of measurement data are used for calculating the maximum amplitude value H of the pulse pressure waveW MaxAnd corresponding pressure value FW Max(ii) a With a sampling rate of 200Hz or more and applying FW MaxAt least 5 continuous and complete pulse pressure wave curves are acquired and provided to the step S10 in real time locally or in the cloud; then continuing to carry out subsequent measurement according to the step value f;
s35, when the q is+nThe secondary measurement observes that the pulse pressure wave amplitude is less than 85% of the maximum amplitude, namely Hq+n≤HWAt 85%, performing subsequent measurements at 1/2f step until the pulse pressure wave disappears;
s36, using the three groups of measurement data of 1 st, 2 nd and 3 rd times to carry out zero point convergence on the longitudinal axis, and calculating the occurrence H of pulse pressure waveA MinPressure value F at → 0A Min(ii) a Using the last three groups of measurement data to carry out zero point convergence of the longitudinal axis and calculate the disappearance H of the pulse pressure waveB MinPressure value F at → 0B Min
S37, completing pulse wave pressure/blood pressure conversion calculation in real time at local or cloud:
P=F/
p is a blood pressure value, F is a pressure value and is a pulse wave pressure/blood pressure conversion factor;
wherein: n + [ (T-T)/a + (D-D)/b ] } c
N is the effective area value when the sensor measures the pulse wave pressure, T is the measured height value T is the system reference height value, D is the measured wrist circumference value D is the system reference wrist circumference value, a, b, c are the system constant value;
by FA Min、FB Min、FW MaxRespectively calculating the diastolic pressure DP, the systolic pressure SP and the mean arterial pressure MAP corresponding to the above steps:
DP=FA Min/
SP=FB Min/
MAP=FW Max/
s38, since the diastolic blood pressure DP, the systolic blood pressure SP and the mean arterial pressure MAP are independent results calculated from independent measurements, they can be used with each other as known:
and carrying out cross check on formulas such as MAP (SP +2 xDP)/3 or MAP (DP + 1/3) (DP-SP), and weighting adjustment by using a least square method to obtain the most probable values of diastolic pressure DP, systolic pressure SP and average arterial pressure MAP.
CN202010557719.0A 2020-06-18 2020-06-18 Construction method of pulse pressure wave mathematical model Active CN111685749B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010557719.0A CN111685749B (en) 2020-06-18 2020-06-18 Construction method of pulse pressure wave mathematical model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010557719.0A CN111685749B (en) 2020-06-18 2020-06-18 Construction method of pulse pressure wave mathematical model

Publications (2)

Publication Number Publication Date
CN111685749A true CN111685749A (en) 2020-09-22
CN111685749B CN111685749B (en) 2022-09-02

Family

ID=72482019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010557719.0A Active CN111685749B (en) 2020-06-18 2020-06-18 Construction method of pulse pressure wave mathematical model

Country Status (1)

Country Link
CN (1) CN111685749B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1121798A (en) * 1994-08-16 1996-05-08 北京工业大学 Cardiovascular function dynamic parameter testing analysis method and apparatus
CN102018507A (en) * 2010-11-25 2011-04-20 北京悦琦创通科技有限公司 Data analysis processing method for measuring blood pressure by oscillography
CN102813507A (en) * 2011-06-08 2012-12-12 中国医学科学院阜外心血管病医院 Device and method for measuring blood pressure value of human body when human body is in stable state
CN104644151A (en) * 2015-02-01 2015-05-27 北京工业大学 Photoelectric volume pulse signal-based propagation prediction method for pulse pressure waveform
US20150320364A1 (en) * 2013-01-25 2015-11-12 Up-Med Gmbh Method of Approximating a Patient's Pulse Wave Based on Non-Invasive Blood Pressure Measurement, A Logic Unit Therefore and a System Therefore
CN108053708A (en) * 2017-12-22 2018-05-18 天津天堰科技股份有限公司 A kind of Dynamic Pulse emulation module
CN108882872A (en) * 2016-04-15 2018-11-23 欧姆龙株式会社 Biont information analytical equipment, Biont information analysis system, program and Biont information analysis method
CN109820498A (en) * 2019-01-28 2019-05-31 重庆邮电大学 A kind of Wrist blood pressure meter that binding identity information and its pulse wave approximating method
CN110461224A (en) * 2016-12-15 2019-11-15 薇心健康有限公司 Wearable pulse wave measuring system and method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1121798A (en) * 1994-08-16 1996-05-08 北京工业大学 Cardiovascular function dynamic parameter testing analysis method and apparatus
CN102018507A (en) * 2010-11-25 2011-04-20 北京悦琦创通科技有限公司 Data analysis processing method for measuring blood pressure by oscillography
CN102813507A (en) * 2011-06-08 2012-12-12 中国医学科学院阜外心血管病医院 Device and method for measuring blood pressure value of human body when human body is in stable state
US20150320364A1 (en) * 2013-01-25 2015-11-12 Up-Med Gmbh Method of Approximating a Patient's Pulse Wave Based on Non-Invasive Blood Pressure Measurement, A Logic Unit Therefore and a System Therefore
CN104644151A (en) * 2015-02-01 2015-05-27 北京工业大学 Photoelectric volume pulse signal-based propagation prediction method for pulse pressure waveform
CN108882872A (en) * 2016-04-15 2018-11-23 欧姆龙株式会社 Biont information analytical equipment, Biont information analysis system, program and Biont information analysis method
CN110461224A (en) * 2016-12-15 2019-11-15 薇心健康有限公司 Wearable pulse wave measuring system and method
CN108053708A (en) * 2017-12-22 2018-05-18 天津天堰科技股份有限公司 A kind of Dynamic Pulse emulation module
CN109820498A (en) * 2019-01-28 2019-05-31 重庆邮电大学 A kind of Wrist blood pressure meter that binding identity information and its pulse wave approximating method

Also Published As

Publication number Publication date
CN111685749B (en) 2022-09-02

Similar Documents

Publication Publication Date Title
Papaioannou et al. Non-invasive 24 hour ambulatory monitoring of aortic wave reflection and arterial stiffness by a novel oscillometric device: the first feasibility and reproducibility study
Liu et al. Patient-specific oscillometric blood pressure measurement
CN103892818B (en) A kind of non-invasive central arterial blood pressure measuring method and equipment
CN101732040B (en) Non-invasive multipath pulse wave detection device, system and analytical system
US20190175042A1 (en) Method and apparatus for determining blood pressure
CN103479343B (en) Central aortic pressure detection system and method based on oscillating sphygmomanometer signals
US10004473B2 (en) Heart rate detection method and device using heart sound acquired from auscultation positions
CN114145724A (en) Method for dynamically monitoring blood pressure based on ECG (electrocardiogram) and PPG (photoplethysmography) multiple physiological characteristic parameters
CN102657525B (en) Nondestructive detection system and method of cardiac function
CN111493855A (en) System and method for non-invasive measurement of individualized cardiac output
Wang et al. The non-invasive and continuous estimation of cardiac output using a photoplethysmogram and electrocardiogram during incremental exercise
CN113160921A (en) Construction method and application of digital human cardiovascular system based on hemodynamics
CN112674736A (en) Monitoring display method and system for automatically evaluating deformation of blood vessel
Choudhary et al. Analyzing seismocardiographic approach for heart rate variability measurement
CN108852307B (en) Non-contact non-invasive arteriosclerosis detection device
CN111685749B (en) Construction method of pulse pressure wave mathematical model
CN210095711U (en) Noninvasive continuous blood pressure measuring equipment
Dong et al. An integrated framework for evaluation on typical ECG-derived respiration waveform extraction and respiration
CN111839487B (en) Blood pressure measuring system without air bag
CN113670516B (en) Compression position positioning and pressure measuring method based on photoplethysmography imaging
CN112842355B (en) Electrocardiosignal heart beat detection and identification method based on deep learning target detection
CN105310678A (en) Detecting method for calculating SV (stroke volume) of heart on basis of pulse wave analysis method
Song et al. Non-invasive hemodynamic diagnosis based on non-linear pulse wave theory applied to four limbs
Hehenkamp et al. Comparison of Portapres® with standard sphygmomanometry in pregnancy
Sharman et al. Measurements of arterial pressure and flow in vivo

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant