基于时域及相空间的脑电放松度识别方法及装置与流程
本发明涉及放松治疗领域，尤其涉及一种基于时域及相空间的脑电放松度识别方法及装置。
背景技术：
放松训练是行为疗法中使用最广的技术之一，是在心理学实验的基础上建立和发展起来的咨询和治疗方法，其在治疗焦虑抑郁症、神经性头痛、失眠、高血压病，减轻更年期综合征和转变不良行为模式等方面取得了较好的疗效。
现有的放松训练主要有录音指导、口头指导和生物反馈指导。其中，录音指导方法僵化、没有变化，无法根据受训者的状态变化内容；口头指导则要求对口头指导的对象要求很高，且受到时间、场地限制；生物反馈指导以脑电反馈为主，能够结合前两种方式的优点，因而受到广泛关注。
进行生物反馈指导需要识别用户的放松度，而计算放松度首先需要从用户的脑电信号中提取各个频段的脑电波(包括delta、theta、alpha、beta、gamma波)，再提取各个脑电波的特征，将这些特征输入到分类器进行分类识别。
现有的特征提取方法一般只能从单一角度提取脑电波的特征，评价方式单一，不能保证分类结果的准确性。而且现有的特征提取算法的计算和处理过程复杂，一方面，增加了对硬件的要求，另一方面，由于计算复杂，也无法及时得到分类结果，进而影响了放松治疗的效果。
技术实现要素：
针对上述问题，本发明的目的在于提供一种基于时域及相空间的脑电放松度识别方法及装置，可全面的提取出各个脑电波的特征。
本发明提供了一种基于时域及相空间的脑电放松度识别方法，包括如下步骤：
对接收的原始脑电信号进行切片，生成具有指定时间长度的待处理脑电序列信号；
从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波；
计算对应于各个脑电波的信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，获取所述待处理脑电序列信号在时域空间的特征量；
根据所述各个脑电波的信号波及所述待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量；
根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度。
优选地，所述计算对应于各个脑电波的信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，获取所述待处理脑电序列信号在时域空间的特征量具体包括：
计算各个信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，计算得到每个信号波的等电位线能量参数变化率基线；其中，所述能量参数至少包括以下其中之一：能量、能量密度；
基于最小二乘法计算每个信号波在各个时刻的能量参数变化率；
统计所述每个信号波的能量参数变化率小于所述等电位线能量参数变化率基线的个数；
根据统计得到的对应于每个信号波的个数与对应的信号波的总采样点数的比值，获取所述待处理脑电序列信号在时域空间的特征量。
优选地，所述根据所述各个脑电波的信号波及待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量具体包括：
根据所述各个脑电波的信号波及所述待处理脑电序列信号形成对应的包含有信号的二维图表；
在各个二维图表上覆盖与所述二维图表等大小的m*m的网格对应，并统计覆盖有信号的网格数；其中，m为大于1的整数；
根据覆盖有信号的网格数与所述网格的全部网格数的比值，生成各个信号波及待处理脑电序列信号的相空间分布密度，获得待处理脑电序列信号在相空间的特征量。
优选地，在所述根据所述各个脑电波的信号波及待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量之后，所述根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度之前，还包括：
基于方差分析和f检验计算每个特征量的显著性指标，并选取显著性指标大于预设阈值的特征量。
优选地，所述基于方差分析和f检验计算每个特征量的显著性指标，并选取显著性指标大于预设阈值的特征量具体包括：
根据所述特征量、由标准设备同步采集得到的脑电放松度及预先拟合的与每个特征量对应的线性拟合曲线进行方差分析，计算得到每个特征量的拟合值与期望的平方和及原值与拟合值的平方和；
根据拟合值与期望的平方和、原值与拟合值的平方和及自由度，计算得到每个特征量的显著性指标；
根据每个特征量的显著性指标，选取出显著性指标大于预设阈值的特征量。
优选地，在所述根据所述各个脑电波的信号波及待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量之后，所述根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度之前还包括：
基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行降维处理，获得降维后的特征量；
其中，所述基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行降维处理，获得降维后的特征量具体包括：
将所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量设置为输入样本空间中的特征量，并对所述输入样本空间进行数据标准化处理；
根据数据标准化处理后的所述输入样本空间得到协方差矩阵；
计算所述协方差矩阵的特征根和与每个特征根对应的特征向量；其中，所述特征根的数量为p个，且所述的p个特征根按照大小排序；
获取所述的p个特征根中，贡献率之和大于预定阈值的前m个特征根；其中，每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和；
根据与所述的前m个特征根对应的特征向量及所述输入样本空间，得到主成分得分矩阵；其中，所述主成分得分矩阵中的特征量为所述降维后的特征量。
优选地，所述根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度具体包括：
基于由不同学习算法训练得到的至少两个神经网络模型对输入的所有特征量进行分类，得到待处理脑电序列信号在各个神经网络模型下的分类；
将出现次数最多的分类设置为所述待处理脑电序列信号的分类；
根据所述分类与脑电放松度的对应关系，识别得到与所述待处理脑电序列信号对应的脑电放松度。
本发明还提供了一种基于时域及相空间的脑电放松度识别装置，包括：
信号切片单元，用于对接收的原始脑电信号进行切片，生成具有指定时间长度的待处理脑电序列信号；
信号提取单元，用于从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波；
第一特征提取单元，用于计算对应于各个脑电波的信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，获取所述待处理脑电序列信号在时域空间的特征量；
第二特征提取单元，用于根据所述各个脑电波的信号波及待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量；以及，
脑电放松度识别单元，用于根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度。
优选地，所述第一特征提取单元具体包括：
变化率基线计算模块，用于计算各个信号波的能量参数，并基于各个信号波的能量参数的最大值及每个脑电波的等电位线能量参数变化率比率，计算得到每个信号波的等电位线能量参数变化率基线；
能量参数变化率计算模块，用于基于最小二乘法计算每个信号波在各个时刻的能量参数变化率；
统计模块，用于统计所述每个信号波的能量参数变化率小于所述等电位线能量参数变化率基线的个数；以及，
比值计算模块，用于根据统计得到的对应于每个信号波的个数与对应的信号波的总采样点数的比值，获取所述待处理脑电序列信号在时域空间的特征量。
优选地，所述第二特征提取单元具体包括：
二维图表生成模块，用于根据所述各个脑电波的信号波及所述待处理脑电序列信号形成对应的包含有信号的二维图表；
网格数统计模块，用于在各个二维图表上覆盖与所述二维图表等大小的m*m的网格对应，并统计覆盖有信号的网格数；其中，m为大于1的整数；以及，
特征量计算模块，用于根据覆盖有信号的网格数与所述网格的全部网格数的比值，生成各个信号波及待处理脑电序列信号的相空间分布密度，获得待处理脑电序列信号在相空间的特征量。
优选地，所述基于时域及相空间的脑电放松度识别装置还包括：
特征选取单元，用于基于方差分析和f检验计算每个特征量的显著性指标，并选取显著性指标大于预设阈值的特征量。
优选地，所述特征选取单元具体包括：
方差分析模块，用于根据所述特征量、由标准设备同步采集得到的脑电放松度及预先拟合的与每个特征量对应的线性拟合曲线进行方差分析，计算得到每个特征量的拟合值与期望的平方和及原值与拟合值的平方和；
显著性指标计算模块，用于根据拟合值与期望的平方和、原值与拟合值的平方和及自由度，计算得到每个特征量的显著性指标；以及，
选取模块，用于根据每个特征量的显著性指标，选取出显著性指标大于预设阈值的特征量。
优选地，所述基于时域及相空间的脑电放松度识别装置还包括：
特征降维单元，用于基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行降维处理，获得降维后的特征量；
其中，所述特征降维单元具体包括：
标准化处理模块，用于将所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量设置为输入样本空间中的特征量，并对所述输入样本空间进行数据标准化处理；
协方差矩阵计算模块，用于根据数据标准化处理后的所述输入样本空间得到协方差矩阵；
特征计算模块，用于计算所述协方差矩阵的特征根和与每个特征根对应的特征向量；其中，所述特征根的数量为p个，且所述的p个特征根按照大小排序；
筛选模块，用于获取所述的p个特征根中，贡献率之和大于预定阈值的前m个特征根；其中，每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和；以及，
降维特征量获得模块，用于根据与所述的前m个特征根对应的特征向量及所述输入样本空间，得到主成分得分矩阵；其中，所述主成分得分矩阵中的特征量为所述降维后的特征量。
优选地，所述脑电放松度识别单元具体包括：
分类模块，用于基于由不同学习算法训练得到的至少两个神经网络模型对输入的所有特征量进行分类，得到待处理脑电序列信号在各个神经网络模型下的分类；
分类统计模块，用于将出现次数最多的分类设置为所述待处理脑电序列信号的分类；以及，
放松度识别模块，用于根据所述分类与脑电放松度的对应关系，识别得到与所述待处理脑电序列信号对应的脑电放松度。
本发明提供的基于时域及相空间的脑电放松度识别方法及装置，同时从时域空间、相空间两个角度提取脑电序列序号的特征量，并基于提取得到的特征量进行分类识别得到最终的脑电放松度。相比于单一角度的特征提取，对信号的评价方式更多样化，能更全面的体现信号的特性，避免了单一角度的特征提取容易导致的特征过于片面的问题而影响最终识别精度的问题。本发明可以提高识别分类的精度和准确度，为放松治疗提供准确、快速的依据。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于时域及相空间的脑电放松度识别方法的流程示意图。
图2是通过切片得到待处理脑电序列信号的示意图。
图3是本发明实施例提供的对对原始脑电序列信号进行加权移动平均计算的原理图。
图4是本发明实施例提供的基于时域及相空间的脑电放松度识别装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本发明实施例提供了一种基于时域及相空间的脑电放松度识别方法，其可包括如下步骤：
s101，对接收的原始脑电信号进行切片，生成具有指定时间长度的待处理脑电序列信号。
具体地：对原始脑电信号进行切片，获得时间长度为30秒的待处理脑电序列信号，并对所述待处理脑电序列信号的各个时刻的脑电信号进行加权移动平均计算，得到去除低频直流信息后的待处理脑电序列信号。
如图2所示，在本发明实施例中，所述待处理脑电序列信号可通过对原始脑电信号进行切片获得。其中，所述原始脑电信号可通过脑电电极采集获得。一般地，脑电电极采集的原始脑电信号的持续时间比较长(如数小时甚至更长)，因此进行需要对原始脑电信号进行切片，例如，每个切片的片段为30秒，即每段所述待处理脑电序列信号的长度为30秒。
在本优选实施例中，为了保证提取和滤波的效率和准确性，还可对脑电序列信号进行相应的预处理，例如，通过预处理去除待处理脑电序列信号中的低频直流信息，以避免这些低频直流信息对脑电波提取的干扰。
在本优选实施例中，为了去除原始脑电序列信号中的低频直流信息，可基于加权移动平均算法对降采样后的原始脑电序列信号的各个时刻的脑电信号进行计算，得到所述待处理脑电序列信号。具体地：
首先，基于当前的第j个时刻的脑电信号，获取所述原始脑电序列信号中的位于第(j-(n-1)/2)个时刻至第(j+(n-1)/2)个时刻之间的n个脑电信号的能量；其中，n为预设的影响数，且n为奇数，j为大于(n+1)/2的整数。
例如，假设当前要预测的脑电信号x(j)的时刻为第10个时刻(即j＝10)，影响数n为5，则对当前要预测的脑电信号有影响的脑电信号为第8个至第12个时刻的脑电信号，即x(8)～x(12)。此时，先获取这5个时刻的脑电信号的能量。
接着，根据预设的权值分布函数为获取的n个脑电信号的能量分配权值；其中，n个脑电信号的能量的权值之和为1。
在本优选实施例中，所述权值分布函数为正态分布函数，如可为：
需要说明的是，在计算得到各个脑电信号的能量的权值后，还需要进行归一化，保证n个脑电信号的能量的权值之和为1。
然后，对所述n个脑电信号的能量根据分配的权值进行加权求和，得到新的第j个时刻的脑电信号的能量。
即：
最后，依次对所述原始脑电序列信号的各个时刻的脑电信号的能量进行加权求和后，根据所有时刻的新的脑电信号的能量，生成待处理脑电序列信号。
在本发明实施例中，一般地，需要将所述待处理脑电序列信号进行再次切片，如切成时间长度为6秒的片段。
本优选实施例中，对脑电信号进行去低频直流信息，避免了这些低频直流信息与脑电波的频率发生重叠而影响提取的效果，而且对待处理脑电序列信号的各个时刻的脑电信号进行加权移动平均计算，可以避免去除低频直流信息时，造成信号失真，保证信号的真实度。
s102，从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波。
在本发明实施例中，所述的各个脑电波可包括delta波、theta波、alpha波、beta波、gamma波。其中，一般地，delta波的频率范围为0.5～3hz，theta波的频率范围为3～7hz、alpha波的频率范围为8～13hz、beta波的频率范围为14～17hz、gamma波的频率范围为34～50hz。
其中，delta波：深度睡眠脑波状态。
当人的大脑频率处于delta波时，为深度睡眠、无意识状态。人的睡眠品质好坏与delta波有非常直接的关系。delta波睡眠是一种很深沉的睡眠状态，如果在辗转难眠时自己召唤出近似delta波状态，就能很快地摆脱失眠而进入深沉睡眠。
theta波：深度放松、无压力的潜意识状态。
当人的大脑频率处于theta波时，人的意识中断，身体深沉放松，对于外界的信息呈现高度的受暗示状态，即被催眠状态。theta波对于触发深沉记忆、强化长期记忆等帮助极大，所以theta波被称为"通往记忆与学习的闸门"。
alpha波：学习与思考的最佳脑波状态。
当人的大脑频率处于alpha波时，人的意识清醒，但身体却是放松的，它提供意识与潜意识的“桥梁”。在这种状态下，身心能量耗费最少，相对的脑部获得的能量较高，运作就会更加快速、顺畅、敏锐。alpha波被认为是人们学习与思考的最佳脑波状态。
beta波：紧张、压力、脑疲劳时的脑波状态。
人们清醒时，大部分时间大脑频率处于beta波状态。随着beta波的增加，身体逐渐呈紧张状态，因而削减了体内免疫系统能力，此时人的能量消耗加剧，容易疲倦，若不充分休息，容易堆积压力。适当的beta波对注意力提升以及认知行为的发展有积极作用。
在本发明实施例中，在获得所述待处理脑电序列信号后，可根据每个脑电波的频率范围通过滤波(如卡尔曼滤波)、小波变换或者自回归模型等提取算法从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波。其中，提取时可仅用一个算法提取得到对应于各个脑电波的信号波，也可同时通过多个算法进行提取，再对不同算法提取得到的结果进行加权求和，得到最终的信号波。利用多个提取算法进行信号波的提取，可避免单一算法提取出现的误差偏大或稳定性不高的问题。
s103，计算对应于各个脑电波的信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，获取所述待处理脑电序列信号在时域空间的特征量。
在本发明实施例中，在获得与各个脑电波对应的信号波后，即可提取其在时域空间的特征。
本发明实施例中，可通过等电位线变化率算法提取各个信号波的等电位线能量变化率和等电位线能量密度变化率，从而获得各个信号波的时域特征。
具体地：
s1031，计算各个信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，计算得到每个信号波的等电位线能量参数变化率基线。
在本发明实施例中，所述能量参数至少包括以下其中之一：能量、能量密度。以所述能量参数为能量为例：
首先，对于一个信号波a(t)，其能量及能量密度的计算公式如下：
在本发明实施例中，等电位线能量变化率比率xlepsrate可通过大量数据参数自学习得到(如通过大量学习各个脑电波的能量变化率基线与能量最大值的比值，求得一个平均的比值作为等电位线能量变化率比率xlepsrate)，且针对delta、theta、alpha、beta、gamma波的等电位线能量变化率比率xlepsrate值都不一样。
然后，计算等电位线能量变化率基线xleps。
xleps＝max(p)*xlepsrate(3)
其中，max(p)为每个信号波的能量的最大值。
s1032，基于最小二乘法计算每个信号波在各个时刻的能量参数变化率。
具体地：
首先，构造关于能量的函数表达式p(t)＝at+b；其中，a为能量参数变化率，b为截距，t为时刻。
然后，计算使得所述函数表达式的残差平方和最小时的能量变化率，获得每个时刻的能量变化率。
在本发明实施例中，要想利用最小二乘算法求得所述函数表达式的参数a，b，就要使残差平方和最小，具体的求解参数a，b的过程如公式4至公式6所示：
s1033，统计所述每个信号波的能量参数变化率小于所述等电位线能量参数变化率基线的个数。
即统计a<xleps的个数xlcount。
s1034，根据统计得到的对应于每个信号波的个数与对应的信号波的总采样点数的比值，获取所述待处理脑电序列信号在时域空间的特征量。
在本发明实施例中，将统计得到的个数xlcount除以对应的信号波的总采样点数(或全部采样时刻)，得到比值xlcountp，即得到各个信号波的等电位线能量变化率特征。
在本发明实施例中，根据上述方法分别计算得到每个信号波(delta、theta、alpha、beta、gamma波)的等电位线能量变化率xlcountp1～xlcountp5，此后，再按照同样的方法计算每个脑电波(delta、theta、alpha、beta、gamma波)的等电位线能量密度变化率xlcounts1～xlcounts5，就获取得到待处理脑电序列信号在时域空间的特征量。
s104，根据所述各个脑电波的信号波及所述待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量。
具体地：
首先，根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表。
其中，二维图表的横坐标为时间，纵坐标为信号的幅值。
其次，在各个二维图表上覆盖与所述二维图表等大小的m*m的网格对应，并统计覆盖有信号的网格数；其中，m为大于1的整数。
其中，m为大于1的整数，且m的取值由信号的长度决定。
最后，根据覆盖有信号的网格数与所述网格的全部网格数，计算各个信号波及待处理脑电序列信号的相空间分布密度，获得待处理脑电序列信号在相空间的特征量。
即，相空间分布密度＝md/m2。其中，md为覆盖有信号的网格数。
在本发明实施例中，在得到各个信号波及所述待处理脑电序列信号的相空间分布密度后，就得到了所述待处理脑电序列信号在相空间的特征量。
s105，根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度。
在本发明实施例中，在得到所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量之后，将其输入到预先设置好的分类器中，就可以得到当前的脑电放松度。
在本发明实施例中，在得到脑电放松度后，就可以根据脑电放松度进行放松治疗，如可以根据脑电放松度进行放松引导内容的选择、标记和播放，能够精准地选取最适合用户的放松引导内容，搭配耳机，播放给用户；同时伴随着基于放松度的放松引导内容播放音量调制，帮助使用者放松身心，缓解焦虑抑郁，陶冶情操、改善个性弱点、消除心理行为障碍、保持心理和躯体健康。
本发明实施例提供的基于时域及相空间的脑电放松度识别方法，同时从时域空间、相空间两个角度提取各个脑电波的信号波的特征，得到待处理脑电序列序号的特征量并基于提取得到的特征量进行分类识别得到最终的脑电放松度。相比于单一角度的特征提取，对信号的评价方式更多样化，能更全面的体现信号的特性，避免了单一角度的特征提取容易导致的特征过于片面的问题而影响最终识别精度的问题。本发明可以大大提高识别分类的精度和准确度，为放松治疗提供准确的依据。此外，本发明实施例提供的特征提取方法计算简单、快速，对硬件要求低，因而能及时输出分类结果，便于实时的放松治疗。
优选地，所述步骤s105具体包括：
s1051，基于由不同学习算法训练得到的至少两个神经网络模型对输入的所有特征量进行分类，得到待处理脑电序列信号在各个神经网络模型下的分类。
s1052，将出现次数最多的分类设置为所述待处理脑电序列信号的分类。
s1053，根据所述分类与脑电放松度的对应关系，识别得到与所述待处理脑电序列信号对应的脑电放松度。
在本优选实施例中，所述的至少两个神经网络模型可采用levenberg-marquart算法、标准bp算法、增加动量项的bp算法、改进的ga-bp算法进行训练得到，其中，每个算法训练得到一个神经网络模型，即对应一个分类器。
训练时，将提取得到的各个特征量作为训练神经网络的输入样本x，将神念设备(或其他标准设备)同步采集得到的“放松度”作为金标准，也就是神经网络的输出y。(x，y)共同组成神经网络的训练样本对，再根据相应的训练算法进行神经网络训练。
在本优选实施例中，利用训练得到的神经网络模型，就可以作为分类器进行脑电放松度识别，识别得到脑电放松度的等级，如可以分别为弱、较弱、中、较强、强五级(即1～5级)。当前也可以分为更其他等级数，如4级，7级等，本发明不做具体限定。
在本优选实施例中，可将提取得到的所有特征量输入到各个分类器，每个分类器都将识别得到一个分类，在确定最终分类时，将出现次数最多的分类设置为所述待处理脑电序列信号的分类，根据所述分类与脑电放松度的对应关系，识别得到与所述待处理脑电序列信号对应的脑电放松度。
本优选实施例中，通过多个分类器同步识别分类的方法，并将出现最多的分类作为最终的分类，可避免单一分类器不稳定或出现偏差而导致的放松度识别结果不准确，进而影响放松治疗的效果。
优选地，在步骤s104之后，步骤s105之前，还可以包括：
s106，基于方差分析和f检验计算每个特征量的显著性指标，并选取显著性指标大于预设阈值的特征量。
具体地：
首先，根据所述特征量、由标准设备同步采集得到的脑电放松度及预先拟合的与每个特征量对应的线性拟合曲线进行方差分析，计算得到每个特征量的拟合值与期望的平方和及原值与拟合值的平方和。
在本发明实施例中，在准备阶段，在采集脑电切片信号的同时，还会利用标准设备(如神念设备)同步采集脑电放松度作为标准的脑电放松度。这样当采集一系列脑电切片信号(a1，a2，a3…an)后，就可以同时得到一系列标准脑电放松度(y1，y2….yn)。其中，在采集脑电切片信号时，可从所述脑电切片信号提取多个特征量，例如，对于脑电切片信号a1，其可提取特征量(x11，x11，x11….)，对于脑电切片信号a2，其可提取特征量(x21，x21，x21….)。这样，对于同一种特征量(例如，这里设x11，x21，…xn1是同一种特征量)，可拟合出其与脑电放松度的函数关系，即对坐标点(x11，y1)、(x21，y2)、…(xn1，yn)进行拟合，这里假设放松度与特征量的关系均是线性关系的，因此每种特征量对应的线性拟合曲线可表示为yi＝xβi+εi。
在选取特征量阶段，在计算得到特征量后，将当前的所述特征量及由标准设备同步采集得到的脑电放松度代入到与每个特征量对应的线性拟合曲线中进行方差分析，即可分别计算得到对应的拟合值与期望的平方和(ssr)，原值与拟合值的平方和(sse)。
然后，根据拟合值与期望的平方和、原值与拟合值的平方和及自由度，计算得到每个特征量的显著性指标。
其中，显著性指标f＝msr/mse，且msr＝ssr/1；mse＝sse/n-2；其中，1和n-2表示自由度。
最后，根据每个特征量的显著性指标，选取出显著性指标大于预设阈值的特征量。
在本发明实施中，设定显著性指标的阈值fmin。选取所述特征量中，f>fmin的特征量作为需要输入的特征量。
本优选实施例中，结合方差分析与f检验对所述特征量进行选取，选出较为显著的特征，从而可以提高分类的精确度和速度。
优选地，在步骤s104之后，步骤s105之前，还可以包括：
s107，基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行降维处理，获得降维后的特征量。
具体地：
s1071，将所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量设置为输入样本空间中的特征量，并对所述输入样本空间进行数据标准化处理。
具体地，将待处理脑电序列信号在时域空间的特征量及在相空间的特征量设置为输入样本空间x中的元素。对该样本空间x进行数据标准化处理具体为：
其中：
其中，x′ij是标准化后的新数据；mj、sj分别表示原始数据某一列的算术平均值和标准(偏)差。
s1072，根据数据标准化处理后的所述输入样本空间得到协方差矩阵。
其中，协方差矩阵d＝xtx，即：
其中：
s1073，计算所述协方差矩阵的特征根和与每个特征根对应的特征向量；其中，所述特征根的数量为p个，且所述的p个特征根按照大小排序。
其中，dp＝pλ(12)
当只考虑第j个特征值时，有dpj＝pjλj，即求解|d-λji|＝0。依次解出各个λ，并使其按大小顺序排列，即λ1≥λ2≥…,≥λp≥0；然后可求出各特征值对应的特征向量p，进而特征方程求解完成。
s1074，获取所述的p个特征根中，贡献率之和大于预定阈值的前m个特征根。
其中，每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和。
首先，计算单个主成分的贡献率并进行累计，根据累计贡献率确定主成分的个数m，从而确定所需要选取的主成分。贡献率的计算公式如公式15所述。累计贡献率即前m个贡献率的累积和，如公式16所示。所述阈值dmax一般取在85％～95％之间。根据上一步骤中的特征根排序可知，λ1≥λ2≥…,≥λp≥0，从前往后(也是从大到小)依次对特征根进行累加，当累计贡献率
s1075，根据与所述的前m个特征根对应的特征向量及所述输入样本空间，得到主成分得分矩阵。
其中，所述主成分得分矩阵中的特征量为所述降维后的特征量。
其中，所述主成分得分矩阵
其中，该主成分得分矩阵t中的各个元素即为经过降维后的特征量。
需要说明的是，在本发明实施例中，还可计算主成分的载荷，其中，所述主成分载荷主要反映主成分得分与原变量xj的关联程度，计算公式为：
在本优选实施例中，在采用主成分分析法筛选出所获得的待处理脑电序列信号在时域空间和相空间的特征量中较为重要的特征量后，即可得到降维后的待处理脑电序列信号的特征量。通过对待处理脑电序列信号在时域空间和相空间的特征量进行降维处理，使得后续进行脑电放松度识别时的数据计算量减小，从而提高脑电放松度识别的速度，实现脑电放松度的实时识别。
此外，需要说明的是，上述的基于方差分析和f检验的降维方法与本优选实施例的主成分分析法可同时使用，本发明不做具体限定。
请参阅图4，本发明还提供一种基于时域及相空间的脑电放松度识别装置100，包括：
信号切片单元10，用于对接收的原始脑电信号进行切片，生成具有指定时间长度的待处理脑电序列信号；
信号提取单元20，用于从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波；
第一特征提取单元30，用于计算对应于各个脑电波的信号波的能量参数，并基于各个信号波的能量参数的最大值及预设的每个脑电波的等电位线能量参数变化率比率，获取所述待处理脑电序列信号在时域空间的特征量；
第二特征提取单元40，用于根据所述各个脑电波的信号波及待处理脑电序列信号，形成对应的包含有信号的二维图表，并通过统计各个二维图表上的信号覆盖密度，获得所述待处理脑电序列信号在相空间的特征量；以及，
脑电放松度识别单元50，用于根据所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行分类识别，得到脑电放松度。
优选地，所述第一特征提取单元30具体包括：
变化率基线计算模块，用于计算各个信号波的能量参数，并基于各个信号波的能量参数的最大值及每个脑电波的等电位线能量参数变化率比率，计算得到每个信号波的等电位线能量参数变化率基线；
能量参数变化率计算模块，用于基于最小二乘法计算每个信号波在各个时刻的能量参数变化率；
统计模块，用于统计所述每个信号波的能量参数变化率小于所述等电位线能量参数变化率基线的个数；以及，
比值计算模块，用于根据统计得到的对应于每个信号波的个数与对应的信号波的总采样点数的比值，获取所述待处理脑电序列信号在时域空间的特征量。
优选地，所述第二特征提取单元40具体包括：
二维图表生成模块，用于根据所述各个脑电波的信号波及所述待处理脑电序列信号形成对应的包含有信号的二维图表；
网格数统计模块，用于在各个二维图表上覆盖与所述二维图表等大小的m*m的网格对应，并统计覆盖有信号的网格数；其中，m为大于1的整数；以及，
特征量计算模块，用于根据覆盖有信号的网格数与所述网格的全部网格数的比值，生成各个信号波及待处理脑电序列信号的相空间分布密度，获得待处理脑电序列信号在相空间的特征量。
优选地，所述基于时域及相空间的脑电放松度识别装置还包括：
特征选取单元，用于基于方差分析和f检验计算每个特征量的显著性指标，并选取显著性指标大于预设阈值的特征量。
优选地，所述特征选取单元具体包括：
方差分析模块，用于根据所述特征量、由标准设备同步采集得到的脑电放松度及预先拟合的与每个特征量对应的线性拟合曲线进行方差分析，计算得到每个特征量的拟合值与期望的平方和及原值与拟合值的平方和；
显著性指标计算模块，用于根据拟合值与期望的平方和、原值与拟合值的平方和及自由度，计算得到每个特征量的显著性指标；以及，
选取模块，用于根据每个特征量的显著性指标，选取出显著性指标大于预设阈值的特征量。
优选地，所述基于时域及相空间的脑电放松度识别装置还包括：
特征降维单元，用于基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量进行降维处理，获得降维后的特征量；
其中，所述特征降维单元具体包括：
标准化处理模块，用于将所述待处理脑电序列信号在时域空间的特征量及在相空间的特征量设置为输入样本空间中的特征量，并对所述输入样本空间进行数据标准化处理；
协方差矩阵计算模块，用于根据数据标准化处理后的所述输入样本空间得到协方差矩阵；
特征计算模块，用于计算所述协方差矩阵的特征根和与每个特征根对应的特征向量；其中，所述特征根的数量为p个，且所述的p个特征根按照大小排序；
筛选模块，用于获取所述的p个特征根中，贡献率之和大于预定阈值的前m个特征根；其中，每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和；以及，
降维特征量获得模块，用于根据与所述的前m个特征根对应的特征向量及所述输入样本空间，得到主成分得分矩阵；其中，所述主成分得分矩阵中的特征量为所述降维后的特征量。
优选地，所述脑电放松度识别单元50具体包括：
分类模块，用于基于由不同学习算法训练得到的至少两个神经网络模型对输入的所有特征量进行分类，得到待处理脑电序列信号在各个神经网络模型下的分类；
分类统计模块，用于将出现次数最多的分类设置为所述待处理脑电序列信号的分类；以及，
放松度识别模块，用于根据所述分类与脑电放松度的对应关系，识别得到与所述待处理脑电序列信号对应的脑电放松度。
以上所揭露的仅为本发明一种较佳实施例而已，当然不能以此来限定本发明之权利范围，本领域普通技术人员可以理解实现上述实施例的全部或部分流程，并依本发明权利要求所作的等同变化，仍属于发明所涵盖的范围。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(randomaccessmemory，ram)等。
- 还没有人留言评论。精彩留言会获得点赞！