Hubery-个人项目经历记录

研究生期间很有幸的进入到了崔老师的组,从此也就进入到了分析人体生理信号的领域,充满挑战的同时也充满了乐趣。借着CSDN整理一下近几年来参与的项目,这里蕴含着我各种美好的回忆,也作为一个展示自己的平台吧。
开始之前,必须真心的感谢崔兴然老师、黄鄂院士、彭仲康教授,遇到这几位老师真的是我一生的荣幸。也要感谢一下读研究生以来实习和工作的爱谛安克X(化名,在职呢哈哈哈)公司,他们提供了太好的平台供我发挥。

音乐与脑电

这个是毕设课题,所有学习的内容比较多,也因此音乐的诸多物理特征和算法都有所了解。因尝试建立音乐和脑电的关系,也不断学习了很多种非线性算法。这里不一一列举,只给出一些好玩的现象。

基音缺失(misssing fundamental)

这里院士提出了Holo的思想,从调制频率方向上解释了为什么确实基频后,听到的还是那个频率的音高,Holo算法也被应用于多种其他时序信号中。如用于分析脑电来判定阿尔兹海默症脑小血管症
在这里插入图片描述

MSE(多尺度样本熵)

彭仲康教授提出的算法,我个人通过研究时发现音乐的时序特征熵越高,对应的脑电的熵也越高。比如听Jazz和Mozart的一首歌,大尺度上的熵有很好的对应关系。
在这里插入图片描述
在这里插入图片描述

DFA(去趋势波动分析)

DFA也是彭仲康教授提出的方法,在研究莫扎特效应时我曾计算过K448的DFA指数,很接近与1,与大自然的声音(粉红噪声)接近。
在这里插入图片描述

音乐时大脑的功能连接

很有趣的结果,听音乐时大脑的连接强度是变小的(部分论文显示是升高的)。
在这里插入图片描述

助眠:ASMR与ASSR

曾通过分析一些有助眠效果音乐,还有ASMR(比如掏耳朵,揉纸团,海浪等声音),曾尝试使用音乐的ASSR听觉稳态诱发响应)来解释,但因数据缺失耽搁了,后续有机会再次收些数据试试。

睡眠与脑电(或PSG)

主要做的任务是睡眠分期呼吸暂停事件检测
具体可见连接:https://hubery-hu.blog.csdn.net/article/details/135231518
目前睡眠5分期效果已有86.5%,同时为集成转化了matlab,python,c++等多种版本。
在这里插入图片描述

情绪与脑电 & 运动想象

这里从运动想象出发,学习了其中最常见的CSP(共空间模式)算法,将其应用于少通道脑电的情绪识别中,结合EMD和SVM常见特征,最终8种情绪的分类效果在80%左右。
在这里插入图片描述

脑机接口(脑控)

最初的版本是在安克做的,使用了眨眼做特征实现脑控,但是考虑到其应用面较少,脑控的算法也不成熟,最终不了了之,目前自己购买了Muse的头环,打算先开发一个脑控的猜拳小游戏。有兴趣的小伙伴也可以联系我一起测试哦~
在这里插入图片描述

脑电与疾病监测

后来逐渐将脑电分析用在了各种疾病的监测上,比如癫痫(癫痫发作,尖棘波识别)麻醉深度的监测爆发抑制的识别等。下方是我模拟的爆发抑制结果(做疾病监测时开发的第一个算法)和两位专家标注结果的比较,很有成就感~
在这里插入图片描述

陀螺仪

欧拉角(用于空间音频)

这里基于TDK、ST等厂商的6轴陀螺仪芯片,应用四元数的方法实现了欧拉角的计算。下图是对Yaw(偏航角)的中间结果。
在这里插入图片描述

运动识别和姿势识别

做这个任务时,每天带着耳机蹦蹦跳跳的(测了自行车,跑步,走路,爬楼梯,跳绳,仰卧起坐,俯卧撑等),还是很怀念(I人被迫吸引别人注意力哈哈哈)~
下面分别是Jump和Walk的6轴值以及计算出的欧拉角,可以看出计步计数)还是很简单的,姿势识别起来也很简单(比如某些角度的规律变化?当然具体识别算法需融合更多信息,之后在其他博客中介绍)
有了欧拉角,一些低头提醒久坐提醒的功能也就很好开发了,根据场景配置即可。
在这里插入图片描述
在这里插入图片描述

音乐信号

除了常见的音乐特征(比如MFCC,色度频率,过零率等)特征的计算外,也做过一些好玩的尝试,比如改变相位和填充包络等,音乐风格真的很好玩,大家也可以试试哦。

改变相位

这里改变了K448的相位,这两首歌的频谱还是一样的。
在这里插入图片描述

填充包络

提取K448的包络后,波形不一致了,但整体节奏感还是一样的,音乐变得很魔性
在这里插入图片描述

其他生理信号

心电和HRV

心电的算法已经很常见了,这里不再赘述,
我曾使用HRV去计算放松度来和BrainCo的结果做比对,整体曲线大趋势一致,可惜量级不一致,且细节上毛刺变化差异很大。但这个模型需建立在大量数据集的基础上,后来也没有BrainCo头环了。。。也耽搁了,后续有兴趣再试试吧。

PPG信号和呼吸率

