Detailed description of the invention
For making the object, 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.
It should be noted that, in accompanying drawing or description description, similar or identical part is all used identical figure number.The implementation that does not illustrate in accompanying drawing or describe is form known to a person of ordinary skill in the art in affiliated technical field.In addition, although the demonstration of the parameter that comprises particular value can be provided herein, should be appreciated that, parameter is without definitely equaling 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, 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, 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, to facilitate tested object to carry out pulse frequency monitoring 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 because for frequency domain algorithm, sample frequency is fixed, 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 situation that requirement of real-time is not high, also can proper extension data acquisition time, for example: 4s, 8s, 10s, 20s, 60s etc.
Step B, the acceleration signal of synchronous acquisition 2s, measures the acceleration of tested object in X, Y and tri-directions of Z: a by 3-axis acceleration sensor
x, a
yand a
z,, the counting as n of collection; 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
x5; a
y1, a
y2, a
y3, a
y4, a
y5; a
z1, a
z2, a
z3, a
z4, a
z5.Wherein a
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 in three directions.It will be apparent to those skilled in the art that and also can adopt the sensor that can obtain one direction acceleration, and needn't strictly distinguish three directions, thereby reduce costs.
In this enforcement, pulse transducer and 3-axis acceleration sensor are integrated on an equipment.
Step C, judge tested object at the acceleration of X, Y and tri-directions of Z the statistical value a at the 5th class
x5, a
y5, a
z5whether be all greater than predetermined threshold value T
min, if so, execution step D, otherwise, re-execute steps A;
Under normal circumstances, this predetermined threshold value T
minget 0.8n, n is the number that once gathers acceleration information on each axle.In this step, a
x5, a
y5, a
z5whether be all greater than predetermined threshold value T, limit exactly acceleration and exceed [20,20] interval ratio data and be less than 0.2n, be also just equivalent to accekeration and be greater than 0.8n in [20,20] interval ratio.
Experiment showed, to exceed [20,20] interval ratio by acceleration more large-signal is poorer, in the time exceeding 0.2n, signal just can not have been differentiated, and the follow-up pulse value drawing will be very inaccurate.Therefore, by very necessary within the specific limits the acceleration restriction of tested object, otherwise will need resampling.
Step D, adopts the band filter that free transmission range is 0.65Hz-3.20Hz to carry out filtering to time domain pulse signal;
In the prior art, conventionally adopt the method for wavelet filtering to carry out filtering to time domain pulse signal, but Wavelets Filtering Algorithm is comparatively complicated, holding time is long.And method processing procedure of the present invention is simple, the time is short, is conducive to Real-time Collection pulse data.
Step e, to time domain pulse signal, identifies each pulse wave by crest and the trough of the main ripple of pulse wave, dicrotic notch and dicrotic wave, obtains the first pulse frequency n according to the time domain pulse signal gathering
1;
As shown in Figure 2, adopt and identify pulse wave by the crest (C) of the main crest of pulse wave (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 the main ripple of Direct Recognition pulse wave to identify pulse wave, doing like this accuracy rate can be slightly poor, but can't affect enforcement of the present invention.
Step F, carries out Fourier (FFT) conversion by time domain pulse signal, obtains frequency domain pulse signal;
Step G, for frequency domain pulse signal, determines its crest respective value in 0.65Hz-3.20Hz frequency range, i.e. pulse frequency value, obtains the second pulse frequency n by this pulse frequency
2;
Experiment showed, that the maximum feature of utilizing frequency domain method to carry out pulses measure is exactly good stability, capacity of resisting disturbance is strong, and still, if simple proportion method can make again real-time variation, therefore the present embodiment adopts the pulse algorithm that time-frequency mixes.
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 tri-directions of Z to 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], then according to three axial acceleration statistical value a on each class
x1, a
x2, a
x3, a
x4, a
x5; a
y1, a
y2, a
y3, a
y4, a
y5; a
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
z1when)>=0.8n, w
1=0.9, w
2=0.1;
When being false above but min (a
x2, a
y2, a
z2when)>=0.8n, w
1=0.7, w
2=0.3;
When being false above but min (a
x3, a
y3, a
z3when)>=0.8n, w
1=0.5, w
2=0.5;
When being false above but min (a
x4, a
y4, a
z4when)>=0.8n, w
1=0.3, w
2=0.7;
When being false above but min (a
x5, a
y5, a
z5when)>=0.8n, w
1=0.1, w
2=0.9.
So far, the present embodiment measurement pulse frequency process finishes.
The invention provides one applicable in people's daily life, especially, under motion conditions, measure the method for pulse frequency.The data that the method utilizes pulse transducer collection to come up, then in conjunction with the acceleration information of a 3-axis acceleration sensor, time domain pulse information and frequency domain pulse information are combined to calculating Pulse Rate.The present invention has substantial degradation to pulse prescription, and stable performance is calculated simply, and the result of calculation of pulse is also more accurate.
It should be noted that, the above-mentioned definition to each element is not limited in various concrete structures or the shape in embodiment, mentioned, and those of ordinary skill in the art can know simply and replace it, 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; object of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any amendment of making, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.