下载

0下载券

加入VIP
  • 专属下载特权
  • 现金文档折扣购买
  • VIP免费专区
  • 千万文档免费下载

上传资料

关闭

关闭

关闭

封号提示

内容

首页 心电信号的计算机分析final

心电信号的计算机分析final.doc

心电信号的计算机分析final

何海得
2019-02-25 0人阅读 举报 0 0 0 暂无简介

简介:本文档为《心电信号的计算机分析finaldoc》,可适用于医药卫生领域

B陆建荣生命科学和生物工程学院心电信号的计算机分析【实验目的】:通过理论结合实际,用C语言编程对MIT心电信号数据进行分析,实现低通滤波、高通滤波、QRS检测、特征提取、心律失常分析,从中了解和掌握数字信号处理的方法和应用。【实验要求】读取数据QRS检测特征参数提取心率失常分析功率谱分析【实验报告】一实验介绍心脏在有节律的活动过程中,能在人体表面产生微弱的电信号,如果我们在人体表面的特定部位安放电极,就能在电极上获得微弱的心电信号,此信号经放大、处理后,描记在记录纸上就是心电图,它能够反映心脏的功能及病情。在获取心电图的过程中,由于心电信号比较微弱,仅为毫伏(mV)级,所以极易受环境的影响。对心电信号引起干扰得主要因数有:工频干扰、电极接触噪声、运动伪迹、呼吸引起的基线漂移和心电幅度变化、信号记录和处理中电子设备产生的干扰、电外科噪声等。为了增强心电信号中的有效成分,抑制噪声和伪迹,提高波形检测准确率,除了对心电记录仪的硬件抗干扰能力有较高的要求外,心电信号AD变换后的处理也至为重要。用于心电信号数字处理的方法主要有:消除电源干扰的工频滤波器,消除采样时间段引起信号失真的汉宁平滑滤波器,消除高频肌电的低通滤波器,消除直流偏移和基线漂移等低频噪声的高通滤波器,以及用于QRS波检测的带通滤波器。本实验利用MIT心电信号数据库,简单设计了对心电信号进行计算机分析的实验,实验主要分成两部分:信号处理和心电参数分析信号处理的方法有低通滤波、高通滤波、微分(查分运算):,对处理后的信号进行如下分析:QRS检测心率失常分析参数提取功率谱分析。本实验的整个过程是:先读取文件数据,将数据显示在计算机屏幕上,并可进行翻页显示,然后对所读心电数据进行低通滤波、高通滤波、微分(查分运算)等处理,同时将处理后的数据显示在屏幕上对心电信号的分析是采用处理后的的数据,先对QRS波进行检测,然后计算特征参数,然分析是采用处理后分析是否心率失常,此实验还对信号进行了功率谱分析,参见图。图:心电信号计算机分析二实验的理论与实现此实验的是用C编程来实现整个处理和分析过程,整个实验的分析过程是:先对信号进行低通滤波、高通滤波和微分然后再进行QRS检测、参数提取和心律失常分析(如图),这样可以滤除干扰,保证分析的可信度。图:心电信号计算机分析流程程序有三个功能模块:数据操作、信号处理和心电信号分析。数据操作主要对数据文件进行读写操作,显示波形并可浏览波形信号处理对I信号进行低通滤波、高通滤波和微分心电信号的分析是对心电信号进行QRS检测、参数提取、心率失常分析,另外还对心电信号进行了功率谱分析,处理和分析程序请参见附页。下面分别介绍各个不同步骤。读取数据并进行显示打开选择的数据文件,将心电数据读入内存中,并安点显示的方式以数据的大小作为幅度,点的先后次序作为时间轴显示在活动窗口通过操作刷新按钮来变化显示数据段。显示结果如图。信号处理()低通滤波:此实验采用二阶低通滤波器,设计的传递函数为:)()()(=zzzHY(Z)X(Z)=z^z^()z^()z^()y(n)y(n)y(n)=x(n)x(n)x(n)y(n)=y(n)y(n)x(n)x(n)x(n)它的差分方程是:)()()()()()(TnTxTnTxnTxTnTyTnTynTy=绝止频率为:Hz。显示结果如图()高通滤波高通滤波的设计方法是将一个全通滤波器减去一个一阶低通滤波器。所采用的低通滤波器是一个整数型滤波器,其传递函数为:)()()(==zzzXzYzHlp其差分方程为:)()()()(TnTxnTxTnTynTy=从原始信号中减去除以直流增益的低通滤波器的输出,即可得到高通滤波。他的传递函数为:)()()()(zHzzXzPzHlphp==其差分方程为:)()()()()(TnTxnTxTnTyTnTxnTp=此高通滤波器的绝止频率为Hz,显示结果如图。()微分对心电信号进行微分可以得到QRS复波的斜率信息。我们采用点微分,它的传递函数为:)()(=zzzzH差分方程为:)()()()()(TnTxTnTxTnTxnTxnTy=结果显示如图。图读取数据显示波形图低通滤波后显示波形图高通滤波后心电波形图差分运算结果心电信号分析心电信号能反映心脏的活动功能,正确分析心电信号对于疾病的诊断和预防具有重要的意义,心脏的活动每时每刻都在进行,对病人实际发病情况研究表明,发病和心电信号异常一般时间歇发生,如要依靠医生在不间断的心电信号中发现异常信号是一件不容易做到的事,随着计算机技术的发展,利用计算机对心电信号进行分析能减轻医生的工作强度,同时业能提高捡出率。本实验中对MIT数据库中的信号进行了简单的检测和分析,工作主要有QRS检测、参数提取、心率失常分析及功率谱分析。下面分别介绍。()QRS检测:QRS波是心电信号中具有明显特征的波形,通常用它来计算心率,同时他也能反映心脏收缩期的特征,检测QRS波在心电信号分析过程中至关重要。QRS扽检测过程大体是:原始心电信号先进型低通滤波、高通滤波,在进行微分,然后把微分后的信号去平方,在用移动积分窗框破进行积分,最后用阈值进行比较,就得到QRS复波峰值点,参见图。实现程序参见所附程序。图QRS检测流程检测结果参见图。()参数提取:对QRS间隔时间进行计算,得到RR间期,由RR间期计算心率,将临近个值进行平均得出平均心率值。另外还计算了QRS间期。结果参见图。()心率失常分析:QRS宽度及RR间期都是由RS检测器所确定的。根据这两个参数就可将ECG信号进行分类。在这个以两个参数为基准的得到映射图中(参见图),算法通过对一套由临床医生确定的个QRS复波进行识别来建立一个称为“正常”的区域。这些QRS复波包括特定病人的“正常”心率与形态。这一识别过程在两组维映射空间内,建立一个初始的“正常”区域的中心。除了“”区域之外,对映射图中其他区域的边界均计算出距离常规区域中心位置(定位)的百分比。“”区域的固定边界是根据生理指标(极限)而定的。所有被映射到“”区域的指示点都被认为是噪声,因为它不在通常所希望的最小的RR间期或QRS宽度的生理限度之内。结果参见图。()功率谱计算:心电信号可以看作为周期信号,假定他以N为周期,则其DFT将在以采样频率为宽度的各区间内重复。有定律可知:)()(∑∑=∞==NknkXNnx如果对一个函数c(t)进行采样,并用FFT法计算他的DFT,则会得到:∑==)(NnNjknkencCπ这时对N个频率点的功率谱的周期图估算定义为:)(CNP=图基于RR间隔和QRS间期映射的心率失常算法,,,,()(==NkCCNkPkNk)(NCNNP=功率普计算结果见图。三讨论本实验使用C语言编程实现了对心电信号的处理和分析,实验中所采用的方法多出自于书中内容,方法简单,易于理解,但实际使用时对于不常见病人的心电信号可能会出现误分析的结果。现在的数学理论在不断发展,神经网络、小波变换、最优化控制等都可用于心电信号的处理和分析,新理论的应用定能极大地提高心电信号的处理效果和分析的正确率。图QRS检测结果显示图心电信号特征参数提取图心率失常分析结果图功率谱图陆建荣年月日

用户评价(0)

关闭

新课改视野下建构高中语文教学实验成果报告(32KB)

抱歉,积分不足下载失败,请稍后再试!

提示

试读已结束,如需要继续阅读或者下载,敬请购买!

评分:

/12

VIP

在线
客服

免费
邮箱

爱问共享资料服务号

扫描关注领取更多福利