专利名称：基于支持向量多元分类的下肢动作模式的识别方法
技术领域：
本发明涉及一种下肢动作模式的识别方法，尤其涉及一种基于支持向量多 元分类的下肢动作模式的识别方法。
背景技术：
为控制自动机械，再现人体真实运动，作为智能机器人发展重点之一的智 能肢体动作辅助系统或交互式、示教式遥操作系统，通常的控制方法需要事先 测量或描述操作者各肢体部位及关节的动作和位置。其中应用较广泛的是利用 多种传感器如角度、速度、加速度以及力传感器等获得被描述肢体的运动学和 动力学数据。这种方法不可避免地涉及到多传感器信息融合问题，同时随传感
器种类和个数的增加，系统的复杂度及成本会迅速上升。
从功能需求分析，智能系统就是要完成如同正常人类肢体所应该完成的任 务。从人体动作的运动学描述来看， 一个复杂肢体动作可以视为是由几种常见 的简单动作片段，根据一定的规则组合而成。因此，只要可以识别出这些简单 的动作片段，就可以给智能系统提供控制指令。
当前国内外已有大量的研究，期望利用肌电信号来判断肢体运动状况，从 研究成果上来看，多集中在人体上肢运动，并取得了较好的识别效果。但是上 肢肌电信号是相对平稳而连续的高斯信号，下肢肌电信号却是典型的时变非平 稳信号，而且经常夹杂较大的下肢运动噪声，因此，如何利用肌电信号快速准 确地识别出下肢的动作模式存在一定的难度。
本发明利用小波变换方法提取肌电信号不同节律的时变信息，克服传统肌 电信号谱分析的不足，为肌电信号的动态分析和下肢动作模式的识别提供鲁棒 性好且易于识别的特征向量。
发明内容
为克服现有技术的上述缺陷，本发明基于支持向量多元分类的下肢动作模 式的识别方法提出 一种基于核聚类简化的支持向量多元分类算法，从而降低因
4数据集合不平衡性？ 1起的分类错误，并在保证分类器泛化性能的同时简化分类 器的支持向量。
为实现上述技术目的，本发明采用的技术方案如下 一种基于支持向量多元 分类的下肢动作模式的识别方法，该方法的主要步骤为首先是采集肌电信号， 接着对采集的信号进行预处理，其后进行信号的时域及频域特征提取，接着对 信号进行PCA主元分析，分析后依据足底压力信号的压力值将压力信号分为支 撑目标集与摆动目标集，随后简化SVM多元分类器对目标集数据进行分类而后 输出识别结果，其特征在于对预处理后的信号进行时域及频域特征提取时利 用一个小的移动窗口对信号进行釆样，同时将计算分散到小的时间段内进行； 肌电信号特征向量空间建立时在时域上选4奪月几电均方才艮值、绝对值平均以及方 差三个分段统计值为特征值，在频域上采用Mallat分解方法获取频域特征值。
本发明基于支持向量多元分类的下肢动作模式的识别方法与现有技术相 比，具有如下有益效果该方法通过在时频平面上表述信号的时变特征，能够 更清晰地反映出信号的频率特性随时间的变化规律，能将原始信号按照时间对 应关系分解到不同频带上或对信号进行不同的精度表达；能同时提供关于信号 时域和频域的两方面信息，从而达到有效的信号处理目的。
图1是本发明基于支持向量多元分类的下肢动作模式的识别方法的框图。 图2是本发明基于数据流分割和移动窗口概念的模式识别策略。 图3是本发明起立过程256点肌电信号的Mailt分解。
具体实施例方式
以下结合附图1至3对本发明基于支持向量多元分类的下肢动作模式的识 别方法作进一步描述。
本发明基于支持向量多元分类的下肢动作模式的识别方法首先是肌电信号 的采集，及对采集的信号进^f亍预处理。
人体在运动过程中，月几电信号的幅值和频率都随时 的改变而改变，而这 种时变性正是识别动作的有力武器，采用有效的识别策略，就可以对人体动作模式进行连续识别。传统的动作模式识别策略通常事先选定一个用于判断的时 间长度、在该时间段内信息采样完全完成后，对这采集到的数据进行特征值提 取，并进而做动作识别判断。由于该方法将特征提取和动作判断统一放到一段
时间内处理，而在采样过程中，处理器只用于控制A/D转换进行采样，因此增
加了系统对动作的反应时。本发明根据"数据流分割"和"移动窗"的概念，利用一 个小的移动窗口对信号进行采样，同时将计算分散到小的时间段内进行，提出
如图2所示的下肢动作才莫式识别方法。
其中，^为进行模式识别判断的时间长度，r为"移动窗"的长度，同时也是 窗口转换的时间。该方法不在每个采样窗内进行模式判断而只是计算该窗口信 息的特征值，在。之后再选择模式识别算法利用包含所有窗口段的特征向量进行 动作识别。 一方面在采样的同时处理器进行特征提:取计算，以充分利用处理器 的计算能力；另外一个方面，通过"移动窗"产生交叠时间段，增加了用于动作判 断的时间长度，从而提高了算法的鲁棒性。
对采集的肌电信号进行预处理后接着是对预处理的信号进行时域及频域特 征提取。人体在完成不同下肢动作任务时，相应肌群各肌肉的兴奋持续时间和 兴奋程度各不相同。因此肌电信号特征值应该包括针对振幅的特征值和针对变 化趋势的特征值。
肌电信号特征向量空间建立时在时域上选择肌电均方根值(RMS)、绝对值 平均(MAV)以及方差(VAR)三个分段统计值为特征值 一是短时统计特征 值计算量小，可以在采样的间隙完成特征提取；二是经过分段的信号统计特征 可以近似反映统计特征随时间的变化规律。在频域上从计算量、占用存储器 情况以及后续动作识别难易程度等各方面考虑，采用Mallat分解方法获取频域 特征值。
各特征量定义如下
(1) RMS: F" ' ( 1 )式中，EMG(t)为第t时刻肌肉表面电信
号的采样值(肌电信号幅值)，T为在一个采样时间段的时间长度。<formula>formula see original document page 7</formula>
式中，N为该段斗间内的采样点数，x;为第i各采样点的EMG信号幅值< 该方法只需要把EMG》文大，整波、滤波即可。
(3) VAR:用于反映信号的功率。其定义为，=jl堂(;c,-"2 (3)
)/7V-13、
为该段时间内信号的平均值。
(4) 时频特征量
利用小波Mallat算法提取肌电信号时频域特征的具体算法简介如下 设》w(W^z是"(i )中的正交小波基，则对任意/(/)ei (i ),有如下展开
<formula>formula see original document page 7</formula>
任取/ei (i )，设^为Z (i )4^的正交投影算子，^为Z (i ) —^的正交投影 算子，则有C = 2X^" 2/ = 2>,, (5)
其中<formula>formula see original document page 7</formula>
分别称为/的尺 度系数和小波系数。
由；,F,e^，我们可以得到
<formula>formula see original document page 7</formula>
<formula>formula see original document page 7</formula>
为简化分析，可假定/， p, ^都是实函数，因此可以忽略积分式中的复共
轭。由尺度函数双尺度方程可以得到
<formula>formula see original document page 7</formula>
从而有"meZ r (10)
同理可得《A = 5>(m-2。。+1,m (11) 这就是Mallat分解7>式。 在特征值提取过程中，小波基的选择和小波分解次数对运算量和分析精度 有一定影响。我们取被试1起立过程的256点数据，经过Mallat分解之后在各 个分解尺度得到的信号长度约为128, 64, 32， 16, 8......如图3所示。
从图3中可以看出，如果进行四层以上的分解，得到的点数将下降到10点
左右，影响了分析的精度，所以在实际处理过程中，Mallat分解的次数定为4次。
在小波基的选择上采用实验的方法。首先选用不同的小波基对长度为256 点的信号进行了小波分解，小波基和分解长度之间的关系如表所示
表1长度为256的信号在Mallat分解后长度与不同小波基之间的关系_
小波基 ~第一层 第二层 第三层 f^l
注Haar表示haar小波，dbN表示N阶的Daubechies小波，coifN表示N阶的coiflets小波，symN表 示N阶的Symlets小波，dmey 4<表Discrete Meyer小波。_^^^^^^^^^^^^_
从表1中可以发现小波函数越复杂，其滤器长度也就越长，得到的Mallat 分解的结果由于受到信号边界条件的影响，与理论长度之间的偏离也就越大。 所以从运算量和特征值提取方便性两个方面考虑，应该选取结构较为简单的小
8
16 32 64 128
16 32 64 128
18 34 66 129
20 36 67 130
20 36 67 130
26 41 72 133
31 46 76 136
18 34 66 129
22 38 69 131
30 45 75 135
110 120 139 178
Llf-C f CCCSSSD
o o o y y y it
c c c s 5 S Q波基。但是小波基如果太简单，就会造成在波形上的不连续，影响Mallat分解 后表现原信号的能力。经过多次试验并参考他人的研究成果，我们选取了 db2 小波作为Mallat分解的小波基。Mallat分解的结果代表了在不同带宽内信号的 强度随时间的变化情况。在每级分解中，采用奇异值分解的方法将高频系数压 缩成一个参数。这样，对于每段肌电信号可以得到4个频域特征参数。
在实际计算过程中，以256点为一个计算单位(即在采用频率为1024Hz条 件下，^取256ms)，每过64点(即r取64ms)求取一个计算单位的特征值，采 用前后两段数据192点重叠的数据采集方法，来保证预测的实时性。
这样，在时频特征量求取中，通过mallat分解得到4个特征值；在统计量 指标求取中，将256点的数据分为四个大段，每段64点，获得共4x3=12个特 征值，则每块肌肉有16个特征值，五路肌电信号组成的特征向量维数共80。由 于过多的向量维数极大增加算法复杂度，因此，在识别计算前先利用主元分析 对这个80维的空间向量降维处理，然后再进行分类。
随后通过主元分析(Principal Component Analysis: PCA)对提取获得的数 据进行处理。主元分析是通过对多元统计观测数据的协方差结构进行分析，求 出能简约表达这些数据依赖关系的主分量。即通过线性变换将原始iV维观测向 量转化为一组新特征向量，而这些新的特征向量能够反映原始向量的绝大多数 信息。因此，主元分析作为一种特征抽取的方法，也可以认为是一种数据压缩 (降维)的方法，处理方法简介如下
给定一个中心化输入向量:^ = 1,...,/)且力x,O (12) 其中，x,维数为w, —般附</。
PCA的目标就是将x(线性转化为& = "、 ( 13 )其中t/为mxw正交矩阵，第z'列t^是样本协方差矩阵的第f个特征向量。
为此，设矩阵C为 (14) 其特征值求解为义,"jC,",，卜l,…，m (15)
义,为C的一个特征值，w,是对应的特征向量。对x,.进4亍正交变化，可以得到 主元V W) = ",、V' = U，w (16)
主元具有以下性质
(1) A(/),hl,…，m互不相关；
(2) ; ,( 1"..,附依次有最大的方差；
(3) 原始输入样本的最小均方误差，第一个主元最小。
称/1,/|>7为主成份X的贡献率，为主成份力,h,...,h的累积贡献
率，表示《个主成份提取原来的特征指标xpx2，...,、信息的多少。
根据前述方法在每段肌电信号中总共提取出80个特征，在提取新特征向量 时，样本的每个统计特征指标需要用均值-方差标准化方法做预处理，以避免某 些变量占据绝对统治地位。
表2前二十个主元贡献率及累积贡献率
主元编号12345678910
各主元贡献率 (100%)30.322.613,89.838.666.10.9480.9330.860.837
累积贡献率 (100%)30.352.966.776.5385.1991.2992.2493.1794.0394.87
主元编号11121314151617181920
各主元贡献率 (100%)0.710.620. 590. 570. 5310. 4340. 3040. 1080. 1060.08
累积贡献率 (100%)95.5896.2096.7997.3697.8998.3298.6398.7498.8498.92
从表2中可以看出，前六个主元累积贡献率达到了 91.29%，可以认为，这 前六个主元在一定程度上可以代表80个指标的总体特征。在后期的分类识别过
10程中，我们将采用6个指标进^f亍下月支运动状态的识别。
接着依据足底压力信号进行动作识别，由于足底压力信息可以简单而准确 的将各摆动动作与支撑动作(含起立动作)区分开来，因此本发明中引入足底
压力信息，通过压力数据将下肢动作识别目标集分解为两个不同集合——摆动 动作集合和支撑动作集合，在此基础上利用肌电信号对每一个集合内的动作进
行进一步的细分。
本发明中基于核聚类筒化的支持向量多元分类方法具体为在建立正类和 负类集合时，首先在特征空间中对正负两类支持向量分别进行聚类，然后用聚 类后所形成的簇的质心代替蔟内的支持向量来简化SVM分类器，从而降低因数 据集合不平衡性引起的分类错误；另外，在基本保持原始分类器泛化性能的同 时简化支持向量，以提高SVM分类器的分类速度。
对于给定集合^"")，"1，…，、A"^'^，2，…，W的M类别分类问题，为 快速有效地对该集合进行分类，可采用以下步骤
(D对j'",…，M进行如下运算把第7类看做正类，把其余M-l类看做负类。
此时可得到正类集合^+，负类集合^-；
(2)简单判断正负集合的不平衡性，选定需要简化的支持向量集合并进行聚
类
假设待聚类的正(负)支持向量集为* =汰《 '"舒(^^。 = 1"..")，聚类半径 设为r ， P为将输入空间中的点映射到特征空间F的非线性映射。
二^(^i"，(^-^^i),C/Wer —m^ = l，Z = {S2"..D 。选择才羊本见eZ ,i十算P(舒禾口
第A个类的质心
之间的距离式中，p(&)， (^1，…，W为特征空间中从属于类G的支持向量。
在此基础上找到^舒距离最近的质心A,即
O =arg — 柳(;<^)
' "1 (18)
如果"("^A)^，则将"舒加入，即。=。1^(^,则类。的质心调整为
O — +^微
广^^71 ^ (19)
其中、为类q所包含的样本数目。
调整"'="'+1, 保}，并重复上述步骤至Z^0。
如果WW^^)"，则增加一个新类
C7wWer — "w附=C7ws&r — wwm +1CC7z^/^r — ww附={^ (^^Oc/w^r—训附=
(3)簇质心原像定位
特征空间中的蔟质心在形式上表示为簇内支持向量的线性组合，由于映射 ^未知，直接将此表达式代入SVM判定函数并不能达到简化SVM分类器的目 的。因此，需要簇c"^"'…，c^敏-""附的质心^在输入空间的原像^，使得 "^) = 0"要确定原像&，首先需要建立输入空间和特征空间之间的距离关系，
虽然目前只能对高斯核函数*"，力=《(1卜_,确立这种距离关系，但考虑到此类 核函数在实际应用中使用最为广泛，因此本简化法仍具有显著的实用性。
在特征空间中，任意一样本点《到质心^的距离计算如下式所示
*(0^(幼="^餘2|> *微》I
"i /vh (20)
定义向量"2=[《《'J"2/ (21)
式中，《(z = l'2'L ，"J为质心a的原像&和它的近邻^在输入空间中的距离。
对于a在特征空间中的个近邻w^^^《，，…w為r^,计算此"*个近邻在输入空间中的原像^^^^…，《be^的均值^(1、)2^1^0，并构建一个新的坐标系
H = j一_Liir
建立矩阵1 =[、， "." ]和中心矩阵—(22)
其中，/是一个 x 的单位矩阵。 则矩阵Jff/是以^为中心的中心矩阵 Z/f = [5夂一x,x^ —x，…，x、 一x] ( 23 )
假设矩阵I//的秩为《，对其进行奇异值分解(Singular Value Decomposition:
SVD)，有
A, 0. 0 0
=五A^f =五!r
(24)
其中S 为由 一 组标准正交列向量e,组成的"x《矩阵；
r二A,《屮,，C2,…^]为一《x 矩阵；列向量c,为向量^-^在A上的投影。此时
i|C "||2,/ = 1"..，""
定义一个 xi的向量"卜[IM2，lk2112，…，11、112『。则经过变换后可以得到质心"
在输入空间中的原^象的近似值[i]: ^ =PlAl—^(《—"2)" (25 )
(4) 简化向量的最佳权值计算
为使IT:,"W(幼和它的近似值A^")尽量相等。定义
"(a;hiap(z4)—iT:iA,(幼i卩，A^,…,c'w"-"""7 (26)
式中，^为类&的质心的原像，A(''=1'2，L ，"J为从属于类G的支持向量，％ 为对应的斥又值。
对(26)式关于A求导以获得使"(A)取最小值的最佳权值A。即令 ▽AW(A)) = o,可得A=g ^",V"",4), hi，…,c/附妙—"， (27)
对于高斯核函数而言，"z"z4) = i,此时，A^2^"J(^，x。。
(5) 简化向量后，对支持向量集合进行分类
13此时，相应的判定函数为户W-sgn(^(x)) (28)
其中 ^ (29)
式中，^为待分类的向量，^为支持向量的数量，6为偏置。"，乂) (''=1，2，L ，Ws )就是所谓的支持向量，它们对应的拉格朗日乘子"'大于零。
从(28)及(29)式可以看出判定一个未知类别的样本所需要的时间和支持向 量的数目成正比。因此，削减支持向量的数量能够有效地提高分类机的分类速 度。
(6)通过上述步骤，求出每一个决策函数户")，判定x的从属类别J。 J是 g'(x), ^W中最大者的上标。
在应用上述算法进行分类求解时，必须注意下述问题
向量削减过程中，通常会设定差异阈值对简化前后分类超平面变化的最大 值作出限制，其目的是限制分类器的泛化性能损失。在利用迭代法增加聚类半 径的大小来尽可能地简化支持向量过程中，每次迭代后，计算所形成的超平面 和原始分类超平面之间的差异如果二者之间的差异超过差异阈值、则迭代简 化过程终止，取上一次迭代所得的简化SVM分类器作为最终的简化SVM分类 器；否则，增大聚类半径以进一步简化支持向量。
在利用支持向量进行分类训练时，核函数选取为高斯核函数。因此需要预 先确定学习参数C、 c7,本文采用交叉校验(Cross Validation)的方法来确定这 两个参数。将训练集分为IO等份，每次留一份验证模型性能，其余90°/。用来训 练模型，最后以模型在20次验证数据上的性能平均值作为这一学习参数下的模 型性能。
将参数C、 c7作为变量，以分类错误率(ErrorRate: ER)为指标衡量SVM分类模型性能，用梯度下降优化方法来极小化SVM分类模型性能的具体算法如
下
(1)设置初始的参数值(C、 c )和参数改变步进量AC、 Act,置迭代次 数/,ew = 1 ，性能没有改善的连续迭代次数M ctowge = 0 ; (2 )计算当前参数值交叉校验的平均性能^ (C,c7);
(3) 将参数值分别更改为(C-AC,c7)和(C + AC,(t),计算相应的平均性能 ￡/ (C - AC, o")和五及(C + AC， o");若￡7 (C - AC, cr)和五i (C + AC， ct)均大于￡i (C, o")，则 ^己= 0 ,若￡i (C — AC, cr) <五i (C + AC, cr) ， i己A￡i c =— AC, o")—五i (C, cr),反 之，记A￡i r = ￡i (C， cr)—五i (C + AC, (t);
(4) 与(3)类似，将参数值分别更改为(C,o"-Ao")和(C,(7 + Ao"),计算扁，; (5 )令M7 = max(flfe(AEi c, A￡i^)),若A￡i = 0 ,令iVoctewge = 7Vocto"ge +1 ，同
时令AC二AC/2, Acr = Ao"/2;否则令C = C +AC*A￡i c/A￡i , o" = ct + Act*AEi^/AEJ , iVoc/z朋ge = 0 ;
(6 ) /fera = +1 ，若/fera大于^见定的最大迭^RP欠数或iVoc/2fl"ge大于允许的 最大无性能改善迭代次数，则退出；否则返回(2)。
从被试1的数据中，每一个动作数据等间隔地取1/2作为测试样本，另外1/2 作为训练样本。为更好地对参数优化前后SVM多元分类方法性能进行分析，同 时与神经网络(NeuralNetworks: NN)中的BP算法进行对照，隐层个数为10。 下肢动作识别的SVM参数实验结果如表3所示
表3支持向量机参数优化前后性能比较
预测方法SVM (优化前)SVM (优化后)BP
训练集 测试集错误率(％) 时间(s) 错误率(％) 时间(s)12.7 0.09 24.2 0.089.5 0.13 12.6 0.1222.4 0.06 22.8 0.07从表3中可以看出，经交叉校验和梯度下降法优化SVM性能后，系统的 错误率性能改善了约48%,其中优化前参数为C = 10, o" = 0.05，经过上述迭代 过程可得到最优的参数为C-156、 CT = 0.5。从表也可以看出，优化前SVM性能 在分类准确性上略差于BP，而优化后分类准确性优于BP。 ^f旦同时我们也可以 观察到，优化后SVM所需时间要长于BP和未优化前的SVM。
在确定SVM的参数C456、 ^7 = 0.5后，为改善SVM在多元分类的实时性， 加快分类速度，我们用核聚类简化的SVM方法对数据集进行测试，考察不同差 异阈值r下，算法的分类性能，测试数据同上，实验结果如表4所示
表4不同差异阈值下的分类器性能
训练
BP
原始 SVM
差
异
阈 值
0.15
0.25
0.3
0.35
0.45
训练集错误率(％)12.579.5810.9
时间(s)0.130.110.12
测试集错误率(％) 时间(s)13.22 0.0813.2 0.1212.4 0.07
训练集错误率(％)9.579.59.54
时间(s)0.210.200.20
测试集错误率(％)10,01.9.69.9
时间(s)0.180.180.18
训练集错误率(％)9.569.59.55
时间(s)0.220.210.22
测试集错误率(％)9.939.779.91
时间(s)0.200.180.19
训练集错误率(％)9.599.399.56
时间(s)0.200.190.20
测试集错误率(％)10,229.69.88
时间(s)0.160.170.17
训练集错误率(％)9.579.512.4
时间(s)0.160.150.16
测试集错误率(％)11.2211.611.4
时间(s)0.160.170.16
训练集错误率(％)10.579.9810.13
时间(s)0.160.160.16
测试集错误率(％) 时间(s)12.22 0.149.89 0.1611.04 0.15
训练集错误率(％) 时间(s)11.57 0.1310.5 0.1411.4 0.14
测试集错误率(％)13.229.610.4
时间(s)0.120.12'0.12
从表4中可以看出，随着差异阈值取值增大，简化法对支持向量的简化率升 高，相应简化后的SVM分类器在测试集上的分类速度也随着提高。虽然筒化SVM分类器的泛化性能随着差异阈值增大而有所降低，但同简化率相比，泛化 性能的损失微乎其微。当差异阈值取0.45时，简化SVM分类器在测试集上的分 类速度是原SVM分类器的2倍，却仅有7%的分类精度损失，这说明本章所提 出的简化方法在极大削减支持向量的同时，基本保持了原SVM分类器的分类精 度，可以较大地提高SVM分类器的分类效率，解决SVM应用于实时分类所存 在的速度瓶颈问题。
简化SVM多元分类器对目标集进行分类后将动作识别结果输出，动作识别 结束，动作识别没有结束则重新进入肌电信号采集步骤， 一次新的动作识别循 环开始。
权利要求
1、一种基于支持向量多元分类的下肢动作模式的识别方法，该方法的主要步骤为首先是采集肌电信号，接着对采集的信号进行预处理，其后进行信号的时域及频域特征提取，接着对信号进行PCA主元分析，分析后依据足底压力信号的压力值将压力信号分为支撑目标集与摆动目标集，随后简化SVM多元分类器对目标集数据进行分类而后输出识别结果，其特征在于对预处理后的信号进行时域及频域特征提取时利用一个小的移动窗口对信号进行采样，同时将计算分散到小的时间段内进行；肌电信号特征向量空间建立时在时域上选择肌电均方根值、绝对值平均以及方差三个分段统计值为特征值，在频域上采用Mallat分解方法获取频域特征值。
2、 如权利要求1所述的一种基于支持向量多元分类的下肢动作模式的识别 方法，其特征在于主元分析处理方法为给定一个中心化输入向量x々=l,...，/) 且1>, =0其中，A维数为附， 一般附</， PCA将^线性转化为 其中"为腦m正交矩阵，第/列",是样本协方差矩阵的第/个特征向量，设矩阵C为 z^ ''其特征值求解为>,=C,",，/ = l,...，m， A为C的一个特征值，"'是对应的特征向量，对、进行正交变化，得到主元^0-"'、，"1，…^。
3、 如权利要求1所述的一种基于支持向量多元分类的下肢动作模式的识别 方法，其特征在于简化SVM多元分类器在建立正类和负类集合时，首先在特 征空间中对正负两类支持向量分别进行聚类，然后用聚类后所形成的簇的质心 代替簇内的支持向量来简化SVM分类器。
4、 如权利要求3所述的一种基于支持向量多元分类的下肢动作模式的识别 方法，其特征在于SVM分类器对于给定集合的分类采用以下步骤(1)集合5 = (^乂),/ = 1"."/， x,e/^,y,e(l,2,…,ika对y-l，…,M进行如下运算把第y类看做正类，把其余M-1类看做负类，此时可得到正类集合^+,负类集合(2) 简单判断正负集合的不平衡性，选定需要简化的支持向量集合并进行聚类；(3) 簇质心原像定位；(4) 筒化向量的最佳权值计算；(5) 简化向量后，对支持向量集合进行分类；相应的判定函数为<formula>formula see original document page 3</formula>其中<formula>formula see original document page 3</formula>式中，x为待分类的向量，^为支持向量的数量，6为偏置。"，X) 。 = 1，2，L ，^ )就是所谓的支持向量，它们对 应的拉格朗日乘子"'大于零(6) 通过上述步骤，求出每一个决策函数尸")，判定x的从属类别J, J是 g1"), g"")中最大者的上标。
5、如权利要求4所述的一种基于支持向量多元分类的下肢动作模式的识别 方法，其特征在于Mallat分解的次数定为4次。
全文摘要
一种基于支持向量多元分类的下肢动作模式的识别方法，该方法的主要步骤为首先是采集肌电信号，接着对采集的信号进行预处理，其后进行信号的时域及频域特征提取，提取时利用一个小的移动窗口对信号进行采样，同时将计算分散到小的时间段内进行；肌电信号特征向量空间建立时在时域上选择肌电均方根值、绝对值平均以及方差三个分段统计值为特征值，在频域上采用Mallat分解方法获取频域特征值，接着对信号进行PCA主元分析，随后依据足底压力信号的压力值将压力信号分为支撑目标集与摆动目标集，随后简化SVM多元分类器对目标集数据进行分类而后输出识别结果。本方法克服传统肌电信号谱分析的不足，为下肢动作模式的识别提供鲁棒性好且易于识别的特征向量。
文档编号G06K9/00GK101587546SQ200910099598
公开日2009年11月25日 申请日期2009年6月12日 优先权日2009年6月12日
发明者群 吴, 吴剑锋, 孙守迁, 曾志强 申请人:浙江大学