CN110755078A - 基于置信度加权的膝关节运动疲劳评估预测系统及方法 - Google Patents

基于置信度加权的膝关节运动疲劳评估预测系统及方法 Download PDF

Info

Publication number
CN110755078A
CN110755078A CN201911040603.3A CN201911040603A CN110755078A CN 110755078 A CN110755078 A CN 110755078A CN 201911040603 A CN201911040603 A CN 201911040603A CN 110755078 A CN110755078 A CN 110755078A
Authority
CN
China
Prior art keywords
knee joint
confidence
bending
fatigue
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
Application number
CN201911040603.3A
Other languages
English (en)
Other versions
CN110755078B (zh
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.)
Fujian Agriculture and Forestry University
Original Assignee
Fujian Agriculture and Forestry University
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 Fujian Agriculture and Forestry University filed Critical Fujian Agriculture and Forestry University
Priority to CN201911040603.3A priority Critical patent/CN110755078B/zh
Publication of CN110755078A publication Critical patent/CN110755078A/zh
Application granted granted Critical
Publication of CN110755078B publication Critical patent/CN110755078B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/1036Measuring load distribution, e.g. podologic studies
    • A61B5/1038Measuring plantar pressure during gait
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1071Measuring physical dimensions, e.g. size of the entire body or parts thereof measuring angles, e.g. using goniometers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00Evaluating a particular growth phase or type of persons or animals
    • A61B2503/10Athletes

Abstract

本发明涉及一种基于置信度加权的膝关节运动疲劳评估预测方法,包括学习阶段和预测阶段,学习阶段包括:1)分别通过弯曲传感器、陀螺仪测量膝关节的弯曲度;通过足底压力传感器测得足底的压力负载;2)计算每个时刻的综合性负载累积值;3)以弯曲传感器和陀螺仪测得值的一致性定义置信度,计算并更新置信度值;4)当膝关节的最大弯曲度每减少一整数度时,计算当时的综合性负载累积值和标准化置信度值,并估算关节疲劳度;5)进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的回归关系;进入预测阶段,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。该系统及方法有利于对膝关节的运动疲劳程度进行准确的评估及预测。

Description

