[go: up one dir, main page]

WO2017086294A1 - Dispositif d'estimation de la pression artérielle, dispositif de traitement d'informations, procédé d'estimation de la pression artérielle et programme d'estimation de la pression artérielle - Google Patents

Dispositif d'estimation de la pression artérielle, dispositif de traitement d'informations, procédé d'estimation de la pression artérielle et programme d'estimation de la pression artérielle Download PDF

Info

Publication number
WO2017086294A1
WO2017086294A1 PCT/JP2016/083770 JP2016083770W WO2017086294A1 WO 2017086294 A1 WO2017086294 A1 WO 2017086294A1 JP 2016083770 W JP2016083770 W JP 2016083770W WO 2017086294 A1 WO2017086294 A1 WO 2017086294A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood pressure
pressure estimation
time
estimation device
parameter
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.)
Ceased
Application number
PCT/JP2016/083770
Other languages
English (en)
Japanese (ja)
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.)
Tohoku University NUC
Original Assignee
Tohoku University NUC
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 Tohoku University NUC filed Critical Tohoku University NUC
Priority to JP2017551876A priority Critical patent/JP6580158B2/ja
Publication of WO2017086294A1 publication Critical patent/WO2017086294A1/fr
Priority to US15/981,075 priority patent/US11284805B2/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02116—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave amplitude
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7278—Artificial waveform generation or derivation, e.g. synthesizing signals from measured signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024—Measuring pulse rate or heart rate
    • A61B5/02416—Measuring pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024—Measuring pulse rate or heart rate
    • A61B5/0245—Measuring pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements 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/6802—Sensor mounted on worn items
    • A61B5/681—Wristwatch-type devices

