CN104644151B - A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal - Google Patents

A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal Download PDF

Info

Publication number
CN104644151B
CN104644151B CN201510051739.XA CN201510051739A CN104644151B CN 104644151 B CN104644151 B CN 104644151B CN 201510051739 A CN201510051739 A CN 201510051739A CN 104644151 B CN104644151 B CN 104644151B
Authority
CN
China
Prior art keywords
waveform
wave
pressure
pressure pulse
pulse signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510051739.XA
Other languages
Chinese (zh)
Other versions
CN104644151A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201510051739.XA priority Critical patent/CN104644151B/en
Publication of CN104644151A publication Critical patent/CN104644151A/en
Application granted granted Critical
Publication of CN104644151B publication Critical patent/CN104644151B/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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy

Landscapes

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

Abstract

A kind of pressure pulse wave wave travel prediction meanss based on photoelectricity volume pulse signal, it is characterised in that:Including Waveform Input module, waveform conditioning module, waveform fitting module, waveform modular converter, and waveform output module.Wherein:Waveform conditioning module includes pretreatment circuit, it is single to clap separator and normalization circuit, waveform fitting module includes fitting function setting apparatus, waveform fitting device and waveform quality arbiter, waveform modular converter includes position setting apparatus, object function setting apparatus, characterizing population group's setting apparatus, parameter converter and waveform synthesizer.The device can have partes corporis humani position pressure pulse wave signal using partes corporis humani's position photoelectricity volume pulse signal according to the prediction of physiology statistical law, and the use scope of the device, prediction effect and stability all improve compared with existing apparatus.

Description

