CN102727185B - Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method - Google Patents

Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method Download PDF

Info

Publication number
CN102727185B
CN102727185B CN 201210249466 CN201210249466A CN102727185B CN 102727185 B CN102727185 B CN 102727185B CN 201210249466 CN201210249466 CN 201210249466 CN 201210249466 A CN201210249466 A CN 201210249466A CN 102727185 B CN102727185 B CN 102727185B
Authority
CN
China
Prior art keywords
span
energy consumption
heart rate
acceleration
axis
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
CN 201210249466
Other languages
Chinese (zh)
Other versions
CN102727185A (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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN 201210249466 priority Critical patent/CN102727185B/en
Publication of CN102727185A publication Critical patent/CN102727185A/en
Application granted granted Critical
Publication of CN102727185B publication Critical patent/CN102727185B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The embodiment of the invention provides a motion energy consumption measuring apparatus. The motion energy consumption measuring apparatus comprises a human body individuation parameter acquisition unit, a heart rate sensor, a basic energy consumption calculating unit, an acceleration sensor, a motion acting calculating unit and a motion energy consumption calculating unit, wherein the human body individuation parameter acquisition unit is used for acquiring and storing human body individuation parameters; the heart rate sensor is used for measuring heart rate HR; the basic energy consumption calculating unit is used for calculating basic energy consumption by utilizing the measured heart rate HR and the human body individuation parameters; the acceleration sensor is used for measuring the motion acceleration; the motion acting calculating unit is used for calculating the external acting by utilizing the motion acceleration and the weight W; and the motion energy consumption calculating unit is used for calculating the motion energy consumption by adding the basic energy consumption BEE and the external acting EEact. The invention also provides a measuring method of the motion energy consumption. According to the measuring apparatus and the measuring method, linear regression equations of the basic energy consumption, the heart rate and the individuation parameters are built, the motion acceleration is measured by using the acceleration sensor, the influence factors of different directions are further considered, the external acting calculation is applied by the weight in the three directions, and the heart rate and the motion acceleration are effectively combined, so that the measurement accuracy is improved.

Description

A kind of sports energy consumption measuring instrument and measuring method based on heart rate and acceleration
Technical field
The present invention relates to motion and can consume fields of measurement, particularly a kind of sports energy consumption measuring instrument and measuring method based on heart rate and acceleration.
Background technology
Along with the development of science and technology, people's living standard is more and more higher at present, and human motion is increasing to the importance of health.The very few purpose that can not arrive sport and body-building or fat-reducing of quantity of motion during human motion, quantity of motion too much can allow again the people feel exhausted and damage the health of human body.Therefore must measure accurately energy expenditure in the human motion process in order to reach scientific and reasonable motion.
The energy that human body consumed in one day roughly is divided into three parts, and the human body Basal Energy Expenditure is basal energy expenditure, and various energy expenditure of physical action are the food heat effect with satisfying specific dynamic action.The human body basal energy expenditure accounts for 65% of one day total power consumption of human body, and the activity energy expenditure accounts for 25% of total power consumption, and the food heat effect accounts for 10% of total power consumption.Wherein activity energy expenditure is to be subjected to human body control and change, and greatly then activity energy expenditure of quantity of motion is large, and the little then activity energy expenditure of quantity of motion is little.In the human motion process, because athletic meeting causes human heart rate's increase, the heart rate increase can cause the increase of the oxygen consumption of human body, and the increase of oxygen consumption can affect the increase of the basal energy expenditure of human body.
Usually kinergety consumption measuring method comprises: direct calorimetry, indirect calorimetry, two mark water law, rhythm of the heart method, acceleration transducer method.Direct calorimetry, indirect calorimetry and two mark water law belong to physiology's category, above-mentioned three kinds of methods have higher certainty of measurement, therefore often with the standard of this several method as evaluation energy consumption measurement method, wherein two mark water laws are regarded as the golden standard in the energy consumption measurement, but its measurement device is complicated, the energy expenditure in a certain motor process of description human body that somewhat expensive and measurement result can not be real-time.Heart rate detection method and the measurement of acceleration transducer method are convenient, the energy expenditure of measurement campaign a period of time that can be real-time, but use separately heart rate and acceleration analysis energy expenditure, its accuracy is not high, and is subject to easily self and extraneous interference.
Summary of the invention
The key issue that will solve required for the present invention is to provide a kind of kinergety consumption models in conjunction with acceleration and heart rate, the catabiotic measuring instrument of institute and measuring method in the measurement human motion process.
A kind of sports energy consumption measuring instrument comprises:
Human body personalizing parameters collecting unit be used for to gather and storage human body personalizing parameters, comprises height H(cm), body weight W (kg), age Y (year), sex G and cut-off heart rate FLEX_HR (inferior/as to divide);
Heart rate sensor is used for measuring heart rate HR;
Energy consumption calculation unit, basis is used for utilizing measurement heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Acceleration transducer is used for measuring acceleration of motion;
Motion acting computing unit is used for utilizing acceleration of motion and described body weight W to calculate externally acting EE Act
The sports energy consumption computing unit is with basic energy consumption BEE and external acting EE ActAddition is calculated sports energy consumption EE, i.e. EE=BEE+EE Act
Heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR more than or equal to cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE, heart rate HR and personalizing parameters:
BEE=G×(α 1×W+β 1×Y+χ 1×HR+δ 1)+(1-G)×(α 2×W+β 2×Y+χ 2×HR+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be verification coefficient, α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19.
When measuring heart rate HR less than cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE and human body personalizing parameters:
BEE=G×(α 1×W+β 1×H+χ 1×Y+δ 1)+(1-G)×(α 2×W+β 2×H+χ 2×Y+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient; δ 1, δ 2Be the verification coefficient; α 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
The acceleration that described utilization is measured and body weight W calculate externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
Figure GDA0000370980470000032
In the following formula, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation, V are got in sym () expression 0Be initial velocity.
The present invention also provides a kind of sports energy consumption measuring method, comprising:
Gather the human body personalizing parameters, comprise height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR;
Measure heart rate HR;
Utilize and measure heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Measure acceleration of motion;
Utilize acceleration of motion and described body weight W to calculate externally acting EE Act
With basic energy consumption BEE and the EE that externally does work ActAddition is calculated sports energy consumption EE, i.e. EE=BEE+EE Act
Heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR more than or equal to cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE, heart rate HR and personalizing parameters:
BEE=G×(α 1×W+β 1×Y+χ 1×HR+δ 1)+(1-G)×(α 2×W+β 2×Y+χ 2×HR+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be the verification coefficient, value records by experiment, α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19;
When measuring heart rate HR less than cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE and human body personalizing parameters:
BEE=G×(α 1×W+β 1×H+χ 1×Y+δ 1)+(1-G)×(α 2×W+β 2×H+χ 2×Y+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient; δ 1, δ 2Be the verification coefficient; α 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
Described acceleration of motion and the described body weight W of utilizing calculates externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
Figure GDA0000370980470000042
In the following formula, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation, V are got in sym () expression 0Be initial velocity.
Compared with prior art, the present invention utilizes heart rate sensor to measure heart rate and human body personalizing parameters Calculating Foundation energy consumption, set up the equation of linear regression of basal energy expenditure, heart rate and personalizing parameters, utilize acceleration transducer to measure acceleration of motion, further consider the influence factor of different directions, implement respectively externally acting in conjunction with body weight W in three directions and calculate, the present invention has promoted certainty of measurement with heart rate and the effective combination of acceleration of motion.
Description of drawings
Fig. 1 is sports energy consumption measuring instrument preferred embodiment structured flowchart of the present invention;
Fig. 2 is Calculating Foundation energy consumption preferred embodiment flow chart of the present invention;
Fig. 3 is that heart rate of the present invention and basic energy consumption concern sketch map;
Fig. 4 calculates the preferred embodiment flow chart that done work that moves by the present invention;
Fig. 5 is sports energy consumption measuring method preferred embodiment flow chart of the present invention.
The specific embodiment
In order to make purpose of the present invention, technical scheme and advantage clearer, below in conjunction with drawings and Examples, the present invention is described in further details.
The invention provides a kind of sports energy consumption measuring instrument, preferred embodiment, structured flowchart as shown in Figure 1 mainly comprises:
Human body personalizing parameters collecting unit be used for to gather and storage human body personalizing parameters, comprise height H (cm), body weight W (kg), age Y (year), sex G and cut-off heart rate FLEX_HR (inferior/as to divide);
Heart rate sensor is used for measuring heart rate HR;
Energy consumption calculation unit, basis is used for utilizing measurement heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Acceleration transducer is used for measuring acceleration of motion;
Motion acting computing unit is used for utilizing the acceleration and the described body weight W that measure to calculate externally acting EE Act
The sports energy consumption computing unit is with basic energy consumption BEE and external acting addition EE Act, calculate sports energy consumption EE, i.e. EE=BEE+EE Act
But as a kind of implementation, described human body personalizing parameters collecting unit can be the input terminal that comprises keyboard, display and memorizer, personalizing parameters (comprising height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR) by keyboard input subjects, input for ease of parameter, can be by the relevant information of display simultaneous display input, and be stored in the memorizer, the parameter of calling in the memorizer is used for calculating energy consumption.
But as the another kind implementation, described human body personalizing parameters collecting unit can be the wireless collection device that comprises a bluetooth receiver module and memorizer, the subjects's that bluetooth receiver module reception bluetooth transmitting terminal sends personalizing parameters (comprising height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR), and be stored in the memorizer, the parameter of calling in the memorizer is used for calculating energy consumption.
Art technology is obviously clear, and described human body personalizing parameters collecting unit can also adopt state of the art to realize, enumerates no longer one by one.
Before test, first test or acquisition subjects's personalizing parameters, personalizing parameters (comprising height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR), the maximum heart rate of per minute when described cut-off heart rate FLEX_HR is the human body rest.
Described heart rate sensor gathers human motion per minute heart rate, and is typical, the heart rate sensor HK-08A that can adopt Hefei China section electronics institute to produce.
Described heart rate sensor can be the heart rate of only measuring a subjects in one minute, as measuring heart rate HR;
Described heart rate sensor can be the heart rate of measuring 4 subjectss in one minute, and calculating mean value is as measuring heart rate HR.
The real time acceleration value of described acceleration transducer Real-time Collection human body three directions in motor process, typical, the capacitance acceleration transducer MMA7260 that can adopt Freescale company to produce.Three directions are X-axis, Y-axis and Z axis; typically; described X-axis is the direction of motion; Y-axis is direction of motion left direction, and the Z axis direction of motion is direction vertically upward, and art technology is obviously clear; definition for direction has relative property; can take other to define with reference to the orientation, taking other to define with reference to the orientation does not affect realization of the present invention, should belong to protection domain of the present invention yet yet.
Described basic energy consumption calculation unit Calculating Foundation energy consumption BEE, utilize the heart rate of heart rate sensor collection and the data of human body personalizing parameters input block to calculate, as shown in Figure 2, flow chart for basal energy expenditure in the human body motor process is measured is divided into and heart rate relevant portion and uncorrelated two parts with heart rate.
As a kind of embodiment, as the heart rate HR that measures during more than or equal to cut-off heart rate FLEX_HR, basal energy expenditure and the heart rate of this moment are linear, as shown in Figure 3, heart rate and energy consumption do not have linear relationship during the motion beginning, when heart rate increases to certain value, present linear relationship, set up basal energy expenditure, the equation of linear regression of heart rate and personalizing parameters
BEE=G.(α 1.W+β 1.Y+χ 1.HR+δ 1)+(1-G)(α 2.W+β 2.Y+χ 2.HR+δ 2)
Wherein, G is 1 during the male, and G is 0 during the women; W is subjects's body weight; Y is subjects's age; HR is subjects's heart rate; α 1, β 1, χ 1, α 2, β 2, χ 2Be the correlation coefficient of BEE and subjects's height, body weight, age and heart rate, δ 1, δ 2Be the verification coefficient of test, value records by experiment, wherein α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α wherein 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19.
As another kind of embodiment, as the heart rate HR that measures during less than cut-off heart rate FLEX_HR, the basal energy expenditure of this moment and subjects's heart rate do not have obvious linear relationship, basal energy expenditure when subjects's basal energy expenditure is the human body rest at this moment, can replace with the rest energy expenditure SMR of human body, basal energy expenditure when SMR is sleep quality, only relevant with the personalizing parameters of human body, do not have linear relationship with heart rate, set up the equation of linear regression of basal energy expenditure and human body personalizing parameters
BEE=G.(α 1.W+β 1.H+χ 1.Y+δ 1)+(1-G).(α 2.W+β 2.H+χ 2.Y+δ 2)
Wherein, G is 1 during the male, and G is 0 during the women; W is subjects's quality in the formula; H is subjects's height; Y is subjects's age; α wherein 1, β 1, χ 1, α 2, β 2, χ 2Be respectively the correlation coefficient at BEE and body weight, height and age; δ 1, δ 2Be respectively the verification coefficient of BEE and test, by calculating wherein correlation coefficient and verification coefficient with the indirect thermal measurement method.α wherein 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α wherein 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
As another kind of embodiment, described basic energy consumption BEE can also calculate according to the oxygen consumption in the human motion process, namely comes the Calculating Foundation energy expenditure by the oxygen consumption that gathers in the motor process, sets up basal energy expenditure and oxygen consumption linear relationship:
BEE=α×OC+β
Wherein, the oxygen consumption when OC is motion (liter/min), α is the fitting coefficient between basal energy expenditure and the oxygen consumption, and span is 2.4 – 2.7, and β is the verification coefficient, and span is 0.30-0.5.
The described motion EE that externally does work ActComputing unit calculates the external institute of motion work, utilizes the body weight W of human body personalizing parameters and the acceleration of acceleration transducer collection to calculate.Acceleration on three directions of acceleration transducer collection campaign (being expressed as X-axis, Y-axis and Z axis), record X by sensor, the Y-axis accekeration is the acceleration on this direction in the human body motor process just, Z axis records the impact that acceleration will be subject to gravity acceleration g, therefore must eliminate the effect of g when calculating.
Preferably, it is the product of power and displacement that the present invention adopts the method for calculating merit in the kinesiology, namely
EE act=F×S×cos(θ)
EE in the formula ActBe power institute work, F is the size of power, and S is the displacement on the force direction, and θ is the angle of power and displacement.Because power is consistent with the direction of displacement, therefore when computing power institute work, use respectively
EE act=F×S
S in the formula must calculate first the displacement on the place direction of exerting oneself for what change when calculating merit.S during motion can be expressed as
S = ∫ 0 T ∫ 0 t adndt + V 0 × T
S in the formula is the displacement on the motoricity direction, and T is motion total time, and t is the real time acceleration value that acceleration transducer gathers for the integration period of motion, a, and n is integration variable, V 0Initial velocity during for motion, preferred initial velocity is 0.
EE act = ∫ 0 T ∫ 0 t W × a 2 × ( V 0 + a × t ) dndt
EE in the formula ActBe motion institute work, T is motion total time, and t is the integration period of motion, and W is subjects's body weight, and t is the real time acceleration value that acceleration transducer gathers for the integration period of motion, a, and n is integration variable, V 0Initial velocity during for motion, preferred initial velocity is 0.
Acceleration in the human motion process contains the almost-periodic function of multiple harmonic for take the step frequency as first-harmonic, must take absolute value when calculating integration with such function, otherwise positive negative acceleration can be cancelled out each other to time integral
EE act = ∫ 0 T ∫ 0 t W × a 2 × ( V 0 + | a | × t ) dndt
EE in the formula ActBe motion institute work, T is motion total time, and t is the integration period of motion, and W is subjects's body weight, and t is the real time acceleration value that acceleration transducer gathers for the integration period of motion, a, and n is integration variable, V 0Be initial velocity, preferred initial velocity is 0.
Because not only with to record real-time acceleration relevant on the Z-direction, also the gravity acceleration g with human body is relevant, so will eliminate the gravity acceleration g function influence in calculating for the acceleration on the Z-direction.If record the real-time acceleration of Z axis greater than 0, show that this moment, the gravity center of human body moved upward, total acceleration is a on Z-direction so Z-g, at this moment on the Z-direction acting as shown in the formula
EE act = ∫ 0 T ∫ 0 t W × ( a z - g ) 2 × ( V 0 + | a z - g | × t ) dndt
G in the formula is acceleration of gravity.
If record the Z axis real time acceleration less than 0, the center of gravity that shows human body this moment is that total acceleration is a on Z-direction so downwards Z+ g, at this moment institute's work is as follows on the Z-direction
EE act = ∫ 0 T ∫ 0 t W × ( a z + g ) 2 × ( V 0 + | a z + g | × t ) dndt
In sum, utilize the acceleration and the body weight W that measure to calculate externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
Figure GDA0000370980470000095
Wherein, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation is got in sym () expression, when doing work on calculating X-axis, Y-axis and three directions of Z axis, adopts respectively following formula to carry out, when calculating externally acting EE ActThe time, with on X-axis, Y-axis and three directions of Z axis the addition of doing work and get final product.For example shown in Figure 4, the external preferred embodiment flow chart that does work when moving for calculating, by following formula, utilize the acceleration of body weight and measurement first, calculate externally acting on X-axis and the Y-axis, when the calculating Z axis is done, whether make a decision the Z axis acceleration measurement greater than 0, accordingly the acceleration of measuring is revised, and then calculate Z axis and externally do work, with three direction acting additions, namely obtain externally acting EE ActWhen calculating integration, because the common frequency of human motion is 20HZ, frequency acquisition should be greater than 40HZ, and frequency acquisition of the present invention is 50HZ, and the time cycle that gathers the motion brief acceleration is 0.02s, so the Δ t in the time of will measuring is taken as 0.02s.
Consider the physiologic effect of human body acting, namely the people can feel tired in walking or after a period of time of running, if but come the above action of repetition can feel same fatigue at pedal with identical time and identical stride.The physiologic effect that on the mechanics brow landing motion is produced can be firmly represents the integration of time, namely
EE act = ∫ 0 T | a | d ( t )
EE wherein ActBe motion institute work; T is total time of motion; A is the acceleration of motion that collects; T is integration variable.It is simpler than said method calculating that firmly long-pending mode is calculated the external acting of motion, but the precision of its calculating is lower than said method.
Described sports energy consumption computing unit with the results added of above basic energy consumption calculation unit and motion acting computing unit, can obtain energy consumption total in the motor process within the unit interval.
A kind of sports energy consumption measuring method, as shown in Figure 5, flow process comprises:
Gather the human body personalizing parameters, comprise height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR, describe the embodiment referring to the record of human body personalizing parameters collecting unit in detail, no longer describe in detail.
Measure heart rate HR;
Utilize and measure heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Measure acceleration of motion;
Utilize acceleration of motion and described body weight W to calculate externally acting EE Act
With basic energy consumption BEE and the EE that externally does work ActAddition is calculated sports energy consumption EE, i.e. EE=BEE+EE Act
Heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR more than or equal to cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE, heart rate HR and personalizing parameters:
BEE=G×(α 1×W+β 1×Y+χ 1×HR+δ 1)+(1-G)×(α 2×W+β 2×Y+χ 2×HR+δ 2)
Wherein, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be the verification coefficient, value records by experiment, α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19;
When measuring heart rate HR less than cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE and human body personalizing parameters:
BEE=G×(α 1×W+β 1×H+χ 1×Y+δ 1)+(1-G)×(α 2×W+β 2×H+χ 2×Y+δ 2)
Wherein, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient; δ 1, δ 2Be the verification coefficient; α 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
Detailed description no longer describes in detail referring to the embodiment of basic energy consumption calculation unit record.
Described acceleration of motion and the described body weight W of utilizing calculates externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
Figure GDA0000370980470000112
Wherein, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation, V are got in sym () expression 0Be initial velocity.
Detailed description no longer describes in detail referring to the embodiment of motion acting computing unit record.
Compared with prior art, the present invention utilizes heart rate sensor to measure heart rate and human body personalizing parameters Calculating Foundation energy consumption, set up the equation of linear regression of basal energy expenditure, heart rate and personalizing parameters, utilize acceleration transducer to measure acceleration of motion, further consider the influence factor of different directions, implement respectively externally acting in conjunction with body weight W in three directions and calculate, the present invention has promoted certainty of measurement with heart rate and the effective combination of acceleration of motion.
The present invention further describes the purpose, technical solutions and advantages of the present invention for embodiment or embodiment; institute is understood that; above only is preferred implementation of the present invention for embodiment or embodiment; not in order to limit the present invention; all any modifications of within the spirit and principles in the present invention the present invention being done, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (8)

1. a sports energy consumption measuring instrument is characterized in that, comprising:
Human body personalizing parameters collecting unit is used for gathering and storage human body personalizing parameters, comprises height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR;
Heart rate sensor is used for measuring heart rate HR;
Energy consumption calculation unit, basis is used for utilizing measurement heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Acceleration transducer is used for measuring acceleration of motion;
Motion acting computing unit is used for utilizing acceleration of motion and described body weight W to calculate externally acting EE Act
The sports energy consumption computing unit is with basic energy consumption BEE and external acting EE ActAddition is calculated sports energy consumption EE, i.e. EE=BEE+EE Act
The acceleration that described utilization is measured and body weight W calculate externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
In the following formula, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation, V are got in sym () expression 0Be initial velocity.
2. sports energy consumption measuring instrument as claimed in claim 1, it is characterized in that, described human body personalizing parameters collecting unit is the input terminal that comprises keyboard, display and memorizer, by keyboard input human body personalizing parameters, the information that display simultaneous display input is relevant, with the Parameter storage of input in memorizer.
3. sports energy consumption measuring instrument as claimed in claim 1, it is characterized in that, described human body personalizing parameters collecting unit is the wireless collection device that comprises bluetooth receiver module and memorizer, the bluetooth receiver module receives the human body personalizing parameters that the bluetooth transmitting terminal sends, with the Parameter storage that receives in memorizer.
4. sports energy consumption measuring instrument as claimed in claim 1 is characterized in that, described heart rate sensor was measured the heart rate of 4 human bodies in one minute, and calculating mean value is as measuring heart rate HR.
5. sports energy consumption measuring instrument as claimed in claim 1 is characterized in that, heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR more than or equal to cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE, heart rate HR and personalizing parameters:
BEE=G×(α 1×W+β 1×Y+χ 1×HR+δ 1)+(1-G)×(α 2×W+β 2×Y+χ 2×HR+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be verification coefficient, α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19.
6. sports energy consumption measuring instrument as claimed in claim 1 is characterized in that, heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR less than cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE and human body personalizing parameters:
BEE=G×(α 1×W+β 1×H+χ 1×Y+δ 1)+(1-G)×(α 2×W+β 2×H+χ 2×Y+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be the verification coefficient; α 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
7. a sports energy consumption measuring method is characterized in that, comprising:
Gather the human body personalizing parameters, comprise height H, body weight W, age Y, sex G and cut-off heart rate FLEX_HR;
Measure heart rate HR;
Utilize and measure heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption BEE;
Measure acceleration of motion;
Utilize acceleration of motion and described body weight W to calculate externally acting EE Act
With basic energy consumption BEE and the EE that externally does work ActAddition is calculated sports energy consumption EE, i.e. EE=BEE+EE Act
Described acceleration of motion and the described body weight W of utilizing calculates externally acting EE ActBy calculating the sum of being done work on X-axis, Y-axis and three directions of Z axis, doing work on described calculating X-axis, Y-axis and three directions of Z axis is respectively:
EE xyz = ∫ 0 T ∫ 0 t W × ACC 2 × ( V 0 + ACC × t ) dndt
Figure FDA0000370980460000032
In the following formula, a XBe the acceleration that X-axis is measured, a YBe the acceleration that Y-axis is measured, a ZBe the acceleration that Z axis is measured, g is acceleration of gravity, and symbolic operation, V are got in sym () expression 0Be initial velocity.
8. sports energy consumption measuring method as claimed in claim 7 is characterized in that, heart rate HR and described human body personalizing parameters Calculating Foundation energy consumption are measured in described utilization, comprising:
When measuring heart rate HR more than or equal to cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE, heart rate HR and personalizing parameters:
BEE=G×(α 1×W+β 1×Y+χ 1×HR+δ 1)+(1-G)×(α 2×W+β 2×Y+χ 2×HR+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient, δ 1, δ 2Be the verification coefficient, value records by experiment, α 1Span be 0.15-0.24, β 1Span be 0.2-0.3, χ 1Span be 0.6-0.7, δ 1Span be 55-60; α 2Span be-0.2--0.1 β 2Span be 0.06-0.075, χ 2Span be 0.44-0.50, δ 2Span be-21--19;
When measuring heart rate HR less than cut-off heart rate FLEX_HR, set up the equation of linear regression of basic energy consumption BEE and human body personalizing parameters:
BEE=G×(α 1×W+β 1×H+χ 1×Y+δ 1)+(1-G)×(α 2×W+β 2×H+χ 2×Y+δ 2)
In the following formula, G is 1 during the male, and G is 0 during the women; α 1, β 1, χ 1, α 2, β 2, χ 2Be correlation coefficient; δ 1, δ 2Be the verification coefficient; α 1Span be 13-14, β 1Span be 4.9-5.1, χ 1Span be-7.0--6.7 δ 1Span be 60-70; α 2Span be 9.5-10, β 2Span be 1.7-2.0, χ 2Span be-4.9--4.5 δ 2Span be 650-670.
CN 201210249466 2012-07-18 2012-07-18 Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method Active CN102727185B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210249466 CN102727185B (en) 2012-07-18 2012-07-18 Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210249466 CN102727185B (en) 2012-07-18 2012-07-18 Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method

Publications (2)

Publication Number Publication Date
CN102727185A CN102727185A (en) 2012-10-17
CN102727185B true CN102727185B (en) 2013-10-30

Family

ID=46984018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210249466 Active CN102727185B (en) 2012-07-18 2012-07-18 Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method

Country Status (1)

Country Link
CN (1) CN102727185B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI559901B (en) * 2014-04-14 2016-12-01 Quanta Comp Inc Method and device of sleep detection
CN104111978B (en) * 2014-06-25 2017-08-29 京东方科技集团股份有限公司 Energy consumption measurement method and energy consumption measurement system
CN106934746A (en) * 2015-12-29 2017-07-07 上海歌志网络科技有限公司 A kind of health and fitness facilities motion consumes the computational methods of calorie
CN106037738A (en) * 2016-07-11 2016-10-26 中国人民解放军第三军医大学第三附属医院 Wearable physical training evaluation system and device
CN106551684A (en) * 2017-02-04 2017-04-05 深圳市哎优智能科技有限公司 A kind of trunk motion exercise efficiency estimation method and device
CN109222945B (en) * 2017-07-10 2023-09-29 博世汽车部件(苏州)有限公司 Systems, methods, devices, and apparatus for fitness monitoring in a vehicle
WO2022221987A1 (en) * 2021-04-19 2022-10-27 中国科学院深圳先进技术研究院 Energy metabolism evaluation method and system based on wearable sensing information fusion
CN117412699A (en) * 2022-05-13 2024-01-16 北京小米移动软件有限公司 Heat consumption estimation method, device and storage medium

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6554776B1 (en) * 2001-11-21 2003-04-29 Medical Graphics Corporation Method for detecting anaerobic threshold and prescribing a training zone to maximize fat utilization or improved cardiovascular fitness
CN1723838A (en) * 2005-07-21 2006-01-25 高春平 Method and device for individualized and three-D type reducing weight
US7643873B2 (en) * 2006-04-28 2010-01-05 Idt Technology Limited Exercise data apparatus
JP4521040B2 (en) * 2008-02-28 2010-08-11 株式会社シマノ Calorie consumption measuring device and calorie consumption measuring method
US8182424B2 (en) * 2008-03-19 2012-05-22 Microsoft Corporation Diary-free calorimeter
CN101518442B (en) * 2008-07-05 2011-04-20 杭州义盛祥通信技术有限公司 Sport quantization watch and sport quantitative analysis method
CN102512174B (en) * 2011-11-21 2014-10-01 中国人民解放军第四军医大学 Portable type physical ability consumption and physiological parameter monitoring instrument

Also Published As

Publication number Publication date
CN102727185A (en) 2012-10-17

Similar Documents

Publication Publication Date Title
CN102727185B (en) Motion energy consumption measuring apparatus based on heart rate and acceleration as well as measuring method
US10098549B2 (en) Local model for calorimetry
TWI393579B (en) The state of the muscle movement state analysis system, methods and computer program products
CN105797350A (en) Intelligent method and system for body building posture recognition, evaluation, early-warning and intensity estimation
CN102292025A (en) Determining energy expenditure of a user
CN102680025A (en) Indoor thermal comfort evaluation system
CN103206962A (en) Electric pedometer based on approaching type sensing
TW201233977A (en) Portable evaluator of amount of exercise and method of evaluating amount of exercise thereof
CN103654798B (en) Method and device for monitoring and recording emotion
CN204496173U (en) A kind of intelligent watch for counting of skipping rope
CN104173052B (en) A kind of human comfort measuring device based on RULA and method thereof
CN103801048A (en) Intelligent skipping rope and control method
CN108981745A (en) A kind of step-recording method, device, equipment and storage medium
US20160058373A1 (en) Running Energy Efficiency
CN105997043A (en) Pulse frequency extraction method based on wrist type wearable device
CN102440781A (en) Method and device for measuring exercise amount
CN105771224B (en) A kind of locomotion evaluation system based on multisensor
CN107966161A (en) Walking detection method based on FFT
Fan et al. Wearable motion attitude detection and data analysis based on Internet of Things
TWM478203U (en) Fitness system with physical ability diagnosis function
CN201683882U (en) Human body physical health detection and fitness tester
CN104138658A (en) Motion measuring and animation feedback system based on smart shoes
CN205120130U (en) FFT calculates equipment of step number and human motion consumption calorie based on improve
CN102243248B (en) Speed detection device and method capable of automatically differentiating walking and running modes
CN107490470A (en) A kind of detection method for upper limbs ectoskeleton power-assisted efficiency

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant