CN102750441B - 一种骑自行车所消耗卡路里的计算方法 - Google Patents
一种骑自行车所消耗卡路里的计算方法 Download PDFInfo
- Publication number
- CN102750441B CN102750441B CN201210179085.5A CN201210179085A CN102750441B CN 102750441 B CN102750441 B CN 102750441B CN 201210179085 A CN201210179085 A CN 201210179085A CN 102750441 B CN102750441 B CN 102750441B
- Authority
- CN
- China
- Prior art keywords
- bicycle
- formula
- bike
- consumed
- calorie
- 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.)
- Expired - Fee Related
Links
Landscapes
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了一种骑自行车所消耗卡路里的计算方法,其特征是通过参数标定的方法获得自行车和骑车人的总重量、自行车和骑车人在骑行中所受风阻的受力面积以及自行车阻力参数,然后计算单位时间内骑自行车时所做的功包括有:克服风的阻力所做的功、克服人和自行车的重力所做的功、改变人和自行车的动能所做的功,以及克服自行车阻力所做的功,从而计算出骑自行车所消耗卡路里。本发明能根据不同运动环境下,实时计算出人体所消耗的卡路里值,计算结果更加精确可靠。
Description
技术领域
本发明涉及一种卡路里计算方法,具体地说是一种骑自行车时所消耗卡路里的计算方法。
背景技术
越来越多的人将自行车作为户外健身设备,为了实现更专业化的健身效果,需要计算骑自行车过程中人体所消耗的卡路里情况。但市场上相关的电子设备关于骑自行车所消耗卡路里量,大多是采用基于行驶的路程的方法或者基于检测人体身体参数的方法来计算。
基于行驶路程来计算骑自行车消耗卡路里量的方法,在计算过程中采用的是估算的方式,即粗略的将消耗的卡路里量和行驶的距离对应起来。这样的计算方法忽略了人的体重的差异、行驶路面不同以及不同环境下风的阻力所带来的影响,也忽略了自行车本身的旧损所带来的影响,计算起来较为简单,但无法满足精确计算的要求。
基于检测人体身体参数来计算骑自行车消耗卡路里量的方法,需要将检测装置穿在身上或者戴在手上。这样会使运动者有不舒服的感觉,影响健身效果和个人运动情绪;而且有些设备对人体会造成一定的辐射伤害,不利于舒适健康的户外健身。
发明内容
本发明是为了克服现有技术的不足之处,提供一种骑自行车所消耗卡路里的计算方法,通过参数标定的方法确定在不同运动环境下人体所消耗的卡路里值,使得计算结果更加精确可靠。
本发明为解决技术问题采用如下技术方案:
本发明一种骑自行车所消耗卡路里的计算方法按如下过程进行:
t0为骑自行车卡路里消耗的初始时刻,tn为骑自行车卡路里消耗的结束时刻,在初始时刻t0到结束时刻tn时间段内取时间间隔Δt为ti到ti+1的单位时间,ti表示一个时间点,ti+1表示下一个时间点,i=0,1,...,n;n为整数,取值为 是检测获得的ti时刻自行车的行驶速度;是检测获得的ti+1时刻自行车的行驶速度;是检测获得的ti时刻行驶中人和自行车所受风阻的压强;是检测获得的ti时刻自行车行驶的路面与水平面的夹角;M是自行车和骑车人的总重量;S是自行车和骑车人在骑行中所受风阻的受力面积;ψ是自行车阻力参数;g是重力加速度;
Wpi为ti到ti+1的单位时间Δt内克服风的阻力所做的功,由式(1)计算获得:
Wgi为ti到ti+1的单位时间Δt内克服人和自行车的重力所做的功,由式(2)计算获得:
Wdi为ti到ti+1的单位时间Δt内改变人和自行车的动能所做的功,由式(3)计算获得:
Wfi为假设自行车在水平面上行驶时,车速保持固定且在没有风阻的情况下,ti到ti+1的单位时间Δt内克服自行车阻力所做的功,由式(4)计算获得:
则在ti到ti+1的单位时间Δt内骑自行车所消耗的卡路里为Wki,按式(5)计算获得:
Wki=Wpi+Wgi+Wdi+Wfi (5)
假设:
由式(1)、式(2)、式(3)、式(4)、式(5)、式(6)、式(7)、式(8)和式(9)获得式(10):
Wki=S×wpi+M×(wgi+wdi)+ψ×wfi (10)
则骑自行车所消耗的卡路里Wk按照式(11)计算获得:
本发明一种骑自行车所消耗卡路里的计算方法的特点也在于:
所述自行车和骑车人的总重量M、所述自行车和骑车人在骑行中所受风阻的受力面积S以及所述自行车阻力参数ψ是通过参数标定的方法按照如下步骤获得:
在骑自行车过程中,骑车人停止向自行车做功,则骑自行车所消耗的卡路里Wkj=0(j=0,1,...,m),m为大于等于3的整数;利用式(10)每隔一段单位时间Δt内得出一个骑自行车所消耗的卡路里等式,获得如下方程组:
S×wp1+M×(wg1+wd1)+ψ×wf1=Wk1=0 方程1
S×wp2+M×(wg2+wd2)+ψ×wf2=Wk2=0 方程2
S×wpm+M×(wgm+wdm)+ψ×wfm=Wkm=0 方程m
利用所述方程组,采用最小二乘法分别计算获得所述自行车和骑车人的总重量M、所述自行车和骑车人在骑行中所受风阻的受力面积S和所述自行车阻力参数ψ。
与已有技术相比,本发明有益效果体现在:
1、本发明从人骑自行车所消耗卡路里的能量转换角度出发,分析出骑自行车时所做的功包括有:克服风的阻力所做的功、克服人和自行车的重力所做的功、改变人和自行车的动能所做的功,以及克服自行车阻力所做的功,通过分别计算出卡路里消耗的各个分量,使卡路里消耗量的计算结果更为合理准确。
2、本发明通过参数标定方法能根据不同运动环境中的不同行车参数,以及运动者本身的体能消耗情况和运动速度的变化情况,实时计算出卡路里消耗的值的大小,从而保证计算的准确性和时效性。
具体实施方式
本实施例中,骑自行车所消耗卡路里的计算方法按如下过程进行:
t0为骑自行车卡路里消耗的初始时刻,tn为骑自行车卡路里消耗的结束时刻,在初始时刻t0到结束时刻tn时间段内取时间间隔Δt为ti到ti+1的单位时间,ti表示一个时间点,ti+1表示下一个时间点,i=0,1,...,n,n为整数,取值为具体实施中单位时间Δt取值为0.1s;
1、本计算方法所用的参数包括:
a、自行车的行驶速度
是检测获得的ti时刻自行车的行驶速度;是检测获得的ti+1时刻自行车的行驶速度;具体实施中采用的是12V/6W的摩擦发电机,摩擦发电机自身体积小,质量轻,安装方便,可以安装在自行车前轮或者后轮的轴上,紧贴自行车轮。摩擦发电机的周长为23cm,每旋转一周能形成3个正弦信号,通过光耦将正弦信号转换脉冲信号,通过计算单位时间Δt内检测到的脉冲信号的个数来计算自行车的行驶速度当单位时间Δt为0.1s时,初始时刻ti到结束时刻ti+1时间段内检测到的脉冲个数为l,则该时间段内的自行车的行驶速度可以按式(1)计算:
b、行驶中人和自行车所受风阻的压强
压强是检测获得的ti时刻自行车行驶中所受风阻力的压强,具体实施中采用气体压力传感器进行检测;
c、自行车行驶的路面与水平面的夹角
是检测获得的ti时刻自行车行驶的路面与水平面的夹角,具体实施中采用三轴陀螺仪传感器进行检测;
d、重力加速度g;
具体实施中取值为9.8(m/s2);
e、自行车和骑车人的总重量M;
f、自行车和骑车人在骑行中所受风阻的受力面积S;
g、自行车阻力参数ψ;
具体实施中参数e、f和g通过参数标定的方式获得。
2.本计算方法所用的计算公式包括:
Wpi为单位时间Δt内克服风阻所做的功,是由所受风的阻力和行驶位移的乘积按照式(2)计算获得:
式(2)中,在单位时间Δt为0.1s时,骑车人所受的风的阻力为行驶的距离为
Wgi为单位时间Δt内克服重力所做的功,是由人和自行车所受的重力乘以人和自行车在竖直方向上行驶的高度按照式(3)计算获得:
式(3)中,人和自行车所受的重力为M×g,在单位时间Δt为0.1s时,人和自行车在竖直方向上行驶的高度为
Wdi为单位时间Δt内改变动能所做的功;是由人和自行车所受的重力乘以人和自行车速度的平方差的一半按照式(4)计算获得:
Wfi为单位时间Δt内克服自行车阻力所做的功,是在理想状态下,假设人骑自行车时不受风的阻力、不改变人和自行车的动能、不改变人和自行车势能(即不克服重力做功)的情况下,单纯克服自行车阻力所做的功,这种功只与路程以及自行车阻力参数ψ有关,其中自行车阻力参数ψ会因人的重量差异、自行车的不同、自行车的磨损情况差异及路面的粗糙程度的不同而不同,但一旦骑自行车的人固定、自行车固定、行驶在固定的路面上,则自行车阻力参数ψ就能确定下来,按照(5)计算获得:
Wki为单位时间Δt内骑自行车所消耗的卡路里,按式(6)计算获得:
Wki=Wpi+Wgi+Wdi+Wfi (6)
式(6)从人骑自行车所消耗卡路里的能量转换角度出发,分析了卡路里消耗主要用于克服风的阻力做过、克服人和自行车的重力做功、改变人和自行车的动能做功以及克服自行车的阻力做功。因此在单位时间Δt内骑自行车所消耗的卡路里为该时间段里各个做功项的求和。
假设:
由式(2)、式(3)、式(4)、式(5)、式(6)、式(7)、式(8)、式(9)和式(10)获得式(11):
Wki=S×wpi+M×(wgi+wdi)+ψ×wfi (11)
式(11)中,单位时间Δt内,wpi、wgi、wdi和wfi可以根据各传感器实时采集到的行驶速度压强和夹角的值分别计算获得;
当在一个固定的环境下,即骑车人固定、自行车固定、行驶在相同的路面上,则行车参数总重量M、受力面积S和自行车阻力参数ψ就可以确定下来;从而计算ti到ti+1的单位时间间隔内的骑自行车消耗的卡路里Wki;
则骑自行车所消耗的卡路里Wk按照式(12)计算获得:
3.本计算方法的参数标定的方法按如下过程进行:
骑自行车过程中,骑车人停止向自行车做功,即让骑自行车所消耗的卡路里Wkj=0(j=0,1,...,m),m为大于等于3的整数;利用式(11)每隔一段单位时间Δt得出一个骑自行车所消耗的卡路里等式,获得如下方程组:
S×wp1+M×(wg1+wd1)+ψ×wf1=Wk1=0 方程1
S×wp2+M×(wg2+wd2)+ψ×wf2=Wk2=0 方程2
S×wpm+M×(wgm+wdm)+ψ×wfm=Wkm=0 方程m
利用上述方程组,采用最小二乘法分别计算获得所述自行车和骑车人的总重量M、所述人骑自行车过程中所受风阻力的受力面积S和所述自行车阻力参数ψ。
参数标定模式是为了检测不同的人在不同骑车环境中的行车参数,如:骑车人的重量差异、骑车人的穿衣差别,自行车的不同、自行车的磨损情况差异以及路面的粗糙程度的不同而导致行车参数各不相同,但在一个固定的环境,即骑车人固定、自行车固定、行驶在相同的路面上,则行车参数总重量M、受力面积S和自行车阻力参数ψ就可以确定下来,从而计算获得单位时间Δt内的卡路里消耗量Wki;
而当骑车环境改变时,例如道路改变等,只要停止向自行车做功,即让卡路里消耗量Wkj=0;则得到m个方程,通过最小二乘法来求解固定的方程组,就可以重新标定出行车参数总重量M、受力面积S和自行车阻力参数ψ,从而保证了不同环境下,计算卡路里消耗的准确性。
Claims (2)
1.一种骑自行车所消耗卡路里的测算方法,其特征按如下步骤进行:
t0为骑自行车卡路里消耗的初始时刻,tn为骑自行车卡路里消耗的结束时刻,在初始时刻t0到结束时刻tn时间段内取时间间隔Δt为ti到ti+1的单位时间,ti表示一个时间点,ti+1表示下一个时间点,i=0,1,...,n;n为整数,取值为 是检测获得的ti时刻自行车的行驶速度;是检测获得的ti+1时刻自行车的行驶速度;是检测获得的ti时刻行驶中人和自行车所受风阻的压强;是检测获得的ti时刻自行车行驶的路面与水平面的夹角;M是自行车和骑车人的总重量;S是自行车和骑车人在骑行中所受风阻的受力面积;ψ是自行车阻力参数;g是重力加速度;
Wpi为ti到ti+1的单位时间Δt内克服风的阻力所做的功,由式(1)计算获得:
Wgi为ti到ti+1的单位时间Δt内克服人和自行车的重力所做的功,由式(2)计算获得:
Wdi为ti到ti+1的单位时间Δt内改变人和自行车的动能所做的功,由式(3)计算获得:
Wfi为假设自行车在水平面上行驶时,车速保持固定且在没有风阻的情况下,ti到ti+1的单位时间Δt内克服自行车阻力所做的功,由式(4)计算获得:
则在ti到ti+1的单位时间Δt内骑自行车所消耗的卡路里为Wki,按式(5)计算获得:
Wki=Wpi+Wgi+Wdi+Wfi (5)
假设:
由式(1)、式(2)、式(3)、式(4)、式(5)、式(6)、式(7)、式(8)和式(9)获得式(10):
Wki=S×wpi+M×(wgi+wdi)+ψ×wfi (10)
则骑自行车所消耗的卡路里Wk按照式(11)计算获得:
2.根据权利要求1所述的骑自行车所消耗卡路里的测算方法,其特征是:
所述自行车和骑车人的总重量M、所述自行车和骑车人在骑行中所受风阻的受力面积S以及所述自行车阻力参数ψ是通过参数标定的方法按照如下步骤获得:
在骑自行车过程中,骑车人停止向自行车做功,则骑自行车所消耗的卡路里Wkj=0(j=0,1,...,m),m为大于等于3的整数;利用式(10)每隔一段单位时间Δt内得出一个骑自行车所消耗的卡路里等式,获得如下方程组:
S×wp1+M×(wg1+wd1)+ψ×wf1=Wk1=0 方程1
S×wp2+M×(wg2+wd2)+ψ×wf2=Wk2=0 方程2
·
·
·
S×wpm+M×(wgm+wdm)+ψ×wfm=Wkm=0 方程m
利用所述方程组,采用最小二乘法分别计算获得所述自行车和骑车人的总重量M、所述自行车和骑车人在骑行中所受风阻的受力面积S和所述自行车阻力参数ψ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210179085.5A CN102750441B (zh) | 2012-06-01 | 2012-06-01 | 一种骑自行车所消耗卡路里的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210179085.5A CN102750441B (zh) | 2012-06-01 | 2012-06-01 | 一种骑自行车所消耗卡路里的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102750441A CN102750441A (zh) | 2012-10-24 |
CN102750441B true CN102750441B (zh) | 2015-07-01 |
Family
ID=47030617
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210179085.5A Expired - Fee Related CN102750441B (zh) | 2012-06-01 | 2012-06-01 | 一种骑自行车所消耗卡路里的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102750441B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3158928B1 (en) * | 2015-05-21 | 2021-07-07 | bOMDIC Inc. | Stamina monitoring method and device |
CN106644208B (zh) * | 2017-01-05 | 2019-06-14 | 武汉齐物科技有限公司 | 一种骑行能力分析***及分析方法 |
CN108446432B (zh) * | 2018-02-06 | 2021-12-17 | 浙江工业大学 | 一种基于模型的虚拟自行车骑手骑速的计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201032488Y (zh) * | 2007-02-02 | 2008-03-05 | 桂盟企业股份有限公司 | 车辆用全球定位装置 |
CN201218740Y (zh) * | 2008-04-02 | 2009-04-08 | 朱培辉 | 自行车运动辅助表 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005054872A2 (en) * | 2003-11-25 | 2005-06-16 | Julian Gordon | Method and device for determination of airspeed in human-powered motion, and methods for accurate calibration thereof |
US20070156335A1 (en) * | 2006-01-03 | 2007-07-05 | Mcbride Sandra Lynn | Computer-Aided Route Selection |
-
2012
- 2012-06-01 CN CN201210179085.5A patent/CN102750441B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201032488Y (zh) * | 2007-02-02 | 2008-03-05 | 桂盟企业股份有限公司 | 车辆用全球定位装置 |
CN201218740Y (zh) * | 2008-04-02 | 2009-04-08 | 朱培辉 | 自行车运动辅助表 |
Also Published As
Publication number | Publication date |
---|---|
CN102750441A (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Martin et al. | Validation of a mathematical model for road cycling power | |
CN106394561B (zh) | 一种车辆的纵向车速的估计方法和装置 | |
GB2493476A (en) | Real-time calculation of total longitudinal force and aerodynamic drag acting on a rider on a vechile | |
CN102951158B (zh) | 车辆质量估算方法 | |
CN103454442B (zh) | 用于铁路车辆的行驶速度补偿装置及其方法 | |
CN103909933B (zh) | 一种分布式电驱动车辆的前轮侧向力估算方法 | |
CN105109490B (zh) | 一种基于三轴加速度传感器判断车辆急转弯的方法 | |
CN110239554A (zh) | 一种车载路面纵向坡度实时识别方法及其装置 | |
CN102750441B (zh) | 一种骑自行车所消耗卡路里的计算方法 | |
Mannion et al. | Analysis of crosswind aerodynamics for competitive hand-cycling | |
CN108510737A (zh) | 一种融合风环境的无人驾驶车辆电源实时监控方法及装置 | |
CN106644208A (zh) | 一种骑行能力分析***及分析方法 | |
D'Hooge et al. | Application of real-world wind conditions for assessing aerodynamic drag for on-road range prediction | |
CN102991505A (zh) | 一种前方道路曲率估计装置及其估计方法 | |
CN204188244U (zh) | 一种车辆动态称重装置 | |
CN105270409A (zh) | 一种路面峰值附着系数测试装置及测试方法 | |
CN102688020B (zh) | 一种骑自行车所消耗卡路里的动态测量*** | |
CN114186362A (zh) | 一种基于运营列车振动响应预测轨道不平顺及其传感器优化配置的分析方法 | |
AU2021102801A4 (en) | Velocity estimation from cycle pedaling by using extended kalman filter algorithm | |
ATE303581T1 (de) | Fahrrad-neigungsmesser mit beschleunigungssensor und raddrehzahl-sensor zur berechnung der fahrbahnneigung, der bergauf erzeugten energie, der höhenlage und der pedaltrittfrequenz | |
Crouch et al. | A phase-averaged analysis of the pedalling cyclist wake | |
CN110260094A (zh) | 行走式管道长度测量装置 | |
CN109990766A (zh) | 基于路面综合阻力系数的道路坡度检测方法与*** | |
US20170138807A1 (en) | Method of measuring motion power and resistance coefficient of bicycle | |
CN103454022A (zh) | 一种骑行自行车做功测量***及方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150701 |