CN108564128A - 一种融合数据结构信息的脑电信号在线识别方法 - Google Patents

Abstract

本发明涉及一种融合数据结构信息的脑电信号在线识别方法，包括S1：利用由少量标记脑电样本组成的小训练集建立基于在线序贯极限学习机(online sequential extreme learning machine,OS‑ELM)算法的分类模型，作为半监督学习中的初始分类模型；S2：利用在线模糊聚类方法建立结构学习模型，基于标记脑电样本的先验信息估计在线采集脑电样本批量增加后数据分布的全局结构；S3：利用分类模型对在线采集的脑电样本进行标记，并采用批量学习的模式，基于由结构学习模型估计的结构信息，选择一批满足一定条件的在线采集脑电样本加入到训练集中，利用更新后的训练集重新训练分类模型；S4：利用更新后的分类模型对采集的脑电信号进行在线识别。

Description

一种融合数据结构信息的脑电信号在线识别方法

技术领域

本发明属于脑电信号处理技术领域，具体涉及一种融合数据结构信息的脑电信号在线识别方法。

背景技术

脑－机接口(brain-computer interface,BCI)作为一种基于大脑功能的电生理测量，实现人与外界信息交流和控制的全新通讯系统，在医疗、康复、航空、军事、认知科学、教育、生活娱乐等多个领域具有极高的应用价值和广阔的应用前景，成为了目前研究的热点之一。其中，基于头皮脑电信号的非侵入式BCI是当今BCI研究中的主流模式。

研究者们对基于头皮脑电信号的BCI已经开展了广泛的研究。但发明人在实现本发明的过程中，发现在BCI分类识别这一关键技术研究中存在以下的不足：

1、采用监督学习的方法来实现模式分类是主流模式，但监督方法需要大量的标记脑电样本才能构建有效的分类器，并且分类器在工作期间是固化的，不具有在线适应性。

2、目前有少数采用半监督分类方法的研究，可以利用未标记脑电样本来训练和更新分类器，但现有方法采用逐一将样本用于分类器学习的方式，存在学习的盲目性、误标记的累积效应、泛化能力较低等影响信号在线识别性能的问题。

发明内容

有鉴于此，本发明的目的在于提供一种融合数据结构信息的脑电信号在线识别方法，以解决现有技术存在的上述问题。

为达到上述目的，本发明提供如下技术方案：

一种融合数据结构信息的脑电信号在线识别方法，包括：

S1：利用由少量标记脑电样本组成的小训练集建立基于OS-ELM算法的分类模型，作为半监督学习中的初始分类模型；

S2：利用在线模糊聚类方法建立结构学习模型，基于标记脑电样本的先验信息估计在线采集脑电样本批量增加后数据分布的全局结构；

S3：利用分类模型对在线采集的脑电样本进行标记，并采用批量学习的模式，基于由结构学习模型估计的结构信息，选择一批满足一定条件的在线采集脑电样本加入到训练集中，利用更新后的训练集重新训练分类模型；

S4：利用更新后的分类模型对采集的脑电信号进行在线识别。

进一步，在步骤S1中，包括：

S11：由标记脑电样本组成的初始训练集样本个数N₀，样本特征维数m，包含K个隐层节点的初始隐含层输出矩阵H₀：

其中，a_j,j＝1,2,…,K为隐层节点的输入权参数，b_j,j＝1,2,…,K为隐层节点的偏移参数，N₀≥K，a_j和b_j这两组参数随机生成，为Sigmoidal激活函数；

计算初始输出权重：

式中，为目标向量，

k＝0，k为加入训练集的脑电样本批次；

S12：利用初始分类器对在线采集未标记脑电样本进行分类识别(标记)：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β₀

式中x代表在线采集的未标记脑电样本。

进一步，在步骤S2中，对于一批在线采集未标记脑电样本X＝{x₁,x₂,…,x_n},x_j∈R^m，每一类的聚类中心为V＝[v_i]_2×m，建立如下结构学习模型目标函数：

式中，u_ij表示第j个未标记样本属于第i类的隶属度，且满足u′_ij表示第j个标记样本属于第i类的隶属度，且满足d_ij表示第j个未标记样本到第i类聚类中心v_i的距离；d′_ij表示第j个标记样本到第i类聚类中心v_i的距离；p为模糊参数，通常取值为2；N为标记样本的个数(N≥N₀)。

通过拉格朗日乘子法最小化目标函数，可以获得数据结构属性，实现对样本的划分，具体步骤为：

1)初始化模糊聚类矩阵且使其满足约束条件；计算初始聚类中心：

2)进行迭代计算：

标记样本隶属度采用如下公式进行更新：

