1 简介
大脑是由亿万个神经元组成的复杂系统,负责身体的各个功能的协调运作,通过大脑皮层上电极记录下
的大脑细胞群的电位活动称为脑电信号。通过对脑电信号进行研究可以获得丰富的心理及生理的疾病信
息,所以脑电信号的分析及去噪算法的研究无论是在临床的诊断和急病治疗上都是十分重要的。一般来
说,脑电信号具有背景噪声强、信号微弱等特点,所以如何消除脑电数据的噪声,更好地获取对大脑有
用的信息是当今研究的热门话题。
近年来小波变换广泛应用于信号处理和图像融合等领域。小波去噪的方法主要有模极大值法、空间滤
波、阈值法,其中阈值法国内外学者广泛应用的方法[2]。阈值的选取是根据信号的长度设定,其设定的
方法具有单一性,无法一次性对信号参数进行选定,由于选取过于零散,无法满足信号复杂等特点。
1 小波变换
小波变换是对基本小波进行伸缩和平移得到的,小波系数是原信号与小波基函数的相似系数。设函数
Φ(t)为可积函数且满足t=L2(R),则傅里叶变换满足如下条件:
2 小波基的选择
小波基的选择对脑电去噪效果至关重要,一般对称性好的小波基不会产生相位畸变[8]。正则性好的小波
基在进行脑电重构时脑电信号更加平滑,紧支撑的小波在脑电去噪时处理的速度更快。
本文选择的小波基为db Nv小波,可以满足以上的条件,db1(简称Haar)小波是小波基当中最简单的
小波形式,局域性较差,因此本文选择db6小波,db6的特点是处理速度更快,MATLAB软件应用方便,
其波形和脑电波形相似,对脑电去噪效果最好。db6小波曲线如图1所示。
图1 db6小波函数
3 小波去噪原理
小波去噪原理是将小波进行分解,再对小波进行多尺度变换,尽可能提取有用的脑电信号。然后再根据
波恩脑电信号特征和噪声特点,选择合适的去噪模型,利用贝叶斯估计后的系数进行小波重构,从而得
到去噪后的脑电信号。其基本原理如图2所示,在进行小波变换前,脑电信号的突变处峰值很高,突变处
集中了噪声,很多平稳有用的脑电信号幅值很小,从而大大增加了去噪难度,所以本文以拉普拉斯为去
噪模型,结合贝叶斯进行小波估计,已达到去噪的目的。