A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal
Technical field
It is more particularly to a kind of that each position pressure is predicted based on photoplethysmographic the present invention relates to technical field of medical equipment The method of power pulse waveform.
Background technology:
Contain abundant haemodynamics information in pulse wave.Always as clinical diagnosis and the foundation for the treatment of, largely Clinical measured result confirm that the feature of pulse wave has close relationship with cardiovascular physiology state.Pulse wave is shown The synthesis of the aspects such as the form (shape of ripple), intensity (amplitude of ripple), speed (speed of ripple) and the rhythm and pace of moving things (cycle of ripple) that come Information reflects many physiology and pathological characters of cardiovascular system of human body to a great extent.
In terms of acquisition principle, the acquisition mode of current pulse wave mainly has to be adopted using the pressure pulse wave of pressure sensor Collection or the photoplethysmographic collection using photoelectric sensor.Intraarterial pressure pulse wave has obtained more comprehensive analysis and has ground Study carefully, fluid mechanic model is more clear and definite, corresponding Hydrodynamic character and cardiovascular physiology meaning application are also relatively broad.So And, the collection of human pressure's pulse wave is highly prone to many-sided interference such as collection position, and required operative skill requirement is higher, and lacks Weary repeatability, is not easy to continuous detection.The collection of photoplethysmographic finger tip has preferable stability and repeatability, but its More perfect pressure pulse wave is studied to the accuracy that cardiovascular function judges not good enough.
The content of the invention:
Existing technical scheme is mainly using gathering position photoelectricity plethysmogram signal and correspondence target site pressure pulse wave is believed Number the non-physiological parameter model such as transmission function set up of power spectrumanalysis, and focus mainly on finger volume pulse signal and oar Arterial pressure pulse signal.Because the individual difference of pulse signal is larger, said apparatus in large-scale application stability compared with Difference, function is more single, and set up model is unfavorable for the amendment and improvement to model in itself without physiological significance.
To solve the above problems, the present invention propose respectively a kind of volume pulsation signal based on pulse wave physiological characteristic and The waveform of pressure pulse signal expression formula containing ginseng, and set up collection position photoplethysmographic expression using priori statistical law Regression equation group between formula parameter and corresponding target site pressure pulse waveform expression argument.So as to the photoelectricity for realizing being input into Volume pulsation signal is converted to the pressure pulse wave signal of target site.The method of waveform fitting has preferable stability, returns Returning equation group has clearer and more definite physiological significance, convenient to be finely adjusted and improvement for different physiological status.It is possible thereby to solve Prediction stability of waveform is poor in the prior art, and the problem that model cannot be corrected.
To reach above-mentioned purpose, the technical solution used in the present invention is:A kind of pressure based on photoelectricity volume pulse signal Power pulse wave propagation forecast device, it is characterised in that:Including signal input conditioning module, waveform fitting module, waveform modulus of conversion Block, Waveform composition module and output module.
The Waveform Input module is received from the time domain photoelectricity volume pulse signal of a certain position actual measurement of human body.
The Signal-regulated kinase, the time domain photoelectricity volume pulse signal to being input into is pre-processed, using prior art It is broken down into the single of correspondence single cardiac cycle and claps signal, place is normalized with wavelength to each single amplitude for clapping signal Reason.
The waveform fitting module receives the single bat pulse signal after normalization, and using given waveform expression formula containing ginseng fIApplication curves fitting algorithm is fitted to it, and the expression formula containing ginseng of waveform by representing the main ripple of pulse waveform, dicrotic pulse respectively Ripple, the expression formula containing ginseng of back wave are added and determine, fitting gained analytical expression parameters vector is II, as wave character Parameter vector.
Photoelectricity volume pulse signal containing ginseng expression formula be:
Wherein, Hin, bin, WinIt is parameter, t is independent variable, represents sampling number.N=1 in expression formula, 2,3 part point Not Dui Ying pulse waveform main ripple, back wave, dicrotic pulse waveform.Wherein, HinRepresent the amplitude of fluctuation, binRepresent fluctuation Center, WinThe width of fluctuation.
Curve fitting process algorithm uses least-squares algorithm, and parameters scope is entered according to each fluctuation physiological significance Row is limited, afterwards setting fitting primary condition Hi1>Hi2>Hi3, bi1<bi2<bi3,Win>0.It is fitted the characteristic parameter vector I for determiningI For:
II=[Hi1,Hi2,Hi3,bi1,bi2,bi3,Wi1,Wi2,Wi3]
Under the conditions of above-mentioned waveform expression formula containing ginseng, fitting effect depends primarily on the interference journey that pulse wave collection is subject to Degree, therefore determine coefficients R to be fitted2As the quantification standard that waveform quality differentiates.Fitting determines coefficients R2It is conventional judgement two The computational methods of curve similarity degree.R in the present invention2Equally as acquisition quality discriminant parameter, correspondence R2Less than certain value Single waveform of clapping is thought acquisition quality difference and is given up.
R2Computing formula is as follows:
Wherein, fi,Represent that measured light Power Capacity Pulse Rate strong point, measured light Power Capacity pulse data are average respectively Value and data point fitting desired value, pl is unicast data point number.
The waveform modular converter, is grouped according to measured's sex, age, mean arterial pressure index, and according to right The priori system between lower photoplethysmographic actual measurement position and required forecast pressure pulse wave position waveform feature parameter should be grouped Meter rule, prediction corresponding position pressure pulse wave waveform feature parameter is calculated using measured light Power Capacity pulse wave characteristic parameters.
Above-mentioned correspondence classification lower photoplethysmographic actual measurement position and required forecast pressure pulse wave position wave character Priori statistical law method for building up between parameter is as follows:
(1) first, the crowd for participating in experiment is grouped according to sex, age, mean arterial pressure (MAP), its middle age Age was starting with 20 years old, was within 5 years old interval.Mean arterial pressure is starting with 70mmHg, and 10mmHg is interval.People to participating in experiment Group is grouped.Detect ear, finger end, photoplethysmographic at toe end simultaneously to above-mentioned each group Subject Population respectively, And using pressure pulse wave signal at pressure sensor detection radial artery, arteria brachialis, arteria carotis.So as to obtain different crowd feature Observed pressure pulse involve photoplethysmographic signal.
(2) set up after between actual measurement photoplethysmographic and actual measurement target site pressure pulse wave waveform feature parameter Statistics relation.It is similar to above-mentioned photoplethysmographic feature extraction mode, to extract pressure pulse wave waveform feature parameter, Observed pressure pulse waveform is fitted using pressure pulse wave expression formula, pressure pulse wave expression formula containing ginseng is:
Expression argument and domain of definition are and fIIt is identical.Fit procedure algorithm uses least-squares algorithm, and according to each ripple Dynamic physiological significance is defined to parameters scope, afterwards setting fitting primary condition Ho1>Ho2>Ho3, bo1<bo2<bo3,Won> 0。
Its characteristic parameter vector is IO=[Ho1,Ho2,Ho3,bo1,bo2,bo3,Wo1,Wo2,Wo3]
(3) measured waveform to different grouping different parts is utilized respectively fIAnd fOTo the photoplethysmographic ripple surveyed Shape and pressure pulse waveform are fitted using waveform fitting module, obtain the I of correspondence each photoelectricity volume pulse signal of actual measurementI And the I of pressure pulse signalOVector.For each site pressure pulse signal IOVectorial each parameter, sets up correspondence different grouping While gather photoelectricity volume pulse signal characteristic parameter vector multiple linear regression equations, i.e.,:
Ho1=TM11×Hi1+TM12×Hi2+......+TM19×Wi3+CM1
Ho2=TM21×Hi1+TM22×Hi2+......+TM29×Wi3+CM2
Ho3=TM31×Hi1+TM32×Hi2+......+TM39×Wi3+CM3
bo1=TM41×Hi1+TM42×Hi2+......+TM49×Wi3+CM4
bo2=TM51×Hi1+TM52×Hi2+......+TM59×Wi3+CM5
bo3=TM61×Hi1+TM62×Hi2+......+TM69×Wi3+CM6
Wo1=TM71×Hi1+TM72×Hi2+......+TM79×Wi3+CM7
Wo2=TM81×Hi1+TM82×Hi2+......+TM89×Wi3+CM8
Wo3=TM91×Hi1+TM92×Hi2+......+TM99×Wi3+CM9
Each term coefficients of TM and CM are determined by multiple linear regression equations.Arrange IOParameters set up equation group, arrange system Matrix number and constant matrices, can obtain TM and CM, wherein, TM is 9 rank square formations, and CM is 9 element column vectors.
It is special using the waveform for gathering position in application process according to above-mentioned utilization prior probability gained TM and CM matrixes Levy parameter vector IICalculate the characteristic parameter vector I of target siteO.According to surveyed photoplethysmographic, its feature I is extractedI, Target site characteristic vector I can then be obtainedOFor:
IO=TM × II+CM
The target site characteristic parameter vector I that will be obtained afterwardsOSubstitute into pressure pulse wave expression formula containing ginseng of corresponding position fO, obtain corresponding pressure pulse wave analytical expression.Complete waveform conversion.
The waveform output module, by result and target site characteristic parameter the vector I of above-mentioned Waveform composition moduleOBy will Form is asked to export.
The beneficial effects of the present invention are:
(1) apparatus of the present invention only need to detect the single position photoelectricity volume pulse signal of human body, you can obtained under certain precision Take the pressure pulse wave signal of different parts.It is simple to operate, and waveform quality stabilization.Normal pressures pulse wave is overcome to adopt Testing process is complicated during collection, using inconvenience, and is difficult to obtain the problem of stabilization pulse wave.In actual application The direct times of collection of pressure pulse wave is effectively reduced, the demand to operative skill is reduced, measured's level of comfort is improve.Through Substantial amounts of confirmatory experiment is crossed, is worked well.
(2) the waveform fitting device based on pulse wave physiological fluctuating characteristic is applied to multiple location pulse wave first, So as to analyze the change of the physiological in its communication process.The adjustment and correction of model are facilitated, is that larger scale clinical application is established Certain basis.
Brief description of the drawings:
Fig. 1 is structured flowchart of the present invention
Fig. 2 is operational flowchart of the present invention
Fig. 3 is waveform fitting and characteristic parameter schematic diagram
Fig. 4 is an actual measurement example, and gives the contrast with observed pressure waveform
Fig. 5 is that the present invention is predicting the experiment effect figure of radial artery pressure pulse wave using finger tip photoplethysmographic. R2Data are that observed pressure pulse list claps waveform and the single cross validation effect clapped between waveform of prediction, and 426 are tested altogether.
Specific embodiment:
A kind of more typical specific embodiment of the invention is described in detail below in conjunction with accompanying drawing.
A kind of typical application scenarios of the invention are to predict radial artery using human body finger tip photoplethysmographic signal Pressure pulse wave signal.It is possible thereby to the mature technology gathered using finger tip and high-quality waveform obtain physiological significance definitely Radial artery pressure pulse wave data.
As shown in figure 1, after prediction process starts, first according to the Sex, Age of measured, blood pressure selects corresponding feature Crowd, and corresponding prediction position is selected as needed.With 56 years old age women, pressure value is tested for 90/130mmHg's As a example by person, sex, age and pressure value are set in step s 201, collection position is finger tip, and target prediction position is radial artery. According to setting in S202, system is the packet situation according to priori statistical law, automatically selects out correspondence crowd and position TM and CM matrixes.
System starts to receive measured waveform in step S203, and step S204 nursed one's health input waveform, filter baseline with Hz noise, by continuous wave according to cardiac cycle separated component vertical wave shape, and carries out wave-shape amplitude and ripple in units of waveform Normalization long, wherein, the normalization of wavelength is realized by the method for interpolation.
In S205 steps, according to the collection waveform fitting expression formula given in S203, to each single waveform of clapping using most Small least square method carries out curve fitting.
Fitting expression is:
Setting primary condition.So as to obtain each single I for clapping waveform of correspondenceIAnd calculate R2
In this example, waveforms amplitude and wavelength are respectively defined as 100 units, then
II=[Hi1,Hi2,Hi3,bi1,bi2,bi3,Wi1,Wi2,Wi3]=[43,69,49,15,27,48,14,25,52]
R2>0.99
According to R in step S2062Numerical value carries out waveform quality judgement, to R2<0.95 waveform is given up, and records house Abandon ratio.If waveform occur to give up, next bat waveform is extracted again and is analyzed.R under normal circumstances2<0.95 waveform ratio< 3%, when giving up ratio more than 5%, it is considered as adjusting acquisition mode.
Is of the step S207 to up-to-standard waveform in S206IChanged, calculated the characteristic vector I of target site waveformO, According to formula:
IO=TM × II+CM
In this example, TM and CM is by the radial artery pressure pulse wave and finger tip photoelectricity volume measured before the patient simultaneously Pulse wave signal is respectively according to fiWith foFitting, and corresponding regression equation acquisition is set up by the parameter for obtaining respectively.
It is computed obtaining:
IO=[Ho1,Ho2,Ho3,bo1,bo2,bo3,Wo1,Wo2,Wo3]=[64,71,35,14,25,49,13,22,48]
By I in step S208OMiddle parameters substitute into given target site pressure pulse wave expression formula, i.e.,
Target site prediction waveform expression formula and respectively corresponding main ripple are obtained, reflection involves dicrotic wave prediction waveform.
Step S209 exports above-mentioned waveform and parameter according to specified format.

