CN105426842B - 基于支持向量机和表面肌电信号的多类手部动作识别方法 - Google Patents

基于支持向量机和表面肌电信号的多类手部动作识别方法 Download PDF

Info

Publication number
CN105426842B
CN105426842B CN201510801198.8A CN201510801198A CN105426842B CN 105426842 B CN105426842 B CN 105426842B CN 201510801198 A CN201510801198 A CN 201510801198A CN 105426842 B CN105426842 B CN 105426842B
Authority
CN
China
Prior art keywords
feature
data
training
collection
test
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.)
Expired - Fee Related
Application number
CN201510801198.8A
Other languages
English (en)
Other versions
CN105426842A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201510801198.8A priority Critical patent/CN105426842B/zh
Publication of CN105426842A publication Critical patent/CN105426842A/zh
Application granted granted Critical
Publication of CN105426842B publication Critical patent/CN105426842B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/20—Movements or behaviour, e.g. gesture recognition
    • G06V40/28—Recognition of hand or arm movements, e.g. recognition of deaf sign language
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing

Abstract

本发明公开了一种基于支持向量机的表面肌电信号多类手部动作识别方法。方法的主要步骤如下：1)获取肌电数据，对信号进行平滑滤波，并通过不同尺度的采样窗口生成数据样本；2)以数据样本为单位，从每个数据样本中提取出一个包含19种时域，频域以及时频域特征的新型多特征特征集，并对特征向量进行归一化和基于最小冗余最大相关性准则的特征选择；3)设计基于皮尔逊VII广义核的支持向量机分类器，使用基于交叉验证的粗网格搜索优化算法，对支持向量机的参数进行优化；4)使用训练集中的数据样本和步骤3)参数优化过程得到的分类器最优参数，训练分类模型，并将测试集中的数据样本输入分类模型进行分类测试。

Description