耳内采集的PPG信号在静息放松的时候还可以同时监测到呼吸率(参考EDR算法)。
从下图可以看出,PPG的信号在上下包络上也有周期性波动,在信号质量比较好的时候,也可以用它来计算呼吸率。
在这里插入图片描述

血压信号的SVV和PPV

具体可见链接:https://hubery-hu.blog.csdn.net/article/details/133696610
具体效果(这里只展示最好看的波形,实际算法考虑了多种异常情况,可看上方链接了解):
在这里插入图片描述
在这里插入图片描述

其他算法

去噪

小波去噪和EMD去噪:具体见链接:
效果如下:https://hubery-hu.blog.csdn.net/article/details/133862472
在这里插入图片描述

Hample算法:https://hubery-hu.blog.csdn.net/article/details/132020843
效果如下:
在这里插入图片描述

ICA

这里复写了Sobi和MARA算法,实现脑电的ICA。
下图展示了分解出的眼电成分(眨眼和眼瞟)。
在这里插入图片描述

各种深度学习方法

CNN、LSTM、Transformer、GNN等都有所涉及,不一一赘述了。

混沌、分形、滤波等多种算法

总结

希望这些介绍能让大家对我有个基础的映像,不定期更新哦~

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/318781.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

初始数字孪生

文章目录 概念定义 架构框图 基本要求 功能要求 服务应用 参考文献 概念定义 数字孪生(digital twin),是指具有保证物理状态和虚拟状态之间以适当速率和精度同步的数据连接的特定目标实体的数字化表达。 架构框图 数字孪生的整体架构&a…

MATLAB逻辑与流程控制

MATLAB逻辑与流程控制 文章目录 MATLAB逻辑与流程控制1、循环结构1.1、for ...end1.2、while ... end 2、分支结构2.1、if ...else2.1、if ...else ... end2.2、switch ... case ... end 1、循环结构 1.1、for …end 格式: 其中,步长的默认值为1,可以省…

STM32H5XX和STM32H7XX选型对比

文章目录 STM32H563/H573STM32H743/753对比内核不同H5独有安全管理器H7的外设资源更丰富 STM32H563/H573 STM32H563和STM32H573微控制器扩展了STM32高性能产品组合。这两款微控制器具有增强的性能和安全性、更高的能效和更多的片内外设。 STM32H563/573产品系列提供1至2 MB的…

查询运行的java程序线程总数

(1)方法一 ps -ef|grep java pstree -p pid|wc -l 安装pstree包,yum install psmisc (2)方法二 top -H -p pid

python 学生信息管理系统

Python项目 学生信息管理系统 带sql文件,适合初学者。 使用MySQL数据库。 送运行视频导入教程视频。 实现功能大概如下: 基本功能增删改查 限制成绩在各种情况下只能为0-100 实现成绩的统计功能(柱状图显示) 筛选各个成绩区间的学生信息 实现排序…

Leecode刷题:024反转单项链表

题目: 是否没看题解:有思路但是没做出来 最初的解答思路: 定义一个pre指向头结点,cur为pre.next,定义last为cur.next,将cur.next指向pre定位到cur位置,然后cur定位至last节点。 最开始的代码: class S…

Spark---RDD算子(单值类型Value)

文章目录 1.RDD算子介绍2.转换算子2.1 Value类型2.1.1 map2.1.2 mapPartitions2.1.3 mapPartitionsWithIndex2.1.4 flatMap2.1.5 glom2.1.6 groupBy2.1.7 filter2.1.8 sample2.1.9 distinct2.1.10 coalesce2.1.11 repartition2.1.12 sortBy 1.RDD算子介绍 RDD算子是用于对RDD进…

Python基础-05(输出输入、if、if else和elif)

文章目录 前言一、输出(print())和输入(input())二、if、if else、elif1.if2.if else3.关于输入input的默认值4.elif 前言 今天复习一些非常基础的内容,以及if、if else和elif语句 一、输出(print()&…

雪花主键_在postgres的实现

一、主键的属性: 1. 全局唯一,不重复 2. 趋势有序,后面的值大于前面的值 3. 高性能, 读写效率高。qps不能太低,否则容易造成线程堵塞 4. 可用性好,业务对ID的可用性高,不存在单点故障 二、常…

漏洞复现-任我行CRM系统SmsDataList接口SQL注入漏洞(附漏洞检测脚本)

免责声明 文章中涉及的漏洞均已修复,敏感信息均已做打码处理,文章仅做经验分享用途,切勿当真,未授权的攻击属于非法行为!文章中敏感信息均已做多层打马处理。传播、利用本文章所提供的信息而造成的任何直接或者间接的…

3DMAX 中的 VR 渲染器如何设置局部渲染?

首先我们要得打开渲染设置,在3damx里按F10,调出渲染设置。选定渲染器为Vary渲染器: 设置VR的局部渲染,需要打开帧缓冲,我们在V-ary项下,打开帧缓冲(点击显示VFB): 然后点击茶壶按钮(…

如何打开wps的备份中心查找备份文件

备份中心在我们使用WPS Office时扮演着重要的角色。经常保存文件的同时,我们也应该学会备份文件,以免意外损失。本文将向您介绍如何使用WPS备份中心来查找并恢复备份文件,方便您在需要时快速找到所需文件。 图片来源于网络,如有侵…