CN108261192B - Continuous blood pressure estimation method, device and equipment - Google Patents

Continuous blood pressure estimation method, device and equipment Download PDF

Info

Publication number
CN108261192B
CN108261192B CN201611270666.4A CN201611270666A CN108261192B CN 108261192 B CN108261192 B CN 108261192B CN 201611270666 A CN201611270666 A CN 201611270666A CN 108261192 B CN108261192 B CN 108261192B
Authority
CN
China
Prior art keywords
pulse
pulse wave
wave
blood pressure
order
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201611270666.4A
Other languages
Chinese (zh)
Other versions
CN108261192A (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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201611270666.4A priority Critical patent/CN108261192B/en
Publication of CN108261192A publication Critical patent/CN108261192A/en
Application granted granted Critical
Publication of CN108261192B publication Critical patent/CN108261192B/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/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • 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/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802Sensor mounted on worn items
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Cardiology (AREA)
  • Physiology (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 continuous blood pressure estimation method, a continuous blood pressure estimation device and continuous blood pressure estimation equipment. The method comprises the following steps: 1) and (3) a calibration link: sampling the pulse wave of the tested object to obtain a sample pulse characteristic parameter and acquiring a reference blood pressure value of the tested object, wherein the first pulse characteristic parameter comprises: the slope of the ascending branch of the second order pulse leading wave, the amplitude difference of the ascending branch of the second order pulse leading wave, the coverage area of the ascending branch of the second order pulse leading wave, the slope of the descending branch of the second order pulse leading wave, the amplitude difference of the descending branch of the second order pulse leading wave and the coverage area of the descending branch of the second order pulse leading wave; calculating to obtain a coefficient parameter of the blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter; measuring the pulse wave of the measured object in real time to obtain a real-time pulse characteristic parameter; and inputting the coefficient parameters and the real-time pulse characteristic parameters into a blood pressure estimation model so as to estimate the real-time pulse continuous blood pressure value of the measured object.

Description

Continuous blood pressure estimation method, device and equipment
Technical Field
The invention relates to the field of medical instruments, in particular to a continuous blood pressure estimation method, a continuous blood pressure estimation device and continuous blood pressure estimation equipment.
Background
The blood pressure is the pressure required to transport the blood of a person to various parts of the body. Blood pressure is the pressure acting on the wall of a blood vessel when blood flows in the blood vessel, and is the driving force for driving the blood to flow in the blood vessel. The ventricles contract and blood flows from the ventricles into the arteries, at which time the pressure of the blood against the arteries is highest, called the Systolic Blood Pressure (SBP). The ventricles relax, the arterial vessels retract elastically, and blood continues to flow forward slowly, but the blood pressure drops, when it drops to a minimum, the pressure at this point is called the Diastolic Blood Pressure (DBP). Blood pressure is an important index of human health state, and in various medical occasions, the blood pressure of a patient needs to be measured in real time. However, the conventional mercury sphygmomanometer has relatively high measurement accuracy, but is not suitable for continuously measuring the blood pressure in real time, and the conventional method for continuously measuring the blood pressure in real time has low measurement accuracy.
Disclosure of Invention
The technical problem to be solved by the embodiments of the present invention is to provide a continuous blood pressure estimation method, device and apparatus, which can effectively and simply realize real-time measurement of continuous blood pressure per stroke.
The invention provides a continuous blood pressure estimation method in a first aspect, which comprises the following steps:
the continuous blood pressure measuring equipment samples the pulse wave of a tested object to obtain a sample pulse characteristic parameter and obtains a reference blood pressure value of the tested object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the tested object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of the slope of the ascending branch of the second-order pulse wave, the amplitude difference of the ascending branch of the second-order pulse wave, the coverage area of the ascending branch of the second-order pulse wave, the slope of the descending branch of the second-order pulse wave, the amplitude difference of the descending branch of the second-order pulse wave and the coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by performing twice derivation on an original pulse wave signal;
calculating by the continuous blood pressure measuring equipment according to the reference blood pressure value and the sample pulse characteristic parameter to obtain a coefficient parameter of a blood pressure estimation model;
the continuous blood pressure measuring equipment measures the pulse wave of the measured object in real time so as to obtain real-time pulse characteristic parameters;
and the continuous blood pressure measuring equipment inputs the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
Wherein the content of the first and second substances,
slope of rising branch of the second order lead pulse wave
Figure GDA0002775627070000021
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1Is the initial wave trough of the second order pulse wave in the heartbeat cycleAmplitude of the dot, TPThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV1The time value corresponding to the initial wave trough point of the second-order pulse wave in the heartbeat period is obtained;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude of the initial valley point of the second order pulse wave in the heartbeat period;
the coverage area of the ascending branch of the second order pulse wave
Figure GDA0002775627070000022
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V1 is the initial wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV1The amplitude of the starting valley point V1 of the second order pulse leading wave in the heartbeat period;
slope of descending branch of the second order lead pulse wave
Figure GDA0002775627070000023
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude value T of the end point valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV2A time value corresponding to a trough point of a terminal point of the second-order pulse wave in the heartbeat period;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude of a terminal wave valley point of the second order pulse wave in the heartbeat period;
the coverage area of the descending branch of the second order pulse wave
Figure GDA0002775627070000031
Wherein P isThe main wave peak point of the second order pulse wave in the heartbeat period, V2 is the end wave valley point of the second order pulse wave in the heartbeat period, I is the free variable, I isiIs the amplitude of point I, IV2Is the amplitude of the end-point valley point V2 of the second-order pulse leading wave in the heartbeat cycle.
Wherein the blood pressure estimation model comprises a linear regression model, a non-linear regression model, or a machine learning model.
The linear regression model is y ═ Ax + B, wherein a and B are coefficient parameters of the linear regression model, x is a pulse characteristic parameter, and y is a blood pressure value.
Wherein, the pulse characteristic parameter still includes second pulse characteristic parameter, second pulse characteristic parameter includes: at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, descending isthmus relative height, counterpulsation wave relative height, lowest value of pulse wave data waveform, pulse cycle, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to descending isthmus, and rising velocity of pulse wave.
Wherein the linear regression model is c ═ k1a+k2b+k3Wherein k is1、k2And k3The parameter values are coefficient parameters of a linear regression model, a is the first pulse characteristic parameter, b is the second pulse characteristic parameter, and c is a blood pressure value.
The invention provides a continuous blood pressure measuring device, which comprises an acquisition module, a calculation module, an actual measurement module and an estimation module,
the acquisition module is used for sampling the pulse wave of the measured object to obtain a sample pulse characteristic parameter and acquiring a reference blood pressure value of the measured object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the measured object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of the slope of the ascending branch of the second-order pulse wave, the amplitude difference of the ascending branch of the second-order pulse wave, the coverage area of the ascending branch of the second-order pulse wave, the slope of the descending branch of the second-order pulse wave, the amplitude difference of the descending branch of the second-order pulse wave and the coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by performing twice derivation on an original pulse wave signal;
the calculation module is used for calculating to obtain a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
the actual measurement module is used for measuring the pulse wave of the measured object in real time so as to obtain real-time pulse characteristic parameters;
the estimation module is used for inputting the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
Wherein the content of the first and second substances,
slope of rising branch of the second order lead pulse wave
Figure GDA0002775627070000041
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude value T of the initial valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV1The time value corresponding to the initial wave trough point of the second-order pulse wave in the heartbeat period is obtained;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude of the initial valley point of the second order pulse wave in the heartbeat period;
the coverage area of the ascending branch of the second order pulse wave
Figure GDA0002775627070000042
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V1 is the initial wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV1The amplitude of the starting valley point V1 of the second order pulse leading wave in the heartbeat period;
slope of descending branch of the second order lead pulse wave
Figure GDA0002775627070000051
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude value T of the end point valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV2A time value corresponding to a trough point of a terminal point of the second-order pulse wave in the heartbeat period;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude of a terminal wave valley point of the second order pulse wave in the heartbeat period;
the coverage area of the descending branch of the second order pulse wave
Figure GDA0002775627070000052
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V2 is the terminal wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV2Is the amplitude of the end-point valley point V2 of the second-order pulse leading wave in the heartbeat cycle.
Wherein the blood pressure estimation model comprises a linear regression model, a non-linear regression model, or a machine learning model.
The linear regression model is y ═ Ax + B, wherein a and B are coefficient parameters of the linear regression model, x is a pulse characteristic parameter, and y is a blood pressure value.
Wherein, the pulse characteristic parameter still includes second pulse characteristic parameter, second pulse characteristic parameter includes: at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, descending isthmus relative height, counterpulsation wave relative height, lowest value of pulse wave data waveform, pulse cycle, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to descending isthmus, and rising velocity of pulse wave.
Wherein the linear regression model is c ═ k1a+k2b+k3Wherein k is1、k2And k3The parameter values are coefficient parameters of a linear regression model, a is the first pulse characteristic parameter, b is the second pulse characteristic parameter, and c is a blood pressure value.
A third aspect of the invention provides a continuous blood pressure measurement device comprising an interface circuit, a memory, and a processor, wherein the memory stores a set of program code therein, and the processor is configured to invoke the program code stored in the memory for performing the following operations:
sampling a pulse wave of a tested object to obtain a sample pulse characteristic parameter, and acquiring a reference blood pressure value of the tested object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the tested object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of a slope of an ascending branch of the second-order pulse wave, an amplitude difference of the ascending branch of the second-order pulse wave, a coverage area of the ascending branch of the second-order pulse wave, a slope of a descending branch of the second-order pulse wave, an amplitude difference of the descending branch of the second-order pulse wave and a coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by carrying out secondary derivation on the pulse wave;
calculating to obtain a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
measuring the pulse wave of the measured object in real time to obtain a real-time pulse characteristic parameter;
and inputting the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
In the above embodiment, the real-time pulse-per-beat continuous blood pressure value of the measured object is calculated by using at least one pulse characteristic parameter selected from the slope of the ascending branch of the second order lead pulse wave, the amplitude difference of the ascending branch of the second order lead pulse wave, the coverage area of the ascending branch of the second order lead pulse wave, the slope of the descending branch of the second order lead pulse wave, the amplitude difference of the descending branch of the second order lead pulse wave, and the coverage area of the descending branch of the second order lead pulse wave, the continuous blood pressure measurement can be realized without a cuff, the electrocardio measurement is also not required, the miniaturized and comfortable wearable design is convenient to realize, and the experimental data proves that the precision of calculating the blood pressure by adopting the pulse characteristic parameters is higher, so that the use requirements of various medical occasions can be met.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of one embodiment of a pulse wave of a heartbeat cycle in accordance with the present disclosure;
FIG. 2 is a schematic diagram of an embodiment of a second order derivative pulse wave of a heartbeat cycle according to the disclosure;
FIG. 3 is a flow chart of a continuous blood pressure estimation method disclosed in the practice of the present invention;
FIG. 4 is a flow chart of a first embodiment of a continuous blood pressure estimation method disclosed in the present application;
FIG. 5 is a flowchart of a second embodiment of a continuous blood pressure estimation method disclosed in the present application;
FIG. 6 is a flowchart of a third embodiment of a continuous blood pressure estimation method disclosed in the present application;
FIGS. 7a to d are graphs of experimental results obtained by a continuous blood pressure estimation method disclosed in the practice of the present invention;
FIG. 8 is a schematic structural diagram of a continuous blood pressure measuring device according to the present disclosure;
fig. 9 is a schematic structural diagram of a continuous blood pressure measuring device disclosed in the implementation of the invention.
Detailed Description
It will be understood that, although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first element may be termed a second element, and, similarly, a second element may be termed a first element, without departing from the scope of example embodiments. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used herein, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The pulse wave form of each heartbeat cycle is substantially similar in a plurality of heartbeat cycles, so the pulse wave form of one heartbeat cycle is described below. As shown in FIG. 1, when the x-axis is time and the y-axis is the amplitude value of the pulse wave, the pulse wave of one heart cycle is shown in FIG. 1. In fig. 1, point a is the main peak point of the pulse wave, point b is the start valley point of the pulse wave, and point c is the end valley point of the pulse wave. The ascending branch is called from point b to point a, and the descending branch is called from point a to point c. The pulse wave shown in fig. 1 is merely an example, and is not particularly limited, and the shape of the pulse wave is generally different from one human body to another, but the pulse wave has features of a main peak point, a start valley point, and an end valley point.
The second order derivation is performed on the pulse wave of one heart cycle to obtain a second order derived pulse wave of one heart cycle as shown in fig. 2. In fig. 2, point a is a main peak point of the second order lead pulse wave, point B is a start valley point of the second order lead pulse wave, and point C is an end valley point of the second order lead pulse wave. The ascending branch is called from the point B to the point A, and the descending branch is called from the point A to the point C.
The pulse characteristic parameters are parameters capable of reflecting the characteristics of the pulse waves, and include both the characteristic parameters of the pulse waves and the characteristic parameters of the second-order pulse waves, for example, the amplitude of the main wave peak of the pulse waves is a pulse characteristic parameter of the pulse waves, which can reflect the maximum value of the pulse waves; the amplitude of the initial valley point of the pulse wave is a pulse characteristic parameter of the pulse wave, which can reflect the minimum value of the pulse wave at the ascending branch, the amplitude of the end valley point of the pulse wave is a pulse characteristic parameter of the pulse wave, which can reflect the minimum value of the pulse wave at the descending branch, and so on. In the embodiment of the present invention, the pulse characteristic parameters are artificially divided into a first pulse characteristic parameter and a second pulse characteristic parameter. The first pulse characteristic parameter is a newly discovered and applied pulse characteristic parameter in the invention, and the second pulse characteristic parameter is a more common pulse characteristic parameter. The first pulse characteristic parameter is one or more combination of the slope of the ascending branch of the second order pulse leading wave, the amplitude difference of the ascending branch of the second order pulse leading wave, the coverage area of the ascending branch of the second order pulse leading wave, the slope of the descending branch of the second order pulse leading wave, the amplitude difference of the descending branch of the second order pulse leading wave and the coverage area of the descending branch of the second order pulse leading wave. The second pulse characteristic parameter is one or more of a combination of pulse wave transmission time, pulse wave peak-to-trough amplitude ratio, main wave height, isthmus relative height, dicrotic wave relative height, minimum value of pulse wave data waveform, pulse period, duration of diastolic period, time ratio of systolic period to diastolic period, main wave rise time, area ratio of systolic period to diastolic period, area of pulse wave data waveform, slope of pulse wave peak to isthmus, and rise speed of pulse wave, etc. It is understood that, for reasons of space, only some of the second pulse characteristic parameters are listed here, and that more of the second pulse characteristic parameters are not listed here.
Since the first pulse characteristic parameters are found and applied first in the present invention, six first pulse characteristic parameters (the slope of the ascending branch of the second order pulse leading wave, the amplitude difference of the ascending branch of the second order pulse leading wave, the coverage area of the ascending branch of the second order pulse leading wave, the slope of the descending branch of the second order pulse leading wave, the amplitude difference of the descending branch of the second order pulse leading wave, and the coverage area of the descending branch of the second order pulse leading wave) will be described below with reference to fig. 2, while the second pulse characteristic parameters existing in the prior art will not be specifically described herein, and refer to the existing documents.
Slope of rising branch of the second order lead pulse wave
Figure GDA0002775627070000091
Wherein, IPIs the amplitude, I, of the main wave peak point A of the second order pulse wave in the heartbeat cycle of FIG. 2V1Is the amplitude, T, of the initial valley point B of the second order derivative pulse wave in the heartbeat cycle of FIG. 2PIs two in the heartbeat cycle of fig. 2The time value T corresponding to the main wave peak point A of the order leading pulse waveV1The time value corresponding to the start valley point B of the second-order pulse wave in the heartbeat cycle of fig. 2;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude, I, of the main wave peak point A of the second order pulse wave in the heartbeat cycle of FIG. 2V1The amplitude of the initial valley point B of the second-order pulse wave in the heartbeat cycle of FIG. 2;
the coverage area of the ascending branch of the second order pulse wave
Figure GDA0002775627070000092
Wherein, P is the main wave peak point a of the second order pulse leading wave in the heartbeat cycle of fig. 2, V1 is the initial wave valley point B of the second order pulse leading wave in the heartbeat cycle of fig. 2, I is the free variable, I is the peak point of the second order pulse leading wave in the heartbeat cycle of fig. 2iIs the amplitude of point I, IV1The amplitude of the start valley point V1 of the second order derivative pulse wave in the heartbeat cycle of fig. 2, and the end valley point is point C in fig. 2; it can be understood that the coverage area of the ascending branch of the second order derivative pulse wave is the "1" area shown in fig. 2;
slope of descending branch of the second order lead pulse wave
Figure GDA0002775627070000101
Wherein, IPIs the amplitude, I, of the main wave peak point A of the second order pulse wave in the heartbeat cycle of FIG. 2V2Is the amplitude, T, of the end-point valley point C of the second order lead pulse wave in the heartbeat cycle of FIG. 2PIs the time value, T, corresponding to the main wave peak point A of the second order pulse wave in the heartbeat cycle of FIG. 2V2The time value corresponding to the end point valley point C of the second-order pulse wave in the heartbeat cycle of fig. 2;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude, I, of the main wave peak point A of the second order pulse wave in the heartbeat cycle of FIG. 2V2The amplitude of the end-point valley point C of the second-order derived pulse wave in the heartbeat cycle of fig. 2;
the second mentionedCoverage area of descending branch of order guide pulse wave
Figure GDA0002775627070000102
Wherein, P is the main wave peak point a of the second order pulse leading wave in the heartbeat cycle of fig. 2, V2 is the end wave valley point C of the second order pulse leading wave in the heartbeat cycle, I is the free variable, I isiIs the amplitude of point I, IV2The amplitude of the end point valley point V2 of the second order pulse wave in the heartbeat cycle2Point C in fig. 2. It can be understood that the coverage area of the descending branch of the second order derivative pulse wave is the "2" area shown in fig. 2.
In order to continuously measure the blood pressure in real time, the invention provides a continuous blood pressure measuring device, which can execute a continuous blood pressure estimation method disclosed by the embodiment of the invention, and the specific steps of the continuous blood pressure measuring method are referred to below and are not described at first. The principle that the continuous blood pressure measuring device can continuously measure the blood pressure in real time is as follows: a blood pressure estimation model is established in the continuous blood pressure measuring equipment, and the blood pressure estimation model can reflect the relation between a first pulse characteristic parameter of a pulse wave of a measured object and a blood pressure value. The continuous blood pressure measuring equipment can continuously estimate the real-time pulse continuous blood pressure value of the measured object in real time by continuously measuring the first pulse characteristic parameter of the pulse wave in real time to obtain the real-time pulse characteristic parameter and inputting the real-time pulse characteristic parameter into the blood pressure estimation model. The blood pressure value may be any one of diastolic pressure and systolic pressure or a combination of both. The reason why the continuous blood pressure measuring device can continuously measure the first pulse characteristic parameter of the pulse wave in real time is that the continuous blood pressure measuring device measures the pulse wave through a photoelectric sensor placed at a measuring part outside the measured object, and the photoelectric sensor sends the measured pulse wave to a processor for calculation and analysis so as to obtain the pulse wave characteristic parameter. The measurement site may be any superficial artery capable of measuring pulse wave, such as fingertip artery, earlobe artery, and wrist artery. The pulse wave can be measured by the photoelectric sensor arranged outside the measured object, so the continuous blood pressure measuring device of the invention is non-invasive (i.e. does not cause wound to human body), and the photoelectric sensor can continuously collect the pulse wave of the measured object in real time, therefore, the continuous blood pressure measuring device can naturally and continuously estimate the blood pressure value of the measured object in real time. More specific process is described in the following text for a continuous blood pressure estimation method.
In the embodiment of the present invention, the six first pulse characteristic parameters may be used individually, or may be used partially or wholly in combination. Moreover, the effect of using the six first pulse feature parameters in combination, in part or in whole, by a linear regression model, a non-linear regression model, or a machine learning model is higher than the estimation accuracy of using any one of the six first pulse feature parameters alone. It is understood that, for purposes of space, only a portion of the models are illustrated, and that more regression prediction models are not illustrated.
As shown in fig. 3, the embodiment of the present invention provides a continuous blood pressure estimation method. The continuous blood pressure measuring method provided by the embodiment of the invention comprises the following steps:
310: the continuous blood pressure measuring equipment samples the pulse wave of the measured object to obtain a sample pulse characteristic parameter and obtains a reference blood pressure value of the measured object.
In an embodiment of the present invention, the pulse characteristic parameter is a characteristic parameter of a pulse wave of the measured object in a heartbeat cycle, and the pulse characteristic parameter includes a first pulse characteristic parameter, and the first pulse characteristic parameter includes: the slope of the ascending branch of the second order pulse leading wave, the amplitude difference of the ascending branch of the second order pulse leading wave, the coverage area of the ascending branch of the second order pulse leading wave, the slope of the descending branch of the second order pulse leading wave, the amplitude difference of the descending branch of the second order pulse leading wave and the coverage area of the descending branch of the second order pulse leading wave. In an optional embodiment, the pulse characteristic parameters further include a second pulse characteristic parameter, and the second pulse characteristic parameter includes at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, isthmus relative height, dicrotic wave relative height, lowest value of pulse wave data waveform, pulse period, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to isthmus, and rise speed of pulse wave.
In the embodiment of the invention, the sampling pulse characteristic parameters of the tested object are obtained by sampling through a continuous blood pressure measuring device. Wherein the continuous blood pressure measuring apparatus is a non-invasive continuous blood pressure measuring apparatus and is a sleeveless measuring apparatus. The continuous blood pressure measuring device is designed to be wearable so that a user can perform continuous blood pressure measurement anytime and anywhere, and the continuous blood pressure measuring device is not separated from a measured object due to movement or movement of the measured object.
In the embodiment of the invention, the reference blood pressure value of the measured object is obtained by measuring through a precision blood pressure measuring instrument and is sent to the continuous blood pressure measuring equipment in a wireless or wired mode. In another embodiment, the reference blood pressure value of the measured object may be manually input to the continuous blood pressure measuring apparatus by medical staff or the like. The precise blood pressure measuring instrument can be an invasive instrument or a non-invasive instrument. For example, when the precision blood pressure measuring instrument is an invasive instrument, the precision blood pressure measuring instrument may be a measuring instrument that measures blood pressure by directly inserting a catheter of a pressure sensor into a blood vessel of a patient, or the like; when The precision blood pressure measuring instrument is a non-invasive instrument, The precision blood pressure measuring instrument may be a mercury sphygmomanometer or a mature commercial device certified by The Medical instrument certification authority such as The american society for The Advancement of Medical Instrumentation (AAMI), and The like. The precise blood pressure measuring instrument may be a cuff-type measuring instrument or a cuff-less measuring instrument.
It can be understood that the blood pressure value of the human body may change with time, and in order to ensure the accuracy of the coefficient parameter of the blood pressure estimation model calculated in step 320, it is required to ensure that the sampling pulse characteristic parameter of the measured object sampled by the continuous blood pressure measuring device and the sampling blood pressure of the measured object measured by the precise blood pressure measuring device are performed simultaneously. For example, at the same time, the blood pressure measuring device samples the sampled pulse characteristic parameters of the measured object at the finger tip of the index finger of the left hand of the measured object, and the precision blood pressure measuring device measures the reference blood pressure of the measured object at the left arm of the measured object. Or, at the same time, the blood pressure measuring device samples the sampling pulse characteristic parameters of the measured object at the left wrist of the measured object, and the precision blood pressure measuring device measures the reference blood pressure of the measured object at the left arm of the measured object.
320: and calculating to obtain a coefficient parameter of the blood pressure estimation model by the continuous blood pressure measuring equipment according to the reference blood pressure value and the sample pulse characteristic parameter.
In the embodiment of the present invention, the coefficient parameters of the blood pressure estimation models of different subjects are usually different for different subjects, so before formally using the continuous blood pressure measurement device to measure the real-time pulse-to-pulse continuous blood pressure value of the subject, the coefficient parameters of the blood pressure estimation models need to be calculated by using the reference blood pressure value obtained in step 310 and the sample pulse characteristic parameters. Wherein the model is derived from a linear regression model, a non-linear regression model, or a machine learning model. It is understood that, for purposes of space, only a portion of the models are illustrated, and that more regression prediction models are not illustrated.
330: and the continuous blood pressure measuring equipment measures the pulse wave of the measured object in real time so as to obtain real-time pulse characteristic parameters.
340: and the continuous blood pressure measuring equipment inputs the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
The continuous blood pressure measuring method of the present invention is described below by way of several specific examples, and it is to be understood that the following examples are for illustrative clarity only and are not to be construed as limiting.
The first embodiment is as follows: the pulse characteristic parameter of the measured object only takes one first pulse characteristic parameter (for example, the slope of the ascending branch of the second-order leading pulse wave), and at this time, the blood pressure estimation model is a linear regression model y ═ Ax + B, where a and B are both coefficient parameters of the linear regression model, x is the slope of the ascending branch, and y is the blood pressure value. As shown in figure 4 of the drawings,
410: the continuous blood pressure measuring equipment samples the slope of the ascending branch of the second-order pulse wave of the tested object to obtain the slope x of the ascending branch of the second-order pulse wave of two samples1,x2Simultaneously, two reference blood pressure values y of the measured object are obtained1,y2
420: the continuous blood pressure measuring equipment is used for measuring the blood pressure according to the reference blood pressure value y1,y2And the slope x of the rising branch of the sample second derivative pulse wave1,x2And calculating to obtain coefficient parameters A and B of the blood pressure estimation model. In particular, the reference blood pressure value y may be1,y2And the slope x of the rising branch of the sample second derivative pulse wave1,x2And respectively substituting the linear regression model into y ═ Ax + B to obtain an equation system (1):
Figure GDA0002775627070000141
the solution of the equation set (1) can be obtained,
Figure GDA0002775627070000142
therefore, the linear regression model can be expressed as
Figure GDA0002775627070000143
It can be understood that, in the above-mentioned two sampling calculations, the values of one coefficient parameter a and one coefficient parameter B are obtained, and in practical application, the values of the coefficient parameter a and the coefficient parameter B can be obtained through multiple sampling calculations by using a rule such as least square fitting, so as to improve the accuracy of the coefficient parameters a and B.
430: the continuous blood pressure measuring equipment measures the pulse wave of the measured object in real time so as to obtain the slope x of the ascending branch of the real-time second-order pulse wavei
440: the continuous blood pressure measuring equipment uses the coefficient parameter
Figure GDA0002775627070000144
And the slope x of the ascending branch of the real-time second order pulse waveiInputting the blood pressure estimation model to estimate real-time stroke continuous blood pressure value of the measured object
Figure GDA0002775627070000145
Example two: the pulse characteristic parameters of the object to be measured are two first pulse characteristic parameters (for example, the slope of the ascending branch of the second-order pulse wave and the slope of the descending branch of the second-order pulse wave), and at this time, the blood pressure estimation model is a linear regression model with y ═ A1xi1+A2xi2+ B, wherein, A1、A2And B are both coefficient parameters of the linear regression model, xi1Is the slope of the ascending branch of the second order leading pulse wave, xi2The slope of the descending branch of the second-order pulse leading wave, y is the blood pressure value, and i is the sampling time. As shown in figure 5 of the drawings,
510: the continuous blood pressure measuring equipment respectively samples the slope of the ascending branch of the second order pulse wave and the slope of the descending branch of the second order pulse wave of the tested object to obtain the slope of the ascending branch of the second order pulse wave and the slope (x) of the descending branch of the second order pulse wave of 3 groups of samples11,x12),(x21,x22) And (x)31,x32) And obtaining 3 reference blood pressure values (y) of the measured object1,y2,y3). Wherein x is11,x21,x31Is the sample value of the slope of the rising branch of the second derivative pulse wave of the sample, x12,x22,x32Is the sample value of the slope of the descending branch of the sample second order lead pulse wave.
520: continuous blood pressureThe measuring equipment is used for measuring the slope (x) of the ascending branch of the sample second-order lead pulse wave and the slope (x) of the descending branch of the sample second-order lead pulse wave according to 3 groups of the ascending branches of the sample second-order lead pulse wave11,x12),(x21,x22),(x31,x32) And 3 reference blood pressure values (y)1,y2,y3) Calculating to obtain a coefficient parameter A of the blood pressure estimation model1、A2And B. Specifically, the slopes of the ascending branches of the sample second-order lead pulse wave and the slopes of the descending branches of the sample second-order lead pulse wave (x) may be divided into 3 groups11,x12),(x21,x22),(x31,x32) And 3 of said reference blood pressure values (y)1,y2,y3) Respectively substituting the linear regression model into y ═ A1xi1+A2xi2+ B, thereby obtaining equation set (2):
Figure GDA0002775627070000151
solving equation set (2) to obtain A1、A2And the value of B.
It will be appreciated that the above calculation of a coefficient parameter a by three samples results in1A coefficient parameter A2And a value of a coefficient parameter B, and in practical application, the coefficient parameter A can be obtained by multiple sampling and calculation by adopting a least square fitting rule and other rules1、A2And the value of the coefficient parameter B to increase the coefficient parameter A1、A2And the accuracy of B.
530: the continuous blood pressure measuring equipment measures the pulse wave of the measured object in real time so as to obtain the slope x of the ascending branch of the real-time second-order pulse wavek1And the slope x of the descending branch of the real-time second-order pulse leading wavek2
540: the continuous blood pressure measuring equipment uses the coefficient parameter A1、A2And B and the slope x of the ascending branch of the real-time second order pulse wavei1And the slope x of the descending branch of the real-time second-order pulse leading wavei2Inputting the blood pressure estimation model to estimate the real-time pulse continuous blood pressure value y of the measured objectk=A1xk1+A2xk2+B。
Example three: the pulse characteristic parameters of the subject are 1 first pulse characteristic parameter (for example, the slope of the ascending branch of the second-order pulse wave) and 1 second pulse characteristic parameter (for example, the pulse wave transmission time), and the blood pressure estimation model is a linear regression model c ═ k at this time1a+k2b+k3Wherein k is1、k2And k3All the parameters are coefficient parameters of a linear regression model, wherein a is the slope of the ascending branch, b is the pulse wave transmission time, and c is the blood pressure value. As shown in figure 6 of the drawings,
610: the continuous blood pressure measuring equipment samples the slope of the rising branch of the second order pulse wave of the measured object and the pulse wave transmission time respectively to obtain the slope of the rising branch of 3 groups of samples and the sample pulse wave transmission time (a)1,b1),(a2,b2) And (a)3,b3) And obtaining 3 reference blood pressure values (c) of the measured object1,c2,c3). Wherein, a1,a2,a3Is the sample value of the slope of the rising branch of the second derivative pulse wave of the sample, b1,b2,b3Is the sampled value of the transmission time of the sample pulse wave.
620: the continuous blood pressure measuring equipment is based on the slope of the rising branch of the 3 groups of the sample second order derivative pulse waves and the sample pulse wave transmission time (a)1,b1),(a2,b2) And (a)3,b3) And 3 reference blood pressure values (c)1,c2,c3) Calculating to obtain a coefficient parameter k of the blood pressure estimation model1,k2And k3. Specifically, the slope of the rising branch of the second-order derivative pulse wave of 3 sets of the sample and the sample pulse wave transmission time (a) may be set1,b1),(a2,b2) And (a)3,b3) And 3 of said reference blood pressure values (c)1,c2,c3) Respectively substituting the linear regression model into c ═ k1a+k2b+k3Thereby obtaining equation set (3):
Figure GDA0002775627070000161
k can be obtained by solving equation set (3)1,k2,k3The value of (c).
It will be appreciated that the above calculation of a coefficient parameter k by three samples yields1A coefficient parameter k2And a coefficient parameter k3In practical application, the coefficient parameter k can be obtained by multiple sampling and calculation by adopting principles such as least square fitting and the like1、k2And k3To increase the coefficient parameter k1,k2,k3To the accuracy of (2).
630: the continuous blood pressure measuring equipment measures the pulse wave of the measured object in real time so as to obtain the slope a of the ascending branch of the real-time second-order pulse waveiAnd real time pulse wave transmission time bi
640: the continuous blood pressure measuring equipment measures the coefficient parameter k1,k2,k3And the slope a of the ascending branch of the real-time second order pulse leading waveiAnd real time pulse wave transmission time biInputting the blood pressure estimation model to estimate the real-time stroke continuous blood pressure value c of the measured objecti=k1ai+k2bi+k3
In order to prove that the accuracy of the real-time continuous blood pressure value measurement by using the first pulse wave characteristic parameter provided by the invention is higher than the accuracy of the real-time continuous blood pressure value measurement by using the second pulse wave characteristic parameter used in the prior art, an experimenter collects the continuous blood pressure data of 22 volunteers for experiment, thereby obtaining the experiment result shown in fig. 7.
As shown in fig. 7a, the area (a) in the figure shows the experimental results of the mean systolic blood pressure error values obtained by real-time continuous blood pressure value measurement respectively using the single first pulse wave feature parameter proposed by the present invention; the area (b) in the figure shows the experimental result of the mean value of systolic blood pressure errors obtained by respectively measuring the real-time continuous blood pressure values by using a single second pulse characteristic parameter used in the prior art; the area (c) in the figure shows the experimental result of the mean systolic blood pressure error value obtained by respectively carrying out real-time continuous blood pressure value measurement by adopting the single first pulse wave characteristic parameter provided by the invention and the pulse wave transmission time in the second pulse wave characteristic parameter used by the prior art. The larger the dark area in the figure, the larger the error and the worse the effect. As can be seen from the figure, the dark region of the region (b) is the most, the dark region of the region (a) is the next to the dark region, and the dark region of the region (c) is the least. The vertical axis represents the number of sampling samples for calibration, from small to large, the larger the number of sampling samples for calibration, the smaller the estimation error, and the better the effect.
As shown in fig. 7b, the area (a) in the figure represents the experimental result of the mean diastolic blood pressure error obtained by real-time continuous blood pressure value measurement respectively using the single first pulse wave feature parameter proposed by the present invention; the area (b) in the figure shows the experimental result of the mean value of diastolic pressure errors obtained by respectively carrying out real-time continuous blood pressure value measurement by using a single second pulse wave characteristic parameter used in the prior art; the area (c) in the figure shows the experimental result of the mean diastolic pressure error value obtained by real-time continuous blood pressure value measurement respectively by using the single first pulse wave characteristic parameter provided by the invention and the pulse wave transmission time in the second pulse wave characteristic parameter used in the prior art. The larger the dark area in the figure, the larger the error and the worse the effect. As can be seen from the figure, the dark region of the region (b) is the most, the dark region of the region (a) is the next to the dark region, and the dark region of the region (c) is the least. The vertical axis represents the number of sampling samples for calibration, from small to large, the larger the number of sampling samples for calibration, the smaller the estimation error, and the better the effect.
As shown in fig. 7c, the area (a) in the figure shows the experimental results of the standard deviation of systolic blood pressure error obtained by real-time continuous blood pressure value measurement respectively using the single first pulse wave feature parameter proposed by the present invention; the area (b) in the figure shows the experimental result of the standard deviation of systolic blood pressure error obtained by real-time continuous blood pressure value measurement by using a single second pulse characteristic parameter used in the prior art; the area (c) in the figure shows the experimental result of the systolic blood pressure error standard deviation obtained by respectively carrying out real-time continuous blood pressure value measurement by adopting the single first pulse wave characteristic parameter provided by the invention and the pulse wave transmission time in the second pulse wave characteristic parameter used by the prior art. The larger the dark area in the figure, the larger the error and the worse the effect. As can be seen from the figure, the dark region of the region (b) is the most, the dark region of the region (a) is the next to the dark region, and the dark region of the region (c) is the least. The vertical axis represents the number of sampling samples for calibration, from small to large, the larger the number of sampling samples for calibration, the smaller the estimation error, and the better the effect.
As shown in fig. 7d, the area (a) in the figure represents the experimental result of the diastolic pressure error standard deviation obtained by respectively performing real-time continuous blood pressure value measurement by using the single first pulse wave feature parameter proposed by the present invention; the area (b) in the figure shows the experimental result of the diastolic pressure error standard deviation obtained by respectively carrying out real-time continuous blood pressure value measurement by using a single second pulse characteristic parameter used in the prior art; the area (c) in the figure shows the experimental result of the diastolic pressure error standard deviation obtained by respectively carrying out real-time continuous blood pressure value measurement by adopting the single first pulse wave characteristic parameter provided by the invention and the pulse wave transmission time in the second pulse wave characteristic parameter used in the prior art. The larger the dark area in the figure, the larger the error and the worse the effect. As can be seen from the figure, the dark region of the region (b) is the most, the dark region of the region (a) is the next to the dark region, and the dark region of the region (c) is the least. The vertical axis represents the number of sampling samples for calibration, from small to large, the larger the number of sampling samples for calibration, the smaller the estimation error, and the better the effect.
While the method of the embodiments of the present invention has been described in detail, in order to better implement the above-described aspects of the embodiments of the present invention, the following also provides the related apparatus for implementing the aspects.
Referring to fig. 8, fig. 8 is a schematic structural diagram of a continuous blood pressure measuring device according to the present invention. The continuous blood pressure measuring apparatus of the present embodiment includes: an acquisition module 810, a calculation module 820, an actual measurement module 830, and an estimation module 840.
The obtaining module 810 is configured to sample a pulse wave of a measured object to obtain a sample pulse characteristic parameter, and obtain a reference blood pressure value of the measured object, where the pulse characteristic parameter is a characteristic parameter of the pulse wave of the measured object in a heartbeat cycle, the pulse characteristic parameter includes a first pulse characteristic parameter, and the first pulse characteristic parameter includes: at least one of a slope of an ascending branch of the second-order pulse wave, an amplitude difference of the ascending branch of the second-order pulse wave, a coverage area of the ascending branch of the second-order pulse wave, a slope of a descending branch of the second-order pulse wave, an amplitude difference of the descending branch of the second-order pulse wave and a coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by carrying out secondary derivation on the pulse wave;
the calculation module 820 is configured to calculate a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
the actual measurement module 830 is configured to measure the pulse wave of the measured object in real time to obtain a real-time pulse characteristic parameter;
the estimation module 840 is configured to input the coefficient parameter and the real-time pulse feature parameter into the blood pressure estimation model, so as to estimate a real-time pulse-to-pulse continuous blood pressure value of the measured object.
Alternatively,
slope of rising branch of the second order lead pulse wave
Figure GDA0002775627070000191
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude value T of the initial valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV1The time value corresponding to the initial wave trough point of the second-order pulse wave in the heartbeat period is obtained;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude of the initial valley point of the second order pulse wave in the heartbeat period;
coverage area of ascending branch of second order lead pulse wave
Figure GDA0002775627070000192
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V1 is the initial wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV1The amplitude of the starting valley point V1 of the second order pulse leading wave in the heartbeat period;
slope of descending branch of the second order lead pulse wave
Figure GDA0002775627070000193
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude value T of the end point valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV2A time value corresponding to a trough point of a terminal point of the second-order pulse wave in the heartbeat period;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude of a terminal wave valley point of the second order pulse wave in the heartbeat period;
the coverage area of the descending branch of the second order pulse wave
Figure GDA0002775627070000201
Wherein P is the second order pulse leading in the heartbeat periodThe main wave peak point of the wave, V2 the end point of the second order pulse wave in the heart cycle, I the free variable, IiIs the amplitude of point I, IV2Is the amplitude of the end-point valley point V2 of the second-order pulse leading wave in the heartbeat cycle.
Optionally, the blood pressure estimation model comprises a linear regression model, a non-linear regression model, or a machine learning model.
Optionally, the linear regression model is y ═ Ax + B, where a and B are both coefficient parameters of the linear regression model, x is a pulse characteristic parameter, and y is a blood pressure value.
Optionally, the pulse characteristic parameters further include a second pulse characteristic parameter, which includes: at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, descending isthmus relative height, counterpulsation wave relative height, lowest value of pulse wave data waveform, pulse cycle, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to descending isthmus, and rising velocity of pulse wave.
Optionally, the linear regression model is c ═ k1a+k2b+k3Wherein k is1、k2And k3The parameter values are coefficient parameters of a linear regression model, a is the first pulse characteristic parameter, b is the second pulse characteristic parameter, and c is a blood pressure value.
The continuous blood pressure measuring device according to the embodiment of the present invention can implement the continuous blood pressure measuring method shown in fig. 3 to 6, and please refer to fig. 3 to 6 and related embodiments specifically, which will not be repeated herein.
Referring to fig. 9, a continuous blood pressure measuring device provided by an embodiment of the present application includes: a storage unit 910, a communication interface 920, and a processor 930 coupled to the storage unit 910 and the communication interface 920. The storage unit 910 is configured to store instructions, the processor 920 is configured to execute the instructions, and the communication interface 920 is configured to communicate with other devices under the control of the processor 930. The processor 930 when executing the instructions may perform any one of the continuous blood pressure estimation methods in the above embodiments of the present application according to the instructions.
Processor 930 may also be referred to as a Central Processing Unit (CPU). The memory unit 910 may include a read-only memory and a random access memory, and provides instructions and data, etc. to the processor 930. A portion of the storage unit 910 may also include a non-volatile random access memory. The components of the wallpaper changing apparatus in a particular application are coupled together, for example, by a bus system. The bus system may include a power bus, a control bus, a status signal bus, and the like, in addition to the data bus. For clarity of illustration, however, the various buses are identified in the figure as the bus system 940. The method disclosed in the above embodiments of the present invention can be applied to the processor 930 or implemented by the processor 930. Processor 930 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 90. The processor 930 may be a general-purpose processor, a digital signal processor, an application specific integrated circuit, an off-the-shelf programmable gate array or other programmable logic device, discrete gate or transistor logic, or discrete hardware components, among others. Processor 930 may implement or perform the methods, steps, and logic blocks disclosed in embodiments of the present invention. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in the storage unit 910, and the processor 930 may read information in the storage unit 910, and implement the steps of the above method in combination with hardware thereof.
Sampling a pulse wave of a tested object to obtain a sample pulse characteristic parameter, and acquiring a reference blood pressure value of the tested object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the tested object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of a slope of an ascending branch of the second-order pulse wave, an amplitude difference of the ascending branch of the second-order pulse wave, a coverage area of the ascending branch of the second-order pulse wave, a slope of a descending branch of the second-order pulse wave, an amplitude difference of the descending branch of the second-order pulse wave and a coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by carrying out secondary derivation on the pulse wave;
calculating to obtain a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
measuring the pulse wave of the measured object in real time to obtain a real-time pulse characteristic parameter;
and inputting the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
Alternatively,
slope of rising branch of the second order lead pulse wave
Figure GDA0002775627070000221
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude value T of the initial valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV1The time value corresponding to the initial wave trough point of the second-order pulse wave in the heartbeat period is obtained;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1Is the heart beatThe amplitude of the initial valley point of the second order lead pulse wave in the period;
the coverage area of the ascending branch of the second order pulse wave
Figure GDA0002775627070000222
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V1 is the initial wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV1The amplitude of the starting valley point V1 of the second order pulse leading wave in the heartbeat period;
slope of descending branch of the second order lead pulse wave
Figure GDA0002775627070000223
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude value T of the end point valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV2A time value corresponding to a trough point of a terminal point of the second-order pulse wave in the heartbeat period;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude of a terminal wave valley point of the second order pulse wave in the heartbeat period;
the coverage area of the descending branch of the second order pulse wave
Figure GDA0002775627070000231
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V2 is the terminal wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV2Is the amplitude of the end-point valley point V2 of the second-order pulse leading wave in the heartbeat cycle.
Optionally, the blood pressure estimation model comprises a linear regression model, a non-linear regression model, or a machine learning model.
Optionally, the linear regression model is y ═ Ax + B, where a and B are both coefficient parameters of the linear regression model, x is a pulse characteristic parameter, and y is a blood pressure value.
Optionally, the pulse characteristic parameters further include a second pulse characteristic parameter, which includes: at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, descending isthmus relative height, counterpulsation wave relative height, lowest value of pulse wave data waveform, pulse cycle, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to descending isthmus, and rising velocity of pulse wave.
Optionally, the linear regression model is c ═ k1a+k2b+k3Wherein k is1、k2And k3The parameter values are coefficient parameters of a linear regression model, a is the first pulse characteristic parameter, b is the second pulse characteristic parameter, and c is a blood pressure value.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (7)