未标记样本隶属度采用如下公式进行更新：

聚类中心更新公式为：

式中，n'为标记样本的数量，n为新增未标记样本的数量；

3)设定计算精度ε，如果J^r-J^r-1<ε,则停止迭代；否则，重复上一步。J^r-1表示第r-1次迭代计算的目标函数值J，J^r表示第r次迭代计算的目标函数值J。

进一步，在步骤S3中，根据由结构学习模型估计的样本结构信息并依据置信度准则，从第k批由分类器在线标记过的脑电样本中选择出N_k个隶属度大于设定阈值的脑电样本，组成集合计算隐含层输出矩阵H_k：

目标向量为：

更新M_k和输出权重β_k：

进一步，在步骤S4中，利用更新后的分类器对在线采集脑电信号进行分类识别(标记)：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β_k

式中，x代表在线采集的未标记脑电样本；

同时，重复步骤S2、S3和S4的在线学习。

本发明的有益效果在于：本发明针对现有BCI研究中脑电信号在线模式分类存在的问题，提出一种基于数据全局结构信息的脑电信号在线识别方法，解决分类器训练的盲目性问题，降低样本错误标记的风险，提高算法的泛化能力和在线学习的安全性。

附图说明

为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：

图1为本发明实施例所述的融合数据结构信息的脑电信号在线识别方法的逻辑框图。

具体实施方式

下面将结合附图，对本发明的优选实施例进行详细的描述。

如图1所示，本发明为了解决所述技术问题，采用了如下的技术方案：首先，采用在线序贯极限学习机(online sequential extreme learning machine,OS-ELM)算法构建能满足在线训练速度要求的半监督分类模型；其次，利用模糊聚类方法对由标记和未标记脑电样本共同构建的信号空间进行模糊划分，允许未标记样本同时属于多个类，对应不同的软隶属度，建立能在线学习的结构学习模型；最后，以分类模型的在线更新为接口，实现结构学习模型与半监督分类模型的融合，构建基于数据结构信息的在线自适应分类模型，提高分类模型的泛化能力和在线学习的安全性。

具体的，为了实现以上目的，本发明方法包括以下步骤：

101、利用由少量标记脑电样本组成的小训练集建立基于OS-ELM算法的分类模型，作为半监督学习中的初始分类模型。

102、利用在线模糊聚类方法建立结构学习模型，基于标记脑电样本的先验信息估计在线采集脑电样本批量增加后数据分布的全局结构。

103、利用分类模型对在线采集的脑电样本进行标记，并采用批量学习的模式，基于由结构学习模型估计的结构信息，选择一批满足一定条件的在线采集脑电样本加入到训练集中，利用更新后的训练集重新训练分类模型。

104、利用更新后的分类模型对采集的脑电信号进行在线识别。

本发明所述脑电信号在线识别方法包括了OS-ELM初始分类模型的建立、脑电数据结构信息的估计、分类模型的在线训练和更新、在线采集脑电信号的识别这几个关键步骤，下面对各个步骤进行详细的说明。

1、建立OS-ELM初始分类模型。设由标记脑电样本组成的初始训练集为N₀为样本个数，m为样本特征维数。计算包含K个隐层节点的初始隐含层输出矩阵H₀：

其中，a_j,j＝1,2,…,K为隐层节点的输入权参数，b_j,j＝1,2,…,K为隐层节点的偏移参数，N₀≥K，a_j和b_j这两组参数随机生成。为Sigmoidal激活函数。

计算初始输出权重：

式中，为目标向量，

设k＝0，k为加入训练集的脑电样本批次。

2、利用初始分类器对在线采集未标记脑电样本进行分类识别(标记)：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β₀

式中x代表在线采集的未标记脑电样本。

3、建立脑电数据结构学习模型，估计数据结构信息。对于一批在线采集未标记样本X＝{x₁,x₂,…,x_n},x_j∈R^m，每一类的聚类中心为V＝[v_i]_2×m，建立如下结构学习模型目标函数：

式中，u_ij表示第j个未标记样本属于第i类的隶属度，且满足u′_ij表示第j个标记样本属于第i类的隶属度，且满足d_ij表示第j个未标记样本到第i类聚类中心v_i的距离；d′_ij表示第j个标记样本到第i类聚类中心v_i的距离；p为模糊参数，通常取值为2；N为标记样本的个数(N≥N₀)。

通过拉格朗日乘子法最小化目标函数，可以获得数据结构属性，实现对样本的划分，具体步骤为：

1)初始化模糊聚类矩阵且使其满足约束条件。计算初始聚类中心：

2)进行迭代计算。

标记样本隶属度采用如下公式进行更新：

