图3信号频谱图

从频谱图中可以看出,低频信号和11Hz左右的信号特别强,25Hz以上的信号几乎为零。由于α波的频率为8~13Hz,由此可以知道,该信号中α波比较多。在脑电采集过程中,被采集者没有处于深度睡眠状态,因此接近0Hz的低频信号可以确定为噪声。15~25Hz频段的信号很微弱,因此可以以判断出,信号中几乎没有β波。

谱估计是数字信号处理的主要内容之一,主要研究信号在频域中的各种特征,目的是根据有限数据在频域内提取被淹没在噪声中的有用信号。调用eeg_psd,m文件,可以绘制信号的功率谱图,如下图4所示。

图4 信号功率谱图

原始的功率谱图是关于f=128Hz对称的,为了便于分析,截取了0~65Hz的一段。信号的功率谱显示,在10Hz作用频率处能量出现峰值,表明在平静状态下采集的脑电信号中α波能量最大,符合生理学的研究结论。

2.3 脑电信号小波分解各层重构波形

实验选取db8小波对前述采集的样本信号进行分解,调用wavelet_dec.m文件,分解得到的各层信号如下图5所示。

图5 脑电信号经小波分解后的各层分量

由于采样频率是256Hz,所以对于自带的频率范围是:

[0,fsfsfsfsfsfsfsfsfsfsfsfsfsfsfs],[,],[,],[,],[,],[,],[,],[,]28282727262625252424232322222

频率范围 0Hz~1Hz 1Hz~2Hz 2Hz~4H 4Hz~8Hz 8Hz~16Hz 16Hz~32Hz 32Hz~64Hz 64Hz~128Hz 128Hz~256Hz

主要信号成分 低频干扰 δ波 δ波 θ波 α波 β波 高频噪声 高频噪声 高频噪声

相应的各子带的频率成分如下表所示:

子带信号 a8 d8 d7 d6 d5 d4 d3 d2 d1

表1 小波分解后脑电信号子带的频率范围

一般情况下,a8子带内是低频干扰,d3、d2、d1子带内是高频噪声,d8、d7、d6、d5、d4子带内是脑电信号,但也可能混有一定的噪声,需要根据实际情况来分析。

2.4 脑电信号节律提取

α波是节律性脑电波中最明显的波。由以上小波分解得到的各层分量分析可知,α波主要集中在d5子带内。所以选取了d5子带内第5500~8000之间的2500个数据点,将信号图绘制出来,如下图6所示。

图6 提取出来的α波节律

从上图中可以明显看出α波波幅由小到大,再由大到小作规律性变化,呈棱状图形。

2.5 小波分解去噪和重构波形

2.5.1

小波变换默认阈值去噪处理

