CN107296593A - A kind of hemodynamic parameter acquisition methods and device - Google Patents

A kind of hemodynamic parameter acquisition methods and device Download PDF

Info

Publication number
CN107296593A
CN107296593A CN201710398262.1A CN201710398262A CN107296593A CN 107296593 A CN107296593 A CN 107296593A CN 201710398262 A CN201710398262 A CN 201710398262A CN 107296593 A CN107296593 A CN 107296593A
Authority
CN
China
Prior art keywords
electrocardiosignal
photoplethysmographic
signal
ripples
point
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.)
Pending
Application number
CN201710398262.1A
Other languages
Chinese (zh)
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.)
Jiangsu Standard Medical Technology Co Ltd
Original Assignee
Jiangsu Standard Medical Technology Co Ltd
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 Jiangsu Standard Medical Technology Co Ltd filed Critical Jiangsu Standard Medical Technology Co Ltd
Priority to CN201710398262.1A priority Critical patent/CN107296593A/en
Publication of CN107296593A publication Critical patent/CN107296593A/en
Pending legal-status Critical Current

Links

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/02028Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/023Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart

Abstract

The invention provides a kind of hemodynamic parameter acquisition methods and device.This method includes:Gather the electrocardiosignal and photoplethysmographic signal of targeted customer;Electrocardiosignal and photoplethysmographic signal are pre-processed;Pretreated electrocardiosignal Q ripples and R ripples are detected, and feature point identification is carried out to pretreated photoplethysmographic signal;Multiple hemodynamic parameters are obtained according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and photoplethysmographic signal.The embodiment of the present invention belongs to continuous non-invasive hemodynamic monitoring, the logic analysis of systematic science is carried out by the electrocardiosignal and photoplethysmographic signal that are produced to targeted customer in the daily lifes such as tranquillization, motion, accurately obtain multiple hemodynamic parameters, compared to existing continuous non-invasive hemodynamic monitoring, measurement is accurate, it is adapted to the specific demands, more safety and comfort such as operating room, Intensive Care Therapy, daily daily life, motion, the limit, antijamming capability is stronger.

Description

