CN105726000B - A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse - Google Patents
A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse Download PDFInfo
- Publication number
- CN105726000B CN105726000B CN201610064419.2A CN201610064419A CN105726000B CN 105726000 B CN105726000 B CN 105726000B CN 201610064419 A CN201610064419 A CN 201610064419A CN 105726000 B CN105726000 B CN 105726000B
- Authority
- CN
- China
- Prior art keywords
- blood
- limbs
- pressure
- pulse wave
- blood flow
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording 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/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
- A61B5/02035—Determining blood viscosity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording 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/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
A kind of calculation method of the cardiovascular functional parameter based on blood pressure of four limbs pulse, belongs to biomethanics technical field.Using the arterial pressure and pressure pulse wave of pulse wave sensor and blood pressure sensor detection four limbs, cardiac function, the noninvasive synchronous detection of rete arteriosum hardenability, heart and the blood flow parameter of vascular function are realized.Its step is first to acquire the blood pressure and pulse wave information of four limbs arteria brachialis and ankle artery, carries out wave-average filtering to a series of pulse waves of acquisition, obtains the average waveform of four limbs arteria brachialis and ankle artery respectively.Based on wave character parameter and non-linear pulse wave theory that waveform analysis obtains, the relationship between four limbs pulse wave and barometric gradient, blood flow is established, then in blood flow and four limbs pulse blood pressure data Information base, calculates and obtains cardiovascular functional parameter.The present invention will provide regular one kind, hurtless measure, easy to operate, fast output, cheap detection means for cardiovascular system.
Description
Technical field
The invention belongs to biomethanics technologies, are related to computer program, circuit design etc., can be used for cardiovascular function ginseng
Several noninvasive calculating.
Background technique
Cardiovascular disease has become the adult health in the harm whole world and China and leads to dead primary illness, but painstaking effort
Pipe disease early clinic symptom is unobvious, until complication deterioration can just diagnose cardiovascular disease.Currently, cardiovascular disease
The diagnostic techniques of disease is broadly divided into direct method and two kinds of indirect method.The direct method of measurement can provide that precision is higher, recordable company
Continuous data, but the acquisition of these supplemental characteristics is obtained by invasive detection means, not only there is certain wind for this method
It is dangerous, and the technology of the clinical level or Medical Devices to doctor has when high requirement.The indirect method of measurement is clinically answered
With it is more be to be diagnosed by cardiovascular instrument, such as Echocardiography, impedance cardiography, spiral CT, nuclear medicine, magnetic resonance
The methods of (MRI)；Although Non-invasive detection, detection environmental requirement is stringent, the disadvantages of cannot carrying, and those who are investigated are only
It can be detected in special medical space, by the healthcare givers of profession.Although the detection technique of current cardiovascular disease exists
The prediction technique applied and provided clinically is had been combined, but still there are many high-risk patients undiscovered.
Domestic and foreign scholars have been devoted to study the Non-invasive detection technology of cardiovascular system based on pulse wave, have built up
Devise Non-invasive detection model, method and apparatus, but there are still some problems, as: evaluation parameter is single, involved in model
The problems such as more hypothesis and conjecture and product technology barrier.So the non-invasive detection methods of cardiovascular system are up for into one
Step is improved and optimization.
Further, since the health perception of modern people greatly improves, for above-mentioned heart function detection method is at high cost, behaviour
Make complicated and be only applicable to special medical institutions, there is also very big differences for the general inspection of health and family health care for being centrifuged vascular diseases
Away from.The invention proposes the cardiovascular noninvasive calculation method based on blood pressure of four limbs pulse, this method will be cardiovascular function
It can detect and regular one kind, hurtless measure, easy to operate, fast output, cheap detection method and means are provided；It is simultaneously also the heart
The non-invasive detection methods of dirty vascular system are provided fundamental basis.
Summary of the invention
The object of the present invention is to provide a kind of calculating of cardiovascular system functional parameter based on blood pressure of four limbs pulse wave
Method, this method can detect 8 tunnel information of the four limbs arteria brachialis of human body and the blood pressure and pulse wave of ankle artery with noninvasive synchronization, according to
Cardiac function, rete arteriosum hardenability, hemorheology according to 8 tunnel information detection human body synchronous with non-linear pulse wave theory refer to
The state parameter of the cardiovasculars functions such as mark.It is intended that improving the accuracy that Non-invasive detection calculates cardiac functional parameter, avoid
The problems such as invasive and high cost of existing invasive detection means.
The four limbs arteria brachialis and ankle arterial pulse wave need to carry out wave-average filtering, and average is that four limbs pressure sensor mentions
The step of pressure pulse wave component taken, wave-average filtering, comprising: reading Wave data first demarcates complete waveform and finds out week
Then phase carries out period and amplitude normalized to waveform, pulse wave unicast data is made to have unified amplitude and length, and
Acquire average waveform.
The non-linear pulse wave theory establishes blood pressure of four limbs, blood flow two basic haemodynamics parameters and four
Relational expression between limb pulse wave, and according to the inherent feature parameter or geometric parameter of human heart vascular system, it derives
Blood pressure of four limbs pulse data, the relational expression of blood flow and cardiovascular functional parameter are simultaneously calculated.
The calculation method mainly utilizes non-linear pulse wave theory and blood pressure of four limbs pulse wave data, establishes pressure
Gradient and blood flow, the relational model of blood flow and cardiovascular functional parameter；Go out a week aroused in interest using numerical analysis
The numerical solution of the four limbs blood flow of different moments, the recurrence of four limbs blood flow is calculated using statistical homing method respectively in phase
Coefficient, amendment obtain human bloodstream amount and mean blood flow；In blood flow and four limbs pulse blood pressure data Information base, calculate
Solve cardiovascular functional parameter parameter.Calculated result will be carried out by error analysis judging whether to reach convergence in algorithm,
If convergence, exports result；Otherwise calculating will be iterated.
The detection device of the 8 tunnel information, including pulse wave sensor, blood pressure sensor, A/D conversion circuit, interface
Circuit and accessory circuit etc. pass through usb hub to the synchronous 8 circuit-switched data information of blood pressure of four limbs and pulse wave for detecting and acquiring
Be transferred to the end PC, by the process control slave computer of host computer, and the blood pressure waveform data detected are handled, store and
Output.The hardware components, it is easy to operate, it easily detects, at low cost, the used time is short.
The detecting state of the 8 tunnel information, including after quiescent condition, quantitative exercise and movement again under quiescent condition
Blood pressure and pulse wave signal.It is intended that indicating the physiological status at individual a certain moment with static parameter, dynamic parameter list is used
The individual self-control made is reacted after showing itself regulatory function and movement when Individual Experience movement.
A kind of calculation method of the cardiovascular functional parameter based on blood pressure of four limbs pulse can be the cardiovascular of human body
System provides regular one kind, hurtless measure, easy to operate, fast output, cheap detection means；It can be adapted for community doctor
Treatment, rehabilitation, sport and body-building, leisure such as are recuperated at the general inspections of health of health care institutions.
It is a feature of the present invention that successively realized according to the following steps in upper computer and lower computer:
Step (1), upper computer and lower computer initialization:
Host computer input: the letter such as personal information of measured, including gender, age, height, weight, shoulder breadth medical history
Breath.
The next machine testing: the pulse wave and blood pressure of measured's arteria brachialis and ankle artery four limbs.
Host computer record: the pulse waveform and blood pressure pressure value of measured's arteria brachialis and ankle artery four limbs.
Step (2), host computer detect automatically, record, saving the pulse waveform figure of four limbs, utilize to obtained waveform diagram
Period and amplitude normalized carry out wave-average filtering, are then calculated as follows and respectively obtain measured in a cardiac cycle
Four limbs mean arterial pressure Pmf, mean arterial pressure is the average value of the corresponding arterial pressure of four limbs in a cardiac cycle.
T be a cardiac cycle, f be four limbs label, f=1,2,3,4 respectively indicate right upper extremity, left upper extremity, right lower extremity,
Left lower extremity, P (t)fFor the corresponding pulse pressure value of different moments each sampled point in a cardiac cycle, (f=1,2,3,4).
Step (3), is calculated as follows shape factor kf, (f=1,2,3,4).
PdfFor diastolic pressure, PsfIt is the actual measured value of four limbs arteria brachialis and ankle artery for systolic pressure, (f=1,2,3,
4)。
The spread speed c (t) of arteria brachialis and ankle arterial pressure pulse wave is calculated as follows in step (4)f:
Wherein:
P(t)fFor the corresponding pulse pressure value of different moments each sampled point in a cardiac cycle, (f=1,2,3,4)；
ρ is density of blood, ρ=1.05*10-3kg/cm3；
bfIt for pulse waveform parameter, is calculated as follows, (f=1,2,3,4):
A is parameter related with velocity profile, a=0.57；
The space pressure gradient F (t) of arteria brachialis and ankle arterial pressure pulse wave is calculated as follows in step (5)f, thus
The barometric gradient of different moments in a cardiac cycle, (f=1,2,3,4).
Wherein:
The derivative of moment Pressure versus Time, c when being different in a cardiac cycle1f(f=1,2,3,4) is four main drives
The spread speed of blood pulse wave in arteries and veins, ωmFor the angular frequency of pulse wave in artery of extremity, bmr, bmIIt is four
Coefficient m=1,2 of limb intra-arterial pulse wave, 3 ..., recursion obtains as the following formula:
Wherein
PmFor mean arterial pressure, P is function of the blood pressure in artery internal coordinate Z and time t, Fourier leaf form are as follows:Amplitude and frequency etc. are obtained by measured waveform.
The blood flow Q of different moments arteria brachialis and ankle artery four limbs in a cardiac cycle is calculated as follows in step (6)
(t)f, (f=1,2,3,4).
Wherein:
λ(t)fThe first power coefficient changed over time for blood flow；ε(t)fThe secondary power system changed over time for blood flow
Number；A(t)fIt is main to describe blood vessel bullet in a cardiac cycle for the function that vessel radius in a cardiac cycle changes over time
The dynamic process of property shrinkage expansion, Q (t)fIt is corresponding for different moments arteria brachialis in a cardiac cycle and each sampled point of ankle artery
Blood flow magnitude (f=1,2,3,4).
λ(t)fIt is indicated with following formula, (f=1,2,3,4):
Wherein:
α is human body correction factor under physiological condition, is to calculate λ (t)fWhen be related in a cardiac cycle blood flow at any time
Non-linear pulse wave propagation coefficient when variation, α=0.57；
β is human body correction factor under physiological condition, is to calculate λ (t)fWhen one cardiac cycle in vessel radius become at any time
Non-linear pulse wave propagation coefficient when change, β=0.1.173 related with vessel radius elastic dilatation；
γ is blood motion viscosity, is the ratio of hemodynamics viscosity and density of blood；
β1For the ratio of the length in body length and blood vessel stress balance of blood vessel under normal physiological conditions, β1=0.57；
β2mWhen for blood vessel stress balance radius and it is undeformed when vessel radius ratio, β2m=1.105；
α2For human body correction factor under physiological condition,
ε(t)fIt is indicated with following formula, (f=1,2,3,4):
Wherein:
For the semi-cone angle in blood vessel under natural conditions,
β0It is to calculate ε (t) for human body correction factor under physiological conditionfWhen one cardiac cycle in vessel radius become at any time
Non-linear pulse wave propagation coefficient when change, β related with vessel radius elastic shrinkage0=-1.5；
A(t)fIt is indicated with following formula, (f=1,2,3,4):
Wherein:
Wherein:
G is measured's age, and h is measured's height, and it is measured value that l, which is measured's shoulder breadth,.
Step (7), as the following formula and using estimate-modified method calculates arteria brachialis and ankle artery four limbs blood flow Q (t)f
Discreet valueWith correction value Q (t)(n+1)f, (f=1,2,3,4) is different in a cardiac cycle to which analysis meter calculates
The corresponding blood flow numerical solution of moment four limbs, (f=1,2,3,4).
N expression time series number, n=1,2 ... ... N,Indicate the four limbs blood flow discreet value at (n+1) moment,
Q(t)(n+1)fIndicating the limb blood flow correction value at (n+1) moment, G (Q (t), τ) is intermediate function, And τn
=(n-1) Δ τ, τ are time dimension, and N is the sampling number in cardiac cycle T, are numerically equal to cardiac cycle T
Interior sampling instant number.
Step (8), is calculated as follows the human bloodstream amount Q (t) of a different moments cardiac cycle,
Q (t)=δ1Q(t)1+δ2Q(t)2+δ3Q(t)3+δ4Q(t)4+ δ is wherein:
δ1, δ2, δ3, δ4It is Q (t)1, Q (t)2, Q (t)3, Q (t)4Corrected parameter, be constant term, δ is constant term, and join
Several numerical solutions is determined by the non-linear regression method in Medical Statistics.
Step (9), iterative calculation judges the error of blood flow as the following formula:
Wherein:
Q(t)0For the reference value of the human bloodstream amount of healthy population.The relative error of iterative calculation blood flow should meet mistake
Rate is less than 5%, if step (2)-(8) should be repeated by not meeting, until error determines less than 5%；If after executing aforesaid operations accidentally
Difference is still greater than 5%, then corrects the regression coefficient of blood flow fitting formula, until
The human body mean blood flow Q in a cardiac cycle is calculated as follows in step (10)m,
Step (11) calculates cardiac function blood flow parameter, cardiac output SV, cardiac output CO, heartbeat index according to the following formula
SI, cardiac index CI and total blood volume BV.
BV=2.65*BSA (L) is wherein:
HR is heart rate, is obtained by actual measurement.BSA is to indicate human body total surface area, is calculated as follows.
BSA(m2)=0.0061H (cm)+0.0128W (kg) -0.1592 is wherein:
H and W is respectively height and weight, is obtained by actual measurement.
The blood vessel parameter of peripheral blood vessel, two parameters of peripheral resistance TPR and compliance C are calculated as follows in step (12).
Wherein:
PSFor right upper extremity arteria brachialis systolic pressure, PdFor right upper extremity arteria brachialis diastolic pressure, has and obtained by actual measurement.PmFor
Mean arterial pressure is calculated as follows
The blood flow parameter of microcirculation, half turnover rate ALK of blood flow, blood flow half renewal time ALT is calculated as follows in step (13)
With the parameters such as blood mean residence time TM.
Heart acting parameter is calculated as follows in step (14), and heart mean power MPP, heart, which are often fought, moves function SKE.
MPP=Pm*Qm(mJ)
Wherein:
Ts is intersystole, and S is aortic valve orifice area 2.5-3.5cm2, ρ is density of blood, ρ=1.05*10-3kg/
cm3, PSFor right upper extremity arteria brachialis systolic pressure.
This project group has been based on the shape information of non-linear pulse wave theory and radial artery, has obtained the nothing of blood viscosity
Detection algorithm is created, resulting result is calculated and meets clinical healthy biological value range.Present invention innovatory algorithm on this basis, and
The health status of cardiovascular system is evaluated and tested based on the blood pressure pulse data of arteria brachialis and ankle artery four limbs.In addition, of the invention
Based on hardware device, easy to operate, easy to detect, the used time is few, spend it is few, it is strong for cardiac function and artery sclerosis early stage
Kang Pu cubing has very strong practicability.
Detailed description of the invention
Fig. 1 is instrumentation plan of the invention；
Fig. 2 is flow diagram of the invention；
Fig. 3 is that waveform of the invention, blood flow and cardiac functional parameter calculate schematic diagram.
Specific embodiment
Below with reference to following embodiment, the present invention will be further described.
Referring to flow diagram of the invention, specific implementation steps are as follows:
Step 1: input, acquisition, the Wave data for saving personal information and blood pressure of four limbs.
Step 2: waveform calibration, amplitude and period normalize, unify the number of collection point in waveform, obtain one completely
Average waveform, calculate the corresponding arterial pressure average value of four limbs, and extract the characteristic point in waveform.
Step 3: relevant waveform parameter is calculated, shape factor kf, pulse wave spread speed c (t)f, pulse waveform
Parameter bf, vessel radius parameter A (t)f, the first power coefficient lambda (t) that changes over time of blood flowf, blood flow changes over time
Secondary power coefficient ε (t)f；Set physiological parameter relevant to waveform and correction factor in human body, density of blood: ρ=1.05*10- 3kg/cm3, blood motion viscosity γ, the semi-cone angle in blood vessel under natural conditions:Human body correction factor: a=0.57,
β=0.1.173,β1=0.57, β2m=1.105, β0=-1.5.
Step 4: calculating the space pressure gradient F (t) of arteria brachialis and ankle arterial pressure pulse wavef。
Step 5: calculating the blood flow Q (t) of different moments arteria brachialis and ankle artery four limbs in a cardiac cyclef, and deposit
Storage.
Wherein:
Step 6: calculating the human bloodstream amount Q (t) of a different moments cardiac cycle
Q (t)=δ1Q(t)1+δ2Q(t)2+δ3Q(t)3+δ4Q(t)4+ δ is wherein:
δ1, δ2, δ3, δ4It is Q (t)1, Q (t)2, Q (t)3, Q (t)4Corrected parameter, be constant term, δ is constant term, and join
Several numerical solutions is determined by the non-linear regression method in Medical Statistics.
Step 7: iterative calculation judges the error of blood flow:
If error is greater than 5%, return step two, and re-executes step 2 to step 6, until error judgement is less than
5%；If error is still greater than 5% after executing aforesaid operations, the regression coefficient of blood flow fitting formula is corrected, until error
Determine to be less than or equal to 5%.
Step 8: calculating the human body mean blood flow Q in a cardiac cyclem,
Step 9: calculating cardiac function blood flow parameter, cardiac output SV, cardiac output CO, heartbeat index SI, cardiac index
CI and total blood volume BV.
Step 10: calculating the blood vessel parameter of peripheral blood vessel, two parameters of peripheral resistance TPR and compliance C.
Step 11: calculating the blood flow parameter of microcirculation, half turnover rate ALK of blood flow, blood flow half renewal time ALT and blood
The parameters such as mean residence time TM.
Step 12: calculating heart acting parameter, heart mean power MPP, heart, which are often fought, moves function SKE.
Claims (4)
1. a kind of device of the calculating cardiovascular functional parameter based on blood pressure of four limbs pulse wave, which is characterized in that synchronous acquisition
Totally 8 road signal divides a series of pulse waves progress wave-average filtering of acquisition to the blood pressure and pulse wave information of arteria brachialis and ankle artery
Not Huo Qu arteria brachialis and ankle artery average waveform；The wave character parameter obtained based on waveform analysis and non-linear pulse wave reason
By foundation is based on the relationship between blood pressure of four limbs pulse wave and blood flow, then according to the intrinsic spy of human heart vascular system
Parameter or geometric parameter are levied, derives the relational expression of blood flow, blood pressure of four limbs pulse data and cardiovascular functional parameter simultaneously
It is calculated；
The device successively obtains cardiovascular functional parameter according to the following steps in upper computer and lower computer:
Step (1), upper computer and lower computer initialization:
Host computer input: the personal information of measured, including gender, age, height, weight, shoulder breadth, medical history information；
The next machine testing: the pulse wave and blood pressure of measured's arteria brachialis and ankle artery four limbs；
Host computer record: the pulse waveform and blood pressure pressure value of measured's arteria brachialis and ankle artery four limbs；
Step (2), host computer detect automatically, record, saving the pulse waveform figure of four limbs, utilize the period to obtained waveform diagram
Wave-average filtering is carried out with amplitude normalized, the four limbs for respectively obtaining measured in a cardiac cycle are then calculated as follows
Mean arterial pressure Pmf, mean arterial pressure is the average value of the corresponding arterial pressure of four limbs in a cardiac cycle；
T is a cardiac cycle, and f is the label of four limbs, and f=1,2,3,4 respectively indicate right upper extremity, left upper extremity, right lower extremity, lower-left
Limb, P (t)fFor the corresponding pulse pressure value of different moments each sampled point in a cardiac cycle, f=1,2,3,4；
Step (3), is calculated as follows shape factor kf, f=1,2,3,4；
PdfFor diastolic pressure, PsfIt is the actual measured value of four limbs arteria brachialis and ankle artery for systolic pressure, f=1,2,3,4；
The spread speed c (t) of arteria brachialis and ankle arterial pressure pulse wave is calculated as follows in step (4)f:
Wherein:
P(t)fFor the corresponding pulse pressure value of different moments each sampled point in a cardiac cycle, f=1,2,3,4；
ρ is density of blood, ρ=1.05*10-3kg/cm3；
bfIt for pulse waveform parameter, is calculated as follows, f=1,2,3,4:
A is parameter related with velocity profile, a=0.57；
The space pressure gradient F (t) of arteria brachialis and ankle arterial pressure pulse wave is calculated as follows in step (5)f, to obtain one
The barometric gradient of different moments in cardiac cycle；
Wherein:
For the derivative of different moments Pressure versus Time in a cardiac cycle, c1f(f=1,2,3,4) is in artery of extremity
The spread speed of blood pulse wave, ωmFor the angular frequency of pulse wave in artery of extremity, bmr, bmIFor artery of extremity
The coefficient of interior pulse wave, m=1,2,3 ..., recursion obtains as the following formula:
Wherein:
PmFor mean arterial pressure, P is function of the blood pressure in artery internal coordinate Z and time t, Fourier leaf form are as follows:Amplitude and frequency are obtained by measured waveform；
The blood flow Q (t) of different moments arteria brachialis and ankle artery four limbs in a cardiac cycle is calculated as follows in step (6)f；
Wherein:
λ(t)fThe first power coefficient changed over time for blood flow；ε(t)fThe secondary power coefficient changed over time for blood flow；A
(t)fIt is main to describe blood vessel elasticity receipts in a cardiac cycle for the function that vessel radius in a cardiac cycle changes over time
The dynamic process of reducing and expansion, Q (t)fFor different moments arteria brachialis in a cardiac cycle and the corresponding blood of each sampled point of ankle artery
Flow value, f=1,2,3,4；
λ(t)fIt is indicated with following formula, f=1,2,3,4:
Wherein:
α is human body correction factor under physiological condition, is to calculate λ (t)fWhen be related to blood flow in a cardiac cycle and change over time
When non-linear pulse wave propagation coefficient, α=0.57；
β is human body correction factor under physiological condition, is to calculate λ (t)fWhen one cardiac cycle in vessel radius when changing over time
Non-linear pulse wave propagation coefficient, β=0.1.173 related with vessel radius elastic dilatation；
γ is blood motion viscosity, is the ratio of hemodynamics viscosity and density of blood；
β1For the ratio of the length in body length and blood vessel stress balance of blood vessel under normal physiological conditions, β1=0.57；
β2mWhen for blood vessel stress balance radius and it is undeformed when vessel radius ratio, β2m=1.105；
α2For human body correction factor under physiological condition,
R2(t)fBy this stepIn obtain
ε(t)fIt is indicated with following formula, (f=1,2,3,4):
Wherein:
For the semi-cone angle in blood vessel under natural conditions,
β0It is to calculate ε (t) for human body correction factor under physiological conditionfWhen one cardiac cycle in vessel radius when changing over time
Non-linear pulse wave propagation coefficient, β related with vessel radius elastic shrinkage0=-1.5；
A(t)fIt is indicated with following formula, f=1,2,3,4:
Wherein:
Wherein:
G is measured's age, and h is measured's height, and it is measured value that l, which is measured's shoulder breadth,；
Step (7), as the following formula and using estimate-modified method calculates arteria brachialis and ankle artery four limbs blood flow Q (t)fEstimate
ValueWith correction value Q (t)(n+1)f, f=1,2,3,4 calculate different moments four limbs in a cardiac cycle to analysis meter
Corresponding blood flow numerical solution；
N expression time series number, n=1,2 ... ... N,Indicate the four limbs blood flow discreet value at (n+1) moment, Q
(t)(n+1)fIndicating the limb blood flow correction value at (n+1) moment, G (Q (t), τ) is intermediate function, And τn
=(n-1) Δ τ, τ are time dimension, and N is the sampling number in cardiac cycle T, are numerically equal to cardiac cycle T
Interior sampling instant number；
Step (8), is calculated as follows the human bloodstream amount Q (t) of different moments in a cardiac cycle
Q (t)=δ1Q(t)1+δ2Q(t)2+δ3Q(t)3+δ4Q(t)4+ δ is wherein:
δ1, δ2, δ3, δ4It is Q (t)1, Q (t)2, Q (t)3, Q (t)4Corrected parameter, be constant term, δ is constant term, and parameter
Numerical solution is determined by the non-linear regression method in Medical Statistics；
The error judgement of step (9), blood flow iterates to calculate error in judgement as the following formula:
Wherein:
Q(t)0For the reference value of the blood flow of healthy population；The relative error of iterative calculation blood flow should meet error rate and be less than
5%, if not meeting should repeat the above steps, until error determines less than 5%；If error is still greater than after executing aforesaid operations
5%, then the regression coefficient of blood flow fitting formula is corrected, until
The human body mean blood flow Q in a cardiac cycle is calculated as follows in step (10)m:
Step (11) calculates cardiac function blood flow parameter: cardiac output SV, cardiac output CO, heartbeat index SI, the heart according to the following formula
Dirty index CI and total blood volume BV；
BV=2.65*BSA (L) is wherein:
HR is heart rate, is obtained by actual measurement；BSA is to indicate human body total surface area, is calculated as follows:
BSA(m2)=0.0061H (cm)+0.0128W (kg) -0.1592 is wherein:
H and W is respectively height and weight, is obtained by actual measurement；
The blood vessel parameter of peripheral blood vessel: two parameters of peripheral resistance TPR and compliance C is calculated as follows in step (12):
Wherein:
PSFor right upper extremity arteria brachialis systolic pressure, PdFor right upper extremity arteria brachialis diastolic pressure, has and obtained by actual measurement；It is averagely dynamic
Pulse pressure is calculated as follows
The blood flow parameter of microcirculation: half turnover rate ALK of blood flow, blood flow half renewal time ALT and blood is calculated as follows in step (13)
Liquid mean residence time TM parameter:
Step (14), heart acting parameter is calculated as follows: heart mean power MPP, heart, which are often fought, moves function SKE；
MPP=Pm*Qm(mJ)
Wherein:
Ts is intersystole, and S is aortic valve orifice area 2.5-3.5cm2, ρ is density of blood, ρ=1.05* 10-3kg/cm3, PS
For right upper extremity arteria brachialis systolic pressure.
2. the device of the calculating cardiovascular functional parameter according to claim 1 based on blood pressure of four limbs pulse wave, special
Sign is, the pulse wave signal, from 4 pressure sensors；Blood pressure signal, from 4 cufves.
3. the device of the calculating cardiovascular functional parameter according to claim 1 based on blood pressure of four limbs pulse wave, special
Sign is that the four limbs arteria brachialis and ankle arterial pulse wave need to carry out wave-average filtering, and average is that four limbs pressure sensor mentions
The step of pressure pulse wave component taken, wave-average filtering, comprising: reading Wave data first demarcates complete waveform and finds out week
Then phase carries out period and amplitude normalized to waveform, pulse wave unicast data is made to have unified amplitude and length, and
Acquire average waveform.
4. the device of the calculating cardiovascular functional parameter according to claim 1 based on blood pressure of four limbs pulse wave, special
Sign is that the non-linear pulse wave theory establishes blood pressure of four limbs, blood flow two basic haemodynamics parameters and four limbs
Relational expression between pulse wave, and according to the inherent feature parameter or geometric parameter of human heart vascular system, derive four
Limb blood pressure pulse data, the relational expression of blood flow and cardiovascular functional parameter are simultaneously calculated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610064419.2A CN105726000B (en) | 2016-01-29 | 2016-01-29 | A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610064419.2A CN105726000B (en) | 2016-01-29 | 2016-01-29 | A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105726000A CN105726000A (en) | 2016-07-06 |
| CN105726000B true CN105726000B (en) | 2019-03-22 |
Family
ID=56247107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610064419.2A Active CN105726000B (en) | 2016-01-29 | 2016-01-29 | A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105726000B (en) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018010107A1 (en) * | 2016-07-13 | 2018-01-18 | 悦享趋势科技（北京）有限责任公司 | Method and device for determining whether pulse wave is valid |
| CN106250680B (en) * | 2016-07-22 | 2019-04-05 | 曾金生 | Health of heart index detection system |
| EP3512431A1 (en) * | 2016-09-16 | 2019-07-24 | Institut National de la Sante et de la Recherche Medicale (INSERM) | Method for imaging a sample with blood and associated devices |
| CN108175387A (en) * | 2017-12-19 | 2018-06-19 | 桂林电子科技大学 | A kind of peripheral vascular resistance detection device and detection method based on electrocardio and pulse wave Morphologic Parameters |
| CN108186000B (en) * | 2018-02-07 | 2024-04-02 | 河北工业大学 | Real-time blood pressure monitoring system and method based on ballistocardiogram signal and photoelectric signal |
| CN110833403A (en) * | 2018-08-16 | 2020-02-25 | 佳纶生技股份有限公司 | Blood pressure measuring system, method and device |
| CN109394187B (en) * | 2018-11-08 | 2021-10-22 | 重庆医科大学 | Wearable cardiovascular health monitoring system based on monomer side signal detection |
| CN111345791B (en) * | 2018-12-24 | 2023-10-24 | 深圳迈瑞生物医疗电子股份有限公司 | Pulse wave measuring device |
| CN110710960B (en) * | 2019-09-12 | 2022-12-09 | 东莞市康助医疗科技有限公司 | Arterial blood vessel age estimation method and device |
| CN111493855B (en) * | 2020-04-21 | 2023-01-06 | 重庆理工大学 | System and method for non-invasive measurement of individualized cardiac output |
| CN112494022A (en) * | 2020-11-26 | 2021-03-16 | 苏州润迈德医疗科技有限公司 | Method for obtaining coronary artery blood vessel evaluation parameter and storage medium |
| CN113160921A (en) * | 2021-05-26 | 2021-07-23 | 南京大学 | Construction method and application of digital human cardiovascular system based on hemodynamics |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1631319A (en) * | 2003-12-23 | 2005-06-29 | 罗志昌 | Blood flow parameter non-invasive detecting method based on finger tip volume pulsation blood flow wave |
| CN103054562A (en) * | 2012-12-25 | 2013-04-24 | 中国科学院合肥物质科学研究院 | Cardiovascular function detection method based on multi-channel pulse wave form analysis and device thereof |
| CN104188642A (en) * | 2014-09-24 | 2014-12-10 | 中国科学院合肥物质科学研究院 | Device for noninvasive acquiring centre artery pressure waveform and relevant parameters and acquiring method thereof |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3184349B2 (en) * | 1992-12-25 | 2001-07-09 | フクダ電子株式会社 | Peripheral arteriosclerosis index measuring device |
-
2016
- 2016-01-29 CN CN201610064419.2A patent/CN105726000B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1631319A (en) * | 2003-12-23 | 2005-06-29 | 罗志昌 | Blood flow parameter non-invasive detecting method based on finger tip volume pulsation blood flow wave |
| CN103054562A (en) * | 2012-12-25 | 2013-04-24 | 中国科学院合肥物质科学研究院 | Cardiovascular function detection method based on multi-channel pulse wave form analysis and device thereof |
| CN104188642A (en) * | 2014-09-24 | 2014-12-10 | 中国科学院合肥物质科学研究院 | Device for noninvasive acquiring centre artery pressure waveform and relevant parameters and acquiring method thereof |
Non-Patent Citations (1)
| Title |
|---|
| 基于脉搏波理论的心血管功能参数无创检测算法及实验验证;吴萍;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20150315(第3期);论文第9-15页、第27-29页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105726000A (en) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105726000B (en) | A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse | |
| EP2393423B1 (en) | Calculating cardiovascular parameters | |
| CN109730663A (en) | Assessment of blood pressure method based on pulse wave conduction speed nonlinear analysis | |
| EP3571984B1 (en) | Systems for creating a model for use in detecting a peripheral arterial pressure decoupling in a subject and for detecting peripheral arterial pressure decoupling using arterial pressure waveform data and said model | |
| Willemet et al. | Computational assessment of hemodynamics-based diagnostic tools using a database of virtual subjects: Application to three case studies | |
| CA3122115C (en) | Method of predicting fluid responsiveness in patients | |
| CN102894982B (en) | Non-invasive detecting method for blood viscosity based on pulse wave | |
| JP2018534023A (en) | A method and apparatus for estimating the propagation time of an aortic pulse (aortic PTT) from a time difference measured between a plurality of reference points of a cardiac elastic diagram (BCG wave). | |
| Hacıhamdioğlu et al. | Preperitoneal fat tissue may be associated with arterial stiffness in obese adolescents | |
| Johnson et al. | Detecting aortic valve-induced abnormal flow with seismocardiography and cardiac MRI | |
| EP2416701B1 (en) | Monitoring peripheral decoupling | |
| CN104027097A (en) | Vascular function noninvasive detecting method and device | |
| Raphael et al. | Feasibility of cardiovascular magnetic resonance derived coronary wave intensity analysis | |
| CN203341726U (en) | Novel intelligent sphygmomanometer | |
| Luo et al. | The construction of unsmooth pulse images in traditional chinese medicine based on wave intensity technology | |
| CN113712524A (en) | Data processing device, system and kit for auxiliary evaluation of cardiovascular risk | |
| Wang et al. | A comparison of three types of pulse signals: Physical meaning and diagnosis performance | |
| Zaid et al. | Cardiovascular sex-differences: insights via physiology-based modeling and potential for noninvasive sensing via ballistocardiography | |
| Bollache et al. | Numerical modeling of arterial pulse wave propagation to characterize aortic hemodynamic: Validation using magnetic resonance data | |
| Hassan et al. | A portable continuous blood pressure monitoring kit | |
| Feng et al. | Continuous Estimation of Left Ventricular Hemodynamic Parameters Based on Heart Sound and PPG Signals Using Deep Neural Network | |
| WO2024031845A2 (en) | Multi-point calibration-based cuff-free dynamic blood pressure measurement method | |
| Havlik et al. | Design and realization of hardware for measurement of hemodynamic parameters | |
| Harana | Non-invasive, MRI-based calculation of the aortic blood pressure waveform by | |
| CN103340614A (en) | Novel intelligent sphygmomanometer |
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 |