matlab脑电信号处理,基于matlab的脑电信号处理相关推荐

  1. Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分、篮板、助攻

    Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分.篮板.助攻 目录 输出结果 预测勒布朗詹姆斯的20 ...

  2. matlab功率谱的峰值,基于Matlab功率谱密度估计方法.doc

    基于Matlab功率谱密度估计方法 基于Matlab功率谱密度估计方法要在实际情况下, 许多平稳信号无法导出数学表达式, 要准确获取这些信号的功率谱密度存在一定的困难.根据维纳辛钦 (Wiener K ...

  3. 如何用matlab进行脑电信号处理,基于matlab的脑电信号处理

    调用noise_reduce.m文件,可以实现脑电信号的小波变换默认阈值去噪处理,原始数据及去噪处理结果对比如下图7所示. 图7 原始信号与小波默认阈值去噪结果图的对比 从原始脑电信号与去噪处理后的效 ...

  4. matlab中对伺服电机,基于MATLAB仿真的伺服电机动态特性分析

    基于MATLAB仿真的伺服电机动态特性分析 摘要:伺服电机(servo motor )是指在伺服系统中控制机械元件运转的发动机,是一种补助马达间接变速装置.伺服电机可使控制速度,位置精度非常准确,可以 ...

  5. 基于matlab的msk仿真,基于MATLAB的MSK系统的仿真

    基于MATLAB的MSK系统的仿真Tag内容描述: 1.基于MATLAB的MSK系统的仿真设计1绪论11数字通信的发展通信按照传统的理解就是信息的传输与交换,为了传递消息,各种消息需要转换成电信号,消 ...

  6. 基于matlab fdma传输系统设计,基于MATLAB的LTE系统仿真研究

    摘  要: 根据LTE系统的原理和模块构成,分析并建立了一个基于MATLAB的系统级仿真平台.针对LTE系统中被广泛认可的轮询.比例公平和最大载干比三种经典调度算法,利用仿真平台产生的数据对三者的性能 ...

  7. matlab自适应逆控制,基于matlab仿真模块的自适应有源噪声逆控制研究

    ATE 农业技术与装备 X2 c(q-1) p(q-1) u2 C( x) p(q-1) LMS M(q-1) p(q-1) w( q-1) w( q-1) ! ! e2 d2 + - - + p(q ...

  8. 基于matlab的信号去噪,基于matlab的信号去噪研究.doc

    基于matlab的信号去噪研究.doc 学 生 毕 业 设 计(论 文) 课题名称 基于 MATLAB 的信号去噪研究 姓 名 陈文 学 号 051220206 院 系 物电系 专 业 电子信息工程 ...

  9. matlab生产计划问题,基于MATLAB的生产计划最优化系统设计

    龙源期刊网 http://www.doczj.com/doc/d0588b0949d7c1c708a1284ac850ad02de8007fa.html 基于MATLAB的生产计划最优化系统设计作者: ...

  10. matlab仿真疏散,276基于matlab的疏散仿真程序简介

    基于matlab的疏散仿真程序简介 朱伟向大海刘方 重庆大学城市建设与环境工程学院 400045 摘要:本文简要介绍了人员疏散仿真现状,作者基于matlab平台开发了疏散仿真模型,该仿真模型基于精细网 ...

最新文章

  1. Android TextVeiw 在java代码中添加空格的方法
  2. $.ajax上传文件或者上传图片
  3. 华中科技大学应用高等工程数学_专业解析【第152期】| 机械电子工程课程设置及研究方向...
  4. FLASK上传时有中文文件名的解决方案
  5. php开发神器,PHP开发神器vim
  6. 蓝桥杯 反置数 数论
  7. Kerberos KDC not reachable
  8. 火狐控制台的html,怎么使用火狐浏览器调试网页
  9. 用计算机算出女儿身高,孩子身高预测计算器准吗
  10. 程序员用python给了女友一个七夕惊喜!
  11. 2019最应该投资什么?是你明年的北大核心
  12. FFmpeg 快速上手:命令行详解、工具、教程、电子书
  13. fastjson使用toJSONString时null值不序列化问题
  14. 企业简介和公司介绍快闪PPT模板
  15. 用python输入 菱形
  16. 赋能IT,创新价值——互联网时代企业IT转型实践
  17. 基于SpringBoot 适合学习的开源社区平台
  18. 【论文汇总】2D目标检测文章汇总,持续更新
  19. Python批量格式化文件名为000xxx.jpg
  20. 优信拍集团php面试题_【优信拍PHP面试】一次不错的面试体验,面试官根据实际出发提问。-看准网...

热门文章

  1. conscience ['kɔnʃəns]
  2. 袁绍困局与张朝阳的雄心——类微博的狐友能让搜狐重回主战场吗
  3. 软件维护费 一般按照多少比例_2021年国家电网招聘:笔试进面一般是多少比例?...
  4. aqs java_Java并发之AQS详解
  5. 根据测试路径自动生成测试用例_自拍教程75Python 根据测试用例选择测试资源
  6. 扩展虚拟机Ubuntu磁盘容量空间
  7. Python定义函数的三种方式
  8. Python :给类或者类的对象添加打印内容 def __repr__(self)
  9. Python中有用的字符串方法
  10. linux 用户进程结束后 malloc申请的内存会自动释放吗,进程退出后malloc的内存是否会被释放?