The specific embodiment
For making the purpose, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.
Need to prove that in accompanying drawing or description description, similar or identical part is all used identical figure number.The implementation that does not illustrate in the accompanying drawing or describe is form known to a person of ordinary skill in the art in the affiliated technical field.In addition, although this paper can provide the demonstration of the parameter that comprises particular value, should be appreciated that parameter need not definitely to equal corresponding value, but can in acceptable error margin or design constraint, be similar to corresponding value.
In one exemplary embodiment of the present invention, provide a kind of method of measuring pulse frequency.As shown in Figure 1, this measurement pulse frequency method comprises:
Steps A is by the pulse transducer time domain pulse signal of collection 2s time;
This pulse transducer is installed with in the position of tested object pulse, such as chest, wrist, finger etc., the pulse of the tested object of Real-time Collection.Preferably, this pulse transducer is attached at tested object chest, carries out the pulse frequency monitoring with convenient tested object under kinestate.
In this step, the sampling duration of time domain pulse signal is 1-2s, and 2s is more quite a lot of, and sample frequency is fixed because for frequency domain algorithm, and sampling number is more, and frequency resolution is higher.Comprehensive requirement of real-time again, 2s is more quite a lot of.Certainly, in the not high situation of requirement of real-time, also can the proper extension data acquisition time, such as: 4s, 8s, 10s, 20s, 60s etc.
Step B, the acceleration signal of synchronous acquisition 2s is measured tested object at the acceleration of X, Y and three directions of Z: a by 3-axis acceleration sensor
x, a
yAnd a
z,, counting of collection is n; Then add up the value of three axial acceleration in [3,3], [6,6], [10,10], [15,15], counting between [20,20] is designated as respectively a
X1, a
X2, a
X3, a
X4, a
X5a
Y1, a
Y2, a
Y3, a
Y4, a
Y5a
Z1, a
Z2, a
Z3, a
Z4, a
Z5A wherein
X1, a
X2, a
X3, a
X4, a
X5For tested object at the acceleration of directions X respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5]; a
Y1, a
Y2, a
Y3, a
Y4, a
Y5For tested object at the acceleration of Y-direction respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5]; a
Z1, a
Z2, a
Z3, a
Z4, a
Z5For tested object at the acceleration of Z direction respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5].
In the present embodiment, the 3-axis acceleration sensor of employing can obtain the acceleration of tested object on three directions.It will be apparent to those skilled in the art that also and can adopt the sensor that can obtain the one direction acceleration, and needn't strictly distinguish three directions, thereby reduce cost.
In this enforcement, pulse transducer and 3-axis acceleration sensor are integrated on the equipment.
Step C judges that tested object is at the acceleration of X, Y and three directions of the Z statistical value a at the 5th class
X5, a
Y5, a
Z5Whether all greater than predetermined threshold value T
Min, if so, execution in step D, otherwise, re-execute steps A;
Generally, this predetermined threshold value T
MinGet 0.8n, n is for once gathering the number of acceleration information on each axle.In this step, a
X5, a
Y5, a
Z5Whether all greater than predetermined threshold value T, limit exactly acceleration and surpass [20,20] interval ratio data less than 0.2n, also just be equivalent to accekeration and be greater than 0.8n in [20,20] interval ratio.
Experiment showed, to surpass [20,20] interval ratio by acceleration more large-signal is poorer that signal just can not have been differentiated when surpassing 0.2n, the follow-up pulse value that draws will be very inaccurate.Therefore, the acceleration restriction of tested object is very necessary within the specific limits, otherwise will need resampling.
Step D, employing free transmission range are that the band filter of 0.65Hz-3.20Hz carries out filtering to the time domain pulse signal;
In the prior art, usually adopt the method for wavelet filtering that the time domain pulse signal is carried out filtering, but Wavelets Filtering Algorithm is comparatively complicated, holding time is long.And method processing procedure of the present invention is simple, and the time is short, is conducive to the Real-time Collection pulse data.
Step e to the time domain pulse signal, is identified each pulse wave by crest and the trough of pulse wave master ripple, dicrotic notch and dicrotic wave, obtains the first pulse frequency n according to the time domain pulse signal that gathers
1
As shown in Figure 2, adopt and identify pulse wave by the crest (C) of pulse wave master crest (A), dicrotic notch crest (B), dicrotic wave and the trough (D) of dicrotic wave, its advantage is to realize accurate identification.It will be apparent to those skilled in the art that and also can take the mode of Direct Recognition pulse wave master ripple to identify pulse wave, doing like this accuracy rate can be slightly poor, but can't affect enforcement of the present invention.
Step F is carried out Fourier (FFT) conversion with the time domain pulse signal, obtains the frequency domain pulse signal;
Step G for the frequency domain pulse signal, determines its crest respective value in the 0.65Hz-3.20Hz frequency range, and namely the pulse frequency value obtains the second pulse frequency n by this pulse frequency
2
Experiment showed, that the maximum characteristics of utilizing frequency domain method to carry out pulses measure are exactly good stability, capacity of resisting disturbance is strong, and still, if simple proportion method can make the real-time variation again, so present embodiment adopts the pulse algorithm of time-frequency mixing.
Step H utilizes following formula, by the first pulse frequency n
1With the second pulse frequency n
2Calculate the final pulse frequency of the unit interval of tested object:
N=n
1×w
1+n
2×w
2
Wherein, w
1And w
2Be respectively the weight of the first pulse frequency and the second pulse frequency, both are relevant at each interval statistical value at the acceleration of X, Y and three directions of Z with tested object, and both sums are 1.
In the present embodiment, the present invention's degree of will speed up is divided into five class statistics, is respectively [3,3], and [6,6], [10,10], [15,15], [20,20] are again according to three axial acceleration statistical value a on each class
X1, a
X2, a
X3, a
X4, a
X5a
Y1, a
Y2, a
Y3, a
Y4, a
Y5a
Z1, a
Z2, a
Z3, a
Z4, a
Z5, and the sampling n that always counts determines the weights of corresponding the first pulse frequency and the second pulse frequency:
As min (a
X1, a
Y1, a
Z1During) 〉=0.8n, w
1=0.9, w
2=0.1;
Be false but min (a when above
X2, a
Y2, a
Z2During) 〉=0.8n, w
1=0.7, w
2=0.3;
Be false but min (a when above
X3, a
Y3, a
Z3During) 〉=0.8n, w
1=0.5, w
2=0.5;
Be false but min (a when above
X4, a
Y4, a
Z4During) 〉=0.8n, w
1=0.3, w
2=0.7;
Be false but min (a when above
X5, a
Y5, a
Z5During) 〉=0.8n, w
1=0.1, w
2=0.9.
So far, present embodiment is measured the end of pulse frequency process.
The invention provides a kind ofly applicable in people's daily life, especially under motion conditions, measure the method for pulse frequency.The data that the method utilizes the pulse transducer collection to come up in conjunction with the acceleration information of a 3-axis acceleration sensor, combine the calculating Pulse Rate with time domain pulse information and frequency domain pulse information again.The present invention has substantial degradation to the pulse prescription, and stable performance is calculated simply, and the result of calculation of pulse is also more accurate.
Need to prove that above-mentioned definition to each element is not limited in various concrete structures or the shape of mentioning in the embodiment, those of ordinary skill in the art can replace simply to it with knowing, for example:
The weighted value of (1) first pulse frequency and the second pulse frequency can be adjusted as required;
(2) number of acceleration class, upper lower limit value wherein can be adjusted as required;
(3) 3-axis acceleration sensor can replace with other multiaxis acceleration transducers.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above only is specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of making, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.