基于置信度加权的膝关节运动疲劳评估预测系统及方法
技术领域
本发明涉及涉及智能穿戴设备技术领域,具体涉及一种基于置信度加权的膝关节运动疲劳评估预测系统及方法。
背景技术
由于运动锻炼时常需要进行屈膝运动,例如跑步、羽毛球、足球等,而过量的重负载屈膝运动超过一定限度后将会损伤膝关节,因此发明一直膝关节运动疲劳评估设备及预测方法就显得十分必要。如果能在运动即将过量之时进行报警,敦促运动员或是运动爱好者及时降低运动负荷或停止运动,将达到保护膝关节的目的。
发明内容
本发明的目的在于提供一种基于置信度加权的膝关节运动疲劳评估预测系统及方法,该系统及方法有利于对膝关节的运动疲劳程度进行准确的评估及预测。
为实现上述目的,本发明采用的技术方案是:一种基于置信度加权的膝关节运动疲劳评估预测系统,包括:
弯曲传感器,设置于膝关节处,用于通过弯曲传感器的电阻变化测得膝关节的弯曲度;
第一陀螺仪和第二陀螺仪,分别设置于大腿和小腿处,用于通过两个陀螺仪之间的角度差测得膝关节的弯曲度;
足底压力传感器,设置于足底,用于测得足底的压力负载;以及评估预测系统,通过无线网络分别与弯曲传感器、第一陀螺仪、第二陀螺仪及足底压力传感器连接以接收各传感器测得的数据,用于进行置信度加权的线性回归分析,以得到综合性负载与膝关节疲劳度之间的回归关系,还用于在得到了综合性负载与膝关节疲劳度之间的回归关系后,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
进一步地,所述弯曲传感器包括分别设置于膝关节内外侧的第一弯曲传感器和第二弯曲传感器,以通过取平均值减少测得的膝关节弯曲度的误差。
本发明还提供了一种基于置信度加权的膝关节运动疲劳评估预测方法,包括学习阶段和预测阶段,所述学习阶段包括以下步骤:
步骤1)在膝关节处设置弯曲传感器,以通过弯曲传感器的电阻变化测得膝关节的弯曲度;在大腿和小腿处分别设置第一陀螺仪和第二陀螺仪,以通过两个陀螺仪之间的角度差测得膝关节的弯曲度;在足底处设置足底压力传感器,以测得足底的压力负载;
步骤2)计算每个时刻的综合性负载累积值at,即从运动开始时刻直到t时刻,由于人体运动带给膝关节的压力效应的累积值;
步骤3)以弯曲传感器和陀螺仪测得的弯曲度的一致性定义置信度,计算每个时刻的置信度值Ct new,并对临近时刻的置信度值进行更新;
步骤4)通过弯曲传感器和陀螺仪测量得到膝关节的最大弯曲度角度值,当膝关节的最大弯曲度每减少一整数度时,计算当时对应的综合性负载累积值ai和标准化置信度值Ci norm,并通过减少的角度估算关节疲劳度fi,直至得到满足数量要求的n个点对(ai,fi),i=1,2,…,n,否则返回步骤2)进行循环计算;
步骤5)利用得到的点对(ai,fi)进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的回归关系;
使用者在学习阶段得到综合性负载与其膝关节疲劳度的个性化回归关系后,进入所述预测阶段,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
进一步地,在膝关节的内外侧分别设置第一弯曲传感器和第二弯曲传感器,通过取平均值减少测得的膝关节弯曲度的误差。
进一步地,步骤2)中,t时刻的综合性负载累积值at的计算方法为:
通过式(1)计算t时刻的综合性负载lt
Figure BDA0002252716500000021
其中,lt表示t时刻的综合性负载;P表示足底压力传感器测得的足底压力;Pthr表示足底压力阈值;Ca是陀螺仪和弯曲传感器测得的弯曲度弧度值的平均值;Cthr是在平时行走情况下采样的弯曲角弧度阈值;β和γ是设定的调整因子;k是用于调整曲线形状的指数函数底数;BMI表示体重指数;δ是调整参数,是正常体重中BMI的最高值;Φ是设定的补充调整参数;b是另一个用于调整曲线形状的指数函数底数;Δtg是上一次膝关节弯曲与当前膝关节弯曲之间的间隔时间,用于表征膝关节运动的频率;Δtthr是在平时行走情况下采样的膝关节弯曲时间间隔阈值;ρ是设定的调整因子;h也是用于调整曲线形状的指数函数底数;ε、θ、ω是设定的调整因子;
得到t时刻的综合性负载lt后,通过式(2)计算0-t时刻的综合性负载累积值at
Figure BDA0002252716500000022
进一步地,步骤3)中,t时刻的置信度值Ct new的计算方法为:
通过式(3)计算t时刻弯曲传感器和陀螺仪测得的弯曲角弧度值的差Δe:
Δe=|Cg-Ct| (3)
其中,Cg、Ct分别表示陀螺仪和弯曲传感器测得的弯曲角弧度值;
获得了Δe后,通过式(4)和式(5)对t时刻及其临近时刻的置信度值进行计算:
Figure BDA0002252716500000031
Figure BDA0002252716500000032
Ct new和Ct pre分别为更新后和更新前的置信度值,置信度初值设为1;g(x)是置信度传播机制的调整函数,其返回值反映传感器误差的影响,是用于调整曲线g(x)形状的系数,增大
Figure BDA00022527165000000310
值将增加置信度值对传感器误差的灵敏度;σ是用于调整速度的加权系数,η是设定的调节常数;
进一步可通过式(7)计算t时刻的标准化置信度值Ci norm
Figure BDA0002252716500000034
进一步地,步骤5)中,由式(8)和式(9)计算置信度加权的线性回归分析:
Figure BDA0002252716500000035
Figure BDA0002252716500000036
其中,
Figure BDA0002252716500000037
表示回归系数;A表示由每个点处的综合负荷累积值ai组成的向量;F表示由膝关节疲劳度fi组成的向量;C表示权重对角矩阵,其以标准置信度值作为其元素。
进一步地,步骤6)中,通过式(10)预测膝关节的疲劳度:
其中,fend表示设定的疲劳度,
Figure BDA0002252716500000039
是回归系数,anow是当前综合负载累积值,lnow表示当前的综合负载;预测时间t用于警告用户在当前锻炼负载下还剩多少时间会导致膝盖极限疲劳,从而使用户降低负载强度或减少运动时间,以保障膝关节健康。
相较于现有技术,本发明具有以下有益效果:提供了一种基于置信度加权的膝关节运动疲劳评估预测系统及方法,该系统及方法设置弯曲传感器和陀螺仪分别测量膝关节的弯曲度,并通过两种传感器测得值的一致性来定义置信度,然后进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的回归关系,在此基础上,即可通过监测综合性负载累积值,对膝关节的疲劳度进行预测,从而实现了对膝关节运动疲劳度的准确评估及预测,能在运动即将过量之时进行报警,达到保护膝关节的目的,因此,具有很强的实用性和广阔的应用前景。
附图说明
图1是本发明实施例的方法中学习阶段的实现流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步说明。
本发明提供了一种基于置信度加权的膝关节运动疲劳评估预测系统,包括:
第一弯曲传感器和第二弯曲传感器,分别设置于膝关节内外侧,用于通过弯曲传感器的电阻变化测得膝关节的弯曲度;设置两个弯曲传感器,是为了取平均值减少测得的膝关节弯曲度的误差;
第一陀螺仪和第二陀螺仪,分别设置于大腿和小腿处,用于通过两个陀螺仪之间的角度差测得膝关节的弯曲度;
足底压力传感器,设置于足底,用于测得足底的压力负载;以及
评估预测系统,通过无线网络分别与弯曲传感器、第一陀螺仪、第二陀螺仪及足底压力传感器连接以接收各传感器测得的数据,用于进行置信度加权的线性回归分析,以得到综合性负载与膝关节疲劳度之间的回归关系,还用于在得到了综合性负载与膝关节疲劳度之间的回归关系后,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
本发明还提供了基于上述系统的一种基于置信度加权的膝关节运动疲劳评估预测方法,包括学习阶段和预测阶段,如图1所示,所述学习阶段包括以下步骤:
步骤1)在膝关节的内外侧处分别设置第一弯曲传感器和第二弯曲传感器,膝关节的弯曲会改变弯曲传感器的电阻,以通过弯曲传感器的电阻变化测得膝关节的弯曲度,通过取平均值减少测得的膝关节弯曲度的误差;在大腿和小腿处分别设置第一陀螺仪和第二陀螺仪,以通过两个陀螺仪之间的角度差测得膝关节的弯曲度;在足底处设置足底压力传感器,以测得足底的压力负载。
步骤2)计算每个时刻的综合性负载累积值at,即从运动开始时刻直到t时刻,由于人体运动带给膝关节的压力效应的累积值,其与膝关节的疲劳程度直接相关,可用于判断及预测膝关节的疲劳程度。其中,t时刻的综合性负载累积值at的计算方法为:
通过式(1)计算t时刻的综合性负载lt
其中,lt表示t时刻的综合性负载;P表示足底压力传感器测得的足底压力;Pthr表示足底压力阈值,其值是在用户每日步行期间采样得到的常规压力值;Ca是陀螺仪和弯曲传感器测得的弯曲度弧度值的平均值;Cthr是在平时行走情况下采样的弯曲角弧度阈值;β和γ是设定的调整因子;k是用于调整曲线形状的指数函数底数;BMI表示体重指数,是目前世界上最常用的测量体重与身高之比的工具,它使用身高和体重之间的比率来衡量一个人是太瘦还是太胖;δ是调整参数,是正常体重中BMI的最高值,世界上通常为25,中国地区为24;Φ是设定的补充调整参数;b是另一个用于调整曲线形状的指数函数底数;Δtg是上一次膝关节弯曲与当前膝关节弯曲之间的间隔时间,用于表征膝关节运动的频率;Δtthr是在平时行走情况下采样的膝关节弯曲时间间隔阈值;ρ是设定的调整因子;h也是用于调整曲线形状的指数函数底数;ε、θ、ω是设定的调整因子;
得到t时刻的综合性负载lt后,通过式(2)计算0-t时刻的综合性负载累积值at
Figure BDA0002252716500000052
步骤3)以弯曲传感器和陀螺仪测得的弯曲度的一致性定义置信度,计算每个时刻的置信度值Ct new,并对临近时刻的置信度值进行更新。其中,t时刻的置信度值Ct new的计算方法为:
通过式(3)计算t时刻弯曲传感器和陀螺仪测得的弯曲角弧度值的差Δe:
Δe=|Cg-Ct| (3)
其中,Cg、Ct分别表示陀螺仪和弯曲传感器测得的弯曲角弧度值;
获得了Δe后,通过式(4)和式(5)对t时刻及其临近时刻的置信度值进行计算:
Figure BDA0002252716500000053
Figure BDA0002252716500000054
Ct new和Ct pre分别为更新后和更新前的置信度值,置信度初值设为1,有误差就被减,误差越大,周围点误差越大,就减得越多,最后剩的值越小,说明越不可靠;g(x)是置信度传播机制的调整函数,其返回值反映传感器误差的影响,是用于调整曲线g(x)形状的系数,增大
Figure BDA0002252716500000061
值将增加置信度值对传感器误差的灵敏度;σ是用于调整速度的加权系数,η是设定的调节常数;
进一步可通过式(7)计算t时刻的标准化置信度值Ci norm
Figure BDA0002252716500000062
步骤4)通过弯曲传感器和陀螺仪测量得到膝关节的最大弯曲度角度值,这里取两种传感器测量所得值的平均值,当膝关节的最大弯曲度每减少一整数度时,计算当时对应的综合性负载累积值ai和标准化置信度值Ci norm,并通过减少的角度估算关节疲劳度fi,直至得到满足数量要求的n个点对(ai,fi),i=1,2,…,n,否则返回步骤2)进行循环计算。
步骤5)利用得到的点对(ai,fi)进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的回归关系。其中,由式(8)和式(9)计算置信度加权的线性回归分析:
Figure BDA0002252716500000063
Figure BDA0002252716500000064
其中,表示回归系数;A表示由每个点处的综合负荷累积值ai组成的向量;F表示由膝关节疲劳度fi组成的向量;C表示权重对角矩阵,其以标准置信度值作为其元素。
使用者在学习阶段得到综合性负载与其膝关节疲劳度的个性化回归关系后,进入所述预测阶段,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
其中,通过式(10)预测膝关节的疲劳度:
Figure BDA0002252716500000066
其中,fend表示设定的疲劳度,是回归系数,anow是当前综合负载累积值,lnow表示当前的综合负载;预测时间t用于警告用户在当前锻炼负载下还剩多少时间会导致膝盖极限疲劳,从而使用户降低负载强度或减少运动时间,以保障膝关节健康。
下面进一步说明本发明的发明思路和技术方案。
本发明的膝关节疲劳度评估及预测基于综合性负载l来实现,综合性负载l反映了人体自身重量、膝关节运动频率及足底压力对膝关节疲劳度进展的影响。时刻t时的综合性负载lt的计算如式(1)所示:
Figure BDA0002252716500000071
其中lt表示时刻t时的综合负荷;P表示足底压力传感器获得的足底压力;Pthr表示足底压力设置阈值,其值是在用户每日步行期间采样得到的常规压力值;Ca是陀螺仪和弯曲传感器获得的弯曲度弧度值的平均值;Cthr是在平时行走情况下采样的弯曲角弧度阈值;β和γ是由用户设定的调整因子;k是用于调整曲线形状的指数函数底数;BMI代表体重指数,它是目前世界上最常用的测量体重与身高之比的工具,它使用身高和体重之间的比率来衡量一个人是太瘦还是太胖;δ是一个调整参数,是正常体重中BMI的最高值,世界上通常为25,中国地区为24;Φ是用户设定的补充调整参数;b是另一个用于调整曲线形状的指数函数底数;Δtg是上一次膝关节弯曲与当前膝关节弯曲之间的间隔时间,用于表征膝关节运动的频率;Δtthr是在平时行走情况下采样的膝关节弯曲时间间隔阈值;ρ是用户设定的调整因子;h也是用于调整曲线形状的指数函数底数;ε,θ,ω是由用户设定的调整因子。
Figure BDA0002252716500000072
在定义了t时刻的综合性负载lt后,通过式(2)可以得到0-t时刻的综合性负载累积值at。综合性负载累积值at表示了从运动开始时刻直到时刻t,由于人体运动带给膝关节的压力效应的累积值,其与膝关节的疲劳程度直接相关,可用于判断及预测膝关节的疲劳程度。
由于本发明硬件设备配置了陀螺仪及弯曲传感器两种可以获得膝关节弯曲角度的传感器,所以该两种传感器的值可以用于相互校验。由于运动中的抖动及其他效应,单一一种传感器的实时值可能有一定误差,导致综合性负载累计值at产生偏差,从而影响关节疲劳度的判断和估计。因此本发明提出用两种传感器检测值的一致性来定义置信度,并使置信度高的时刻在进行线性规划寻找规律时的参与度高于置信度低的时刻,从而削弱噪声及测量误差对最终结果的影响。式(3)可计算出t时刻两种不同传感器测量得到弧度值的差Δe,其中Cg及Ct分别表示由陀螺仪和弯曲传感器测得的弯曲角弧度值。
Δe=|Cg-Ct| (3)
Figure BDA0002252716500000073
Figure BDA0002252716500000081
获得了Δe后,本发明通过式(4)与式(5)对t时刻及其临近时刻的置信度值进行更新。Ct new及Ct pre分别为更新前和更新后的置信度值,由于抖动在时域上具有连续性,且剧烈的程度随距离抖动中心的远离而逐步衰减,因此式(5)采用1/2的次方数来体现逐渐削减的效果。g(x)是置信度传播机制的调整函数,其返回值决定了传感器误差的影响,而
Figure BDA0002252716500000082
用于调整曲线g(x)的形状。较大的
Figure BDA0002252716500000083
值会增加置信度值对传感器误差的灵敏度。由于强烈运动引起的传感误差是连续的,如果时间t的传感误差超过阈值ethr,不仅t时刻的置信度值需减小,而且相邻时间片的置信度值也应减小以反映这种影响。此外,随着离t时刻距离的增加,减小的幅度逐渐减小。σ是用于调整速度的加权系数,η是由设计者决定的调节常数。式(6)给出了t-2时刻的置信度更新公式示例。由于距时间片t的距离为2,因此削减幅度为1/22
Figure BDA0002252716500000084
Figure BDA0002252716500000085
在置信度值更新后,通过式(7)可取得t时刻的标准化置信度值Ct norm,其值小于等于1且囊括了时刻0-t的全部置信度情况。
在获得了综合性负载累积值at后,可以通过分析其与膝关节疲劳度之间的关系,来判断当前的膝关节疲劳度值ft以及预测膝关节将来还能承受多长时间此等负载。为了获得综合性负载累积值at与膝关节疲劳度ft直接的关系,本发明将其假设为大致线性关系,于是需要通过线性回归分析来获得两者间的关系。于是该关系分析问题就变成了一种有监督的模式识别问题,可通过线性规划来解决,由于上述提到的传感器感应误差的问题,可将上述设置的标准化置信度值Ct norm纳入分析,从而形成置信度加权的线性回归分析。由本领域相关研究得知,膝关节疲劳会使其最大弯曲度降低,因此通过检测最大弯曲度可以初步判断当前的膝关节疲劳度ft,但却无法做出预测。从膝关节初始阶段到最疲劳阶段,最大弯曲角度减少约有10度,故以每减少1度作为节点,分别采样当时的综合性负载累积值ai及标准置信度累计值Ci norm,并通过角度减少估测疲劳度fi(其中i为序号),就可以形成点对(ai,fi)以进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的关系,从而实现对未来的预测。整个分析学习流程如图1所示,置信度加权的线性回归分析由式(8)及式(9)计算所得。其中代表回归系数;A表示由每个关键点处的综合负荷累积值ai组成的向量;F是由膝关节疲劳度fi组成的向量;C是如式(9)所示的权重对角矩阵,其以标准置信度值作为其元素。
Figure BDA0002252716500000091
Figure BDA0002252716500000092
在获得了回归系数后,就可以在日常的运动中,通过监测综合性负载累积值at来对膝关节的疲劳程度进行预测,从而提醒用户及时休息。预测公式如式(10)所示,其中fend表示用户设定的疲劳度,是回归系数,anow是当前综合负载累积值,lnow表示当前的综合负载。预测时间t可用于警告用户在当前锻炼负载下还剩多少时间会导致膝盖极限疲劳,从而使用户降低负载强度或是减少运动时间以保障膝关节健康。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。

