matlab eeg信号处理,基于MATLAB的脑电信号处理.pdf

南京航空航天大学

基于Matlab 的脑电信号处理

姓名 陆想想

专业领域 生物医学工程

课程名称 数字信号处理

二О一三年四月

摘 要：脑电信号属于非平稳随机信号，且易受到各种噪声干扰。本文基于Matlab仿真系统，

主要研究了小波变换在脑电信号处理方面的应用，包括小波变换自动阈值去噪处理、强制去

噪处理，以α波为例，提取小波分解得到的各层频率段的信号，并做了一定的分析和评价。

关键词：脑电信号；小波变换；去噪重构；频谱分析

0 引言

脑电信号EEG(Electroencephalograph)是人体一种基本生理信号，蕴涵着丰富的生理、

心理及病理信息，脑电信号的分析及处理无论是在临床上对一些脑疾病的诊断和治疗，还是

在脑认知科学研究领域都是十分重要的。由于脑电信号的非平稳性且极易受到各种噪声干扰，

特别是工频干扰。因此消除原始脑电数据中的噪声，更好地获取反映大脑活动和状态的有用

信息是进行脑电分析的一个重要前提。本文的研究目的是利用脑电采集仪器获得的脑电信号，

利用Fourier 变换、小波变换等方法对脑电信号进行分析处理，以提取脑电信号α波的“梭形”

节律，并对脑电信号进行功率谱分析和去噪重构。

1 实验原理和方法

1.1实验原理

1.1.1 脑电信号

[1]

根据频率和振幅的不同，可以将脑电波分为4 种基本类型 ，即δ波、θ波、α波、β波。

4 种波形的起源和功能也不相同，如图1所示。

图1 脑电图的四种基本波形

α波的频率为8~13Hz，振幅为为20~100µV，它是节律性脑电波中最明显的波，整个皮

层均可产生α波。正常成人在清醒、安静、闭目时，波幅呈现有小变大，再由大变小，如此

反复进行，形成所谓α节律的“梭形”。每一“梭形”持续时间约为1~2s。当被试者睁眼、警觉、

思考问题或接受其他刺激时，α波立即消失而代之以快波，这种现象称之为“α波阻断”。一般

认为，α波是大脑皮质处于清醒安静状态时电活动的主要表现。

β波的频率是18~30Hz，振幅为5~20µV，是一种快波。β波的出现以为着大脑比较兴

奋。

θ波的频率是4~7Hz，振幅为10~50µV，它是在困倦时，神经系统处于抑制状态时所记

录的波形。

δ波在睡眠、深度麻醉、缺氧或大脑有器质性病变时出现，频率是 1~3.5Hz，振幅为

20~200µV。

1.1.2 小波变换

小波变换的概念是由从事石油信号处理的法国工程师J.Morlet 在 1974 年首先提出的。

与Fourier 变换相比，小波变换是时间 (空间)频率的局部化分析，它通过伸缩平移运算，

对信号逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频

信号分析的要求，从而可以聚焦到信号的任意细节，解决了Fourier 变换的困难问题。

在噪声中如何准确地检测到信号一直是信号处理领域所关心的内容。小波变换是一种信

号的时间一尺度分析方法，由于具有多分辨率分析的特点，良好的时频局部化特性，能够对

各种时变信号进行有效的分解，从而较好地将信号与噪声加以分离，获得满意的去噪效果。

小波分析去噪处理的方法一般有三种：

默认阈值去噪处理。该方法是利用ddencmp 函数生成信号的默认阈值，然后利用

wdencmp 函数进行去噪处理。

给定阈值进行去噪处理。在实际的去噪过程中，阈值往往可以通过经验公式获得，并且

这种阈值比默认阈值可信度高。在进行量化处理时可利用wthresh 函数。

强制去噪处理。该方法是将小波分解结构中的高频或者低频系数设置为0，即滤掉所有

高频部分或低频部分。这种方法比较简单，且去噪得到的信号比较平滑，但是容易丢失信号

[2]

中的有用成分 。

本文采用了两种去噪方法，并分析比较了他们的去噪效果。

1.2实验方法与步骤

1.2.1 脑电信号的读取

本文使用的脑电数据是使用南航生物医学光子学实验室的脑电采集系统采集获得的，原

始数据格式为.eeg。为了方便在Matlab 环境下对数据进行分析，将文件转换为.tx