一种肌电信号特征提取方法技术

本发明专利技术公开了一种肌电信号特征提取方法，首先对肌电信号作傅里叶变换取傅里叶系数，选取预设频段的傅里叶系数并归一化；构建堆叠自动编码器，设置每一层网络参数进行训练更新参数；构建堆叠的受限玻尔兹曼机，设置每一层网络参数进行训练更新参数；在预训练自动编码器和受限玻尔兹曼机完成后，将两部分预训练的参数前后连接起来，即自动编码器的输出为受限玻尔兹曼机的输入，形成一个贯通的网络，实现肌电信号特征提取，解决现有肌电信号提取方法时面临的时域特征变化大、频域特征提取不充分的问题。

A Feature Extraction Method for EMG Signals

The invention discloses a method for feature extraction of EMG signals. Firstly, the Fourier coefficients of EMG signals are extracted by Fourier transform, and the Fourier coefficients of preset frequency bands are selected and normalized; stacked automatic encoders are constructed to set the network parameters of each layer for training and updating parameters; stacked restricted Boltzmann machines are constructed to set the network parameters of each layer for training and updating parameters. After the completion of pre-training automatic encoder and limited Boltzmann machine, the two parts of pre-training parameters are connected, that is, the output of automatic encoder is the input of limited Boltzmann machine, forming a through network, realizing feature extraction of EMG signal, and solving the problems of time-domain feature change and insufficient feature extraction in frequency-domain faced by existing methods of extracting EMG signal.

