CN106413534A - Blood-pressure continuous-measurement device, measurement model establishment method, and system - Google Patents

Blood-pressure continuous-measurement device, measurement model establishment method, and system Download PDF

Info

Publication number
CN106413534A
CN106413534A CN201580000417.0A CN201580000417A CN106413534A CN 106413534 A CN106413534 A CN 106413534A CN 201580000417 A CN201580000417 A CN 201580000417A CN 106413534 A CN106413534 A CN 106413534A
Authority
CN
China
Prior art keywords
pulse wave
characteristic parameter
blood pressure
feature point
pressure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201580000417.0A
Other languages
Chinese (zh)
Other versions
CN106413534B (en
Inventor
李烨
沈子骁
苗芬
何青云
尹丽妍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Publication of CN106413534A publication Critical patent/CN106413534A/en
Application granted granted Critical
Publication of CN106413534B publication Critical patent/CN106413534B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/021Measuring pressure in heart or blood vessels
    • 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

Abstract

The invention provides a blood-pressure continuous-measurement device, a measurement model establishment method, and a system. The device comprises: a signal collector (100), used for collecting electrocardiosignals and pulse wave signals of a measured object for multiple consecutive cardiac cycles, so as to obtain two sets of synchronous electrocardiosignals and pulse wave signals; a signal reception module (210), used for receiving the two sets of synchronous electrocardiosignals and pulse wave signals; and a blood pressure computing module (220), used for extracting a first characteristic point of the pulse wave signals according to the two sets of synchronous electrocardiosignals and pulse wave signals, extracting a second characteristic point of the electrocardiosignals in a same cardiac cycle according to the first characteristic point, acquiring corresponding characteristic values according to the first characteristic point and the second characteristic point, selecting a characteristic parameter from the characteristic values according to an evaluation parameter, and acquiring a blood pressure value according to a pre-established blood pressure computing model of characteristic parameters. The device, the method and the system have the characteristics of high precision, low complexity, low cost, convenience in measurement, and high real-time quality.

Description

Continuous blood pressure measurer, measurement model method for building up and system
Technical field
The present invention relates to blood pressure measurement technology, more particularly to a kind of continuous blood pressure measurer, measurement model foundation side Method and system.
Background technology
Compared to traditional discontinuous blood pressure measuring method, continuous BP measurement method, can in each cardiac cycle Provide pressure value by shooting.This, for understanding blood pressure rule, prevents cardiovascular and cerebrovascular disease, has great importance.Meanwhile, Due to continuous BP measurement method can real-time tracing by the pressure value of person, not only can reduce " white coat hypertension " and " hidden Property hypertension " hidden danger, contribute to the daily monitoring of hyperpietic additionally it is possible to the indirect reflection effect to human body for the antihypertensive drugs Really, to avoid the abuse of medicine.In addition, continuous BP measurement, additionally aid the analysis of blood pressure variability, this is for further Solve tested health, the occurrence risk of assessment cardiovascular and cerebrovascular disease, there is important reference value.
However, current noninvasive continuous BP measurement equipment, it is former based on volume-compensation method and angiosthenia method mostly Reason design.Measuring apparatus using angiosthenia method can reach higher precision, substantially can reach long period noninvasive company The requirement of continuous blood pressure measurement.But due to the high sensitivity to displacement for the sensor, long-time sensor measurement position to be kept It is relatively fixed relatively difficult, air bag pressue device also affects measured's comfort level, therefore base in long-time measurement process simultaneously Blood pressure measurement device in angiosthenia still has many problems demand to change in motion measurement, long-term measurement and simplification operating aspect Enter.Measuring apparatus using volume-compensation method can continuously measure blood pressure of often fighting, and can measure blood pressure waveform without distortion.But Due to the effect of gasbag pressure, long-time measurement can lead to venous congestion to affect certainty of measurement, can bring to measured simultaneously Uncomfortable.In addition, the noninvasive continuous BP measurement equipment of principle design based on volume-compensation method and angiosthenia method is heavy, price Expensive, carry and gather equal inconvenience, be unfavorable for the daily physiology monitoring of family.
Content of the invention
Based on this it is necessary to be directed to deficiency of the prior art, a kind of continuous blood pressure measurer, measurement model is provided to build Cube method and system.
A kind of continuous blood pressure measurer, described device includes:
Signal picker, including Electrocardial signal acquisition device and pulse wave signal collector, for synchronous acquisition measurement object Multiple continuous cardiac cycle in electrocardiosignal and pulse wave signal, obtain two groups of synchronous electrocardiosignals and pulse wave letter Number;
Signal receiving module, for receiving two groups of synchronous electrocardiosignals and pulse wave signal;
Blood pressure computing module, for the electrocardiosignal synchronous according to two groups and pulse wave signal, extracts described pulse wave letter Number fisrt feature point, and according to described fisrt feature point obtain same cardiac cycle in electrocardiosignal second feature point, Corresponding characteristic value is obtained according to described fisrt feature point and described second feature point, using assessment parameter from described characteristic value Selected characteristic parameter, according to the blood pressure computation model of the described characteristic parameter pre-building, obtains pressure value.
Wherein in some embodiments, described blood pressure computing module includes:
Characteristic value acquiring unit, for the electrocardiosignal synchronous according to two groups and pulse wave signal, extracts described pulse wave The fisrt feature point of signal, and the second feature of the electrocardiosignal in same cardiac cycle is obtained according to described fisrt feature point Point;
Characteristic parameter acquiring unit, special accordingly for according to described fisrt feature point and described second feature point, obtaining Levy parameter;With
Blood pressure calculation unit, for the blood pressure computation model according to the described characteristic parameter pre-building, obtains pressure value.
Wherein in some embodiments, described blood pressure computing module also includes:
Filter unit, for being filtered processing to described two groups synchronous electrocardiosignals and pulse wave signal respectively.
Wherein in some embodiments, described characteristic parameter includes:Pulse wave propagation time, cardiac cycle, pulse wave are received Rise time contracting phase and pulse wave K value.
Wherein in some embodiments, the blood pressure computation model of described characteristic parameter includes multivariate linear model, described many First linear model is:
SBP=P D G,
MAP=P E G,
DBP=P F G,
Wherein, SBP is systolic pressure, and MAP is mean blood pressure, and DBP is diastolic pressure, P=[1 PWTT EcgPeriod PpgSystolicTime ppg_K], PWTT is described pulse wave propagation time, and EcgPeriod is described cardiac cycle, PpgSystolicTime is the described rise time in pulse wave systole phase, and ppg_K is described pulse wave K value, and D, E and F are respectively institute State systolic pressure, described diastolic pressure and described mean blood pressure corresponding bp coefficient matrix, G is the physiological parameter square of measurement object Battle array,Gender is sex, and age is the age, and height is height, and weight is body weight, and BMI is physique ArmLength is brachium to ratio.
Wherein in some embodiments, the blood pressure computation model of described characteristic parameter includes nonlinear model, described non-thread Property model includes the nonlinear regression model (NLRM) based on SVMs, the described blood pressure computation model setting up described characteristic parameter Process includes:
Obtain the training set between described characteristic parameter and described systolic pressure;
Construct decision function using the described nonlinear regression model (NLRM) based on SVMs and calculate optimum solution;
Obtain systolic pressure and the diastolic pressure coefficient matrix based on described data set according to described decision function;
Final pressure value is obtained according to described systolic pressure and diastolic pressure coefficient matrix.
Wherein in some embodiments, described characteristic parameter is chosen unit and is additionally operable to, and is initialized special according to described characteristic value Levy parameter set and characteristic parameter subset;Choose multiple candidate feature subset;Special to the plurality of candidate using described valuation functions Levy subset to be estimated, obtain the characteristic parameter conforming to a predetermined condition.
Wherein in some embodiments, described device also includes:
Memory module, for storing described characteristic value, described characteristic parameter and its corresponding pressure value;
Display module, for showing described characteristic value, described characteristic parameter and its corresponding pressure value.
Wherein in some embodiments, described device also includes:
Communication module, for receiving two groups of synchronous electrocardiosignals and the pulse wave signal of described signal picker collection, And the various data is activations storing memory module are to external device (ED);Or,
It is respectively provided with wireless signal Transmit-Receive Unit, signal picker and signal in signal picker and signal receiving module Carried out data transmission by wireless signal between receiver module and communicate.
System set up by a kind of continuous BP measurement model, and described system includes:
Signal acquisition module, for receiving electrocardiosignal and pulse wave in two groups of synchronous multiple continuous cardiac cycles Signal;
Feature point extraction module, for extracting the fisrt feature of the described pulse wave signal that described signal acquisition module obtains Point, and the second feature point of the electrocardiosignal in same cardiac cycle is obtained according to described fisrt feature point;
Characteristic value acquisition module, for according to described fisrt feature point and described second feature point, obtaining corresponding feature Parameter;
Model building module, for setting up the blood pressure computation model of described characteristic parameter.
A kind of continuous BP measurement method for establishing model, methods described includes:
Obtain the electrocardiosignal in synchronous multiple continuous cardiac cycle and pulse wave signal;
Extract the fisrt feature point of described pulse wave signal, and obtained according to the fisrt feature point of described pulse wave signal same The second feature point of the electrocardiosignal in one cardiac cycle;
According to described fisrt feature point and described second feature point, obtain corresponding characteristic parameter;
Based on corresponding characteristic parameter, set up the blood pressure computation model of described characteristic parameter.
Wherein in some embodiments, described characteristic parameter includes:Pulse wave propagation time, cardiac cycle, pulse wave are received Rise time contracting phase and pulse wave K value.
Wherein in some embodiments, the blood pressure computation model of described characteristic parameter includes multivariate linear model, described many First linear model is:
SBP=P D G,
MAP=P E G,
DBP=P F G,
Wherein, SBP is described systolic pressure, and MAP is described mean blood pressure, and DBP is described diastolic pressure, P=[1 PWTT EcgPeriod PpgSystolicTime ppg_K], PWTT is described pulse wave propagation time, and EcgPeriod is described aroused in interest In the cycle, PpgSystolicTime is the described rise time in pulse wave systole phase, and ppg_K is described pulse wave K value, and D, E and F divide Not Wei described systolic pressure, described diastolic pressure and described mean blood pressure corresponding bp coefficient matrix, Gender is described sex, and age is the described age, and height is described height, and weight is described body weight, and BMI is described body Matter ratio, armLength is described brachium.
Wherein in some embodiments, the blood pressure computation model of described characteristic parameter includes nonlinear model, described non-thread Property model includes the nonlinear regression model (NLRM) based on SVMs, the non-thread based on SVMs for the described blood pressure calculation unit Property regression model obtain pressure value process include:
Wherein in some embodiments, characteristic parameter obtaining step includes:
According to characteristic value initialization feature parameter set and characteristic parameter subset;
Select evaluation function, and multiple characteristic parameter subsets are generated according to described evaluation function;
Optimal characteristics subset of parameters is chosen from the plurality of characteristic parameter subset according to described evaluation function.
The blood pressure measuring device of above-described embodiment and blood pressure calculating method and system just can be real based on several characteristic parameters Now the unperturbed formula of blood pressure is continuously measured.It has high precisely, low complex degree, low cost, measurement is convenient, real-time is high spy Point.
Brief description
Fig. 1 is the structured flowchart of the continuous blood pressure measurer of one embodiment of the invention;
Fig. 2 is the structured flowchart of the continuous blood pressure measurer of another embodiment of the present invention;
Fig. 3 is the structured flowchart of the signal picker of one embodiment of the invention;
Fig. 4 is the flow chart of the continuous BP measurement method for establishing model of one embodiment of the invention;
Fig. 5 is filtered electrocardiosignal and the pulse wave signal schematic diagram of one embodiment of the invention;
Fig. 6 is the schematic diagram of the fisrt feature point of the pulse wave signal of extraction of one embodiment of the invention;
Fig. 7 is the schematic diagram of the second feature point of the electrocardiosignal of extraction of one embodiment of the invention;
Fig. 8 is the schematic diagram of the second feature point of the electrocardiosignal of extraction of another embodiment of the present invention;
Fig. 9 is the rise time in systole phase of one embodiment of the invention and the schematic diagram of fall time diastole;
Figure 10 is the pulse wave K value schematic diagram of one embodiment of the invention;
Figure 11 is the pulse wave propagation time schematic diagram of one embodiment of the invention;
Figure 12 is that the structured flowchart of system set up by the continuous BP measurement model of one embodiment of the invention.
Specific embodiment
In order that the objects, technical solutions and advantages of the present invention become more apparent, below in conjunction with drawings and Examples, right The present invention is further elaborated.It should be appreciated that specific embodiment described herein is only in order to explain the present invention, and It is not used in the restriction present invention.
In the exemplary embodiment, when blood pressure is that blood flows in intravascular, the pressure of vasoactive wall, it is to promote The power that blood flows in intravascular.Ventricular contraction, blood from ventricle flow into artery, the now pressure highest to artery for the blood, Referred to as systolic pressure (systolic blood pressure, SBP).Ventricular diastole, arteries elastical retraction, blood is still slowly Continuation flow forward, but drop in blood pressure, pressure now is referred to as diastolic pressure (diastolic blood pressure, DBP).The heart Dynamic cycle (cardiac cycle) refers to initiate the initial of heartbeat next time from a heartbeat, and cardiovascular system is experienced Process.During diastole internal pressure reduce, vena cave blood flow back enter the heart, during heart contraction internal pressure raise, by blood pump arrive move Arteries and veins.Heart often shrinks and once constitutes a cardiac cycle with diastole.In one cardiac cycle, the time domain average value of arterial pressure is Mean arterial pressure (mean arterial blood pressure, MAP).
Sphygmomanometry substantially can be divided into the direct method of measurement and the indirect method of measurement two big class.The direct method of measurement is will to connect pressure The conduit of force snesor is percutaneously inserted directly into main artery or heart detection blood pressure signal, can continuously be measured.Due to this Method directly records blood pressure, and data is the most accurate, is therefore regarded as the goldstandard of blood pressure detecting by the world, but its technical requirements is relatively Height, and have necessarily traumatic, so being only applicable to rescue and the major operation patient of urgent patient.The indirect method of measurement is by detection The parameters such as the beating of ductus arteriosus wall, capacity of blood vessel change obtain blood pressure indirectly.Method is simple due to this, therefore in clinic On be used widely.
The indirect method of measurement can be divided into batch (-type) mensuration and continous way mensuration two big class again.With Korotkoff's Sound method and oscillographic method The general principle of the batch (-type) mensuration for representing is with gas sleeve interruption artery blood flow, then by detection during venting In characteristic point moment gas sleeve, pressure is determining pressure value.The shortcoming of this method be all with blood pressure phase of often fighting in measurement process before Carry, record pressure value that blood pressure is a certain particular moment (only not 2 values of SBP and DBP within same cardiac cycle) and differ Surely it is the representational pressure value of measured.Continous way mensuration is noninvasive continuous measurement blood pressure in a certain amount of time, can Detect often fight blood pressure and continuous arterial pressure waveform, be that clinical diagnosis provides more fully foundation with treatment, particularly exist Clinical monitoring and in particular cases observation blood pressure consecutive variations aspect have the incomparable advantage of conventional method.Hereinafter will have Body is described with reference to the accompanying drawings each way of example of the present invention.
Terminology used in the present invention " module ", " unit " are made up of assembly, and " assembly " refers to and present system Related entity, or the combination of software, hardware, hardware and software.For example, assembly can be entering of running on a processor Journey, object, program and computer.As an example, the application program (application, APP) operating in mobile terminal can be Assembly.In addition, assembly can include one or more assemblies.
As shown in figure 1, the continuous blood pressure measurer 10 of the embodiment of the present invention includes:Signal picker 100, control module 200th, memory module 300 and display module 400.
Signal picker 100 collection is related to each item data of the raw body activity of measurement object.Here, it is related to raw body activity Various types of data (biological information) can be related to electrocardio, pulse frequency, cardiac rate, etc. data (information).Signals collecting Device 100 has such as Electrocardial signal acquisition device 110 and pulse wave signal collector 120.Using Electrocardial signal acquisition device 110 from survey The ECG detecting position of amount object obtains electrocardiosignal.Pulse wave signal collector 120 is from the pulse wave test section of measurement object Position obtains pulse wave signal.In some of them embodiment of the present invention, electrocardiosignal includes electrocardiogram, and pulse wave signal includes arteries and veins Fight waveform figure.By with reference to Fig. 3 in Electrocardial signal acquisition device 110 described in detail below and pulse wave signal collector 120 Configuration.
Control module 200 integrally control and measure device 10, and process the measurement for example being gathered by signal picker 100 The electrocardiosignal of object and pulse wave signal.Specifically, the measurement object that control module 200 is gathered based on signal picker 100 Electrocardiosignal and pulse wave signal carry out the blood pressure of computation and measurement object.Additionally, control module 200 is based on acquisition control module 240 carry out control signal collector 100 carry out signals collecting.
Describe function and the configuration of control module 200 in more detail below.Control module 200 has signal receiving module 210th, blood pressure computing module 220, display control module 230 and acquisition control module 240.
Signal receiving module 210 is used for the electrocardiosignal of measurement object and the pulse wave of receipt signal collector 100 collection Signal.The electrocardiosignal of reception and pulse wave signal are sent to blood pressure computing module 220 to enter promoting circulation of blood by signal receiving module 210 Pressure calculates.
The electrocardiosignal that blood pressure computing module 220 is received based on signal receiving module 210 and pulse wave signal are calculated, Obtain corresponding characteristic parameter and be calculated the blood of test object according to the blood pressure computation model of the characteristic parameter prestoring Pressure value.Specifically, blood pressure computing module 220 is used for, according to two groups of synchronous electrocardiosignals and pulse wave signal, extracting pulse wave The fisrt feature point of signal, and the second feature point of the electrocardiosignal in same cardiac cycle, root is obtained according to fisrt feature point Obtain corresponding characteristic value according to fisrt feature point and second feature point, using assessment parameter selected characteristic parameter from characteristic value, According to the blood pressure computation model of the characteristic parameter pre-building, obtain pressure value.Here characteristic parameter includes:Pulse wave is propagated Time, cardiac cycle, rise time in pulse wave systole phase and pulse wave K value etc..Blood pressure computing module 220 is about pressure value meter The detailed process calculated can be found in the related description of hereinafter step 404, step 406, step 408.
And in some of them embodiment of the present invention, blood pressure computing module 220 is also by calculated pressure value etc. Information sends to display control module 230.
In some embodiments of the invention, as shown in Fig. 2 blood pressure computing module 220 includes filter unit 222, feature Value acquiring unit 224, characteristic parameter acquiring unit 226 and blood pressure calculation unit 228.
Filter unit 222, for signal receiving module 210 is received two groups synchronous electrocardiosignals and pulse wave signal It is filtered, to obtain filtered electrocardiosignal and pulse wave signal.Detection process in electrocardiosignal and pulse wave signal In, due to real work situation nonideality. often contain very strong ambient noise in detection signal.These ambient noises The various random noises causing including industrial frequency noise, myoelectricity noise, respiratory wave noise, human action and electrode polarization etc..
Characteristic value acquiring unit 224, for the filtered two groups of synchronous electrocardiosignals of filtered unit 222 and pulse wave Signal, extracts the fisrt feature point of filtered pulse wave signal, then obtains same cardiac cycle according to fisrt feature point The second feature point of interior electrocardiosignal.In the present embodiment, fisrt feature point includes crest and the trough of pulse wave signal.The Two characteristic points include ecg-r wave.
Characteristic parameter acquiring unit 226, for according to fisrt feature point and described second feature point, obtaining corresponding feature Parameter.Here characteristic parameter includes cardiac electrophysiology parameter, pulse wave physiological parameter and electrocardiograph pulse ripple physiological parameter.
Blood pressure calculation unit 228, for the blood pressure computation model according to the characteristic parameter pre-building, obtains pressure value. Here blood pressure includes systolic pressure, diastolic pressure and mean blood pressure.
Display control module 230 controls and shows the various types of numbers being processed by control module 200 on display module 400 It is believed that breath and the result by processing these acquisition of informations.For example, display control module 230 control display module 400 show by The calculated pressure value of blood pressure computing module 220.In addition, display control module 230 controls display module 400 to show by information The information of electrocardiosignal that receiver module 210 receives and pulse wave signal etc..
Acquisition control module 240 control signal collector 100 carries out signals collecting.For example, acquisition control module 240 controls Electrocardial signal acquisition device 110 gathers electrocardiosignal, and acquisition control module 240 controls pulse wave signal collector 120 synchronous simultaneously Ground gathers the pulse wave signal in multiple continuous cardiac cycles.
Schematically joining according to the control module 200 in some embodiments disclosed by the invention described in detail above Put, the configuration of control module 200 is not limited to the example shown in Fig. 1.As long as configuration meets above-mentioned functions it is possible to pass through any Functional block carrys out configuration control module 200.
Next, the memory module 300 of some of them embodiment disclosed by the invention will be described.Storage there is not module 300 Store various types of data of the expression raw body activity gathering with signal picker 100 and/or each of control module 200 process The data of type.Memory module 300 stores such as electrocardio related data and pulse wave related data.Additionally, memory module 300 Can also the calculated pressure value of data storage processing module.Memory module 300 can also be stored in following continuous blood pressures and survey Characteristic point, characteristic parameter and physiological parameter used in amount method, and for information such as the blood pressure models about characteristic parameter. Data processing module is referred to the characteristic point of storage, characteristic parameter and physiological parameter in memory module 300, and for relevant The information such as the blood pressure model of characteristic parameter are calculating blood pressure.
Additionally, in FIG, although giving the example that memory module 300 is arranged in continuous blood pressure measurer 10, It is embodiment not limited to this disclosed by the invention.For example, measurement apparatus 10 can also include the connectivity port of external device (ED) connection (not shown), and it is connected to the external memory unit (not shown) being arranged on outside via this connectivity port.If measurement dress Put 10 and be connected to external memory unit, then the above-mentioned various types of data storages that can will be stored in memory module 300 exist In external memory unit.Here, connectivity port can be card connector, and external memory unit can be storage card (SD Card).
Display module 400 show under the control of display control module 200 by control module 200 process various types of Information.For example, display module 400 can be display.If display module 400 is display, display module 400 can example As in the form of numeral or figure come the electrocardiogram (ECG) data measured by showing and pulse wave data and the pressure value that calculated etc..
Additionally, measurement apparatus 10 may further include the sound of the (not shown) such as composition such as loudspeaker, headphone Frequency output module.If measurement apparatus 10 include dio Output Modules, above-mentioned by display module 400 show various types of Information can be realized by dio Output Modules output alarm tone.Additionally, according to embodiment disclosed by the invention, measurement apparatus 10 Can also include for sending information and the communication module (not shown) from external device (ED) receive information to various external device (ED)s.Example As, the data of various raw body activities of the measurement object that signal picker 100 can be gathered by communication module and control module obtain To the information such as various types of data send to external device (ED).Specifically, signal picker 100 can be gathered by communication module The electrocardiogram (ECG) data of measurement object and the information such as the pressure value that obtains of pulse wave data data processing module send to outside dress Put.Here, the communication module time of various types of data is activations to external device (ED) can be that each measurement is related to raw body work The real time of dynamic various types of data, or various types of data can terminate a series of blood pressure measurement process it Jointly send afterwards.
Additionally, communication module can also receive characteristic point used in following continuous BP measurement process methods, feature ginseng Number and physiological parameter, and for information such as the blood pressure models about characteristic parameter, and by the information Store being received upper In the memory module 300 stated or external memory unit.Therefore, it can to update the blood pressure of these characteristic parameters via communication unit The information such as model.
Additionally, communication module can be using wired or wireless communication mode.If communication module has wireless transmitting function, The radio transmission method then using can be the bluetooth of such as near field communication system or the body area network being standardized IEEE802.15.6.
Additionally, the external device (ED) via communication module and measurement apparatus 10 communication can be PC (personal computer) or intelligence Terminal.These external device (ED)s can be realized and control module 200 identical function.Because external device (ED) has and control module 200 identical functions, then this external device (ED) the various types of data sending from measurement apparatus 10 can be carried out with control mould Block 200 identical is processed.
Below, it is described with reference to Figure 3 Electrocardial signal acquisition device 110 and the pulse wave signal that signal picker 100 includes The illustrative arrangement of collector 120.Fig. 3 shows the structural frames of the illustrative arrangement of signal picker 100 that figure 1 illustrates Figure.Electrocardial signal acquisition device 110 includes electrocardioelectrode 112, ecg signal acquiring module 114 and electrocardiosignal analog-to-digital conversion module 116.Electrocardioelectrode 112 needs direct and human contact, can be made up of the copper sheet having plated one layer of silver chlorate, Main Function is handle Human ecg signal is transmitted to electrocardiogram acquisition module up well.Electrocardiogram acquisition module 114 can adopt electrocardiogram acquisition chip, This electrocardiogram acquisition chip includes preposition amplification, right leg drive, filtering and rearmounted amplifying circuit, at electrocardiogram acquisition module 114 Useful electrocardiosignal after reason, after analog-digital converter, is transferred to signal receiving module 210.Pulse wave signal collector 120 Including pulse wave sensor 122, pulse wave signal acquisition module 124 and pulse signal analog-digital converter 126.Pulse wave sensor 122 can be made up of light emitting diode and photodiode, can adopt two kinds of metering systems of transmission or reflection.Light emitting diode The illumination (green glow/infrared light) of special spectrum can be launched, this illumination is mapped on human vas, be received thoroughly by photodiode It is emitted through blood vessel or the light returned from vasoreflex, the change of the intensity of this transmission/reflection/angle light just reflects human pulse Change, thus collect the pulse wave of human body.Pulse wave signal acquisition module 124 is come in pulse wave sensor 122 transmission Pulse wave signal is sampled, and pulse signal analog-digital converter 126 carries out transmission after analog-to-digital conversion to the pulse wave signal of sampling To signal receiving module 210.
Blood pressure measurement mould in measurement apparatus 10 in above-described embodiment disclosed by the invention to be described below with reference to Fig. 4 Type method for building up.
Step 402, obtains the electrocardiosignal in synchronous multiple continuous cardiac cycle and pulse wave signal.
In some of them embodiment of the present invention, the electrocardiosignal receiving and pulse wave signal are filtered locate Reason, to remove the interference of the noise in electrocardiosignal and pulse wave signal.
Specifically, the characteristic according to electrocardiosignal and pulse wave signal, the method that can adopt wavelet transformation is to acquisition Electrocardiosignal and pulse wave signal filtering, filtered signal is as shown in Figure 5.
Step 404, extracts the fisrt feature point of pulse wave signal, and is obtained same according to the fisrt feature point of pulse wave signal The second feature point of the electrocardiosignal in one cardiac cycle.In some of them embodiment of the present invention, difference threshold can be adopted Value method extracts the fisrt feature point of pulse wave signal, and this fisrt feature point includes crest and the trough of pulse wave signal.Specifically real Existing process is as follows:
(1) obtain the radio-frequency component of pulse wave signal.
By differential process are carried out to pulse wave signal, extract its high-frequency signal.For example, according to the following formula to pulse Ripple signal asks first derivative, dPPG (i)=PPG (i+1)-PPG (i), that is, be equivalent to and carry out high-pass filtering to it.
Wherein, PPG (i) represents the amplitude of i-th point of signal in PPG waveform, and dPPG represents the first-order difference of PPG signal.
(2) segmentation is carried out to the radio-frequency component of pulse wave signal.
Pulse frequency due to human body changes in certain scope, and in general, pulse frequency all can be more than 30 beats/min.Thus, right Pulse wave signal waveform carries out segmentation, such as with 2s for one section, thus all can have a pulse cycle in ensureing every section.
(3) signal segment that above-mentioned steps (2) are obtained carries out threshold test.
Search the maximum of every section of signal amplitude;
Acquisition amplitude exceedes the position of all data points (BigNumber) of maximum prearranged multiple, in the present embodiment, Prearranged multiple can be 0.3 times.
The data point of above-mentioned acquisition is traveled through, the difference searching adjacent two data point position exceedes scheduled time length Data point.Because if the position of two data points is more than first scheduled time it becomes possible to determine that this two data points are located at On tandem two different pulse waves.In the present embodiment, first scheduled time can be 0.2s.
(4) search the crest of pulse wave signal.
The position of the previous data point that above-mentioned (3) are detected is defined as Label such that it is able to determine pulse crest value Position just near Label.Shown by statistics, pulse wave amplitude exceedes the data point of maximum amplitude prearranged multiple, its time Persistence length is not over second scheduled time.Therefore, traveled through in second scheduled time after Label, found out Maximum during this, finally determines the crest of pulse wave.In certain embodiments, second scheduled time can be 0.4s.
(5) search the trough of pulse wave signal.
Some characteristics being had in itself due to pulse wave, enabling substantially determine pulse wave trough position crest it In the range of front second scheduled time.Then, search the minimum of a value in this interval, so that it is determined that the trough of pulse wave.
Fisrt feature point extraction is carried out by said process to pulse wave signal, the result obtaining is as shown in Figure 6.
The second feature of the electrocardiosignal in same cardiac cycle after determining fisrt feature point, is obtained according to fisrt feature point Point.Second feature point includes cardiac electrical R wave characteristic point.The method of the acquisition of second feature point of electrocardiosignal can be realized relatively Many, in some embodiments of the invention, it is worth most as a example method and Wavelet Transform by interval and illustrates.
(1) interval is worth method most
Because electrocardio and pulse wave are to be produced by same cardiac cycle, therefore when the peak-to-valley value of pulse wave all determines Afterwards, cardiac electrical R ripple can just be positioned.It concretely comprises the following steps:
After the Wave crest and wave trough of the PPG that has good positioning, between two adjacent PPG peak time, traversal corresponding ECG letter Number, the R ripple of maximizing point wherein, as ECG.Its effect is as shown in Figure 7.
(2) the R ripple detecting step based on Wavelet Transform is as follows:
Multi-resolution decomposition process:Choose quadratic spline B small echo, for morther wavelet, little multi-resolution decomposition is carried out to electrocardiosignal, In the present embodiment, multi-resolution decomposition is 24 Scale Decompositions;
Maximum detection process:Because the extreme portions energy of R ripple concentrates on yardstick 3, the present invention will enter in yardstick 3 The detection of row maximum;The wavelet coefficient d3 of yardstick 3 is scanned, coefficient value is more than 0 and big with left adjoint point line slope Point in 0 constitutes positive maximum Candidate Set PosC.Simultaneously by coefficient value be less than 0 and with left adjoint point line slope be less than 0 point structure Become negative maximum Candidate Set NegC;
Binarization:First threshold (th1) and Second Threshold (th2) are determined based on Candidate Set mean value method, and will just In maximum point Candidate Set PosC, every point more than first threshold th1 is stored in positive maximum collection Pos, and its value is set to 1. Point less than Second Threshold th2 every in negative maximum point Candidate Set NegC is stored in negative maximum collection Neg, and its value is put For -1;
The pairing process of positive and negative maximum pair:It is contemplated that at R crest value point after generating positive maximum collection and negative maximum collection In maximum between interval in.So needing to carry out the pairing work of positive and negative maximum pair.Defining very big value set is Loca, all data that positive maximum collection and negative maximum are concentrated are stored in loca.Find at adjacent 2 points in loca set N1, N2, and N1 value is 1 for -1, N2 value, i.e. adjacent two points negative maximum points, one is negative maximum point.Define modulus maxima Value is combined into Pair to collection.If the actual range of N1, N2 is less than at 80 points, then the physical location of N1 is stored in Pair.
Candidate point acquisition process:According to wavelet conversion characteristics, the corresponding data point of R crest value is through multi-resolution decomposition Afterwards, zero crossing that can be corresponding to the maximum on a certain yardstick to line.But signal can produce displacement in wavelet transformation, that is, Point corresponding to R wave crest point reality has the zero crossing that certain probability is not maximum pair, has certain distance between 2 points.But R Wave crest point is necessarily in interval interior between them.So, in some of them embodiment of the present invention, R wave crest point is carried out Positioning, using analysis maximum between interval range method.In interval range, the point of true amplitude maximum is set to R crest Candidate point.Definition R wave crest point Candidate Set is R_C, and the candidate point (i.e. second feature point) meeting threshold condition is stored in Candidate Set For in R_C.
With Wavelet Transform, electrocardiosignal is positioned, obtain waveform as shown in Figure 8.
Step 406, according to fisrt feature point and described second feature point, obtains corresponding characteristic parameter.
According to above-mentioned steps, the fisrt feature point obtaining and second feature point calculate corresponding characteristic value respectively, and these are special Value indicative includes ecg characteristics, pulse wave characteristic, and electrocardiograph pulse ripple fusion feature.These ecg characteristics, pulse wave characteristic and the heart Electric pulse wave fusion feature includes cardiac electrophysiology parameter, pulse wave physiological parameter and electrocardiograph pulse ripple physiological parameter.
Wherein, cardiac electrophysiology parameter mainly includes but are not limited to:
Heart rate (heart rate, HR), refers to the heart number of times beated per minute;
Cardiac electrical cycle (EcgPeriod), represents the time span of each cardiac cycle;
SDNN:The standard deviation of phase (i.e. cardiac cycle) sequence between NN;
pNN50:The number that between NN, the difference of interim two neighboring phase is more than 50 milliseconds, accounts for phase sum between all of NN Percentage;
HRV triangle index:Between NN between total heart rate of phase and NN phase histogram height business;
LF:The integration of low-frequency range (0.04Hz-0.15Hz) in phase power spectral density plot between NN;
HF:The integration of high band (0.15Hz-0.4Hz) in phase power spectral density plot between NN;
LF/HF:The ratio of low-frequency range energy and high band energy in phase power spectral density plot between NN.
Pulse wave physiological parameter mainly includes but are not limited to:
Rise time in systole phase (ppgSystolicTime):In same cardiac cycle, the trough of pulse wave to crest when Between be spaced (as shown in t1 in Fig. 9);
Fall time diastole (ppgDiastolicTime):In same period, the crest of pulse wave is to next trough Time interval (as shown in t2 in Fig. 9);
Pulse wave K value:The mean value (Pm-Pd) of pulse wave pressure fluctuation component, in flutter component maximum (Ps-Pd) Shared percentage (as shown in Figure 10);
The first derivative of pulse wave (PPG), second dervative characteristic value;
Electrocardiograph pulse ripple physiological parameter:Pulse wave propagation time (Pulse Wave Transit Time, PWTT):It is ECG R wave crest point to PPG initiate time interval (as Figure 11).
Further, using assessment parameter selected characteristic parameter from the characteristic value of above-mentioned acquisition.
For realizing the low complex degree matching of continuous BP measurement, uncorrelated or redundancy characteristic value need to be rejected, thus reaching Reduce Characteristic Number, improve measuring accuracy, reduce the purpose running complexity.
In some of which embodiment of the present invention, characteristic parameter selection process is as follows:
(1) according to characteristic value initialization feature parameter set and characteristic parameter subset.
Specifically, characteristic parameter collection F={ f1,f2,...,fNInclude characteristic value acquiring unit 24 can obtain all of Cardiac electrophysiology parameter, pulse wave physiological parameter and electrocardiograph pulse ripple physiological parameter.Characteristic parameter subset F ' it is initially set to empty set.
(2) select evaluation function, and multiple characteristic parameter subsets are generated according to this evaluation function.
In some of which embodiment of the present invention, using the evaluation function based on information gain.Wherein, information gain Computational methods are as follows:
Assume existing characteristics subset A and character subset B, classified variable is C, then comentropy H (C) of classified variable C can It is expressed as:
By feature FjConditional information entropy H for classified variable C (C | Fj) be expressed as:
Then select feature FjThe change of the comentropy of C in front and back becomes the information gain (Information Gain) of C, can It is expressed as:IG(C|Fj)=H (C)-H (C | Fj).
In some embodiments of the invention, by the search to proper subspace, the feature selecting side based on decision tree Method, is maximized based on information gain and selects a characteristic parameter as the Split Attribute of decision tree, recursively give birth to from top to bottom Become child node, until data set is inseparable, stop decision tree and stop growing.Multiple characteristic parameter are generated by said method Collection, above-mentioned search procedure can be using search, heuristic search or random search completely.
(3) decision tree beta pruning
Using the evaluation function based on information gain defined above, beta pruning is carried out to set up decision tree, gained is finally determined Feature at each branch node of plan tree is optimal characteristics collection.That is, the multiple characteristic parameter subset choosings obtaining from above-mentioned steps (2) Take optimal characteristics subset of parameters.
In some of them embodiment of the present invention, rear beta pruning is carried out to decision tree based on evaluation function, and return optimal Fork attribute is as optimal characteristics subset of parameters.In addition, evaluation function can be selected based on screening washer and wrapper.Conventional evaluation Function includes:Correlation, distance, information gain, uniformity, grader error rate.
In some of them embodiment of the present invention, the final characteristic parameter (i.e. optimal characteristics subset of parameters) determining includes: Pulse wave propagation time, cardiac cycle, rise time in pulse wave systole phase and pulse wave K value.
Step 408, based on corresponding characteristic parameter, sets up the blood pressure computation model of characteristic parameter.
In some embodiments of the invention, pressure value at least includes:Systolic pressure, diastolic pressure and mean blood pressure.
The blood pressure computation model of the characteristic parameter pre-building includes multivariate linear model and nonlinear model.
In some embodiments of the invention, multivariate linear model is the multinomial physiological parameter based on human body and characteristic parameter Set up, it can obtain the corresponding result of characteristic ginseng value.This multinomial physiological parameter includes:Sex, the age, height, body weight, Physique than with brachium at least one.
For example, on training sample set, based on pulse wave translation time, cardiac electrical cycle, time in pulse wave systole phase, pulse Wave characteristic K value, sets up systolic pressure (SBP), the multiple linear regression model of mean blood pressure (MAP) and diastolic pressure (DBP) is such as respectively Under:
SBP=a1+a2×PWTT+a3×EcgPeriod+a4×PpgSystolicTime+a5× ppg_K, (1)
MAP=b1+b2×PWTT+b3×EcgPeriod+b4×PpgSystolicTime+b5× ppg_K, (2)
DBP=c1+c2×PWTT+c3×EcgPeriod+c4×PpgSystolicTime+c5× ppg_K, (3)
For different subjects, the coefficient a between themi,bi,ciDifferent.Therefore will for different subjects, To its regression coefficient ai,bi,ciIt is analyzed.Excavated respectively by the method for the data separate statistical analysis of substantial amounts of subject Relevance between individual regression coefficient and the tested body physiological parameter of itself, from every physiological parameter of subject, picks out Sex (gender), age (age), height (height), body weight (weight), physique ratio (BMI) and brachium (armLength) physiological parameter such as, sets up multilinear fitting regression model.Regression coefficient ai,bi,ciWith various physiological parameters Relation can be expressed as:
Above-mentioned model formation (1), (2) and (3) can develop into following formula (4), (5), (6):
SBP=P D G, (4)
MAP=P E G, (5)
DBP=P F G, (6)
Wherein, SBP is systolic pressure, and MAP is mean blood pressure, and DBP is diastolic pressure, P=[1 PWTT EcgPeriod PpgSystolicTime ppg_K], PWTT is pulse wave propagation time, and EcgPeriod is cardiac cycle, PpgSystolicTime is the rise time in pulse wave systole phase, and ppg_K is pulse wave K value, and D, E and F are respectively systolic pressure, relax Open pressure and mean blood pressure corresponding bp coefficient matrix,
Gender is sex, and age is the age, and height is height, and weight is body weight, and BMI is physique ratio, ArmLength is brachium.
On training dataset, corresponding bp coefficient matrix D, E and F can be obtained, and set up corresponding characteristic parameter Blood pressure computation model.
In some embodiments of the invention, the model of characteristic parameter includes nonlinear model, and nonlinear model includes base In the nonlinear regression model (NLRM) of SVMs, the step setting up the blood pressure computation model of characteristic parameter includes:
1) obtain the training set between characteristic parameter and systolic pressure.
If known training set { (x1,SBP1),(x2,SBP2),...,(xN,SBPN), wherein
Xi={ PWTTi,EcgPeriodi,PpgSystolicTimei,PPG_Ki};
2) SVMs commonly use kernel function include linearly, multinomial inner product, Sigmoid inner product, Radial basis kernel function etc., Kernel function can by former feature space linearly inseparable Projection Character to another space to realize effectively classifying.In the present invention In, RBF is selected as kernel function and slack variable (ε, C) based on extensive statistics result;
3) construct and solve optimization problem
Obtain optimal solution
4) construct decision function
Wherein, b can obtain in the following manner,
5) on training set, corresponding systolic pressure and diastolic pressure coefficient matrix can be obtainedWith b.The computing formula of final blood pressure is represented by:
By processing to the electrocardio of input and pulse wave signal, extract corresponding characteristic value, and be updated to above-mentioned The blood pressure computation model of the characteristic parameter obtaining, can obtain the pressure value of measurement object, such as systolic pressure, diastolic pressure and average Blood pressure.
As shown in figure 12, also propose a kind of continuous BP measurement model in the present invention and set up system 500, this system 500 is wrapped Include:Signal acquisition module 502, feature point extraction module 504, characteristic value acquisition module 506 and model building module 508.
Signal acquisition module 502, for receiving electrocardiosignal and arteries and veins in two groups of synchronous multiple continuous cardiac cycles Fight ripple signal;
Feature point extraction module 504, for extracting the fisrt feature point of the pulse wave signal of signal acquisition module acquisition, and Obtain the second feature point of the electrocardiosignal in same cardiac cycle according to fisrt feature point;
Characteristic value acquisition module 506, for according to fisrt feature point and second feature point, obtaining corresponding characteristic parameter;
Model building module 508, for setting up the blood pressure computation model of characteristic parameter.
The continuous BP measurement model of the present embodiment is set up system 500 and is used for realizing aforesaid continuous BP measurement model building Cube method, therefore continuous BP measurement model are set up being embodied as in system and be can be found in continuous BP measurement model foundation above The embodiment part of method, for example, signal acquisition module 502, feature point extraction module 504, characteristic value acquisition module 506 and mould Type set up module 508 be respectively used to realize step 402 in above-mentioned continuous BP measurement method for establishing model, 404,406 and 408, So, its specific implementation can refer to the description of hereinbefore each embodiment about step 402,404,406 and 408, This is not repeated.
Fig. 4 is the schematic flow sheet of the continuous BP measurement method for establishing model of one embodiment of the invention.It should be appreciated that Be although each step in the flow chart of Fig. 4 shows successively according to the instruction of arrow, but these steps are not certainty Order according to arrow instruction executes successively.Unless expressly stated otherwise herein, the execution of these steps is not strict Order limits, and it can execute in the other order.And, at least a portion step in Fig. 4 can include many sub-steps Or multiple stages, these sub-steps or stage are not necessarily to complete in synchronization execution, but can be different Moment executes, and its execution sequence is also not necessarily and carries out successively, but can be with the sub-step of other steps or other steps Or at least a portion executed in parallel in stage or alternately execute.
Above each embodiment only implementation just for corresponding steps in illustrating is set forth, Ran Hou In the case of logic is not conflicting, each embodiment above-mentioned be can be mutually combined and form new technical scheme, and be somebody's turn to do New technical scheme is still in the open scope of this specific embodiment.
Through the above description of the embodiments, those skilled in the art can be understood that above-described embodiment side Method can be realized by the mode of software plus necessary general hardware platform naturally it is also possible to pass through hardware, but in many cases The former is more preferably embodiment.Based on such understanding, technical scheme is substantially done to prior art in other words Go out partly can embodying in the form of software product of contribution, this computer software product is carried on a non-volatile meter In calculation machine readable storage medium (as ROM, magnetic disc, CD, server cloud space), including some instructions with so that a station terminal Equipment (can be mobile phone, computer, server, or network equipment etc.) executes the method described in each embodiment of the present invention.
Each technical characteristic of embodiment described above can arbitrarily be combined, for making description succinct, not to above-mentioned reality The all possible combination of each technical characteristic applied in example is all described, as long as however, the combination of these technical characteristics is not deposited In contradiction, all it is considered to be the scope of this specification record.
Embodiment described above only have expressed the several embodiments of the present invention, and its description is more concrete and detailed, but simultaneously Can not therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art Say, without departing from the inventive concept of the premise, some deformation can also be made and improve, these broadly fall into the protection of the present invention Scope.Therefore, the protection domain of patent of the present invention should be defined by claims.

