基于MATLAB的心电信号分析系统的设计与仿真设计matlab.docxVIP

据库共有48个病例，每个病例数据长30min，总计约有116000多个心拍，包含有正常心拍和各种异常itle(高通滤波后信号频谱{9.997})subplot(3,1,3)plot(f2,abs(t(2,2,3)5/20plot(t,F)title(差值后信号时域波形 据库共有48个病例，每个病例数据长30min，总计约有116000多个心拍，包含有正常心拍和各种异常 itle(高通滤波后信号频谱{9.997})subplot(3,1,3)plot(f2,abs( t(2,2,3)5/20plot(t,F)title(差值后信号时域波形)axis([02.5- butt(OmegaP1,OmegaS1,Rp1,As1);[b1,a1]=impinvr(cs1, 摘要： 本文是利用 MATLA件对美国麻省理工学院提供的 MIT-BIH数据库的 122 号心电信 号病例进行分析，利用 MATLAB软件与 simulink 平台对 122 号心电信号的病例进行读取、插 值、高通滤波、低通滤波等的处理。将心电信号中的高频和低频的杂波进行滤除后对插值前 后滤波前后的时域波形与频谱进行分析。同时也将滤波器的系统函数进行读取，分析，画出 滤波的信号流程图，也分析各个系统与级联后的系统的冲击响应、幅频响应、相位响应和零 极点图来判断系统的稳定性，并用 MATLA件将图形画出，以便于以后的对系统进行分析。 关键词： MATLA，Bsimulink ，心电信号，滤波器 1. 课程设计的目的、意义： 本设计课题主要研究数字心电信号的初步分析与滤波器的应用。通过完成本课题的设计，拟 主要达到以下几个目的： (1) 了解 MATLA件的特点和使用方法， 熟悉基于 Simulink 的动态建模和仿真的步骤和 过程； (2) 了解人体心电信号的时域特征和频谱特征； (3) 进一步了解数字信号的分析方法； (4) 通过应用具体的滤波器进一步加深对滤波器理解； (5) 通过本课题的设计， 培养学生运用所学知识分析和解决实际问题的能力。 2 设计任务与技术指标： 设计一个简单的心电信号分析系统。对输入的原始心电信号，进行一定的数字信号处理，进行频谱分 析。 采用 Matlab 语言设计，要求分别采用两种方式进行仿真，即直接采用 Matlab 语言编程的静态仿真方 式、采用 Simulink 进行动态建模和仿真的方式。根据具体设计要求完成系统的程序编写、调试与功能测 试。 试 2.1 必做部分： 2.1.1 读取原始心电信号 美国麻省理工学院提供的 MIT-BIH数据库是一个权威性的国际心电图检测标准库， 近年 来应用广泛，为我国的医学工程界所重视。MIT-BIH数据库共有 48 个病例， 每个病例数据长 30min ，总计约有 116000 多个心拍，包含有正常心拍和各种异常心拍，容丰富完整。 为了读取简单方便， 采用其txt 格式的数据文件作为我们的原心电信号数据。利用 Matlab 提供的文件 textread 或 textscan 函数，读取 txt 数据文件中的信号，并且还原实际波形。 2.1.2 对原始心电信号做线性插值 由于原始心电信号数据不是通过等间隔采样得到的，也就是说原始的心电数据并不是均匀的， 而用 1 / 20 bplot(3,1,1);plot(H1);title(subplot(3,1,2);plot(H2M=length(F0);m=1:M;fs=1000;f2=m*fs/M;F01=fft(F0);Y.2596B=1.0e+028*Columns1through10000000000005.3949系统首先应对原始心电信号做线性插值处理，使其变为等间隔的数字信号，否则直接处理后会出现偏差，根据心电i i 1 ，i i 1 ， j bplot(3,1,1);plot(H1);title(subplot(3,1,2);plot(H2 M=length(F0);m=1:M;fs=1000;f2=m*fs/M;F01=fft(F0);Y .2596B=1.0e+028*Columns1through10000000000005.3949 系统首先应对原始心电信号做线性插值处理，使其变为等间隔的数字信号，否则直接处理后会出现偏差，根据心电 i i 1 ， i i 1 ， j j 1 j 1 i 1 j 1 i 1 j j 分隔成 0.001s 。添加的幅值点采用一次线性插值。对二维数据进行插值, 相连幅值间数据的插值根据时间 进行, 运算公式如下: t t t N t