本发明涉及医疗设备领域，尤其涉及一种基于多源信息融合的步频检测方法及设备。
背景技术：
随着我国社会老龄化程度的加剧，由于脑卒中引起的下肢运动功能障碍患者正不断地增加。此外，由于工伤、交通事故、疾病等原因造成精神或肢体损伤的患者也显著增加，严重影响患者的生活质量，并给患者家庭及社会带来极大负担。研究表明，及时进行一定程度的集中和重复训练可以修复受损的中枢神经，有助于下肢运动功能的恢复。传统的康复训练存在医师不足、设备功能单一且成本较高、康复周期长、过程乏味主动性差以及无法进行准确的康复状态评价等缺陷，因此研制能够克服上述缺点的下肢康复训练及评价设备显得尤为必要。
下肢康复机器人作为一种康复医疗设备，通过辅助患者进行科学地、有效地康复训练，从而达到恢复患者运动功能的目的。它在康复训练中为患者提供安全保护的同时，加强患者的主动运动意图对患者运动功能康复具有强化和促进作用，有利于保持患者精神紧张度，加强神经肌肉运动通路的修复。
步频的准确获取是进行康复训练过程中康复评估的重要环节，医生需要根据评估结果实时了解患者状态并进行康复计划调整，患者需要根据康复评估结果了解自身状态，步频作为康复评估重要评估参数之一，其检测准确性就因此显得尤为重要。国内外已经出现不少成熟的产品，例如北京度量科技有限公司的nokov的康复步态系统等，常见的步频步速检测方法包括基于视觉的动作捕捉检测，以及基于姿态传感器、力传感器等数据进行步态的检测，基于视觉的动作捕捉检测是通过计算人体运动特征点的变化频率获取人体步频，基于其他类传感器的步频检测方法多是将传感信号进行低通滤波后，进行傅里叶变换或者短时傅里叶变换，从频域分析人体步频。
上述步频检测技术应用于康复训练的传统方法对下肢康复机器人，尤其是针对于平衡及步行训练的下肢康复机器人具有以下几点不足。
1.基于视觉的动作检测，具有针对性，对硬件要求高，成本高，数据处理复杂，步频准确性以及时间不易保证。
2.基于单传感信号的傅里叶变换或者短时傅里叶变换等的传感信号处理，提取患者步速步频具有非实时性及准确性不足。
技术实现要素：
有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是提供一种具有改进的实时性及准确性的基于多源信息融合的步频检测方法及设备。
为实现上述目的，本发明提供了一种基于多源信息融合的步频检测方法，包括：获取具有随用户人体上下运动的骨盆机构和供用户步行的跑台的系统的骨盆机构高度数据以及跑台电机电流数据；将获取的所述骨盆机构高度数据和所述跑台电机电流数据分别进行卡尔曼滤波，将通过所述卡尔曼滤波获取的数据采用方差权值平均法进行数据融合处理，实时提取出用户的实时步频。
优选地，所述系统还包括传感器，通过所述传感器获取所述骨盆高度数据以及所述跑台电机电流数据。
优选地，所述数据融合处理采用方差权值平均法，包括：
基于数据的瞬态均值μi(k)以及方差
其中，di(k)表示第i个信息源在第k个时刻的值；n表示信息源数据系列的维度；
采用下述公式通过前一时刻的值进行在线更新数据，
采用下述公式执行信息融合：
m表示信息源数量；ωi(k)表示第i个传感器的权值；d(k)表示融合后的信息；
其中根据当前时刻的计算权值ωi(k)：
式中
优选地，针对融合后的信息，选择幅值最大的频率作为实时步频。
为实现上述目的，本发明还提供了一种基于多源信息融合的步频检测设备，包括：随用户人体上下运动的骨盆机构、供用户步行的跑台以及处理装置；其中，处理装置用于执行下述步骤：获取具有随用户人体上下运动的骨盆机构和供用户步行的跑台的系统的骨盆机构高度数据以及跑台电机电流数据；将获取的骨盆机构高度数据和所述跑台电机电流数据分别进行卡尔曼滤波，将通过所述卡尔曼滤波获取的数据采用方差权值平均法进行数据融合处理，实时提取出用户的实时步频。
优选地，利用所述系统本身的传感器获取骨盆机构高度数据以及跑台电机电流数据。
优选地，所述数据融合处理采用方差权值平均法，包括：
基于数据的瞬态均值μi(k)以及方差
其中，di(k)表示第i个信息源在第k个时刻的值；n表示信息源数据系列的维度；
采用下述公式通过前一时刻的值进行在线更新数据，
采用下述公式执行信息融合：
m表示信息源数量；ωi(k)表示第i个传感器的权值；d(k)表示融合后的信息；
其中根据当前时刻的计算权值ωi(k)：
式中
优选地，针对融合后的信息，选择幅值最大的频率作为实时步频。
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。
附图说明
结合附图，并通过参考下面的详细描述，将会更容易地对本发明有更完整的理解并且更容易地理解其伴随的优点和特征，其中：
图1是根据本发明优选实施例的基于多源信息融合的步频检测方法的流程图。
图2是根据本发明优选实施例的基于多源信息融合的步频检测设备的示意图。
需要说明的是，附图用于说明本发明，而非限制本发明。注意，表示结构的附图可能并非按比例绘制。并且，附图中，相同或者类似的元件标有相同或者类似的标号。
具体实施方式
针对康复机器人系统步频检测，现有技术多需要专门配置传感器进行步频的检测，将使得系统成本增加，同时基于傅里叶变换的检测算法具有非实时的特点。而本发明所采用的方法是基于多源信息融合的人体步频检测方法。
图1是根据本发明优选实施例的基于多源信息融合的步频检测方法的流程图。本方法的实施例可基于下肢康复机器人系统，如图2所示，康复机器人系统包括：随人体上下运动的骨盆机构10、跑台20、机器人电控系统30以及显示设备40等。
步频检测作为康复机器人系统功能之一，本方法无需外部传感器，直接依据人机交互过程中，骨盆机构随人体重心上下运动的高度信息，以及人在跑台上行走过程，跑台驱动系统电流输出变化的信息，通过基于信息融合的方法进行步频分析。
如图1所示，根据本发明优选实施例的基于多源信息融合的步频检测方法包括：获取具有随用户人体上下运动的骨盆机构10和供用户步行的跑台20的系统的骨盆机构高度数据以及跑台电机电流数据；将获取的骨盆机构高度数据和所述跑台电机电流数据分别进行卡尔曼(kalman)滤波，将通过所述卡尔曼滤波获取的数据采用方差权值平均法进行数据融合处理，实时提取出用户的实时步频。
例如，利用所述系统本身的传感器获取骨盆机构高度数据以及跑台电机电流数据。
具体地，例如，vma融合方法如下所示：
vma权值参数依赖于数据的瞬态均值μi(k)以及方差
di(k)——第i个信息源在第k个时刻的值；n——信息源数据系列的维度。
所用到的方程与均值通过前一时刻的值进行在线更新
vma融合算法将不同传感器的信息通过加权的形式融合在一起，
m——信息源数量；ωi(k)——第i个传感器的权值；d(k)——融合后的信息
权值ωi(k)根据当前时刻的计算得到
式中
当然，ωi(k)的更新方式有很多种，以上仅为其中一种。
小波变换的含义是把某一被称为基本小波(motherwavelet)的函数位移τ后，再在不同尺度α下，与待分析信号x(t)作内积，即
式中α>0，称为尺度因子，其作用是对基本小波
小波函数具有不唯一性，可选择常用的小波基haar小波、mexcianhat小波、morlet小波。
通过在尺度方向加和变换因子，即可得小波频谱。最终选择幅值最大的频率即为步频。
本发明的方法不需要配置外部传感器，根据人体运动过程中，机器人系统动力输出的变化，进行信息融合，通过卡尔曼滤波及小波分析，实现了实时准确地分析人体运动步频信息。本发明不用配置外部传感器，直接利用机器人驱动系统电流信息作为信号源，分析患者运动步频，降低了设备的生产制造成本；本发明通过骨盆机构和跑台两路信号源的输入，进行信息融合,获得更为精确可靠的信息，能更好的保证了结果的可靠性与准确性；本发明基于卡尔曼滤波和小波分析的方法，能实时分析提取出患者的实时步频，方便治疗师把控治疗过程以及康复疗效评估。
在具体应用时，本发明可适用于一种基于骨盆机构支撑机构的下肢康复机器人，可检测患者运动步频，作为医生对患者康复状态评估的直接参考。本发明属于但不局限于康复机器人技术，而是也可以应用至下肢物理治疗和步态训练等领域。
例如，可以利用诸如下肢康复机器人中的处理系统执行上述步骤。由此，可以看出，本发明还公开了与上述方法对应的设备及系统。
上述说明示出并描述了本发明的优选实施例，如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。