Summary of the invention
The present invention is directed to above-mentioned deficiency provide a kind of can continuous measurement and the system for detecting blood flow parameters of analyzing pulse wave.
Principle of the present invention is through the pulse wave curves to being gathered, and adopts the area barycenter displacement analytic process that its wave character is measured, and is quantized into characteristic index.System can also be that calculated factor is done further calculating and generated a series of quantized and hearts, blood vessel, hemodynamic index that blood flow is relevant with above-mentioned characteristic index.Above-mentioned serial characteristic index and the further hemodynamic index that generates are carried out continuous monitoring, record, observation, storage, playback, statistics, analysis and research, promptly realized noinvasive, in real time, the pulse wave characteristic of continuously monitoring human body and the process of hemodynamic state and its dynamic change of tracing observation.
For realizing above-mentioned purpose, system of the present invention comprises:
Data acquisition module is used to gather pulse wave signal, gathers finger tip volumetric blood flow pulse wave signal by the finger-clipped photoelectric sensor, to signal amplify, processing such as filtering and analog digital conversion; System supports the blood pressure acquisition module as optional module, when needed the data analysis module of connecting system;
Data analysis module is used to analyze the pulse wave digital signal from data acquisition module, and the pulse waveform characteristic quantity is changed into characteristic index, can be that calculated factor is done further calculating and generated a series of quantized hemodynamic indexs in addition with the characteristic index.Above-mentioned characteristic index and/or hemodynamic index are carried out time and frequency domain analysis, generate parameter list and statistical graph.
Specifically, data analysis module calculates at least a following characteristic index according to pulse wave signal:
1) integral body under the pulse wave curves and local area and area coefficient thereof comprise whole pulse wave-wave area of pictural surface A ', the ripple area of pictural surface Ac ' of the ripple area of pictural surface Ad ' of upstroke counterpart and area coefficient Cd thereof and decent counterpart and area coefficient Cc thereof;
2) position of centre of gravity of area under the pulse wave curves comprises the ripple area of pictural surface center of gravity G1 of whole pulse wave-wave area of pictural surface center of gravity G and upstroke counterpart and the ripple area of pictural surface center of gravity G2 of decent counterpart;
3) the pairing pulsation period T of pulse wave-wave figure, comprise the upstroke counterpart the time phase Td and decent counterpart time Tc mutually;
4) to 1), 2), 3) described in index carry out the new characteristic index that computing produces.
Wherein 4) the new index described in includes but not limited to: Ad '/Ac ', Cd/Cc, Cc-Cd, Ad '/A ', Ac '/A ', X2-X1, Y1-Y2, X-X1, Y1-Y, X2-X, Y-Y2, (tCc-tc)/(tc-tDc), (qc-qCc)/(qDc-qc), (X2-X)/(X-X1), (Y-Y2)/(Y1-Y), X-Td/T, (X2-Td/T), (Td/T-X1), Td/Tc, (Tc-Td)/T or the like, can form new index according to these index computings in addition.
Referring to Fig. 2, area under the described pulse wave curves comprises whole pulse wave-wave area of pictural surface A ' and calculates the ripple area of pictural surface Ad ' of upstroke counterpart respectively and the ripple area of pictural surface Ac ' of decent counterpart.
Phase T when pulse wave-wave figure is pairing: pulsation period;
The time phase Td of upstroke counterpart: angioplerosis phase;
The time phase Tc of decent counterpart: blood vessel retraction phase.
Referring to Fig. 3 and Fig. 8 a, b, the barycentric coodinates position of area under the described pulse wave curves comprises whole pulse wave-wave area of pictural surface center of gravity G coordinate (tc, qc) and relative coordinate (X, Y) and the ripple area of pictural surface center of gravity G1 coordinate (t of upstroke counterpartDc, qDc) and relative coordinate (X1 is Y1) with the ripple area of pictural surface center of gravity G2 coordinate (t of decent counterpartCc, qCc) and suitable coordinate (X2, Y2);
Data analysis module is further analyzed and is calculated blood pressure and pulse wave number word signal, generates hemodynamic index.Described hemodynamic index comprises a kind of in " blood vessel elasticity, vascular compliance, Peripheral resistance, mean arterial pressure, cardiac output, cardiac index, SV, SI etc. " index at least.
Said result of calculation and parameter list comprise a kind of in the following numerical value at least: value of calculation, period average and overall average in real time; Said statistical graph comprises a kind of in the following form at least: normal distribution, Lorenz scatterplot and trendgram.
In view of finger tip volumetric blood flow pulse wave is gathered with respect to radial artery, brachial artery, carotid artery etc. more easily; Especially in real time, dynamically, the continuous acquisition aspect has advantage; Yet, for the pressure wave that radial artery, brachial artery, carotid artery etc. are located, finger tip volumetric blood flow pulse wave signal a little less than; Waveform is also relatively slicker and more sly; Its rising and descend all compare slowly, and dicrotic wave is low flat even do not have a crest under a lot of situation, and this is that the filter action to waveform generation causes because pressure pulse wave is through the repeatedly branch leakage of microvascular resistance of tip and capillary network and blood capillary surrounding tissue etc.Because discernible intuitionistic feature point is less, increased suitable difficulty just for the quantification and the analysis of waveform.Brand-new characteristic point and quantizating index have been generated through system of the present invention to pulse wave; And it is different with characteristic point on the pulse wave curves shown in Figure 1; No matter how individual variation, kinestate, physiology and pathological condition change; These characteristic points and quantizating index all can not disappear, and are fit to continue tracing observation, thereby efficiently solve this difficult problem of pulse wave continuous monitoring.
Adopting system of the present invention to carry out that pulse wave measurement and hemodynamics detect is to be based upon on the continuous monitoring basis of pulse wave.Finger tip volumetric blood flow pulse wave signal includes the unify information of microcirculation aspect of a large amount of cardiovascular systems; Article one, the morphological characteristic of pulse wave curves is blood flow and the blood vessel results of interaction by particular individual to be detected, has reflected human body microvascular hemodynamics characteristic under special time and state.Under the unobstructed situation of trunk and blood capillary blood flow, wherein include abundant human bloodstream dynamic information.The variation tendency of continuous monitoring pulse wave characteristic index just can be followed the tracks of, observes, the slight change and the rule of the hemodynamics characteristic of analysis, researching human body.
Area under the finger tip volumetric blood flow pulse wave curves Q (t) can be expressed as:
Qmax is a finger volume pulse blood flow maximum; Qmin is a finger volume pulse blood flow minima; T is the pulsation period.See Fig. 4.
For further wave character being quantized, with reflection from blood vessel begin full up to maximum process be the upstroke of pulse wave corresponding the time be defined as filling period mutually, represent with Td.With reflected blood vessel begin under the arterial wall elastic reaction to bounce back up to venous backflow make the full process of heart be the decent of pulse wave corresponding the time be defined as the retraction phase mutually, represent with Tc.
Whole oscillogram cartographic represenation of area under the pulse wave curves Q (t) is:
Pulse wave curves Q (t) during following filling period corresponding area can be expressed as:
Cd=Ad'/A ' is the ratio of pulse wave curves corresponding area and whole oscillogram area during following filling period, is defined as pulse wave filling period area coefficient.See Fig. 5.
Pulse wave curves Q (t) down retraction during the phase corresponding area can be expressed as:
Cc=Ac'/A ' is the ratio of retraction corresponding area and whole oscillogram area during the phase under the pulse wave curves, is defined as the pulse wave phase area coefficient that bounces back.See Fig. 5.
Area barycentric coodinates under the pulse wave curves Q (t) are seen Fig. 6, can be expressed as:
Abscissa:Vertical coordinate:
The morphological characteristic that quantizes pulse wave also can realize through the relative coordinate of reference area center of gravity, sees Fig. 7.If T is a long measure, Qmax-Qmin is a width unit, we can obtain waveform area center of gravity G (X, relative coordinate value Y) is:
With quadrat method can obtain respectively Q (t) down with the barycentric coodinates of filling period and corresponding area of retraction phase, same, the morphological characteristic that quantizes pulse wave also can be seen Fig. 8 through the relative coordinate realization of the corresponding area center of gravity of calculating.
Pulse wave curves Q (t) is the area barycentric coodinates G1 (t of upstroke part downDc, qDc) see and can be expressed as Fig. 8 a:
Abscissa:Vertical coordinate:
If T is a long measure, Qmax-Qmin is a width unit, we can obtain this portion waveshape area center of gravity relative coordinate G1 (X1 Y1), sees Fig. 8 b, can be expressed as:
Pulse wave curves Q (t) is the area barycentric coodinates G2 (t of decent part downCc, qCc) see and can be expressed as Fig. 8 a:
Abscissa:Vertical coordinate:
If T is a long measure, Qmax-Qmin is a width unit, we can obtain this portion waveshape area center of gravity relative coordinate G2 (X2 Y2), sees Fig. 8 b, can be expressed as:
Area A ', Ad ', Ac ' and phase T, Td, Tc when corresponding, comprise area coefficient Cd, Cc and barycentric coodinates (X, Y), (X1, Y1), (X2 Y2) is the index of the whole and local feature of reflection pulse wave.Whole and local area or area coefficient variation and corresponding displacement of center of gravity thereof through the continuous detecting pulse wave curves; Can observe the slight change of pulse wave; Above-mentioned characteristic index is carried out time and frequency domain analysis, can dynamic observe its change procedure, trend and rule.
Through the continuous detecting hemodynamic index, can observe the slight change of cardiovascular system, hemodynamic index is carried out time and frequency domain analysis, can dynamic observe its change procedure, trend and rule.
System of the present invention has generated brand-new characteristic point and quantizating index for the pulse wave signal of being gathered; And no matter how individual variation, kinestate, physiology and pathological condition change; These characteristic points and quantizating index all can not disappear; Be fit to continue tracing observation, thereby efficiently solve this difficult problem of pulse wave continuous monitoring.Above-mentioned characteristic index is that life science provides effective technical means; Through a series of quantifications to the pulse wave characteristic; Help to distinguish the difference between interindividual difference and same individual different conditions more carefully, in depth study the life scientific phenomena.
The data acquisition module of system of the present invention can comprise finger-clipped photoelectric sensor, amplifier, wave filter and mould/number converter, the pulse wave signal that pick off is caught through signal amplification, filtering, change into digital signal and send the data analysis processing module to.System supports the blood pressure acquisition module as optional module, and the data analysis module of connecting system is measured blood pressure automatically when needed; System also supports the pressure value that other sphygomanometers are measured is manually imported native system.
Data processing module can be a high-performance processor, is sampled by processor and traces, and writes down each finger tip volumetric blood flow pulse waveform figure.The oscillogram characteristic quantity that is obtained is changed into characteristic index, comprise area and area coefficient and when corresponding mutually with barycentric coodinates etc., also comprise by These parameters calculating must all characteristic indexs and/or a series of hemodynamic indexs that generated.Testing result can be selected the storage playback, output to other computers, server, display or printer according to user's needs, can also set up sample database when needing, and generates the assessment template, is connected into central monitoring system network or remote medical service system.
System of the present invention have signal characteristic clear, do not receive waveform restriction, be easy to observe and characteristics such as analysis, especially be fit in real time, dynamically, tracing observation pulse wave continuously.Can be used for fields such as physical stamina test, motion monitoring, training examination and judging, physical fitness research, disorder in screening, healthy classification and health control towards sports professional person and ordinary populace.
The specific embodiment
Come further to set forth the present invention below in conjunction with concrete embodiment.Should be appreciated that these embodiment only are used to explain the present invention, and can not limit protection scope of the present invention.
Embodiment 1
As shown in Figure 9, in this example, system for detecting blood flow parameters comprises: data acquisition module, be used to gather pulse wave signal, gather finger tip volumetric blood flow pulse wave signal by the finger-clipped photoelectric sensor, to signal amplify, processing such as filtering and analog digital conversion; Data analysis module is used to analyze the pulse wave digital signal from data acquisition module, and the pulse waveform characteristic quantity is changed into characteristic index, and system also is that calculated factor is done further calculating and generated a series of quantized hemodynamic indexs with the characteristic index.These parameters is carried out time and frequency domain analysis, generate parameter list and statistical graph.
These characteristic indexs comprise:
1) integral body under the pulse wave curves and local area and area coefficient thereof comprise whole pulse wave-wave area of pictural surface A ', the ripple area of pictural surface Ac ' of the ripple area of pictural surface Ad ' of upstroke counterpart and area coefficient Cd thereof and decent counterpart and area coefficient Cc thereof;
2) position of centre of gravity of area under the pulse wave curves comprises the ripple area of pictural surface center of gravity G1 of whole pulse wave-wave area of pictural surface center of gravity G and upstroke counterpart and the ripple area of pictural surface center of gravity G2 of decent counterpart;
3) the pairing pulsation period T of pulse wave-wave figure, comprise the upstroke counterpart the time phase Td and decent counterpart time Tc mutually;
4) to 1), 2), 3) described in index carry out the new characteristic index that computing produces.
Wherein 4) the new index described in includes but not limited to: Ad '/Ac ', Cd/Cc, Cc-Cd, Ad '/A ', Ac '/A ', X2-X1, Y1-Y2, X-X1, Y1-Y, X2-X, Y-Y2, (tCc-tc)/(tc-tDc), (qc-qCc)/(qDc-qc), (X2-X)/(X-X1), (Y-Y2)/(Y1-Y), X-Td/T, (X2-Td/T), (Td/T-X1), Td/Tc, (Tc-Td)/T or the like, can form new index according to these index computings in addition.
Data analysis module is also done further calculating as calculated factor and is generated a series of quantized hemodynamic indexs, and data analysis module is also carried out time domain and/or frequency-domain analysis to characteristic index, and generates result of calculation, parameter list and/or statistical graph.Said result of calculation and parameter list comprise: real-time value of calculation, period average and overall average; Said statistics icon comprises: normal distribution, Lorenz scatterplot and trendgram.
Data acquisition module comprises finger-clipped photoelectric sensor, amplifier, wave filter and mould/number converter, the pulse wave signal that pick off is caught through signal amplification, filtering, change into digital signal and send the data analysis processing module to.But blood pressure acquisition module and cuff are arrangement, both can select to take the built-in blood pressure module of system to measure blood pressure automatically, and system also supports the pressure value that other sphygomanometers are measured is manually imported native system.
This example uses the finger-clipped photoelectric sensor to measure pulse wave, has or not blood pressure module that system is generated the pulse waveform characteristic index and has no influence.Seeing that the finger-clipped photoelectric sensor is for the pressure pulse pick off; Require not strict to the detection position; Sensor probe is easy to location and convenient fixing; The adaptability of acquired signal, stability and repeatability are strong, can satisfy Static Detection and the dynamically requirement of continuous monitoring, especially are fit to carry out time and frequency domain analysis.
Data processing module is a high-performance processor, is sampled by processor and traces, and writes down each finger tip volumetric blood flow pulse waveform figure.The oscillogram characteristic quantity that is obtained is changed into characteristic index, comprise area and area coefficient and when corresponding mutually with barycentric coodinates etc., also comprise by These parameters calculating must all characteristic indexs, generate a series of hemodynamic indexs.Testing result can be selected the storage playback, output to other computers, server or printer according to user's needs, can also set up sample database when needing, and generates the assessment template, is connected into central monitoring system network or remote medical service system.
Embodiment 2 clinical practices
Use system of the present invention through the finger tip image data, and analyze.The common healthy population sample of six all ages and classes sections that table 1 demonstration is randomly drawed is totally 2000 examples; Age, span scope was at 18~82 years old; Gather the finger tip volumetric blood flow pulse wave signal under each individual rest state; Adopt system of the present invention that the pulse wave characteristic index that is obtained is calculated and added up, the result is following:
The pulse wave form feature analysis of the common healthy population of table 1 all ages and classes section
| Age bracket | Mean age | Sample number | (X2-X)/(X-X1) age characteristics index average | Standard deviation |
| 18~30 years old | 24.10 | 380 | 0.345 | 0.061 |
| 31~40 years old | 34.00 | 310 | 0.385 | 0.055 |
| 41~50 years old | 45.73 | 310 | 0.456 | 0.053 |
| 51~60 years old | 54.02 | 310 | 0.508 | 0.048 |
| 61~70 years old | 65.61 | 310 | 0.558 | 0.047 |
| 71~82 years old | 75.63 | 380 | 0.586 | 0.045 |
Table 1 and Figure 10 show that the finger tip volumetric blood flow pulse wave characteristic index that six groups of crowds of all ages and classes section are obtained demonstrates variation clocklike along with the increase at age under rest state; Wherein (X2-X)/(X-X1) is the trend that rises gradually with age growth; Here we use (X2-X)/(X-X1) as " age characteristics index ", can find the Changing Pattern that this index of age groups rose with the age.
Randomly draw each 400 example of sample of adhering to sports for a long time and participating in sports hardly; The range of age was at 18~82 years old; Be divided into " motion is arranged " and " non-motion " two groups, gather the finger tip volumetric blood flow pulse wave under the rest state, extract (X2-X)/(X-X1) index; Be the age characteristics index, the result sees table 2 and table 3.Further will " motion be arranged " and " non-motion " two groups compares with 2000 routine healthy population Figure 10, the result sees Figure 13.
Table 2 all ages and classes section has the pulse wave form feature analysis of exercise group
| Age bracket | Mean age | Sample number | (X2-X)/(X-X1) age characteristics index average | Standard deviation |
| 18~30 years old | 25.20 | 68 | 0.303 | 0.062 |
| 31~40 years old | 34.08 | 65 | 0.343 | 0.059 |
| 41~50 years old | 45.98 | 66 | 0.411 | 0.056 |
| 51~60 years old | 54.87 | 65 | 0.465 | 0.053 |
| 61~70 years old | 66.01 | 67 | 0.517 | 0.049 |
| 71~82 years old | 74.99 | 69 | 0.547 | 0.048 |
Table 2 and Figure 11 show the Changing Pattern that the age rises with the age at 18~82 years old " age characteristics index " that the exercise group crowd is arranged.
The pulse wave form feature analysis of the non-exercise group of table 3 all ages and classes section
| Age bracket | Mean age | Sample number | (X2-X)/(X-X1) age characteristics index average | Standard deviation |
| 18~30 years old | 24.90 | 68 | 0.416 | 0.071 |
| 31~40 years old | 34.82 | 66 | 0.451 | 0.067 |
| 41~50 years old | 46.03 | 65 | 0.513 | 0.061 |
| 51~60 years old | 55.01 | 66 | 0.560 | 0.051 |
| 61~70 years old | 64.98 | 67 | 0.598 | 0.045 |
| 71~82 years old | 75.61 | 68 | 0.625 | 0.039 |
Table 3 and Figure 12 show the Changing Pattern that the age rises with the age at 18~82 years old non-exercise group crowd " age characteristics index ".
Figure 13 has the average of exercise group, non-exercise group and healthy subjects group to compare; The result shows that the sample of adhering to sports for a long time the characteristic that becomes younger occurred with respect to the hemodynamics characteristic of common healthy sample; The sample of on the contrary, participating in sports hardly the aging characteristic then occurred with respect to the hemodynamics characteristic of common healthy sample.Show that system of the present invention uses this index to can be used for the healthy classification to different crowd.
Figure 14 is the Lorenz scatterplot of " the age characteristics index " of a certain sample of employing system of the present invention output.
Figure 15 is normal distribution (Nogata) figure of " the age characteristics index " of a certain sample of employing system of the present invention output.
Obviously; " the age characteristics index " of one of a series of pulse wave characteristic indexs that generated by system of the present invention can reflect the hemodynamics characteristic well and increase old and feeble corresponding relation in age; And adhere to sports for a long time and participate in sports hardly can obtaining wanting useful quantitative to the exponential influence degree of age characteristics, can observe old and feeble track continuously through this index.And coming more simple and conveniently to hemodynamics The Characteristic Study, statistics and analysis through the measurement of volumetric blood flow pulse wave nuance is also feasible, the result is also more clear obviously.
Above-mentioned application result shows, system of the present invention can be with a series of integral body and the local characteristic quantification of pulse wave, and can further generate more pulse wave characteristic index through the combination calculation to these eigenvalues, is applied in the real-time continuous monitoring.As be applied in stress and dynamic changing process thereof the monitoring of human body under the sports load, and therefrom find traditional detection method detect less than pulse wave change and difference.Can be used for the noinvasive fast detecting of hemodynamic state and the typing of different crowd; Through the tendency of track record finger tip volumetric blood flow pulse wave characteristic index, the information such as beginning and ending time point, amplitude of variation and cycle of appearance differentiation, the hemodynamic situation of quantitative analysis human body, the research body is fitted the relation of ability and motion risk.This system is simple to operate, easy to use, characteristics be particularly useful for moving monitoring, physical stamina test, training guidance and motion intervention etc. efficiently and effectively.
Simultaneously; As a kind of system for detecting blood flow parameters based on pulse wave measurement and analysis; Through quantizing the series of features index of finger tip volumetric blood flow pulse wave, we can find pulse waveform attitude and hemodynamics characteristic between different crowd nuance and and physiology and pathological state between relation.Through the collection and the statistical analysis of large sample amount, set up crowd's sample database, for digitized template Evaluation Method and intelligent medical logic calculation function module provide scientific basis.