arma模型谱估计matlab_基于机器学习的心律失常分类(四)——心电信号特征提取[MATLAB]...

36fc3682414b9f36f931a1804e670392.png

目前比较常用的特征提取是提取心电信号的各波形间期长度、波峰高度等,本文是使用ARMA模型对心电信号进行处理,使用其系数来作为特征。

一、心拍划分

大多伴有异常波形的心律失常信号,通常都会具体表现在单个心拍中,需要对信号进行心拍划分。心拍信号截取的常见方法之一是以R波位置作为定位点,前后分别截取估计的采样点数。

491a4f9a9650bdf017567a8fdb05f57e.png
R波位置
上图是用matlab找出的R波位置,下面附上代码。但是我直接用的注释文件中标注的R波位置。
%ecg去噪后标记QRS波、R点。

啦啦啦

本文采用的MIT-BIH心律失常数据库以及其中由权威专家标注的R波位置。

e6dcda67afe33b5a38d49d0eb815315a.png
100.dat前十秒R波标注

心律失常的波形异常主要出现在PR间期和QRS波群,本文以R波位置为中心定位点,为了能基本包含所有的波形特性,向前向后各截取150个样本点共计301个样本点作为一个样本心拍。

%去噪后信号&标注信息

本文从MIT-BIH数据库的48条记录种的MLII通道数据中进行心拍划分,总计获取96871个样本心拍,其中正常心拍(N)有74466个,左束支传导阻滞心拍(L)有8051个,右束支传导阻滞心拍(R)有7239个,室性早搏心拍(V)有7115个。

二、ARMA模型

744a8ba4ca7a809e3a0da0ade7e34f13.png

三、特征提取

在开始设计和使用ARMA模型前首先需要通过定阶来确定阶数,常见的进行定阶的方法有:通过一个序列的自相关函数(ACF)和偏自相关函数(PACF)进行定阶、通过赤池信息量准则(AIC)进行定阶。

根据葛丁飞等人的利用MDL和 MDL敏感性函数作为准则进行模型定阶,或者根据毛雪岷等人所提出的按照信噪比和赤池信息量准则进行定阶,结果都可以表明当 n = 4,m = 2 时,模型拟合得到的信号与原始心电信号最为相似。因此本文选用ARMA(4,2)模型,通过R语言的arima()函数来建立模型,输出拟合的模型系数。

RStudio:

library(readxl)
X100 <- read_excel("D:/DATA/test1.xlsx", col_names = FALSE)
X=X100
N=nrow(X)ARMAnum=matrix(nrow = N,ncol = 6,byrow = TRUE)for(i in 1:N){TE=X[i,1:301]TE=t(TE)TE_AR=arima(TE,order=c(4,2,2))TE_CO=TE_AR$coeffor (n in 1:6) {ARMAnum[i,n]=TE_CO[n]}
}View(ARMAnum)write.csv (ARMAnum, file ="C:UsersDesktop数据2_处理后数据3_ARMA系数test1.csv")

本文从划分出的心拍中选取四类心拍各1700个,共计6800个样本心拍数据。

6d0d90b09762704f54e8688e5a545e5b.png
部分正常心电信号样本心拍的ARMA模型拟合系数

然后就可以用ARMA系数作为分类模型的输入进行训练了。

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

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

相关文章

抖音创作规范_抖音短视频运营小规则

抖音短视频规则一、抖音的定义1、抖音是个去中心化的平台&#xff0c;根本没有什么权重&#xff0c;主要是看你账号的活跃度&#xff0c;互动关系&#xff0c;和行为习惯标签&#xff0c;抖音后台技术有去重机制&#xff0c;机器会自动判断你的视频重复率。重复率达到一定的比例…

目睹DNS怪现状 本地对DNS的设定故障分析

分享一下我老师大神的人工智能教程。零基础&#xff01;通俗易懂&#xff01;风趣幽默&#xff01;还带黄段子&#xff01;希望你也加入到我们人工智能的队伍中来&#xff01;https://blog.csdn.net/jiangjunshowDNS&#xff08;Domain Name System&#xff09;在tcp/ip网络中&…

P1242 新汉诺塔

题目描述 设有n个大小不等的中空圆盘&#xff0c;按从小到大的顺序从1到n编号。将这n个圆盘任意的迭套在三根立柱上&#xff0c;立柱的编号分别为A、B、C&#xff0c;这个状态称为初始状态。 现在要求找到一种步数最少的移动方案&#xff0c;使得从初始状态转变为目标状态。 移…

python文件输入符_python文件IO与file操作