Claims (15)

1. a kind of continuous blood pressure measurer, described device includes:
Signal picker is including Electrocardial signal acquisition device and pulse wave signal collector, many for synchronous acquisition measurement object Electrocardiosignal in individual continuous cardiac cycle and pulse wave signal, obtain two groups of synchronous electrocardiosignals and pulse wave signal;
Signal receiving module, for receiving two groups of synchronous electrocardiosignals and pulse wave signal;
Blood pressure computing module, for the electrocardiosignal synchronous according to two groups and pulse wave signal, extracts described pulse wave signal Fisrt feature point, and the second feature point of the electrocardiosignal in same cardiac cycle is obtained according to described fisrt feature point, according to Described fisrt feature point and described second feature point obtain corresponding characteristic value, are chosen from described characteristic value using assessment parameter Characteristic parameter, according to the blood pressure computation model of the described characteristic parameter pre-building, obtains pressure value.
2. device according to claim 1 is it is characterised in that described blood pressure computing module includes:
Characteristic value acquiring unit, for the electrocardiosignal synchronous according to two groups and pulse wave signal, extracts described pulse wave signal Fisrt feature point, and according to described fisrt feature point obtain same cardiac cycle in electrocardiosignal second feature point;
Characteristic parameter acquiring unit, for according to described fisrt feature point and described second feature point, obtaining corresponding feature ginseng Number;With
Blood pressure calculation unit, for the blood pressure computation model according to the described characteristic parameter pre-building, obtains pressure value.
3. device according to claim 1 is it is characterised in that described blood pressure computing module also includes:
Filter unit, for being filtered processing to described two groups synchronous electrocardiosignals and pulse wave signal respectively.
4. device according to claim 1 is it is characterised in that described characteristic parameter includes:Pulse wave propagation time, aroused in interest Cycle, rise time in pulse wave systole phase and pulse wave K value.
5. the device according to claim 2 or 4 it is characterised in that the blood pressure computation model of described characteristic parameter include many First linear model, described multivariate linear model is:
SBP=P D G,
MAP=P E G,
DBP=P F G,
Wherein, SBP is systolic pressure, and MAP is mean blood pressure, and DBP is diastolic pressure, P=[1 PWTT EcgPeriod PpgSystolicTime ppg_K], PWTT is described pulse wave propagation time, and EcgPeriod is described cardiac cycle, PpgSystolicTime is the described rise time in pulse wave systole phase, and ppg_K is described pulse wave K value, and D, E and F are respectively institute State systolic pressure, described diastolic pressure and described mean blood pressure corresponding bp coefficient matrix, G is the physiological parameter square of measurement object Battle array, G = 1 g e n d e r a g e h e i g h t w e i g h t B M I a r m L e n g t h , Gender is sex, and age is the age, and height is height, and weight is body weight, and BMI is physique ArmLength is brachium to ratio.
6. the device according to claim 2 or 4 it is characterised in that the blood pressure computation model of described characteristic parameter include non- Linear model, described nonlinear model includes the nonlinear regression model (NLRM) based on SVMs, described blood pressure calculation unit base Process in the nonlinear regression model (NLRM) acquisition pressure value of SVMs includes:
Obtain the training set between described characteristic parameter and described systolic pressure;
Construct decision function using the described nonlinear regression model (NLRM) based on SVMs and calculate optimum solution;
Obtain systolic pressure and the diastolic pressure coefficient matrix based on described data set according to described decision function;
Final pressure value is obtained according to described systolic pressure and diastolic pressure coefficient matrix.
7. device according to claim 2 is it is characterised in that described characteristic parameter selection unit is additionally operable to according to described spy Value indicative initialization feature parameter set and characteristic parameter subset, choose multiple candidate feature subset, using described valuation functions to institute State multiple candidate feature subset to be estimated, obtain the characteristic parameter conforming to a predetermined condition.
8. device according to claim 1 is it is characterised in that described device also includes:
Memory module, for storing described characteristic value, described characteristic parameter and its corresponding pressure value;
Display module, for showing described characteristic value, described characteristic parameter and its corresponding pressure value.
9. the device described in 8 is wanted it is characterised in that described device also includes according to right:
Communication module, for receiving two groups of synchronous electrocardiosignals and the pulse wave signal of described signal picker collection, and will The various data is activations of memory module storage are to external device (ED);Or,
It is respectively provided with wireless signal Transmit-Receive Unit in signal picker and signal receiving module, signal picker and signal receive Carried out data transmission by wireless signal between module and communicate.
10. a kind of continuous BP measurement system, described system includes:
Signal acquisition module, for receiving the electrocardiosignal in two groups of synchronous multiple continuous cardiac cycles and pulse wave letter Number;
Feature point extraction module, for extracting the fisrt feature point of the described pulse wave signal that described signal acquisition module obtains, And the second feature point of the electrocardiosignal in same cardiac cycle is obtained according to described fisrt feature point;
Characteristic value acquisition module, for according to described fisrt feature point and described second feature point, obtaining corresponding characteristic parameter;
Model building module, for setting up the blood pressure computation model of described characteristic parameter.
A kind of 11. method for building up of continuous BP measurement model, methods described includes:
Obtain the electrocardiosignal in synchronous multiple continuous cardiac cycle and pulse wave signal;
Extract the fisrt feature point of described pulse wave signal, and obtained with wholeheartedly according to the fisrt feature point of described pulse wave signal The second feature point of the electrocardiosignal in the dynamic cycle;
According to described fisrt feature point and described second feature point, obtain corresponding characteristic parameter;
Based on corresponding characteristic parameter, set up the blood pressure computation model of described characteristic parameter.
12. methods according to claim 11 are it is characterised in that described characteristic parameter includes:Pulse wave propagation time, the heart Dynamic cycle, rise time in pulse wave systole phase and pulse wave K value.
13. methods according to claim 11 or 12 are it is characterised in that the blood pressure computation model of described characteristic parameter includes Multivariate linear model, described multivariate linear model is:
SBP=P D G,
MAP=P E G,
DBP=P F G,
Wherein, SBP is described systolic pressure, and MAP is described mean blood pressure, and DBP is described diastolic pressure, P=[1 PWTT EcgPeriod PpgSystolicTime ppg_K], PWTT is described pulse wave propagation time, and EcgPeriod is described aroused in interest In the cycle, PpgSystolicTime is the described rise time in pulse wave systole phase, and ppg_K is described pulse wave K value, and D, E and F divide Not Wei described systolic pressure, described diastolic pressure and described mean blood pressure corresponding bp coefficient matrix, G = 1 g e n d e r a g e h e i g h t w e i g h t B M I a r m L e n g t h , Gender is described sex, and age is the described age, and height is described height, and weight is described body weight, and BMI is described body Matter ratio, armLength is described brachium.
14. methods according to claim 11 or 12 are it is characterised in that the blood pressure computation model of described characteristic parameter includes Nonlinear model, described nonlinear model includes the nonlinear regression model (NLRM) based on SVMs, described sets up described feature The step of the blood pressure computation model of parameter includes:
Obtain the training set between described characteristic parameter and described systolic pressure;
Construct decision function using the described nonlinear regression model (NLRM) based on SVMs and calculate optimum solution;
Obtain systolic pressure and the diastolic pressure coefficient matrix based on described data set according to described decision function;
Final pressure value is obtained according to described systolic pressure and diastolic pressure coefficient matrix.
15. methods according to claim 11 are it is characterised in that characteristic parameter obtaining step includes:
According to characteristic value initialization feature parameter set and characteristic parameter subset;
Select evaluation function, and multiple characteristic parameter subsets are generated according to described evaluation function;
Optimal characteristics subset of parameters is chosen from the plurality of characteristic parameter subset according to described evaluation function.
CN201580000417.0A 2015-08-08 2015-08-08 Continuous blood pressure measuring device, measuring model establishing method and system Active CN106413534B (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2015/086429 WO2017024457A1 (en) 2015-08-08 2015-08-08 Blood-pressure continuous-measurement device, measurement model establishment method, and system

Publications (2)

Publication Number Publication Date
CN106413534A true CN106413534A (en) 2017-02-15
CN106413534B CN106413534B (en) 2020-06-30

Family

ID=57982845

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580000417.0A Active CN106413534B (en) 2015-08-08 2015-08-08 Continuous blood pressure measuring device, measuring model establishing method and system

Country Status (2)

Country Link
CN (1) CN106413534B (en)
WO (1) WO2017024457A1 (en)

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106974631A (en) * 2017-03-17 2017-07-25 深圳大学 A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal
CN107361753A (en) * 2017-08-29 2017-11-21 哈尔滨理工大学 Health state monitoring method based on pulse wave characteristic point
CN107550481A (en) * 2017-08-24 2018-01-09 京东方科技集团股份有限公司 A kind of portable equipment and blood pressure measuring method
CN107692993A (en) * 2017-11-08 2018-02-16 深圳市玉成创新科技有限公司 A kind of blood pressure instrument
CN107854123A (en) * 2017-12-13 2018-03-30 清华大学深圳研究生院 One kind is without oversleeve continuous blood pressure monitoring method and apparatus
WO2018152713A1 (en) * 2017-02-22 2018-08-30 清华大学深圳研究生院 Data processing method for use with blood pressure measuring devices
CN108670231A (en) * 2018-03-14 2018-10-19 深圳竹信科技有限公司 Blood pressure measuring method, terminal and computer readable storage medium
CN108720821A (en) * 2018-05-28 2018-11-02 京东方科技集团股份有限公司 A kind of blood pressure estimation device and blood pressure estimating system
CN108926334A (en) * 2017-05-26 2018-12-04 深圳市玉成创新科技有限公司 Blood pressure acquisition methods and its system and device based on pulse wave
CN108937896A (en) * 2017-05-23 2018-12-07 深圳先进技术研究院 A kind of central aortic blood pressure measuring method and equipment
CN108992049A (en) * 2018-06-14 2018-12-14 深圳鑫想科技有限责任公司 A kind of smart phone blood pressure test method and system
CN109009048A (en) * 2018-06-12 2018-12-18 四川斐讯信息技术有限公司 A kind of PPG peak detection method and system
CN109259745A (en) * 2018-10-25 2019-01-25 贵州医科大学附属医院 A kind of wearable cardiovascular and cerebrovascular disease intelligent monitor system and method
CN109350031A (en) * 2018-12-07 2019-02-19 深圳先进技术研究院 A kind of photoplethysmographic signal feature point detecting method and device
WO2019105105A1 (en) * 2017-11-28 2019-06-06 深圳市岩尚科技有限公司 Classification prediction data processing method for photoplethysmography-based blood pressure measurement device
CN109890276A (en) * 2017-04-01 2019-06-14 华为技术有限公司 Monitoring of blood pressure method, apparatus and equipment
CN110251105A (en) * 2019-06-12 2019-09-20 广州视源电子科技股份有限公司 A kind of non-invasive blood pressure measuring method, device, equipment and system
CN110367957A (en) * 2018-04-12 2019-10-25 三星电子株式会社 Biological information measurement equipment and method for measuring biological information
CN110464345A (en) * 2019-08-22 2019-11-19 北京航空航天大学 A kind of separate head bioelectrical power signal interference elimination method and system
CN110664386A (en) * 2019-09-04 2020-01-10 深圳先进技术研究院 Acquisition device and method for pulse wave signals
CN110974170A (en) * 2019-12-13 2020-04-10 中国人民解放军63919部队 Physiological parameter detection system
CN111065321A (en) * 2017-09-12 2020-04-24 欧姆龙株式会社 Pulse wave measuring device, blood pressure measuring apparatus, pulse wave measuring method, and blood pressure measuring method
CN111297341A (en) * 2020-02-20 2020-06-19 京东方科技集团股份有限公司 Developments blood pressure check out test set and pulse wave characteristic extraction equipment
CN111631733A (en) * 2020-06-19 2020-09-08 浙江澍源智能技术有限公司 Arterial blood spectrum detection method and device
CN111904404A (en) * 2020-08-06 2020-11-10 苏州国科医工科技发展(集团)有限公司 Blood pressure regulation and control equipment containing closed-loop monitoring and based on ear vagus nerve stimulation
CN111973165A (en) * 2020-08-14 2020-11-24 北京航空航天大学 Linear and nonlinear mixed non-invasive continuous blood pressure measuring system based on PPG
CN112656393A (en) * 2020-12-08 2021-04-16 山东中科先进技术研究院有限公司 Heart rate variability detection method and system
CN113208573A (en) * 2021-04-21 2021-08-06 北京雪扬科技有限公司 Support wearable equipment of PPG + ECG function
CN113520350A (en) * 2021-07-27 2021-10-22 香港心脑血管健康工程研究中心有限公司 Processing method and device for acquiring relevant characteristic parameters and index information of blood pressure map signals
CN113854985A (en) * 2021-08-27 2021-12-31 联卫医疗科技(上海)有限公司 Method and device for obtaining machine learning model samples for blood pressure prediction
CN113855007A (en) * 2021-08-27 2021-12-31 联卫医疗科技(上海)有限公司 Method and device for obtaining machine learning model samples for blood glucose prediction
CN114145724A (en) * 2021-12-08 2022-03-08 四川北易信息技术有限公司 Method for dynamically monitoring blood pressure based on ECG (electrocardiogram) and PPG (photoplethysmography) multiple physiological characteristic parameters
CN114305364A (en) * 2022-01-05 2022-04-12 北京科技大学 Blood pressure detection method, system and equipment based on millimeter wave radar
CN115670408A (en) * 2022-12-28 2023-02-03 可孚医疗科技股份有限公司 Blood pressure measuring device, linear model coefficient self-correction method and system thereof, and measuring method
CN113520350B (en) * 2021-07-27 2024-05-17 香港心脑血管健康工程研究中心有限公司 Processing method and device for acquiring relevant characteristic parameters and index information of blood pressure chart signals

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106343976B (en) * 2016-09-14 2018-09-07 京东方科技集团股份有限公司 The method and apparatus established the method and apparatus of blood pressure model and determine blood pressure
CN107184194B (en) * 2017-07-06 2023-03-10 中国科学院合肥物质科学研究院 Automatic blood pressure measuring device and method based on digital control
TWI729387B (en) 2019-05-03 2021-06-01 華碩電腦股份有限公司 Wearable device for physiological data detection
CN111248881B (en) * 2020-02-21 2022-08-02 乐普(北京)医疗器械股份有限公司 Method and device for predicting blood pressure
CN111419205A (en) * 2020-03-12 2020-07-17 天津大学 Three-element cuff-free continuous blood pressure detection system based on artificial neural network
CN112932429B (en) * 2021-02-02 2022-08-26 成都泰盟软件有限公司 Automatic real-time calibration method for continuous blood pressure measurement
CN113288090B (en) * 2021-05-06 2022-04-22 广东工业大学 Blood pressure prediction method, system, device and storage medium
CN113303776B (en) * 2021-05-28 2024-02-09 西安理工大学 Non-contact blood pressure measurement method based on circulating neural network
CN113647924A (en) * 2021-08-20 2021-11-16 刘明 Arm-wearing type blood pressure meter based on hongmeng operating system
CN114515147B (en) * 2021-12-20 2023-10-03 之江实验室 Physiological monitoring system based on BCG signal and PPG signal fusion

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101193588A (en) * 2005-03-21 2008-06-04 海尔思-斯玛特有限公司 System for continuous blood pressure monitoring
CN102008296A (en) * 2010-12-24 2011-04-13 吉林大学 Device and method for measuring arterial blood pressures based on pulse wave signals and electrocardiosignals
CN102397064A (en) * 2011-12-14 2012-04-04 中国航天员科研训练中心 Continuous blood pressure measuring device
CN103519794A (en) * 2012-07-04 2014-01-22 索尼公司 Measurement apparatus, measurement method, program, storage medium, and measurement system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4869262A (en) * 1987-02-02 1989-09-26 Pulse Time Products Limited Device for displaying blood pressure
JP3318727B2 (en) * 1994-06-06 2002-08-26 日本光電工業株式会社 Pulse wave transit time sphygmomanometer
CN103385702B (en) * 2013-07-26 2015-08-26 中国科学院深圳先进技术研究院 A kind of non-invasive blood pressure continuous detection apparatus and method
JP6435128B2 (en) * 2013-07-26 2018-12-05 タタ・コンサルタンシー・サーヴィシズ・リミテッド Physiological parameter monitoring

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101193588A (en) * 2005-03-21 2008-06-04 海尔思-斯玛特有限公司 System for continuous blood pressure monitoring
CN102008296A (en) * 2010-12-24 2011-04-13 吉林大学 Device and method for measuring arterial blood pressures based on pulse wave signals and electrocardiosignals
CN102397064A (en) * 2011-12-14 2012-04-04 中国航天员科研训练中心 Continuous blood pressure measuring device
CN103519794A (en) * 2012-07-04 2014-01-22 索尼公司 Measurement apparatus, measurement method, program, storage medium, and measurement system

Cited By (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018152713A1 (en) * 2017-02-22 2018-08-30 清华大学深圳研究生院 Data processing method for use with blood pressure measuring devices
CN106974631B (en) * 2017-03-17 2020-01-07 深圳大学 Blood pressure measuring method and device based on pulse wave waveform and electrocardiosignal
CN106974631A (en) * 2017-03-17 2017-07-25 深圳大学 A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal
CN109890276A (en) * 2017-04-01 2019-06-14 华为技术有限公司 Monitoring of blood pressure method, apparatus and equipment
CN109890276B (en) * 2017-04-01 2021-05-18 华为技术有限公司 Blood pressure monitoring method, device and equipment
CN108937896A (en) * 2017-05-23 2018-12-07 深圳先进技术研究院 A kind of central aortic blood pressure measuring method and equipment
CN108926334A (en) * 2017-05-26 2018-12-04 深圳市玉成创新科技有限公司 Blood pressure acquisition methods and its system and device based on pulse wave
CN107550481A (en) * 2017-08-24 2018-01-09 京东方科技集团股份有限公司 A kind of portable equipment and blood pressure measuring method
CN107361753A (en) * 2017-08-29 2017-11-21 哈尔滨理工大学 Health state monitoring method based on pulse wave characteristic point
CN111065321A (en) * 2017-09-12 2020-04-24 欧姆龙株式会社 Pulse wave measuring device, blood pressure measuring apparatus, pulse wave measuring method, and blood pressure measuring method
CN107692993A (en) * 2017-11-08 2018-02-16 深圳市玉成创新科技有限公司 A kind of blood pressure instrument
WO2019105105A1 (en) * 2017-11-28 2019-06-06 深圳市岩尚科技有限公司 Classification prediction data processing method for photoplethysmography-based blood pressure measurement device
CN107854123A (en) * 2017-12-13 2018-03-30 清华大学深圳研究生院 One kind is without oversleeve continuous blood pressure monitoring method and apparatus
CN107854123B (en) * 2017-12-13 2021-02-02 清华大学深圳研究生院 Sleeveless continuous blood pressure monitoring method and device
CN108670231A (en) * 2018-03-14 2018-10-19 深圳竹信科技有限公司 Blood pressure measuring method, terminal and computer readable storage medium
CN110367957A (en) * 2018-04-12 2019-10-25 三星电子株式会社 Biological information measurement equipment and method for measuring biological information
CN108720821A (en) * 2018-05-28 2018-11-02 京东方科技集团股份有限公司 A kind of blood pressure estimation device and blood pressure estimating system
CN109009048A (en) * 2018-06-12 2018-12-18 四川斐讯信息技术有限公司 A kind of PPG peak detection method and system
CN108992049A (en) * 2018-06-14 2018-12-14 深圳鑫想科技有限责任公司 A kind of smart phone blood pressure test method and system
CN109259745A (en) * 2018-10-25 2019-01-25 贵州医科大学附属医院 A kind of wearable cardiovascular and cerebrovascular disease intelligent monitor system and method
CN109350031A (en) * 2018-12-07 2019-02-19 深圳先进技术研究院 A kind of photoplethysmographic signal feature point detecting method and device
WO2020114448A1 (en) * 2018-12-07 2020-06-11 深圳先进技术研究院 Photoplethysmography signal feature point detecting method and device
CN110251105A (en) * 2019-06-12 2019-09-20 广州视源电子科技股份有限公司 A kind of non-invasive blood pressure measuring method, device, equipment and system
CN110251105B (en) * 2019-06-12 2022-06-21 广州视源电子科技股份有限公司 Noninvasive blood pressure measuring method, device, equipment and system
CN110464345A (en) * 2019-08-22 2019-11-19 北京航空航天大学 A kind of separate head bioelectrical power signal interference elimination method and system
CN110664386A (en) * 2019-09-04 2020-01-10 深圳先进技术研究院 Acquisition device and method for pulse wave signals
CN110974170A (en) * 2019-12-13 2020-04-10 中国人民解放军63919部队 Physiological parameter detection system
CN111297341A (en) * 2020-02-20 2020-06-19 京东方科技集团股份有限公司 Developments blood pressure check out test set and pulse wave characteristic extraction equipment
CN111631733A (en) * 2020-06-19 2020-09-08 浙江澍源智能技术有限公司 Arterial blood spectrum detection method and device
CN111631733B (en) * 2020-06-19 2024-01-26 浙江澍源智能技术有限公司 Arterial blood spectrum detection method and device
CN111904404A (en) * 2020-08-06 2020-11-10 苏州国科医工科技发展(集团)有限公司 Blood pressure regulation and control equipment containing closed-loop monitoring and based on ear vagus nerve stimulation
CN111973165A (en) * 2020-08-14 2020-11-24 北京航空航天大学 Linear and nonlinear mixed non-invasive continuous blood pressure measuring system based on PPG
CN112656393A (en) * 2020-12-08 2021-04-16 山东中科先进技术研究院有限公司 Heart rate variability detection method and system
CN112656393B (en) * 2020-12-08 2023-06-27 山东中科先进技术有限公司 Method and system for detecting heart rate variability
CN113208573A (en) * 2021-04-21 2021-08-06 北京雪扬科技有限公司 Support wearable equipment of PPG + ECG function
CN113208573B (en) * 2021-04-21 2022-07-12 北京雪扬科技有限公司 Support wearable equipment of PPG + ECG function
CN113520350A (en) * 2021-07-27 2021-10-22 香港心脑血管健康工程研究中心有限公司 Processing method and device for acquiring relevant characteristic parameters and index information of blood pressure map signals
CN113520350B (en) * 2021-07-27 2024-05-17 香港心脑血管健康工程研究中心有限公司 Processing method and device for acquiring relevant characteristic parameters and index information of blood pressure chart signals
CN113855007A (en) * 2021-08-27 2021-12-31 联卫医疗科技(上海)有限公司 Method and device for obtaining machine learning model samples for blood glucose prediction
CN113854985A (en) * 2021-08-27 2021-12-31 联卫医疗科技(上海)有限公司 Method and device for obtaining machine learning model samples for blood pressure prediction
CN114145724A (en) * 2021-12-08 2022-03-08 四川北易信息技术有限公司 Method for dynamically monitoring blood pressure based on ECG (electrocardiogram) and PPG (photoplethysmography) multiple physiological characteristic parameters
CN114305364A (en) * 2022-01-05 2022-04-12 北京科技大学 Blood pressure detection method, system and equipment based on millimeter wave radar
CN114305364B (en) * 2022-01-05 2023-08-11 北京科技大学 Blood pressure detection method, system and equipment based on millimeter wave radar
CN115670408A (en) * 2022-12-28 2023-02-03 可孚医疗科技股份有限公司 Blood pressure measuring device, linear model coefficient self-correction method and system thereof, and measuring method

Also Published As

Publication number Publication date
CN106413534B (en) 2020-06-30
WO2017024457A1 (en) 2017-02-16

Similar Documents

Publication Publication Date Title
CN106413534A (en) Blood-pressure continuous-measurement device, measurement model establishment method, and system
CN102429649B (en) Continuous blood pressure measuring device
JP6721153B2 (en) Biological information analysis device, system, and program
CN104873186B (en) A kind of wearable artery detection device and its data processing method
CN102397064B (en) Continuous blood pressure measuring device
CN102488503B (en) Continuous blood pressure measurer
WO2019161609A1 (en) Method for analyzing multi-parameter monitoring data and multi-parameter monitor
CN108577830B (en) User-oriented physical sign information dynamic monitoring method and system
US8755871B2 (en) Systems and methods for detecting arrhythmia from a physiological signal
WO2019161608A1 (en) Multi-parameter monitoring data analysis method and multi-parameter monitoring system
WO2020007041A1 (en) Blood pressure measurement apparatus and method, electronic device, and computer readable storage medium
CN110251105A (en) A kind of non-invasive blood pressure measuring method, device, equipment and system
US20150088003A1 (en) Data integrity
JP2018508273A (en) Method and apparatus for deriving mean arterial pressure in a subject
US11517223B2 (en) Apparatus and method for estimating blood glucose
CN204708828U (en) A kind of wearable noinvasive arterial health checkout gear
CN110236508A (en) A kind of non-invasive blood pressure continuous monitoring method
KR20200021208A (en) Apparatus and method for estimating blood pressure
WO2021184803A1 (en) Blood pressure classification method and apparatus
CN111248881B (en) Method and device for predicting blood pressure
CN104739394A (en) Portable human body physiological signal monitoring and alarming system
CN108124419B (en) Blood pressure measuring method, intelligent heart rate earphone and system
WO2020160058A1 (en) Systems and methods for computationally efficient non-invasive blood quality measurement
CN113040738B (en) Blood pressure detecting device
KR20220053439A (en) Continuous blood pressure monitoring system based on photoplethysmography by using convolutionalbidirectional long short-term memory neural networks

Legal Events

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