【技术实现步骤摘要】
一种肌电信号特征提取方法
本专利技术涉及信号特征提取
，具体涉及一种肌电信号特征提取方法。
技术介绍
肌电信号能在一定程度上反映神经肌肉的活动，在临床医学、人机功效学、康复医学以及体育科学等方面均有重要的实用价值。特征提取能发现更具意义的潜在变量，有助于更加深入地了解肌电信号数据，减少数据存储和输入数据带宽，减少冗余。传统的肌电信号处理的方法是将肌电信号看成均值为零、方差随信号强度变化而变化的随机信号。时域特征的提取相对比较简单，故时域分析方法在肌电信号应用领域得到了比较广泛的应用。然而，时域特性虽然容易提取，但大量研究表明，当肌肉收缩力大小稍有变化时，表面肌电信号的时域特征变化较大，并不稳定，而通过快速傅里叶变换将信号变换成频域中的频谱或功率谱后，频谱或功率谱的波形变化较小，即肌电信号的频域描述相对比较稳定。肌电信号的频谱稳定性使得频域方法成为肌电信号处理技术的主流，得到了广泛的应用。但是肌电信号的复杂性决定了单独使用某一类方法进行肌电信号处理都可能无法充分利用所得信息。
技术实现思路
为解决上述问题，本专利技术提出一种肌电信号特征提取方法，实现肌电信号特征提取，解决现有肌电信号提取方法时面临的时域特征变化大、频域特征提取不充分的问题。本专利技术采用如下技术方案，一种肌电信号特征提取方法，具体步骤如下：1)对肌电信号作傅里叶变换取傅里叶系数，选取预设频段的傅里叶系数并归一化；2)构建堆叠自动编码器，设置每一层网络参数进行训练更新参数；3)构建堆叠的受限玻尔兹曼机，设置每一层网络参数进行训练更新参数；4)在预训练自动编码器和受限玻尔兹曼机完成后，将两部分预训练的参数前后连接起来，即自动编码器的输出为受限玻尔兹曼机的输入，形成一个贯通的网络，实现肌电信号特征提取。优选地，所述步骤1)中所述肌电信号为肌电信号通过滤波进行除噪后的肌电信号。优选地，所述步骤1)中所述预设频段为11HZ～50HZ。优选地，所述步骤2)对自动编码器进行训练采用随机梯度下降法，具体步骤为：21)设置自动编码器的网络参数；22)将训练数据分成含多个样本的小批量数据进行分批次训练；23)在前馈网络中根据权重和偏置计算自动编码器的输出；24)使用反向传播算法反向传播，计算损失函数，损失函数对权重和偏置求偏导，再根据随机梯度下降法更新每个连接的权重或偏置，损失函数为：其中m是每次小批量数据的样本个数，y(x)是训练数据对应的真实标签，a(x)是当前权重或偏置下的输出值；损失函数分别对权重或偏置求偏导得到梯度值，使用随机梯度下降法更新权重；25)堆叠多个自动编码器，上一个自动编码器的隐藏层作为下一个自动编码器的输入，循环步骤21)-24)进行训练。优选地，所述网络参数包括：学习速率、动量项和权值惩罚项。优选地，所述步骤3)中采用梯度下降法进行训练更新参数，具体步骤为：31)定义一个能量函数E(v，h)，E(v，h)＝-aTv-cTh-hTuv其中v为可视层输入，h为隐藏层，a、c分别表示可视层和隐藏层的偏置，u表示可视层与隐藏层的权重。根据能量函数E(v，h)得到可视层与隐藏层的联合概率分布32)基于贝叶斯公式由联合概率分布得到条件概率分布，使用似然函数计算损失函数；33)设置受限玻尔兹曼机的网络参数；34)训练的目标为最大化步骤32)中的似然函数，似然函数对参数取偏导来，利用梯度下降法更新参数。专利技术所达到的有益效果：本专利技术是一种肌电信号特征提取方法，实现肌电信号特征提取，解决现有肌电信号提取方法时面临的时域特征变化大、频域特征提取不充分的问题。本专利技术将自动编码器应用于肌电信号分析领域，确保输出的低维度特征能够最大限度地表征输入的高维度数据，抓住了输入数据的特点，使其特征保持不变，减少数据存储和输入数据带宽，减少冗余，充分提取肌电信号。附图说明图1是本专利技术实施例的一种肌电信号特征提取方法的流程示意图。具体实施方式下面根据附图并结合实施例对本专利技术的技术方案作进一步阐述。本实施例所用的肌电信号数据是由四个人做十个不同的手势动作获取的肌电信号数据，每个人的不同的手势表示一个类别，每个手势都是一个8维数据，每一维数据有10000个样本点。一种肌电信号特征提取方法，构建堆叠自动编码器网络架构，即堆叠多个自动编码器，每一个自动编码器都是一个三层神经网络，设置每一层神经网络的参数且利用随机梯度下降法(Stochasticgradientdescent，SGD)进行训练更新权重；建立只有可视层和隐藏层的受限玻尔兹曼机(restrictedBoltzmannmachine，RBM)，使用梯度下降法进行训练更新权重。训练完自动编码器和受限玻尔兹曼机后，将两部分预训练的权重参数前后连接起来形成一个贯通的全连接网络，实现肌电信号特征提取。自动编码器将神经网络的隐含层视为一个编码器和解码器，输入数据经过隐含层的编码和解码，到达输出层时，确保输出的低维度特征能够最大限度地表征输入的高维度数据。受限玻尔兹曼机是一种随机神经网络，用于学习一组数据的内在表示，通过使用数据集设置可视层上的神经元来训练受限玻尔兹曼机以匹配数据集中的数据点，训练后的受限玻尔兹曼机可以用于对新的数据进行分类。一种肌电信号特征提取方法，如图1所示，具体步骤如下：1)对肌电信号作傅里叶变换取傅里叶系数，选取预设频段的傅里叶系数并归一化；作为一种较佳的实施例，所述步骤1)中所述肌电信号为肌电信号通过滤波进行除噪后的肌电信号。所述步骤1)中所述预设频段为11HZ～50HZ。将归一化后的傅里叶变换系数随机划分为训练数据和测试数据，标记部分数据作为理想的输出数据，即标签数据。归一化可以加快梯度下降，加快最优解的速度，提高精度。2)构建堆叠自动编码器，设置每一层网络参数进行训练更新参数；对自动编码器进行训练采用随机梯度下降法，具体步骤为：21)设置自动编码器的网络参数；所述网络参数包括：学习速率、动量项和权值惩罚项；学习速率设为0.05，用以控制学习速度；设置动量项为0.01，调节在小批量训练网络时本次更新值与上一次更新值的相对比例，在梯度陡峭的地方就下降地快，梯度小的地方下降地慢，避免错过梯度下降法的最优解；设置权值惩罚项防止参数(权重或偏置)过大，避免过拟合现象；将激活函数设置为sigmoid函数；22)将训练数据分成含多个样本的小批量数据进行分批次训练；既能加快训练速度，又能使模型收敛到较优的解。本实施例中每个小批量数据中样本的个数为32，设置样本训练次数为200，即epoch迭代次数，迭代200次；23)在前馈网络中根据权重和偏置计算自动编码器的输出，输出等于权重与输入的乘积，结果再加上一个偏置，最后用激活函数对结果进行映射。具体为：y′＝f(wx+b)，其中y’表示输出，x为前馈网络输入数据，w为前馈网络的权重，f表示激活函数，本实施例中采用sigmoid函数，参数b为偏置，在计算过程中，需要“复制”b维数，复制后的b和wx乘积结果可以相加；但是这样每次复制b，会增大计算量，本实施例中将偏置项添加到权重w中变为w’，对输入数据x前一列添加一列1元素得到x’，此时的w′x′即是(wx+b)；24)使用反向传播算法反向传播，计算损失函数，损失函数对权重和偏置求偏导，再根据随机梯度下降法更新每个连接的权本文档来自技高网...

