一种用于可穿戴式设备的人体髋关节位移测量方法与流程

1.本发明属于可穿戴传感器领域，具体涉及一种用于可穿戴式设备的人体髋关节位移测量方法。


背景技术：

2.髋关节运动是反映人类运动能力的重要指标之一，拥有较大的应用价值，因此髋关节的位移测量被广泛研究。大型实验室测量设备如光学式运动捕捉系统、测力台等等常常用于临床上的髋关节位移测量，但其存在价格昂贵，受空间限制，使用不便等缺点。可穿戴传感器作为包括惯性测量单元，超声波传感器，微型摄像头等等，相比于大型实验室测量步态参数的设备如光学式运动捕捉系统、测力台等等，以其小巧、廉价，不受时间、空间限制、易于推广等优点作为新兴技术被广泛应用于步态测量领域。然而，目前仅有较少研究使用放置在人体下肢部位的可穿戴传感器测量髋关节位移。


技术实现要素：

3.本发明的目的在于提供一种用于可穿戴式设备的人体髋关节位移测量方法。
4.本发明中所涉及的部分名词含义如下：
5.摆动期及支撑期为人走路步态周期的两个组成部分，人的走路步态周期如图1所示，以右腿为例，右脚落地是该步态周期的开始，右腿支撑人体，左腿向前迈出，左脚落地，随后左腿开始支撑人体，右脚离地，这段时间为右腿的支撑期；右脚离地后，右腿向前迈出并落地，为右腿的摆动期；右脚落地后完成一个步态周期。
6.平足期是指支撑期内脚面与地面贴合，踝关节几乎静止的一段时间。
7.跨步长是指人行走过程中同一只脚相邻着地点之间的距离。
8.梯度下降法是常用的优化方法之一，属于迭代法的一种，可以用于求解最小二乘问题。
9.为方便描述，定义矢状面，如图2所示，矢状面是垂直于水平面并将人体分为左右两部分的面，人在行走时，主要运动都发生在矢状面内。
10.为方便描述，定义x轴、y轴，如图3所示：x轴为矢状面内水平向前方向，y轴为矢状面内竖直向上方向。
11.为方便描述，定义小腿姿态角度为矢状面内小腿与y轴的夹角，定义大腿姿态角度为矢状面内大腿与y轴的夹角，定义对侧腿姿态角为身体对侧的腿(将身体对侧的小腿和大腿视为一个整体)与y轴的夹角，如图5所示。
12.本发明为解决技术问题，所采用的具体技术方案如下：
13.一种用于可穿戴式设备的人体髋关节位移测量方法，其包括以下步骤：
14.s1：利用佩戴于两小腿脚踝外侧上方的两个可穿戴式设备识别用户行走过程中的支撑期、摆动期、平足期和步态周期，并测量当前用户的小腿姿态角度、跨步长、踝关节x轴位移和踝关节y轴位移；
15.s2：通过所述小腿姿态角度、踝关节x轴位移、踝关节y轴位移计算支撑期内各时刻髋关节x轴位移和髋关节y轴位移；
16.s3：基于所述小腿姿态角度、跨步长、踝关节x轴位移和踝关节y轴位移，首先用待定的多项式来表征人体摆动期过程中大腿姿态角度随时间的变化曲线，利用处于摆动期内的当前腿对应的大腿姿态角度、小腿姿态角度、大腿长度和小腿长度，计算出摆动期内髋关节的x轴和y轴位移；利用人体在所述当前腿处于摆动期时身体对侧腿处于支撑期且身体对侧膝关节无弯曲的假设，计算所述当前腿的对侧腿姿态角度；再利用对侧腿姿态角度和对侧腿的长度，计算对侧腿的髋关节运动x轴和y轴位移，将两腿的髋关节位移做差，利用人体运动能量总是取最小原理，迭代优化所述多项式的待定系数，直至两腿髋关节位移差最小，利用优化迭代后的大腿姿态角度随时间的变化曲线，求得当前用户摆动期内的髋关节x轴位移和髋关节y轴位移。
17.作为优选，所述s2中，支撑期内各时刻髋关节x轴位移h
x
和髋关节y轴位移h
y
：
18.h
x
＝a
x
+l
·
sinθ
s
19.h
y
＝a
y
+l
·
cosθ
s
20.式中：a
x
、a
y
分别为踝关节x轴位移、踝关节y轴位移，l为用户的大腿长度和小腿长度之和，θ
s
为小腿姿态角度。
21.作为优选，所述s3中，摆动期内的髋关节x轴位移和髋关节y轴位移的计算过程包括以下步骤：
22.s31：基于三次多项式定义人体摆动期过程中大腿姿态角度随时间的变化曲线为：
23.θ
t
＝p1t3+p2t2+p3t p424.式中：θ
t
为大腿姿态角度，p1、p2、p3、p4为大腿姿态角度变化曲线的参数；
25.s32：设定p1、p2、p3、p4的初始值；
26.s33：基于当前的大腿姿态角度随时间的变化曲线计算大腿姿态角度θ
t
，再计算摆动期内各时刻髋关节x轴位移h
x
和髋关节y轴位移h
y
：
27.h
x
＝a
x
+l
s
·
sinθ
s
+l
t
·
sinθ
t
28.h
y
＝a
y
+l
s
·
cosθ
s
+l
t
·
cosθ
t
29.式中：l
s
为用户小腿长度，l
t
为用户大腿长度；
30.s34：计算对侧腿姿态角度θ
f
：
[0031][0032]
式中：s为跨步长；
[0033]
s35：计算对侧髋关节位移：
[0034][0035]
h
oy
＝l
·
cosθ
f
[0036]
式中：h
ox
、h
oy
分别为对侧髋关节x轴位移、侧髋关节y轴位移；
[0037]
s36：计算两侧髋关节位移差：
[0038]
[0039]
式中：δh为两侧髋关节位移差，max为求最大值的函数，其输出值为摆动期内不同时刻的最大值；
[0040]
s37：调整p1、p2、p3、p4的值，并重复步骤s33至s36，使用优化算法求得满足δh最小的p1、p2、p3、p4的最优值并代入所述三次多项式中，得到最终的大腿姿态角度随时间的变化曲线，重复步骤s33计算得到髋关节x轴位移、髋关节y轴位移作为最终解。
[0041]
作为优选，所述优化算法为梯度下降法。
[0042]
作为优选，所述用于可穿戴式设备的人体髋关节位移测量方法包括相连的单片机和惯性传感器，惯性传感器包含三维加速度计以及三维角速度计。
[0043]
作为优选，所述惯性传感器为基于mpu6050芯片的惯性传感器。
[0044]
作为优选，所述两个可穿戴式设备佩戴于两小腿脚踝外侧上方2～4cm。
[0045]
作为优选，所述两个可穿戴式设备的采样频率不低于100hz。
[0046]
上述各优选方式中的技术特征在不相互冲突的前提下，均可进行相互组合，不构成限制。
[0047]
本发明相对于现有技术而言，其有益效果是：
[0048]
1)使用本发明进行髋关节位移测量，不受场地限制，便捷，易于推广。
[0049]
2)本发明适应多种病态的步态，拥有较好的应用价值以及广泛的应用范围。
[0050]
3)通过可穿戴式设备进行测量，可以方便应用于日常生活。
附图说明
[0051]
图1为本发明中人体走路步态周期示意图；
[0052]
图2为本发明中矢状面示意图；
[0053]
图3为本发明中x轴、y轴、传感器放置位置示意图；
[0054]
图4为本发明中支撑期内髋关节位移计算示意图；
[0055]
图5为本发明中摆动期内髋关节位移计算示意图；
[0056]
图6为本发明中髋关节位移计算结果示意图；
[0057]
上述图1中，a为右腿脚落地动作，b为左腿脚离地动作，c为左腿脚落地动作，d为右腿脚离地动作；
[0058]
上述图3、图4、图5中，标记1为放置在右小腿上的惯性传感器；
[0059]
上述图4、图5中，标记2为小腿、标记3为大腿，标记4为对侧腿。
具体实施方式
[0060]
下面结合附图对本发明进行进一步说明，因便于更好地理解。本发明中的技术特征在不相互冲突的前提下，均可进行相互组合，不构成限制。
[0061]
本发明使用可穿戴式设备对人体髋关节位移进行测量。以某一用户为例，本发明中用于可穿戴式设备的人体髋关节位移测量方法具体实施过程如下：
[0062]
(1)准备工作：
[0063]
本实施例中，目标用户穿戴该可穿戴式设备在右侧小腿脚踝外侧上方3cm处，用于测量用户右侧髋关节位移，如图3所示。同样的，在用户左侧小腿的脚踝外侧上方3cm处也佩
戴有另一个可穿戴式设备，用于测量用户左侧髋关节位移，两侧的髋关节位移测量方法完全相同，下面主要以右腿为例进行说明。
[0064]
该可穿戴式设备包含用于处理信号数据的单片机以及基于mpu6050芯片的惯性传感器，惯性传感器包含三维加速度计以及三维角速度计，采样频率为100hz。
[0065]
用户在平坦的地面上行走，惯性传感器采集加速度、角速度数据；单片机基于内置算法通过加速度、角速度数据识别用户行走过程中的右腿的支撑期、摆动期、平足期和步态周期，并测量右侧小腿姿态角度、跨步长、踝关节x轴位移、踝关节y轴位移。下面以右腿其中一个步态周期为例进行接下来的计算步骤，同时踝关节x轴位移、踝关节y轴位移分别以该步态周期平足期时踝关节的x轴位移、踝关节y轴位移为零点。
[0066]
(2)支撑期内髋关节位移计算：
[0067]
假设用户的腿处于支撑期时膝关节无弯曲，小腿姿态角度等于大腿姿态角度，如图4所示，由此即可通过小腿姿态角度、踝关节x轴位移、踝关节y轴位移计算支撑期内各时刻髋关节x轴位移和髋关节y轴位移：
[0068]
h
x
＝a
x
+l
·
sinθ
s
[0069]
h
y
＝a
y
+l
·
cosθ
s
[0070]
式中：h
x
、h
y
分别为髋关节x轴位移、髋关节y轴位移，a
x
、a
y
分别为踝关节x轴位移、踝关节y轴位移，l为用户腿长(小腿与大腿长度之和)，θ
s
为小腿姿态角度。
[0071]
(3)摆动期内髋关节位移计算：
[0072]
摆动期内髋关节位移计算的原理如下：基于小腿姿态角度、跨步长、踝关节x轴位移和踝关节y轴位移，首先用待定的多项式来表征人体摆动期过程中大腿姿态角度随时间的变化曲线，利用处于摆动期内的当前腿对应的大腿姿态角度、小腿姿态角度、大腿长度和小腿长度，计算出摆动期内髋关节的x轴和y轴位移；利用人体在所述当前腿处于摆动期时身体对侧腿处于支撑期且身体对侧膝关节无弯曲的假设，计算所述当前腿的对侧腿姿态角度；再利用对侧腿姿态角度和对侧腿的长度，计算对侧腿的髋关节运动x轴和y轴位移，将两腿的髋关节位移做差，利用人体运动能量总是取最小原理，迭代优化所述多项式的待定系数，直至两腿髋关节位移差最小，利用优化迭代后的大腿姿态角度随时间的变化曲线，求得当前用户摆动期内的髋关节x轴位移和髋关节y轴位移。
[0073]
在本实施例中，摆动期内髋关节位移计算过程具体如下：
[0074]
1)基于三次多项式定义人体摆动期过程中大腿姿态角度随时间的变化曲线为：
[0075]
θ
t
＝p1t3+p2t2+p3t+p4[0076]
式中：θ
t
为大腿姿态角度，p1、p2、p3、p4为大腿姿态角度变化曲线的参数。
[0077]
2)设定p1、p2、p3、p4的初始值，本实施例中选定p1、p2、p3、p4的初始值均为0。
[0078]
3)基于当前的大腿姿态角度随时间的变化曲线计算大腿姿态角度θ
t
，再计算摆动期内各时刻当前腿侧的髋关节x轴位移h
x
和髋关节y轴位移h
y
：
[0079]
h
x
＝a
x
+l
s
·
sinθ
s
+l
t
·
sinθ
t
[0080]
h
y
＝a
y
+l
s
·
cosθ
s
+l
t
·
cosθ
t
[0081]
式中：l
s
为用户小腿长度，l
t
为用户大腿长度；
[0082]
4)在当前腿处于摆动期时，身体对侧腿处于支撑期，因此假设身体对侧膝关节无弯曲，计算对侧腿姿态角度θ
f
：
[0083][0084]
式中：s为跨步长。
[0085]
5)根据对侧腿姿态角度计算身体对侧的髋关节位移：
[0086][0087]
h
oy
＝l
·
cosθ
f
[0088]
式中：h
ox
、h
oy
分别为对侧髋关节x轴位移、侧髋关节y轴位移；
[0089]
6)计算两侧髋关节位移差：
[0090][0091]
式中：δh为两侧髋关节位移差；max为求最大值的函数，由于中的参数均是随时间t变化的，因此函数的输出值为摆动期内不同时刻的最大值。
[0092]
7)调整p1、p2、p3、p4的值，并重复步骤3)至6)，使用优化算法求得满足δh最小的p1、p2、p3、p4的最优值并代入前述的三次多项式中，得到最终的大腿姿态角度随时间的变化曲线，然后再次重复步骤3)计算即可得到髋关节x轴位移、髋关节y轴位移作为最终解。本实施例中优化算法为梯度下降法。
[0093]
(3)髋关节位移计算效果：
[0094]
本实施例中，上述方法可以准确测量得到用户的髋关节位移。举例而言，该用户小腿长为0.45米，大腿长为0.45米，右腿步态周期内右侧髋关节位移的计算结果如图6所示，由此可见，本发明的装置和方法相对于现有技术而言，可实现对用户的髋关节位移进行测量。
[0095]
以上所述的实施例只是本发明的一些较佳的方案，然而其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。例如，上述实施例也可以使用其他算法或者使用其他传感器来计算所述小腿姿态角度、大腿姿态角度、跨步长以及踝关节x轴、踝关节y轴位移，如使用毫米波雷达、红外传感器、激光雷达等传感器。上述可穿戴设备也可以同时使用两套，以达到同时测量左右两侧髋关节位移的目的。上述可穿戴设备也可以采用现有技术中的其他结构或方式进行改动，如使用别的惯性传感器芯片、采用更高的采样频率等等。
[0096]
由此可见，凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。