A kind of hemodynamic parameter acquisition methods and device
Technical field
The present invention relates to technical field of data processing, and in particular to a kind of hemodynamic parameter acquisition methods and device.
Background technology
Hemodynamics is using blood flow and heart, blood vessel distortion as the subject of research object.Existing hemodynamics Hygienic monitoring on hands of childhood method is broadly divided into invasive monitoring and the major class of non-invasive monitoring two.
Existing invasive monitoring can obtain multiple hemodynamic parameters, but inconvenient for use, expensive, exist certain Infection risk, bring pain to testee.It is main in existing non-invasive monitoring to use impedance cardiography, it is necessary to once Property electrode transmission high frequency electric passes through thoracic cavity by a narrow margin, device structure is complicated, has the electromagnetic radiation of doses to being monitored person, damages The health of monitored person is done harm to, comfortableness, anti-interference are poor, be not suitable for long-term quiescent condition monitoring and in daily life, fortune Monitored in dynamic, limit the application and popularizations of technology.
The content of the invention
The embodiment of the present invention provides a kind of hemodynamic parameter acquisition methods and device, for solving how to provide accurate The problem of obtaining the method for multiple hemodynamic parameters.
The embodiments of the invention provide a kind of hemodynamic parameter acquisition methods, including:
Gather the electrocardiosignal and photoplethysmographic signal of targeted customer;
The electrocardiosignal and the photoplethysmographic signal are pre-processed;
Pretreated electrocardiosignal Q ripples and R ripples are detected, and pretreated photoplethysmographic signal is carried out special Levy point identification;
Obtain multiple according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Hemodynamic parameter.
Alternatively, the electrocardiosignal and the photoplethysmographic signal are pre-processed, including:
The baseline and noise of the electrocardiosignal are removed, it is average smooth to removing the electrocardiosignal progress after baseline and noise Processing, to removing baseline and noise, carrying out the electrocardiosignal after average smooth treatment, carries out first-order difference processing, to a jump Electrocardiosignal after the reason of office does average smooth treatment, obtains the smooth signal of first-order difference of the electrocardiosignal;
Average smooth treatment is carried out to the photoplethysmographic signal, to the photoelectricity volume arteries and veins after average smooth treatment Ripple signal of fighting carries out first-order difference processing, and the photoplethysmographic signal after average smooth treatment and first-order difference processing is done Average smooth treatment, obtains the smooth signal of first-order difference of the photoplethysmographic, by the photoplethysmographic The smooth signal of first-order difference is moved to the left preset time, and the preset time is propagated for the pulse wave of the targeted customer from heart To the time needed for the photoplethysmographic signal collection point.
Alternatively, the characteristic point of the pretreated photoplethysmographic signal of mark includes:
Photoplethysmographic maximum slope point.
Alternatively, obtained according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Multiple hemodynamic parameters are taken, including:
According to the R ripples crest location of the electrocardiosignal and the position acquisition of the photoplethysmographic maximum slope point Blood systolic pressure and blood diastolic pressure;
The stroke output and external perihaemal canal resistance of heart are obtained according to the blood systolic pressure and the blood diastolic pressure;
Cardiac output is obtained according to the stroke output of the heart;
The left heart acting is obtained according to the blood systolic pressure, the blood diastolic pressure and the cardiac output;
Arterial compliance is obtained according to the blood systolic pressure, the blood diastolic pressure and the stroke output.
Alternatively,
Blood systolic pressure and blood diastolic pressure are obtained according to equation below:
SBP=a1+a2·PWV
DBP=a3+a4·PWV
Wherein, SBP represents blood systolic pressure, and DBP represents blood diastolic pressure;PWV represents pulse wave velocity, according to such as Lower formula is obtained:
Height is the height of the targeted customer;PWTT represents the R ripples crest and the photocapacitance of the electrocardiosignal Time interval between product pulse wave maximum slope point;a1、a2、a3And a4It is the constant related to the targeted customer;According to such as Lower formula obtains the stroke output of heart:
Wherein, SV represents the stroke output of heart;Characteristic quantity K=1/3, HR represent heart rate;a5And a6It is and the target The related constant of user;
External perihaemal canal resistance is obtained according to equation below:
SVR=a10+a11·(SBP-DBP)+a12·HR+a13·Age+a14·K
Wherein, SVR represents external perihaemal canal resistance;Age represents the age of the targeted customer;a10、a11、a12、a13And a14 It is the constant related to the targeted customer;Cardiac output is obtained according to equation below:
Wherein, CO represents cardiac output;
Left heart acting is obtained according to equation below:
LCW=a15·(MAP-a16)·CO
Wherein, LCW represents left heart acting;a15、a16It is the constant related to the targeted customer;MAP represents average blood Pressure, MAP=(2DSP+SBP)/3;
Arterial compliance is obtained according to equation below:
C=SV (SBP-DBP)
Wherein, C represents arterial compliance.
Alternatively, the characteristic point of the pretreated photoplethysmographic signal of mark also includes:
Point that slope before sequences of ventricular depolarization point, photoplethysmographic maximum slope point is zero, photoplethysmographic Trailing edge slope smallest point.
Alternatively, obtained according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Multiple hemodynamic parameters are taken, including:
The sequences of ventricular depolarization point and photoplethysmographic slope determined according to the original position of the Q ripples of the electrocardiosignal of detection Time interval between the point that slope before maximum point is zero determines pre-ejection;
According to the slope before the photoplethysmographic maximum slope point be zero point and photoplethysmographic under Drop determines left heart ejection time along the time interval between slope smallest point.
The embodiment of the present invention provides a kind of hemodynamic parameter acquisition device, including:
Signal gathering unit, electrocardiosignal and photoplethysmographic signal for gathering targeted customer;
Pretreatment unit, for being pre-processed to the electrocardiosignal and the photoplethysmographic signal;
Processing unit, for detecting pretreated electrocardiosignal Q ripples and R ripples, and to pretreated photoelectricity volume arteries and veins Ripple signal of fighting carries out feature point identification;
Hemodynamic parameter acquiring unit, for according to the pretreated electrocardiosignal and the photoelectricity volume arteries and veins The characteristic point of ripple signal of fighting obtains multiple hemodynamic parameters.
The embodiment of the present invention provides a kind of electronic equipment, including:Processor, memory and bus;Wherein,
Processor and memory complete mutual communication by bus;
Processor is used to call the programmed instruction in memory, to perform above-mentioned hemodynamic parameter acquisition methods.
The embodiment of the present invention provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage Medium storing computer is instructed, and the computer instruction makes the computer perform above-mentioned hemodynamic parameter acquisition side Method.
Hemodynamic parameter acquisition methods and device provided in an embodiment of the present invention, gather the electrocardiosignal of targeted customer And photoplethysmographic signal;The electrocardiosignal and the photoplethysmographic signal are pre-processed;Detection is pre- Electrocardiosignal Q ripples and R ripples after processing, and feature point identification is carried out to pretreated photoplethysmographic signal;According to The electrocardiosignal Q ripples of detection and the characteristic point of R ripples and the photoplethysmographic signal obtain multiple Hemodynamics ginsengs Number.The embodiment of the present invention belongs to movable type continuous non-invasive hemodynamic monitoring, proposes that photoplethysmographic signal is included There are VPV and volume change information, believed by the electrocardio of the generation to targeted customer in the daily lifes such as tranquillization, motion Number and photoplethysmographic signal carry out systematic science logic analysis, accurately obtain multiple hemodynamic parameters, compare Existing noninvasive hemodynamics monitoring, it is to avoid the electromagnetic radiation to targeted customer, reduces equipment power dissipation, more safety is relaxed Suitable, antijamming capability is stronger.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing There is the accompanying drawing used required in technology description to do simple introduction, it should be apparent that, drawings in the following description are the present invention Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can also basis These accompanying drawings obtain other accompanying drawings.
Fig. 1 is the schematic flow sheet of the hemodynamic parameter acquisition methods of one embodiment of the invention;
Fig. 2 is the schematic diagram of the hemodynamic parameter acquisition methods of one embodiment of the invention;
Fig. 3 is the structural representation of the movable many life sign monitor systems of type continuous non-invasive of one embodiment of the invention Figure;
Fig. 4 is the structural representation of the hemodynamic parameter acquisition device of one embodiment of the invention;
Fig. 5 is the entity structure schematic diagram of the electronic equipment of one embodiment of the invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention In accompanying drawing, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Fig. 1 is the schematic flow sheet of the hemodynamic parameter acquisition methods of one embodiment of the invention.As shown in figure 1, The method of the embodiment includes:
S11:Gather the electrocardiosignal and photoplethysmographic signal of targeted customer;
Containing abundant human body in daily lifes such as tranquillization, motions in electrocardiosignal (Electrocardiograph, ECG) Physiological data in activity.Electrocardiosignal Q ripples and R ripples (Q and R represent two home positions of ecg wave form respectively) detection are can Occupy critical role in many life sign monitor systems of sports type continuous non-invasive.In actual applications, the targeted customer of collection Electrocardiosignal is the three lead electrocardiosignals containing powerful motor electrical noise.
Photoplethysmographic (Photoplethysmogram of a Pulse Wave, PPW) contain Hemodynamics, All multi informations such as heartbeat function, vegetative nerve function, refer to the light intensity after two-beam reflects or reflected through tissue blood with the time Change corresponding relation, cardiac output continuous monitoring can be carried out by extracting blood flow information therein.
In actual applications, (wavelength is respectively the electrocardiosignal and two-beam Power Capacity pulse wave of collection targeted customer 660nm feux rouges and 940nm infrared light) synchronizing signal.
S12:The electrocardiosignal and the photoplethysmographic signal are pre-processed;
In actual applications, the pretreatment to electrocardiosignal includes removing the baseline and motion artifacts of electrocardiosignal, is averaged Smooth treatment, calculating first-order difference simultaneously do average smooth treatment to first-order difference signal;To the pre- of photoplethysmographic signal Processing includes average smooth treatment and calculates first-order difference and do average smooth treatment to first-order difference signal.
S13:Pretreated electrocardiosignal Q ripples and R ripples are detected, and pretreated photoplethysmographic signal is entered Row feature point identification;
S14:Obtained according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Multiple hemodynamic parameters;
It should be noted that impedance electrocardio is that high frequency electric is passed by thoracic aorta according to resistance minimal path by a narrow margin Lead, each heartbeat, due to VPV in sustainer and volume change, measure impedance.Photoplethysmographic is in tested portion When the intravascular VPV in position and volume change, transmission or reflection luminous intensity is measured.Impedance electrocardio and photoplethysmographic are all Include VPV and volume change information;The embodiment of the present invention passes through the analysis to electrocardiosignal and photoplethysmographic Hemodynamic parameter is obtained, without high frequency electric, by thoracic cavity, enhances interference free performance by a narrow margin in the prior art, it is to avoid Electromagnetic radiation to targeted customer, reduces equipment power dissipation.
Hemodynamic parameter acquisition methods provided in an embodiment of the present invention belong to movable type continuous non-invasive hemodynamics Hygienic monitoring on hands of childhood, propose photoplethysmographic signal include VPV and volume change information, by targeted customer quiet The electrocardiosignal and photoplethysmographic signal of generation in the daily lifes such as breath, motion carry out the logic point of systematic science Analysis, accurately obtains multiple hemodynamic parameters, compared to existing continuous non-invasive hemodynamic monitoring more safety and comfort, resists Interference performance is stronger.
Further, the electrocardiosignal and the photoplethysmographic signal are pre-processed, including:
The baseline and noise of the electrocardiosignal are removed, it is average smooth to removing the electrocardiosignal progress after baseline and noise Processing, to removing baseline and noise, carrying out the electrocardiosignal after average smooth treatment, carries out first-order difference processing, to described one Electrocardiosignal after order difference processing does average smooth treatment, obtains the smooth signal of first-order difference of the electrocardiosignal;
Average smooth treatment is carried out to the photoplethysmographic signal, to the photoelectricity volume arteries and veins after average smooth treatment Ripple signal of fighting carries out first-order difference processing, and the photoplethysmographic signal after average smooth treatment and first-order difference processing is entered The average smooth treatment of row, obtains the smooth signal of first-order difference of the photoplethysmographic, by the photoplethysmographic The smooth signal of first-order difference be moved to the left preset time, the preset time passes for the pulse wave of the targeted customer from heart The time cast to needed for the photoplethysmographic signal collection point.
It should be noted that the process to the removal noise of the electrocardiosignal of targeted customer is individual including removal baseline drift, M The average smooth treatment of point, the electrocardiosignal after processing is x (i), and i is sample data sequence number.
In actual applications, can use has limit for length's unit impulse response (Finite Impulse Response, FIR) filter Ripple device removes baseline drift, it would however also be possible to employ other modes remove baseline drift, the invention is not limited in this regard.
It should be noted that the number M of the point of average smooth treatment, need to integrate R wave widths TRWith sample rate fsTo determine. Usual M point at least needs to cover TR, i.e.,
M≥fs·TR
For example, fs=500Hz, TR=0.030 second, then M >=15.
First-order difference processing and average smooth treatment are carried out to removing the electrocardiosignal after noise and the smooth treatment that is averaged, Obtain the smooth signal of first-order difference of the electrocardiosignal;First-order difference is asked to x (i):
D (i)=x (i+1)-x (i)
N number of point is done to d (i) again average smooth, obtain the smooth signal y (i) of first-order difference of electrocardiosignal.
It should be noted that pulse wave is generated in heart, by preset time TpTravel to photoplethysmographic collection Point.When the embodiment of the present invention by the photoplethysmographic after first-order difference and the smooth treatment that is averaged by being moved to the left default Between, so as to combine electrocardiosignal and photoplethysmographic signal acquisition hemodynamic parameter exactly.
Further, the characteristic point of the pretreated photoplethysmographic signal of mark includes:
Photoplethysmographic maximum slope point.
Further, according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Multiple hemodynamic parameters are obtained, including:
According to the R ripples crest location of the electrocardiosignal and the position acquisition of the photoplethysmographic maximum slope point Blood systolic pressure and blood diastolic pressure;
The stroke output and external perihaemal canal resistance of heart are obtained according to the blood systolic pressure and the blood diastolic pressure;
Cardiac output is obtained according to the stroke output of the heart;
The left heart acting is obtained according to the blood systolic pressure, the blood diastolic pressure and the cardiac output;
Arterial compliance is obtained according to the blood systolic pressure, the blood diastolic pressure and the stroke output.
Specifically, blood systolic pressure and blood diastolic pressure are obtained according to equation below:
SBP=a1+a2·PWV
DBP=a3+a4·PWV
Wherein, SBP represents blood systolic pressure, and DBP represents blood diastolic pressure;PWV represents pulse wave velocity, according to such as Lower formula is obtained:
Height is the height of the targeted customer;PWTT represents the R ripples crest and the photocapacitance of the electrocardiosignal Time interval between product pulse wave maximum slope point;a1、a2、a3And a4It is the constant related to the targeted customer, it is necessary to school It is accurate;The unit of blood pressure is mmHg.
The blood oxygen saturation SpO of blood2Calculated with following equation:
Wherein, SpO2Represent blood oxygen saturation;ACFeux rouges、DCFeux rouges、ACInfrared light、DCInfrared lightIt is feux rouges and infrared in emergent light respectively The amplitude of the AC and DC part of light.
The stroke output of heart is obtained according to equation below:
Wherein, SV represents the stroke output of heart;Characteristic quantity K=1/3, HR represent heart rate;a5And a6It is and the target The related constant of user;Usual a5=14.390, a6When=0.0014, more accurate stroke output can be obtained.
Heart rate HR is the number of times of heartbeat per minute, is calculated according to equation below:
HR=60/T
Wherein, T represents the time interval of adjacent R crest peak position in electrocardiosignal.
External perihaemal canal resistance is obtained according to equation below:
SVR=a10+a11·(SBP-DBP)+a12·HR+a13·Age+a14·K
Wherein, SVR represents external perihaemal canal resistance;Age represents the age of the targeted customer;a10、a11、a12、a13And a14 It is the constant related to the targeted customer;K is constant, generally takes 1/3;Generally take a10=0.0034, a11=-0.015, a12 =0.007, a13=-0.003, a14When=4.244, preferable measurement result can be obtained;External perihaemal canal resistance SVR unit It is dynesscm-5
Body surface area is calculated according to equation below
Wherein, a7、a8、a9It is the constant coefficient relevant with crowd or individual, it is necessary to calibrate.Usual a7=0.007184, a8 =0.725, a9Result is more accurate when=0.425;Wherein, height, body weight, BSA unit are respectively:Centimetre, kilogram, square Rice.
External perihaemal canal drag index SI is obtained according to equation below
SI=SVRBSA
Wherein, SI represents external perihaemal canal drag index, and unit is dynesscm-5/m2
Cardiac output is obtained according to equation below:
Wherein, CO represents cardiac output;Cardiac output CO is due to heart contraction, per minute that the blood come is pumped out from left ventricle Liquid total amount, unit is to rise.
Stroke volume index SVI and cardiac output index CI is obtained according to equation below:
SVI=SV/BSA
CI=CO/BSA
Left heart acting is obtained according to equation below:
LCW=a15·(MAP-a16)·CO
Wherein, LCW represents left heart acting;a15、a16It is the constant related to the targeted customer;MAP represents average blood Pressure, MAP=(2DSP+SBP)/3;Left heart acting LCW unit is kgm.
The left heart is obtained according to equation below and is work index LCWI
LCWI=LCW/BSA
Wherein, the unit that the left heart is work index LCWI is kgm/m2
Arterial compliance is obtained according to equation below:
C=SV (SBP-DBP)
Wherein, C represents arterial compliance, and unit is cm5/dynes。
Damping time constant RC is obtained according to equation below
RC=SVRC/ (SBP-DBP)
Wherein, damping time constant RC unit is s.
Further, the characteristic point of the pretreated photoplethysmographic signal of mark also includes:
Point that slope before sequences of ventricular depolarization point, photoplethysmographic maximum slope point is zero, photoplethysmographic Trailing edge slope smallest point.
Further, according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Multiple hemodynamic parameters are obtained, including:
The sequences of ventricular depolarization point and photoplethysmographic slope determined according to the original position of the Q ripples of the electrocardiosignal of detection Time interval between the point that slope before maximum point is zero determines pre-ejection;
According to the slope before the photoplethysmographic maximum slope point be zero point and photoplethysmographic under Drop determines left heart ejection time along the time interval between slope smallest point.
It should be noted that early stage ejection time PEP is since being electrically excited occurs in ventricle, until aorta petal opens it Between time interval;Left heart ejection time LVET be from aorta petal be opened to aortic valve closing between time interval.
As shown in Fig. 2 in Q ripple original positions, sequences of ventricular depolarization starts.Slope before pulse wave gradient maxima is zero B points, aorta petal and pulmonary valve are opened, in the X points of pulse wave trailing edge slope minimum value, aortic valve closing.By pre- If time TpPulse wave has traveled to photoplethysmographic collection point.
Thus, the B points that the slope before pre-ejection PEP is from Q ripples original position to pulse wave gradient maxima is zero it Between time interval;Left heart ejection time LVET is minimum from B points on photoplethysmographic figure to pulse wave trailing edge slope Time interval between the X points of value.
In actual applications,, can be by T with the method for statistics to certain group of people or individualpAs undetermined constant, by with Lower formula is calculated and obtained:
Tp=a17+a18·PWTT
PWTT represents the time between the R ripples crest of the electrocardiosignal and the photoplethysmographic maximum slope point Interval.After early stage ejection time PEP and left heart ejection time LVET is obtained, systolic time ratio is obtained according to equation below STR:
STR=PEP/LVET.
Rich row variability SVV is that the maximum in passing 30 stroke output SV subtracts minimum value, then divided by average value:
The embodiment of the present invention is moved according to ecg wave form, photoplethysmographic and its difference waveform by the blood of innovation Mechanics algorithm, accurately obtains multiple hemodynamic parameters, including stroke output, cardiac output, external perihaemal canal resistance, cardiac muscle Contraction of indices etc..
In actual applications, the hemodynamic parameter acquisition methods of the embodiment of the present invention can be applied to as shown in Figure 3 In many life sign monitor systems of sports type continuous non-invasive.The many life sign monitor systems of sports type continuous non-invasive, are after application The various invasive and noninvasive life sign monitor system used when operating room, hospital, family are being monitored static person and rest Afterwards, it is of future generation to carry out health and the core technology of physiological parameter monitoring in long-range tranquillization and motion.
As shown in figure 3, many life sign monitor systems of sports type continuous non-invasive include monitoring host machine part 31 and wearing portion Divide 32;Monitor host machine part 31 include the peripheral data such as system host, keyboard, mouse, display screen, printer, network, figure and The input such as energy, output and communication apparatus;Wearing portion 32 includes electrocardiogram acquisition electrode 321, photoplethysmographic and gathered Device 322 and monitor 323.Monitor 323 gathers electrocardiogram acquisition electrode 321 and photoplethysmographic collector 322 Electrocardiosignal and pulse signal are sent to monitoring host machine part 31;Monitor the method pair according to embodiments of the present invention of host machine part 31 Electrocardiosignal and photoplethysmographic signal are handled;And by the electrocardiosignal after processing, photoplethysmographic signal, Hemodynamic parameter and relevant information are sent to monitor 323;Electrocardiosignal light after the processing that 323 pairs of monitor is received Power Capacity pulse wave signal and monitoring relevant parameter are shown with data.In actual applications, monitoring host machine part 31 and prison Surveying instrument 323 can be communicated by USB serial ports, bluetooth, netting twine and WiFi, can also shared processor, memory and bus.Should With the software of the hemodynamic parameter acquisition methods of the invention provided, the PC for meeting hardware configuration requirement can also be loaded into Computer, notebook computer, server, single-chip microcomputer, iPad, tablet personal computer, on the smart mobile phone of iOS and Android operation system, make The said equipment turns into many life sign monitor system host machine parts of sports type continuous non-invasive.
Fig. 4 is the structural representation of the hemodynamic parameter acquisition device of one embodiment of the invention.As shown in figure 4, The device of the embodiment of the present invention includes signal gathering unit 41, pretreatment unit 42, processing unit 43 and hemodynamic parameter Acquiring unit 44, specifically:
Signal gathering unit 41, electrocardiosignal and photoplethysmographic signal for gathering targeted customer;
Pretreatment unit 42, for being pre-processed to the electrocardiosignal and the photoplethysmographic signal;
Processing unit 43, for detecting pretreated electrocardiosignal Q ripples and R ripples, and to pretreated photoelectricity volume Pulse wave signal carries out feature point identification;
Hemodynamic parameter acquiring unit 44, for according to the pretreated electrocardiosignal and the photoelectricity volume The characteristic point of pulse wave signal obtains multiple hemodynamic parameters.
The hemodynamic parameter acquisition device of the embodiment of the present invention can be used for performing above method embodiment, its principle Similar with technique effect, here is omitted.
Fig. 5 is the structural representation of the electronic equipment of one embodiment of the invention.
Reference picture 5, electronic equipment includes:Processor (processor) 51, memory (memory) 52 and bus (bus) 53;Wherein,
Processor 51 and memory 52 complete mutual communication by bus 53;
Processor 51 is used to call the programmed instruction in memory 52, to perform the blood that above-mentioned each method embodiment is provided Hemodynamics parameter acquiring method.
In addition, the logical order in above-mentioned memory 52 can be realized by the form of SFU software functional unit and is used as solely Vertical production marketing in use, can be stored in a computer read/write memory medium.Understood based on such, this hair The part or the part of the technical scheme that bright technical scheme substantially contributes to prior art in other words can be with soft The form of part product is embodied, and the computer software product is stored in a storage medium, including some instructions are to make It (can be personal computer, server, smart mobile phone, single-chip microcomputer, tablet personal computer or the network equipment to obtain a computer equipment Deng) perform all or part of step of each of the invention embodiment methods described.And foregoing storage medium includes:USB flash disk, movement Hard disk, hard disk, read-only storage (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
The present embodiment provides a kind of computer program product, and the computer program product includes being stored in non-transient calculating Computer program on machine readable storage medium storing program for executing, the computer program includes programmed instruction, when described program instruction is calculated When machine is performed, computer is able to carry out the hemodynamic parameter acquisition methods that above-mentioned each method embodiment is provided.
The present embodiment provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage medium storing program for executing Computer instruction is stored, the computer instruction makes the computer perform the hemodynamics that above-mentioned each method embodiment is provided Learn parameter acquiring method.
Hemodynamic parameter acquisition methods and device provided in an embodiment of the present invention, belong to noninvasive hemodynamics prison Survey, propose photoplethysmographic signal include VPV and volume change information, by targeted customer tranquillization, fortune The dynamic logic analysis for waiting the electrocardiosignal and photoplethysmographic signal of the generation in daily life to carry out systematic science, accurately Multiple hemodynamic parameters are obtained, compared to existing noninvasive hemodynamics monitoring more safety and comfort, antijamming capability is more By force.
It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Apply the form of example.Moreover, the present invention can be used in one or more computers for wherein including computer usable program code The computer program production that usable storage medium is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.) The form of product.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product Figure and/or block diagram are described.Being interpreted as can be by each in computer program instructions implementation process figure and/or block diagram Flow and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer journeys can be provided Sequence instruction to all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices processor with Produce a machine so that being produced by the instruction of computer or the computing device of other programmable data processing devices is used for Realize the dress for the function of being specified in one flow of flow chart or multiple flows and/or one square frame of block diagram or multiple square frames Put.
It should be noted that term " comprising ", "comprising" or its any other variant are intended to the bag of nonexcludability Contain, so that process, method, article or equipment including a series of key elements are not only including those key elements, but also including Other key elements being not expressly set out, or also include for this process, method, article or the intrinsic key element of equipment. In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that including the key element Process, method, article or equipment in also there is other identical element.
In the specification of the present invention, numerous specific details are set forth.Although it is understood that, embodiments of the invention can To be put into practice in the case of these no details.In some instances, known method, structure and skill is not been shown in detail Art, so as not to obscure the understanding of this description.Similarly, it will be appreciated that disclose in order to simplify the present invention and helps to understand respectively One or more of individual inventive aspect, above in the description of the exemplary embodiment of the present invention, each of the invention is special Levy and be grouped together into sometimes in single embodiment, figure or descriptions thereof.However, should not be by the method solution of the disclosure It is interpreted into the following intention of reflection:I.e. the present invention for required protection requirement is than the feature that is expressly recited in each claim more Many features.More precisely, as the following claims reflect, inventive aspect is to be less than single reality disclosed above Apply all features of example.Therefore, it then follows thus claims of embodiment are expressly incorporated in the embodiment, Wherein each claim is in itself as the separate embodiments of the present invention.
Above example is merely to illustrate technical scheme, rather than its limitations;Although with reference to the foregoing embodiments The present invention is described in detail, it will be understood by those within the art that:It still can be to foregoing each implementation Technical scheme described in example is modified, or carries out equivalent substitution to which part technical characteristic;And these are changed or replaced Change, the essence of appropriate technical solution is departed from the spirit and scope of various embodiments of the present invention technical scheme.