Claims (1)

1. a kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal, it is characterised in that:Including ripple Shape input module, Signal-regulated kinase, waveform fitting module, waveform modular converter, and waveform output module;
The Waveform Input module is received from the time domain photoelectricity volume pulse signal of a certain position actual measurement of human body;
The Signal-regulated kinase, the time domain photoelectricity volume pulse signal to being input into is pre-processed, and is broken down into correspondence single The single of one cardiac cycle claps pulse signal, and each single amplitude for clapping pulse signal is normalized with wavelength;
The waveform fitting module receives the single bat pulse signal after normalization, and using given waveform expression formula containing ginseng fIShould It is fitted with curve fitting algorithm, the expression formula containing ginseng of waveform by representing the main ripple of pulse waveform, dicrotic wave, anti-respectively The expression formula containing ginseng of ejected wave is added and determines, fitting gained analytical expression parameters vector is II, as waveform feature parameter Vector;
The waveform of photoelectricity volume pulse signal contains ginseng expression formula:
Wherein, Hin, bin, WinIt is parameter, t is independent variable, represents sampling number;N=1 in expression formula, 2,3 part is right respectively Answer main ripple, back wave, the dicrotic pulse waveform of pulse waveform;Wherein, HinRepresent the amplitude of fluctuation, binRepresent the center of fluctuation Position, WinThe width of fluctuation;
Curve fitting process algorithm uses least-squares algorithm, and parameters scope is limited according to each fluctuation physiological significance Fixed, setting afterwards is fitted primary condition Hi1>Hi2>Hi3, bi1<bi2<bi3,Win>0;It is fitted the characteristic parameter vector I for determiningIFor:
II=[Hi1,Hi2,Hi3,bi1,bi2,bi3,Wi1,Wi2,Wi3]
Under the conditions of above-mentioned waveform expression formula containing ginseng, fitting effect depends primarily on the annoyance level that pulse wave collection is subject to, therefore Determine coefficients R to be fitted2As the quantification standard that waveform quality differentiates;Fitting determines coefficients R2It is the conventional curve of judgement two The computational methods of similarity degree;R2As acquisition quality discriminant parameter, correspondence R2Single waveform of clapping less than certain value is thought to gather matter Amount difference is simultaneously given up;
R2Computing formula is as follows:
Wherein,Measured light Power Capacity Pulse Rate strong point, measured light Power Capacity pulse data average value and number are represented respectively Strong point is fitted desired value, and pl claps pulse signal data point number for single;
The waveform modular converter, is grouped according to measured's sex, age, mean arterial pressure index, and according to correspondence point Priori statistics rule between the lower photoplethysmographic actual measurement position of group and required forecast pressure pulse wave position waveform feature parameter Rule, prediction corresponding position pressure pulse wave waveform feature parameter is calculated using measured light Power Capacity pulse wave characteristic parameters;
Above-mentioned correspondence packet lower photoplethysmographic actual measurement position and required forecast pressure pulse wave position waveform feature parameter Between priori statistical law method for building up it is as follows:
(1) first, the crowd for participating in experiment is grouped according to sex, age, mean arterial pressure, the wherein age was with 20 years old Starting, 5 years old is interval;Mean arterial pressure is starting with 70mmHg, and 10mmHg is interval;Crowd to participating in experiment is divided Group;Detect ear, finger end, photoplethysmographic at toe end simultaneously to above-mentioned each group Subject Population respectively, and using pressure Pressure pulse wave signal at force snesor detection radial artery, arteria brachialis, arteria carotis;So as to obtain the actual measurement pressure of different crowd feature Power pulse involves photoplethysmographic signal;
(2) statistics between actual measurement photoplethysmographic and actual measurement target site pressure pulse wave waveform feature parameter is set up after Relation;It is similar to above-mentioned photoplethysmographic waveform feature parameter vector extracting mode, to extract pressure pulse waveform Characteristic parameter, is fitted, pressure pulse wave using pressure pulse waveform expression formula containing ginseng to observed pressure pulse waveform Waveform contains ginseng expression formula:
Expression argument and domain of definition are and fIIt is identical;Fit procedure algorithm uses least-squares algorithm, and according to each fluctuation physiology Meaning is defined to parameters scope, afterwards setting fitting primary condition Ho1>Ho2>Ho3, bo1<bo2<bo3,Won>0;
Its characteristic parameter vector is IO=[Ho1,Ho2,Ho3,bo1,bo2,bo3,Wo1,Wo2,Wo3]
(3) measured waveform to different grouping different parts is utilized respectively fIAnd fOTo survey photoplethysmographic waveform and Pressure pulse waveform is fitted using waveform fitting module, obtains the I of correspondence each photoelectricity volume pulse signal of actual measurementIAnd pressure The I of power pulse signalOVector;For each site pressure pulse signal IOVectorial each parameter, sets up the same of correspondence different grouping When gather photoelectricity volume pulse signal characteristic parameter vector multiple linear regression equations, i.e.,:
Ho1=TM11×Hi1+TM12×Hi2+......+TM19×Wi3+CM1
Ho2=TM21×Hi1+TM22×Hi2+......+TM29×Wi3+CM2
Ho3=TM31×Hi1+TM32×Hi2+......+TM39×Wi3+CM3
bo1=TM41×Hi1+TM42×Hi2+......+TM49×Wi3+CM4
bo2=TM51×Hi1+TM52×Hi2+......+TM59×Wi3+CM5
bo3=TM61×Hi1+TM62×Hi2+......+TM69×Wi3+CM6
Wo1=TM71×Hi1+TM72×Hi2+......+TM79×Wi3+CM7
Wo2=TM81×Hi1+TM82×Hi2+......+TM89×Wi3+CM8
Wo3=TM91×Hi1+TM92×Hi2+......+TM99×Wi3+CM9
Each term coefficients of TM and CM are determined by multiple linear regression equations;Arrange IOParameters set up equation group, arrange coefficient matrix And constant matrices, TM and CM is obtained, wherein, TM is 9 rank square formations, and CM is 9 element column vectors;
According to above-mentioned utilization prior probability gained TM and CM matrixes, in application process, joined using the wave character for gathering position Number vector IICalculate the characteristic parameter vector I of target siteO;According to surveyed photoplethysmographic, its feature I is extractedI, then Target site characteristic parameter vector IOFor:
IO=TM × II+CM
The target site characteristic parameter vector I that will be obtained afterwardsOSubstitute into pressure pulse waveform expression formula containing ginseng of corresponding position fO, obtain corresponding pressure pulse wave analytical expression;Complete waveform conversion;
The waveform output module, by result and target site characteristic parameter the vector I of above-mentioned waveform modular converterOShape on request Formula is exported.
CN201510051739.XA 2015-02-01 2015-02-01 A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal Active CN104644151B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510051739.XA CN104644151B (en) 2015-02-01 2015-02-01 A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510051739.XA CN104644151B (en) 2015-02-01 2015-02-01 A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal

