基于子空间正交向量的峰电位检测方法专利检索-信噪比信号处理专利检索查询-专利查询网

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

阅读：0发布：2024-02-22

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

权利要求

1. 一种基于子空间正交向量的峰电位检测方法，包括： (1) 根据输入信号X e Rixm的标准峰波，构造不完备冗余字典D E ，其中m表示输入脑电信号的长度，1表示标准峰波的长度； (2) 通过步骤（1)中构建的不完备冗余字典D填充子空间Vi，V2，…，Vj，…Vi，其中Vj e RDXm 表示第j个子空间，其中^。 (3) 将输入脑电信号分别投影到各个子空间中，得至Ijl个系数向量:Ci，C2,...，Ct，...打， CtERnX嗦示第t个系数向量，; (4) 设定阔值T为1.25-1.45,保留Ct中大于阔值T的值并将小于T的值全部置0,将Cl, C2,...，Ct，.. .Cl依次拼接得到重构系数向量LeRfm-Wxi。 (5) 将字典D的伪逆与重构系数向量L相乘得到重构信号X，此时;^中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即为峰电位。
2. 根据权利要求1所述的方法，其中步骤（1)中根据输入信号X的标准峰波构造字典D，按如下步骤进行： (la) 设矩阵A(ai，a2. . .，as, . . .anOTeRfm-i+Ux®中所有元素为0, 将A中的第S个行向量as的第i个到i+1个元素置为标准峰波，对as进行归一化处理，使其模值Mas II =1; (lb) 令as = di，其中di表示字典D中的第i个原子，得到字典D(di，d2. . .，di， ...cU)T，di 的模值为 Mdi||=l。
3. 根据权利要求1所述的方法，其中步骤(2)中填充子空间Vi，V2,…，ν^，···νι，是先将子空间Vj(vji，vj2,. . .，vjk,. . .V化)Τ中的原子Vjk初始化为0,再将V化中的第j+(k-l) XI菌Jj+kX 1元素置为标准峰波，1，Vj中的每个原子都是标准正交的。

说明书全文

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

技术领域

[0001] 本发明属于信号处理技术领域，特别涉及对脑电信号的峰电位的检测，可用于对胞外单电极峰电位的分类。

背景技术

[0002] 大脑是生命活动最复杂和最重要的器官之一，是神经系统最重要的组成部分。神经元是大脑信号加工传递和信息处理的基本单元，大脑中传递的信息都通过神经元产生的电信号表现出来，因而神经电信号是研究大脑信息传递的重要方面，其中包括单个神经元细胞外动作电位，又称作峰电位。而胞外单电极记录作为一种非常成熟的神经信号采集的传统手段，在目前的神经电生理学领域仍然被广泛使用。由于胞外单电极记录所获得的信号含有比较高的背景噪音，而且依赖于不同的电极尖端阻抗和形状，单根电极可能记录到多个神经元的信号，因此从单根电极电压信号提取和分类动作电位一直是研究点。通过研究峰电位，知晓峰电位序列，可以了解神经元编码机制和神经元之间的相互连接关系，能够为大脑神经机制的研究提供重要信息。大脑神经元胞外单细胞峰电位的检测与分类，是许多脑功能研究的前提。但由于高强度的背景噪声的存在，以及同一脑区内的神经元发放的峰电位在形态和尺度上很相似，通常很难将单个神经元发放的峰电位有效地提取出来。然而，使用神经元峰电位检测与分类算法，可以大大地提高对于单个神经元发放的峰电位的提取效果，从而更好的满足神经生理学实验的要求。

[0003] 目前为止，从峰电信号中检测并提取出峰电位的各种算法中，应用最广泛的仍是幅度阈值检测法及其改进算法。该算法在对峰电位检测之前，首先对其进行降噪，可以不同程度的降低信号中噪声对检测效果的影响。阈值法是设定一个阈值，当峰电信号超出此阈值时，认为是有峰电位出现。关于阈值的设定，起初通常根据实验经验设置定阈值，改进的阈值设定方法通常根据采集信号的标准差或者灰度值等。其中阈值的选定对检测结果也有很大的影响。阈值设置过大，可能无法将幅值较小的峰电位检测出，形成漏检；阈值设置过小，可能将幅值较大的噪声当做峰电位检出，形成误检。

[0004] 通过输入脑电信号的标准峰波构建冗余字典的检测方法目前也常常被采用，由上海交通大学提出的专利申请"基于稀疏重表示的脑电信号检测方法"，专利申请号为 CN201110051448.2,公开了一种基于稀疏重表示的用于检测驾驶员是否疲劳驾驶的脑电信号检测方法。该方法主要采用小波变换来检测与驾驶员的疲劳驾驶状态相关的脑电信号的频率特征信息，构建一个模板矩阵，对于任意的外来脑电信号序列，通过求解模板矩阵用来构建外来序列的稀疏解来实现对驾驶员脑电信号的检测。然而该方法需要对输入脑电信号进行较复杂的预处理工作，且对脑电信号做小波变换后，利用压缩感知理论，转化为求解一范数凸优化问题，该方法计算复杂度高，检测效率低下。

[0005] 综上所述，仅仅根据输入脑电信号的几何特征设定相应阈值的阈值检测法虽然实现简单且易于操作，但该算法不具有良好的鲁棒性，容易形成漏检或者误检;而采用直接将输入信号投影到冗余的即稀疏重表示的字典中，然后转化为优化问题的求解方式虽然在一定程度上提高了检测成功率，但是计算复杂度太高，导致检测效率低下。

发明内容

[0006] 本发明的目的在于针对上述已有技术中的不足，提出一种基于子空间正交向量的脑电信号峰电位检测方法，以降低计算复杂度，提高检测效率。

[0007] 本发明的技术方案是:根据输入脑电信号的特征构建一个不完备的冗余字典，再通过字典构建子空间，其中子空间中的向量是标准正交的，然后通过将输入脑电信号投影到各个子空间中，设定阈值，将保留的系数用反变换的方法重构信号，完成对峰电位的检测。其实现步骤包括如下：

[0008] (1)根据输入信号X e Rlxm的标准峰波，构造不完备冗余字典D e R(-1+1)x'其中m表不输入脑电信号的长度，1表不标准峰波的长度；

[0009] (2)通过步骤（1)中构建的不完备冗余字典D填充子空间V1J2, ···*，··，，其中Vj 表示第j个子空间，1彡j彡1，其中。