1 标准输入输出IO - (1) 打印到屏幕 print() print(self, *args, sep’ ‘, end’n’, fileNone): 把传递的表达式 转换成一个 字符串表达式 known special case of print """ print(value, ..., sep , endn, filesys.stdout, flushFalse) Prints the values to…

和县机电工程学校工业机器人_【校企合作】学校举行工业机器人教学系统捐赠仪式...

点击上方蓝字关注我们12月22日下午&#xff0c;我校与广东诚飞智能科技有限公司在公共实训中心会议室举行了“工业机器人教学系统”捐赠仪式。广东诚飞智能科技有限公司总经理杨成芳、区域经理杨少杰&#xff0c;学校副校长罗红金、教务科科长冯华以及机电专业组教师代表参加了…

此地钱多人傻速来

分享一下我老师大神的人工智能教程。零基础&#xff01;通俗易懂&#xff01;风趣幽默&#xff01;还带黄段子&#xff01;希望你也加入到我们人工智能的队伍中来&#xff01;https://blog.csdn.net/jiangjunshow原文&#xff1a;http://www.wulujia.com/“此地钱多人傻速来”这…

TCP三次握手(待细研究)

xu言&#xff1a; 看到一张不错清晰的Tcp三次握手图&#xff0c;收藏 Initiator 发起人 Receiver 接收者 LISTENING 状态xx服务启动后首先处于侦听&#xff08;LISTENING&#xff09;状态。 ESTABLISHED 状态ESTABLISHED的意思是建立连接。表示两台机器正在通信。 CLOSE_WAIT…

ubuntu20.04下安装vmshare tools_Go语言的安装与配置-vscode

vscode是比较好用的&#xff0c;虽然goland自带插件但是收费&#xff0c;所以配置vscode虽然麻烦但是也值得。Step1.下载安装vscode和golangDownload Visual Studio Code - Mac, Linux, Windows​code.visualstudio.comThe Go Programming Language​golang.google.cnnext step…

小说中人称转换作用_长篇儿童小说《合欢》:一树红绒落马缨

浙江少年儿童出版社出版肖复兴的第二部长篇儿童小说《合欢》&#xff0c;因对儿童独立人格的自觉尊重与强调&#xff0c;以及对人物复杂内心世界的深度开掘与勘探&#xff0c;回应了“儿童的被发现”这一重要课题。《合欢》不仅是作者对“远逝的童年”的一次深情回眸&#xff0…

内置函数之sorted,filter,map

# 4&#xff0c;用map来处理字符串列表,把列表中所有人都变成sb,比方alex_sb # name[oldboy,alex,wusir] # print(list(map(lambda i:i"_sb",name))) #5&#xff0c;用map来处理下述l&#xff0c;然后用list得到一个新的列表&#xff0c;列表中每个人的名字都是sb结尾…

asp按时间自动递增编号_约束力最强的手铐——美国ASP钢性手铐

美国ASP公司生产的钢性手铐&#xff0c;是目前为止世界上约束力最大的手铐。这种手铐的最初设计&#xff0c;是借鉴了中国古代的刑具——枷。钢性的手铐链接&#xff0c;非常好的约束罪犯的双手&#xff0c;使其没有任何活动的空间。而且最大的好处就是&#xff0c;执法人员只要…

python按行读取文件取消空白行_在Python中读取文件时忽略空行的最简单方法

我将堆栈生成器表达式&#xff1a;with open(filename) as f_in: lines (line.rstrip() for line in f_in) # All lines including the blank ones lines (line for line in lines if line) # Non-blank lines 现在&#xff0c;lines都是非空行。这样你就不用在电话线上打两次…

SQL Server2012如何导出sql脚本并且还原数据库

一 备份数据库 1 选择某一个数据库&#xff0c;右键依次选择&#xff1a;任务》生成脚本&#xff1b; 2 选择要编写脚本的数据库对象&#xff0c;注意此处可以选择特定的数据库对象&#xff0c;我们可以选择我们需要备份的数据表。 3 在当前页面下面第一步选择高级选项&am…

脚本启动慢_Linux 常用运维脚本,建议收藏

同步本机时间ntpdate 210.72.145.44清除系统缓存&#xff0c;空出更多内存free && sync && echo 3 > /proc/sys/vm/drop_caches && free杀掉僵尸进程kill $(ps -A -ostat,ppid | awk /[zZ]/ && !a[$2] {print $2})显示全部arp解析tcpdump a…

eclipse maven项目 class类部署不到tomcat下_Servlet tomcat部署

网上搜到的很多利用eclipse结合tomcat开发servlet的教程都要修改server.xml感觉这种方式太粗暴了&#xff0c;server.xml最好是与应用无关的&#xff0c;这里比较推荐export war的方式进行部署先记录一下环境配置过程我的系统是ubuntu 14.04&#xff1b;下载eclipse for javaee…

JS数组遍历

### 三 数组的遍历( 1 ) for 循环 let ary [a, b, c, d, 12, 3, 4, f]; for (let i 0; i < ary.length; i) { console.log(ary[i]); }( 2 ) for of &#xff08;ES6 &#xff09; let ary [a, b, c, d, 12, 3, 4, f]; for(let item of ary){ console.log(…

python多线程并发写入_Python多线程并发的简单测试

之前也写了一些简单的Python程序&#xff0c;对于多线程的并发一直没有涉及&#xff0c;今天决定先突破一下&#xff0c;把这个部分的内容先快速的掌握&#xff0c;然后在这个基础上细化改进。我的好友冰川擅长Python技术&#xff0c;所以就拿来主义&#xff0c;参考了他的文章…

tensorflow 指定cpu 但是还会运行 gpu_PyTorch VS TensorFlow 谁最强?这是标星 15000+ Transformers 库的运行结果...

点击上方“AI有道”&#xff0c;选择“星标”公众号重磅干货&#xff0c;第一时间送达作者 | Lysandre Debut译者 | 陆离出品 | AI科技大本营(ID: rgznai100)【导语】自然语言处理预训练模型库 Transformers 实现了几种用于 NLP 任务的最先进的 Transformer 架构&#xff0c;如…

CSS 组合选择符

CSS 组合选择符 一、简介 组合选择符说明了两个选择器直接的关系. CSS组合选择符包括各种简单选择符的组合方式。 在 CSS3 中包含了四种组合方式: 后代选择器(以空格分隔)子元素选择器(以大于号分隔&#xff09;相邻兄弟选择器&#xff08;以加号分隔&#xff09;普通兄弟选择器…

python爬取知乎live_python爬虫——xml数据爬取

回想去年写毕业论文时候需要数据&#xff0c;虽然数据量不大&#xff0c;但是手工录入&#xff0c;也是颇费眼神。现在自学了Python&#xff0c;就拿曾经的需求当做案例吧&#xff01; 爬取目标&#xff1a;最近20年每个月的仔猪价格&#xff0c;生猪价格和猪肉价格 爬取软件&a…