一种基于肌电信号的双向递归肌电特征选择方法技术

本发明专利技术是一种基于肌电信号的双向递归肌电特征选择方法，所述方法包括基于最大相关的前向动态递归特征选择和基于最小冗余的后向动态过滤特征选择。对于前向选择，首先构建原始肌电特征集，并通过使用互信息准则去除对手势识别贡献最低的特征产生肌电特征总集；其次，结合本文所提特征搜索策略和特定机器学习模型对肌电特征总集进行相关性分析，得到前向最佳候选肌电特征子集。对于后向选择，从原始肌电特征集开始，通过本文所提评估方式，每次迭代去除一个最不相关的肌电特征，分析特征集的冗余性，直至得到后向最佳候选肌电特征子集。最后取前向与后向最佳候选肌电特征子集的交集，并进行结果验证，输出最佳肌电特征集。

A bi-directional recursive feature selection method based on EMG signal

【技术实现步骤摘要】
一种基于肌电信号的双向递归肌电特征选择方法
本专利技术涉及信号处理与人机交互领域，特别涉及一种基于肌电信号的双向递归肌电特征选择方法。
技术介绍
随着科技的发展，人机交互的研究与应用受到越来越多的学者关注。其中，手势识别作为人机交互的一个主要分支，是人机交互领域的重要研究课题之一。人机交互中的手势对使用者来说具有便捷、交互性强、表达内容丰富等优点，一直是人机交互领域的首要选择。肌电信号是在人体肌肉活动过程中产生的一种生物电信号，是不同的运动单位动作电势在肌肉纤维和细胞组织中传导叠加的综合结果。肌电信号主要是通过提取肌电特征进行手势识别，选择合适的肌电特征进行组合可以获得更高的手势识别率。鉴于此需要一种肌电特征选择方法，通过意大利经济学家帕累托的关键少数法则，利用20％的最佳肌电特征代替余下80％的肌电特征所做的贡献，在取得最佳效果的同时减少资源损耗，并提高识别效率。
技术实现思路
本专利技术提供了一种基于肌电信号的双向递归肌电特征选择方法，可以从多种肌电特征中，通过前向动态递归特征选择(FDRFS)和后向动态过滤特征选择(BDFFS)分别得到前向和后向最佳肌电特征组合，然后将它们取交集得到最佳肌电特征组合。为解决上述技术问题，根据
技术实现思路
，采用如下技术方案为一种基于肌电信号的双向递归肌电特征选择方法，具体包括以下步骤：步骤1：通过n个肌电特征向量构建肌电特征集；步骤2：计算每个肌电特征向量的方差以及互信息，剔除肌电特征集中无贡献的机电特征向量或者评分最低的前m1(m1<0.1n)个肌电特征向量，得到第一次迭代肌电特征集合；步骤3：从过滤后肌电特征总集中随机挑选一定数量的肌电特征向量构建随机肌电特征集，计算随机肌电特征集的组合数和识别率；步骤4：剔除候选肌电特征集中相对冗余特征；步骤5：更新动态递归系数，确定候选肌电特征集；步骤6：重复步骤3–步骤5，直至动态递归系数达到最大迭代次数，输出前向最佳肌电特征组合；步骤7：输入n个肌电特征构建的肌电特征集，确定反向肌电特征组合；步骤8：计算肌电特征子集的识别率，确定后向候选肌电特征集；步骤9：重复步骤8，直至达到迭代条件，输出后向最佳肌电特征组合；步骤10：通过取前向最佳肌电特征组合和后向最佳肌电特征组合的交集得到最佳肌电特征集。作为优选，步骤1中所述初始肌电特征集为迭代次数为h∈[0,stepmax]，b为输入肌电特征向量个数，b∈[0,n]，为初始肌电特征集；所述初始肌电特征集为：其中，Xi，i∈[1,n]为第i个肌电特征向量，为l维度向量，具体定义为：其中，Xi,j，i∈[1,n]，j∈[1,l]为第i个肌电特征向量中第j个肌电特征。作为优选，步骤2中所述每个肌电特征向量的方差为：其中，Xi表示第i组特征向量，l表示特征向量中的样本数，Xi,j表示特征向量中第i个样本值，μ(Xi)表示该特征向量的平均值；步骤2中所述计算每个原始肌电特征的互信息为：其中，Y表示手势的标签向量，p(xi,j,y)表示Xi与Y两个变量之间的联合概率密度分布函数，p(xi)和p(y)分别表示Xi和Y两者的边缘概率分布；步骤2中所述剔除肌电特征集中无贡献或贡献小的肌电特征向量，具体为：Var(X)＜α或I(X；Y)＜β认为是贡献小的肌电特征向量，其中α,β分别是方差阈值和互信息阈值；从步骤1的特征集合中筛选出贡献小的肌电特征向量，然后剔除贡献小较低的前m1(m1<0.1n)个肌电特征向量；步骤2中所述第一次迭代肌电特征集合为：T1k＝{X1,X2,…,Xk,k<n}作为优选，步骤3中所述随机挑选一定数量的过滤后肌电特征向量构建随机肌电特征集为：从T1k中任意挑选组(p1·k)种组合肌电特征特征向量，其中，p1为第一动态递归系数，q为降采样系数；利用支持向量机(SVM)f(·)计算所有特征向量组合的识别率具体计算过程为：将训练集输入到SVM中，得到决策边界：其中，为拉格朗日算子对偶最优解，K(Xi,X)为高斯核函数，b*为最优截距，yi∈Y；再将测试集输入到SVM中，可计算出为前向第一次迭代中第一个特征向量组合的识别率；选择识别率最高特征向量组合作为该批次下的次优特征向量组合，并放入候选肌电特征集中：作为优选，步骤4中所述剔除候选肌电特征集中相对冗余特征，具体过程为：将步骤3得到的S1中挑选p1·p2·kl种肌电特征向量剔除，使得余下特征组合的f(·)相对最优SAcc1；其中，p2为第二动态递归系数；从T1k减去过滤后的S1作为第二次迭代的肌电特征集合供下一步选择；作为优选，步骤5中所述更新动态递归系数为：p1＝p1-Δ，p2保持不变，Δ为系数步长；步骤5中所述确定候选肌电特征集为：从步骤4中得到第二次迭代特征集合中按照步骤3中方法挑选特征的组合，每种组合均与步骤4中所述候选特征子集S1取并集，并利用f(·)计算此时所有特征组合的将最高的组合作为最优特征组合放入候选肌电特征集S2中；作为优选，步骤6中所述动态递归系数达到最大迭代次数为：p1·k＜1达到最大迭代次数即Stepmax时；步骤6中所述前向最佳肌电特征组合为作为优先，步骤7所述输入n个特征组合而成的肌电特征集为RTsn，确定反向肌电特征组合，为后向初始肌电特征集；步骤7所述确定肌电特征组合，具体为：所述肌电特征集进行每一种特征的依次剔除，余下所有特征形成肌电特征子集，具体如下：作为优选，通过f(·)计算步骤8中所有后向肌电特征子集的识别率为s为迭代次数，具体计算过程如下：将训练集输入到SVM中，得到决策边界：其中，为拉格朗日算子对偶最优解，K(Xi,X)为高斯核函数，b*为最优截距，yi∈Y；为后向第一次迭代中第一个肌电特征组合的识别率；利用f(·)计算每一种肌电特征子集的识别率选出识别率最高MRAcc1的特征组合作为本轮候选肌电特征子集RS1，此时认为特征Xi相对余下所有特征是最不相关的，将其从RT1n中剔除得到供下一步选择；作为优选，步骤9中所述迭代条件为：或者时停止，η为识别率阈值，为第s次迭代后的特征数；后向最佳肌电特征组合为本专利技术优点在于，可以提供一种更好的特征选择方法，通过前向动态递归和反向动态过滤两个方向对肌电总集进行特征选择，将两个方向各自得到的最优特征组合取交集，得到最终的最优特征组合，该特征组合满足最大相关和最小冗余的要求。附图说明图1：为本专利技术方法的流程图；具体实施方式下面结合附图对本专利技术的较佳实施例进行详细阐述，以使本专利技术的优点和特征能更易于被本领域技术人员理解，从而对本专利技术的保护范围做出更为清楚本文档来自技高网...