1. A continuous blood pressure measuring device is characterized by comprising an acquisition module, a calculation module, an actual measurement module and an estimation module,
the acquisition module is used for sampling the pulse wave of the measured object to obtain a sample pulse characteristic parameter and acquiring a reference blood pressure value of the measured object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the measured object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of the slope of the ascending branch of the second-order pulse wave, the amplitude difference of the ascending branch of the second-order pulse wave, the coverage area of the ascending branch of the second-order pulse wave, the slope of the descending branch of the second-order pulse wave, the amplitude difference of the descending branch of the second-order pulse wave and the coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by performing twice derivation on an original pulse wave signal;
the calculation module is used for calculating to obtain a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
the actual measurement module is used for measuring the pulse wave of the measured object in real time so as to obtain real-time pulse characteristic parameters;
the estimation module is used for inputting the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
2. The apparatus of claim 1,
slope of rising branch of the second order lead pulse wave
Figure FDA0002775627060000011
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude value T of the initial valley point of the second order pulse wave in the heartbeat cyclePIs the time value corresponding to the main wave crest point of the second order pulse wave in the heartbeat period,TV1the time value corresponding to the initial wave trough point of the second-order pulse wave in the heartbeat period is obtained;
the amplitude difference PPG _ AID ═ I of the ascending branch of the second order pulse waveP-IV1I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV1The amplitude of the initial valley point of the second order pulse wave in the heartbeat period;
the coverage area of the ascending branch of the second order pulse wave
Figure FDA0002775627060000021
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, V1 is the initial wave valley point of the second order pulse wave in the heartbeat cycle, I is the free variable, IiIs the amplitude of point I, IV1The amplitude of the starting valley point V1 of the second order pulse leading wave in the heartbeat period;
slope of descending branch of the second order lead pulse wave
Figure FDA0002775627060000022
Wherein, IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude value T of the end point valley point of the second order pulse wave in the heartbeat cyclePThe time value T corresponding to the main wave crest point of the second order pulse wave in the heartbeat cycleV2A time value corresponding to a trough point of a terminal point of the second-order pulse wave in the heartbeat period;
the amplitude difference PPG _ DID ═ I of the descending branch of the second order pulse leading waveP-IV2I, wherein IPIs the amplitude of the main wave peak point of the second order pulse wave in the heartbeat cycle, IV2The amplitude of a terminal wave valley point of the second order pulse wave in the heartbeat period;
the coverage area of the descending branch of the second order pulse wave
Figure FDA0002775627060000023
Wherein, P is the main wave peak point of the second order pulse wave in the heartbeat cycle, and V2 is the main wave peak point in the heartbeat cycleThe second order lead pulse wave has the end point valley point, I is free variable, IiIs the amplitude of point I, IV2Is the amplitude of the end-point valley point V2 of the second-order pulse leading wave in the heartbeat cycle.
3. The apparatus of claim 1, wherein the blood pressure estimation model comprises a linear regression model, a non-linear regression model, or a machine learning model.
4. The apparatus according to claim 3, wherein the linear regression model is y-Ax + B, where A and B are both coefficient parameters of the linear regression model, x is a pulse feature parameter, and y is a blood pressure value.
5. The apparatus of claim 3, wherein the pulse feature parameters further comprise a second pulse feature parameter, the second pulse feature parameter comprising: at least one of pulse wave transmission time, pulse wave peak to trough amplitude ratio, dominant wave height, descending isthmus relative height, counterpulsation wave relative height, lowest value of pulse wave data waveform, pulse cycle, duration of diastole, time ratio of systole to diastole, dominant wave rise time, area ratio of systole to diastole, area of pulse wave data waveform, slope of pulse wave peak to descending isthmus, and rising velocity of pulse wave.
6. The apparatus of claim 5, wherein the linear regression model is c-k1a+k2b+k3Wherein k is1、k2And k3The parameter values are coefficient parameters of a linear regression model, a is the first pulse characteristic parameter, b is the second pulse characteristic parameter, and c is a blood pressure value.
7. A continuous blood pressure measurement device comprising an interface circuit, a memory, and a processor, wherein the memory stores a set of program code therein, and wherein the processor is configured to invoke the program code stored in the memory for performing the following operations:
sampling a pulse wave of a tested object to obtain a sample pulse characteristic parameter, and acquiring a reference blood pressure value of the tested object, wherein the pulse characteristic parameter is a characteristic parameter of the pulse wave of the tested object in a heartbeat period, the pulse characteristic parameter comprises a first pulse characteristic parameter, and the first pulse characteristic parameter comprises: at least one of the slope of the ascending branch of the second-order pulse wave, the amplitude difference of the ascending branch of the second-order pulse wave, the coverage area of the ascending branch of the second-order pulse wave, the slope of the descending branch of the second-order pulse wave, the amplitude difference of the descending branch of the second-order pulse wave and the coverage area of the descending branch of the second-order pulse wave, wherein the reference blood pressure value of the measured object is measured by a precision blood pressure measuring instrument, and the second-order pulse wave is obtained by performing twice derivation on an original pulse wave signal;
calculating to obtain a coefficient parameter of a blood pressure estimation model according to the reference blood pressure value and the sample pulse characteristic parameter;
measuring the pulse wave of the measured object in real time to obtain a real-time pulse characteristic parameter;
and inputting the coefficient parameters and the real-time pulse characteristic parameters into the blood pressure estimation model so as to estimate the real-time pulse-to-pulse continuous blood pressure value of the measured object.
CN201611270666.4A 2016-12-30 2016-12-30 Continuous blood pressure estimation method, device and equipment Active CN108261192B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611270666.4A CN108261192B (en) 2016-12-30 2016-12-30 Continuous blood pressure estimation method, device and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611270666.4A CN108261192B (en) 2016-12-30 2016-12-30 Continuous blood pressure estimation method, device and equipment