Claims (8)

1.一种基于置信度加权的膝关节运动疲劳评估预测系统,其特征在于,包括:
弯曲传感器,设置于膝关节处,用于通过弯曲传感器的电阻变化测得膝关节的弯曲度;
第一陀螺仪和第二陀螺仪,分别设置于大腿和小腿处,用于通过两个陀螺仪之间的角度差测得膝关节的弯曲度;
足底压力传感器,设置于足底,用于测得足底的压力负载;以及
评估预测系统,通过无线网络分别与弯曲传感器、第一陀螺仪、第二陀螺仪及足底压力传感器连接以接收各传感器测得的数据,用于进行置信度加权的线性回归分析,以得到综合性负载与膝关节疲劳度之间的回归关系,还用于在得到了综合性负载与膝关节疲劳度之间的回归关系后,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
2.根据权利要求1所述的基于置信度加权的膝关节运动疲劳评估预测系统,其特征在于,所述弯曲传感器包括分别设置于膝关节内外侧的第一弯曲传感器和第二弯曲传感器,以通过取平均值减少测得的膝关节弯曲度的误差。
3.一种基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,包括学习阶段和预测阶段,所述学习阶段包括以下步骤:
步骤1)在膝关节处设置弯曲传感器,以通过弯曲传感器的电阻变化测得膝关节的弯曲度;在大腿和小腿处分别设置第一陀螺仪和第二陀螺仪,以通过两个陀螺仪之间的角度差测得膝关节的弯曲度;在足底处设置足底压力传感器,以测得足底的压力负载;
步骤2)计算每个时刻的综合性负载累积值at,即从运动开始时刻直到t时刻,由于人体运动带给膝关节的压力效应的累积值;
步骤3)以弯曲传感器和陀螺仪测得的弯曲度的一致性定义置信度,计算每个时刻的置信度值Ct new,并对临近时刻的置信度值进行更新;
步骤4)通过弯曲传感器和陀螺仪测量得到膝关节的最大弯曲度角度值,当膝关节的最大弯曲度每减少一整数度时,计算当时对应的综合性负载累积值ai和标准化置信度值Ci norm,并通过减少的角度估算关节疲劳度fi,直至得到满足数量要求的n个点对(ai,fi),i=1,2,…,n,否则返回步骤2)进行循环计算;
步骤5)利用得到的点对(ai,fi)进行置信度加权的线性回归分析,得到综合性负载与膝关节疲劳度之间的回归关系;
使用者在学习阶段得到综合性负载与其膝关节疲劳度的个性化回归关系后,进入所述预测阶段,通过监测综合性负载累积值,对膝关节的疲劳度进行预测。
4.根据权利要求3所述的基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,在膝关节的内外侧分别设置第一弯曲传感器和第二弯曲传感器,通过取平均值减少测得的膝关节弯曲度的误差。
5.根据权利要求3所述的基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,步骤2)中,t时刻的综合性负载累积值at的计算方法为:
通过式(1)计算t时刻的综合性负载lt
其中,lt表示t时刻的综合性负载;P表示足底压力传感器测得的足底压力;Pthr表示足底压力阈值;Ca是陀螺仪和弯曲传感器测得的弯曲度弧度值的平均值;Cthr是在平时行走情况下采样的弯曲角弧度阈值;β和γ是设定的调整因子;k是用于调整曲线形状的指数函数底数;BMI表示体重指数;δ是调整参数,是正常体重中BMI的最高值;Φ是设定的补充调整参数;b是另一个用于调整曲线形状的指数函数底数;Δtg是上一次膝关节弯曲与当前膝关节弯曲之间的间隔时间,用于表征膝关节运动的频率;Δtthr是在平时行走情况下采样的膝关节弯曲时间间隔阈值;ρ是设定的调整因子;h也是用于调整曲线形状的指数函数底数;ε、θ、ω是设定的调整因子;
得到t时刻的综合性负载lt后,通过式(2)计算0-t时刻的综合性负载累积值at
Figure FDA0002252716490000022
6.根据权利要求5所述的基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,步骤3)中,t时刻的置信度值Ct new的计算方法为:
通过式(3)计算t时刻弯曲传感器和陀螺仪测得的弯曲角弧度值的差Δe:
Δe=|Cg-Ct| (3)
其中,Cg、Ct分别表示陀螺仪和弯曲传感器测得的弯曲角弧度值;
获得了Δe后,通过式(4)和式(5)对t时刻及其临近时刻的置信度值进行计算:
Figure FDA0002252716490000023
Figure FDA0002252716490000024
Ct new和Ct pre分别为更新后和更新前的置信度值,置信度初值设为1;g(x)是置信度传播机制的调整函数,其返回值反映传感器误差的影响,
Figure FDA0002252716490000025
是用于调整曲线g(x)形状的系数,增大值将增加置信度值对传感器误差的灵敏度;σ是用于调整速度的加权系数,η是设定的调节常数;
进一步可通过式(7)计算t时刻的标准化置信度值Ci norm
Figure FDA0002252716490000032
7.根据权利要求6所述的基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,步骤5)中,由式(8)和式(9)计算置信度加权的线性回归分析:
Figure FDA0002252716490000033
Figure FDA0002252716490000034
其中,
Figure FDA0002252716490000035
表示回归系数;A表示由每个点处的综合负荷累积值ai组成的向量;F表示由膝关节疲劳度fi组成的向量;C表示权重对角矩阵,其以标准置信度值作为其元素。
8.根据权利要求7所述的基于置信度加权的膝关节运动疲劳评估预测方法,其特征在于,步骤6)中,通过式(10)预测膝关节的疲劳度:
Figure FDA0002252716490000036
其中,fend表示设定的疲劳度,
Figure FDA0002252716490000037
是回归系数,anow是当前综合负载累积值,lnow表示当前的综合负载;预测时间t用于警告用户在当前锻炼负载下还剩多少时间会导致膝盖极限疲劳,从而使用户降低负载强度或减少运动时间,以保障膝关节健康。
CN201911040603.3A 2019-10-29 2019-10-29 基于置信度加权的膝关节运动疲劳评估预测系统及方法 Active CN110755078B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911040603.3A CN110755078B (zh) 2019-10-29 2019-10-29 基于置信度加权的膝关节运动疲劳评估预测系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911040603.3A CN110755078B (zh) 2019-10-29 2019-10-29 基于置信度加权的膝关节运动疲劳评估预测系统及方法

