[发明专利]一种脑电信号多尺度定量分析方法

说明书

技术领域

本发明涉及一种脑电信号分析方法。特别是涉及一种脑电信号多尺度定量分析方法。

背景技术

脑电图是大脑特定部位神经细胞群在生理过程中自发或诱发的电活动在头皮表面或颅内特定部位引起的电位信号，它反映了大脑生物电的节律性活动规律。脑电信号是一种高度复杂的非线性信号，目前临床工作中仍依靠医生对视觉脑电图的经验性判断，所以对于脑电图的分析基本停留于视觉判断、经验性的感性认识，几乎处于定性阶段，带有明显的主观经验依赖性，缺乏客观的定量指标。应用新的物理概念、新的数字信号分析方法对脑电信号进行定量分析，提取有助于临床诊断和发病机理研究的新的参数的方法具有重要的意义。

子波分析是一种自适应的信号分析方法，被称为“数字显微镜”。子波分析应用于脑电信号的分析具有独特的优势，可以精确观察不同频率(尺度)脑电信号的细节特征及随时间的演变规律。通过比较不同生理病理状态下脑电信号的多尺度定量特征参数，提取有助于疾病诊断和分型的电生理参数。

子波分析(wavelet analysis)是近几年新发展起来的一种数学方法，通过信号与一个被称为子波的解析函数进行卷积，将含有多尺度成分的复杂信号进行分解。其中子波基函数是由一个子波母函数经过平移和伸缩变换而得到的一族外形相似但尺度不同的局部波动函数。子波函数在物理空间(时间、空间)代表局部确定位置和有限尺度范围内的具有特定物理过程或形状的“事件”或者“扰动小波”。子波变换就是将含有多尺度成分的复杂信号与一定位置一定尺度下代表特定的物理过程或形状的子波基函数进行局部互相关分析，如果其相关性高，表明该信号在该处含有该子波成分，使得信号在该处于该子波具有较好的相似性。利用子波变换可以将含有多尺度成分的复杂信号中的不同成分分解为多种单一尺度的成分，每一种单一尺度的成分在物理空间(时间、空间)局部确定位置和有限尺度范围内具有特定物理过程或形状。

发明内容

本发明所要解决的技术问题是，提供一种能够使EEG技术在临床诊断和基础研究中发挥更大作用的脑电信号多尺度定量分析方法。

本发明所采用的技术方案是：一种脑电信号多尺度定量分析方法，包括如下步骤：

1)采用多导数字脑电图仪和16位A/D采集卡，采集不同导联的数字脑电信号，将采集到的原始脑电数据存入计算机，并根据分析内容设定采样频率和分析数据的时间长度；

2)对原始脑电数据分析进行连续子波变换

设一维脑电信号s(t)在子波函数下的子波分析定义为：

其中，W_s(a,b)为子波系数，W_ab(t)为子波函数族，是由子波的母波函数W(t)经过平移b和伸缩变换a而来：

其中a表示时间尺度参数，b表示时刻参数；

3)根据多尺度子波系数,信号的总能量按照尺度进行分解，各尺度信号占有的能量的总和等于信号的总能量。

步骤3)是根据子波系数W_s(a,b)，一维脑电信号s(t)的能量进行分解：

其中各尺度子波功率谱密度随尺度参数的变化：

其中，C_W为子波母函数的允许常数，与每种子波母函数的傅里叶谱有关。

按照下式计算不同尺度脑电信号功率的百分比，观察不同频率脑电信号多尺度功率谱随尺度即频率的分布图：

本发明的一种脑电信号多尺度定量分析方法，将子波分析应用于定量EEG的分析，具有如下效果：

1，它不需要信号是平稳的，可以自动分解和提取信号瞬时变化的局部多尺度特征；

2，可以同时提取信号时域、频域的信息及瞬时演变过程。

采用本发明的方法可以将看似复杂无序的脑电信号分解在不同的尺度上，实现脑电信号在每一尺度上的精确表达。从而使EEG技术在临床诊断和基础研究中发挥更大的作用。

附图说明

图1是正常成人视觉脑电图；

图2a是正常成人脑电信号子波系数等值线图；

图2b是正常儿童脑电信号子波系数等值线图；

图2c是儿童失神癫痫发作期脑电信号子波系数等值线图；

图2d是青少年肌阵挛癫痫发作期脑电信号子波系数等值线图；

图3a是正常成人脑电信号分尺度功率随频率的分布图；

图3b是正常儿童脑电信号分尺度功率随频率的分布图；

图3c是儿童失神癫痫脑电信号分尺度功率随频率的分布图；