基于支持向量机和表面肌电信号的多类手部动作识别方法
技术领域
本发明涉及一种基于支持向量机和表面肌电信号的多类手部动作识别方法，属于模式识别技术领域。
背景技术
特征提取对于肌电信号模式识别问题的最终识别率有着至关重要的影响。肌电信号特征可以分为时域，频域，时频域三种，其中时域特征包括一些基于信号幅值的特征，频域特征包括一些基于信号功率谱的特征，而时频域特征包括一些经过小波分析技术提取出的特征。时域特征存在对非静态信号不鲁棒，对信号幅值变化较为敏感等缺陷，频域特征存在对经过一些预处理步骤(如全波校正)后的信号不鲁棒的缺陷，此外一些文献表明单一频域特征性能较差。本文通过构建包括时域，频域和时频域特征的高性能混合多特征集，来克服不同域单一特征的固有缺陷。
多通道肌电信号提取出的多特征集会存在高维度的特性，维度过高会导致模型推广能力的下降，同时极大地增加计算负载。一般通过特征降维来解决数据维度过高的问题，常见的数据降维方法主要分为三类：特征选择，特征映射和特征聚类。特征选择试图从原特征集中选取一个更具有代表性的特征子集，特征映射通过将特征从原始高维空间映射到一个低维空间来降低特征维度，特征聚类从原始特征产生多个聚类，用聚类质心的来取代原始特征集，作为一种新的低维度表述形式。本文通过一种基于互信息的最小冗余最大相关性特征选择方法来完成特征空间降维。
模式识别方法的选择同样在很大程度上影响肌电模式识别问题的最终准确率。已有工作表明支持向量机在肌电模式识别方面有着较好的性能。在对Atzori等人发布的NinaPro数据集基准测试结果中，非线性支持向量机是多种线性和非线性分类器中唯一一种能够在变换特征情况下依然能够持续获得高识别率的分类器。国内已有涉及支持向量机进行肌电模式识别的发明，大多使用较为经典的线性核、径向基核、多项式核、sigmoid核等等核函数。等人提出了一种基于皮尔逊VII函数的广义核函数，本发明通过使用皮尔逊VII广义核函数，在一定程度上提高了识别准确率。
发明内容
本发明的目的在于针对现有技术的不足，提供一种基于支持向量机和表面肌电信号的多类手部动作识别方法，通过构建一种新型特征集和引入皮尔逊VII广义核函数，提高了手部动作识别的准确性。
本发明的目的是通过以下技术方案来实现的：一种基于支持向量机和表面肌电信号的多类手部动作识别方法，包括如下步骤：
(1)数据获取，预处理，分割，训练集测试集的划分，包括以下子步骤：
(1.1)获取多通道肌电数据；
(1.2)对数据进行平滑滤波；
(1.3)使用100ms，150ms，200ms，250ms四种长度的采样窗口对数据进行采样分割，其中采样窗口的移动步长为窗口长度的25％；
(1.4)根据评估方法，将数据样本划分为训练集和测试集；
(2)特征提取，归一化和自动特征选择，包括以下子步骤：
(2.1)对每个采样窗口内的数据，提取一种多特征的特征集；
(2.2)对特征进行归一化处理；
(2.3)应用基于最小冗余最大相关性(Minimum Redundancy Maximum Relevance,MRMR)准则的自动特征选择算法进行特征选择；
(3)基于皮尔逊VII广义核支持向量机的多类手部动作识别，包括以下子步骤：
(3.1)设计基于皮尔逊VII广义核(Pearson VII Universal Kernel，PUK)的支持向量机分类器；
(3.2)使用训练集数据样本，展开网格搜索参数优化过程，搜索最优的分类器参数；
(3.3)使用步骤(3.2)得到的最优分类器参数和所有的训练集数据样本，训练分类器；
(3.4)将测试集数据样本依次输入分类器，输出分类结果。
进一步地，步骤(1.1)中，数据来源为NinaPro开源数据集1内的10通道肌电数据。
进一步地，步骤(1.2)中，使用50毫秒长度的窗口对数据进行均值平滑滤波。
进一步地，步骤(2.1)中的特征集包含19种特征，它们是：信号幅值的绝对值之和(IEMG)，信号幅值的绝对平均值(MAV)，改进信号幅值的绝对平均值1(MMAV1)，改进信号幅值的绝对平均值2(MMAV2)，信号均方根(RMS)，估计肌肉收缩力的非线性检测器(v-order)，估测肌肉收缩力的对数检测器(LOG)，波形长度(WL)，幅值平均变化值(AAC)，波形长度标准差(WL-DASDV)，Willison幅值(WAMP)，自回归系数(ARC)，绝对平均值斜率(MAVSLP)，功率谱的平均功率(MNP)，Daubechies1小波变换边界(MDWT-DB1)，基于Daubechies1小波的多分辨率小波分析(MRWA-DB1)，小波包变换后信号的均值(DWPT-MEAN)，小波包变换后信号标准差(DWPT-SD)，自回归余数残差29种统计量(ARR-29)。
进一步地，步骤(2.2)中对训练特征集和测试特征集采用不同的特征归一化方法：
对训练特征集的归一化方法为，用训练特征集每一列减去其自身的均值并除以其方差；
对测试特征集的归一化方法为，用测试特征集样本每一列减去训练特征集对应列的均值并除以训练特征集对应列的方差；
其具体公式为：
其中train为归一化前的训练特征集，train′为归一化后的训练特征集，test为归一化前的测试特征集，test′为归一化后的测试特征集。
进一步地，步骤(2.3)中所述基于最小冗余最大相关性准则的自动特征选择算法，其具体实现方式为：
设存在一个特征子集S,S中每一列特征xi到类标签c的依赖性D定义为：
S中特征两两之间的冗余R定义为：
上述公式中，I为互信息，两个离散变量X和Y之间互信息的计算公式为：
最小冗余最大相关性准则的自动特征选择算法试图找到一个特征子集S，其满足拥有最大的依赖性D，最小的冗余R，即满足：
maxφ(D,R),φ(D,R)＝D-R
实际的手势识别实验中，根据“在整个分类器的训练过程中测试集信息始终保持未知”的模式识别原则，对训练特征集和测试特征集的特征选择使用不同的策略，具体步骤为：
对训练特征集的特征选择，主要是应用基于最小冗余最大相关性准则的自动特征选择算法，选取最优的N维子集S。将子集中每一维特征对应于原始特征集中的序号保存下来；
对测试特征集的特征选择，直接应用对训练集特征选择时保存下来的序号。
进一步地，步骤(3.1)具体为：采用皮尔逊VII广义核函数作为支持向量机的核函数，对于支持向量机的具体原理描述如下：
对于n个样本的训练样本集其中xi为一个训练样本，yi为样本xi的标签；
设变换φ(·)将输入空间映射到特征空间上，将数据分成两类的线性超平面可以定义为：
ω×φ(x)+b＝0 ω∈Rp,b∈R (1)
支持向量机试图找到一个最优的超平面，可以最好地划分训练数据集中的两类数据，这个优化过程等价于求解下列二次规划问题：
其中α为拉格朗日乘子，C∈[0,∞)为惩罚因子，xi和xj为训练集中两个不同的样本，带有非零α的xi被定义为支持向量，设tj(j＝1,...,s)为s个支持向量的序号，ω可以由求得，通过引入核函数k(xi,xj)＝φ(xi)T·φ(xj)，公式(2)可以由以下形式表述：
公式(3)中的核函数k(xi,xj)有多种形式，线性核k(xi,xj)＝φ(xi)T·φ(xj)可以等价于k(xi,xj)＝xi·xj
另一种常用核函数为径向基核，公式为：
k(xi,xj)＝exp(-γ|xi-xj|2)，其中|·|为欧氏距离，γ为径向基核函数参数。
所述皮尔逊VII广义核函数的定义如下：
其中，σ和λ为皮尔逊VII广义核函数的两个参数。
进一步地，步骤(3.2)具体为：首先确定σ和λ的候选值，以及SVM惩罚因子C的候选值。假设皮尔逊VII广义核参数σ有M个候选值，参数λ有N个候选值，首先产生一个M×N的网格，网格中每个点都是一个候选(σ,λ)组合。随后在网格每个节点上，通过五折交叉验证，即将训练特征集数据分成5等分，每次依次取1份作为测试集，剩余4份作为训练集，依次重复5次，每次都用该节点上的(σ,λ)值作为分类器参数，训练分类器并计算识别结果，最终得到5次测试的识别率，统计5次平均识别率，确定获得最高5折交叉验证识别率的(σopt,λopt)组合。接下来对惩罚因子C进行调优，假设惩罚因子C有K个候选值，则对每个候选值，设定之前获得最优的(σopt,λopt)组合参数，并使用5折交叉验证步骤训练分类器，得到K个候选惩罚因子C各自的5折交叉验证识别率，统计获得最高识别率的惩罚因子Copt，作为最优惩罚因子。最终，输出并且保存上述调优过程获得的最优参数(σopt,λopt,Copt)。
本发明的有益效果是：通过引入新型皮尔逊VII广义核函数，提升了支持向量机对多类手部运动肌电信号的分类性能。对比测试53种时域、频域、时频域特征，构建了四种新型多特征集。通过组合时域、频域、时频域特征构建多特征特征集，克服了单一特征由于其本身性质导致的各种缺陷。对于基于表面肌电信号的假肢控制，人机交互等领域有着重要意义。
附图说明
图1为本发明涉及实验所选取的NinaPro数据集3个手势集，(a)为5个腕部运动的手势集；(b)为8个手部姿态的手势集；(c)为12个手指运动的手势集；
图2为本发明所述方法流程图；
图3为150ms长度采样窗口情况下，使用皮尔逊VII广义核函数进行分类，在不同手势集上的识别结果，(a)为150ms长度采样窗口情况下，使用皮尔逊VII广义核函数进行分类，在不同手势集上的不投票识别准确率；(b)为150ms长度采样窗口情况下，使用皮尔逊VII广义核函数进行分类，在不同手势集上的投票识别准确率，投票窗口长度单位为一段完整的动作；
图4为使用特征集3情况下，不同支持向量机核函数的分类性能对比，(a)为使用特征集3情况下，在识别5个腕部运动时，不同支持向量机核函数分类准确率和采样窗口的关系；(b)为使用特征集3情况下，在识别8个手部姿态时，不同支持向量机核函数分类准确率和采样窗口的关系；(c)为使用特征集3情况下，在识别12个手指运动时，不同支持向量机核函数分类准确率和采样窗口的关系。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
如图2所示，本发明一种基于支持向量机和表面肌电信号的多类手部动作识别方法，包括如下步骤：
步骤(1)下载NinaPro数据集1的10通道肌电数据，并对数据进行平滑滤波，滤波使用50毫秒长度的窗口，平滑滤波具体公式为：
其中为平滑滤波后第i帧的信号值，xi为平滑滤波前第i帧的信号值。
步骤(2)使用滑动采样窗口对数据进行采样，生成数据样本。把数据样本按照具体评估方法划分为训练集和测试集。
步骤(3)以数据样本为单位，从每个数据样本中提取时域，频域以及时频域肌电特征，构建一个新型多特征特征集，具体步骤为：
(3.1)从已有肌电信号相关文献中总结53种时域，频域和时频域特征；
(3.2)将这53种单一特征按照其属性分为四类，这四类分别是时域类，频域类，时频域类和其他类；
(3.3)将时域类特征进一步按照其散点图和数学属性分成四类，这样总计有七个特征类；
(3.4)选取NinaPro数据集1中的52种手势和数据集2中的40种手势，获取肌电信号并按照被试内十折交叉验证方法划分训练集和测试集；
(3.5)使用线性判别分析分类器，统计53种单一特征各自的被试内十折交叉验证准确率；
(3.6)对于七个特征类中的每一类，选取识别率最高的单一特征。同时，将识别率大于50％的单一特征也选出来。由于本发明在NinaPro的两个数据集上分别统计识别率，统计两个优选特征集合，它们分别是NinaPro1上最高识别率的特征集和NinaPro1上识别率大于50％的特征集合的并集(下文中记作NinaPro1best+50％)，NinaPro2上最高识别率的特征集和NinaPro2上识别率大于50％的特征集合的并集(下文中记作NinaPro2best+50％)。随后，取NinaPro1best+50％和NinaPro2best+50％的交集，作为我们提出的新型特征集。该新型特征集共包含19种不同的肌电信号特征；
(3.7)实际测试时，对训练数据集和输入测试数据集的每个样本，均提取3.6)中总结的新型特征集，以构建训练特征集和测试特征集；
步骤(4)对特征向量进行归一化，具体步骤为：
对训练特征集的归一化方法为：用训练特征集每一列减去其均值并除以其方差。
对测试特征集的归一化方法为：测试特征集样本每一列的元素减去训练特征集对应列的均值并除以训练特征集对应列的方差。
其具体公式为：
其中train为归一化前的训练特征集，train′为归一化后的训练特征集，test为归一化前的测试特征集，test′为归一化后的测试特征集。
对训练集的归一化在进行分类测试之前进行。离线测试中，本发明在进行分类测试之前完成对测试集的归一化，实际应用中，本发明推荐在分类测试时，对每个新输入分类器的测试样本分别进行归一化操作。
步骤(5)对归一化后的训练集特征向量进行基于最小冗余最大相关性准则的自动特征选择，最小冗余最大相关性准则为：
设存在一个特征子集S,S中每一列特征xi到类标签c的依赖性D定义为
S中特征两两之间的冗余R定义为：
上述公式中，I为互信息，两个离散变量X和Y之间互信息的计算公式为：
满足最小冗余最大相关性准则的特征子集Sopt，其满足拥有最大的依赖性D，最小的冗余R，即满足：
maxφ(D,R),φ(D,R)＝D-R
从训练集中挑选出最优特征子集Sopt后，将Sopt中每一列特征的序号保存起来，实际测试中，对测试集每一个测试样本，依照Sopt中的序号来选择参与测试的特征向量。
步骤(6)设计基于皮尔逊VII广义核的支持向量机分类器，对支持向量机的定义为：
对于n个样本的训练样本集其中xi为一个训练样本，yi为样本xi的标签，设变换φ(·)将输入空间映射到特征空间上，将数据分成两类的线性超平面可以定义为：
ω×φ(x)+b＝0 ω∈Rp,b∈R (1)
支持向量机试图找到一个最优的超平面，可以最好地划分训练数据集中的两类数据，这个优化过程等价于求解下列二次规划问题：
其中α为拉格朗日乘子，C∈[0,∞)为惩罚因子，带有非零α的xi被定义为支持向量，设tj(j＝1,...,s)为s个支持向量的序号，ω可以由求得，通过引入核函数k(xi,xj)＝φ(xi)T·φ(xj)，公式(2)可以由以下形式表述：
公式(3)为支持向量机的基本形式，其中的核函数k(xi,xj)可以有多种形式，线性核k(xi,xj)＝φ(xi)T·φ(xj)可以等价于：
k(xi,xj)＝xi·xj (4)
本发明采用的新型皮尔逊VII广义核函数，其定义为：
此外，作为和皮尔逊VII广义核函数的对比，本发明涉及的实验中还测试了径向基核函数和线性核函数，线性核函数公式参考公式(4)
径向基核函数公式为：
k(xi,xj)＝exp(-γ|xi-xj|2)，其中|·|为欧氏距离。
本发明使用序列最小最优化支持向量机分类器，使用一对一(OAO)策略解决多分类问题。
步骤(7)使用基于交叉验证的粗网格搜索优化算法，对支持向量机分类器的参数进行优化，粗网格搜索优化具体步骤为：
(7.1)确定σ和λ的候选值，以及SVM惩罚因子C的候选值。假设皮尔逊VII广义核参数σ有M个候选值，参数λ有N个候选值，首先产生一个M×N的网格，网格中每个点都是一个候选(σ,λ)组合；
(7.2)在网格每个节点上，通过五折交叉验证，即将训练特征集数据分成5等分，每次依次取1份作为测试集，剩余4份作为训练集，依次重复5次，每次都用该节点上的(σ,λ)值作为分类器参数，训练分类器并计算识别结果，最终得到5次测试的识别率，统计5次平均识别率，确定获得最高5折交叉验证识别率的(σopt,λopt)组合；
(7.3)对惩罚因子C进行调优过程。假设惩罚因子C有K个候选值，则对每个候选值，设定7.2)中获得最优的(σopt,λopt)组合参数，并使用8.2)中相同的5折交叉验证步骤训练分类器，得到K个候选惩罚因子C各自的5折交叉验证识别率，统计获得最高识别率的惩罚因子Copt，作为最优惩罚因子；
(7.4)输出并且保存上述调优过程获得的最优参数(σopt,λopt,Copt)；
步骤(8)使用训练集中的数据样本和步骤7)得到的最优参数(σopt,λopt,Copt)，训练分类器，并将测试集中的数据样本输入训练好的分类器进行分类测试。
步骤(9)使用多数同意规则投票确定最终测试结果。设每N个数据样本为一个投票窗口，对每个投票窗口内的N个数据样本，统计出现次数最多的标签，作为该窗口的标签。最终识别率计算公式为：
本发明测试不投票的识别结果，和以整段动作为长度的投票窗口的识别结果。
实施例
步骤(1)本发明使用开源NinaPro数据集作为肌电数据来源，选取NinaPro数据集1中5个腕部运动，8个手部姿态，12个手指运动总共25类手势的数据。本发明涉及的手势动作参考图1。
步骤(2)对原始数据采用50ms的窗口进行平滑滤波，并且按照100ms，150ms，200ms，250ms四种长度的滑动窗口来进行采样，滑动窗口的移动步长均为窗口长度的25％。
步骤(3)按照被试内十折交叉验证方法，划分离线算法性能评估的训练集和测试集。NinaPro数据集1总共涉及27名被试，选取每个被试的所有动作数据，分成10份，每次依次取1份作为测试集，剩余9份作为训练集，这样每名被试都会产生10折相互独立的分类测试，总共有270折测试。
步骤(4)从每个采样窗口的数据中提取本发明提出的新型特征集，同时提取4种已知文献中总结的经典特征集作为识别率参考，这4种经典特征集分别为：Phinyomark特征集，改进Phinyomark特征集，Hudgin特征集，Du特征集。
步骤(5)对训练集和测试集分别进行归一化。
步骤(6)对归一化后的特征向量进行基于最小冗余最大相关性准则的自动特征选择，自动特征选择仅对训练集数据进行，在测试时，对测试数据直接应用训练集数据的特征选择结果。自动特征选择算法的输出为训练集特征每一列按照最小冗余最大相关性得分由高到低重新排列的列序号，本发明选取得分最高的D维特征向量，其中D的范围为30,60,100,200,400,600，将选中的特征向量序号保存下来。
步骤(7)设计基于皮尔逊VII广义核函数的支持向量机分类器，同时作为皮尔逊VII广义核函数的性能对比，测试还使用了传统的线性核函数和径向基核函数。
步骤(8)使用基于交叉验证的粗网格搜索优化算法，对支持向量机各个核函数分类器的参数进行优化，其中皮尔逊VII广义核函数参数ω的候选参数集为{27，26，25，...，2-1，2-2，2-3}，参数σ的候选参数集为{27，26，25，...，22，21，20}，惩罚因子C的候选参数集为{25，24，23，22，21，20}。径向基核函数参数γ的候选参数集为{2-3，2-5，2-7，...，2-11，2-13，2-15}，径向基核函数惩罚因子C的候选参数集为{29，27，25，...，21，2-1，2-3}。线性核函数惩罚因子C的候选参数集为{25，24，23，...，2-13，2-14，2-15}。
步骤(9)使用训练集中的数据样本和步骤(8)得到的最优参数，训练分类器，并将测试集中的数据样本输入分类器进行测试。
步骤(10)使用多数同意规则投票确定最终测试结果。本文测试不投票的识别结果，和以整段动作为单位的投票窗口的识别结果，这里整段动作指的是一个包含起势，稳定态和收势的完整手势动作。
步骤(11)设上述步骤所涉及的不同采样窗口长度，不同特征集，不同支持向量机核函数等设置为不同的实验变量，对每种实验变量组合，均进行一次完整的步骤(1)-步骤(10)的实验流程，统计所有实验变量组合的识别率。
150ms采样窗口，三个手势集的实验结果如下表所示：
5个腕部运动
8个手部姿态
12个手指运动
从表中结果可以看出，本发明提出的新型特征集的不投票准确率均超过了Phinyomark，Hudgin和Du三种传统特征集，改进Phinyomark特征集也是一种新型特征集，由Doswald等人在2014年提出，本发明提出的新型特征集在识别8个手部姿态时，有部分识别率超过了改进Phinyomark特征集。经过投票后，本发明提出的新特征集在大部分情况下和其他特征集相比均有一定的优势。
图3(a)、(b)分别表示的是在150ms长度窗口下，使用皮尔逊VII广义核支持向量机作为分类器，在使用不同特征集时所获得的不投票准确率和投票准确率柱状图。
图4(a)、(b)、(c)分别表示的是使用本发明提出的新型特征集时，不同支持向量机核函数以及不同采样窗口长度组合，在识别5个腕部运动，8个手部姿态，12个手指运动时的准确率。可以看出，皮尔逊VII广义核支持向量机的分类性能在所有情况下都超过了线性核函数支持向量机，在大部分情况下超过了径向基核函数支持向量机。
在使用150ms长度采样窗口时，皮尔逊VII广义核支持向量机可以在识别5个腕部动作时获得98.74％的投票准确率，在识别8个手部姿态时可以获得96.48％的投票准确率，在识别12个手指运动时可以获得96.64％的投票准确率。

