课题二基于MATLAB平台的心电信号分析系统设计与仿真.pdfVIP

课题二 基 MATLABDE 的心电信号分析系统的设计与仿真 一、本课题的目的 本设计课题主要研究数字心电信号的初步分析及滤波器的应用。通过完成本课题的设 计，拟主要达到以下几个目的： (1)了解MATLAB 软件的特点和使用方法，熟悉基 Simulink 的动态建模和仿真的步骤 和过程； (2)了解人体心电信号的时域特征和频谱特征； (3)进一步了解数字信号的分析方法； (4)通过应用具体的滤波器进一步加深对滤波器理解； (5)通过本课题的设计，培养学生运用所学知识分析和解决实际问题的能力。 二、课题任务 设计一个简单的心电信号分析系统。对输入的原始心电信号，进行一定的数字信号处理， 进行频谱分析。采用Matlab语言设计，要求分别采用两种方式进行仿真，即直接采用Matlab 语言编程的静态仿真方式、采用 Simulink 进行动态建模和仿真的方式。根据具体设计要求 完成系统的程序编写、调试及功能测试。 (1)对原始数字心电信号进行读取，由数字信号数据绘制出其时域波形。 (2)对数字信号数据做一次线性插值，使其成为均匀数字信号，以便后面的信号分析。 (3)根据心电信号的频域特征（自己查阅相关资料），设计相应的低通和高通滤波器。 (4)编程绘制实现信号处理前后的频谱，做频谱分析，得出相关结论。 (5)对系统进行综合测试，整理数据，撰写设计报告。 三、主要设备和软件 (1)PC机一台。 (2)MATLAB6.5以上版本软件，一套。 四、设计内容、步骤和要求 4.1必做部分 4.1.1利用Matlab对MIT-BIH数据库提供的数字心电信号进行读取，并还原实际波形 美国麻省理工学院提供的MIT-BIH数据库是一个权威性的国际心电图检测标准库，近 年来应用广泛，为我国的医学工程界所重视。MIT-BIH数据库共有48个病例，每个病例数 据长30min，总计约有116000多个心拍，包含有正常心拍和各种异常心拍，内容丰富完整。 为了读取简单方便，采用其tx 格式的数据文件作为我们的原心电信号数据。利用Matlab 提供的文 textread或textscan函数，读取tx 数据文件中的信号，并且还原实际波形。 4.1.2对原始心电信号做线性插值 由于原始心电信号数据不是通过等间隔采样得到的，也就是说原始的心电数据并不是均 匀的，而用Matlab 中提供的数字滤波器处理数据时，要求数据是等间隔的。因此设计的系 统首先应对原始心电信号做线性插值处理，使其变为等间隔的数字信号，否则直接处理后会 出现偏差，根据心电信号的特点, 把时间分隔成0.001s。添加的幅值点采用一次线性插值。 对二维数据进行插值,相连幅值间数据的插值根据时间进行,运算公式如下: tt t Nt/0.001 AA A t t 0.001 A A A/N i i1， ， i i1 ， j j1 ， j j1 t A 其中 是第 i个数据时间点，A 是与之对应的数据，N 是两数据之间需要的插值数， 是i i 需要插值的两点数据差，i1，2，3，arraysize，j 1，2，3，N 1，t t ，A A ， j1 i1 j1 i1 i，j 1时数组t ,A 依次排列,即得到了插值后等间隔的新数据。 j j 将插值后的幅值数据另存为txt文件，以供后续程序调用。 4.1.3根据心电信号的频域特征，设计相应的低通和带通滤波器  一般正常人的心电信号频率在0.7～100HZ范围内，幅度为10 V (胎儿)～5mV(成人)。 人体心电信号微弱，信噪比小，因此，在采集心电信号时，易受到仪器、人体活动等因素的 影响，而且所采集的心电信号常伴有干扰。采集心电数据时，由于人的说话呼吸，常常会混 有约为0.1Hz到0.25Hz频段的干扰，对于这些低频干