一种行人步频计算系统及其方法
【技术领域】
[0001] 本发明涉及到行人步态测量领域，特指一种行人步频计算系统及其方法。
【背景技术】
[0002] 行人步频是基于行人的步态过程，通过各类传感器计算行人所完成的步态周期个 数，进而计算得到的行人步行频率。行人步频是行人步态分析、行人航迹推算等应用中的基 础要素，在康复医疗、体育科学以及行人导航定位等领域有着广泛的应用。
[0003] 行人步频计算最重要的是准确计算出行人步行过程中的每个步态周期，对步态周 期的计算有很多常规的方法：如论文"基于数字视频和数字图像处理的步态分析在卒中康 复训练中的应用"《2009年中国脑血管病杂志，胡雪艳，456页-460页》一文中，基于图像对 步频进行计算，在图像采集阶段存在标记点标定及摄像机定位等问题，操作较为复杂，而 且视频信号数据量较大很难实现实时步频计算；如论文"基于GPS和自包含传感器的行人 室内外无缝定位算法研宄"《中国博士学位论文全文数据库，陈伟，61页-69页》一文中基 于小腿EMG肌电信号计算行人步频，如论文"基于惯性传感器和WiFi的室内定位系统设计 和实现"《中国优秀硕士学位论文全文数据库，张世哲，17页-26页》一文中基于加速度传感 器计算行人步频，这类步频计算系统中一般基于传感器信号采用峰值计算法、过零计算法、 平带计算法等计算行人步频，在测量信号具有噪声情况下，行人步频计算误差较大，而且需 要对测量信号连续不断的进行搜索和计算，算法效率偏低。
[0004] 常规的行人步频算法在测量信号具有噪声情况下，针对信号的峰值计算、过零计 算以及平带计算等误差较大，在测量信号具有噪声情况下存在行人步频计算误差较大的问 题；而且常规的行人步频算法需要对测量信号连续不断的进行计算，存在算法效率偏低的 问题。
【发明内容】
[0005] 本发明要解决的技术问题就在于：针对现有技术存在的技术问题，本发明提供一 种无需连续计算、在测量信号具有噪声情况下也能实现准确、高效的行人步频计算的系统 及其方法。
[0006] 为解决上述技术问题，本发明提出的解决方案为：
[0007]本发明提出一种行人步频计算系统，它包括腰部测量装置1、大腿部测量装置2、 小腿部测量装置3、足面测量装置4以及网关5,其中，腰部量装置1和大腿部测量装置2用 于测量行人髋关节角度，大腿部测量装置2和小腿部测量装置3用于测量行人膝关节角度， 小腿部测量装置3和足面测量装置4用于测量行人踝关节角度，通过无线通信将各个测量 数据发送至网关5,在网关5中进行行人步频计算和显示，其特征在于：
[0008] 所述的腰部测量装置1安装于行人腰部位置，包括腰部惯性测量模块11、腰部嵌 入式微处理器12、腰部无线数据收发模块13,其中，
[0009] 腰部惯性测量模块11用于读取测量行人腰部的旋转角速度数据和加速度数据， 并传送给腰部嵌入式微处理器12;
[0010] 腰部嵌入式微处理器12用于接收腰部惯性测量模块11测量的旋转角速度数据以 及加速度数据；
[0011] 腰部无线数据收发模块13接收或发送旋转角速度数据以及加速度数据，将其发 送给网关5;
[0012] 所述的大腿部测量装置2安装于行人大腿部位置，包括大腿部惯性测量模块21、 大腿部嵌入式微处理器22、大腿部无线数据收发模块23,其中，
[0013] 大腿部惯性测量模块21用于读取测量行人大腿部的旋转角速度数据和加速度数 据并传送给大腿部嵌入式微处理器22;
[0014]大腿部嵌入式微处理器22用于接收大腿部惯性测量模块21测量的旋转角速度数 据以及加速度数据；
[0015] 大腿部无线数据收发模块23接收或发送旋转角速度数据以及加速度数据，将其 发送给网关5;
[0016] 所述的小腿部测量装置3安装于行人小腿部位置，包括小腿部惯性测量模块31、 小腿部嵌入式微处理器32、小腿部无线数据收发模块33,其中，
[0017] 小腿部惯性测量模块31用于读取测量行人小腿部的旋转角速度数据和加速度数 据并传送给小腿部嵌入式微处理器32;
[0018] 小腿部嵌入式微处理器32用于接收小腿部惯性测量模块31测量的旋转角速度数 据以及加速度数据；
[0019] 小腿部无线数据收发模块33接收或发送旋转角速度数据以及加速度数据，将其 发送给网关5;
[0020] 所述的足面测量装置4安装于行人足面位置，包括足面惯性测量模块41、足面嵌 入式微处理器42、足面无线数据收发模块43,其中，
[0021] 足面惯性测量模块41用于读取测量行人足面的旋转角速度数据和加速度数据并 传送给足面嵌入式微处理器42;
[0022] 足面嵌入式微处理器42用于接收足面惯性测量模块41测量的旋转角速度数据以 及加速度数据；
[0023] 足面无线数据收发模块43接收或发送旋转角速度数据以及加速度数据，将其发 送给网关5;
[0024]所述的网关5用于计算行人的步频F(j)，并显示行人的步频F(j)计算结果，网关 5包括网关无线数据收发模块51、网关嵌入式微处理器52、网关显示模块53,其中：
[0025] 所述网关无线数据收发模块51用于接收腰部测量装置1、大腿部测量装置2、小腿 部测量装置3、足面测量装置4发来的各旋转角速度数据以及加速度数据，并发送给网关嵌 入式微处理器模块52;
[0026]所述网关嵌入式微处理器52用于计算行人的步频F(j)，并将行人的步频F(j)发 送给网关显示模块53;
[0027] 所述网关显不t旲块53用于显不彳丁人步频F(j);
[0028] 一种行人步频计算方法，该方法采用行人步频测量系统测量行人步频，其步骤 为：
[0029] (I).将腰部惯性测量模块11测量腰部的竖直方向加速度数据的峰值标识表示为 Gf(J),
[0030] 若Gf (j) = 1，则确定腰部惯性测量模块11测量腰部的竖直方向加速度数据为峰 值，
[0031] 若Gf (j) = 0,则确定腰部惯性测量模块11测量腰部的竖直方向加速度数据为不 是峰值；
[0032] Gf(j) = 0 ；
[0033] 将大腿部惯性测量模块21测量大腿部的前进方向上的加速度数据正负标识表示 SGz(j)，
[0034] gGz(j) = 1，则确定大腿部惯性测量模块21的前进方向上的加速度数据为正值，
[0035] 若Gz (j)=0,则确定大腿部惯性测量模块21的前进方向上的加速度数据为负值， j表不彳丁人步态周期序；
[0036] 初女台化Gz(j)=0;
[0037] (2).分别将腰部测量装置1、大腿部测量装置2、小腿部测量装置3、足面测量装置 4测量的旋转角速度数据和加速度数据发送至网关5;
[0038] (3).采用网关5中的网关嵌入式微处理器52的测量数据，分别计算腰部惯性测量 模块12的竖直方向加速度数据的峰值标识Gf (j)和大腿部惯性测量模块12的前进方向上 的加速度数据正负标KGz (j)，当腰部惯性测量模块12的竖直方向加速度数据达到峰值，即 Gf (j) = 1，贝IJ转步骤（4)，否贝IJ转步骤（2);
[0039] (4).分别设置行人步态双脚支撑期后脚目标函数和前脚目标函数，计算行人步态 周期内后脚的计数值和前脚的计数值，具体方法如下：
[0040] (4-1).当Gf(j) = 1 且Gz(j) =0时，即Gf (j) ? (Gz(j)+1) = 1 时，计算行人步态 周期计时器的计数值Cstart，具体步骤如下：
[0041] (4-1-1).设置行人步态双脚支撑期后脚目标函数，记为Z1 (t)，计算行人步态双脚 支撑期的后脚的关节角度：
[0043] 式（1)中，0Jt)表示利用腰部惯性测量模块11和大腿部惯性测量模块21的旋 转角速度数据和加速度数据计算得到的髋关节角度；0 3(t)表示利用小腿部惯性测量模块 31和足面惯性测量模块41的旋转角速度数据和加速度数据计算得到的踝关节角度；t为行 人步态双脚支撑周期内的时间；
[0044] (4-1-2) ?当[Z1 (t+nk) -Z1 (t-k) ] [Z1 (t-k) -Z1 (t-nk) ] < 0,则认为行人的步态未到 达双脚支撑期内，返回步骤（2)，
[0045]当[Z1 (t+nk)-Z1 (t-k)] [Z1 (t-k)-Z1 (t-nk)]>0,其中，n为小于10的常数，k表示 行人步态双脚支撑周期，则认为行人的步态已到达双脚支撑期内，装在网关嵌入式微处理 器52中的计时器的计数，其计数值记为Cstart;
[0046] (4-2) ?当Gf(j) = 1 且Gz(j) = 1 时，即Gf (j) ? (Gz(j)+1) = 2 时，计算步态周期 计时器的计数值Cmd，具体步骤如下：
[0047](4-2-1).设置行人步态双脚支撑期前脚目标函数，记为Z2 (t)，计算行人步态双脚 支撑期的前脚的关节角度：
[0048] Z2 (t) = (J)J(t)-<i> 2 (t) + <i>3(t) ⑵
[0049] 式⑵中，Mt)表示利用腰部惯