Claims (10)

1. a kind of hemodynamic parameter acquisition methods, it is characterised in that including:
Gather the electrocardiosignal and photoplethysmographic signal of targeted customer;
The electrocardiosignal and the photoplethysmographic signal are pre-processed;
Pretreated electrocardiosignal Q ripples and R ripples are detected, and characteristic point is carried out to pretreated photoplethysmographic signal Mark;
Multiple blood are obtained according to the characteristic point of the electrocardiosignal Q ripples of detection and R ripples and the photoplethysmographic signal Kinetic parameter.
2. according to the method described in claim 1, it is characterised in that the electrocardiosignal and the photoplethysmographic are believed Number pre-processed, including:
The baseline and noise of the electrocardiosignal are removed, average smooth place is carried out to removing the electrocardiosignal after baseline and noise Reason, to removing baseline and noise, carrying out the electrocardiosignal after average smooth treatment, carries out first-order difference processing, to first-order difference Electrocardiosignal after processing does average smooth treatment, obtains the smooth signal of first-order difference of the electrocardiosignal;
Average smooth treatment is carried out to the photoplethysmographic signal, to the photoplethysmographic after average smooth treatment Signal carries out first-order difference processing, and the photoplethysmographic signal after average smooth treatment and first-order difference processing is put down Equal smooth treatment, obtains the smooth signal of first-order difference of the photoplethysmographic signal;By the photoplethysmographic The smooth signal of first-order difference of signal is moved to the left preset time, the preset time for the targeted customer pulse wave from the heart The dirty time propagated to needed for the photoplethysmographic signal collection point.
3. according to the method described in claim 1, it is characterised in that the pretreated photoplethysmographic signal of mark Characteristic point includes:
Photoplethysmographic maximum slope point.
4. method according to claim 3, it is characterised in that according to the electrocardiosignal Q ripples and R ripples of detection and the light The characteristic point of Power Capacity pulse wave signal obtains multiple hemodynamic parameters, including:
According to the R ripples crest location of the electrocardiosignal and the position acquisition blood of the photoplethysmographic maximum slope point Systolic pressure and blood diastolic pressure;
The stroke output and external perihaemal canal resistance of heart are obtained according to the blood systolic pressure and the blood diastolic pressure;
Cardiac output is obtained according to the stroke output of the heart;
The left heart acting is obtained according to the blood systolic pressure, the blood diastolic pressure and the cardiac output;
Arterial compliance is obtained according to the blood systolic pressure, the blood diastolic pressure and the stroke output.
5. method according to claim 4, it is characterised in that
Blood systolic pressure and blood diastolic pressure are obtained according to equation below:
SBP=a1+a2·PWV
DBP=a3+a4·PWV
Wherein, SBP represents blood systolic pressure, and DBP represents blood diastolic pressure;PWV represents pulse wave velocity, according to following public affairs Formula is obtained:
<mrow> <mi>P</mi> <mi>W</mi> <mi>V</mi> <mo>=</mo> <mfrac> <mrow> <mi>H</mi> <mi>e</mi> <mi>i</mi> <mi>g</mi> <mi>h</mi> <mi>t</mi> </mrow> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>P</mi> <mi>W</mi> <mi>T</mi> <mi>T</mi> </mrow> </mfrac> </mrow>
Height is the height of the targeted customer;PWTT represents the R ripples crest and the photoelectricity volume arteries and veins of the electrocardiosignal The time interval fought between ripple maximum slope point;a1、a2、a3And a4It is the constant related to the targeted customer;According to following public affairs Formula obtains the stroke output of heart:
<mrow> <mi>S</mi> <mi>V</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>a</mi> <mn>5</mn> </msub> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>S</mi> <mi>B</mi> <mi>P</mi> <mo>-</mo> <mi>D</mi> <mi>B</mi> <mi>P</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <mi>K</mi> <mn>2</mn> </msup> <mo>&amp;CenterDot;</mo> <mi>H</mi> <mi>R</mi> </mrow> </mfrac> <mo>+</mo> <msub> <mi>a</mi> <mn>6</mn> </msub> </mrow>
Wherein, SV represents the stroke output of heart;Characteristic quantity K=1/3, HR represent heart rate;a5And a6It is and the targeted customer Related constant;
External perihaemal canal resistance is obtained according to equation below:
SVR=a10+a11·(SBP-DBP)+a12·HR+a13·Age+a14·K
Wherein, SVR represents external perihaemal canal resistance;Age represents the age of the targeted customer;a10、a11、a12、a13And a14It is and institute State the related constant of targeted customer;Cardiac output is obtained according to equation below:
<mrow> <mi>C</mi> <mi>O</mi> <mo>=</mo> <mfrac> <mrow> <mi>S</mi> <mi>V</mi> <mo>&amp;CenterDot;</mo> <mi>H</mi> <mi>R</mi> </mrow> <mn>1000</mn> </mfrac> </mrow>
Wherein, CO represents cardiac output;
Left heart acting is obtained according to equation below:
LCW=a15·(MAP-a16)·CO
Wherein, LCW represents left heart acting;a15、a16It is the constant related to the targeted customer;MAP represents mean blood pressure, MAP =(2DSP+SBP)/3;
Arterial compliance is obtained according to equation below:
C=SV (SBP-DBP)
Wherein, C represents arterial compliance.
6. method according to claim 3, it is characterised in that the pretreated photoplethysmographic signal of mark Characteristic point also includes:
Point that slope before sequences of ventricular depolarization point, photoplethysmographic maximum slope point is zero, photoplethysmographic decline Along slope smallest point.
7. method according to claim 6, it is characterised in that according to the electrocardiosignal Q ripples and R ripples of detection and the light The characteristic point of Power Capacity pulse wave signal obtains multiple hemodynamic parameters, including:
The sequences of ventricular depolarization point and photoplethysmographic maximum slope determined according to the original position of the Q ripples of the electrocardiosignal of detection Time interval between the point that slope before point is zero determines pre-ejection;
The point and photoplethysmographic trailing edge for being zero according to the slope before the photoplethysmographic maximum slope point Time interval between slope smallest point determines left heart ejection time.
8. a kind of hemodynamic parameter acquisition device, it is characterised in that including:
Signal gathering unit, electrocardiosignal and photoplethysmographic signal for gathering targeted customer;
Pretreatment unit, for being pre-processed to the electrocardiosignal and the photoplethysmographic signal;
Processing unit, for detecting pretreated electrocardiosignal Q ripples and R ripples, and to pretreated photoplethysmographic Signal carries out feature point identification;
Hemodynamic parameter acquiring unit, for according to the pretreated electrocardiosignal and the photoplethysmographic The characteristic point of signal obtains multiple hemodynamic parameters.
9. a kind of electronic equipment, it is characterised in that including:Processor, memory and bus;Wherein,
Processor and memory complete mutual communication by bus;
Processor is used to call the programmed instruction in memory, and the Hemodynamics described in 1~7 any one is required with perform claim Parameter acquiring method.
10. a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage medium storing program for executing storage computer refers to Order, the computer instruction makes the computer perform claim require the hemodynamic parameter acquisition side described in 1~7 any one Method.
CN201710398262.1A 2017-05-31 2017-05-31 A kind of hemodynamic parameter acquisition methods and device Pending CN107296593A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710398262.1A CN107296593A (en) 2017-05-31 2017-05-31 A kind of hemodynamic parameter acquisition methods and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710398262.1A CN107296593A (en) 2017-05-31 2017-05-31 A kind of hemodynamic parameter acquisition methods and device

