基于子空间正交向量的峰电位检测方法技术

本发明专利技术公开了一种基于子空间正交向量的峰电位检测方法，主要用于解决现有同类技术计算复杂度高，且对高噪脑电信号检测准确率过低的问题。其实现步骤为：1)根据输入信号构建冗余字典D；2)将字典D分解到l个子空间中；3)将输入信号投影到各子空间得到l个系数矩阵；4)设定阈值T，取出各系数矩阵中大于T的值构成重构系数向量；5)计算D的伪逆与系数向量的乘积得到检测结果。本发明专利技术计算复杂度低，检测准确率高，实验表明，本发明专利技术对信噪比为+30dB至‑11.6dB的低噪脑电信号峰电位检测准确率为100％，对信噪比为‑11.6dB至‑16.2的高噪脑电信号峰电位的检测准确率为80％以上，可用于对胞外单电极峰电位的分类。

【技术实现步骤摘要】

本专利技术属于信号处理
，特别涉及对脑电信号的峰电位的检测，可用于对胞外单电极峰电位的分类。
技术介绍
大脑是生命活动最复杂和最重要的器官之一，是神经系统最重要的组成部分。神经元是大脑信号加工传递和信息处理的基本单元，大脑中传递的信息都通过神经元产生的电信号表现出来，因而神经电信号是研究大脑信息传递的重要方面，其中包括单个神经元细胞外动作电位，又称作峰电位。而胞外单电极记录作为一种非常成熟的神经信号采集的传统手段，在目前的神经电生理学领域仍然被广泛使用。由于胞外单电极记录所获得的信号含有比较高的背景噪音，而且依赖于不同的电极尖端阻抗和形状，单根电极可能记录到多个神经元的信号，因此从单根电极电压信号提取和分类动作电位一直是研究点。通过研究峰电位，知晓峰电位序列，可以了解神经元编码机制和神经元之间的相互连接关系，能够为大脑神经机制的研究提供重要信息。大脑神经元胞外单细胞峰电位的检测与分类，是许多脑功能研究的前提。但由于高强度的背景噪声的存在，以及同一脑区内的神经元发放的峰电位在形态和尺度上很相似，通常很难将单个神经元发放的峰电位有效地提取出来。然而，使用神经元峰电位检测与分类算法，可以大大地提高对于单个神经元发放的峰电位的提取效果，从而更好的满足神经生理学实验的要求。目前为止，从峰电信号中检测并提取出峰电位的各种算法中，应用最广泛的仍是幅度阈值检测法及其改进算法。该算法在对峰电位检测之前，首先对其进行降噪，可以不同程度的降低信号中噪声对检测效果的影响。阈值法是设定一个阈值，当峰电信号超出此阈值时，认为是有峰电位出现。关于阈值的设定，起初通常根据实验经验设置定阈值，改进的阈值设定方法通常根据采集信号的标准差或者灰度值等。其中阈值的选定对检测结果也有很大的影响。阈值设置过大，可能无法将幅值较小的峰电位检测出，形成漏检；阈值设置过小，可能将幅值较大的噪声当做峰电位检出，形成误检。通过输入脑电信号的标准峰波构建冗余字典的检测方法目前也常常被采用，由上海
交通大学提出的专利申请“基于稀疏重表示的脑电信号检测方法”，专利申请号为CN201110051448.2，公开了一种基于稀疏重表示的用于检测驾驶员是否疲劳驾驶的脑电信号检测方法。该方法主要采用小波变换来检测与驾驶员的疲劳驾驶状态相关的脑电信号的频率特征信息，构建一个模板矩阵，对于任意的外来脑电信号序列，通过求解模板矩阵用来构建外来序列的稀疏解来实现对驾驶员脑电信号的检测。然而该方法需要对输入脑电信号进行较复杂的预处理工作，且对脑电信号做小波变换后，利用压缩感知理论，转化为求解一范数凸优化问题，该方法计算复杂度高，检测效率低下。综上所述，仅仅根据输入脑电信号的几何特征设定相应阈值的阈值检测法虽然实现简单且易于操作，但该算法不具有良好的鲁棒性，容易形成漏检或者误检；而采用直接将输入信号投影到冗余的即稀疏重表示的字典中，然后转化为优化问题的求解方式虽然在一定程度上提高了检测成功率，但是计算复杂度太高，导致检测效率低下。
技术实现思路
本专利技术的目的在于针对上述已有技术中的不足，提出一种基于子空间正交向量的脑电信号峰电位检测方法，以降低计算复杂度，提高检测效率。本专利技术的技术方案是：根据输入脑电信号的特征构建一个不完备的冗余字典，再通过字典构建子空间，其中子空间中的向量是标准正交的，然后通过将输入脑电信号投影到各个子空间中，设定阈值，将保留的系数用反变换的方法重构信号，完成对峰电位的检测。其实现步骤包括如下：(1)根据输入信号x∈R1×m的标准峰波，构造不完备冗余字典D∈R(m-l+1)×m，其中m表示输入脑电信号的长度，l表示标准峰波的长度；(2)通过步骤(1)中构建的不完备冗余字典D填充子空间V1，V2，…，Vj，…Vl，其中Vj∈Rn×m表示第j个子空间，1≤j≤l，其中(3)将输入脑电信号分别投影到各个子空间中，得到l个系数向量：C1,C2,...,Ct,...Cl，Ct∈Rn×1表示第t个系数向量，1≤t≤l；(4)设定阈值T为1.25-1.45，保留Ct中大于阈值T的值并将小于T的值全部置0，将C1,C2,...,Ct,...Cl依次拼接得到重构系数向量L∈R(m-l+1)×1。(5)将字典D的伪逆与重构系数向量L相乘得到重构信号此时中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即为峰电位。本专利技术与现有技术相比具有以下优点：1.计算复杂度低本专利技术将脑电信号投影到子空间中，且子空间中的原子是相互正交的，通过设定阈值再反变换得到检测结果，极大的降低了计算复杂度。2.鲁棒性好实验结果表明，本专利技术对噪声方差为0.05—0.55，信噪比为+30dB—-11.6dB的低噪脑电信号峰电位检测准确率可达100％，同时对噪声方差为0.55—0.7，信噪比为-11.6dB—-16.2的高噪脑电信号峰电位的检测准确率到达80％以上，表现出良好的鲁棒性。附图说明图1为本专利技术的实现总流程图；图2为本专利技术中冗余字典构建示意图；图3为本专利技术中子空间构建示意图。具体实施方式下面结合附图和实例对本专利技术进行详细说明。参照图1，本专利技术的实现步骤如下：步骤1，根据输入信号构建冗余字典D。参考图2，本步骤的具体实现如下：(1a)设矩阵A(a1,a2...,as,...am)T∈R(m-l+1)×m中所有元素为0，1≤s≤m，将A中的第s个行向量as的第i个到i+l个元素置为标准峰波，对as进行归一化处理，使其模值||as||＝1；(1b)令as＝di，其中di表示字典D中的第i个原子，1≤i≤m，得到字典D(d1,d2...,di,...dm)T，di的模值为||di||＝1，字典D中的第i行di是由标准峰波与m-l个0元素构成，标准峰波起始于di的第i个元素。步骤2，构建子空间V1，V2，…，Vj，…Vl。参考图3，本步骤的具体实现如下：(2a)由标准峰波与m-l个0元素构成第一子空间V1的原子，其中标准峰波起始于原
子v11的第1个元素、原子v12的第1+l个元素、原子v1k的第k+l个元素和v1n的第n+l个元素，得到第1个子空间V1(v11,v12,...,v1k,...v1n)T∈Rn×m，其中1≤k≤n，(2b)由标准峰波与m-l个0元素构成第二子空间V2的原子，其中标准峰波起始于原子v21的第2个元素、原子v22的第2+l个元素、原子v2k的第2+(k-1)×l个元素和v2n的第2+(n-1)×l个元素，得到第2个子空间V2(v21,v22,...,v2k,...v2n)T∈Rn×m；(2c)由标准峰波与m-l个0元素构成第j个子空间Vj的原子，1≤j≤l，标准峰波起始于原子vj1的第j个元素、原子vj2的第j+l个元素、原子vjk的第j+(k-1)×l个元素和vjn的第j+(n-1)×l个元素，得到第j个子空间Vj(vj1,vj2,...,vjk,...vjn)T∈Rn×m；(2d)由标准峰波与m-l个0元素构成第l个子空间Vl的原子,标准峰波起始于原子vl1的第l个元素、原子vl2的第2×l个元素、原子vlk的第k×l个元素和vln的第n×l个元素，得到第l个子空间Vl(vl1,vl2,...,vlk,...vln)T∈Rn×m。步骤3，将输入脑本文档来自技高网...

