《数字信号处理》课程设计项目2 心电信号分析与处理.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》 课程设计项目2 心电信号分析与处理8.1 心电信号中的工频干扰估计与消除8.2 心电信号中的QRS波检测8.3 孕妇心电图中胎儿心电信号的提取 8.1 心电信号中的工频干扰估计与消除一、课程设计的研究背景 数字信号处理技术在生物医学工程领域的应用十分广泛。用于诊断心脏病的主要技术之一是心电图（Electrocadiogram,简称ECG），心电图机通过测量在体表的电位来推测心脏的缺陷。利用数字信号处理技术可对心电信号自动进行滤波、分析和识别的工作，以辅助医生进行疾病诊断并减轻工作强度。人体心电信号常常受到电源引起的工频干扰（频率为50Hz或60Hz），其时域及频域波形如图8.1所示，对心电图诊断产生严重影响，必须对其频率进行估计从而将其消除。图8.1 心电信号中的工频干扰 二、课程设计目标要求 本课题研究的目标是：基于Matlab软件编程，采用FIR滤波器及LMS自适应滤波算法对心电信号中的工频干扰进行估计与消除，总结各种算法的优缺点，进行实验结果分析并撰写4000~5000字的课题论文。课题论文包含以下内容：1)题目、摘要、关键词、引言2)内容：心电信号中的工频干扰估计与消除3)包括：理论、程序(及注释)、图形、结果分析4)结论、参考文献三、课程设计内容与参考设计MATLAB程序模块，实现以下功能：(一) 心电信号的读取与显示 本课题的心电数据来自美国威斯康星大学UW DigiScope软件所产生的仿真心电信号，文件名为：ECG_X1.txt。设计要求: 读取该心电信号,并叠加工频干扰噪声,适当设置其幅度A与初相位；显示含有工频噪声的心电波形（设取样频率fs=150Hz,取样点数N=512）。参考程序：见\课程设计项目参考源码\kcsj2\课程设计题2.1\kcsj211.m (在空格中填入注释) 图8.2 心电信号的读取及显示(二)工频干扰的估计与消除1.FIR滤波器消除工频干扰(1)海宁滤波器 由于工频干扰在心电图中是高频噪声，故可通过平滑滤波器进行低通滤波消噪。海宁滤波器是一个简单的FIR型平滑滤波器，其系统函数为：(8.1) 设计要求: [1]编写程序(程序名：kcsj212.m),利用海宁滤波器对含有工频的心电信号ECGY.txt进行低通滤波消噪，列出程序清单，并画出海宁滤波器的幅频特性、零极点分布图及消噪后的心电波形;观察此波形，工频干扰是否都消除了？参考程序：见\课程设计项目参考源码\kcsj2\课程设计题2.1\kcsj212.m(在空格中填入注释)[2]对原有的程序进行修改，设计3级海宁滤波器进一步改善滤波效果，列出程序清单，并画出消噪后的心电波形。图8.3 海宁滤波器消除心电工频干扰 (2)窗函数法设计FIR带阻滤波器 由于海宁滤波器阶数较低，过渡带较宽，虽然消除了工频噪声，但也抑制了过渡带中的有用频率成分，从而导致了心电信号的衰减。可采用窗函数法设计一个FIR带阻滤波器，在消除工频干扰的同时尽量保留心电信号中的有用频率成分。设计要求：编写程序（程序名：kcsj213.m），利用窗函数法设计一个FIR带阻滤波器，设置适当的滤波参数，消除工频干扰，列出程序清单，并画出该滤波器的单位脉冲响应、幅频特性及消噪后的心电波形及其幅频谱。参考程序：见\课程设计项目参考源码\kcsj2\课程设计题2.1\kcsj213.m(在空格中填入注释)图8.4 FIR带阻滤波器消除心电工频干扰 2.LMS自适应滤波算法消除工频干扰 上述FIR滤波器当工频干扰的中心频率有漂移时，滤波效果将受到一定的影响，可设计一个自适应噪声抵消系统（图8.5）来消除，自适应滤波器于1967年由威德诺（B.Widrow）等人提出，如图8.5所示，输入信号为有用S(n)和噪声信号N(n)之和，参考输入信号可以是仅和N(n)相关的噪声信号R(n), 通过自适应滤波器调整其输出Y(n),在最小均方误差准则下得到N(n)的一个最佳估计，将此估计值和主输入通道相减从而较好地消除主输入通道噪声。输出误差e(n)为： e(n)=S(n)+N(n)-Y(n)均方误差的估计为： E[e2(n)]=E[S2(n)+E{[N(n)-Y(n)]2}+2E{S(n)[N(n)-Y(n)]}由于S(n)与R(n)及N(n)不相关，故期望值最小值为： E[e2(n)]= E[S2(n)]+ E{[N(n)-Y(n)]2}通过调整自适应滤波器，当Y(n)→N(n)时输出的均方误差最小，此时系统输出十分接近于输入信号中的有用成分S(n)。这种通过均方误差反