Publications (2)

Publication Number Publication Date
CN108261192A CN108261192A (en) 2018-07-10
CN108261192B true CN108261192B (en) 2021-01-01

Family

ID=62771270

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611270666.4A Active CN108261192B (en) 2016-12-30 2016-12-30 Continuous blood pressure estimation method, device and equipment

Country Status (1)

Country Link
CN (1) CN108261192B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109157204A (en) * 2018-08-07 2019-01-08 四川智琢科技有限责任公司 A kind of no cuff type wrist artery blood pressure measuring method and system
CN111214218B (en) * 2020-01-13 2024-02-09 京东方科技集团股份有限公司 Multi-physiological parameter detection equipment
CN113499045A (en) * 2021-06-30 2021-10-15 歌尔科技有限公司 Blood pressure monitoring method, system, device, master and slave wearable devices and storage medium
CN117476241B (en) * 2023-12-28 2024-04-19 柏意慧心(杭州)网络科技有限公司 Method, computing device and medium for determining a blood flow of a blood vessel

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008079688A (en) * 2006-09-26 2008-04-10 Citizen Holdings Co Ltd Blood pressure measuring instrument
CN102429649A (en) * 2011-12-14 2012-05-02 中国航天员科研训练中心 Continuous blood pressure measuring device
CN102488503A (en) * 2011-12-14 2012-06-13 中国航天员科研训练中心 Continuous blood pressure measurer
CN103976721A (en) * 2014-04-22 2014-08-13 辛勤 Blood pressure measuring method and embedded device for realizing method
CN104757955A (en) * 2015-03-25 2015-07-08 华中科技大学 Human body blood pressure prediction method based on pulse wave

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20160146393A (en) * 2015-06-12 2016-12-21 삼성전자주식회사 Apparatus for measuring blood pressure based on multiprocessing

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008079688A (en) * 2006-09-26 2008-04-10 Citizen Holdings Co Ltd Blood pressure measuring instrument
CN102429649A (en) * 2011-12-14 2012-05-02 中国航天员科研训练中心 Continuous blood pressure measuring device
CN102488503A (en) * 2011-12-14 2012-06-13 中国航天员科研训练中心 Continuous blood pressure measurer
CN103976721A (en) * 2014-04-22 2014-08-13 辛勤 Blood pressure measuring method and embedded device for realizing method
CN104757955A (en) * 2015-03-25 2015-07-08 华中科技大学 Human body blood pressure prediction method based on pulse wave