【技术保护点】
1.一种基于肌电信号的双向递归肌电特征选择方法，其特征在于，包括以下步骤：/n步骤1：通过n个肌电特征向量构建肌电特征集；/n步骤2：计算每个肌电特征向量的方差以及互信息，剔除肌电特征集中无贡献的机电特征向量或者评分最低的前m

【技术特征摘要】
1.一种基于肌电信号的双向递归肌电特征选择方法，其特征在于，包括以下步骤：
步骤1：通过n个肌电特征向量构建肌电特征集；
步骤2：计算每个肌电特征向量的方差以及互信息，剔除肌电特征集中无贡献的机电特征向量或者评分最低的前m1(m1<0.1n)个肌电特征向量，得到第一次迭代肌电特征集合；
步骤3：从过滤后肌电特征总集中随机挑选一定数量的肌电特征向量构建随机肌电特征集，计算随机肌电特征集的组合数和识别率；
步骤4：剔除候选肌电特征集中相对冗余特征；
步骤5：更新动态递归系数，确定候选肌电特征集；
步骤6：重复步骤3–步骤5，直至动态递归系数达到最大迭代次数，输出前向最佳肌电特征组合；
步骤7：输入n个肌电特征构建的肌电特征集，确定反向肌电特征组合；
步骤8：计算肌电特征子集的识别率，确定后向候选肌电特征集；
步骤9：重复步骤8，直至达到迭代条件，输出后向最佳肌电特征组合；
步骤10：通过取前向最佳肌电特征组合和后向最佳肌电特征组合的交集得到最佳肌电特征集。