【技术保护点】
一种基于子空间正交向量的峰电位检测方法，包括：(1)根据输入信号x∈R1×m的标准峰波，构造不完备冗余字典D∈R(m‑l+1)×m，其中m表示输入脑电信号的长度，l表示标准峰波的长度；(2)通过步骤(1)中构建的不完备冗余字典D填充子空间V1，V2，…，Vj，…Vl，其中Vj∈Rn×m表示第j个子空间，1≤j≤l，其中(3)将输入脑电信号分别投影到各个子空间中，得到l个系数向量：C1,C2,...,Ct,...Cl，Ct∈Rn×1表示第t个系数向量，1≤t≤l；(4)设定阈值T为1.25‑1.45，保留Ct中大于阈值T的值并将小于T的值全部置0，将C1,C2,...,Ct,...Cl依次拼接得到重构系数向量L∈R(m‑l+1)×1。(5)将字典D的伪逆与重构系数向量L相乘得到重构信号此时中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即为峰电位。

【技术特征摘要】
1.一种基于子空间正交向量的峰电位检测方法，包括：(1)根据输入信号x∈R1×m的标准峰波，构造不完备冗余字典D∈R(m-l+1)×m，其中m表示输入脑电信号的长度，l表示标准峰波的长度；(2)通过步骤(1)中构建的不完备冗余字典D填充子空间V1，V2，…，Vj，…Vl，其中Vj∈Rn×m表示第j个子空间，1≤j≤l，其中(3)将输入脑电信号分别投影到各个子空间中，得到l个系数向量：C1,C2,...,Ct,...Cl，Ct∈Rn×1表示第t个系数向量，1≤t≤l；(4)设定阈值T为1.25-1.45，保留Ct中大于阈值T的值并将小于T的值全部置0，将C1,C2,...,Ct,...Cl依次拼接得到重构系数向量L∈R(m-l+1)×1。(5)将字典D的伪逆与重构系数向量L相乘得到重构信号此时中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即...

【专利技术属性】
技术研发人员：谢雪梅，柳林，翁昕，焦阳，石光明，
申请(专利权)人：西安电子科技大学，西安中电科西电科大雷达技术协同创新研究院有限公司，
类型：发明
国别省市：陕西;61