《一种基于奇异谱熵的神经元动作电位特征提取方法.pdf》由会员分享，可在线阅读，更多相关《一种基于奇异谱熵的神经元动作电位特征提取方法.pdf（6页珍藏版）》请在专利查询网上搜索。
1、(10)授权公告号 CN 101828911 B (45)授权公告日 2011.07.20 CN 101828911 B *CN101828911B* (21)申请号 201010159031.3 (22)申请日 2010.04.27 A61B 5/04(2006.01) (73)专利权人 杭州电子科技大学 地址 310018 浙江省杭州市下沙高教园区 2 号大街 (72)发明人 范影乐 钟华 丁颖 (74)专利代理机构 杭州求是专利事务所有限公 司 33200 代理人 杜军 CN 101221554 A,2008.07.16, 全文 . CN 1744073 A,2006.03.08, 全文。
2、 . 张小鹏， 范影乐， 杨勇 . 基于奇异谱熵的脑 电意识任务识别方法的研究 .计算机工程与科 学 .2009, 第 31 卷 ( 第 12 期 ),117-120. (54) 发明名称 一种基于奇异谱熵的神经元动作电位特征提 取方法 (57) 摘要 本发明涉及一种基于奇异谱熵的神经元动作 电位特征提取方法。 现有的方法信息提取不完成、 计算效率低。 本发明方法首先采样动作电位信号， 并在设定的时间窗宽动作电位信号形成时间序列 X ； 其次设定嵌入维数 K 和延迟时间 的取值， 对 该时间序列 X 进行多维相空间重构得到相空间矩 阵 Y ； 然后对相空间矩阵 Y 进行奇异值分解， 得到 相空。
3、间矩阵 Y 的奇异值， 分别计算各非零奇异值 与所有奇异值总和之比得到概率， 并利用信息熵 计算方法计算时间序列 X 的奇异谱熵值 ； 最后移 动时间窗， 依上述方法计算所有奇异谱熵值， 即得 到神经元动作电位特征。本发明方法在计算过程 中的信息量完整且计算速度快。 (51)Int.Cl. (56)对比文件 审查员 张宇 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书 1 页 说明书 3 页 附图 1 页 CN 101828911 B1/1 页 2 1. 一种基于奇异谱熵的神经元动作电位特征提取方法， 其特征在于该方法包括如下步 骤 ： 步骤 (1) 通过动作电位采集系统，。
4、 以 40KHz 采样频率采集得到多个时长为 T 个采样点 的动作电位信号 ； 步骤 (2) 设定时间窗宽 n 的取值， 使 n T， 且 n 个采样点的时长应包括上升支下降支 神经元放电的重要时段， 并记录任一动作电位信号在该时间窗下的各采样点幅值， 形成时 间序列 X ； 步骤 (3) 设定嵌入维数 K 和延迟时间 的取值， 对该时间序列 X 进行多维相空间重构 得到相空间矩阵 Y ； 所述的多维相空间重构的具体方法为 ： 首先在时间序列 X 中从第一个动作电位采样 点为起始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第一个子序列作为 相空间的第一个行向量 ； 其次从第 +。
5、1 个动作电位采样点为起始点， 提取间隔为 的连 续 K 个采样点的动作电位幅值， 组成第二个子序列作为相空间的第二个行向量 ； 再次从第 2+1 个动作电位采样点为起始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第三个子序列作为相空间的第三个行向量， 然后从第 3+1 个动作电位采样点为起 始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第四个子序列作为相空间 的第四个行向量， 以此类推， 从时间序列 X 中共可提取个不同的子序列， 记 则得到一个 NK 阶相空间矩阵 Y ； 步骤 (4) 对相空间矩阵 Y 进行奇异值分解， 得到相空间矩阵 Y 的奇异值， 即为。
6、非负对角 阵 的对角线上的所有特征值， 其中 UTYV ， U 为 NN 阶的正交矩阵， V 为 KK 阶的 正交矩阵 ； 步骤 (5) 分别计算各非零奇异值与所有奇异值总和之比得到概率， 并利用信息熵计算 方法计算时间序列 X 的奇异谱熵值 ； 步骤 (6) 移动时间窗， 记录得到下一个时间序列， 重复步骤 (3) 至步骤 (5) 计算该时间 序列的奇异谱熵值 ； 步骤 (7) 重复步骤 (6)， 计算出该动作电位信号可提取到的所有奇异谱熵值， 即为神经 元动作电位特征。 权 利 要 求 书 CN 101828911 B1/3 页 3 一种基于奇异谱熵的神经元动作电位特征提取方法 技术领域 。
7、0001 本发明属于生物医学工程领域， 涉及一种神经元动作电位特征提取方法， 具体涉 及一种基于奇异谱熵的神经元动作电位特征提取方法。 背景技术 0002 神经元动作电位的特征提取技术是动作电位模式分类技术以及动作电位序列解 码等神经信息学研究的前期基础， 且由于植入式多电极阵列技术的发展， 使得大量神经元 动作电位的获取成为可能， 因此提取动作电位的有效特征， 并根据特征将其划分至各自对 应的神经元， 对于后续理解神经元响应与外在刺激之间的关联性， 起着非常关键的作用。 0003 目前神经元动作电位特征提取及模式分类方法主要包括模板法和信号处理方法。 前者需要一定的先验知识， 要求能够获取信。
8、号中所隐含的动作电位类别数和各自的动作电 位模板特征， 因此分类性能的稳定性通常不够理想 ； 后者则将采样后的动作电位信号视为 多维的点处理时间序列， 通过对其特征的刻画， 来实现信号的自动分类。 目前被普遍采用的 特征提取方法有主成分分析和小波分析方法。 主成分分析通过样本协方差矩阵求取主要特 征值， 但由于动作电位信号具有一定的非线性时变性， 因此主成分分析可能无法反映动作 电位信号的完整信息。小波分析方法通过对动作电位信号进行多层分解， 在时频域上对动 作电位的动态特性进行有效描述， 但是在小波分析中所使用的小波函数具有多样性， 因此 分析得到的小波分量和小波谱只相对所选择的小波基有意义。
9、。 发明内容 0004 本发明针对现有技术的不足， 提供了一种基于奇异谱熵的神经元动作电位特征提 取方法。 该方法提取的动作电位奇异谱熵特征， 能有效地用于后期动作电位的分类， 可解决 传统方法所提取的特征无法反映动作电位完整信息或是无法反映动作电位非线性非平稳 的动态特性等问题。 0005 本发明方法包括以下步骤 ： 0006 步骤 (1) 通过动作电位采集系统， 以 40KHz 采样频率采集得到多个时长为 T 个采 样点的动作电位信号 ； 0007 步骤 (2) 设定时间窗宽 n 的取值， 并记录任一动作电位信号在该时间窗下的各采 样点幅值， 形成时间序列 X ； 0008 步骤 (3) 。
10、设定嵌入维数 K 和延迟时间 的取值， 对该时间序列 X 进行多维相空间 重构得到相空间矩阵 Y ； 0009 所述的多维相空间重构的具体方法为 ： 首先在时间序列 X 中从第一个动作电位 采样点为起始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第一个子序列 作为相空间的第一个行向量 ； 其次从第 +1 个动作电位采样点为起始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第二个子序列作为相空间的第二个行向量 ； 再次 从第 2+1 个动作电位采样点为起始点， 提取间隔为 的连续 K 个采样点的动作电位幅 说 明 书 CN 101828911 B2/3 页 4 值，。
11、 组成第三个子序列作为相空间的第三个行向量， 然后从第 3+1 个动作电位采样点为 起始点， 提取间隔为 的连续 K 个采样点的动作电位幅值， 组成第四个子序列作为相空间 的第四个行向量， 以此类推， 从时间序列 X 中共可提取个不同的子序列， 记 则得到一个 NK 阶相空间矩阵 Y ； 0010 步骤 (4) 对相空间矩阵 Y 进行奇异值分解， 得到相空间矩阵 Y 的奇异值， 即为非负 对角阵 的对角线上的所有特征值， 其中 UTYV ， U 为 NN 阶的正交矩阵， V 为 KK 阶 的正交矩阵 ； 0011 步骤 (5) 分别计算各非零奇异值与所有奇异值总和之比得到概率， 并利用信息熵 。
12、计算方法计算时间序列 X 的奇异谱熵值 ； 0012 步骤 (6) 移动时间窗， 记录得到下一个时间序列， 重复步骤 (3) 至步骤 (5) 计算该 时间序列的奇异谱熵值 ； 0013 步骤 (7) 重复步骤 (6)， 计算出该动作电位信号可提取到的所有奇异谱熵值， 即为 神经元动作电位特征。 0014 本发明方法所具有的有益效果为 ： 0015 1、 奇异谱熵是一种衡量一维时间序列复杂度的平均熵参数， 适用于随机信号或者 确定信号以及随机和确定性信号混合的信号。将奇异谱熵应用于动作电位的特征提取中， 其特征表征了该动作电位在奇异值分解下的信息量， 可有效实现动作电位的分类。 0016 2、 。
13、基于奇异谱熵的神经元动作电位特征提取方法， 对动作电位信号进行多维相空 间重构， 并在重构的相空间内进行奇异值分解， 其分解后的奇异值大小反映了该分量在动 作电位信号中的重要程度 ； 且多维相空间重构方法， 也一定程度上解决了传统方法在信号 分解重构过程中丢失信息的问题。 0017 3、 考虑到动作电位具有较强的非平稳性和非线性， 因此基于奇异谱熵的动作电位 特征提取方法， 在奇异值分解基础上采用熵值评价， 即奇异谱熵来描述在奇异值分解下动 作电位的特征， 反映了动作电位的整体信息以及在奇异谱划分下的不确定性， 是对动作电 位的一种非线性度量。 0018 4、 奇异谱熵算法， 只需较短的序列长。
14、度就能估计出较为稳定的统计值， 能较好地 适用于动作电位时长较短的特点。 0019 5、 基于奇异谱熵的神经元动作电位特征提取方法， 在相空间重构过程中， 时间窗 宽、 嵌入维数和延迟时间的取值将影响相空间的质量以至于影响动作电位的奇异谱熵特征 的质量。 附图说明 0020 图 1 是神经元动作电位波形图 ； 0021 图 2 是神经元动作电位奇异谱熵图。 具体实施方式 0022 以下将结合附图对本发明作进一步说明， 本发明方法利用一维的动作电位信号重 构多维相空间， 在该相空间内进行奇异值分解， 得到动作电位的奇异谱， 并采用熵值描述奇 异值分解下的动作电位特征。 奇异谱熵特征反映了动作电位。
15、的整体信息以及在奇异谱划分 说 明 书 CN 101828911 B3/3 页 5 下的不确定性， 可以作为区分不同动作电位的有效依据。本发明方法的具体步骤为 ： 0023 步骤 (1) 通过动作电位采集系统可采集到时长为 T 个采样点的动作电位信号。为 使每个动作电位都能提取一定数量的奇异谱熵特征， 以便于后期分类时特征的选择， 且所 提取的特征能反映该动作电位的特性， 所以设定合适的时间窗宽 n 使 n T， 且 n 个采样点 的时长应包括上升支下降支等神经元放电的重要时段。如图 1 所示动作电位 T 64， 通过 观察动作电位波形可取 n 12。 0024 步骤 (2) 记录任一动作电位。
16、信号在该时间窗下的各采样点幅值， 形成第一个时间 序列 X x1， x2， x3， .， xn， 选择嵌入维数 K 和延迟时间 重构其相空间， 则相空间行向 量为 Yi xi， xi+， xi+2， .， xi+(K-1)， 其中嵌入维数和延迟时 间的取值也会影响最终提取的动作电位奇异谱熵特征质量。由于动作电位时长较短， 所以 两者取值一般都较小， 但太小会影响重构相空间的质量。如图 1 所示， 本实施例中 K 6， 1。 0025 步骤 (3) 令则重构的多维相空间可表示 为其中 T 表示矩阵转置。Y 为一个 76 阶矩阵。对 Y 进行奇异值 分解， 使得 UTYV ， 其中 U 为一个 7。
17、7 阶正交矩阵， V 为一个 66 阶正交矩阵， 为 76 阶的非负对角阵，S diag(1， 2， 3， .， r)。非零的特征值 1， 2， 3， .， r以及为零的 r+1， r+2， .， 6为多维相空间 Y 的全部奇异值， 其中 0 r 6。 0026 步骤(4)定义非零奇异值i与所有非零奇异值总和之比为其贡献率pi， 如式(1) 所示， 计算所有非零奇异值的贡献率。 0027 0028 0029 步骤(5)通过式(2)求得时间序列X的奇异谱熵值SSE， 移动时间窗截取不同时间 序列， 重复步骤(2)至步骤(4)可计算出该动作电位的所有时间序列的奇异谱熵值， 作为该 动作电位可提取的特征， 进而可求得所有待分类动作电位的不同奇异谱熵值， 作为后期动 作电位分类的特征依据。对于图 1 所示的动作电位， 其奇异谱熵曲线如图 2 所示。 说 明 书 CN 101828911 B1/1 页 6 图 1 图 2 说 明 书 附 图 。