Publications (2)

Publication Number Publication Date
CN104644151A CN104644151A (en) 2015-05-27
CN104644151B true CN104644151B (en) 2017-07-07

Family

ID=53236165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510051739.XA Active CN104644151B (en) 2015-02-01 2015-02-01 A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal

Country Status (1)

Country Link
CN (1) CN104644151B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105615845B (en) * 2016-02-25 2020-05-19 广州视源电子科技股份有限公司 Method and system for detecting interference pulse signal
CN109662702A (en) * 2018-05-23 2019-04-23 李芝宏 Condenser type pulse detection system and method
CN111685749B (en) * 2020-06-18 2022-09-02 郑昕 Construction method of pulse pressure wave mathematical model
CN116981398A (en) * 2020-12-04 2023-10-31 华为技术有限公司 Blood pressure prediction method, blood pressure prediction device, and computer program
CN112826459B (en) * 2021-01-08 2022-11-29 北京工业大学 Pulse wave waveform reconstruction method and system based on convolution self-encoder

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS586691A (en) * 1981-07-06 1983-01-14 Hitachi Denshi Ltd Simultaneous monitoring system for plural waveform
CN1121798A (en) * 1994-08-16 1996-05-08 北京工业大学 Cardiovascular function dynamic parameter testing analysis method and apparatus
CN102894982A (en) * 2012-09-28 2013-01-30 北京工业大学 Non-invasive detecting method for blood viscosity based on pulse wave
CN104138253A (en) * 2013-05-11 2014-11-12 吴健康 Noninvasive continuous arterial blood pressure measuring method and equipment

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070287923A1 (en) * 2006-05-15 2007-12-13 Charles Adkins Wrist plethysmograph
US20150031971A1 (en) * 2013-07-26 2015-01-29 Covidien Lp Methods and systems for using an estimate signal to improve signal resolution in a physiological monitor

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS586691A (en) * 1981-07-06 1983-01-14 Hitachi Denshi Ltd Simultaneous monitoring system for plural waveform
CN1121798A (en) * 1994-08-16 1996-05-08 北京工业大学 Cardiovascular function dynamic parameter testing analysis method and apparatus
CN102894982A (en) * 2012-09-28 2013-01-30 北京工业大学 Non-invasive detecting method for blood viscosity based on pulse wave
CN104138253A (en) * 2013-05-11 2014-11-12 吴健康 Noninvasive continuous arterial blood pressure measuring method and equipment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《基于光电容积脉搏波描记法的无创连续血压测量》;李章俊 等;《中国生物医学工程学报》;20120831;第31卷(第4期);607-614 *

