本发明专利技术涉及一种基于双密度小波邻域相关阈值处理的脑电信号消噪方法。当前大都采用经典离散小波变换结合传统阈值法对脑电信号进行消噪，而已有的基于经典小波变换结合传统阈值法的消噪方法都存在不足之处。本发明专利技术首先从大脑皮层上采集脑电信号，然后用双密度小波正变换对脑电信号进行分解，得到多层的信号高频系数，根据小波系数的局部统计依赖性，运用邻域相关阈值处理算法进行收缩，最后将收缩后的小波系数进行重构得到消噪后的信号。本发明专利技术根据脑电信号特性和干扰噪声特点，以信噪比为目标函数，采用网格寻优法对邻域相关阈值处理算法中的三个可调参数进行寻优，进而有效地平滑噪声并保留EEG的细节特征。
【技术实现步骤摘要】
本专利技术属于生物医学信号消噪领域，涉及一种。
技术介绍
脑电信号(Electroencephalogram, EEG)是中枢神经系统产生的生物电活动而衍射到脑皮的募集信号，人在主动思维或受到不同的感觉刺激时，EEG的特征具有明显差异。对EEG进行分析，可以获得大量的生理、心理及病理信息。EEG是一种非线性非平稳信号，而且非常微弱，其幅值仅为UV级，极易受到噪声干扰，如脉冲干扰、工频干扰、呼吸干扰、眼动干扰、心电干扰、肌电干扰、头皮电极的抖动等，这些干扰产生的噪声混杂在EEG中，使得EEG不能正确反映大脑的真实生物电活动特征,给医生诊断带来误诊，也同样给研究人员对EEG的进一步深入研究，例如基于EEG的脑机接口(BCI)中运动意识任务的模式识别带来困难。因此，消除EEG中的噪声就显得非常必要。小波变换因其具有良好的时频局部化特性、多分辨率特性、去相关性等特点，在生物医学信号消噪领域得到了广泛应用，也是传统的脑电信号消噪方法。经典离散小波变换虽功能强大，但它采用的是临界采样，导致经典离散小波变换具有数据敏感性，输入数据在时域上微小的变化会对离散小波变换系数产生不可预测的结果，即使输入有很小的移位都可能导致在各个尺度上的能量分布有明显的变化，这将降低BCI系统的EEG控制信号鲁棒性。小波系数阈值处理策略的选取是影响基于小波变换阈值消噪方法降噪效果的主要因素之一。信号经小波分解后，每个子空间内的小波系数具有一定的统计依赖性，利用这种局部依赖性，可以有效地提高信号的降噪效果。但传统软、硬阈值算法只是单独对每一个小波系数进行阈值处理，并未考虑小波系数的局部依赖性，造成系数估计值偏差过大，进而使得传统小波阈值消噪方法对EEG的消噪性能降低。综上所述，已有的基于经典小波变换和传统软阈值法的脑电信号消噪方法存在不足之处。
技术实现思路
本专利技术的目的就是针对现有技术的不足，提出一种基于双密度小波邻域相关阈值处理的脑电信号消噪新方法。双密度小波变换(Double-Density Discrete WaveletTransform, DD-DffT)是一种新的小波分析方法,它的滤波器组是由一个尺度函数和两个彼此之间偏移0.5个单位的小波函数构成，同时使用过采样代替临界采样，因而它很好的客服了经典离散小波变换的局限性，具有运算速度快、近似的平移不变性、完美的重构性和有限的冗余性等特点，且与经典小波变换相比能够更准确地描述信号的细节特征，进而提高信号特征的鲁棒性；邻域相关阈值处理算法在小波系数的估计计算中融入了邻域内小波系数的信息，可以减小系数估计值的偏差，进而有效地平滑噪声并保留EEG的细节特征。实现本专利技术方法的主要思路是:利用双密度小波对采集到的EEG分解，得到多层的信号高频系数；根据小波系数的局部统计依赖性，运用邻域相关阈值处理算法进行收缩，将收缩后的小波系数进行重构得到消噪后的信号。为了实现以上目的，本专利技术方法主要包括以下步骤:步骤一.获取人大脑运动想象脑电信号样本数据。步骤二.将步骤一获取的运动想象脑电信号进行双密度小波分解，得到低频小波系数和各层高频小波系数，分别记为G(I)其中J表示最大分解尺度；j表示分解尺度；h1、h2表示两个高通滤波器；1表示小波系数序列的标号，其中(I)、Wih2(I)分别表示子空间(Ih1)和子空间(j，h2)内的小波系数序列。步骤三.计算各分解子空间的Donoho阈值，记为入j, hk:本文档来自技高网...
【技术保护点】
基于双密度小波邻域相关阈值处理的脑电信号消噪方法，其特征在于该方法包括如下步骤：步骤一.获取大脑运动想象脑电信号样本数据；步骤二.将步骤一获取的运动想象脑电信号进行双密度小波分解，得到低频小波系数和各层高频小波系数，分别记为CJ(l)和wj,h1(l)、wj,h2(l)；其中J表示最大分解尺度；j表示分解尺度；h1、h2表示两个高通滤波器；l表示小波系数序列的标号，其中wj,h1(l)、wj,h2(l)分别表示子空间(j,h1)和子空间(j,h2)内的小波系数序列；步骤三.计算各分解子空间的Donoho阈值，记为λj,hk：λj,hk=median(|wj,hk|)*2InM/0.6745其中，M为子空间(j,hk)内的小波系数数目；步骤四.定义加权阈值缩放因子λj，计算式为：λj＝β/(1+In(j))其中，j为分解尺度；步骤五.对Donoho阈值进行缩放得到各分解子空间的阈值，记为Tj,hk：Tj,hk＝λj*λj,hk根据分解尺度对频率相对高子空间的阈值进行放大、频率相对低子空间的阈值进行缩小，进而增强对EEG中高频噪声的抑制，同时保留低频有用信息；步骤六.计算以wj,hk(l),k＝1,2为中心，大小为2m+l的操作邻域窗口内小波系数的均值,记为：w‾j,hk(l)=12m+1Σm1=-mmwj,hk(l+m1)步骤七.使用邻域相关阈值处理函数对高频系数wj,hk(l)进行收缩，其函数构造为:以信噪比为目标函数，使用网格搜索法对m、β、K三个可调参数进行了 寻优，当m＝2、β＝1.15、K＝4时，算法的消噪效果最为理想；参数寻优的具体步骤如下：(1)根据真实EEG波形特点，生成采样率Fs、时间t=0:1/Fs:4、频率范围为2～30Hz的标准信号s(t)=Σm=2132-rand(1)sin(2πnt+2*rand(1)*π)+Σn=1430rand(1)5cos(2πnt+2*rand(1)*π),与EEG的主要节律对应；(2)对标准信号加噪声x(t)＝awgn(s(t),n,“measured“)，即标准信号加噪后的信噪比为n（dB）；(3)将x(t)进行双密度小波分解，得到各层高频小波系数，记为wj,h1(l)、wj,h2(l)；计算各分解子空间的Donoho阈值，记为λj,hk；(4)设定m∈[0，(M?1)/2]，M为最大分解子空间内小波系数数目，步长Δm=1；设定β∈[0，20]，步长Δβ=0.05；设定K∈[1，60]，步长ΔK=0.05；然后，根据步骤四到步骤七计算每种参数组合(m，β，K)的消噪信号以及信噪比SNR（m,β,K）；信噪比定义为：式中，N是信号长度，s(i)是原始标准信号，为去噪后的信号；通常信噪比越大信号的整体消噪效果越好，因此，只要找出SNR（m,β,K）的最大值就可以寻得最优的(m，β，K)参数组合；步骤八.将低频系数和经收缩后的高频系数进行双密度小波逆变换，得到消噪后的EEG。FDA0000397655720000015.jpg,FDA0000397655720000014.jpg,FDA0000397655720000013.jpg,FDA0000397655720000023.jpg,FDA0000397655720000022.jpg,FDA0000397655720000024.jpg...
【技术特征摘要】
1.基于双密度小波邻域相关阈值处理的脑电信号消噪方法，其特征在于该方法包括如下步骤: 步骤一.获取大脑运动想象脑电信号样本数据； 步骤二.将步骤一获取的运动想象脑电信号进行双密度小波分解，得到低频小波系数和各层高频小波系数，分别记为C1(I)和^hl ...
【专利技术属性】
技术研发人员：罗志增，周瑛，席旭刚，高云园，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。