Claims (7)

1.一种基于支持向量机和表面肌电信号的多类手部动作识别方法，其特征在于，包括如下步骤：
(1)数据获取，预处理，分割，训练集测试集的划分，包括以下子步骤：
(1.1)获取多通道肌电数据；
(1.2)对数据进行平滑滤波；
(1.3)使用100ms，150ms，200ms，250ms四种长度的采样窗口对数据进行采样分割，其中采样窗口的移动步长为窗口长度的25％；
(1.4)根据被试内十折交叉验证评估方法将数据样本划分为训练集和测试集；
(2)特征提取，归一化和自动特征选择，包括以下子步骤：
(2.1)对每个采样窗口内的数据，提取一种多特征的特征集；
(2.2)对特征进行归一化处理；
(2.3)应用基于最小冗余最大相关性准则的自动特征选择算法进行特征选择；
(3)基于皮尔逊VII广义核支持向量机的多类手部动作识别，包括以下子步骤：
(3.1)设计基于皮尔逊VII广义核的支持向量机分类器；具体为：采用皮尔逊VII广义核函数作为支持向量机的核函数，对于支持向量机的具体实现过程如下：
对于n个样本的训练样本集其中xi为一个训练样本，yi为样本xi的标签；
设变换φ(·)将输入空间映射到特征空间上，将数据分成两类的线性超平面定义为：
ω×φ(x)+b＝0ω∈Rp,b∈R (1)
其中ω是垂直于这个线性超平面的向量，b是为了防止超平面过原点而引入的位移；
支持向量机试图找到一个最优的超平面，通过该超平面最好地划分训练数据集中的两类数据，这个优化过程等价于求解下列二次规划问题：
其中α为拉格朗日乘子，C∈[0,∞)为惩罚因子，xi和xj为训练集中两个不同的样本，带有非零α的xi被定义为支持向量，设tj为s个支持向量的序号，j＝1,...,s，ω由求得，通过引入核函数k(xi,xj)＝φ(xi)T·φ(xj)，公式(2)由以下形式表述：
公式(3)中的核函数k(xi,xj)采用皮尔逊VII广义核函数，定义如下：
其中，σ和λ为皮尔逊VII广义核函数的两个参数；
(3.2)使用训练集数据样本，展开网格搜索参数优化过程，搜索最优的分类器参数；
(3.3)使用步骤(3.2)得到的最优分类器参数和所有的训练集数据样本，训练分类器；
(3.4)将测试集数据样本依次输入分类器，输出分类结果。
2.根据权利要求1所述的识别方法，其特征在于，步骤(1.1)中，数据来源为NinaPro开源数据集1内的10通道肌电数据。
3.根据权利要求1所述的识别方法，其特征在于，步骤(1.2)中，使用50ms长度的窗口对数据进行均值平滑滤波。
4.根据权利要求1所述的识别方法，其特征在于，步骤(2.1)中的特征集包含19种特征，它们是：信号幅值的绝对值之和(IEMG)，信号幅值的绝对平均值(MAV)，改进信号幅值的绝对平均值1(MMAV1)，改进信号幅值的绝对平均值2(MMAV2)，信号均方根(RMS)，估计肌肉收缩力的非线性检测器(v-order)，估测肌肉收缩力的对数检测器(LOG)，波形长度(WL)，幅值平均变化值(AAC)，波形长度标准差(WL-DASDV)，Willison幅值(WAMP)，自回归系数(ARC)，绝对平均值斜率(MAVSLP)，功率谱的平均功率(MNP)，Daubechies1小波变换边界(MDWT-DB1)，基于Daubechies1小波的多分辨率小波分析(MRWA-DB1)，小波包变换后信号的均值(DWPT-MEAN)，小波包变换后信号标准差(DWPT-SD)，自回归余数残差29种统计量(ARR-29)。
5.根据权利要求1所述的识别方法，其特征在于，步骤(2.2)中对训练特征集和测试特征集采用不同的特征归一化方法：
对训练特征集的归一化方法为，用训练特征集每一列减去其自身的均值并除以其方差；
对测试特征集的归一化方法为，用测试特征集样本每一列减去训练特征集对应列的均值并除以训练特征集对应列的方差；
其具体公式为：
其中train为归一化前的训练特征集，train′为归一化后的训练特征集，test为归一化前的测试特征集，test′为归一化后的测试特征集。
6.根据权利要求1所述的识别方法，其特征在于，步骤(2.3)中所述基于最小冗余最大相关性准则的自动特征选择算法，其具体实现方式为：
设存在一个特征子集S,S中每一列特征xi到类标签c的依赖性D定义为：
S中特征两两之间的冗余R定义为：
上述公式中，I为互信息，两个离散变量X和Y之间互信息的计算公式为：
最小冗余最大相关性准则的自动特征选择算法试图找到一个特征子集S，其满足拥有最大的依赖性D，最小的冗余R，即满足：
maxφ(D,R),φ(D,R)＝D-R
在手势识别时，根据“在整个分类器的训练过程中测试集信息始终保持未知”的模式识别原则，对训练特征集和测试特征集的特征选择使用不同的策略，具体步骤为：
对训练特征集的特征选择，应用基于最小冗余最大相关性准则的自动特征选择算法，选取最优的N维子集S；将子集中每一维特征对应于原始特征集中的序号保存下来；
对测试特征集的特征选择，直接应用对训练集特征选择时保存下来的序号。
7.根据权利要求1所述的识别方法，其特征在于，步骤(3.2)具体为：首先确定σ和λ的候选值，以及惩罚因子C的候选值；假设皮尔逊VII广义核参数σ有M个候选值，参数λ有N个候选值，首先产生一个M×N的网格，网格中每个点都是一个候选(σ,λ)组合；随后在网格每个节点上，通过五折交叉验证，即将训练特征集数据分成5等分，每次依次取1份作为测试集，剩余4份作为训练集，依次重复5次，每次都用该节点上的(σ,λ)值作为分类器参数，训练分类器并计算识别结果，最终得到5次测试的识别率，统计5次平均识别率，确定获得最高5折交叉验证识别率的(σopt,λopt)组合；接下来对惩罚因子C进行调优，假设惩罚因子C有K个候选值，则对每个候选值，设定之前获得最优的(σopt,λopt)组合参数，并使用5折交叉验证步骤训练分类器，得到K个候选惩罚因子C各自的5折交叉验证识别率，统计获得最高识别率的惩罚因子Copt作为最优惩罚因子；最终，输出并且保存上述调优过程获得的最优参数(σopt,λopt,Copt)。
CN201510801198.8A 2015-11-19 2015-11-19 基于支持向量机和表面肌电信号的多类手部动作识别方法 Expired - Fee Related CN105426842B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510801198.8A CN105426842B (zh) 2015-11-19 2015-11-19 基于支持向量机和表面肌电信号的多类手部动作识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510801198.8A CN105426842B (zh) 2015-11-19 2015-11-19 基于支持向量机和表面肌电信号的多类手部动作识别方法