Publications (1)

Publication Number Publication Date
CN107296593A true CN107296593A (en) 2017-10-27

Family

ID=60137355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710398262.1A Pending CN107296593A (en) 2017-05-31 2017-05-31 A kind of hemodynamic parameter acquisition methods and device

Country Status (1)

Country Link
CN (1) CN107296593A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107854123A (en) * 2017-12-13 2018-03-30 清华大学深圳研究生院 One kind is without oversleeve continuous blood pressure monitoring method and apparatus
CN109171677A (en) * 2018-09-20 2019-01-11 芯海科技(深圳)股份有限公司 A kind of Portable type measurement unit and the hemodynamic parameter measurement method based on the device
CN110811592A (en) * 2018-08-10 2020-02-21 株式会社东芝 Blood circulation detection device and blood circulation detection method
CN110974172A (en) * 2019-12-13 2020-04-10 北京理工大学 Real-time physiological parameter measuring system
CN111214218A (en) * 2020-01-13 2020-06-02 京东方科技集团股份有限公司 Detection equipment for multiple physiological parameters
CN111248894A (en) * 2018-12-02 2020-06-09 西诺嘉医药有限公司 Apparatus and method for recovering pulse waveform and other hemodynamics using ECG
CN113288103A (en) * 2021-04-01 2021-08-24 安徽通灵仿生科技有限公司 Non-invasive cardiac output monitoring system and method
CN113616167A (en) * 2021-07-21 2021-11-09 深圳市千帆电子有限公司 Hemodynamic monitoring method and system
ES2875673A1 (en) * 2020-05-07 2021-11-10 Rgb Medical Devices Sa Method, device and system to determine the trend of cardiac output (Machine-translation by Google Translate, not legally binding)
CN113633290A (en) * 2021-08-18 2021-11-12 北京航空航天大学 Method, apparatus, electronic device and readable storage medium for blood perfusion monitoring

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102159132A (en) * 2008-09-22 2011-08-17 奇塔医疗公司 System and method for determining blood flow
CN103796576A (en) * 2011-07-25 2014-05-14 奇塔医疗公司 Method and system for monitoring hemodynamics

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102159132A (en) * 2008-09-22 2011-08-17 奇塔医疗公司 System and method for determining blood flow
CN103796576A (en) * 2011-07-25 2014-05-14 奇塔医疗公司 Method and system for monitoring hemodynamics

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
姜春华: "《活血化瘀研究》", 31 December 1981 *
成志霞: "基于多生理参数的动脉硬化检测技术研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 *
易继锴,侯媛彬: "《智能控制技术》", 31 July 2007 *
王灿: "无创连续心输出量分析算法研究与软件设计", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
董骁: "可穿戴式多生理参数监护系统的研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 *
齐保申等: "1695名正常人群外周血管阻力异常及其相关因素分析", 《中国预防医学杂志》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107854123A (en) * 2017-12-13 2018-03-30 清华大学深圳研究生院 One kind is without oversleeve continuous blood pressure monitoring method and apparatus
CN110811592A (en) * 2018-08-10 2020-02-21 株式会社东芝 Blood circulation detection device and blood circulation detection method
CN109171677A (en) * 2018-09-20 2019-01-11 芯海科技(深圳)股份有限公司 A kind of Portable type measurement unit and the hemodynamic parameter measurement method based on the device
CN111248894B (en) * 2018-12-02 2023-05-26 西诺嘉医药有限公司 Apparatus and method for recovering pulse waveforms and other hemodynamics using ECG
CN111248894A (en) * 2018-12-02 2020-06-09 西诺嘉医药有限公司 Apparatus and method for recovering pulse waveform and other hemodynamics using ECG
CN110974172A (en) * 2019-12-13 2020-04-10 北京理工大学 Real-time physiological parameter measuring system
CN111214218A (en) * 2020-01-13 2020-06-02 京东方科技集团股份有限公司 Detection equipment for multiple physiological parameters
CN111214218B (en) * 2020-01-13 2024-02-09 京东方科技集团股份有限公司 Multi-physiological parameter detection equipment
ES2875673A1 (en) * 2020-05-07 2021-11-10 Rgb Medical Devices Sa Method, device and system to determine the trend of cardiac output (Machine-translation by Google Translate, not legally binding)
CN113288103A (en) * 2021-04-01 2021-08-24 安徽通灵仿生科技有限公司 Non-invasive cardiac output monitoring system and method
CN113616167A (en) * 2021-07-21 2021-11-09 深圳市千帆电子有限公司 Hemodynamic monitoring method and system
CN113633290B (en) * 2021-08-18 2022-04-01 北京航空航天大学 Method, apparatus, electronic device and readable storage medium for blood perfusion monitoring
CN113633290A (en) * 2021-08-18 2021-11-12 北京航空航天大学 Method, apparatus, electronic device and readable storage medium for blood perfusion monitoring