[0010] (3)将输入脑电信号分别投影到各个子空间中，得到1个系数向量：C^C2,...， C t, · · .Ci，CteRnX1表示第t个系数向量，1彡t彡1;

[0011] (4)设定阈值T为1.25-1.45,保留Ct中大于阈值T的值并将小于T的值全部置0,将 &，&，...，C t，.. .Ci依次拼接得到重构系数向量LeR(-1+nxl。

[0012] (5)将字典D的伪逆与重构系数向量L相乘得到重构信号Jc，此时；中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即为峰电位。

[0013] 本发明与现有技术相比具有以下优点：

[0014] 1.计算复杂度低

[0015]本发明将脑电信号投影到子空间中，且子空间中的原子是相互正交的，通过设定阈值再反变换得到检测结果，极大的降低了计算复杂度。

[0016] 2.鲁棒性好

[0017] 实验结果表明，本发明对噪声方差为0.05-0.55,信噪比为+30dB--11.6dB的低噪脑电信号峰电位检测准确率可达100%，同时对噪声方差为0.55-0.7,信噪比为-11.6dB一-16.2的高噪脑电信号峰电位的检测准确率到达80%以上，表现出良好的鲁棒性。附图说明

[0018] 图1为本发明的实现总流程图；

[0019] 图2为本发明中冗余字典构建示意图；

[0020] 图3为本发明中子空间构建示意图。

具体实施方式

[0021] 下面结合附图和实例对本发明进行详细说明。

[0022] 参照图1，本发明的实现步骤如下：

[0023]步骤1，根据输入信号构建冗余字典D。

[0024]参考图2,本步骤的具体实现如下：

[0025] (la)设矩阵A(ai，a2. . ·，as，· · .am)TeR(m-1+1)Xm中所有元素为〇，1彡s彡 m，将A中的第3个行向量~的第i个到i+1个元素置为标准峰波，对~进行归一化处理，使其模值I |as I =1 ；

[0026] (Ib)令as = di，其中di表示字典D中的第i个原子，1彡i彡m，得到字典Ddd2...， (11，...')\(1 1的模值为||(11||=1，字典0中的第1行(11是由标准峰波与111-1个0元素构成，标准峰波起始于Cl 1的第i个元素。

[0027] 步骤2，构建子空间V1，V2，…，Vj，."Vi。

[0028]参考图3,本步骤的具体实现如下：

[0029] (2a)由标准峰波与m-1个0元素构成第一子空间V1的原子，其中标准峰波起始于原子Vn的第1个元素、原子V 12的第1+1个元素、原子m的第k+Ι个元素和Vln的第n+1个元素，得到第 1 个子空间Vi(vn，vi2,...，vlk，.. .Vln)TeRnXm，其4

[0030] (2b)由标准峰波与m-1个0元素构成第二子空间V2的原子，其中标准峰波起始于原子V21的第2个元素、原子V22的第2+1个元素、原子V2k的第2+(k-l) Xl个元素和V2n的第2+(n-1) X 1个元素，得到第2个子空间V2(v21，v22, · · ·，V2k，· · .v2n)TeRnXm;

[0031] (2c)由标准峰波与m-1个0元素构成第j个子空间％的原子，1彡j彡I，标准峰波起始于原子Vji的第j个元素、原子Vj2的第j+Ι个元素、原子Vjk的第j+(k-l) Xl个元素和Vjn的第 j+(n-l) X 1 个元素，得到第j个子空间L(Vjl，Vj2,...，vjk，.. .Vjn)TeRnXm;

[0032] (2d)由标准峰波与m-1个0元素构成第1个子空间¥1的原子，标准峰波起始于原子 Vn的第1个元素、原子Vi 2的第2X 1个元素、原子m的第kXl个元素和Vln的第nXl个元素，得到第 1 个子空间Vi(vn，vi2, · · ·，vik，· · .vin)TeRnXm。

[0033] 步骤3,将输入脑电信号投影到子空间。

[0034] 将输入脑电信号分别投影到各个子空间中，得到1个系数向量：C1, C2, ...，Ct， • . .Ci，其中CtGRnxl表示第t个系数向量，1彡t彡1。

[0035] 步骤4,设定阈值，得到重构系数向量。

[0036] 4a)设定阈值T为1.25-1.45,将Ct中的各元素依次与阈值T进行比较，如果C t中元素大于T，则保留该元素，如果Ct元素小于或等于T，则将该元素置为0;

[0037] 4b)将经过3a)步骤处理后的CllC2, ...，Ct，.. .C1按列依次拼接，得到重构系数向量LGR(m-1+1)X1。

