一种基于机器学习的触觉信息分类计算建模方法
【专利摘要】一种基于机器学习的触觉信息分类计算建模方法，获取训练集样本的触觉序列，使用线性动态系统模型建模提取子触觉序列的动态特征，使用马丁距离计算子触觉序列的动态特征之间的距离，使用K?中心点算法对马丁矩阵进行聚类构建码书，使用码书对每组触觉序列进行表征，得到系统包模型，将训练集样本的系统包模型和训练集样本标签一起送入极限学习机中训练分类器，将待分类样本的系统包模型送入分类器得到物体类型的标签。本发明解决了机器人对非合作目标稳定、柔顺抓取的实际需求，为精细操作任务的完成提供数据基础，并可与其它传感结果融合计算，从而通过多源深度感知增强对不同目标物的描述和辨识能力，为智能操控的实施奠定技术基础。
【专利说明】
-种基于机器学习的触觉信息分类计算建模方法
技术领域
[0001] 本发明设及机器人触觉建模领域，尤其设及一种基于机器学习的触觉信息分类计 算建模方法。
【背景技术】
[0002] 目前机器人为实现精细操作通常配有多种传感器，但各类传感器一般只对不同模 态传感器采用独立的应用方式认识周围的环境，割断了信息之间的内在联系，因此会严重 降低感知动作的智能化程度。为了能够准确提供操作装置本身的状态、操作对象的位置、属 性等信息，需要研究视觉、触觉等多模态融合的理论与方法，从不同的角度对操作对象的材 质、形变、位置、距离进行测量，从而为精细操作任务的完成提供数据基础。而视触觉的模型 建立首先需要建立触觉模型，然而，经专利检索，目前尚无对触觉信息进行分类计算的建模 方法。为满足机器人对非合作目标稳定、柔顺抓取的实际需求，迫切需要研究复杂作业环境 下视触觉信息的感知、高效表示和融合机理，其中，触觉模型的建立尤为重要。
【发明内容】
[0003] 本发明提供一种基于机器学习的触觉信息分类计算建模方法，解决了机器人对非 合作目标稳定、柔顺抓取的实际需求，为精细操作任务的完成提供数据基础，并可与其它传 感结果融合计算，从而通过多源深度感知增强对不同目标物的描述和辨识能力，为智能操 控的实施奠定技术基础。
[0004] 为了达到上述目的，本发明提供一种基于机器学习的触觉信息分类计算建模方 法，包含W下步骤：
[0005] 步骤Sl、利用触觉传感器获取训练集样本的触觉序列；
[0006] 步骤S2、将每组触觉序列合并为子触觉序列，使用线性动态系统模型对每组子触 觉序列进行建模，分别提取子触觉序列的动态特征；
[0007] 步骤S3、使用马下距离作为量度计算子触觉序列的动态特征之间的距离，获得马 下矩阵；
[000引步骤S4、使用K-中屯、点算法对马下矩阵进行聚类，构建码书；
[0009] 步骤S5、使用码书对每组触觉序列进行表征，得到外部特征为直方图的系统包模 型；
[0010] 步骤S6、将训练集样本的系统包模型和训练集样本标签一起送入极限学习机中训 练分类器；
[0011] 步骤S7、利用触觉传感器获取待分类样本的触觉序列，对待分类样本的触觉序列 按照步骤S2~步骤S5的方法进行处理，获得待分类样本的系统包模型；
[0012] 步骤S8、将待分类样本的系统包模型送入步骤S6获得的分类器中，即可得到物体 类型的标签。
[0013] 所述的步骤Sl中，所述的触觉阵列传感器设置在机器人的机械臂末端的灵巧手 上，所述的灵巧手具有=个手指，每个手指末端分别配置一个触觉阵列传感器，每个触觉阵 列传感器具有mXn个传感器单元，所述的触觉序列为通过触觉阵列传感器获得的给定N个 训练时间序列样本SiGRdXT,其中，d是空间维度，T是序列的时间长度，i = l，2,…，N。
[0014]所述的步骤S2中，首先通过降低维数对触觉序列的主要成分进行提取，去除冗杂 的成份。
[001引所述的步骤S2中，
[0016] 通过重塑m X n矩阵作为一个m X n维向量，分别获得灵巧手上=个手指的(m X n)维 动态触觉序列S(l> G RfmXnW，s(2) E RfmXrOXT, s(3) E RUXrOXT ;
[0017] 虫聪3个丰指的m Xn维触觉序列，合并为一个单一的子触觉序列：
[001 引
（1)
[0019] 其中，d=mXnX3;
[0020] 假设噪声为w(t)GRdxT，将子触觉序列S分解为:S(t)=D(t)+w(t)，其中D(t)GRd XT为不含噪声部分；
[0021] 假设存在一个P维滤波器0。:巧一於/XT, a= 1.….化伸得前)ER, 0为隐 层特征空间的维数，有公的)，使得
，其中，
[0022] Ai和Ci为满足上式的特征参数，t是S的长度自变量，t = 1，2，…，d X T，并且v( t) G
RdXT独立戸A如羊口品击姑E盛密賠。、珪蜡结果女日下：
[0023] (2)
[0024] 令Mi = (Ai，Ci)，即为子触觉序列的动态特征，即某一频率范围内的触觉序列。
[00巧]所述的步骤S3中，使用马下距离来进行两个LDS的特征化=(心，打），12=(42，〔2)之 间的距离衡量，得到马下距离矩阵，md表征输入样本的触觉子序列个数；
[0026] 伽常:库列的IDS据佈，间的細则值由下击讲斤宙公，
[0027]
(3)
[0028] 对于任意两个模型之间规则角的计算，使用李雅普诺夫方程求解，保证系统稳定 的重要条件是存在正定的对称矩阵P满足下式：
[0029] A Tdh-D - _rV { a\
[003C
[0031
[003;
[0034] (6)
[003引景终得到两个LDS据佈Ml与]间的马T距离dM(Mi .Ms):
[00巧] 骄寺征间的马下距离：
[0036] 讯。
[0037] 所述的步骤S4中，使用K-中屯、点算法对马下矩阵进行聚类的步骤包含：
[0038] 步骤S4.1、随机选取k个马下距离矩阵中的数据作为聚类中屯、点；
[0039] 步骤S4.2、将矩阵中剩余的数据按照所选取的聚类中屯、点进行分组，分组W每组 中每个数据点与每组中聚类中屯、点的距离数据值的平方和最小为原则；
[0040] 步骤S4.3、重复执行步骤S4.1、步骤S4.2和步骤S4.3，更新每组数据样本的聚类中 屯、点；
[0041] 步骤S4.4、直至聚类中屯、点稳定不变，完成数据聚类；
[0042] 聚类完成后得到由k组LDS特征组成的码书(Codebook) 为表示LDS 特征的元组，对于码书中的任意一组元组(A，C)，称其为码词(Codewords)。
[0043] 所述的步骤S5中，使用码书对触觉序列进行表征后可W得到系统包模型 {BoS巧;，m为样本中触觉序列的个数，运种表征方式的外部体现为直方图模型 /1 =的，f G W"'，可由特征词频率算法得到；
[0044] 假设在第i组触觉序列中，第j组码词出现的次数为CU次，则有下式：
[0045]
(8)
[0046] 式中，hij表示在第i组触觉序列中，第j组码词出现的频率，m为样本中触觉序列个 数，对于hi则有hi=比ir，，hik]，为一组特征向量。
[0047] 所述的步骤S6中，
[004引假定存在一组训练触觉序列则该触觉序列可W使用具有成个隐层结 点且扩展到RBF核函数的单隐层核函数网络进行建模，函数模型如下式所示：
[0049]
(9).
[0050] 式中，是连接第i个核函数与输出神经网络的输出权重向量，(61 (rj)为输出函数，r=[rr-rm]T为相应的特征向量，化为函数中屯、点，Oi为影响因子；
[0051] 高斯核函数可W被定义为：
[005^
(10)
[0053]使极限学习机具有多个输出结点可W完成对数据的多分类，也就是使进行M种类 别识别的分类器具有M组隐含层结点，极限学习机算法对于属于该类别的样本使用+1进行 标记，对于不属于该类的样本使用0进行标记，极限学习机算法在使训练样本误差最小的同
[0化4] 时，也使样本的输出权重最小，即最小化Mhco-tM2与M CO M的值，式中H为神经网络隐含 层输出矩阵，它可由下式进行表示：
(11)
[005引 （12)
[0化6] )
[0057]式中，矩阵巧f为神经网络隐含层输出矩阵H的广义逆矩阵，在训练具有多输出结 点的极限学习机分类器时遵循如下条件：
[0化引 （14)
[0059] U5)
[0060] 上式中，C亦为损耗参数，可由交叉验证获得点=[Cm,…点,k]T为训练集误差向 量，此时，训练极限学习机分类器等效于解决下式的最优化问题： 幽]
(16)
[0062] 式中，为训练样本相应的拉格朗日系数；[0063] 当训练集样本规横不大时，可W得到下式进行描述：
[0064； (17)
[00化：
[0066： (18)
[0067： SJW表示为
[0068； (19)
[0069] 使用Qelm = HHT对极限学习机的核函数矩阵进行定义，并且满足W下核函数公式 Q ELMa,b= 4 (raH (rb)=K(ra，rb)，使用核函数的极限学习机分类器输出函数可W按照下 式进行?^?公
[0070]
(20)
[0071] 通过上述算法可完成对具有多输出结点的极限学习机分类器的训练，并且可W对 测试集样本的物理标签进行预测，所预测的标签即为输出的结点所对应的索引，此索引对 应上式中最大的输出值，且有ltest = arg maxfjjGU，.. .m}，式中fj是第j个输出结点的输 出函数。
[0072] 所述的步骤S8中，样本对应的标签为Ii，其中i = 1，2，…，N，标签是物体的种类信 息，可用数字、文字或字母表示。
[0073] 本发明解决了机器人对非合作目标稳定、柔顺抓取的实际需求，为精细操作任务 的完成提供数据基础，并可与其它传感结果融合计算，从而通过多源深度感知增强对不同 目标物的描述和辨识能力，为智能操控的实施奠定技术基础。
【附图说明】
[0074] 图1是本发明提供的一种基于机器学习的触觉信息分类计算建模方法的流程图。
[0075] 图2是多指灵巧手的示意图。
[0076] 图3是触觉阵列传感器的示意图。
[0077] 图4是马下距离矩阵定量示意图。
[0078] 图5是系统包直方图模型外观与内部组成图。
[0079] 图6是分类准确率曲线。
【具体实施方式】
[0080] W下根据图1-图6,具体说明本发明的较佳实施例。
[0081] 如图1所示，本发明提供一种基于机器学习的触觉信息分类计算建模方法，包含W 下步骤：
[0082] 步骤Sl、利用触觉阵列传感器获取训练集样本的触觉序列；
[0083] 本实施例中，所述的触觉阵列传感器设置在机器人的机械臂末端的灵巧手上，如 图2所示，所述的灵巧手具有=个手指，分别标记成手指1，手指2,手指3,每个手指末端分别 配置一个触觉阵列传感器，如图3所示，每个触觉阵列传感器具有mXn个传感器单元，在每 个时间段，从每个手指上的触觉阵列传感器的mXn个传感器单元中收集触觉数据；
[0084] 所述的触觉序列为通过触觉阵列传感器获得的给定N个训练时间序列样本SiGRd XT，其中，d是空间维度，T是序列的时间长度，i = l，2,…，N;
[0085] 步骤S2、将每组触觉序列合并为子触觉序列，使用线性动态系统模型化inear Dynamic Systems，LDS)对每组子触觉序列进行建模，分别提取子触觉序列的动态特征；
[0086] 由于数据集中触觉序列的维数较高，会使得后续计算时间相对较长，而大部分触 觉序列中都包含了较多杂乱和冗余的信息；因此，在运里通过降低维数对触觉序列的主要 成分进行提取，去除相对冗杂的成份；
[0087] 维数的降低可W通过从高维到低维的函数映射来实现，并且在降低维数的过程 中，依据函数映射的存在与否，划分为线性维数降低算法和非线性维数降低算法两类;线性 维数降低算法为特征选择算法，它是通过选择实验样本的整体特征向量中的一组作为特征 来描述该实验样本;非线性维数降低算法为特征提取算法，是将已经存在的特征向量组合 成新的样本特征向量来进行表示;为进一步融入触觉序列的时空信息，利用LDS来对触觉序 列进行建模；
[0088] 通过重塑mXn矩阵作为一个(mXn)维向量，分别获得手指1、手指2、手指3的(mX n)维动态触觉序列 S(l)ER(mXn)XT，s(2)ER(mXn)XT，s(3)ER(mXn)XT;
[0089] 串联3个手指的(mXn)维触觉序列，合并为一个单一的子触觉序列：
[0090]
(1)
[0091] 其中，d=mXnX3;
[0092] 假设噪声为w(t)GRdxT，将子触觉序列S分解为:S(t)=D(t)+w(t)，其中D(t)GRd XT为不含噪声部分；
[0093] 假设存在一个P维滤波器0 ，〇 = 1，…，Q，使得x(t)GR，Q为隐层特征空 间的维数，韦
，其中，Al和Cl为满足上式的 特征参数，t是S的长度自变量，t=l，2,…，dXT，并且v(t)GRdxT独立同分布于噪声概率密 度q( ?)，建模结果如下；
[0094]
[0095] 令Mi = (Ai，Ci)，即为子触觉序列的动态特征，即某一频率范围内的触觉序列；
[0096] 步骤S3、使用马下距离(Madin Distance)作为量度计算子触觉序列的动态特征 之间的距离，获得马下矩阵；
[0097] 如图4所示，对于提取特征存在于欧式空间的分类算法，其特征之间的距离使用欧 式距离来衡量，然而，对于两个LDS的特征化=(心，打），12=^2，〔2)，因为其存在于非欧式空 间中，所W使用马下距离来进行LDS特征之间的距离衡量，得到马下距离矩阵，md表征输入 样本的触觉子序列个数；
[0098] 马下距离是基于两个系统之间的空间角定义的，运种空间角又称为观测子序列模 型参数的规则角(Principle Angles)，在本发明构建的系统中，即触觉序列的LDS特征之间 的规则角，该规则角可W由下式进行定义：
[0099]
(3)
[0100] 对于任意两个模型之间规则角的计算，使用李雅普诺夫化yapunov)方程求解，保 证系统稳定的重要条件是存在正定的对称矩阵P满足下式：
[0101] ATpQ-P = -CTc (4)
[0102] 式中
,
[0103] 诚巨/击田了味^+锋辛台|、田於T金3走估.
[0104] (5)
[0105] 景终得到两个LDS据佈Ml与]间的马下距离dM(Mi，M2):
[0106]
(：6)
[0107] 使用上述步骤对所提取出的LDS特征数据进行计算后，得到特征间的马下距离：
[010 引
（7)
[0109] 步骤S4、便用K-中屯、点（K-Medoid)算法对与J矩阵进行聚类，构建码书 (Codebook)；
[0110] K-均值化-Means)算法和K-Medoid算法均可用于聚类，相较于K-均值化-Means)算 法，K-Medoid算法是一种基于数据划分的聚类算法，使用簇中最靠近样本中屯、的数据点来 代表该簇，而不是使用数据点中屯、自身，W此避免了噪声和冗余数据对于算法的影响，但是 由于两个LDS特征之间的距离使用马下距离进行衡量，而马下距离存在于非欧式空间，K- Means需将上面进行计算的马下距离转换到欧式空间中，其过程较繁琐，并且K-Means算法 自身也存在对于原始条件过度依赖的缺点，例如聚类中屯、点数量k等都会对聚类效果造成 影响，综合上述因素，使用K-Medoid算法完成聚类，K-Medoid算法还有一个优点是使用简单 方便且易实现；
[0111] 使用K-中屯、点算法对马下矩阵进行聚类的步骤包含：
[0112] 步骤S4.1、随机选取k个马下距离矩阵中的数据作为聚类中屯、点；
[0113] 步骤S4.2、将矩阵中剩余的数据按照所选取的聚类中屯、点进行分组，分组W每组 中每个数据点与每组中聚类中屯、点的距离数据值的平方和最小为原则；
[0114] 步骤S4.3、重复执行步骤S4.1、步骤S4.2和步骤S4.3,更新每组数据样本的聚类中 屯、点；
[0115] 步骤S4.4、直至聚类中屯、点稳定不变，完成数据聚类；
[0116] 聚类完成后得到由4组0)5特征组成的码书((：〇(166〇〇1〇，^[(；4，（：,.^；^为表示0)8 特征的元组，对于码书中的任意一组元组(A，C)，称其为码词(Codewords);
[0117] 步骤S5、使用码书对每组触觉序列进行表征，得到外部特征为直方图的系统包 (Bag-〇f-System，BoS)模型；
[011引使用码书对触觉序列进行表征后可W得到系统包模型{BoS ，m为样本中触觉
序列的个数，如图5所示，运种表征方式的外部体现为直方图模聖 可由特征词频率(Term化equency，TF)算法得到；
[0119] 假设在第i组触觉序列中，第j组码词出现的次数为CU次，则有下式：
[0120]
(8)
[0121] 巧中，hu表不在第i組触觉序列中，第j組妈词出现的频率，m为样本中触觉序列个 数，对于hi则有hi=比ir，Aik]，为一组特征向量；
[0122] 运种算法过程简单且严谨，并且运种算法的归一化处理有助于提高样本的分类精 度，图5中也对于系统包模型的内部数据形式进行了显示；
[0123] 经过上述过程，得到触觉序列系统包，它由m组特征向量组成；
[0124] 步骤S6、将训练集样本的系统包模型和训练集样本标签一起送入极限学习机 化XTREME LEARNING MACHI肥，EL…由祉I链A苯号（即分类函数f);
[01巧]假定存在一组触觉序巧 ，则该触觉序列可W使用具有办个隐层结 点且扩展到RBF核函数的单隐层核函数网络(Single-hidden Layer Fee壯ordward化ural Network，SLFN)进行建模，函数模型如下式所示：
[0126]
(说
[0127] 式中，CO =[ COr-?"]T是连接第i个核函数与输出神经网络的输出权重向量，(61 (rj)为输出函数，r=[rr-rm]T为相应的特征向量，化为函数中屯、点，Oi为影响因子；
[012引高斯核巧数可W被定父为：
[0129] (10)
[0130] 1义攸I化子 >」WbLM具巧歩/|、拖!巧结点可W完成对数据的多分类，也就是使进行M种 类别识别的分类器具有M组隐含层结点，ELM算法对于属于该类别的样本使用+1进行标记， 对于不属于该类的样本使用0进行标记，与一般的分类算法不同的是，ELM算法在使训练样 本误差最小的同时，也使样本的输出权重最小，即最小化MHco-TM2与M CO M的值，式中H 为神经网络隐含民输出巧阵，它可由下式讲行表示；
[0131] (11)
[0側 《巧)
[0133] 3)
[0134] 式中，矩阵分t为神经网络隐含层输出矩阵H的广义逆矩阵，在训练具有多输出结 点的ELM分类器时遵循如下条件：
[0135] 最小 i (14)
[0136] 约束： U5)上式 中，C亦为损耗参数，可由交叉验证获得;Ci = Km,…，Cl,k]T为训练集误差向量，此时，训练 ELM分类器等效于解决下式的最优化问题：
[0137]
,)
[0140] (17)
[0138] 式中，为训练样本相应的拉格朗日系数；[0139] 当训练集样本规模不大时，可W得到下式进行描述：
[0141]
[01创 （献
[0143] ;为
[0144] r 19)
[0145] 使用Qelm = HHT对极限学习机的核函数矩阵进行定义，并且满足W下核函数公式 QELMa,b=(Hra)(Hrb)=K(ra，rb)，使用核函数的ELM分类器输出函数可W按照下式进行定 义
[0146] (撕）
[0147] 通过上述算法可完成对具有多输出结点的ELM分类器的训练，并且可W对测试集 样本的物理标签进行预测，所预测的标签即为输出的结点所对应的索引，此索引对应上式 中最大的输出值，且有ltest = arg maxfjjGU，. . .m}，式中。是第j个输出结点的输出函 数；
[0148] 步骤S7、利用触觉传感器获取待分类样本的触觉序列，对待分类样本的触觉序列 按照步骤S2~步骤S5的方法进行处理，获得待分类样本的系统包模型；
[0149] 步骤S8、将待分类样本的系统包模型送入步骤S6获得的分类器中，即可得到物体 类型的标签。
[0150] 样本对应的标签为Ii,其中i = l，2,…，N，标签是物体的种类信息，可用数字、文字 或字母表示，例如：16种样本，用数字表示分类，那么样本对应的标签为Il = I，12 = 2,13 = 3,...，1i6=16;
[0151] 对于物体样本的分类准确率见图6,横轴和纵轴分别为式（15)所述损耗参数C和式 (9)所述隐层节点# :，.根据不同的C和#分类准确率均高于80%。
[0152] 本发明中的触觉传感器安装在机器人的机械臂末端的灵巧手上，在操作过程中， 当机械臂系统到达指定位置时，安装在机械臂末端的灵巧手执行抓取操作。由于待抓取的 目标物体形状、材质等区别较大，因此，利用安装在灵巧手上的触觉传感器阵列准确获得抓 取过程中的触觉信息，自适应调整"握力"大小，确保可W稳定的抓取不同的物体。
[0153] 物体材质是机器人精细操作过程中一个非常重要的环节。如果能够通过触觉传感 器有效地判别物体的材质，就可W根据其特点来采取有针对性的抓取和操作策略。W往对 材质的分析大都是基于图像纹理的，运一信息源只能刻画单一角度和距离的纹理信息。在 机械臂末端安装的摄像机可W在机械臂末端接近物体过程中不断采集物体的图像，从而形 成关于物体的视频序列，运一数据源可W从不同角度和距离感知物体的纹理信息。尽管如 此，运两类信息源都是在非接触式的情形下采集的，对于一些难W通过视觉方式鉴别的物 体(如用材质做出来的外观相似的物体)效果不佳。本发明结合机器人精细操作平台引入了 另一个触觉数据源，即在机械臂正式抓取物体之前，利用手指轻触物体，利用运种接触式的 触觉序列信息，可W有效地弥补视觉信息在材质分析任务上的不足。
[0154] 本发明解决了机器人对非合作目标稳定、柔顺抓取的实际需求，为精细操作任务 的完成提供数据基础，并可与其它传感结果融合计算，从而通过多源深度感知增强对不同 目标物的描述和辨识能力，为智能操控的实施奠定技术基础。
[0155] 尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的 描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的 多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
【主权项】
1. 一种基于机器学习的触觉信息分类计算建模方法，其特征在于，包含以下步骤： 步骤S1、利用触觉传感器获取训练集样本的触觉序列； 步骤S2、将每组触觉序列合并为子触觉序列，使用线性动态系统模型对每组子触觉序 列进行建模，分别提取子触觉序列的动态特征； 步骤S3、使用马丁距离作为量度计算子触觉序列的动态特征之间的距离，获得马丁矩 阵； 步骤S4、使用K-中心点算法对马丁矩阵进行聚类，构建码书； 步骤S5、使用码书对每组触觉序列进行表征，得到外部特征为直方图的系统包模型； 步骤S6、将训练集样本的系统包模型和训练集样本标签一起送入极限学习机中训练分 类器； 步骤S7、利用触觉传感器获取待分类样本的触觉序列，对待分类样本的触觉序列按照 步骤S2~步骤S5的方法进行处理，获得待分类样本的系统包模型； 步骤S8、将待分类样本的系统包模型送入步骤S6获得的分类器中，即可得到物体类型 的标签。2. 如权利要求1所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S1中，所述的触觉阵列传感器设置在机器人的机械臂末端的灵巧手上，所述的灵巧 手具有三个手指，每个手指末端分别配置一个触觉阵列传感器，每个触觉阵列传感器具有m Xn个传感器单元，所述的触觉序列为通过触觉阵列传感器获得的给定N个训练时间序列样 本 SieRdxT，其中，d是空间维度，T是序列的时间长度，i = l，2,…，N。3. 如权利要求1所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S2中，首先通过降低维数对触觉序列的主要成分进行提取，去除冗杂的成份。4. 如权利要求1或2所述的基于机器学习的触觉信息分类计算建模方法，其特征在于， 所述的步骤S2中， 通过重塑m X η矩阵作为一个m X η维向量，分别获得灵巧手上三个手指的(m X η)维动态 触觉序列 S(1) e R(mXn)XT，s⑵ e R(mXn)XT，s⑶ e R(mXn)XT; 串联3个手指的mXn维触觉序列，合并为一个单一的子触觉序列：(1) 其中，d=mXnX3; 假设噪声为w(t)eRdxT，将子触觉序列S分解为:S(t)=D(t)+w(t)，其中D(t)eRdxT为不 含噪声部分； 假设存在一个p维滤波器如:允->贫〃7，使得冰)￡凡0为 隐层特征空间的维数，有D(t)=^〇<0)，使得其中，Ai和Ci为满足上式的特征参数，t是S的长度自变量，t = 1，2，…，d X T，并且v(t) e (2) RdXT独立同分布于噪声概率密度q( ·)，建模结果如下： 令Mi = (Ai，Ci)，即为子触觉序列的动态特征，即某一频率范围内的触觉序列。5. 如权利要求4所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S3中，使用马丁距离来进行两个LDS的特征吣=(六 1，(：1)，12=(如，(：2)之间的距离衡 量，得到马丁距离矩阵，md表征输入样本的触觉子序列个数； 触觉序列的LDS特征之间的规则角由下式进行定义：(3) 对于任意两个模型之间规则角的计算，使用李雅普诺夫方程求解，保证系统稳定的重 要条件是存在正定的对称矩阵P满足下式： ATPQ-P = -CTC (4) 式牛 V21 U \7 然后使用下式计算子空间角{汉^^的余弦值：(5) 最终得到两个LDS特征施与此之间的马丁距离(66) 使用上述步骤对所提取出的LDS特征数据进行计算后，得到特征间的马丁距离：(7)〇6. 如权利要求5所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S4中，使用K-中心点算法对马丁矩阵进行聚类的步骤包含： 步骤S4.1、随机选取k个马丁距离矩阵中的数据作为聚类中心点； 步骤S4.2、将矩阵中剩余的数据按照所选取的聚类中心点进行分组，分组以每组中每 个数据点与每组中聚类中心点的距离数据值的平方和最小为原则； 步骤S4.3、重复执行步骤S4.1、步骤S4.2和步骤S4.3，更新每组数据样本的聚类中心 占 . 4\\\ ， 步骤S4.4、直至聚类中心点稳定不变，完成数据聚类； 聚类完成后得到由k组LDS特征组成的码书((：〇(1必〇〇1〇，|(4，€^0=1为表示LDS特征的 元组，对于码书中的任意一组元组(A，C)，称其为码词(Codewords)。7. 如权利要求6所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S5中，使用码书对触觉序列进行表征后可以得到系统包模型丨B〇S ，m为样本中 触觉序列的个数，这种表征方式的外部体现为直方图模型為=[/*ρ/ν··，/?,"]Τ eUr，可由特 (8) 征词频率算法得到； 假设在第i组触觉序列中，第j组码词出现的次数为(?次，则有下式： 式中，hij表示在第i组触觉序列中，第j组码词出现的频率，m为样本中触觉序列个数，对 于hi则有hi= [hn. . .hik]，为一组特征向量。8.如权利要求7所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S6中， 假定存在一组训练触觉序列，则该触觉序列可以使用具有#个隐层结点且 扩展到RBF核函数的单隐层核函数网络进行建模，函数模型如下式所示：(9) 式中，ω = [ ω . . com]T是连接第i个核函数与输出神经网络的输出权重向量，φ i(rj) 为输出函数，r=[n.. .rm]T为相应的特征向量，yi为函数中心点，〇i为影响因子； 高斯核函数可以被定义为：(10) 使极限学习机具有多个输出结点可以完成对数据的多分类，也就是使进行Μ种类别识 别的分类器具有Μ组隐含层结点，极限学习机算法对于属于该类别的样本使用+1进行标记， 对于不属于该类的样本使用〇进行标记，极限学习机算法在使训练样本误差最小的同时，也 使样本的输出权重最小，即最小化I |Ηω-Τ| |2与| | ω | |的值，式中Η为神经网络隐含层输出 矩阵，它可由下式进行表示：式中，矩阵为神经网络隐含层输出矩阵Η的广义逆矩阵，在训练具有多输出结点的极 限学习机分类器时遵循如下条件： 最小化（14) ' X i=1 约束条饵(15) 上式中，C亦为损耗参数，可由交叉验证获得. . .，|1>k]T为训练集误差向量， 此时，训练极限学习机分类器等效于解决下式的最优化问题：式中，为训练样本相应的拉格朗日系数； 当训练集样本规模不大时，可以得到下式进行描述：(17) 并可以由此得到ω :(18) 此时，极限学习机算法的输出函数可以表示为(19 使用ΩΜ=ΗΗ4?极限学习机的核函数矩阵进行定义，并且满足以下核函数公式QELMa, b =Φ (ra) Φ (n) =K(ra，rb)，使用核函数的极限学习机分类器输出函数可以按照下式进行 定义(20) 通过上述算法可完成对具有多输出结点的极限学习机分类器的训练，并且可以对测试 集样本的物理标签进行预测，所预测的标签即为输出的结点所对应的索引，此索引对应上 式中最大的输出值，且有ltest = arg max fj，je {1，. . .m}，式中fj是第j个输出结点的输出 函数。9.如权利要求8所述的基于机器学习的触觉信息分类计算建模方法，其特征在于，所述 的步骤S8中，样本对应的标签为1 i，其中i = 1，2，…，N，标签是物体的种类信息，可用数字、 文字或字母表示。
【文档编号】G06F17/50GK105956351SQ201610522439
【公开日】2016年9月21日
【申请日】2016年7月5日
【发明人】侯月阳, 卢山, 田路路, 王奉文, 于学文
【申请人】上海航天控制技术研究所