CN103942445A - 一种健身跑运动方案生成方法 - Google Patents
一种健身跑运动方案生成方法 Download PDFInfo
- Publication number
- CN103942445A CN103942445A CN201410178014.2A CN201410178014A CN103942445A CN 103942445 A CN103942445 A CN 103942445A CN 201410178014 A CN201410178014 A CN 201410178014A CN 103942445 A CN103942445 A CN 103942445A
- Authority
- CN
- China
- Prior art keywords
- jog
- heart rate
- population
- fitness
- value
- 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.)
- Granted
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明涉及一种健身跑运动方案生成方法,包括以下步骤:S1.采集运动者健身跑过程中的速度和心率数据;S2.基于步骤S1采集的速度和心率数据,运用遗传算法对运动者的健身跑运动模型进行参数整定,建立运动者的健身跑运动模型;S3.基于步骤S2建立的健身跑运动模型,以及给定运动时间和安全有效心率区间参数,以在给定运动时间的健身跑过程中,心率落在安全有效心率区间内的时长最长为优化目标,建立运动者的健身跑最优运动方案。该方法可以为运动者建立个性化的健身跑最优运动方案,从而为指导运动者进行安全、有效的健身跑提供了有用的技术手段。
Description
技术领域
本发明涉及健身跑运动优化技术领域,特别是一种健身跑运动方案生成方法。
背景技术
健身跑是一种最为普及的健身和体疗方法,健身跑运动过程建模、优化等相关课题一直是学者们的研究热点。但是,目前定量研案健身跑方面的成果还比较少。Cheng(2008)和Scalzi(2012)提出了健身跑的运动过程模型,在此基础上,Scalzi(2013)等还给出自动生成最优跑步方案的方法。然而Cheng和Scalzi模型中的参数只是针对23–38岁的英国人,Scalzi等人的最优方案也只能用于赛跑的情景,即在指定运动距离下,得到最短时间。这些成果难以直接用于指导健身跑的教学和平时锻炼,这使得健身跑的体育教学和平时锻炼缺乏有效的技术手段和工具。
发明内容
本发明的目的在于提供一种健身跑运动方案生成方法,该方法可以为运动者建立个性化的健身跑最优运动方案,从而为指导运动者进行安全、有效的健身跑提供了有用的技术手段。
为实现上述目的,本发明的技术方案是:一种健身跑运动方案生成方法,包括以下步骤:
S1. 采集运动者健身跑过程中的速度和心率数据;
S2. 基于步骤S1采集的速度和心率数据,运用遗传算法对运动者的健身跑运动模型进行参数整定,建立运动者的健身跑运动模型;
S3. 基于步骤S2建立的健身跑运动模型,以及给定运动时间和安全有效心率区间参数,以在给定运动时间的健身跑过程中,心率落在安全有效心率区间内的时长最长为优化目标,建立运动者的健身跑最优运动方案。
进一步的,在步骤S2中,对于如式(1)的健身跑运动模型:
(1)
其中,x 1(t)、x 2(t)和u(t)分别表示运动者t时刻较0时刻的心率增量、t时刻较0时刻的疲劳程度增量和t时刻的速度,、分别表示x 1(t)、x 2(t)的微分,a 1、a 2、a 3、a 4和a 5为5个待估参数;按如下步骤建立运动者的健身跑运动模型:
S201. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群P(t),t=0,其中初始化参数a 1、a 2、a 3、a 4和a 5,个体编码采用实数编码,编码长度为5;
S202. 适应度计算:以心率的估计值与实测值的残差平方和作为适应度函数,根据式(1)和采集的速度实测值u(t),计算心率估计值x 1(t),将计算结果和采集的心率实测值输入适应度函数,计算各个体的适应度;
S203. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S204. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S203得到的下一代种群,返回步骤S202进行下一次迭代;
将求得的参数a 1、a 2、a 3、a 4和a 5的值代入式(1),得到运动者的健身跑运动模型。
进一步的,在步骤S3中,基于步骤S2建立的健身跑运动模型,将求解运动者的健身跑最优运动方案抽象为一个带约束的优化问题:在给定运动时间Ts和采样周样T下,在速度值u(t)中找出最优的速度值u*(t),使其满足健身跑过程中心率不超过预设的最大值V max,且心率达到安全有效心率区间[a, b]的次数N(t)最多;u(t)和N(t)的定义如式(2)和式(3):
(2)
(3)
其中,表示对T S 除于T的值向下取整,对应T S 时间内的采样次数,s i (t)表示在第i个采样周期内其速度值为1,而在其它采样周期内其速度值为0;c i 为待估参数,i=1,2,…,;按如下步骤估计个参数c i 的值:
S301. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群Q(t),t=0,初始化参数,个体编码采用实数编码,编码长度为;
S302. 适应度计算:以式(3)作为适应度函数,首先估计速度值u(t),然后将该值代入健身跑运动模型,得到心率值x 1(t),最后将计算结果输入适应度函数,计算各个体的适应度,对于x 1(t)的值超过V max的个体,将其适应度赋值为0;
S303. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S304. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S303得到的下一代种群,返回步骤S302进行下一次迭代;
将求得的参数的值代入式(2),得到运动者的健身跑最优运动方案。
本发明的有益效果是:
(1)本发明基于运动者健身跑运动过程采集的速度与心率数据,构建了健身跑的个性化运动过程模型,从而准确刻画健身跑运动过程中个体的各种生理反映与其心率及速度之间的非线性关系,为制订个性化的健身跑运动方案奠定基础。
(2)本发明基于健身跑的个性化运动过程模型,并在指定运动时间和保障安全有效运动的条件下,建立了运动者的健身跑最优运动方案,该最优运动方案可以指导运动者进行科学、安全和有效的健身跑运动,对于体质和心理健康将起着积极促进作用,具有较好的应用和推广价值。
附图说明
图1是本发明实施例的实现流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步说明。
本发明的健身跑运动方案生成方法,如图1所示,包括以下步骤:
S1. 采集运动者在健身跑运动过程中的速度和心率数据。
S2. 基于步骤S1采集的速度和心率数据,运用遗传算法对Cheng和Scalzi等定义的健身跑运动模型进行参数整定,建立运动者的健身跑运动模型。具体为:
对于如式(1)的Cheng和Scalzi定义的健身跑运动模型(微分方程组):
(1)
其中,x 1(t)、x 2(t)和u(t)分别表示运动者t时刻较0时刻(平静状态下)的心率增量、t时刻较0时刻的疲劳程度增量(模拟复杂的生理过程)和t时刻的速度,、分别表示x 1(t)、x 2(t)的微分,a 1、a 2、a 3、a 4和a 5为5个待估参数;具体按如下步骤建立运动者的健身跑运动模型:
S201. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群P(t),t=0,其中初始化参数a 1、a 2、a 3、a 4和a 5,个体编码采用实数编码,编码长度为5;
S202. 适应度计算:以心率的估计值与实测值的残差平方和作为适应度函数(适应值越小,个体越优),根据式(1)和采集的速度实测值u(t),采用四阶龙格-库塔法计算心率估计值x 1(t),将计算结果和采集的心率实测值输入适应度函数,计算各个体的适应度;
S203. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S204. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S203得到的下一代种群,返回步骤S202进行下一次迭代;
将求得的参数a 1、a 2、a 3、a 4和a 5的值代入式(1),得到运动者的健身跑运动模型。
S3. 基于步骤S2建立的健身跑运动模型,以及给定运动时间和安全有效心率区间参数,以在给定运动时间的健身跑过程中,心率落在安全有效心率区间内的时长最长为优化目标,建立运动者的健身跑最优运动方案。具体为:
基于步骤S2建立的健身跑运动模型,将求解运动者的健身跑最优运动方案抽象为一个带约束的优化问题:在给定运动时间Ts和采样周样T下,在速度值u(t)中找出最优的速度值u*(t),使其满足健身跑过程中心率不超过预设的最大值V max,且心率达到安全有效心率区间[a, b]的次数N(t)最多;u(t)和N(t)的定义如式(2)和式(3):
(2)
(3)
其中,表示对T S 除于T的值向下取整,对应T S 时间内的采样次数,u(t)为分段阶跃函数,它反映了在一段足够小的采样周期中跑步的速度是一个定值,它由一些特殊的分段阶跃函数s i (t)进行线性组合得到;其中,s i (t)表示在第i个采样周期内其速度值为1,而在其它采样周期内其速度值为0;c i 为待估参数,有个参数c i 需要估计,即i=1,2,…,;具体按如下步骤估计个参数c i 的值:
S301. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群Q(t),t=0,初始化参数,个体编码采用实数编码,编码长度为;
S302. 适应度计算:以式(3)作为适应度函数(适应值越大,个体越优),首先估计速度值u(t),然后将该值代入到运动者的健身跑运动模型中,得到心率值x 1(t),最后将计算结果输入适应度函数,计算各个体的适应度,对于x 1(t)的值超过V max的个体,将其适应度赋值为0;
S303. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S304. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S303得到的下一代种群,返回步骤S302进行下一次迭代;
将求得的参数的值代入式(2),得到运动者的健身跑最优运动方案。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (3)
1.一种健身跑运动方案生成方法,其特征在于,包括以下步骤:
S1. 采集运动者健身跑过程中的速度和心率数据;
S2. 基于步骤S1采集的速度和心率数据,运用遗传算法对运动者的健身跑运动模型进行参数整定,建立运动者的健身跑运动模型;
S3. 基于步骤S2建立的健身跑运动模型,以及给定运动时间和安全有效心率区间参数,以在给定运动时间的健身跑过程中,心率落在安全有效心率区间内的时长最长为优化目标,建立运动者的健身跑最优运动方案。
2.根据权利要求1所述的一种健身跑运动方案生成方法,其特征在于,在步骤S2中,对于如式(1)的健身跑运动模型:
(1)
其中,x 1(t)、x 2(t)和u(t)分别表示运动者t时刻较0时刻的心率增量、t时刻较0时刻的疲劳程度增量和t时刻的速度,、分别表示x 1(t)、x 2(t)的微分,a 1、a 2、a 3、a 4和a 5为5个待估参数;按如下步骤建立运动者的健身跑运动模型:
S201. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群P(t),t=0,其中初始化参数a 1、a 2、a 3、a 4和a 5,个体编码采用实数编码,编码长度为5;
S202. 适应度计算:以心率的估计值与实测值的残差平方和作为适应度函数,根据式(1)和采集的速度实测值u(t),计算心率估计值x 1(t),将计算结果和采集的心率实测值输入适应度函数,计算各个体的适应度;
S203. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S204. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S203得到的下一代种群,返回步骤S202进行下一次迭代;
将求得的参数a 1、a 2、a 3、a 4和a 5的值代入式(1),得到运动者的健身跑运动模型。
3.根据权利要求1所述的一种健身跑运动方案生成方法,其特征在于,在步骤S3中,基于步骤S2建立的健身跑运动模型,将求解运动者的健身跑最优运动方案抽象为一个带约束的优化问题:在给定运动时间Ts和采样周样T下,在速度值u(t)中找出最优的速度值u*(t),使其满足健身跑过程中心率不超过预设的最大值V max,且心率达到安全有效心率区间[a, b]的次数N(t)最多;u(t)和N(t)的定义如式(2)和式(3):
(2)
(3)
其中,表示对T S 除于T的值向下取整,对应T S 时间内的采样次数,s i (t)表示在第i个采样周期内其速度值为1,而在其它采样周期内其速度值为0;c i 为待估参数,i=1,2,…,;按如下步骤估计个参数c i 的值:
S301. 种群初始化:设置种群规模、交叉概率、变异概率、最大进化代数,随机生成初始种群Q(t),t=0,初始化参数,个体编码采用实数编码,编码长度为;
S302. 适应度计算:以式(3)作为适应度函数,首先估计速度值u(t),然后将该值代入健身跑运动模型,得到心率值x 1(t),最后将计算结果输入适应度函数,计算各个体的适应度,对于x 1(t)的值超过V max的个体,将其适应度赋值为0;
S303. 采用整体算术杂交算子对种群进行交叉操作,然后采用自适应变异算子对种群进行变异操作,再采用带精英保留策略的轮盘赌方法选择下一代种群;
S304. 判断是否满足终止条件:最优个体的适应度达到给定的阈值,或者经过多次迭代后最优个体的适应度不再发生变化,或者迭代次数达到最大进化代数,是则输出最优个体的优化结果,计算结束,否则基于步骤S303得到的下一代种群,返回步骤S302进行下一次迭代;
将求得的参数的值代入式(2),得到运动者的健身跑最优运动方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410178014.2A CN103942445B (zh) | 2014-04-30 | 2014-04-30 | 一种健身跑运动方案生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410178014.2A CN103942445B (zh) | 2014-04-30 | 2014-04-30 | 一种健身跑运动方案生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103942445A true CN103942445A (zh) | 2014-07-23 |
CN103942445B CN103942445B (zh) | 2017-02-15 |
Family
ID=51190113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410178014.2A Active CN103942445B (zh) | 2014-04-30 | 2014-04-30 | 一种健身跑运动方案生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103942445B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467826A (zh) * | 2016-01-06 | 2016-04-06 | 上海应用技术学院 | 一种运动智能手表和最佳跑速获取方法 |
CN107456751A (zh) * | 2017-07-17 | 2017-12-12 | 南京宽乐健康产业有限公司 | 一种共享体适能健身系统智能控制方法 |
CN107705833A (zh) * | 2017-11-22 | 2018-02-16 | 中山大学 | 一种基于人体生理调节模型的运动推荐方法 |
CN108898795A (zh) * | 2018-06-12 | 2018-11-27 | 广州小鹏汽车科技有限公司 | 一种基于移动设备的行人状态实时分析方法及系统 |
CN109966122A (zh) * | 2017-12-28 | 2019-07-05 | 沈阳新松机器人自动化股份有限公司 | 一种智能判断疲劳度的行走辅助设备 |
CN110610233A (zh) * | 2019-09-19 | 2019-12-24 | 福建宜准信息科技有限公司 | 基于领域知识和数据驱动的健身跑心率预测方法 |
CN110634549A (zh) * | 2019-08-13 | 2019-12-31 | 刘芷含 | 针对超重人群的智能化运动管理装置及方法 |
CN111859624A (zh) * | 2020-06-24 | 2020-10-30 | 西安理工大学 | 一种基于nsga-ii算法的运动处方参数管理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279072B (zh) * | 2011-05-12 | 2013-07-10 | 天津大学 | 步行器测力系统中的蚁群遗传融合标定方法 |
-
2014
- 2014-04-30 CN CN201410178014.2A patent/CN103942445B/zh active Active
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467826A (zh) * | 2016-01-06 | 2016-04-06 | 上海应用技术学院 | 一种运动智能手表和最佳跑速获取方法 |
CN107456751A (zh) * | 2017-07-17 | 2017-12-12 | 南京宽乐健康产业有限公司 | 一种共享体适能健身系统智能控制方法 |
CN107705833A (zh) * | 2017-11-22 | 2018-02-16 | 中山大学 | 一种基于人体生理调节模型的运动推荐方法 |
CN109966122A (zh) * | 2017-12-28 | 2019-07-05 | 沈阳新松机器人自动化股份有限公司 | 一种智能判断疲劳度的行走辅助设备 |
CN108898795A (zh) * | 2018-06-12 | 2018-11-27 | 广州小鹏汽车科技有限公司 | 一种基于移动设备的行人状态实时分析方法及系统 |
CN110634549A (zh) * | 2019-08-13 | 2019-12-31 | 刘芷含 | 针对超重人群的智能化运动管理装置及方法 |
CN110610233A (zh) * | 2019-09-19 | 2019-12-24 | 福建宜准信息科技有限公司 | 基于领域知识和数据驱动的健身跑心率预测方法 |
CN111859624A (zh) * | 2020-06-24 | 2020-10-30 | 西安理工大学 | 一种基于nsga-ii算法的运动处方参数管理方法 |
CN111859624B (zh) * | 2020-06-24 | 2024-04-26 | 西安理工大学 | 一种基于nsga-ii算法的运动处方参数管理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103942445B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103942445A (zh) | 一种健身跑运动方案生成方法 | |
Holmgård et al. | Evolving personas for player decision modeling | |
WO2019231729A1 (en) | Physical activity training assistant | |
CN102626557A (zh) | 基于ga-bp算法的分子蒸馏工艺参数优化方法 | |
CN104899432A (zh) | 基于核函数组合的pso-lssvm脉动风速预测方法 | |
CN107563509A (zh) | 一种基于特征返回的条件dcgan模型的动态调整算法 | |
CN113663312B (zh) | 一种基于微惯性的无器械健身动作质量评估方法 | |
CN110610233B (zh) | 基于领域知识和数据驱动的健身跑心率预测方法 | |
CN106651001A (zh) | 一种改进神经网络的金针菇产量预测方法及实施系统 | |
CN109858798A (zh) | 关联改造措施与电压指标的电网投资决策建模方法及装置 | |
CN104217260B (zh) | 一种风场邻近多台风电机测量风速缺损值的组合填充系统 | |
Szabó et al. | Estimating task difficulty in educational games | |
Stankovic et al. | Relation and influences of sports climbers' specific strength on the success in sports climbing | |
Yin et al. | Effects of wave nonlinearity on submerged flexible vegetation dynamics and wave attenuation | |
Sannicandro et al. | Strength ability, endurance and anthropometric parameters in youth football: descriptive analysis and functional relationships | |
CN106178466A (zh) | 一种健身用度量分析方法及系统 | |
CN105044531B (zh) | 一种基于ekf和fsa的动态信号参数辨识方法 | |
CN106408087A (zh) | 面向水下弱信号检测的强化学习自适应随机共振方法 | |
CN103123668A (zh) | 一种基于混合单元法的空间绳系机器人系统的仿真方法 | |
CN103781108A (zh) | 一种基于神经网络的无线传感器网络数据预测方法 | |
CN103473812A (zh) | 一种基于tps的等深线追踪算法 | |
CN104898689B (zh) | 一种区间信息条件下的多无人机目标分配决策方法 | |
CN101884832A (zh) | 虚拟角色健身互动系统及其方法 | |
CN109117462A (zh) | 一种游戏活跃玩家人数的预测方法 | |
CN105973266A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20181031 Address after: 350007 Jinshan Industrial Park, Fuzhou, Cangshan, Fujian, 15 Patentee after: Fujian Ezon Information Technology Co., Ltd. Address before: 350007 No. 8, Shang three road, Cangshan District, Fuzhou, Fujian. Patentee before: Fujian Normal University |