Definitions

  • the present invention relates to a blood pressure estimation device, an information processing device, a blood pressure estimation method, and a blood pressure estimation program.
  • a blood pressure estimation device that estimates the pressure (in other words, blood pressure) of blood flowing through a living body's circulatory system is known (see, for example, Patent Document 1).
  • the blood pressure estimation device detects the arterial diameter at two positions on the surface of the finger and also detects the electrical impedance of the living body between the two positions. Furthermore, the blood pressure estimation device estimates blood pressure based on the detected artery diameter and electrical impedance and the Kalman filter.
  • blood pressure may not be estimated with high accuracy.
  • One of the objects of the present invention is to estimate blood pressure with high accuracy.
  • the blood pressure estimation device includes: A detection unit for detecting a first parameter representing a cycle length that is a length of a heartbeat cycle of a living body; In a mathematical model representing blood flowing in the circulatory system of the living body by a fluid flowing through a flow path in which a plurality of elastically deforming containers are connected in an annular shape, the volume of at least one of the plurality of containers is with respect to time; A processor that determines a change based on the detected first parameter, and that estimates the pressure of the blood based on the determined change and the mathematical model; Is provided.
  • the information processing apparatus is Change in volume of at least one of the plurality of containers with respect to time in a mathematical model representing blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers that are elastically deformed are annularly connected
  • Determining unit based on the first parameter representing the detected cycle length, which is the length of the cycle of the heartbeat of the living body,
  • An estimation unit that estimates the pressure of the blood based on the determined change and the mathematical model; Is provided.
  • the blood pressure estimation method comprises: Change in volume of at least one of the plurality of containers with respect to time in a mathematical model representing blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers that are elastically deformed are annularly connected Is determined based on the first parameter representing the detected cycle length, which is the length of the cycle of the heartbeat of the living body, Based on the determined change and the mathematical model, the pressure of the blood is estimated.
  • the blood pressure estimation program includes: Change in volume of at least one of the plurality of containers with respect to time in a mathematical model representing blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers that are elastically deformed are annularly connected Is determined based on the first parameter representing the detected cycle length, which is the length of the cycle of the heartbeat of the living body, Based on the determined change and the mathematical model, the computer is caused to execute a process for estimating the pressure of the blood.
  • FIGS. 1 to 41 embodiments relating to a blood pressure estimation device, an information processing device, a blood pressure estimation method, and a blood pressure estimation program according to the present invention will be described with reference to FIGS. 1 to 41.
  • the blood pressure estimation device 1 includes a detection unit 10 and a processing unit 20.
  • the blood pressure estimation device 1 is a wristwatch type.
  • the blood pressure estimation device 1 may be a type different from the wristwatch type (for example, a bandage type).
  • the detection unit 10 detects a first parameter representing a cycle length that is a cycle length of a heartbeat (in other words, a pulse) of a living body.
  • the living body is a human living body.
  • the living body may be a living body of an animal other than a human being.
  • the first parameter is the pulse rate.
  • the pulse rate is the number of pulses of the living body per predetermined unit time (in this example, 1 minute).
  • the pulse rate is calculated by dividing the unit time by the period length of one pulse.
  • the period length may be calculated by dividing the unit time by the pulse rate.
  • the first parameter may be a cycle length.
  • the detection unit 10 detects a second parameter indicating the magnitude of the heartbeat of the living body.
  • the second parameter is a difference between the minimum value and the maximum value in one pulse of the signal representing the amount of blood in the artery.
  • the second parameter is a single pulse of the width of the artery, the cross-sectional area of the artery, the amount of blood in the artery, the flow rate of blood in the artery, the flow rate of blood in the artery, or a signal representing at least one of these. It may be the difference between the minimum value and the maximum value.
  • the second parameter may be expressed as a pulse amplitude.
  • the detection unit 10 irradiates light on the artery of the living body, detects the intensity of the light reflected by the living body, and based on the change of the detected intensity with respect to time, the pulse rate and Detect pulse amplitude.
  • the detector 10 detects the pulse rate based on the time between two consecutive peaks in the detected intensity change.
  • the detection unit 10 detects the pulse amplitude based on the difference between the minimum intensity value and the maximum intensity value between two successive peaks in the detected intensity change.
  • the detection unit 10 includes a member that presses the surface of the living body in the vicinity of the artery (for example, the wrist) and detects the pulse rate and the pulse amplitude by detecting the pressure that the member receives from the surface. May be.
  • the detection unit 10 may detect the pressure by using a piezoelectric element.
  • the detection unit 10 includes an electrode attached to the surface of the living body in the vicinity of the heart, and detects the first parameter and the second parameter by detecting the potential on the surface via the electrode. Good.
  • the detection unit 10 may be a bandage type.
  • the processing unit 20 includes a processing device 21, a storage device 22, an input device 23, and an output device 24 that are connected to each other via a bus BU.
  • the processing unit 20 is an example of an information processing device.
  • the processing device 21 controls each element constituting the processing unit 20 by executing a program stored in the storage device 22. Thereby, the processing unit 20 realizes functions to be described later.
  • the processing device 21 includes a CPU (Central Processing Unit).
  • the processing device 21 may include an MPU (Micro Processing Unit) or a DSP (Digital Signal Processor) instead of or in addition to the CPU. Further, the processing device 21 may be configured by an LSI (Large Scale Integration).
  • the storage device 22 stores information in a readable / writable manner.
  • the storage device 22 includes at least one of a RAM (Random Access Memory), a semiconductor memory, and an organic memory.
  • the storage device 22 may include at least one of an HDD (Hard Disk Drive) and an SSD (Solid State Drive).
  • the storage device 22 may include a recording medium such as a flexible disk, an optical disk, a magneto-optical disk, or a semiconductor memory, and a reading device that can read information from the recording medium.
  • the input device 23 inputs information from outside the blood pressure estimation device 1.
  • the input device 23 includes a key type button.
  • the input device 23 may include a microphone.
  • the output device 24 outputs information to the outside of the blood pressure estimation device 1.
  • the output device 24 includes a display.
  • the output device 24 may include a speaker.
  • the processing unit 20 may include a touch panel display that forms both the input device 23 and the output device 24.
  • the processing unit 20 estimates the pressure (in other words, blood pressure) of blood flowing through the circulatory system of the living body.
  • the mathematical model will be described.
  • the mathematical model represents blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers that are elastically deformed are connected in an annular shape.
  • the flow path in the mathematical model is formed by the first to eighth containers FV1 to FV8 and the first to eighth communication pipes FC1 to FC8.
  • the flow path may be formed by nine or more containers.
  • the flow path may be formed by nine or more containers and the same number of communication pipes as the containers.
  • the first to eighth containers FV1 to FV8 are connected in a ring shape.
  • the first to seventh containers FV1 to FV7 and the second to eighth containers FV2 to FV8 are connected by the second to eighth communication pipes FC2 to FC8, respectively.
  • the eighth container FV8 and the first container FV1 are connected by the first communication pipe FC1.
  • the fluid flows from the first to seventh containers FV1 to FV7 to the second to eighth containers FV2 to FV8, respectively. Further, in the flow path, the fluid flows from the eighth container FV8 to the first container FV1.
  • the first to eighth containers FV1 to FV8 include the left atrium, the left ventricle, the artery downstream of the aorta and the aorta, the vena cava and the vein upstream of the vena cava, the right atrium, the right ventricle, Each represents a pulmonary artery and a pulmonary vein.
  • Each of the containers FV1 to FV8 is a spherical shell that is elastically deformed.
  • the first to eighth differential dP i / dt with respect to time t of the pressure P i of the fluid in the container FVi of the i is represented by Equation 1.
  • i represents each integer of 1 to 8.
  • Q i represents the flow rate flowing into the i-th container FVi.
  • Q i + 1 represents the flow rate flowing out from the i-th container FVi.
  • V i is (in other words, no-load volume) pressure P i of the fluid in the container FVi i th volume of the container FVi of the i in the case of 0 represents a.
  • E i represents a predetermined coefficient associated with the i th container FVi.
  • the coefficient E i may be regarded as a parameter representing a ratio of a change in the pressure P i of the fluid in the i-th container FVi with respect to a change in the volume of the i-th container FVi with respect to time.
  • the coefficient E i represents the amount of fluid flowing out from the i-th container FVi per unit time and the change of the no-load volume Vi of the i- th container FVi with respect to time, in the i-th container FVi. It may be regarded as a parameter that represents a ratio of a change in the pressure P i of the fluid in the i-th container FVi with respect to an amount subtracted from the amount of fluid flowing in.
  • Unloaded volumes V 1 , V 3 to V 5 , VV of the containers FV1, FV3 to FV5, FV7, and FV8 other than the second and sixth containers FV2 and FV6 among the first to eighth containers FV1 to FV8. 7, V 8 does not change with respect to time.
  • the unloaded volumes V 2 and V 6 of the second and sixth containers FV2 and FV6 change with time as expressed in Equation 2. Changes over time of the unloaded volumes V 2 , V 6 of the second and sixth containers FV2, FV6 may be taken to represent pulsations of the left and right ventricles.
  • A represents the pulse amplitude detected by the detector 10.
  • b represents the pulse rate detected by the detection unit 10.
  • ⁇ represents the time (in other words, the time within the cycle) from the time when the cycle starts in one pulse cycle.
  • no-load volume f i represents the no-load volume V i .
  • No-load volume f i has, pulse amplitudes a, pulse rate b, and the predetermined value in accordance with the time tau.
  • no-load volume f i is expressed by Equation 3.
  • the unloaded volumes f 2 and f 6 of the second and sixth containers FV2 and FV6 are determined based on the pulse amplitude a and the pulse rate b.
  • f i, 0 represents (in other words, reference unloaded volume) pulse amplitude a reference value of the no-load volume f i when the reference value a 0 a.
  • the reference unloaded volume f i, 0 is a ratio (in other words, a pulse rate ratio) of time ⁇ from a time point when the cycle starts to a reference value b 0 of the pulse rate b in a single pulse cycle (in other words, a pulse rate ratio) b /. b 0 and a predetermined value.
  • the reference unloaded volumes f 2,0 , f 6,0 for the second and sixth containers FV2, FV6 are represented by the curves VL, VR of FIG.
  • the flow rate Q 9 exiting from the container FV8 eighth is equal to the flow rate Q 1 that flows into the first container FV1. Further, the flow rates Q i flowing into the first to eighth containers FV1 to FV8 may be regarded as the flow rates in the first to eighth communication pipes FC1 to FC8, respectively.
  • Equation 5 The flow rate Q i in each of the communication pipes FC1 to FC8 is expressed by Equation 5.
  • the function C is obtained when the pressure P i ⁇ 1 of the upstream container connected to the i th communication pipe FCi is equal to the pressure of the downstream container connected to the i th communication pipe FCi. If less than the pressure P i, it represents 0. Further, as expressed by Expression 6, the function C is obtained when the pressure P i-1 of the upstream container connected to the i- th communication pipe FCi is reduced to the downstream side connected to the i-th communication pipe FCi. When the pressure is not less than the pressure P i of the container, 1 is represented.
  • the function C may be taken to represent a valve that prevents the fluid from flowing back (in other words, a check valve).
  • the pressure P 0 of the upstream container connected to the first communication pipe FC1 is equal to the pressure P 8 of the downstream container connected to the eighth communication pipe FC8.
  • R i represents the resistance to fluid flow from the upstream container connected to the i th communication pipe FCi to the downstream container connected to the i th communication pipe FCi.
  • the resistance R i is regarded as a ratio of the difference between the pressures of the fluid in the pair of containers to the flow rate of the fluid between the pair of containers connected to each other among the plurality of containers FV1 to FV8. It's okay.
  • g i represents a target value of the resistance R i (in other words, target resistance).
  • the resistance R i may be grasped as approaching the target resistance g i with a delay.
  • Target resistance g i have, pulse amplitude a, and the predetermined value in accordance with the pulse rate b.
  • the target resistance g i is expressed by Equation 9.
  • g i, 0 represents (in other words, the reference target resistance) pulse amplitude a reference value of the target resistance g i when the reference value a 0 a.
  • Reference target resistance g i, 0 has a value predetermined according to the pulse rate ratio b / b 0.
  • the reference target resistance g i, 0, as represented in FIG. 5, decreases as the pulse rate ratio b / b 0 increases.
  • gamma i represents (in other words, the reference target resistance ratio) the target resistance g i ratio reference target resistance g i, 0 of the.
  • the reference target resistance ratio ⁇ i has a predetermined value according to the ratio of the pulse amplitude a to the reference value a 0 (in other words, the pulse amplitude ratio) a / a 0 .
  • the reference target resistance ratio gamma i this example, as represented in FIG. 6, with decreases as the pulse amplitude ratio a / a 0 becomes large and the pulse amplitude ratio a / a 0 is 1 1 It is.
  • T i represents the time constant of the change of the resistance R i .
  • the time constant T i is the change in resistance R i may be seen as indicative of the degree of delay with respect to the change in target resistance g i.
  • resistors R 4, R 8 for the fourth and eighth communicating pipe FC4, FC8 of the pulse amplitude a and is determined based on the pulse rate b.
  • the function of the processing unit 20 includes a storage unit 201, a determination unit 202, and an estimation unit 203.
  • the processing unit 20 calculates the pressures P 1 (t) to P 8 (t) of the fluid in the first to eighth containers FV1 to FV8 at time (in other words, time point t) based on Expressions 1 to 9. From the time t, a process for estimating the pressures P 1 (t + ⁇ t) to P 8 (t + ⁇ t) of the fluid in the first to eighth containers FV1 to FV8 at a time t + ⁇ t after a predetermined step time ⁇ t has elapsed. Run repeatedly.
  • the storage unit 201 stores a reference no-load volume f i, 0 and a reference target resistance g i, 0 in advance.
  • the processing unit 20, instead of storing the reference no-load volume f i, 0, reference no-load volume f i, 0 may store function calculated.
  • the processing unit 20, instead of storing the reference target resistance g i, 0, may store a function for calculating a reference target resistance g i, 0.
  • the storage unit 201 stores a reference target resistance ratio ⁇ i in advance.
  • the processing unit 20, instead of storing the reference target resistance ratio gamma i may store a function for calculating a reference target resistance ratio gamma i.
  • the estimation unit 203 sets the time t, the pressures P 1 to P 8 , and the resistances R 1 to R 8 to an initial value t ini , an initial value P 1, ini to P 8, in i, and an initial value R 1, ini. To R8 , ini , respectively.
  • the estimation unit 203 acquires the pulse rate b (t) and the pulse amplitude a (t) at time t based on the pulse rate and pulse amplitude detected by the detection unit 10.
  • the estimation unit 203 performs interpolation (for example, based on the time when the pulse rate and pulse amplitude are detected by the detection unit 10, the time t, and the pulse rate and pulse amplitude detected by the detection unit 10. , Linear interpolation), the pulse rate b (t) and the pulse amplitude a (t) at time t are obtained.
  • the estimation unit 203 calculates the pulse rate and the pulse amplitude detected by the detection unit 10 at the time closest to the time t without performing interpolation, and calculates the pulse rate b (t) and the pulse amplitude a (t at the time t. ).
  • the determination unit 202 is based on the reference target resistance gi , 0 stored in the storage unit 201, Formula 9, and the pulse rate b (t) and the pulse amplitude a (t) acquired by the estimation unit 203. Te, to calculate a target resistance g i.
  • the determination unit 202 uses the reference target resistance g i, 0 stored in the storage unit 201, the reference target resistance ratio ⁇ i stored in the storage unit 201, Formula 9, and the estimation unit 203. and obtained pulse rate b (t) and pulse amplitude a (t), on the basis, to calculate a target resistance g i.
  • Calculating the target resistances g i is an example of determining a target resistance g i.
  • the estimation unit 203 calculates the cycle length ⁇ e at time t based on the acquired pulse rate b (t) and Equation 10.
  • Estimating unit 203 sets the periodic start time t 0 to time t.
  • Period start time t 0 is the time at which the period of the pulse for each pulse begins.
  • the estimation unit 203 calculates the in-period time ⁇ by subtracting the period start time t 0 from the time t.
  • the equation 11 is a differential equation for Equation 8, the target resistance g i calculated by the determination unit 202, based on, at time t + Delta] t after a lapse of step time Delta] t from the time t, 4 And resistances R 4 (t + ⁇ t) and R 8 (t + ⁇ t) for the eighth communication pipes FC4 and FC8 are calculated.
  • Calculation of the resistance R i is an example of a determination of the resistance R i.
  • the determination unit 202 includes the reference no-load volume fi , 0 stored in the storage unit 201, Equation 3, Equation 12, and the pulse rate b (t) and the pulse amplitude a (t acquired by the estimation unit 203. ) And the in-period time ⁇ (t) calculated by the estimation unit 203, the time differentiation ⁇ of the unloaded volumes V 2 and V 6 of the second and sixth containers FV2 and FV6 at the time t. 2 ( ⁇ ) and ⁇ 6 ( ⁇ ) are calculated.
  • Time differential [Phi i unloaded volume V i is an example of a change with time of the no-load volume V i.
  • Calculating the time derivative [Phi i is an example of a determination of the time derivative [Phi i.
  • the storage unit 201 may store a time derivative of the reference no-load volume fi , 0 instead of the reference no-load volume fi , 0 or in addition to the reference no-load volume fi , 0. Good.
  • the determination unit 202 calculates the time differential of the reference unloaded volume fi , 0 stored in the storage unit 201, the pulse rate b (t) and the pulse amplitude a (t) acquired by the estimation unit 203. Based on the in-period time ⁇ (t) calculated by the estimation unit 203, the time differentials ⁇ 2 ( 2) of the unloaded volumes V 2 and V 6 of the second and sixth containers FV2 and FV6 at time t. ⁇ ), ⁇ 6 ( ⁇ ) may be calculated.
  • the processing unit 20 instead of the storage of the time derivative of the reference no-load volume f i, 0, may store a function for calculating a time derivative of the reference no-load volume f i, 0.
  • the unloaded volume V of the containers FV1, FV3 to FV5, FV7, and FV8 other than the second and sixth containers FV2 and FV6 among the first to eighth containers FV1 to FV8. 1 , V 3 to V 5 , V 7 and V 8 do not change with time. Accordingly, the unloaded volumes V 1 , V 3 to V 5 of the containers FV1, FV3 to FV5, FV7, and FV8 other than the second and sixth containers FV2 and FV6 among the first to eighth containers FV1 to FV8. , V 7 , V 8 time derivatives ⁇ 1 ( ⁇ ), ⁇ 3 ( ⁇ ) to ⁇ 5 ( ⁇ ), ⁇ 7 ( ⁇ ), ⁇ 8 ( ⁇ ) are zero.
  • the estimation unit 203 is obtained by applying a fourth-order Runge-Kutta method to the differential equation represented by Equation 1, and after the step time ⁇ t has elapsed from time t, based on Equation 13 to Equation 21.
  • the pressures P 1 (t + ⁇ t) to P 8 (t + ⁇ t) of the fluid in the first to eighth containers FV1 to FV8 at time t + ⁇ t are calculated.
  • Calculation of the pressure P 1 (t + ⁇ t) ⁇ P 8 (t + ⁇ t) are an example of estimation of the pressure P 1 (t + ⁇ t) ⁇ P 8 (t + ⁇ t).
  • the estimation unit 203 After estimating the pressures P 1 (t + ⁇ t) to P 8 (t + ⁇ t), the estimation unit 203 updates the time t to the time t + ⁇ t after the step time ⁇ t has elapsed from the time t.
  • the estimation unit 203 includes a period ⁇ , resistances R 4 (t + ⁇ t), R 8 (t + ⁇ t), unloaded volumes V 2 , V 6 time derivatives ⁇ 2 ( ⁇ ), ⁇ 6 ( ⁇ ), and pressure P
  • the intra-period processing including the calculation of 1 (t + ⁇ t) to P 8 (t + ⁇ t) and the update of the time t is repeatedly executed while the intra-period time ⁇ is equal to or shorter than the period length ⁇ e .
  • Estimation unit 203 if the period within the time tau is greater than the period length tau e, acquisition of pulse rate b (t) and pulse amplitude a (t), calculating the target resistances g i and cycle length tau e, and Then, the cycle start time t 0 is set again, and then the intra-cycle processing is performed again.
  • the processing unit 20 stores the time t and the pressures P 1 (t) to P 8 (t) at the time t in association with each other.
  • the processing unit 20 outputs (for example, displays on the display) the latest blood pressure calculated via the output device 24 every time a predetermined display cycle (for example, 5 seconds) elapses.
  • the processing unit 20 outputs the pressure P 3 of the fluid in the third container FV3 as blood pressure.
  • the estimation unit 203 calculates the unloaded volume of the second and sixth containers FV2 and FV6 in the period between the in-cycle time ⁇ and the time ⁇ + ⁇ t after the step time ⁇ t has elapsed from the in-cycle time ⁇ .
  • a constant value for example, time derivative ⁇ i ( ⁇ ) at the time ⁇ in the cycle
  • V 2 and V 6 may be used as the time derivative of V 2 and V 6 .
  • the estimation unit 203 may use Formulas 22 to 25 instead of Formula 15, Formula 17, Formula 19, and Formula 21, respectively.
  • the estimator 203 determines the resistance to the fourth and eighth communication pipes FC4 and FC8 at time t + ⁇ t / 2 after the lapse of half ⁇ t / 2 of the step time ⁇ t from time t, based on Expression 26.
  • R 4 (t + ⁇ t / 2) and R 8 (t + ⁇ t / 2) may be calculated.
  • the detection unit 10 detects the pulse rate and the pulse amplitude every time a predetermined detection period (for example, 1 second) elapses. Further, the processing unit 20 executes the processing shown in FIG. Hereinafter, the processing of FIG. 8 will be described.
  • the processing unit 20 sets the time t, the pressures P 1 to P 8 , and the resistances R 1 to R 8 to the initial value t ini , the initial values P 1, ini to P 8, ini , and the initial value R 1, ini. To R8 and ini , respectively (step S101 in FIG. 8).
  • the processing unit 20 acquires the pulse rate b (t) and the pulse amplitude a (t) at time t based on the pulse rate and pulse amplitude detected by the detection unit 10 (step S102 in FIG. 8).
  • the process part 20 is based on the reference target resistance gi , 0 memorize
  • the processing unit 20 includes the reference target resistance g i, 0 stored in the storage unit 201, the reference target resistance ratio ⁇ i stored in the storage unit 201, and the pulse acquired in step S102. the number b (t) and pulse amplitude a (t), on the basis, to calculate a target resistance g i.
  • the processing unit 20 calculates the cycle length ⁇ e at time t based on the pulse rate b (t) acquired in step S102, and sets the cycle start time t 0 to time t (FIG. 8). Step S104).
  • step S105 the processing unit 20 from the time t, the value obtained by subtracting the periodic start time t 0 which is set in step S104, is calculated as the period in time tau (step S105 in FIG. 8).
  • the processing unit 20 determines whether or not the in-period time ⁇ calculated in step S105 is equal to or shorter than the period length ⁇ e calculated in step S104 (step S106 in FIG. 8).
  • the processing unit 20 determines “Yes” and proceeds to step S107 in FIG. Then, the processing unit 20 is one other than the fourth and eighth communication pipes FC4 and FC8 among the first to eighth communication pipes FC1 to FC8 at the time t + ⁇ t after the lapse of the step time ⁇ t from the time t. Resistances R 1 (t + ⁇ t) to R 3 (t + ⁇ t) and R 5 (t + ⁇ t) to R 7 (t + ⁇ t) for the communication pipes FC1 to FC3, FC5 to FC7 are changed to resistances R 1 (t) to R 3 (t ), R 5 (t) to R 7 (t), respectively.
  • the processing unit 20 sets the resistances R 4 (t + ⁇ t) and R 8 (t + ⁇ t) to the fourth and eighth communication pipes FC4 and FC8 at time t + ⁇ t after the lapse of step time ⁇ t from time t, at step S103. Is calculated based on the target resistance g i calculated in (1).
  • the processing unit 20 calculates the pressures P 1 (t + ⁇ t) to P 8 (t + ⁇ t) of the fluid in the first to eighth containers FV1 to FV8 at the time t + ⁇ t after the step time ⁇ t has elapsed from the time t. (Step S108 in FIG. 8). The process of step S108 will be described later.
  • the processing unit 20 updates the time t to a time t + ⁇ t after the step time ⁇ t has elapsed from the time t (step S109 in FIG. 8). Then, the processing unit 20 returns to step S105 of FIG. 8, the period in time tau until greater than the period length tau e, it repeats the process from step S105 to step S109.
  • processor 20 determines "No" at step S106 in FIG. 8, the flow returns to step S102 in FIG. 8. And the process part 20 performs the process after step S102 again with respect to the time t updated in step S109.
  • the processing unit 20 proceeds to step S102 in FIG. 8, the latest time among the times when the pulse rate and the pulse amplitude are detected by the detection unit 10 is updated at time t109. Wait until a future (future) time.
  • the processing unit 20 executes the processes shown in FIGS. 9 and 10 as the process of step S108 in FIG. Hereinafter, the processing of FIGS. 9 and 10 will be described.
  • the processing unit 20 uses the first provisional pressures P 1,0 to P 8,0 , which are provisional values of the pressures of the fluid in the first to eighth containers FV1 to FV8, at the time t.
  • the pressures P 1 (t) to P 8 (t) of the fluid in the containers FV1 to FV8 are set.
  • the processing unit 20 sets the first provisional pressures P 0,0 , P 9,0 to be used as the first provisional pressures P 8,0 , P 1,0 for convenience.
  • the processing unit 20 includes the first provisional pressures P 0,0 to P 9,0 , the resistances R 1 (t + ⁇ t) to R 9 (t + ⁇ t) calculated in step S107 in FIG. Based on the time derivatives ⁇ 1 ( ⁇ ) to ⁇ 8 ( ⁇ ) at ⁇ , the first time change rates K 1,1 to K 8,1 are calculated (step S201 in FIG. 9).
  • the processing unit 20 applies the second provisional pressures P 1,1 to P 8,1 , which are provisional values of the pressures of the fluid in the first to eighth containers FV1 to FV8, at the time t. 8 to the fluid pressures P 1 (t) to P 8 (t) in the containers FV1 to FV8, the first time rate of change K 1,1 to K 8,1 calculated in step S201 of FIG. A value obtained by adding a product of a value obtained by dividing the step time ⁇ t by 2 is set.
  • the processing unit 20 sets the second provisional pressures P 0,1 , P 9,1 to be used as the second provisional pressures P 8,1 , P 1,1 for convenience.
  • the processing unit 20 includes the second provisional pressures P 0,1 to P 9,1 , the resistances R 1 (t + ⁇ t) to R 9 (t + ⁇ t) calculated in step S107 in FIG.
  • the second time change rate K 1,2 to K 8,2 is calculated (step S202 in FIG. 9).
  • a third interim pressure P 1,2 ⁇ P 8,2 is a provisional value of the pressure of the fluid in the container FV1 ⁇ FV8 first to eighth, at time t, first to 8 to the fluid pressures P 1 (t) to P 8 (t) in the containers FV1 to FV8, the second time rate of change K 1,2 to K 8,2 calculated in step S202 of FIG.
  • a value obtained by adding a product of a value obtained by dividing the step time ⁇ t by 2 is set.
  • the processing unit 20 sets the third provisional pressures P 0,2 , P 9,2 used as the third provisional pressures P 8,2 , P 1,2 respectively.
  • the processing unit 20 includes the third provisional pressures P 0,2 to P 9,2 , the resistances R 1 (t + ⁇ t) to R 9 (t + ⁇ t) calculated in step S107 in FIG. Based on the time derivatives ⁇ 1 ( ⁇ + ⁇ t / 2) to ⁇ 8 ( ⁇ + ⁇ t / 2) at the time ⁇ + ⁇ t / 2 after half the step time ⁇ t from ⁇ , the third rate of change K 1,3 to K 8,3 is calculated (step S203 in FIG. 9).
  • the processing unit 20 applies the fourth provisional pressures P 1,3 to P 8,3 , which are provisional values of the pressures of the fluid in the first to eighth containers FV1 to FV8, at the time t. 8 to the fluid pressures P 1 (t) to P 8 (t) in the containers FV1 to FV8, the third rate of change K 1,3 to K 8,3 calculated in step S203 of FIG. Each is set to a value obtained by adding the product of the step time ⁇ t.
  • the processing unit 20 sets the fourth provisional pressures P 0,3 , P 9,3 to be used as the fourth provisional pressures P 8,3 , P 1,3 for convenience .
  • the processing unit 20 includes the fourth provisional pressures P 0,3 to P 9,3 , the resistances R 1 (t + ⁇ t) to R 9 (t + ⁇ t) calculated in step S107 in FIG. Based on the time differentials ⁇ 1 ( ⁇ + ⁇ t) to ⁇ 8 ( ⁇ + ⁇ t) at time ⁇ + ⁇ t after ⁇ by step time ⁇ t, fourth time change rates K 1,4 to K 8,4 are calculated (FIG. 10). Step S204).
  • the processing unit 20 calculates the pressures P 1 (t) to P 8 (t) of the fluid in the first to eighth containers FV1 to FV8 at time t and Steps S201 to S204, respectively.
  • Step S205 And the process part 20 complete
  • the blood pressure estimation device 1 detects the first parameter (in this example, the pulse rate b) representing the cycle length that is the length of the cycle of the heartbeat of the living body. Furthermore, the blood pressure estimation device 1 uses a plurality of containers FV1 to FV8 in a mathematical model that represents blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers FV1 to FV8 that are elastically deformed are connected in an annular shape.
  • the first parameter in this example, the pulse rate b
  • the blood pressure estimation device 1 uses a plurality of containers FV1 to FV8 in a mathematical model that represents blood flowing through a circulatory system of a living body by a fluid flowing through a flow path in which a plurality of containers FV1 to FV8 that are elastically deformed are connected in an annular shape.
  • the blood pressure estimation device 1 estimates the pressure of blood flowing through the circulatory system of the living body based on the determined change and the mathematical model.
  • the first parameter representing the cycle length is easily detected with higher accuracy than the artery diameter. Furthermore, by using the change with respect to time of the unloaded volumes V 2 and V 6 of the second and sixth containers FV2 and FV6 determined based on the first parameter, the mathematical model can be obtained from the left ventricle and the living body in the living body. The behavior of the right ventricle can be expressed with high accuracy. As a result, the mathematical model can express the behavior of blood flowing through the living body's circulatory system with high accuracy. Thus, the blood pressure estimation apparatus 1 can estimate the pressure P i which the blood has with high accuracy.
  • the blood pressure estimation device 1 has a pressure between containers of the pressure of the fluid in the pair of containers with respect to the flow rate of the fluid between the pair of containers connected to each other among the plurality of containers FV1 to FV8.
  • the resistor R 4 for the fourth and eighth communicating pipe FC4, FC8 of, R 8) resistance is the ratio of the difference to be determined based on the first parameter detected.
  • the blood pressure estimation device 1 estimates the pressure of blood flowing through the circulatory system of the living body based on the determined resistance.
  • the mathematical model can represent the behavior of the peripheral blood vessels in the living body with higher accuracy than when the resistances R 4 and R 8 are kept constant. Thereby, the mathematical model can express the behavior of blood flowing through the circulatory system of a living body with high accuracy. As a result, the blood pressure estimation apparatus 1 can estimate the pressure P i with blood at high accuracy.
  • the blood pressure estimation apparatus 1 of the first embodiment the resistance (in this example, resistors R 4, R 8 for the fourth and eighth communicating pipe FC4, FC8) of the target value (in this example, the target resistance g 4 , G 8 ) are determined based on the detected first parameter.
  • the blood pressure estimation device 1 determines the resistance so that the resistance approaches the determined target value with a delay.
  • the mathematical model can represent the behavior of the peripheral blood vessels in the living body with higher accuracy.
  • the mathematical model can express the behavior of blood flowing through the circulatory system of a living body with high accuracy.
  • the blood pressure estimation apparatus 1 the pressure P i with the blood flowing through the circulatory system of a living body can be estimated with high accuracy.
  • the blood pressure estimation device 1 of the first embodiment detects a second parameter (pulse amplitude a in this example) that represents the magnitude of the heartbeat.
  • a second parameter pulse amplitude a in this example
  • the blood pressure estimation apparatus 1 based on the second parameter detected, the resistance (in this example, resistors R 4, R 8 for the fourth and eighth communicating pipe FC4, FC8) of the target value of the (in this example Then, the target resistances g 4 and g 8 ) are determined.
  • the mathematical model can represent the behavior of the peripheral blood vessels in the living body with higher accuracy.
  • the mathematical model can express the behavior of blood flowing through the circulatory system of a living body with high accuracy.
  • the blood pressure estimation apparatus 1 the pressure P i with the blood flowing through the circulatory system of a living body can be estimated with high accuracy.
  • the blood pressure estimation device 1 of the first embodiment changes the unloaded volume with respect to time of at least one of the plurality of containers FV1 to FV8 based on the detected second parameter (in this example, the first parameter 2 and the time differentials ⁇ 2 and ⁇ 6 of the unloaded volumes V 2 and V 6 of the sixth containers FV2 and FV6).
  • the mathematical model can represent the behavior of the left ventricle and the right ventricle in the living body with higher accuracy.
  • the mathematical model can express the behavior of blood flowing through the circulatory system of a living body with high accuracy.
  • the blood pressure estimation apparatus 1 the pressure P i with the blood flowing through the circulatory system of a living body can be estimated with high accuracy.
  • the plurality of containers FV1 to FV8 include the left atrium, the left ventricle, the aorta and the artery downstream of the aorta, the vena cava and the vein upstream of the vena cava, the right Represents the atrium, right ventricle, pulmonary artery, and pulmonary vein, respectively.
  • the mathematical model can express the behavior of blood flowing through the circulatory system of a living body with high accuracy.
  • the blood pressure estimation apparatus 1 the pressure P i with the blood flowing through the circulatory system of a living body can be estimated with high accuracy.
  • a curve D11 in FIG. 11 represents an example of the pulse rate detected by the blood pressure estimation device 1 when the living body is at rest.
  • 12 to 14 show an example of an estimation result by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • Circles D12, D13, and D14 in FIG. 12 are the maximum blood pressure (in other words, the highest value of blood pressure), the lowest blood pressure (in other words, the lowest value of blood pressure), and the average blood pressure measured by the blood pressure measurement device as a comparative example. (In other words, the average value of blood pressure).
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • the average blood pressure is a value obtained by averaging the maximum blood pressure and the minimum blood pressure.
  • Curves E12, E13, and E14 in FIG. 12 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively. As shown in FIG. 12, according to the blood pressure estimation device 1 of the first embodiment, blood pressure can be estimated with high accuracy.
  • Curve E15, E16 in FIG. 13 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D18 in FIG. 14 represents the maximum blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E18 in FIG. 14 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 14, according to the blood pressure estimation device 1 of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • 15 represents an example of the pulse rate detected by the blood pressure estimation device 1 during the movement of the living body.
  • 16 to 18 show an example of an estimation result by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D22, D23, and D24 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • FIG. 16 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation apparatus 1, respectively. As shown in FIG. 16, according to the blood pressure estimation device 1 of the first embodiment, blood pressure can be estimated with high accuracy.
  • Curve E25, E26 in FIG. 17 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D28 in FIG. 18 represents the systolic blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E28 in FIG. 18 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 18, according to the blood pressure estimation device 1 of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • the blood pressure estimation device 1 may include a first device that constitutes the detection unit 10 and a second device that constitutes the processing unit 20.
  • the first device and the second device are communicably connected to each other.
  • the second device may be a mobile phone, a smartphone, a personal computer, or a server device.
  • the blood pressure estimation device 1 performs the detection of the pulse rate and the estimation of the blood pressure at the time when the pulse rate is detected in parallel. By the way, after the detection of the pulse rate in the predetermined period is completed, the blood pressure estimation apparatus 1 may estimate the blood pressure in the period.
  • the blood pressure estimation device 1 stores each of a plurality of different times and the pressure at the time in association with each other.
  • the blood pressure estimation device 1 may output each of a plurality of stored times in association with the stored pressure at that time.
  • the blood pressure estimation device 1 may output a graph representing changes in pressure with respect to time.
  • the processing unit 20 sets the resistances R 4 and R 8 for the fourth and eighth communication pipes FC 4 and FC 8 to initial values R 4, ini , R 8, and ini , respectively, and changes them with time. Instead, the blood pressure may be estimated.
  • the blood pressure estimation device 1 receives information representing a pulse rate and blood pressure measured in advance with respect to the living body, and based on the input information, the reference target resistance g i, 0 , the initial value of the resistance At least one of R i, ini and coefficient E i may be determined.
  • a blood pressure estimation device according to a first modification of the first embodiment will be described.
  • the blood pressure estimation device according to the first modification of the first embodiment is different from the blood pressure estimation device according to the first embodiment in that a second parameter representing the magnitude of a heartbeat of a living body is estimated based on the first parameter. is doing.
  • the difference will be mainly described.
  • 1st modification of 1st Embodiment what attached
  • the processing unit 20 estimates the pulse amplitude a based on the pulse rate b detected by the detection unit 10 and the mathematical expression 27. Further, the processing unit 20 determines the unloaded volumes f 2 and f 6 of the second and sixth containers FV2 and FV6 based on the estimated pulse amplitude a, and the fourth and eighth communication pipes FC4. , to calculate the target resistance g 4, g 8 for FC8.
  • T a represents the time constant of the change of the pulse amplitude a.
  • the time constant Ta may be understood as representing the degree to which the change in the pulse amplitude a is delayed with respect to the change in the target amplitude h.
  • the pulse amplitude a may be grasped as approaching the target amplitude h with a delay.
  • the target amplitude h has a value determined in advance according to the pulse rate ratio b / b 0 . In this example, the target amplitude h increases as the pulse rate ratio b / b 0 increases, as represented by FIG.
  • the blood pressure estimation device 1 may determine the target amplitude h based on the information indicating the pulse rate and blood pressure measured in advance with respect to the living body and the input information.
  • the estimation unit 203 acquires the pulse rate b (t) at the time t based on the pulse rate detected by the detection unit 10. Furthermore, the estimation unit 203 calculates the pulse amplitude a at time t + ⁇ t after the step time ⁇ t has elapsed from time t, based on Equation 28, which is a difference equation with respect to Equation 27, and the acquired pulse rate b (t). (T + ⁇ t) is calculated.
  • the calculation of the pulse amplitude a is an example of the estimation of the pulse amplitude a.
  • the detection unit 10 does not have to detect the second parameter indicating the magnitude of the heartbeat of the living body.
  • the same operations and effects as the blood pressure estimation device 1 of the first embodiment are exhibited.
  • 20 to 22 show an example of an estimation result obtained by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D32, D33, and D34 respectively represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure measured by the blood pressure measuring device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E32, E33, and E34 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E35, E36 in FIG. 21 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D38 in FIG. 22 represents the systolic blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E38 in FIG. 22 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 22, according to the blood pressure estimation device 1 of the first modification of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • 23 to 25 show examples of estimation results by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D42, D43, and D44 respectively represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure measured by the blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E42, E43, and E44 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively. As shown in FIG. 23, according to the blood pressure estimation device 1 of the first modification of the first embodiment, blood pressure can be estimated with high accuracy.
  • Curve E45, E46 in FIG. 24 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D48 in FIG. 25 represents the systolic blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E48 in FIG. 25 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 25, according to the blood pressure estimation device 1 of the first modified example of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • the blood pressure estimation apparatus uses a predetermined function to compare the blood pressure estimation apparatus according to the first modification of the first embodiment with a reference target resistance g i, 0 and a reference target resistance.
  • the difference is that the ratio ⁇ i , the resistance R 4,0 to the fourth communication pipe FC4 when the living body is at rest, and the coefficient E 4 associated with the fourth container FV4 are calculated.
  • the difference will be mainly described.
  • symbol used in the 1st modification of 1st Embodiment is the same or substantially the same.
  • Processing unit 20 instead of the storage of the reference target resistance ratio gamma i, and stores a function for calculating a reference target resistance ratio gamma i.
  • the processing unit 20 stores a function represented by Expression 29 as a function for calculating the reference target resistance ratio ⁇ i .
  • the processing unit 20 calculates the reference target resistance ratio ⁇ i based on Expression 29.
  • the processing unit 20 instead of storing the reference target resistance g i, 0, and stores the function for calculating the reference target resistance g i, 0.
  • the processing unit 20 stores a function represented by Expression 30 as a function for calculating the reference target resistance gi , 0 .
  • the processing unit 20 calculates the reference target resistance g i, 0 based on Equation 30.
  • R i, 0 represents the resistance to the i th communication pipe FCi when the living body is at rest.
  • the processing unit 20 sets the resistance R 4,0 to the fourth communication pipe FC4 when the living body is resting, the reference value b 0 of the pulse rate b, and the maximum blood pressure P smax measured when the living body is resting. And based on the minimum blood pressure P smin measured when the living body is resting, the predetermined coefficient ⁇ , and Equation 31.
  • the reference value b 0 of the pulse rate b is a pulse rate measured when the living body is at rest.
  • the coefficient ⁇ may be changed every time blood pressure estimation is executed.
  • the coefficient ⁇ may be set to a value according to the living body.
  • the processing unit 20 uses the coefficient E 4 associated with the fourth container FV4 as the maximum blood pressure P smax measured when the living body is resting and the minimum blood pressure P smin measured when the living body is resting. , Based on a predetermined coefficient ⁇ and Expression 32.
  • the coefficient ⁇ may be changed every time blood pressure estimation is executed. Also, the coefficient ⁇ may be set to a value according to the living body.
  • the processing unit 20 uses a fixed value (in other words, a value common to a plurality of different living bodies) as the resistance R 8,0 for the eighth communication pipe FC8.
  • the resistance R 8,0 to the eighth communication pipe FC8 may be changed every time blood pressure estimation is executed. Further, the resistance R 8,0 for the eighth communication pipe FC8 may be set to a value corresponding to the living body.
  • the processing unit 20 uses a fixed value as the coefficient E j associated with the jth container FVj.
  • j represents each integer other than 4 among 1 to 8 integers.
  • the coefficient E j associated with the jth container FVj may be changed every time blood pressure estimation is executed.
  • the coefficient E j associated with the jth container FVj may be set to a value corresponding to the living body.
  • the same operations and effects as the blood pressure estimation device 1 of the first modification of the first embodiment are exhibited.
  • 26 to 29 show an example of an estimation result by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D92, D93, and D94 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E92, E93, and E94 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E95, E96 in FIG. 27 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D98 in FIG. 28 represents the maximum blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E98 in FIG. 28 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 28, according to the blood pressure estimation device 1 of the second modification of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • the blood pressure measured by the blood pressure measuring device is represented as measured blood pressure
  • the blood pressure estimated by the blood pressure estimating device 1 is represented as estimated blood pressure.
  • a square F91 in FIG. 29 represents the relationship between the maximum blood pressure measured by the blood pressure measurement device as a comparative example and the maximum blood pressure estimated by the blood pressure estimation device 1.
  • a triangle F92 in FIG. 29 represents the relationship between the minimum blood pressure measured by the blood pressure measurement device as a comparative example and the minimum blood pressure estimated by the blood pressure estimation device 1.
  • a circle F93 in FIG. 29 represents the relationship between the average blood pressure measured by the blood pressure measurement device as a comparative example and the average blood pressure estimated by the blood pressure estimation device 1.
  • the maximum blood pressure, the minimum blood pressure, and the average blood pressure can be estimated with high accuracy.
  • a blood pressure estimation device according to a third modification of the first embodiment will be described.
  • the blood pressure estimation device according to the third modification of the first embodiment is different from the blood pressure estimation device according to the first embodiment in that the blood pressure is estimated without being based on the second parameter representing the magnitude of the heartbeat of the living body. ing.
  • the difference will be mainly described.
  • the 3rd modification of 1st Embodiment what attached
  • the target value g 4, g 8 of resistors R 4, R 8 for the fourth and eighth communicating pipe FC4, FC8 of is represented by equation 34 in place of Equation 9.
  • the detection unit 10 does not have to detect the second parameter indicating the magnitude of the heartbeat of the living body.
  • the same operations and effects as the blood pressure estimation device 1 of the first embodiment are exhibited.
  • 30 to 32 show an example of an estimation result by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D52, D53, and D54 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E52, E53, and E54 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E55, E56 in FIG. 31 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D58 in FIG. 32 represents the systolic blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E58 in FIG. 32 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 32, according to the blood pressure estimation device 1 of the third modification of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • 33 to 35 show an example of an estimation result obtained by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D62, D63, and D64 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E62, E63, and E64 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E65, E66 in FIG. 34 represents a fourth target value g 4 of the resistor R 4 for communicating conduit FC4, and the resistance R 4 for the fourth communicating pipe FC4, respectively.
  • a square D68 in FIG. 35 represents the maximum blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E68 in FIG. 35 represents the systolic blood pressure estimated by the blood pressure estimation device 1. As shown in FIG. 35, according to the blood pressure estimation device 1 of the third modification of the first embodiment, the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • a blood pressure estimation device according to a fourth modification of the first embodiment will be described.
  • the blood pressure estimation apparatus according to the fourth modification of the first embodiment is different from the blood pressure estimation apparatus according to the third modification of the first embodiment in that blood pressure is estimated without changing the resistance with respect to time. ing.
  • the difference will be mainly described.
  • symbol used in 1st Embodiment is the same or substantially the same.
  • the processing unit 20 sets the resistances R 4 and R 8 for the fourth and eighth communication pipes FC 4 and FC 8 to the initial values R 4, ini , R 8, and ini , respectively, and does not change them with respect to time. According to the blood pressure estimation device 1 of the fourth modified example of the first embodiment, the same operations and effects as the blood pressure estimation device 1 of the first embodiment are exhibited.
  • 36 to 38 show an example of an estimation result by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D72, D73, and D74 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E72, E73, and E74 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E76 in FIG. 37 represents the resistance R 4 for the fourth communicating pipe FC4.
  • a square D78 in FIG. 38 represents the maximum blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E78 in FIG. 38 represents the systolic blood pressure estimated by the blood pressure estimation device 1.
  • the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.
  • 39 to 41 show examples of estimation results by the blood pressure estimation device 1 when the pulse rate detected by the blood pressure estimation device 1 changes as shown in FIG.
  • circles D82, D83, and D84 represent the systolic blood pressure, the diastolic blood pressure, and the average blood pressure, respectively, measured by a blood pressure measurement device as a comparative example.
  • the blood pressure measurement device measures blood pressure by using a cuff.
  • curves E82, E83, and E84 represent the maximum blood pressure, the minimum blood pressure, and the average blood pressure estimated by the blood pressure estimation device 1, respectively.
  • blood pressure can be estimated with high accuracy.
  • Curve E86 in FIG. 40 represents a resistance R 4 for the fourth communicating pipe FC4.
  • a square D88 in FIG. 41 represents the systolic blood pressure measured by the blood pressure measurement device as a comparative example.
  • a circle E88 in FIG. 41 represents the systolic blood pressure estimated by the blood pressure estimation device 1.
  • the relationship between the maximum blood pressure and the pulse rate can be estimated with high accuracy.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Cardiology (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Hematology (AREA)

Abstract

La présente invention concerne un dispositif d'estimation de la pression artérielle (1) comprenant : une unité de détection (10) qui détecte un premier paramètre indiquant une longueur de période qui représente une période de battement cardiaque d'un corps vivant; et une unité de traitement (20) qui, dans un modèle mathématique où le sang circulant dans le système circulatoire du corps vivant est représenté par un fluide s'écoulant dans un passage d'écoulement formé par la liaison d'une pluralité de vaisseaux à déformation élastique dans une configuration annulaire, détermine une modification temporelle du volume d'au moins un des vaisseaux sur la base du premier paramètre détecté et estime ensuite la pression du sang circulant sur la base du modèle mathématique et de la modification déterminée.
PCT/JP2016/083770 2015-11-17 2016-11-15 Dispositif d'estimation de la pression artérielle, dispositif de traitement d'informations, procédé d'estimation de la pression artérielle et programme d'estimation de la pression artérielle Ceased WO2017086294A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2017551876A JP6580158B2 (ja) 2015-11-17 2016-11-15 血圧推定装置、情報処理装置、血圧推定方法、及び、血圧推定プログラム
US15/981,075 US11284805B2 (en) 2015-11-17 2018-05-16 Blood pressure estimating apparatus, method for estimating blood pressure, and non-transitory computer-readable recording medium having stored therein program for estimating blood pressure

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015-224729 2015-11-17
JP2015224729 2015-11-17

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/981,075 Continuation US11284805B2 (en) 2015-11-17 2018-05-16 Blood pressure estimating apparatus, method for estimating blood pressure, and non-transitory computer-readable recording medium having stored therein program for estimating blood pressure

Publications (1)

Publication Number Publication Date
WO2017086294A1 true WO2017086294A1 (fr) 2017-05-26

Family

ID=58718955

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/083770 Ceased WO2017086294A1 (fr) 2015-11-17 2016-11-15 Dispositif d'estimation de la pression artérielle, dispositif de traitement d'informations, procédé d'estimation de la pression artérielle et programme d'estimation de la pression artérielle

Country Status (3)

Country Link
US (1) US11284805B2 (fr)
JP (1) JP6580158B2 (fr)
WO (1) WO2017086294A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6847489B1 (ja) * 2019-11-05 2021-03-24 国立大学法人東北大学 血圧推定装置、血圧推定方法、及び、血圧推定プログラム
WO2021090779A1 (fr) * 2019-11-05 2021-05-14 国立大学法人東北大学 Dispositif d'estimation de la tension artérielle, procédé d'estimation de la tension artérielle, et programme d'estimation de la tension artérielle

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021044571A1 (fr) * 2019-09-05 2021-03-11 株式会社日立国際電気 Procédé, dispositif et système d'estimation de la fréquence du pouls
DE102019214343A1 (de) * 2019-09-20 2021-03-25 Robert Bosch Gmbh Verfahren zum Ausgleichen von Ladezuständen eines elektrischen Energiespeichers

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007029137A (ja) * 2005-07-22 2007-02-08 Seiko Instruments Inc 血液レオロジー測定装置
JP2008512177A (ja) * 2004-09-13 2008-04-24 ガンブロ・ルンディア・エービー 急激な血圧変化の検知
WO2013005320A1 (fr) * 2011-07-06 2013-01-10 パイオニア株式会社 Dispositif de prédiction d'hypotension
JP2013132407A (ja) * 2011-12-27 2013-07-08 Tohoku Univ 連続血圧測定システム、連続血圧測定方法、連続血圧測定プログラムおよび、連続血圧測定プログラムを記録したコンピュータ読み取り可能な記録媒体

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3893452A (en) * 1974-07-09 1975-07-08 Medtronic Inc Blood pressure monitoring system
US5423322A (en) * 1988-12-29 1995-06-13 Medical Physics, Inc. Total compliance method and apparatus for noninvasive arterial blood pressure measurement
US5140990A (en) * 1990-09-06 1992-08-25 Spacelabs, Inc. Method of measuring blood pressure with a photoplethysmograph
JP3465133B2 (ja) * 1997-08-05 2003-11-10 日本光電工業株式会社 患者監視装置
US6413223B1 (en) 1999-06-01 2002-07-02 Massachussetts Institute Of Technology Cuffless continuous blood pressure monitor
CA2388033A1 (fr) * 1999-10-07 2001-04-12 Alexander K. Mills Procede non invasif de determination en continu de caracteristiques physiologiques et dispositif a cet effet
JP4206218B2 (ja) * 2002-04-03 2009-01-07 セイコーインスツル株式会社 循環動態測定装置
JP4581050B2 (ja) * 2006-02-16 2010-11-17 財団法人ヒューマンサイエンス振興財団 心臓酸素消費量自動最小化システムおよびこれを用いた心疾患治療システム
US7654964B1 (en) * 2006-03-16 2010-02-02 Pacesetter, Inc. System and method for detecting arterial blood pressure based on aortic electrical resistance using an implantable medical device
WO2008144404A1 (fr) * 2007-05-16 2008-11-27 Massachusetts Instutute Of Technology Systèmes et méthodes d'estimation basée sur des modèles du débit cardiaque et de la résistance vasculaire périphérique totale
EP2306896B1 (fr) * 2008-07-08 2021-04-21 The Johns Hopkins University Systèmes non invasifs pour évaluer la pression de remplissage cardiaque
GB201119982D0 (en) * 2011-11-18 2012-01-04 Maarek Albert Medical device and measurement of cardiac output
US9522317B2 (en) * 2011-08-19 2016-12-20 Pulson, Inc. Systems and methods for coordinating musculoskeletal and cardiovascular or cerebrovascular hemodynamics

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008512177A (ja) * 2004-09-13 2008-04-24 ガンブロ・ルンディア・エービー 急激な血圧変化の検知
JP2007029137A (ja) * 2005-07-22 2007-02-08 Seiko Instruments Inc 血液レオロジー測定装置
WO2013005320A1 (fr) * 2011-07-06 2013-01-10 パイオニア株式会社 Dispositif de prédiction d'hypotension
JP2013132407A (ja) * 2011-12-27 2013-07-08 Tohoku Univ 連続血圧測定システム、連続血圧測定方法、連続血圧測定プログラムおよび、連続血圧測定プログラムを記録したコンピュータ読み取り可能な記録媒体

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6847489B1 (ja) * 2019-11-05 2021-03-24 国立大学法人東北大学 血圧推定装置、血圧推定方法、及び、血圧推定プログラム
WO2021090779A1 (fr) * 2019-11-05 2021-05-14 国立大学法人東北大学 Dispositif d'estimation de la tension artérielle, procédé d'estimation de la tension artérielle, et programme d'estimation de la tension artérielle

Also Published As

Publication number Publication date
US20180263514A1 (en) 2018-09-20
JP6580158B2 (ja) 2019-09-25
JPWO2017086294A1 (ja) 2018-09-06
US11284805B2 (en) 2022-03-29

Similar Documents

Publication Publication Date Title
JP6580158B2 (ja) 血圧推定装置、情報処理装置、血圧推定方法、及び、血圧推定プログラム
JP5904660B2 (ja) 動脈圧波形データを用いた血管症状の検出
JP6494664B2 (ja) 血圧測定デバイスの連続キャリブレーション
CN104116503B (zh) 一种无创连续血压的测量装置
CN108245148B (zh) 自适应生物信号特征组合设备及方法
JP2002253519A5 (fr)
JP2017521115A5 (fr)
EP2393423A2 (fr) Calcul de paramètres cardiovasculaires
JP2016508386A5 (fr)
CN107847649A (zh) 用于处理心脏传感器输出的方法和计算机系统
CN109872820A (zh) 一种无袖带血压测量方法、装置、设备及存储介质
CN107898443A (zh) 重搏波检测方法、装置及计算机存储介质
EP3019075B1 (fr) Détermination d'un paramètre hémodynamique
JP2020505177A (ja) 向上した心拍出量測定を提供する血流力学モニター
WO2015006191A1 (fr) Détermination de paramètres cardiovasculaires sur la base de la pression artérielle
JP6835226B2 (ja) 特定装置、特定方法、及び、特定プログラム
CA2756048A1 (fr) Surveillance du decouplage peripherique
EP4329599B1 (fr) Détermination de la pression sanguine
JP5940652B2 (ja) 体調監視装置及び方法
JP6847489B1 (ja) 血圧推定装置、血圧推定方法、及び、血圧推定プログラム
WO2016039241A1 (fr) Procédé d'estimation d'intensité de course favorable pour un dispositif d'assistance d'entrainement et de course
JP6806254B2 (ja) 特定装置、特定方法、及び、特定プログラム
WO2021090779A1 (fr) Dispositif d'estimation de la tension artérielle, procédé d'estimation de la tension artérielle, et programme d'estimation de la tension artérielle
CN118267604B (zh) 基于左心室辅助设备的剩余心脏收缩力估计方法及系统
US20210145298A1 (en) Methods, apparatus, computer programs, systems for calculating a pulse wave velocity of a subject

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16866296

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017551876

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16866296

Country of ref document: EP

Kind code of ref document: A1