[发明专利]一种自动识别和提取睡眠脑电中K复合波的方法

摘要

一种自动识别和提取睡眠脑电中K复合波的方法，包括以下步骤：对脑电信号进行小波分解和重构；对重构的数据采用Teager energy能量算子的运算并取绝对值；对得到的Teager能量曲线平滑并0/1粗粒化；对粗粒化的数据进行阈值检测；对满足阈值的重构信号进行形态学检测，满足形态学检测条件则该位置的信号为K复合波并记录下起点和终点位置以及波峰波谷的值和位置；本发明对小波分解重构后的信号用Teager能量算子进行分析，提取Teager能量算子的绝对值序列并进行平滑、粗粒化处理，该方法实现容易，抗噪能力强，并能比较准确地确定K复合波的位置以及波峰波谷的值和位置，为睡眠分期中非眼动2期的识别及K复合波的研究奠定了基础。

主权项

一种自动识别和提取睡眠脑电中K复合波的方法，其特征在于，具体实现步骤如下：步骤一、读取睡眠脑电信号x(t)，对x(t)分段，段长为30s，对每一段的数据利用db5小波进行八层的小波分解，依据K复合波的主要分布范围，选择第7、8层的低频系数进行重构，得到预处理后脑电信号y(t)；步骤二、对预处理后脑电信号y(t)的最大值max(y(t))和最小值min(y(t))进行对比,若最大值和最小值的差的绝对值大于100μv时，则对预处理后脑电信号y(t)进行Teager能量算子(TEO)的计算并取其绝对值，得到Teager能量算子(TEO)的绝对值序列z(n)，若不满足上述条件，则不处理，进行下一段信号的处理；步骤三、对Teager能量算子(TEO)的绝对值序列z(n)进行平滑处理，方法为加矩形窗求平均值,窗长为125，步长为1，得到新的序列o(n)；步骤四、将平滑后的序列o(n)进行粗粒化转化为0/1序列p(n)；步骤五、检测0/1序列p(n)中连续1的持续时间，若持续时间大于0.5s，则记录其起始和终止的位置，以此为据，去查找重构信号相对应的位置，然后对该位置的重构信号进行形态学检测，满足形态学检测条件则该位置的信号为K复合波，并记录下起点和终点位置以及波峰波谷的值和其位置。