2.沈阳大学 信息工程学院, 辽宁 沈阳 110044
2.School of Information Engineering, Shenyang University, Shenyang 110044, China.
羽毛球运动是一项技术主导型体育竞技运动,具有对抗性强、运动姿态灵活度高等特点.挥拍动作具有灵活性高、重复性弱的特点,从而导致运动员无法客观衡量自己每次挥拍动作的技术特点,只能靠大规模重复动作,自己领悟技术要领.
现有羽毛球挥拍轨迹识别主要采用计算机视觉技术.比如,使用三维高速摄像机[1]记录运动员具体挥拍动作,并根据视频内容计算相关技术参数.或通过普通摄像机获得运动员动作,并采用马尔科夫链随机场图对运动员动作进行识别与预测[2].因此,本文采用一种价格低廉、算法简单、安装方便的智能设备(由加速度计、陀螺仪和磁力计组成)辅助日常羽毛球训练.
使用智能设备进行挥拍轨迹识别,需要解决球拍空中姿态解算和球拍位移计算的问题.羽毛球拍挥拍动作属于三维空间姿态,而陀螺仪漂移、加速度计和磁力计的低频响应[3]、地球重力加速度的存在影响姿态解算精度.为了消除重力干扰,Faber等[4]使用一种剩余矩阵的方法在平稳运动过程中,将偏航角误差减少到2°以内.而对于挥拍轨迹识别的另一重点位移计算来说,有学者通过在鞋尖内置惯性传感器,对步行位移进行2D轨迹估计[5].现有研究结果表明,二维空间内,轨迹误差可以控制在10%以内.通过多个位置布置传感器,垂直位移误差可以达到30%[6].在羽毛球运动中,通过在身体不同部位放置传感器,可以识别多种羽毛球挥拍动作,但并无挥拍轨迹识别计算.
羽毛球运动中的挥拍动作具有姿态转动的变加速运动特点,导致偏航角动态解算误差增大,从而对后续挥拍轨迹识别计算造成更大累积误差.因此,本文提出一种多传感器融合的挥拍轨迹识别方法TRM.该方法使用具有多传感器的智能设备采集挥拍过程中的传感器数据,实现高精度挥拍轨迹识别.通过多次实物实验,验证了本文算法的有效性和高轨迹识别精度.
1 挥拍轨迹的姿态角估计本文通过增加陀螺仪误差作为系统的状态变量和增加加速度计与磁力计计算的姿态角作为观察量输入,使得输出的姿态角更为精确和稳定,得到改进基于卡尔曼滤波器的四元数姿态解算算法.
把陀螺仪观测量和姿态角作为观测量,更新状态变量,建立系统观测方程:
(1) |
(2) |
其中: wx,wy,wz为由陀螺仪观测量; R为测量的噪声; φroll表示横滚角; θpitch表示俯仰角; ψyaw表示偏航角,它们组成一组姿态角.采用椭球拟合校正加速度计和磁力计数据,得到系统观测量输入.
在四元数解算的基础上增加了陀螺仪的漂移误差b,误差估计与修正均由卡尔曼滤波器[7]给出.系统建立状态方程:
(3) |
(4) |
其中: A为由陀螺仪偏导数求得的系统系数矩阵; Qk为系统的过程协方差,即系统过程噪声对应的协方差; q0,q1,q2,q3为四元数; bx,by,bz为陀螺仪三个轴向误差.
在状态方程和观测方程中,实现对陀螺仪误差的修正.采用卡尔曼滤波器来估计球拍的姿态变化,并通过卡尔曼滤波器解算出修正四元数,最终求解出球拍的姿态角.由于陀螺仪的采样数据在球拍变加速挥拍过程中产生抖动误差,导致状态变量X中的陀螺仪误差失真,本文采用加速度计与磁力计对陀螺仪的误差进行估计和修正,经过卡尔曼滤波处理后,取此时的四元数解算球拍的姿态角.
2 频域数值积分轨迹计算 2.1 频域数值积分频域的积分原理就是将原始信号做傅里叶变换后,对频域信号进行积分处理,从而得到相应频谱,再经过傅里叶逆变换进行时域信息还原.
在加速度计数据采样过程中,球拍的持续运动和空中反转,会在加速度计的3个轴向上产生重力分量.为了避免重力分量对后续数值积分的影响[8],本文通过旋转矩阵方法进行重力分量去除.
将提纯后加速度信号进行傅里叶变换,再对变换后的加速度信号在频域内做两次积分可得速度v和位移s:
(5) |
(6) |
式中: w为傅里叶分量所对应的频率; fmin和fmax为最小截止频率和最大截止频率; H(f)为频率响应函数,f为频率,当f∈(fmin,fmax)时,H(f)=1,否则H(f)=0. 对速度v和位移s做傅里叶逆变换即可得到时域内的速度与位移信号.频域积分的误差主要来源于低频,本文结合MPU6050器件手册和多次实验分析,选取最小截止频率为100Hz.
2.2 积分累积误差处理本文采用最小二乘法拟合多项式算法对位移计算的累积误差进行处理.由于每个位移数据都是根据前一次积分得到的,趋势项的存在产生积分后误差,而积分的累积误差导致后续轨迹失真变大.
设积分后的速度为vi(i=0,1,2,…,n),时间为t,假设存在多项式f(t)与vi之间存在均方差之和最小值,则有
(7) |
其中: pk为多项式系数; Ψ为函数集合,它包含所有次数不超过m(m <n)的多项式.由上述假设关系可得
(8) |
根据最小二乘法,I能否取得最小值完全取决于多项式系数pk的取值,由元函数的极值条件可知
(9) |
此时,I能够取得最小值,通过矩阵运算可以求解系数pk,进而通过拟合多项式f(t)消除趋势项,如式(10) 所示:
(10) |
同理,在经过第二次数值积分后所产生的趋势项也可以通过该方法去除.
3 多传感器轨迹识别方法本文提出的多传感器融合挥拍轨迹识别方法TRM,采用空间坐标的形式进行三维轨迹计算,结合拍体姿态角进行挥拍轨迹空间描述.选定初始点为直角坐标系原点,分别对x,y,z三轴的加速度( ax,ay,az)做累积积分即可得到x,y,z三轴的位移量(sx,sy,sz),每个采样点的位移量即可构成一组空间坐标,这样挥拍的三维轨迹就以坐标的形式呈现在空间中.
羽毛球拍挥拍轨迹识别需要5个要素,包括球拍空中运动3个方向的角度,即姿态角、球拍挥拍速度和挥拍位移.
首先,将通过传感器采集图 1中的1~3模块的陀螺仪信号、磁力计信号和加速度信号,并提供给改进基于卡尔曼滤波器的四元素姿态解算算法进行姿态解算,得到球拍空中3个姿态角φroll,θpitch,ψyaw.其中,姿态解算过程中使用陀螺仪观测量与加速度计和磁力计估计得到的姿态角共同作为观测方程参数,使用四元数和陀螺仪误差共同作为状态方程参数.再使用卡尔曼滤波器进行姿态角的计算,从而减少挥拍过程中产生的陀螺仪误差,达到误差修正的目的.
其次,如图 1所示,对加速度信号进行频域数值积分轨迹计算.包括去除3个轴向的重力分量的加速度信号提纯; 对提纯后加速度信号进行频域积分,一次积分得到挥拍速度参数,二次积分得到挥拍位移参数; 对频域积分过程中的累积误差使用最小二乘多项式的方法进行处理,以减少轨迹误差.最后,通过球拍姿态角、速度和位移得到羽毛球拍挥拍动作的三维空间轨迹.
为了不破坏球拍原有平衡点,将智能设备固定于球拍拍柄底部,智能设备如图 2虚线方框中所示.选用全碳钢球拍,质量为81.5g(包括智能设备).智能设备包括集成加速度计和陀螺仪的MPU6050芯片模块及HMC5883磁力计芯片模块.智能设备的采样频率为100Hz.
如图 2所示,球拍坐标与世界坐标关系图.图 2中,世界坐标用X,Y,Z表示.沿拍底指向拍头方向为球拍x轴,以拍体为中心,y轴垂直向内,z轴垂直向外.智能设备放置在球拍尾部,通过蓝牙将数据传输给电脑,在电脑端进行数据存储和球拍轨迹识别的计算.
4.2 轨迹识别有效性验证测试者使用4.1节中的实验设备,分别做正手高远球和反手挑球动作(正反手经典动作),以验证TRM方法的有效性.每种动作分为两种速度,一种为匀速,一种为变加速.如图 3所示,从左至右分别为正手高远球和反手挑球动作.
轨迹识别的有效性验证共进行4组实验,每组实验重复25次.本文采用有效率对每组实验进行有效性评价,有效率η公式为
(11) |
其中:Ai为识别出动作与测试者所做动作一致的次数; B=25为总测试次数.
表 1为有效性验证对照表.该表按照不同挥拍类型和不同挥拍速度共分为4种条件验证.根据表 1所示,匀速条件下的正手高远球和反手挑球有效率均为96.00%,变加速条件下的正手高远球和反手挑球有效率分别为92.00%和90.00%.变加速条件下挥拍动作容易失真,智能设备识别出的轨迹与测试者实际挥拍轨迹有一定差距,特别是反手动作,陀螺仪在变加速转向时计算不够精准.而在匀速挥拍过程中,轨迹估计相对稳定,有效率均高达96.00%.其中,4种条件下的100次挥拍实验有效率均值为93.00%,从而说明本文提出的TRM对羽毛球拍挥拍轨迹可以进行有效识别.
轨迹识别的准确性可以有效地评估轨迹识别算法好坏,而轨迹识别误差率又是对准确性的定量描述.因此,本文通过计算实际挥拍动作与TRM识别出动作之间的3个轴向轨迹识别误差率均值,来评估TRM的准确性.轨迹识别误差率ω均值公式如下:
(12) |
(13) |
式中: ai ,bi,ci分别为球拍X轴,Y轴和Z轴的轨迹误差; Kin为各轴轨迹误差,当n=1时,Ki1=ai;当n=2时,Ki2=bi;当n=3时,Ki3=ci; Ci为TRM计算值; mi为激光测距仪测量值.
在4.1节中,测试者进行了两种速度条件下的正手高远球和反手挑球动作测试.每种条件进行25次,共计100次实验.本节中,针对这100次实验结果,进行位移误差率计算.以人身体后白墙垂直向上方向为Z轴正方向,人竖直站立正前方为X轴正方向,平行墙面右手方向为Y轴正方向.采用激光测距仪测量挥拍动作过程中的X,Y,Z三个方向的位移.以正手高远球挥拍动作为例,如图 4所示,位移分别为X轴移动61mm,Y轴移动201mm,Z轴移动51mm.
为了评价本文提出的TRM算法轨迹识别准确性,本实验中,以激光测距仪测量实际挥拍位移为基准,分别采用时域积分方法[9-10]与TRM进行挥拍位移估计,并计算3个轴向误差率.其中,时域数值积分算法选取复化辛普森公式进行积分,在经多项式拟合消除累积误差后,提取球拍的运动轨迹.分别采用TRM和时域积分方法对4.1节中的100次实验进行挥拍轨迹估计.通过式(12) 计算轨迹识别误差率,TRM和时域积分轨迹识别方法的3个轴向轨迹误差如表 2所示.
实验结果表明,本文提出的TRM轨迹识别方法无论在单个轴向上的位移误差,还是3个轴向位移误差的均值,误差率都比使用时域积分的轨迹识别方法小,并且TRM轨迹识别误差均值在10%以内,低于现有方法均值水平.
5 结论本文针对羽毛球拍具体应用,结合羽毛球拍挥拍动作运动特点,提出一种多传感器融合的羽毛球拍挥拍轨迹识别方法.该方法通过加速度计和磁力计对陀螺仪进行修正,改进现有姿态解算方法.通过频域积分计算挥拍位移,经过典型挥拍动作测试结果表明,挥拍偏航角姿态角精度在1°以内,轨迹误差均值控制在10%以内.羽毛球拍挥拍轨迹识别精度得到明显的提高,为羽毛球动作识别技术的改进提供了高效的解决方案.
[1] |
张博, 杨志平.
羽毛球男子头顶突击杀球技术动作的运动学分析[J]. 沈阳体育学院学报 , 2013, 32 (2) : 102–106.
( Zhang bo, Yang Zhi-ping. Kinematical analysis on the skills of surprising scores overhead in men’s badminton[J]. Journal of Shenyang Sport University , 2013, 32 (2) : 102–106. ) |
[2] | Wang Z,Shi Q,Shen C,et al.Bilinear programming for human activity recognition with unknown MRF graphs[C]// Proceeding of IEEE Conference on Computer Vision and Pattern Recognition (CVPR).Portland,2013:1690-1697. |
[3] | Ahmadi M, Khayatian A, Karimaghaee P. Attitude estimation by divided difference filter in quaternion space[J]. Acta Astronautica , 2012, 75 (6/7) : 95–107. |
[4] | Faber G S, Chang C C, Rizun P, et al. A novel method for assessing the 3-D orientation accuracy of inertial/magnetic sensors[J]. Journal of Biomechanics , 2013, 46 (15) : 2745–2751. DOI:10.1016/j.jbiomech.2013.07.029 |
[5] | Stenum J, Bruijn S M, Jensen B R. The effect of walking speed on local dynamic stability is sensitive to calculation methods[J]. Journal of Biomechanics , 2014, 47 (15) : 3776–3779. DOI:10.1016/j.jbiomech.2014.09.020 |
[6] | Ibata Y,Kitamura S,Motoi K,et al.Measurement of three dimensional posture and trajectory of lower body during standing long jump utilizing body-mounted sensors[C]//Proceeding of 35th Annual International Conference of the IEEE EMBS.Osaka,2013:4891-4894. |
[7] | Jahanchahi C, Mandic D P. A class of quaternion Kalman filters[J]. IEEE Transactions on Neural Networks and Learning Systems , 2014, 25 (3) : 533–544. DOI:10.1109/TNNLS.2013.2277540 |
[8] | Rodriguez-Martin D, Samà A, Perez-Lopez C, et al. SVM-based posture identification with a single waist-located triaxial accelerometer[J]. Expert Systems with Applications , 2013, 40 (18) : 7203–7211. DOI:10.1016/j.eswa.2013.07.028 |
[9] | Rezvanian S, Lockhart T E. Towards real-time detection of freezing of gait using wavelet transform on wireless accelerometer data[J]. Sensors , 2016, 16 (4) : 475–480. DOI:10.3390/s16040475 |
[10] | Gangsig S, Ohseop S. A time-domain method to generate artificial time history from a given reference response spectrum[J]. Nuclear Engineering and Technology , 2016, 48 (3) : 831–839. DOI:10.1016/j.net.2016.01.023 |