Publications (2)

Publication Number Publication Date
CN105426842A CN105426842A (zh) 2016-03-23
CN105426842B true CN105426842B (zh) 2018-08-14

Family

ID=55505042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510801198.8A Expired - Fee Related CN105426842B (zh) 2015-11-19 2015-11-19 基于支持向量机和表面肌电信号的多类手部动作识别方法

Country Status (1)

Country Link
CN (1) CN105426842B (zh)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107037878A (zh) * 2016-12-14 2017-08-11 中国科学院沈阳自动化研究所 一种基于手势的人机交互方法
CN106778610B (zh) * 2016-12-16 2020-04-07 哈尔滨工程大学 一种基于时频图像特征的脉内调制识别方法
CN107126303A (zh) * 2017-02-15 2017-09-05 上海术理智能科技有限公司 一种基于手机app的上下肢运动辅助方法
CN106991386A (zh) * 2017-02-27 2017-07-28 杭州电子科技大学 一种基于深度残差网络的手势识别方法
CN107157450B (zh) * 2017-06-19 2020-03-31 中国科学院计算技术研究所 用于对帕金森病人的手部运动能力进行量化评估方法和系统
CN107292292A (zh) * 2017-07-20 2017-10-24 浙江好络维医疗技术有限公司 一种基于支持向量机的qrs波群分类方法
CN109412882B (zh) * 2017-08-18 2023-01-10 中兴通讯股份有限公司 特征提取方法、装置、服务器及计算机可读存储介质
CN108564105A (zh) * 2018-02-28 2018-09-21 浙江工业大学 一种针对肌电个体差异问题的在线手势识别方法
CN108875799A (zh) * 2018-05-29 2018-11-23 齐鲁工业大学 一种基于改进s变换的运动想象分类识别方法
CN109001702B (zh) * 2018-06-04 2022-05-13 桂林电子科技大学 无载波超宽带雷达人体动作识别方法
CN108985157A (zh) * 2018-06-07 2018-12-11 北京邮电大学 一种手势识别方法和装置
CN109543644B (zh) * 2018-06-28 2022-10-04 济南大学 一种多模态手势的识别方法
CN109492546B (zh) * 2018-10-24 2022-08-12 广东工业大学 一种融合小波包与互信息的生物信号特征提取方法
CN109480871A (zh) * 2018-10-30 2019-03-19 北京机械设备研究所 一种面向rsvp脑-机接口的疲劳检测方法
CN109766939B (zh) * 2018-12-29 2023-09-26 佛山科学技术学院 一种基于照片的镀锌钢和低碳钢分类方法及装置
CN109710081A (zh) * 2019-02-20 2019-05-03 杭州老板电器股份有限公司 手势感应测试装置及方法
CN109948465B (zh) * 2019-02-28 2020-11-13 浙江工业大学 一种基于多目标的表面肌电信号分类方法
CN110288048B (zh) * 2019-07-02 2022-11-29 东北大学 一种svm有向无环图的海底管道风险评估方法
CN110547806B (zh) * 2019-09-11 2022-05-31 湖北工业大学 一种基于表面肌电信号的手势动作在线识别方法及系统
CN110703915B (zh) * 2019-09-27 2022-12-27 杭州电子科技大学 一种用于同步识别手势和抓取力的粒子化方法及粒子化分析法
CN110516762B (zh) * 2019-10-10 2022-11-15 深圳大学 一种肌肉状态量化评定方法、装置、存储介质及智能终端
CN111122973A (zh) * 2019-11-26 2020-05-08 北京无线电测量研究所 一种提高自动测试覆盖率的方法、系统及电子设备
CN111209857B (zh) * 2020-01-06 2023-04-21 武汉科技大学 一种基于肌电信号的双向递归肌电特征选择方法
CN111522435A (zh) * 2020-02-21 2020-08-11 浙江工业大学 一种基于表面肌电信号的机械臂交互方法
CN111709299B (zh) * 2020-05-19 2022-04-22 哈尔滨工程大学 一种基于加权支持向量机的水声目标识别方法
CN111709314B (zh) * 2020-05-27 2023-08-22 杭州电子科技大学 基于面部表面肌电的情感分布识别方法
CN111505632B (zh) * 2020-06-08 2023-03-03 北京富奥星电子技术有限公司 基于功率谱和多普勒特征的超宽带雷达动作姿态识别方法
CN112132192A (zh) * 2020-09-07 2020-12-25 北京海益同展信息科技有限公司 一种模型训练方法、装置、电子设备及存储介质
CN112257953B (zh) * 2020-11-03 2021-09-24 上海电力大学 一种基于极地新能源发电功率预测的数据处理方法
CN113505827B (zh) * 2021-07-08 2024-01-12 西藏大学 一种机器学习分类方法
CN113730190A (zh) * 2021-09-18 2021-12-03 上海交通大学 一种三维空间运动的上肢康复机器人系统
CN113876339A (zh) * 2021-10-11 2022-01-04 浙江工业大学 一种睡眠状态脑电特征信号特征集的构建方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833671A (zh) * 2010-03-30 2010-09-15 西安理工大学 一种基于支持向量机的表面肌电信号多类别模式识别方法
CN103345641A (zh) * 2013-07-16 2013-10-09 杭州电子科技大学 基于小波包熵和支持向量机的手部肌电信号动作识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8447704B2 (en) * 2008-06-26 2013-05-21 Microsoft Corporation Recognizing gestures from forearm EMG signals

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833671A (zh) * 2010-03-30 2010-09-15 西安理工大学 一种基于支持向量机的表面肌电信号多类别模式识别方法
CN103345641A (zh) * 2013-07-16 2013-10-09 杭州电子科技大学 基于小波包熵和支持向量机的手部肌电信号动作识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Determining the Optimal Window Length for Pattern Recognition-Based Myoelectric Control: Balancing the Competing Effects of Classification Error and Controller Delay;L. H. Smith 等;《IEEE Transaction on Neural Systems and Rehabilitation Engineering》;20110430;第19卷(第2期);第186-192页 *
Facilitating the Application of Support Vector Regression by Using a Universal Pearson VII Function Based Kernel;B. Ustun等;《Chemometrics and Intelligent Laboratory Systems》;20051017;第29-40页 *
基于表面肌电信号的人手抓取动作模式识别技术研究;高允领;《中国优秀硕士学位论文全文数据库信息科技辑》;20150115;正文第16-38页 *

