我的中期报告1.docVIP

燕 山 大 学 本科毕业设计（论文）中期报告 课题名称： 学院（系）： 电气工程学院 年级专业： 08 学生姓名： 指导教师： 完成日期： 2012-5-1 毕业设计进展情况 针对普拉提运动脑电动力学分析的研究，我已完成开题报告的检查，包括选题依据、研究内容、研究步骤及工作进度的制定。目前处于毕业设计中期工作准备中，即将接受中期进度的检查，现将这阶段的学习情况总结如下： 在中期的准备中，为了实现普拉提运动脑电动力学分析，首先了解了脑电信号方面的基础知识。重点掌握了脑电信号去噪的各种方法。最常用的方法有伪迹减法、回归方法、自适应滤波、小波变换、主成分分析与独立分量分析等。其中独立分量分析具有良好的去噪能力，可一次性去除多种噪声并而不影响有用信号成分，最终选择该方法作为主要的去噪方法。本研究使用Matlab中的EEGLAB工具箱，对采集到的原始脑电信号进行预处理。 经过一个月的学习，基本掌握了Matlab的简单编程。对经过预处理得到的干净脑电信号进行分析之前，做了一个简单的频谱分析。然后使用Matlab小波工具包对信号进行分解以得到四种脑电基本节律即δ、θ、α、β。为实现分解过程编写了一个简单的程序。接着完成了四种节律波的功率谱分析。 今后还有很多的工作要做，工作重点主要有：应用简单的算法计算普拉提运动脑电指标，包括各节律信号能量计算以及对各通道信号进行小波相干性分析。继续熟悉Matlabk软件的使用并编写相关程序。分析采样的脑电信号结果，总结普拉提运动对脑电信号的有何种影响。撰写毕业论文，准备终期答辩。 二、毕业设计具体方案 1、研究方法 本实验共有5位实验者，他们身体状况相似，让他们同时练习普拉提，取他们六次不同练习阶段的脑电图进行对比研究。 2、脑电信号预处理流程图 ① 运行matlab中的EEGLAB的图形界面。 ② 数据导入 EEGLAB中。 ③ 改变采样频率为250HZ。 ④ 过滤数据，设置带通滤波器为0.05-200HZ 。 ⑤ 转换数据为平均参考电极。从Neuroscan采集系统得到的数据改变参考电极为两乳突处。 ⑥ 提取出想要分析的数据段。 ⑦ 删除基线值。 ⑧ 运行ICA ，滤除各种伪迹，得到干净信号。 多通道的一致性分析 ？？ 多通道的相关性分析 ？？ 不同脑区的功率谱变化 ？？ 得到结论性的结果。普拉提对脑功能的影响 3、小波分解与重构 对采样率为250HZ的脑电信号进行小波分解。本研究对脑电数据进行连续小波分解，小波函数选择Daubechchies4,分解层数为5 ，则小波分解后各层的频段范围显示如下： D1 62.5-125HZ D2 31.25-62.5HZ D3 15.625-31.25HZ D4 7.8125-15.3625HZ D5 3.90625-7.8125HZ A5 0-3.90625HZ 其中 A5近似为δ波，D5近似为θ波，D4近似为α波，D3近似为β波。 4、功率谱分析 对提取出来的四种脑电基本节律信号分别进行功率谱估计。功率谱估计有多种方法，一般可分为参数化方法和非参数化方法。本研究采用Welch法中的pwelch函数。 clear all; close all; load s Fs = 250; nfft=1024; t = 0:1/Fs:1; window=boxcar(100); noverlap=20; [Pxx,f]= pwelch(s,window,noverlap,nfft, Fs); PXX=10*log10(Pxx); Plot(f,PXX)； 4、为得到普拉提对脑电活动的影响需要计算各个节律波的总能量，并对比各个不同练习阶段的变化情况。 sum(s.^2) 其中s为各个节律的幅度信息 三、设计过程中遇到的问题及解决办法 1．Matlab的使用 MATLAB有大量的函数，因为之前对它还不太熟悉，所以在应用初期经常会遇到问题。但通过查阅书籍，使用matlab帮助功能函数，逐渐能应用自如了。 2.小波分解层次的确定 小波分解只能将信号分解到特定频率，无法得到所要求的精确的频段。 四、毕业设计时间进度表 完成计划 完成日期 MATLAB软件2012-5-10 对程序进行调试 2012-5-15 加载数据进行分析 2012-5-18 总结分析结果 2012-5-23 撰写毕业论文 201-6-17 准备答辩 2