基于matlab的心电信号分析与处理小论文

基于MATLAB的心电信号分析与处理摘要本课题设计了一个简单的心电信号分析系统。直接采用MATLAB语言编程对输入的原始心电信号进行处理，并通过MATLAB语言编程设计对其进行时域和频域的波形频谱分析，根据具体设计要求完成系统的程序编写、调试及功能测试,得出一定的结论。THISTOPICHASDESIGNEDASIMPLEECGANALYSISSYSTEMDIRECTUSEOFMATLABPROGRAMMINGLANGUAGEORIGINALECGSIGNALINPUTISPROCESSED,ANDITSWAVEFORMSPECTRUMANALYSISOFTHETIMEDOMAINANDFREQUENCYDOMAINMATLABLANGUAGEPROGRAMMINGTHROUGHDESIGN,PREPAREDINACCORDANCEWITHSPECIFICDESIGNREQUIREMENTSTOCOMPLETETHESYSTEMOFPROCEDURES,DEBUGGINGANDFUNCTIONALTESTING,TOOACERTAINCONCLUSION关键字MATLAB、心电信号、滤波一、课题目的及意义心电信号是人类最早研究并应用于医学临床的生物信号之一，它比其它生物电信号更易于检测，并且具有较直观的规律性，因而心电图分析技术促进了医学的发展。然而，心电图自动诊断还未广泛应用于临床，从国内外的心电图机检测分析来看，自动分析精度还达不到可以替代医生的水平，仅可以为临床医生提供辅助信息。其主要原因是心电波形的识别不准，并且心电图诊断标准不统一。因此，探索新的方法以提高波形识别的准确率，寻找适合计算机实现又具诊断价值的诊断标准，是改进心电图自动诊断效果，扩大其应用范围的根本途径。如何把心电信号的特征更加精确的提取出来进行自动分析，判断出其异常的类型成了亟待解决的焦点问题。本课题通过MATLAB语言编程，对原始心电信号进行一定的分析处理。ECGISTHEFIRSTHUMANSTUDYANDONEBIOLOGICALSIGNALAPPLIEDTOCLINICALMEDICINE,ITISEASIERTODETECTTHANOTHERBIOLOGICALSIGNALS,ANDHASAMOREINTUITIVEREGULARITY,THUSECGANALYSISTECHNOLOGYFORTHEDEVELOPMENTOFMEDICALSCIENCEHOWEVER,ECGAUTOMATICDIAGNOSISHASNOTBEENWIDELYUSEDINCLINICAL,ECGMACHINEDETECTIONANALYSISFROMHOMEANDABROAD,THEACCURACYOFTHEAUTOMATICANALYSISCANREPLACETHEDOCTORHASNOTYETREACHEDTHELEVELOFAIDCANONLYPROVIDEINFORMATIONTOTHECLINICIANTHEMAINREASONISNOTALLOWEDTOIDENTIFYTHEECGWAVEFORMANDECGDIAGNOSTICCRITERIAARENOTUNIFORMTHEREFORE,TOEXPLORENEWWAYSTOIMPROVETHEACCURACYOFWAVEFORMRECOGNITION,SEARCHINGFORCOMPUTERIMPLEMENTEDBUTALSOWITHTHEDIAGNOSTICVALUEOFTHEDIAGNOSTICCRITERIA,THEEFFECTOFIMPROVINGTHEELECTROCARDIOGRAMAUTOMATICDIAGNOSIS,THEFUNDAMENTALWAYTOEXPANDITSRANGEOFAPPLICATIONSHOWTOECGEXTRACTMOREACCURATEFEATUREAUTOMATICALLYANALYZE,DETERMINETHETYPEOFEXCEPTIONHASBECOMETHEFOCALPOINTOFTHEIRPROBLEMSTOBESOLVED二、设计思路1根据给定的一段MITBIH心电信号，画出心电信号的时域波形和频谱图（幅频和相频）；ACCORDINGTOAGIVENPERIODOFMITBIHECG,ECGDRAWATIMEDOMAINWAVEFORMANDSPECTRUMAMPLITUDEANDPHASEFREQUENCY2根据心电信号频率范围及其噪声的频率范围设计2个滤波器（一个IIR，一个FIR）实现对心电信号滤波。滤波器的种类（高通，低通，带通，带阻），滤波器性能指标（通阻带截止频率，衰减系数），滤波器的设计方法（IIR有冲击响应不变法和双线性变换法，FIR有窗函数法及频率抽样法）等自行设计。要求输出所设计的滤波器的系统函数，画出滤波器的频率响应（幅频响应和相频响应）曲线；DESIGNTWOFILTERSAIIR,AFIRFILTERTOACHIEVEECGECGACCORDINGTOTHEFREQUENCYRANGEANDNOISEFREQUENCYRANGEFILTERTYPEHIGHPASS,LOWPASS,BANDPASS,BANDSTOP,FILTERPERFORMANCEINDICATORSPASSSTOPBANDCUTOFFFREQUENCY,ATTENUATIONCOEFFICIENT,FILTERDESIGNMETHODSIIRHASIMPULSERESPONSEINVARIANCEANDBILINEARTRANSFORM,FIRWITHWINDOWFUNCTIONMETHODANDFREQUENCYSAMPLINGMETHODOWNDESIGNSYSTEMFUNCTIONSREQUIREDOUTPUTFILTERDESIGN,DRAWTHEFILTERFREQUENCYRESPONSEAMPLITUDEFREQUENCYRESPONSEANDPHASEFREQUENCYRESPONSECURVE3用该滤波器对心电信号进行滤波，画出滤波以后心电信号的时域波形和频谱（幅频）；分析信号滤波前后心电信号的时域和频域的变化；WITHTHISFILTERECGSIGNALFILTERING,FILTERINGAFTERTHEDRAWECGWAVEFORMANDSPECTRUMAMPLITUDEANDFREQUENCYCHANGESBEFOREANDAFTERTHETIMEDOMAINANALYSISOFECGSIGNALFILTERINGANDFREQUENCYDOMAIN4两个滤波器滤波效果异同分析SIMILARITIESANDDIFFERENCESBETWEENTHETWOFILTERSFILTERINGEFFECTANALYSIS三、设计过程及分析31原始心电信号分析用LOAD函数将原心电信号导入BLOADCUSERSTOSHIBADESKTOP信号处理102TXT，并画出心电信号的时域波形和频谱图（幅频和相频），ECGWITHALOADFUNCTIONWILLIMPORTBLOADCUSERSTOSHIBADESKTOPSIGNALPROCESSING102TXT,ANDDRAWTHEECGWAVEFORMANDSPECTRUMAMPLITUDEANDPHASEFREQUENCY如图1所示图1原始心电信号的时域波形图及频谱图心电信号由于受到人体诸多因素的影响，因而有着一般信号所没有的特点1信号弱。心电信号是体表的电生理信号，一般比较微弱，幅度在10PV5MV，频率为005100HZ。例如从母体腹部收取到的胎儿心电信号仅10ZV50IV。2噪声强。由于人体自身信号弱，加之人体又是一个复杂的系统，因此信号容易受到噪声干扰。3随机性强。心电信号不仅是随机的，而且是非平稳的。同时，在心电图检测过程中极易受到各种噪声源的干扰，从而使图像质量变差，使均匀和连续变化的心电数值产生突变，在心电图上形成一些毛刺。使原本很微弱的信号很难和噪声进行分解。可能出现的噪声有如下的种类1工频干扰工频干扰是由电力系统和人体的分布电容引起的，其频率包括50HZMITBIH数据库数据工频因为是美国标准，所以是60HZ的基波及其各次谐波，其幅值成分在ECG峰一峰值的050范围内变化。2引起基线漂移的干扰心电信号有时候会出现信号基线起伏不平的现象，造成这样的现象有很多原因，主要的有呼吸运动人体呼吸时胸腔内器官和组织会发生一定程度的变化，会对在体表记录到的心电图波形的幅度和形态有所影响，表现为基线随呼吸产生周期性或非周期性漂移，从而导致心电波形的幅度随呼气和吸气而分别上抬和下移。呼吸运动是引起心电基线漂移的主要原因。运动伪迹运动伪迹是由于人体轻微运动造成电极与入体的接触电阻发生变化而引入的一种干扰，它的产生原因仅仅是接触电阻的变化，而不是接触的断续。这种干扰同样导致信号基线的变化，但不是基线的跃变。信号记录和处理中电子设备引起的干扰这种干扰对信号影响很大，严重时可完全淹没心电信号或使得基线剧烈漂移，其中导联开路和放大器的热移是主要因素。这种干扰往往无法通过心电分析算法来校正。由于心电波形已经完全畸变，此时对这些数据分析已无太大意义。所以一般跳过此段数据。3高频噪声心电信号中的高频噪声主要是肌电噪声。肌肉收缩会产生MV级的肌电干扰，表现为心电图上不规则的细小波纹，使心电图模糊不清或产生失真。肌电噪声的特点是频率范围较广，频谱分布非常复杂。32设计滤波器（1）IIR滤波器的设计1）IIR滤波器的设计过程按照技术要求设计一个模拟滤波器，得到模拟低通滤波器的传输函数H（S），再按一定的转换关系将HS转换成数字低通滤波器的系数函数HZ。这样设计的关键问题就是找到这样的转换关系，将S平面上的HS转换成Z平面上的HZ。2）巴特沃斯滤波器分母多项式的因式表示，如表1所示表1巴特沃斯滤波器分母多项式的因式表示3）巴特沃斯低通滤波器的阶数公式NLOG1010AS/101/10RP/101/2LOG10WS/WP）4）巴特沃斯低通滤波器函数由巴特沃斯低通滤波器的阶数公式和巴特沃斯滤波器分母多项式的因式表示求出归一化巴特沃斯低通滤波器HASSN7则HASS1/S1S204450S1S21247S1S218022S15）巴特沃斯低通滤波器的频域特性，如图2所示图2巴特沃斯低通滤波器的相频和幅频特性6）经过巴特沃斯低通滤波器器后心电信号的时域波形和频谱图，AFTERBUTTERWORTHLOWPASSFILTERISWHENTHEECGWAVEFORMANDSPECTRUMDOMAIN如图3所示图3经过巴特沃斯低通滤波器器后心电信号的时域波形图和频谱图对比原始信号的时域波形图和频谱图可得通过低通滤波器后的心电信号波形图可以明显看出波形变得平滑，由工频干扰产生的毛刺被低通滤波器成功滤除。COMPARISONOFTHEORIGINALSIGNALANDTHETIMEDOMAINWAVEFORMSPECTRUMCANBEOBTAINEDAPPARENTWAVEFORMBECOMESSMOOTH,THEBURRFREQUENCYINTERFERENCEGENERATEDBYTHESUCCESSOFTHELOWPASSFILTERTOFILTEROUTECGWAVEFORMAFTERTHELOWPASSFILTER（2）FIR滤波器的设计1）FIR滤波器的设计过程给定理想的频率响应函数HDEJW及技术指标,W求出理想的单位抽样响应HDN根据阻带衰减选择窗函数WN根据过渡带宽度确定N值NA/W求所设计的FIR滤波器的单位脉冲响应HNHDNWN计算频率响应HDEJW,验算指标是否满足要求。2）布拉克曼窗低通滤波器的频域特性，如图4所示图4布拉克曼窗低通滤波器的相频和幅频特性3）经过布拉克曼窗低通滤波器器后心电信号的时域波形和频谱图，AFTERBLACKMANWINDOWLOWPASSFILTERISECGWAVEFORMANDSPECTRUM,如图5所示图5经过布拉克曼窗低通滤波器器后心电信号的时域波形和频谱图对比原始信号的时域波形图和频谱图可得通过低通滤波器后的心电信号波形图可以明显看出波形变得平滑，由工频干扰产生的毛刺被低通滤波器成功滤除。4、结果分析通过用MATLAB设计的IIR低通滤波器和FIR低通滤波器对带有噪声的心电信号进行处理和分析，消除了原心电信号中由工频干扰产生的毛刺现象，虽然仍存在一些误差，但基本上达到了预期效果。通过做课程设计熟悉并掌握了相关的MATLAB操作环境；进一步加深了对滤波器设计相关知识的理解，熟练了滤波器设计的方法和过程。BYIIRANDFIRLOWPASSFILTERWITHALOWPASSFILTERDESIGNMATLABECGWITHNOISEPROCESSINGANDANALYSIS,ELIMINATINGTHEORIGINALECGFREQUENCYINTERFERENCEGENERATEDBYTHEWORKERSBURRPHENOMENON,ALTHOUGHTHEREARESTILLSOMEERRORSBUTBASICALLYACHIEVEDTHEEXPECTEDRESULTSBYDOINGCURRICULUMDESIGNANDFAMILIARWITHTHERELEVANTOPERATINGENVIRONMENTMATLABFURTHERDEEPENEDTHEUNDERSTANDINGOFFILTERDESIGNKNOWLEDGE,SKILLEDFILTERDESIGNMETHODSANDPROCESSES巴特沃斯数字低通滤波器巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。在振幅的对数对角频率的波得图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。窗函数法窗函数法是设计FIR滤波器的最主要方法之一，实际中遇到的离散时间信号总是有限长的，因此不可避免的要遇到数据截短的问题，在信号处理中，对离散序列的截短是通过序列与窗函数相乘来实现的。在信号处理中，窗函数是一种除在给定区间之外取值均为0的实函数。譬如在给定区间内为常数而在区间外为0的窗函数被形象地称为矩形窗。任何函数与窗函数之积仍为窗函数，所以相乘的结果就像透过窗口“看”其他函数一样。窗函数在光谱分析、滤波器设计以及音频数据压缩等方面有广泛的应用。程序流程图CLEARALLCLOSEALLCLC时域波形图以及频谱图BLOADCUSERSTOSHIBADESKTOP信号处理102TXTTB,1时间CB,2幅值FIGURE1SUBPLOT311PLOTT,CTITLE原始心电信号的时域波形图XLABEL时间T/SYLABEL幅值/AN3600MFFTC,NFS360采样频率FFS/N0N1定位坐标SUBPLOT312PLOTF,ABSMTITLE原始心电信号的频谱图XLABEL频率/HZYLABEL幅值/DBAXIS0,360,0,150SUBPLOT313PLOTF,ANGLEMTITLE原始心电信号的相频图XLABEL频率/HZYLABEL相角/RADAXIS0,360,5,5低通IIR滤波器WP2PI10/FS通带数字频率转换成模拟频率WS2PI15/FS通带截至数字频率转换成模拟频率RP2通带最大衰减RS20阻带最小衰减N,WCBUTTORDWP,WS,RP,RS,S确定最小阶数N和频率参数WCBZ,AZBUTTERN,WC得巴特沃斯归一化低通原型H,WFREQZBZ,AZ生成频率响应参数F1W/PIFS/2采样频率转换成模拟采样频率Y1FILTERBZ,AZ,C使