【技术保护点】
1.一种肌电信号特征提取方法，其特征在于，包括以下步骤：1)对肌电信号作傅里叶变换取傅里叶系数，选取预设频段的傅里叶系数并归一化；2)构建堆叠自动编码器，设置每一层网络参数进行训练更新参数；3)构建堆叠的受限玻尔兹曼机，设置每一层网络参数进行训练更新参数；4)在预训练自动编码器和受限玻尔兹曼机完成后，将两部分预训练的参数前后连接起来，即自动编码器的输出为受限玻尔兹曼机的输入，形成一个贯通的网络，实现肌电信号特征提取。

【技术特征摘要】
1.一种肌电信号特征提取方法，其特征在于，包括以下步骤：1)对肌电信号作傅里叶变换取傅里叶系数，选取预设频段的傅里叶系数并归一化；2)构建堆叠自动编码器，设置每一层网络参数进行训练更新参数；3)构建堆叠的受限玻尔兹曼机，设置每一层网络参数进行训练更新参数；4)在预训练自动编码器和受限玻尔兹曼机完成后，将两部分预训练的参数前后连接起来，即自动编码器的输出为受限玻尔兹曼机的输入，形成一个贯通的网络，实现肌电信号特征提取。2.根据权利要求1所述的一种肌电信号特征提取方法，其特征在于，所述步骤1)中所述肌电信号为肌电信号通过滤波进行除噪后的肌电信号。3.根据权利要求1所述的一种肌电信号特征提取方法，其特征在于，所述步骤1)中所述预设频段为11HZ～50HZ。4.根据权利要求1所述的一种肌电信号特征提取方法，其特征在于，所述步骤2)对自动编码器进行训练采用随机梯度下降法，具体步骤为：21)设置自动编码器的网络参数；22)将训练数据分成含多个样本的小批量数据进行分批次训练；23)在前馈网络中根据权重和偏置计算自动编码器的输出；24)使用反向传播算法反向传播，计算损失函数，损失函数对权重和偏置求偏导，再根...

【专利技术属性】
技术研发人员：徐宁，刘妍妍，倪亚南，潘安顺，刘小峰，
申请(专利权)人：河海大学常州校区，
类型：发明
国别省市：江苏,32