Also Published As

Publication number Publication date
CN104644151A (en) 2015-05-27

Similar Documents

Publication Publication Date Title
CN104644151B (en) A kind of pressure pulse wave wave travel Forecasting Methodology based on photoelectricity volume pulse signal
CN102488503B (en) Continuous blood pressure measurer
US10825569B2 (en) Universal non-invasive blood glucose estimation method based on time series analysis
CN102429649B (en) Continuous blood pressure measuring device
CN102397064B (en) Continuous blood pressure measuring device
CN110251105A (en) A kind of non-invasive blood pressure measuring method, device, equipment and system
CN105943005B (en) The non-invasive blood pressure detection device mixed based on photoelectricity green light pulse with electrocardiogram
WO2017024457A1 (en) Blood-pressure continuous-measurement device, measurement model establishment method, and system
CN107736880A (en) A kind of pulse analysis method and system
CN110037668B (en) System for judging age, health state and malignant arrhythmia identification by combining pulse signal time-space domain with model
CN103876723A (en) Method of obtaining blood pressure value by noninvasive radial artery wave calculating pulse wave transmission time
CN109833034A (en) The method and device of blood pressure data is extracted in a kind of pulse wave signal
CN101703396A (en) Radial artery pulse wave based cardiovascular function parameter detection and analysis method and detection device
CN106419937A (en) Mental stress analysis system based on heart sound HRV theory
Chen et al. Machine learning method for continuous noninvasive blood pressure detection based on random forest
WO2023226824A1 (en) Blood pressure detection apparatus based on arteriolar photoplethysmography
CN113160921A (en) Construction method and application of digital human cardiovascular system based on hemodynamics
CN109872820A (en) A kind of no cuff blood pressure measuring method, device, equipment and storage medium
CN113520333A (en) Method, device and equipment for determining core body temperature and readable medium
CN1275179C (en) Computer evaluating method for human body sub-health status
CN107995981B (en) Data processing method for blood pressure measuring device
CN102894982A (en) Non-invasive detecting method for blood viscosity based on pulse wave
CN113040738B (en) Blood pressure detecting device
CN110200642A (en) A kind of measurement method and terminal of cognitive load and psychological pressure
CN113069091A (en) Pulse condition classification device and method for PPG (photoplethysmography) signals

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