Publications (2)

Publication Number Publication Date
CN110755078A true CN110755078A (zh) 2020-02-07
CN110755078B CN110755078B (zh) 2022-08-09

Family

ID=69334516

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911040603.3A Active CN110755078B (zh) 2019-10-29 2019-10-29 基于置信度加权的膝关节运动疲劳评估预测系统及方法

Country Status (1)

Country Link
CN (1) CN110755078B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5197488A (en) * 1991-04-05 1993-03-30 N. K. Biotechnical Engineering Co. Knee joint load measuring instrument and joint prosthesis
CN102438519A (zh) * 2009-08-28 2012-05-02 艾伦·约瑟夫·泽尔纳 通过运动分析来表征身体能力
US20140364962A1 (en) * 2013-05-23 2014-12-11 Rehabilitation Institute Of Chicago Effective shape controller for lower limb
US20160310342A1 (en) * 2013-12-27 2016-10-27 Korea Polytechnic University Industry Academic Cooperation Foundation System and method for controlling joint angle of knee-joint type walking training robot
CN108665956A (zh) * 2018-06-05 2018-10-16 陈燕 体能和肌肉能力的综合评估系统及方法
CN108670264A (zh) * 2018-05-25 2018-10-19 福建农林大学 基于关节弯曲与足底压力的膝关节运动疲劳度监测系统
CN109793500A (zh) * 2019-01-24 2019-05-24 河南省人民医院 膝关节负荷力学分析装置
CN109953761A (zh) * 2017-12-22 2019-07-02 浙江大学 一种下肢康复机器人感知系统及运动意图推理方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5197488A (en) * 1991-04-05 1993-03-30 N. K. Biotechnical Engineering Co. Knee joint load measuring instrument and joint prosthesis
CN102438519A (zh) * 2009-08-28 2012-05-02 艾伦·约瑟夫·泽尔纳 通过运动分析来表征身体能力
US20140364962A1 (en) * 2013-05-23 2014-12-11 Rehabilitation Institute Of Chicago Effective shape controller for lower limb
US20160310342A1 (en) * 2013-12-27 2016-10-27 Korea Polytechnic University Industry Academic Cooperation Foundation System and method for controlling joint angle of knee-joint type walking training robot
CN109953761A (zh) * 2017-12-22 2019-07-02 浙江大学 一种下肢康复机器人感知系统及运动意图推理方法
CN108670264A (zh) * 2018-05-25 2018-10-19 福建农林大学 基于关节弯曲与足底压力的膝关节运动疲劳度监测系统
CN108665956A (zh) * 2018-06-05 2018-10-16 陈燕 体能和肌肉能力的综合评估系统及方法
CN109793500A (zh) * 2019-01-24 2019-05-24 河南省人民医院 膝关节负荷力学分析装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TAIAN VIEIRA等: "Design and Test of a Biomechanical Model for the Estimation of Knee Joint Angle During Indoor Rowing: Implications for FES-Rowing Protocols in Paraplegia", 《IEEE TRANSACTIONS ON NEURAL SYSTEMS AND REHABILITATION ENGINEERING》 *

