基于MATLAB的心电信号的分析与处理设计Word下载.docx
《基于MATLAB的心电信号的分析与处理设计Word下载.docx》由会员分享，可在线阅读，更多相关《基于MATLAB的心电信号的分析与处理设计Word下载.docx（13页珍藏版）》请在冰豆网上搜索。
3．掌握序列离散傅里叶变换的MATLAB实现，并进行频谱分析
4．熟练掌握使用MATLAB设计IIR或FIR数字滤波器
5．学会用MATLAB对信号进行分析和处理
二、设计内容、技术条件和要求
一）设计内容与技术条件
1．根据给定的一段MIT-BIH心电信号（101号），画出心电信号的时域波形和频谱图（幅频和相频）；
2．根据心电信号频率范围及其噪声的频率范围设计2个滤波器（一个IIR，一个FIR）实现对心电信号滤波。
滤波器的种类（高通，低通，带通，带阻），滤波器性能指标（通阻带截止频率，衰减系数），滤波器的设计方法（IIR有冲击响应不变法和双线性变换法，FIR有窗函数法及频率抽样法）等自行设计。
要求输出所设计的滤波器的系统函数，画出滤波器的频率响应（幅频响应和相频响应）曲线；
3．用该滤波器对心电信号进行滤波，画出滤波以后心电信号的时域波形和频谱（幅频）；
分析信号滤波前后心电信号的时域和频域的变化；
4．两个滤波器滤波效果异同分析；
5．运用GUI设计一个心电信号处理系统界面。
（选作）
二）设计要求
1．根据滤波器的性能指标要求，设计数字滤波器；
2．程序中按照IIR滤波器的步骤一步步完成设计；
尽可能的少调用MATLAB自带
的函数文件；
3．设计程序要通用性好，整齐易懂，并要求主要语句有注释；
4．设计结果中的图示要美观，整齐，有标题，有纵横坐标标示；
5．课程设计报告要有理论依据、设计过程，结果分析。
报告要求实事求是、文理
通顺、字迹端正。
三、时间进度安排
9月1日：
查阅资料，复习医学信号处理课程的相关内容，熟悉MATLAB的工作环境，
掌握运算编程方法。
9月2日~5日：
按照课程设计要求，完成课程设计的程序设计部分。
9月8日~10日：
调试并完善程序，撰写课程设计说明书。
9月11~12日：
答辩。
四、主要参考文献
参考资料
（1）数字信号处理——理论、算法与实现胡广书清华大学出版社
（2）数字信号处理及MATLAB实现余成波清华大学出版社
（3）MATLAB7.0在数字信号处理中的应用罗军辉机械工业出版社
（4）数字信号处理原理及其MATLAB实现从玉良电子工业出版社
（5）MATLAB帮助文件
相关网站
（1）MATLAB中国论坛
（2）MATLAB学习网
（3）MATLAB中文论坛
指导教师签字：
侯海燕2014年8月25日
1、设计目的意义·
·
1
2、设计内容·
2.1设计原理·
2.2要求·
2.3内容·
3、计过程及分析·
3.1原始心电信号分析·
3.2设计滤波器·
3
4、结果分析·
7
5、参考文献·
附录·
8
1、设计目的意义
本次课程主要达到以下几个目的：
（1）熟练掌握使用MATLAB程序设计方法
（2）掌握序列离散傅里叶变换的MATLAB实现，并进行频谱分析
（3）熟练掌握使用MATLAB设计IIR或FIR数字滤波器
2、设计内容
2.1设计原理
（1）频谱分析（傅里叶变换）
（2）滤波器设计
IIR（冲激响应不变法，双线性变换法）详解见附录
FIR（窗函数法，频率抽样法）详解见附录
（3）信号滤波（时域卷积，频域相乘）
2.2要求
要求设计出心电数据处理的处理与分析程序。
（1）处理对象：
心电数据；
（2）内容：
心电数据仿真，心电数据处理（仿真数据，真实数据）；
（3）结果：
得到处理结果。
2.3内容
（1）原始心电信号的时域波形和频谱图（幅频和相频）。
（2）IIR滤波器的频率响应（幅频响应和相频响应）曲线；
原始心电信号经IIR滤波器滤波后的心电信号的时域波形和频谱图（幅频和相频）；
（3）FIR滤波器的频率响应（幅频响应和相频响应）曲线；
原始心电信号经FIR滤波器滤波后的心电信号的时域波形和频谱图（幅频和相频）；
3、设计过程及分析
3.1原始心电信号分析
用load函数将原心电信号导入b=load（'
C:
\Users\Administrator\Desktop\课程设计\101ecg.txt'
），并画出心电信号的时域波形和频谱图（幅频和相频），如图1所示:
图1原始心电信号的时域波形图及频谱图
心电信号由于受到人体诸多因素的影响，因而有着一般信号所没有的特点：
（1）信号弱。
心电信号是体表的电生理信号，一般比较微弱，幅度在10pV～5mV，频率为0．05～100Hz。
例如从母体腹部收取到的胎儿心电信号仅10／zV～50／IV。
（2）噪声强。
由于人体自身信号弱，加之人体又是一个复杂的系统，因此信号容易受到噪声干扰。
（3）随机性强。
心电信号不仅是随机的，而且是非平稳的。
同时，在心电图检测过程中极易受到各种噪声源的干扰，从而使图像质量变差，使均匀和连续变化的心电数值产生突变，在心电图上形成一些毛刺。
使原本很微弱的信号很难和噪声进行分解。
可能出现的噪声有如下的种类：
1）工频干扰
工频干扰是由电力系统和人体的分布电容引起的，其频率包括50Hz（MIT-BIH数据库数据工频因为是美国标准，所以是60Hz）的基波及其各次谐波，其幅值成分在ECG峰一峰值的0—50％范围内变化。
2）引起基线漂移的干扰
心电信号有时候会出现信号基线起伏不平的现象，造成这样的现象有很多原因，主要的有：
①呼吸运动人体呼吸时胸腔内器官和组织会发生一定程度的变化，会对在体表记录到的心电图波形的幅度和形态有所影响，表现为基线随呼吸产生周期性或非周期性漂移，从而导致心电波形的幅度随呼气和吸气而分别上抬和下移。
呼吸运动是引起心电基线漂移的主要原因。
②运动伪迹运动伪迹是由于人体轻微运动造成电极与入体的接触电阻发生变化而引入的一种干扰，它的产生原因仅仅是接触电阻的变化，而不是接触的断续。
这种干扰同样导致信号基线的变化，但不是基线的跃变。
③信号记录和处理中电子设备引起的干扰这种干扰对信号影响很大，严重时可完全淹没心电信号或使得基线剧烈漂移，其中导联开路和放大器的热移是主要因素。
这种干扰往往无法通过心电分析算法来校正。
由于心电波形已经完全畸变，此时对这些数据分析已无太大意义。
所以一般跳过此段数据。
3）高频噪声
心电信号中的高频噪声主要是肌电噪声。
肌肉收缩会产生mV级的肌电干扰，表现为心电图上不规则的细小波纹，使心电图模糊不清或产生失真。
肌电噪声的特点是频率范围较广，频谱分布非常复杂。
3.2设计滤波器
（1）IIR滤波器的设计
1）IIR滤波器的设计过程：
按照技术要求设计一个模拟滤波器，得到模拟低通滤波器的传输函数H（s），再按一定的转换关系将H（s）转换成数字低通滤波器的系数函数H（z）。
这样设计的关键问题就是找到这样的转换关系，将s平面上的H（s）转换成z平面上的H（z）。
2）巴特沃斯滤波器分母多项式的因式表示，如表1所示：
表1巴特沃斯滤波器分母多项式的因式表示
3）巴特沃斯低通滤波器的阶数公式
N=log10（（10^（As/10）-1）/（10^（Rp/10）-1））/（2*log10（ws/wp）））
4）巴特沃斯低通滤波器函数
由巴特沃斯低通滤波器的阶数公式和巴特沃斯滤波器分母多项式的因式表示求出归一化巴特沃斯低通滤波器Has（s）
N=7
则Has（s）=1/（（s+1）*（s^2+0.4450s+1）*（s^2+1.247s+1）*（s^2+1.8022s+1））
5）巴特沃斯低通滤波器的频域特性，，如图2所示：
图2巴特沃斯低通滤波器的相频和幅频特性
6）经过巴特沃斯低通滤波器器后心电信号的时域波形和频谱图，如图3所示：
图3经过巴特沃斯低通滤波器器后心电信号的时域波形图和频谱图
对比原始信号的时域波形图和频谱图可得通过低通滤波器后的心电信号波形图可以明显看出波形变得平滑，由工频干扰产生的毛刺被低通滤波器成功滤除。
（2）FIR滤波器的设计
1）FIR滤波器的设计过程：
1给定理想的频率响应函数Hd（e^jw）及技术指标δ,Δw;
2求出理想的单位抽样响应hd（n）;
3根据阻带衰减选择窗函数w（n）;
4根据过渡带宽度确定N值N=A/Δw;
5求所设计的FIR滤波器的单位脉冲响应h（n）=hd（n）*w（n）;
6计算频率响应Hd（e^jw）,验算指标是否满足要求。
2）布拉克曼窗低通滤波器的频域特性，如图4所示：
图4布拉克曼窗低通滤波器的相频和幅频特性
3）经过布拉克曼窗低通滤波器器后心电信号的时域波形和频谱图，如图5所示
图5经过布拉克曼窗低通滤波器器后心电信号的时域波形和频谱图
4、结果分析
通过用matlab设计的IIR低通滤波器和FIR低通滤波器对带有噪声的心电信号进行处理和分析，消除了原心电信号中由工频干扰产生的毛刺现象，虽然仍存在一些误差，但基本上达到了预期效果。
通过做课程设计熟悉并掌握了相关的matlab操作环境；
进一步加深了对滤波器设计相关知识的理解，熟练了滤波器设计的方法和过程。
5、参考文献：
[1]陈天华.数字图像处理[M].北京.清华大学出版社，2009
[2]刘卫国.MALTAB程序设计与应用[M].北京.高等教育出版社，2008
[3]程正兴.小波分析算法与应用.西安.西安交通大学出版社，1998
[4]程佩青.数字信号处理.北京.清华大学出版社，2007
[5]苏金明，王永利.MALTAB应用指南[M].上册.北京电子工业出版社，2004
[6]夏良正.数字图像处理（修订版）[M].南京.东南大学出版社，1999
[7]霍红涛，林小竹，何薇.数字图像处理[M].北京.北京理工大学出版社，2003
[8]张开滋，郭继
升级会员 