心电信号分析系统

信号与线性系统

课程设计报告

课题名称：基于Labview的心电信号分析班级:通信092班

姓名：孙星星

学号：092307

成绩：

指导教师：王宝珠

日期：2011年12月28号

目录

一绪论

1.1 摘要

关键词

1.2 课题简介

1.3 课题设计目的及意义

二设计任务及技术指标

三设计总体方案流程图

四方案论证与比较

4.1 心电信号读取设计方案论证与选择

4.2 线性插值的设计方案论证与选择

4.3 滤波器及谱分析的设计方案论证与选择

五课题设计具体内容

5.1心电信号的读取

5.2 对原始心电信号线性插值

5.3设计相应的低通高通滤波器及其谱分析

5.4 总体的程序框图和前面板

六课题设计遇到的问题及其解决办法

6.1心电信号的读取部分

6.2对原始心电信号线性插值部分

七结论

八课程设计的收获、心得体会与建议

参考文献

一绪论

1.1 提要:

心电信号分析系统是对心电信号文件拾取输入后，对其做一定的数字信号处理，进行频谱分析。

Labview是一种带有图形控制流结构的数据流模式（Date Flow Mode),程序执行是数据驱动，同时也是一种图形化的编程语言--G语言。本设计采用Labview 综合运用其丰富的VI库来实现心电信号的拾取、线性插值、滤波、谱分析。该课题利用VI库中索引数组、数组子集、字符串--数值转换、While循环、chebyshev 滤波器、单边FFT，双边FFT等，得到了简单的拾取、插值、滤波、谱分析等功能。，对心电信号做简单的数字信号处理。

关键词：Labview,心电信号，VI库，数字信号处理

1.2 课题设计目的及意义

本课题主要研究基于Labview的数字心电信号初步分析及其各种滤波器的应用。通过完成本课题的设计，了解基于LabVIEW虚拟仪器（一种软件化的仪器系统）的特点和使用方法，熟悉掌握LabVIEW的使用及练习使用其不同的功能，了解人体心电信号的时域特征和频谱特征，通过对心电信号的滤波处理、频谱分析，进一步了解数字信号的分析方法，进一步加深对各种滤波器（chebyshev,inverse chebyshev,butterworth)的理解。此外，通过本课题的设计，培养运用所学知识分析和解决实际问题的能力。

心电信号分析是一门比较实用的电子工程的专业课程。当今社会，心血管疾病是发病率和死亡率最高、对人类生命威胁最大的疾病。心电信号预处理就是对心电信号的时域特征、频域特征进行了解，以便以后对心电信号的自动识别起到一定的基础作用。另外，Labview具有强大的虚拟仪器功能和软件开发功能，运行速度快、兼容性和移植性好、方便易用，适合于课程设计短期内完成。

二设计任务及技术指标

课题所用信号是美国麻省理工学院提供的MIT-BIH数据库（一个权威性的国际心电图检测标准库），近年来应用广泛，为我国的医学工程界所重视。MIT-BIH 数据库共有48个病例，每个病例数据长30min，总计约有116000多个心拍，包含有正常心拍和各种异常心拍，内容丰富完整。

为了读取简单方便，采用其txt格式的数据文件作为我们的源心电信号数据。利用labvIEW提供的文件I/O函数，读取txt数据文件中的信号，并且还原实际

波形。

设计任务：设计一个基于虚拟仪器的简单的心电信号分析系统，对原始心电信号做输入后，对其做一定的数字信号处理，进行频谱分析，根据具体设计要求完成系统的程序编写及调试。

（1）对原始数字心电信号进行读取，由数字信号数据绘制出其时域波形。（2）对数字信号数据做一次线性插值，使其成为均匀数字信号，以便后面的信号分析。

（3）根据心电信号的频域特征（自己查阅相关资料），设计相应的低通和带通滤波器。

（4）编程绘制实现信号处理前后的频谱，做频谱分析，得出相关结论。（5）对系统进行综合测试，整理数据。

技术指标：一个高频滤波器，低截止频率设置为0.25,来滤波低频信号；对于高频信号干扰，可以让信号再通过一个低频滤波器，其中截止频率设置为99Hz。

三设计方案

开始

对心电信号进行拾取

While循环线性插值

chebyshev高通滤波器、

低通滤波器滤波处理

双边FFT、单边fft谱分析

四 方案论证与比较

如前所述，总体设计方案比较显而易见，现将各部分分述一下。 4.1 心电信号读取设计方案论证与选择

课题所用信号是美国麻省理工学院提供的MIT-BIH 数据库，此数据库的数据文件的前两行为解释说明文字，不是真正的信号数据，读取信号程序要能够自动忽略前两行文字，只读取真正的数字信号数据。另外labvIEW 默认的从文本文件中读取的数据都是字符串，因此在使用心电信号数据前需要将其转换为数值才可以。第一列时间数据均为0:00.007这种格式，因此需要将字符串0:00.007先转化为字符串0.007，即去除字符串中冒号（：）以前的部分，然后再将其转为数值。最后利用已经转为数值的分别代表心电信号时间和幅值的两个一维数组，图形化还原原始心电信号波形，在此推荐利用labvIEW 中XY Graph 。 4.2 线性插值的设计方案论证与选择

由于原始心电信号数据并不是由等间隔采样而得到的，也就是说原始的心电数据并不是均匀的，但是后面肯定要用到数字滤波器对心电信号做处理，而labvIEW 中用数字滤波器被处理的数据应该都是等间隔的，否则处理后会出现偏差，因此要对原始心电信号要做一个线性插值处理，使其变为等间隔的数据。根据原文件心电信号的特点, 把时间分隔成0.001s 。添加的幅值点采用一次线性插值。对二维数据进行插值,相连幅值间数据的插值根据时间进行,运算公式如下:

1--=∆i i t t t

001.0/t N ∆=

1--=∆i i A A A 001.01+=-j j t t

N A A A j j /1∆+=-

其中1,,,,1,3,2,1,,3,2,11111===-==----j i A A t t N j a r r a y s i z e i i j i j 时数组

j j A t ,依次排列,即得到了插值后等间隔的新数据。虽然Labview 提供的插值法有

很多，但是由于自身对Labview 不熟悉，所以采用教师提供的线性插值法（即上