Also Published As

Publication number Publication date
CN108261192A (en) 2018-07-10

Similar Documents

Publication Publication Date Title
CN108261190B (en) Continuous blood pressure estimation method, device and equipment
Kurylyak et al. A Neural Network-based method for continuous blood pressure estimation from a PPG signal
CN108261192B (en) Continuous blood pressure estimation method, device and equipment
JP4571317B2 (en) Pressureless continuous blood pressure monitoring device
Ghosh et al. Continuous blood pressure prediction from pulse transit time using ECG and PPG signals
Sahoo et al. Wavelet based pulse rate and Blood pressure estimation system from ECG and PPG signals
WO2017066149A1 (en) Ambulatory blood pressure and vital sign monitoring apparatus, system and method
CN101990445B (en) Blood pressure estimating device
CN108261191B (en) Continuous blood pressure estimation method, device and equipment
CN109872820B (en) Method, device, equipment and storage medium for measuring blood pressure without cuff
CN112089405B (en) Pulse wave characteristic parameter measuring and displaying device
CN114652351B (en) Continuous blood pressure measuring method and device based on ultrasonic Doppler and electronic equipment
Bote et al. Evaluation of blood pressure estimation models based on pulse arrival time
Mishra et al. Cuffless blood pressure monitoring using PTT and PWV methods
Mamun et al. Using photoplethysmography & ECG towards a non-invasive cuff less blood pressure measurement technique
WO2020114190A1 (en) Recursion method and system for confidence degree of blood oxygen saturation measurement, and storage medium
CN109247929B (en) Blood pressure determination device, method, apparatus, and storage medium
Yen et al. Blood Pressure and Heart Rate Measurements Using Photoplethysmography with Modified LRCN.
Bhattacharjee Cuff-less blood pressure estimation from electrocardiogram and photoplethysmography based on VGG19-LSTM network
CN105726001B (en) A kind of blood pressure measuring method and device
JP2020110443A (en) Blood circulation information calculation apparatus, blood circulation information calculation method, and program
CN106175720A (en) The monitoring of a kind of Physiological And Biochemical Parameters and the method and device recorded
Laurenson et al. A 180 nm CMOS analog adaptive sampler for blood pressure feature extraction
Xing et al. Temporal complexity in photoplethysmography and its influence on blood pressure
TWI795219B (en) Method and system of detecting specific physiological syndrome related to blood circulation and deep sleep based on hemodynamic analysis

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