2.根据权利要求1所述的基于肌电信号的双向递归肌电特征选择方法，其特征在于：步骤1中所述初始肌电特征集为迭代次数为h∈[0,stepmax]，b为输入肌电特征向量个数，b∈[0,n]，为初始肌电特征集；
所述初始肌电特征集为：



其中，Xi，i∈[1,n]为第i个肌电特征向量，为l维度向量，具体定义为：



其中，Xi,j，i∈[1,n]，j∈[1,l]为第i个肌电特征向量中第j个肌电特征。


3.根据权利要求1所述的基于肌电信号的双向递归肌电特征选择方法，其特征在于：步骤2中所述每个肌电特征向量的方差为：



其中，Xi表示第i组离散目标特征向量，l表示特征向量中的样本数，Xi,j表示特征向量中第i个样本值，μ(Xi)表示该特征向量的平均值；
步骤2中所述计算每个原始肌电特征的互信息为：



其中，Y表示手势的标签向量，p(xi,j,y)表示Xi与Y两个变量之间的联合概率密度分布函数，p(xi)和p(y)分别表示Xi和Y两者的边缘概率分布；
步骤2中所述剔除肌电特征集中无贡献或贡献小的肌电特征向量，具体为：
Var(X)＜α或I(X；Y)＜β认为是贡献小的肌电特征向量，其中α,β分别是方差阈值和互信息阈值；
从步骤1的特征集合中筛选出贡献小的肌电特征向量，然后剔除贡献小较低的前m1(m1<0.1n)个肌电特征向量；
步骤2中所述第一次迭代肌电特征集合为：
T1k＝{X1,X2,…,Xk,k<n}。


4.根据权利要求1所述的基于肌电信号的双向递归肌电特征选择方法，其特征在于：步骤3中所述随机挑选一定数量的过滤后肌电特征向量构建随机肌电特征集为：



从T1k中任意挑选组(p1·k)种组合肌电特征特征向量，其中，p1为第一动态递归系数，q为降采...

【专利技术属性】
技术研发人员：程扬威，李公法，蒋国璋，陶波，余明超，徐超，
申请(专利权)人：武汉科技大学，
类型：发明
国别省市：湖北;42