《虚拟仪器大作业实验报告》由会员分享，可在线阅读，更多相关《虚拟仪器大作业实验报告（17页珍藏版）》请在金锄头文库上搜索。
1、东南大学生物科学与医学工程学院 虚 拟 仪 器 实验报告大作业实验名称：基于MIT-BIH心率失常数据库的心电信号系统的设计专 业： 生 物 医 学 工 程 姓 名： 学 号： 同组人员： 学 号： 实 验 室: 综合楼716 实验时间： 2013/11/28 评定成绩： 审阅教师： 目录一实验目的二实验内容基于MIT-BIH心率失常数据库的心电信号系统的设计1.实验要求和说明2.程序设计流程图3.程序各版块介绍说明4.前面板的设计5.调试过程6.结果及分析三实验收获及小结四参考文献一实验目的现代医学表明，心电信号（ECG）含有临床诊断心血管疾病的大量信息，ECG的检测与分析在临床诊断中具有重
2、要价值，是了解心脏的功能与状况、辅助诊断心血管疾病、评估各种治疗方法有效性的重要手段。 本次大作业利用具有直观图形化编程和强大数字信号处理功能的虚拟仪器编程语言LabVIEW作为开发平台，设计一个基于虚拟仪器的简单心电信号分析系统，该系统具有心电信号的读取，处理分析，波形显示、心率显示及报警，波形存储和回放等功能。二实验内容1.实验内容及要求基于MIT-BIH心率失常数据库的心电信号系统的设计1. 本次大作业所用原始信号是从MIT-BIH（Massachusettes Institute of and Beth Israel Hospital，美国麻省理工学院和波士顿贝丝以色列医院）心率 数据
3、库（http:/www.physionet.org/physiobank/database/mitdb/）中选取心电信 号作为实验分析的数据。设计的系统要求对原始心电信号进行读取、绘制出 其时域波形，利用原始心电数据中的时间数据控制显示时间，并具有保存回 放功能，同时具有心率过快或过缓报警提示功能。2. 心电信号是微弱低频生理电信号，通常频率在0.05Hz100Hz，幅值不超过 4mV，它通过安装在皮肤表面的电极来拾取。由于实际检测工况的非理想， 在ECG 信号的采集过程中往往会受到工频噪声及电极极化等各种随机噪声 的影响。噪声的存在降低了诊断的准确性。其中影响最大的是工频干扰和基 线漂移噪声
4、。因此，在ECG 信号检测过程中，如何抑制工频干扰和基线漂 移等是必须解决的问题。要求选择并设计合适的滤波器，除去所给心电信号 的工频干扰和基线漂移。3. 检测心率：检测信号心电的R波，计算平均心率和实时心率（R-R波时间间隔 的倒数），并显示实时心率和平均心率。4. 对任一路心电信号滤波前后的信号进行时域分析和频谱分析，分别显示出结 果。提高部分：1. 由于原始信号数据并不是等间隔采样而得到的，而Labview中用数字滤波器处 理的数据要求等间隔的，由此需要对原始数据做一次线性插值处理，使其成 为均匀数字信号，以便后面的信号进行频域谱分析。提示：根据原文件心电 信号的特点，时间间隔设置成0.
5、001s，采用一次线性插值。2. 对所给的心电信号进行特征点的提取。各波时距的测量，P-R间期 Q-T间期的 测量。能显示心电信号特征参数。3. 使用Report Generation Toolkit for Microsoft Office 设计一个报表生成程序， 要求报表中显示标题栏，测试数据、和心电波形图，包含监测开始时间、结 束时间、监测总时间长度、监测期间平均心率、心率最大值、心率最小值、 监 测期间实时心率的波形图。4. 为了保护用户信息安全，登陆系统时要求设置用户登录模块。当用户运行程 序时，程序会自动弹出需要输入用户名和密码的口令信息对话框，输入不正 确则不能进入运行状态，此外
6、，还有相应的注册用户模块，删除用户模块， 修改密码模块等。2.程序流程图输入用户名和密码进行登录 否 口令错误，重新输入口令登录 口令正确进行数据回放是否进行数据回放 是 否从网站上下载的dat文件中读取数据读取数据之后进行数据处理，将两个通道的数据分开并且整理成12位否将数据转化成十进制进行处理 未处理完毕 数据处理完毕数据处理完之后用波形图和滤波前频谱分析进行显示进行滤波处理，去除基线漂移和工频干扰进行滤波后的频谱分析和显示 存储数据进行滤波后的波形显示 进行数据存储记录计算和显示平均心率，瞬时心率吧判断心率是否有异常3.程序各版块介绍说明数据的读取和处理模块1.首先是读取数据文件部分，先
7、添加一个文件对话框，并对文件类型和文件打开提醒进行设置，然后加入一个打开文件函数，设置其权限和操作分别为和。前期工作完成之后，就可以用二进制文件来读取数据了。2.关于数据的处理部分，第一道导联12位数据的低8位存在第一字节中, 高4位存在第二字节的前4位中; 第二道导联12位数据的低8位, 存在第三字节中, 高4位存在第二字节的后4位中。用索引函数将三个数据分别读取，然后对第二个数据进行处理，第二个数据与 进行操作后得到低四位，然后用拼接函数与第三个数据进行拼接就得到了第二道导联数据；第二个数据进行右移四位操作后保留高四位，然后用拼接函数与第一个数据拼接就得到了第一=导联数据。然后对得到的数据
8、变成十进制数进行处理。3.在读取数据时，由于是三个数据一处理，所以建立的for循环的循环次数应该是数据长度/3.滤波模块1.因为基线漂移部分都是频率非常小的波，我们选择了Butterworth的带通滤波器，高截止频率为100HZ，低截止频率为0.5HZ,滤除了原本波形中0.05HZ0.5HZ的部分，这样可以较好的去除基线漂移。2.为了去除工频干扰部分，在带通滤波器后面加入一个等波纹带阻滤波器，根据实验要求设置其高通带截止频率和低通带截止频率分别为100和0.05，然后根据实验滤波部分的要求将低阻带截止频率和高阻带截止频率分别设置为59和61，这样就可以滤去工频干扰。3.在已经滤去了基线漂移和工
9、频干扰的情况下，为了达到更好地滤波效果在等波纹带阻滤波器后面加入一个3阶切比雪夫滤波器，将波形滤除的更加完整。4.因为实验说明中采样频率为360.所以我们的滤波器的采样频率一率采取了360，这样滤波效果更明显一点。频谱分析模块 1.由于要进行滤波前后的频谱分析对比，所以频谱分析模块要在滤波前后都加入。2.首先要把处理后的数据表现成波形形式，用波形生成函数进行波形生成，并对其接线端子进行相应设置。3.在波形生成之后加入两个频谱函数幅度谱相位谱函数和FFT实部虚部函数，并在后面进行相应的显示波形图连接。4.滤波后的频谱分析与上面一样，波形图生成之后，将两个频谱图放在一起进行对比分析。数据记录与回放
10、模块1.首先在整体框图外面建立一个case结构，case结构的条件判断为。条件为真代表回放，条件为假代表不回放。执行程序时先选择条件假，将数据写入文件；然后再执行条件真，回放数据波形。2.将写入文件的程序框图放入条件为假的结构中，将读取文件的程序框图放入条件为真的结构中。3.写入文件时，首先将动态数据转换成数组，然后将X和Y分别写入一个一维的电子表格文件，读取时只要用电子表格文件读取再转换成动态数据直接连接到波形图上显示就可以了。4.需要特别注意的是，在提醒保存文件位置的对话框出来之后一定要按顺序分别存入x1,y1,x2,y2。读取的时候提醒读取文件位置的对话框出来之后也一定要按顺序选择，否则
11、波形会错乱。显示心率和判断心率异常模块1.首先用波峰检测函数来找出波峰个数，然后将每个波峰的位置传给索引数组函数，然后将找到的波峰个数减一后送入索引位置，这样就可以输出最后一个波峰的横坐标（即所有RR间期的总长度），然后除以波峰的个数就可以的算出一个RR间期的长度。2.用RR间期的长度除以采样周期就可以的得到一个RR间期所需时间，然后求出平均一分钟内出现了几个RR周期（即平均心率）。3.求瞬时心率时，建立一个输入变量，选择所需时刻的心率。然后求出这一时刻离哪一个波峰最近，求出这一波峰的位置和上一个波峰的位置进行相减姐可以得出瞬时RR间期的长度，然后同上面一样求出瞬时心率。4.用输入控件选择心率
12、上限和心率下限，然后将平均心率进行比较，如果在上下限范围之内则显示心率正常，否则显示心率异常。登陆（用户名和密码）模块1.首先建立一个while结构，接下来在while结构里面设置一个条件结构的嵌套。2.外层的条件结构用来判断是否登录，内层的条件结构用来判断用户名和密码是否正确。3.建立局部变量字符串和访问确认，为后面的用户名口令正误做准备。4.建立两个字符串输入变量，分别用作用户名和口令的输入，在外层的条件结构中加入两个字符串常量，设置为你自己的用户名和口令，如果输入变量中的用户名和口令与字符串常量设置的用户名口令相匹配，那么内层条件为真，可以继续执行下面的程序；如果不匹配则外层条件结构中的
13、字符串输出变量则会显示“口令错误”。5.这里的注意事项是，每启动一次程序就必须验证一次用户名和口令，并不是只要打开程序执行时输入用户名口令登陆就可以了，只要程序停止即使未关闭下次运行就必须再次登陆（但是用户名和口令不必再次输入，只要点登录就可以了）。4.前面板的设计说明1.首先因为本次实验前面板要显示的图形比较多，所以首先在前面板上建立几个选项卡。2.第一个选项卡放登录界面、数据显示和第一通道波形显示及滤波后波形；第二选项卡放入第二通道波形及滤波后波形；第三选项卡放入第二通道波形滤波前的幅频分析显示；第四选项卡放入第二通道波形滤波后的幅频分析显示；最后一个选项卡放入回放的第一通道和第二通道波形。3.登陆界面：建立两个字符串输入控件，分别用来输入用户名和口令；在建立一个字符串显示控件：当
