基于MATLAB的心电信号的分析与处理设计医学信号处理.docx

一.设计内容 根据给定的一段 MIT-BIH 心电信号（122 号），画出心电信号的时域波形和频 谱图（幅频和相频）； 根据心电频率范围及其噪声的频率范围设计 2 个滤波器（一个 IIR,一个 FIR） 实现对心电信号滤波。滤波器的种类（高通，低通，带通，带阻），滤波器性 能指标（通阻带截止频率，衰减系数），滤波器的设计方法（有 IIR 冲击响应 不变法和双线性变换法， FIR 有窗函数法及频率抽样法）等自行设计。要求 输出所设计的滤波器的系统函数，画出滤波器的频率响应（幅频响应和相频 响应）曲线。 用该滤波器对心电信号进行滤波，画出滤波以后心电信号的时域波形和频谱 （幅频）；分析信号滤波前后心电信号的时域和频域的变化； 两个滤波器滤波效果分析； 目录 一.设计内容................................................................................................1 二．设计流程.............................................................................................1 2.1 设计步骤…………………………………………………………………………..1 2.2 设计流程图………………………………………………………………………..1 三．设计过程及分析…………………………………………………....3 3.1 心电信号的特点………………………………………………….………………3 3.2 原心电信号导入及分析……………………………………….…………………4 3.3.滤波器设计…………………………………………………….…………………5 3.4.设计低通滤波器……………………………………………….…………………5 3.5.高通滤波器设计………………………………………….………………………7 3.6.滤波前后心电信号对比分析………………………………………………….…7 四．课程设计总结……………………………………………………..10 附录………………………………………………………………….….11 二．设计流程 2.1 设计步骤如下： 将给定的 MIT-BIH 心电信号（122 号）导入 matlab，并画出心电信号的 时域波形和频谱图（幅频和相频）； 用双线性 Z 变换法设计一低通数字滤波器； 用 IIR 冲击响应不变法设计一高通数字滤波器； 对原心电信号进行低通滤波，并画出滤波后的心电信号的时域波形和频 谱图（幅频和相频）； 对原心电信号进行高通滤波，并画出滤波后的心电信号的时域波形和频 谱图（幅频和相频）； 2.2 流程图如图（a）所示： 设计低通滤波器  设计高通滤波器 导入原心电信号  行 低 通 滤 波  对 原 心 电 信 号 进  行 高 通 滤 波  对 原 心 电 信 号 进  通 滤 波 后 波 形 图  滤 波 连 续 通 过 低 ， 高 低 通 滤 波 后 波 形 图  高通滤波后波形 图 图（a） 三．设计过程及分析 3.1 心电信号的特点 ECG 信号是心脏活动过程中心脏的肌肉和神经电活动的综合，它的信源是 心脏。由于生命体自身的机理以及信号源的不可触及性， ECG 信号象其他生命 体系统产生的信号一样，具有以下几个主要特点： 是随机性较强：即信号无法用确定的函数式来描述； 是噪声背景强：预测的有用信号往往淹没在许多无用信号中； 是信号频率低：心电信号基本是周期的信号，同样，它又有着非平稳的性 质，ECG 信号随时都处于微小的变动之中，而不是固定不变的。ECG 节律会随 运动而变，甚至会随某种疾病(如糖尿病)而变。虽然不同的人之间的心电波形差 异很大，但所有的正常的心电波形的都可以划分为 P 波、QRS 波、T 波和 U 波 几个部分。心电信号是心电电压振幅随时间变化的信号，即是时域信号。采用体 表电极作为传感器时，信号的幅度范围是 10uV-～4mV，典型值是 lmV。一个典 型的导联电压对时间的函数关系如图（b）所示： 图（b） 1．P 波：P 波反映左右两心房的电激动过程，也叫心房除极波。其起点表示窦 房结开始激动，终点表示两心房激动结束。正常 P 波具有以下特点： 方向：I、II、aVF、V4^ Ⅳ6 各导联应直立向上，aVR 导联波形应倒置； 时程：P 波时程即波的宽度应小于 O．11 秒； 电压：电压表示了 P 波的振幅，对肢体导联应小于 0．25mV，对心前导联应小 于 0．20mV。 2．PR 间期：PR 问期是 P 波起点至 QRS 波起点的时间间隔，它表示激动从窦房 结发出，经结间传导束房室