基于matlab心电信号分析.docVIP

燕山大学 课 程 设 计 说 明 书 题目：基于matlab的心电信号QRS波检测与分析 燕山大学课程设计（论文）任务书 院（系）：电气工程学院 基层教学单位：生物医学工程系 学 号 学生姓名 专业（班级） 08医疗仪器 设计题目 基于matlab的心电信号QRS波检测与分析 设 计 技 术 参 数 阈值=0.6相对幅值的差 RR间期1mv 设 计 要 求 能够准确的提取与分析QRS波 绘制图表来说明检测结果 工 作 量 学会使用matlab软件 了解人体生理信号的特征 对心电信号有个初步的认识与掌握 搜集一组正常人体心电信号数据，编写matlab程序，并对它进行分析与检测 根据设计结果写出报告 工 作 计 划 查资料，搜集关于此次课设的相关内容 对题目进行分析，确定合理的解决方案 编程并且调试 写报告 参 考 资 料 周辉 数字信号处理基础及Matlab实现 中国林业出版社 2005 肖伟 刘忠 Matlab程序设计与应用 清华大学出版社 2005 指导教师签字 基层教学单位主任签字 说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份 2011年11 月10 日 目录 第1章 概述……………………………………………………………………4 第2章 ECG特征参数及分析…………………………………………………5 2.1 心电信号的特点…………………………………………………5 2.2 心电信号的特征参数……………………………………………5 第3章 QRS波得检测与分析…………………………………………………6 3.1 以软件为主的方法实现QRS波的检测应用计算机分析心电信号，已经越来越广泛的用于心脏功能检查(Holter系统)、心电监护等方面，而心电分析中的首要的关键问题是QRS波的检测可靠的检测不仅是诊断心律失常的重要依据，而且只有在QRS波确定之后，有可能计算心率并进行心率变异分析，才能检测ST段的参数和分析心电的其它细节信息。进行全面综合分析，才能对心脏的功能结构做出正确的判断。 图2.1 如图2.1为完整的心电信号波形图，分别由P波、P—R段、P—R间期、QRS复合波、S—T段、T波和U波组成。本文将重点讲诉QRS波。并且对 它进行相关分析。 QRS复波。代表两个心室兴奋传播过程的电位变化。由窦房结发生的兴奋波经传导系统首先到达室间隔的左侧面，以后按一定路线和方向，并由内层向外层依次传播。随着心室各部位先后去极化形成多个瞬间综合心电向量，在额面的导联轴上的投影，便是心电图肢体导联的QRS复合波。典型的QRS复合波包括三个相连的波动。第一个向下的波为Q波，继Q波后一个狭高向上的波为R波，与R波相连接的又一个向下的波为S波。由于这三个波紧密相连且总时间不超过0.10秒，故合称QRS复合波。QRS复合波所占时间代表心室肌兴奋传播所需时间，正常人在0.06～0.10秒之间 第三章 QRS波得检测与分析 3.1 以软件为主的方法实现QRS波的检测以软件为主的方 法实现QRS波的检测滤波之后的信号一般经过一些变换以提高QRS波的份量，进而采用一系列阈值进行判别，这些阈值有固定阈值法，也有可变阈值法。前者由于可能的干扰或高P、高T波的存在，若其滤波后超过其阈值便会产生假阳性(FP，falsepositive)结果；另外，当心律失常或QRS波幅度变小，阈值设置过高，会导致漏检产生假阴性(FN，falsenegative)结果。由于固定阈值的这些缺点，有研究者提出了用可变阈值检测，以提高检测的精确率，所采用的可变阈值包括幅度阈值、斜率阈值和时间间隔阈值等。clear all; clc; load testsig1; input=testsig1(1:2000); rate=testsig1(2001); sig=input; lensig=length(sig); wtsig1=cwt(sig,6,mexh); lenwtsig1=length(wtsig1); wtsig1(1:20)=0; wtsig1(lenwtsig1-20:lenwtsig1)=0; y=wtsig1; yabs=abs(y); %? sigtemp=y; siglen=length(y); sigmax=[]; for i=1:siglen-2 if (y(i+1)y(i)y(i+1)y(i+2))|(y(i+1)y(i)y(i+1)y(i+2)) sigmax=[sigmax;abs(sigtemp(i+1)),i+1]; end;