Similar Documents

Publication Publication Date Title
CN107296593A (en) A kind of hemodynamic parameter acquisition methods and device
Hernandez et al. Biowatch: estimation of heart and breathing rates from wrist motions
Guidoboni et al. Cardiovascular function and ballistocardiogram: A relationship interpreted via mathematical modeling
Siddiqui et al. A pulse rate estimation algorithm using PPG and smartphone camera
Couceiro et al. Detection of motion artifact patterns in photoplethysmographic signals based on time and period domain analysis
Silva et al. Signal quality estimation with multichannel adaptive filtering in intensive care settings
Visvanathan et al. Estimation of blood pressure levels from reflective photoplethysmograph using smart phones
Lim et al. Adaptive template matching of photoplethysmogram pulses to detect motion artefact
Ibrahim et al. Continuous blood pressure monitoring using wrist-worn bio-impedance sensors with wet electrodes
Ghamari et al. Design and prototyping of a wristband-type wireless photoplethysmographic device for heart rate variability signal analysis
Marques et al. A real time, wearable ECG and blood pressure monitoring system
Kanawade et al. Photoplethysmography based arrhythmia detection and classification
Zhang et al. Study of cuffless blood pressure estimation method based on multiple physiological parameters
Cao et al. Crisp-BP: Continuous wrist PPG-based blood pressure measurement
US20220031220A1 (en) Model-based sensor technology for detection of cardiovascular status
Şentürk et al. Towards wearable blood pressure measurement systems from biosignals: a review
Wang et al. Ballistocardiogram heart rate detection: Improved methodology based on a three-layer filter
Ahmed et al. Feasibility analysis for estimation of blood pressure and heart rate using a smart eye wear
Athaya et al. An efficient fingertip photoplethysmographic signal artifact detection method: A machine learning approach
Krizea et al. Accurate detection of heart rate and blood oxygen saturation in reflective photoplethysmography
Szakacs-Simon et al. Pulse oximeter based monitoring system for people at risk
Janjua et al. A low-cost tonometer alternative: A comparison between photoplethysmogram and finger ballistocardiogram and validation against tonometric waveform
Klum et al. Impedance plethysmography for respiratory flow and rate estimation using multilayer perceptrons
Zainal et al. Real Time PPG Data Acquisition with GUI based Application for HRV Measurement
Villata Cuffless blood pressure measurement

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20171027

RJ01 Rejection of invention patent application after publication