Also Published As

Publication number Publication date
CN105426842A (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
CN105426842B (zh) 基于支持向量机和表面肌电信号的多类手部动作识别方法
CN108491077B (zh) 一种基于多流分治卷积神经网络的表面肌电信号手势识别方法
CN110349676B (zh) 时序生理数据分类方法、装置、存储介质和处理器
CN103728551B (zh) 一种基于级联集成分类器的模拟电路故障诊断方法
CN105446484A (zh) 一种基于隐马尔科夫模型的肌电信号手势识别方法
CN106909784A (zh) 基于二维时频图像深度卷积神经网络的癫痫脑电识别方法
Chauhan et al. An efficient data mining classification approach for detecting lung cancer disease
CN109726751A (zh) 基于深度卷积神经网络识别脑电成像图的方法
KR101232945B1 (ko) 2클래스 분류 예측 모델의 작성 방법, 분류 예측 모델 작성을 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체 및 2클래스 분류 예측 모델의 작성 장치
CN109948647A (zh) 一种基于深度残差网络的心电图分类方法及系统
CN104537382A (zh) 基于遗传算法优化支持向量机的肌电信号步态识别方法
CN109497990B (zh) 一种基于典型相关分析的心电信号身份识别方法及系统
Hajian et al. New automated power quality recognition system for online/offline monitoring
CN106909946A (zh) 一种多模态融合的商品分类系统
CN101833671A (zh) 一种基于支持向量机的表面肌电信号多类别模式识别方法
CN104156690B (zh) 一种基于图像空间金字塔特征包的手势识别方法
CN110399805A (zh) 半监督学习优化svm的运动想象脑电信号分类方法
CN111476158A (zh) 一种基于pso-pca-svm的多通道生理信号体感手势识别方法
CN107491729A (zh) 基于余弦相似度激活的卷积神经网络的手写数字识别方法
Li et al. Dating ancient paintings of Mogao Grottoes using deeply learnt visual codes
CN109145971A (zh) 基于改进匹配网络模型的单样本学习方法
CN106548136A (zh) 一种无线信道场景分类方法
CN103310027B (zh) 用于地图模板匹配的规则提取方法
Banumathi et al. A novel approach for upgrading Indian education by using data mining techniques
Ramakrishnan et al. Epileptic eeg signal classification using multi-class convolutional neural network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180814

Termination date: 20201119