基于MATLAB的心电信号分析系统的设计与仿真VIP免费

信号与线性系统课程设计报告课题二基于MATLAB的心电信号分析系统的设计与仿真摘要：本文是利用MATLAB软件对美国麻省理工学院提供的MIT-BIH数据库的122号心电信号病例进行分析，利用MATLAB软件及simulink平台对122号心电信号的病例进行读取、插值、高通滤波、低通滤波等的处理。将心电信号中的高频和低频的杂波进行滤除后对插值前后滤波前后的时域波形及频谱进行分析。同时也将滤波器的系统函数进行读取，分析，画出滤波的信号流程图，也分析各个系统及级联后的系统的冲击响应、幅频响应、相位响应和零极点图来判断系统的稳定性，并用MATLAB软件将图形画出，以便于以后的对系统进行分析。关键词：MATLAB，simulink，心电信号，滤波器1.课程设计的目的、意义：本设计课题主要研究数字心电信号的初步分析及滤波器的应用。通过完成本课题的设计，拟主要达到以下几个目的：(1)了解MATLAB软件的特点和使用方法，熟悉基于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~信号与线性系统课程设计报告Matlab中提供的数字滤波器处理数据时，要求数据是等间隔的。因此设计的系统首先应对原始心电信号做线性插值处理，使其变为等间隔的数字信号，否则直接处理后会出现偏差，根据心电信号的特点,把时间分隔成0.001s。添加的幅值点采用一次线性插值。对二维数据进行插值,相连幅值间数据的插值根据时间进行,运算公式如下:1iittt，001.0/tN，1iiAAA，001.01jjtt，NAAAjj/1其中it是第i个数据时间点，Ai是与之对应的数据，N是两数据之间需要的插值数，A是需要插值的两点数据差，，，，，，，，，1321321Njarraysizei1111ijijAAtt，，1ji，时数组jjAt,依次排列,即得到了插值后等间隔的新数据。2.1.3根据心电信号的频域特征，设计相应的低通和带通滤波器一般正常人的心电信号频率在0.7～100HZ范围内，幅度为V10(胎儿)～5mV(成人)。人体心电信号微弱，信噪比小，因此，在采集心电信号时，易受到仪器、人体活动等因素的影响，而且所采集的心电信号常伴有干扰。采集心电数据时，由于人的说话呼吸，常常会混有约为0.1Hz到0.25Hz频段的干扰，对于这些低频干扰，可以让信号通过一个高频滤波器，低截止频率设置为0.25,来滤除低频信号，对于高频信号干扰，可以让信号再通过一个低频滤波器，其中截止频率设置为99Hz。也可以直接应用带通滤波器设计。（1）根据以上指标，设计模拟巴特沃斯（切比雪夫）低通、高通或带通滤波器，画出幅频特性（模拟滤波器幅频特性freqs）。（2）根据心电信号频谱范围设计一个3阶以上模拟滤波器对心电信号进行预滤波；（3）采用直接、级联或并联方式，实现该系统，并画出系统的信号流图；（4）分析系统的时域特性（阶跃响应、冲击响应等），并用Matlab绘出相关波形；（5）用Matlab分析幅频特性，并绘出相关波形；（6）分析系统函数零...