未标记样本隶属度采用如下公式进行更新：

聚类中心更新公式为：

式中，n'为标记样本的数量，n为新增未标记样本的数量。

3)设定计算精度ε，如果J^r-J^r-1<ε,则停止迭代；否则，重复上一步。J^r-1表示第r-1次迭代计算的目标函数值J，J^r表示第r次迭代计算的目标函数值J。

4、采用批量学习模式对分类模型进行在线训练和更新。根据由结构学习模型估计的样本结构信息并依据置信度准则，从第k批由分类器在线标记过的脑电样本中选择出N_k个隶属度大于设定阈值的脑电样本，组成集合计算隐含层输出矩阵H_k：

目标向量为：

更新M_k和输出权重β_k：

5、利用更新后的分类器对在线采集脑电信号进行分类识别(标记)：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β_k

式中，x代表在线采集的未标记脑电样本。

同时，重复第3步、第4步和第5步的在线学习步骤。

最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。

Claims (5)

1.一种融合数据结构信息的脑电信号在线识别方法，其特征在于：包括：

S1：利用由少量标记脑电样本组成的小训练集建立基于OS-ELM算法的分类模型，作为半监督学习中的初始分类模型；

S2：利用在线模糊聚类方法建立结构学习模型，基于标记脑电样本的先验信息估计在线采集脑电样本批量增加后数据分布的全局结构；

S3：利用分类模型对在线采集的脑电样本进行标记，并采用批量学习的模式，基于由结构学习模型估计的结构信息，选择一批满足一定条件的在线采集脑电样本加入到训练集中，利用更新后的训练集重新训练分类模型；

S4：利用更新后的分类模型对采集的脑电信号进行在线识别。

2.根据权利要求1所述的融合数据结构信息的脑电信号在线识别方法，其特征在于：在步骤S1中，包括：

S11：由标记脑电样本组成的初始训练集样本个数N₀，样本特征维数m，包含K个隐层节点的初始隐含层输出矩阵H₀：

其中，a_j,j＝1,2,…,K为隐层节点的输入权参数，b_j,j＝1,2,…,K为隐层节点的偏移参数，N₀≥K，a_j和b_j这两组参数随机生成，为Sigmoidal激活函数；

计算初始输出权重：

式中，为目标向量，

k＝0，k为加入训练集的脑电样本批次；

S12：利用初始分类器对在线采集未标记脑电样本进行分类识别：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β₀

式中x代表在线采集的未标记脑电样本。

3.根据权利要求2所述的融合数据结构信息的脑电信号在线识别方法，其特征在于：在步骤S2中，对于一批在线采集未标记脑电样本X＝{x₁,x₂,…,x_n},x_j∈R^m，每一类的聚类中心为V＝[v_i]_2×m，建立如下结构学习模型目标函数：

式中，u_ij表示第j个未标记样本属于第i类的隶属度，且满足u'_ij表示第j个标记样本属于第i类的隶属度，且满足u'_ij≥u'_kj,d_ij表示第j个未标记样本到第i类聚类中心v_i的距离；d'_ij表示第j个标记样本到第i类聚类中心v_i的距离；p为模糊参数，通常取值为2；N为标记样本的个数(N≥N₀)。

通过拉格朗日乘子法最小化目标函数，可以获得数据结构属性，实现对样本的划分，具体步骤为：

1)初始化模糊聚类矩阵且使其满足约束条件；计算初始聚类中心：

2)进行迭代计算：

标记样本隶属度采用如下公式进行更新：

未标记样本隶属度采用如下公式进行更新：

聚类中心更新公式为：

式中，n'为标记样本的数量，n为新增未标记样本的数量；

3)设定计算精度ε，如果J^r-J^r-1<ε,则停止迭代；否则，重复上一步。J^r-1表示第r-1次迭代计算的目标函数值J，J^r表示第r次迭代计算的目标函数值J。

4.根据权利要求3所述的融合数据结构信息的脑电信号在线识别方法，其特征在于：在步骤S3中，根据由结构学习模型估计的样本结构信息并依据置信度准则，从第k批由分类器在线标记过的脑电样本中选择出N_k个隶属度大于设定阈值的脑电样本，组成集合计算隐含层输出矩阵H_k：

目标向量为：

更新M_k和输出权重β_k：

5.根据权利要求4所述的融合数据结构信息的脑电信号在线识别方法，其特征在于：在步骤S4中，利用更新后的分类器对在线采集脑电信号进行分类识别：

f_K(x)＝[G(a₁,x,b₁),…,G(a_K,x,b_K)]β_k

式中，x代表在线采集的未标记脑电样本；

同时，重复步骤S2、S3和S4的在线学习。