[0038] 步骤5，计算重构信号，完成检测。

[0039] 将字典D的伪逆与重构系数向量L相乘得到重构信号χ，此时λ中非零坐标就是检测到的输入信号中的峰波，峰波中的最大值即为峰电位。

[0040] 以上描述仅是本发明的一个具体实例，并不构成对本发明的任何限制。显然对于本领域的专业人员来说，在了解了本发明内容和原理后，都可能在不背离本发明原理、结构的情况下，进行形式和细节上的各种修正和改变，但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。

标题	发布/更新时间	阅读量
一种消除邻近信道干扰的方法、电子设备和可读存储介质	2020-05-08	895
一种快速响应的高精度光散射浊度计测量装置	2020-05-08	575
图像的质量评价方法、电子设备和计算机可读存储介质	2020-05-08	177
基于理想边缘外推的遥感影像超分辨率复原方法	2020-05-08	348
窄带无线通信终端中的同步装置	2020-05-08	260
语音识别方法、装置、存储介质以及电子装置	2020-05-08	251
一种基于图像质量的人脸防伪方法	2020-05-11	174
一种空间光学遥感器有效载荷光机电热一体化分析与优化方法	2020-05-11	290
一种保偏光纤分布式偏振耦合检测弱耦合点提取方法	2020-05-08	543
一种感光基板及显示面板	2020-05-08	662

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

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

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：