Also Published As

Publication number Publication date
CN110755078B (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
CN105210067B (zh) 计算用户的与体育锻炼有关的生理状态
CN102292025B (zh) 确定用户的能量消耗
TWI393579B (zh) The state of the muscle movement state analysis system, methods and computer program products
FI125723B (en) Portable activity tracking device and associated method
US11779226B2 (en) Method and a system for determining the maximum heart rate of a user in a freely performed physical exercise
Li et al. An adaptive and on-line IMU-based locomotion activity classification method using a triplet Markov model
KR102107379B1 (ko) 3축 동작측정기를 이용한 노쇠 예측 방법, 3축 동작측정기를 이용한 노쇠 예측 시스템 및 웨어러블 노쇠 예측 장치
CN106667493A (zh) 人体平衡性评估系统及评估方法
US10892055B2 (en) Motor function estimation information generation apparatus, motor function estimation system, motor function estimation information generation method, motor function estimation method and storage medium
Pande et al. Energy expenditure estimation with smartphone body sensors
KR20140124590A (ko) 비정상 보행 여부를 측정하는 방법
EP3043875A1 (en) Method to determine body's physiological response to physical exercise for assessing readiness and to provide feedback, and system for implementing the method
CN108836344B (zh) 步长步频估算方法和装置及步态检测仪
Liu et al. A wearable gait analysis and recognition method for Parkinson’s disease based on error state Kalman filter
CN110755078B (zh) 基于置信度加权的膝关节运动疲劳评估预测系统及方法
CN113331829B (zh) 一种足底信息监测方法及智能鞋垫装置
US10542933B2 (en) Exercise test evaluation system, exercise test evaluation apparatus, exercise test evaluation method, and non-transitory computer readable recording medium
JP6586173B2 (ja) 歩行状態検出装置及び歩行状態検出プログラム
US20220110587A1 (en) Body motion analysis data treatment
US20220280090A1 (en) Body condition measurement tape and method for providing body condition monitoring service by using same
KR102268445B1 (ko) 관성 정보 기반 보행 안정성 평가 장치 및 보행 안정성 평가 방법
US10670450B2 (en) Apparatus and method for weight measurement during motion using force waveforms
Mazumder et al. Development of an adaptive gait characterizer
JP6738249B2 (ja) 歩行分析方法及び歩行分析装置
KR101441815B1 (ko) 운동량 측정 방법 및 장치 및 이를 이용하는 운동기구

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant