基于gui心电监护仪的设计与实现

南阳理工学院本科毕业设计（论文）基于GUI心电监护仪的设计与实现DESIGNANDIMPLEMENTATIONOFGUIBASEDECG学院计算机与信息工程学院专业通信工程学生姓名学号指导教师（职称）评阅教师完成日期基于GUI心电监护仪的设计与实现摘要目前心脏病的死亡率远远高于其他疾病，已成为威胁人类生命安全的“第一杀手”。由于一些患者在发病时感到不适，而到医院检查时症状消失，进行心电图检查时得不到明显异常的心电信号，最终延误了诊断和治疗。针对这些问题,开发一套携带方便、低成本的远程心电监护系统具有深远的理论研究意义和实际应用价值。通过将采集的数据存储到TXT文本文件中，在利用MATLAB软件读取数据、滤波、并进行一系列的计算后，可以得出检测的结果，并将结果以文件形式保存下来，以便随时查阅。设计充分利用MATLAB的强大分析仿真能力，使系统的完整性和准确性得以很好的实现。关键词MATLAB；GUI；心电信号；监护仪DESIGNANDIMPLEMENTATIONOFGUIBASEDECGABSTRACTCURRENTLYHEARTDISEASEMORTALITYRATEISMUCHHIGHERTHANOTHERDISEASES,HASBECOMEATHREATTOTHESAFETYOFHUMANLIFE,“THEFIRSTKILLER“SINCESOMEPATIENTSFEELUNWELLINTHEDISEASE,BUTNOTSIGNIFICANTLYABNORMALECGTOTHEHOSPITALWHENTHESYMPTOMSDISAPPEAR,ECGEXAMINATION,EVENTUALLYDELAYEDDIAGNOSISANDTREATMENTTOSOLVETHESEPROBLEMS,FACILITATETHEDEVELOPMENTOFAPORTABLE,LOWCOSTREMOTEECGMONITORINGSYSTEMHASFARREACHINGTHEORETICALSIGNIFICANCEANDPRACTICALVALUEBYSTORINGTHEDATAWILLBECOLLECTEDTXTTEXTFILE,USINGMATLABSOFTWARETOREADTHEDATA,FILTERING,ANDASERIESOFCALCULATIONS,THETESTRESULTSCANBEOBTAINED,ANDTHERESULTSSTOREDINDOCUMENTFORMDOWNTOINSPECTIONATANYTIMEDESIGNEDTOTAKEFULLADVANTAGEOFMATLABSPOWERFULANALYSISANDSIMULATIONCAPABILITIES,THECOMPLETENESSANDACCURACYOFTHESYSTEMTOACHIEVEVERYGOODKEYWORDSMATLABGRAPHICALUSERINTERFACEELECTROCARDIOGRAMMONITOR目录1引言111研究背景和意义112研究现状22心电监护仪的医学基础321人体心电信号的产生机理322心电图及心电信号的特征分析4221心脏电传导过程分析4222心电信号时域特征分析4223心电信号频域特征分析5224心电信号的特点63基于MATLAB的心电信号处理631MATLAB处理ECG的优势632心电信号噪声的来源及特点733MATLAB对ECG信号的处理过程7331MITBIH简介7332MATLAB读取心电信号8333心电信号的滤波8334心电信号频谱及功率谱的计算11335心率的计算12336信号的保存和回放124基于MATLAB的GUI界面设计1241GUI所要实现的功能1342GUI界面的布局设计1343GUI对应函数的编写14431信号的读取模块14432信号的滤波模块15433信号的频谱显示17434信号的心率计算模块18435信号保存模块18436信号回放模块1944GUI界面运行结果19总结21参考文献22致谢231引言目前市面上常见的便携式心电仪多数是采用了前后端的实现方式，前端是以单片机为核心的心电信号采集系统，后端多数采用的是处理性能较高的嵌入式微处理器。这种处理器性能强大，它使得心电仪在心电数据采集、处理、存储和显示等功能的基础上，还能够实现对心电数据的分析。1然而，这种心电仪在实现多种功能的同时，也有一些缺点结构比较复杂、功耗较大、成本也较高。另一方面，在导联个数上，在相当长的一段时间内，心电导联系统一般仅仅具有单导或三导联同步记录功能，市场上现在也还有很多这种产品。该类产品因为只支持少数的导联，因而它的液晶屏幕比较小，用户观察起来很不方便，只能通过自带的打印机将心电图打印出来之后才能较好的观察分析。另外，这种产品往往不适合复杂心脏疾病的诊断。目前很多厂商也在竞相开发支持多导联的心电仪产品。深圳迈瑞电子就是其中之一，它在便携式监护仪领域做出了带头作用，典型的产品如PM9000EXPRESS、PMS000等等1。11研究背景和意义心脏是人体血液循环的动力泵，心脏搏动是生命存在的重要标志，心脏搏动的节律也是人体生理状态的重要标志之一。心脏的基本活动包括电活动和机械活动，每个心动周期都是电活动在前，机械活动在后。心电信号是心脏电活动的一种客观表示方式，是一种典型的生物电信号，具有频率、振幅、相位、时间差等特征要素，比其他生物电信号更易于检测1，并具有一定的规律性。由于心电信号从不同方面和层次上反映了心脏的工作状态，因此在心脏疾病的临床诊断和治疗过程中具有非常重要的参考价值。对心电信号的采集和分析一直是生物医学工程领域研究的一个热点，是一项复杂的工程，涉及到降低噪声和抗干扰技术，信号分析和处理技术等不同领域，也依赖于生命科学和临床医学的研究进展。自1903年心电图引入医学临床以来，无论是在生物医学方面，还是在工程学方面，心电信号的记录、处理与诊断技术均得到了飞速的发展，并积累了相当丰富的资料。当前，心电信号的检测、处理仍然是生物医学工程界的重要研究对象之一。伴随着人们生活水平的提高，营养过剩和运动减少，生活和工作节奏的加快，社会老龄化的加剧，心脏病等心血管疾病的发病率明显上升。目前心脏病的死亡率远远高于其他疾病，是威胁人类生命安全的主要疾病，心脏病己成为威胁人类生命安全的“第一杀手“。由于心脏病具有病情隐蔽、发展缓慢、发病危险性高等特点，因此它对心脏病患者、特别是中老年心脏病患者的危害性极大，因此心脏系统疾病的防治和诊断己成为当今医学界面临的首要问题。另外，一些特殊的心脏病患者，在正常工作生活时发病，而到医院检查时症状消失，导致在医院无法检测到异常心电图，无法对病情做出诊断，耽误了治疗的最佳时机，所以心电监护是有重要意义的。随着社会老龄化的加剧，解决长期慢性病的监护目前已经是重要的社会问题。怎样才能使病人在家庭中得到更好医疗保健，同时又减少病人家属及社会的负担，是现在摆在有关研究人员及医生面前的一个重要课题。心电信号是人类最早研究并应用于临床医学的生物电信号之一，与其它生物电信号相比心电信号更易于检测并具有较直观的规律性。在医学发达的今天，心电信号相比心电信号更易于检测并具有较直观的规律性。在医学发达的今天，心电信号的监测与处理仍是生物医学领域重要的研究对象之一，是医学上诊断心脏系统疾病十分重要的检测依据。因此及时通过心电信号来发现并预防心脏病的发生是减少心脏病人发病和死亡率的一个有效途径。心电学自上世纪初到现在已经发展了一百多年。在这期间心电图检测技术本身不断发展完善，各种心电检测方法不断问世，到目前可以说心电检测已趋于成熟。以往的检测手段基本上都需要到医院在庞大的心电图仪器上进行检测，而心脏病人不可能每次发病时都具备到医院检测的条件，还有一点就是一些患者在发病时感到不适，而到医院检查时症状消失，进行心电图检查时得不到明显异常的心电信号，这将影响对患者的诊断和治疗。因此导致了许多心脏病人病情不能得到及时诊断和治疗，延误了病情甚至导致死亡等严重后果。因此，开发一套携带方便、低成本的远程心电监护系统具有深远的理论研究意义和实践应用价值2。其具体意义如下（1）使患者或医护人员减少了路途奔波，节省了时间和社会的医疗资源；（2）把患者的监护在家中完成，既为患者节省了开支，又为医院节省了床位；（3）使患者在熟悉的环境中进行检测，减少了患者的心理压力，提高诊断的准确性；（4）对于自理能力较差的老年人和行动不便的病人的实施远程监护，可以随时了解监护对象的健康状况在患者病情突变恶化时报警，为患者提供及时的救助。12研究现状随着电子与信息技术的不断发展及其在医疗系统中应用的深入，世界各地尤其是欧美国家相继提出了心电检测设备的小型化、家用化要求和建立远程医疗体系的设想。从1980年代开始，国外开始建立以电话线路传输心电信号的心电图监测中心，随后又出现了以数字式电话传输心电图信号的研究。英国牛津大学的JOHNSON教授采用远程监护的方法让孕妇和胎儿在放松的状态下在家中检测血压、血氧、心电图等重要生理指标；德国的一个研究小组则通过宽带视频通信远程监护家中老人的各种生理参数，以便在必要的时刻提供救治和帮助。进入21世纪后，美国和欧盟在20002005年期间各投入150亿美元和175亿美元用来进行远程医疗的研究工作，与此同时，国外各大公司也纷纷跟进，进行心电监护产品的研究开发工作，亚洲的日本在这方面也做了较大的投入，其中SONY，东芝已有类似的监护设备上市，但都价格不菲。国内在这方面的研究晚于西方国家，一个总的特点是起步晚，起点高。但随着中国经济的快速发展，人们对健康的重视程度越来越高，对健康监护产品的需求量也稳步提升，产品的应用范围从危重病人监护，发展到如今普通病房的监护，目前，很多家庭对此也提出了一定的应用需求。国内早期在此方面研究的一个比较典型的案例是清华大学在1994年研制成功的家庭心电/血压监护网系统，该系统在病人不适时具有手动按键报警功能和类似HOLTER的心电图长时间记录发送功能3。2005年6月，山东大学齐鲁医院建成了国内首家心脏远程监护中心，该中心实行24小时监护，只要患者携带的微型发射机处于工作状态，就会将心电的异常变化传输到该中心，监护中心便可以进行相应处理和预警。目前，国内生产便携式心电监护设备的厂家有很多，产品也进入了实用化，但是大多数是以OEM方式进行组装的，具有自主开发能力的较少。总的来说，目前国内心电监护产品主要特点为市场需求越来越大；技术水平和产品质量在不断提高；生产厂家多，但核心技术掌握不足。1随着中国经济水平的不断发展及与国际社会融入程度的不断加深，在这面有着巨大的发展潜力。综上所述，无论国内还是国外都对心电监护设备的研究投入了巨大的人力与物力。伴随着电子技术的飞速发展，其前景必定相当广阔。2心电监护仪的医学基础心电信号有着医学上的重要意义，我们对心电信号进行处理分析以及对心电监护仪的设计模拟都必须在了解心电信号的医学基础上进行，它主要包括心电信号的产生机理，心电信号的时域和频域特征等6。21人体心电信号的产生机理心电是心脏的无数心肌细胞电活动的综合反映6，心电的产生与心肌细胞的除极和复极过程密不可分。心肌细胞在静息状态下，细胞膜外带有正电荷，细胞膜内带有同等数量的负电荷，此种分布状态称为极化状态，这种静息状态下细胞内外的电位差称为静息电位，其值保持相对的恒定。当心肌细胞一端的细胞膜受到一定程度的刺激或阈刺激时，对钾、钠、氯、钙等离子的通透性发生改变，引起膜内外的阴阳离子产生流动，使心肌细胞除极化和复极化，并在此过程中与尚处于静止状态的邻近细胞膜构成一对电偶，此变化过程可用置于体表的一定检测出来。由心脏内部产生的一系列非常协调的电刺激脉冲，分别使心房、心室的肌肉细胞兴奋，使之有节律地舒张和收缩，从而实现“血液泵”的功能，维持人体循环系统的正常运转。心电信号从宏观上记录心脏细胞的除极和复极过程，在一定程度上客观反映了心脏各部位的生理状况，因而在临床医学中有重要意义。22心电图及心电信号的特征分析221心脏电传导过程分析心电生理学资料表明9，心脏不断的进行有节奏的收缩和舒张运动。由心肌激动产生的生物电变化通过心脏周围的导电组织和体液，反映到身体表面上来，使身体各部位在每一心动周期中也都发生有规律的电变化活动。在每个心动周期中，窦房结是心脏的最高起博点也叫一级起搏点，它发出的激动命令经结间束首先传给房室结也称第二级起搏点。房室结向下发出一条传导路，称房室束，它位于室间隔内。房室束往下又不断发左右两个束支，越分越细，最后分别形成互相交织得像网一样的结构，称普肯耶纤维，终止于心肌内。此生物电传递变化十分复杂，呈混沌态，其有序结果通过周围组织传遍全身，使身体各部位出现有规律而各向异性的电变化。将测量电极放置在人体表面的一定部位记录出来的心电信号变化曲线，就是目前临床上常规记录的心电图。222心电信号时域特征分析图21典型的心电信号如图21所示的正常心电图由一系列波群组成，各段波群反映不同阶段的心电信号变化，由于QRS波变化比较集中，所以给出了分解图。下面对每个波形点作详细的介绍1P波最初产生的偏离的波被命名为P波，它反映心房除极过程的电位变化，代表了两个心房的去极。心脏的兴奋发源于窦房结，最先传至心房，故心电图各波中最先出现的是代表左右两心房兴奋过程的P波。兴奋在向两心房传播过程中，其心电去极化的综合向量先指向左下肢，然后逐渐转向左上肢。如将各瞬间心房去极的综合向量连结起来，便形成一个代表心房去极的空间向量环，简称P环。P环在各导联轴上的投影即得出各导联上不同的P波。P波形小而圆钝，随各导联而稍有不同。P波的宽度一般不超过011秒，电压（高度）不超过025毫伏。2QRS波群心室激活产生的最大波，它反映心室肌除极过程的电位变化。正常间隔008012秒。典型的QRS波群是指三个紧密相连的波；第一个向下的波为Q波，这波不一定总是出现。QRS波的第一个向上的波为R波，继R波后第一个向下的波为S波，发生在S波后的向上的波称为R。QRS是广义的代表心室肌的除极波，并不是每一个QRS波群都具有Q、R、S三个波，一个单相的负QRS复合波被称为QS波。3PR间期从P波开始到QRS复合波开始，它代表心房肌开始除极到心室肌开始除极的时限。正常间期是012020秒，测量是从P波的起点到QRS复合波的起点，不管初始波是Q波还是R波。它是房室传导时间的一种度量，由于这个原因，它在临床诊断上很有用。基线是由波的TP段建立的T波末端到下一个P波开始。4ST间期是在QRS波群以后，T波以前的一段平线。代表左、右心室全部除极完毕到复极开始以前的一段时间。该段在确定病理学上比如心肌梗塞升高和局部缺血降低上是很重要的。在正常情况下，它用作测量其它波形幅度的等电势线。5T波代表心室肌复极过程引起的电位变化。6QT间期代表整个心室肌自开始除极至复极完毕的总时间。QT间期代表体现了心室肌肉激活间期和恢复。这个持续时间和心率的变化相反。但通常不采用QT，而采用修正QT，称为QTCQTCQT175心室率60。体表心电图反映的是心电信号的时域特性，经分析可以看出ECG信号的特征段的分界处是波形上的拐点。223心电信号频域特征分析心电图是心电信号的一种时域描述8。心电信号经傅立叶变换可用域表示。心电图对多种心脏病有诊断意义，据而推论频域分析对心功能了解和疾病的诊断可能更有价值，但必须寻找心电信号时域和频域间的对应关系和进行更精确的分析。用数字信号处理方法,将心电信号通过不同频率的滤波器，分析滤波后图形的变异,可以明了心电图各主波的频率组成,时域和频域参数的对应关系。图22正常人心电图自功率谱图22为心电信号的自功率谱图,约有26,甚至30个谱峰,自低频起,最初5个谱峰占能量5060,其中第1谱峰幅度常最大,约占总能量1/3,以后与其后依次谱峰的关系,约为FNF0N,例如第1谱峰为0976HZ,第8谱峰则约为2930HZ,依此类推。心电图波形可用一时间函数表示，周期的时间函数可分解成许多不同频率正弦信号的叠加，各正弦波都有其振幅和相位，按频率高低，表示此不同频率正弦波振幅和相位大小的图形即为信号的频谱9。1965年后出现快速傅里叶变换FFT的高效算法，以后又有多种具FFT分析功能的硬件问世，开始了生物信号与图象的数字处理新阶段。SHER等用扫频仪进行心电图频率分析，指出高于100HZ的频率对电压振幅无大作用。GOLDEN等10用FFT技术对正常人静止心电图进行谱分析，指出谱的最大值发生于基频。过去进行的频率分析工作大多着眼于与心电图仪器设计有关的频带宽度参数问题，近年虽开始注意到其可能的临床意义，但对心电信号时域与频域参数间的对应关系尚未给以应有的注意。常规心电图在临床应用上有一定的限制性6，因此寻找更敏感和更精细反映心脏信号改变的方法则很有意义，将心电信号变换为频域描述提供了这种前景。224心电信号的特点心电信号属于强噪声背景下的低频微弱信号，它是由复杂的生命体发出的不稳定的自然信号，由于受到人体诸多因素的影响，因而有着一般信号所没有的特点。1属低频信号，且能量主要在几百赫兹以下。2信号弱。例如从母体腹部取到的胎儿心电信号仅为10V成人的心电信号范围也仅为5MV。3噪声强。由于人体自身信号弱，加之人体又是一个复杂的整体，因此信号易受噪声的干扰。如胎儿心电混有很强噪声，它一方面来自肌电、工频等干扰，另一方面，在胎儿心电中不可避免地含有母亲心电，母亲心电相对我们要提取的胎儿心电则变成了噪声。4随机性强。心电信号信号不但是随机的，而且是非平稳的。正是因为生物医学信号的这些特点，使得心电信号处理成为当代信号处理技术最可发挥其威力的一个重要领域。3基于MATLAB的心电信号处理心电信号属于强噪声背景下的低频微弱信号，要对其进行分析，就必须经过去噪处理，而MATLAB作为数字信号处理的专业软件能够帮助我们很好的完成去噪分析的功能。31MATLAB处理ECG的优势1友好的工作平台和编程环境。MATLAB由一系列工具组成4。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。2简单易用的程序语言。MATLAB拥有一个高级的矩阵/阵列语言3，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行。3出色的图形处理功能。MATLAB自产生之日起就具有方便的数据可视化功能7，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。另外新版本的MATLAB还着重在图形用户（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。4应用广泛的模块集合工具箱。MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱4。一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。32心电信号噪声的来源及特点心电信号在经过采集、数模转换过程中，不可避免的受到各种类型的噪声干扰，这些干扰使得得到的心电信号的信噪比较低,甚至淹没了心电信号。通常心电信号中主要包括以下3种噪声1工频干扰主要包括50HZ电源线干扰及高次谐波干扰。由于人体分布电容的存在使入体具有天线效应以及较长的导联线暴露在外，50HZ的工频干扰在心电信号中是常见的，依情况不同，其干扰幅度达心电信号峰一峰值的050。2肌电干扰由于病人的紧张或寒冷刺激，以及因某些疾病如甲状腺机能亢进等，都会产生高频肌电噪声，其产生是众多肌纤维分时随机收缩时引起的，频率范围很广DC1000V,谱特性接近白噪声，其频率一般在5HZ2KHZ之间。3基线漂移这种噪声是因呼吸、肢体活动或运动心电图测试所引起的。稍微剧烈的肢体运动将引起心电信号波形发生改变，严重地破坏了心电信号分析的准确性。上下波动和扭曲的心电图也令医师眼花缭乱，影响诊断，其频率一般在005HZ2HZ之间。33MATLAB对ECG信号的处理过程331MITBIH简介MITBIH是由美国麻省理工学院提供的研究心律失常的数据库。目前国际上公认的可作为标准的心电数据库有三个，分别是美国麻省理工学院提供的MITBIH数据库，美国心脏学会的AHA数据库以及欧洲ATT心电数据库。其中MITBIH数据库近年来应用比较广泛6。MITBIH的数据可以从HTTP/WWWPHYSIONETORG/CGIBIN/CHART直接下载，我们可以选择文档格式的数据进行下载。332MATLAB读取心电信号MATLAB中有直接读取TXT文件的函数TEXTREAD,它的基本语法是A,B,CTEXTREAD（FILENAME,FORMAT）A,B,CTEXTREAD（FILENAME,FORMAT,N）其中FILENAME就是文件名，FORMAT就是要读取的格式，A,B,C就是从文件中读取到的数据。从MITBIH中下载的文档直接是数字格式的直接用语句ATEXTREADFILENAME就可以读取信号，并将信号存变量A中。MITBIH中的数据采样频率为360HZ，10秒钟对应3600个点，数据读取之后会形成一个36002的矩阵，第一列保存的是时间，第二列保存的是信号值。文件读取后就可以在MATLAB中进行一系列的处理工作了。333心电信号的滤波心电信号作为一种人体的基本生理信号，是心脏电活动在人体体表的表现，信号一般比较微弱，频率在005HZ100HZ范围内，幅度为10V（胎儿）5MV（成人），心电信号信噪比和频率都较低，在心电的采集、放大、检测等过程中，易受到外界的各种干扰。常见的噪声干扰有第一是基线漂移，一般是由人体呼吸和心肌兴奋所引起的，它的频率低于05HZ，属于低频干扰；其次是肌电干扰，它是由人体肌肉颤动所致，它的发生频率具有随机性，范围在5HZ2000HZ之间；第三是工频干扰，它是由室内照明及动力设备影响到人体的分布电容所引起的，频率为50HZ。6消除或减少这些干扰时识别心电信号特征和参数的前提。心电信号噪声来源不同，频率也存在差异，正是由于这些差异，对不同的信号干扰其滤波方法也不同，滤波可以用硬件实现，但实现过程相对困难，也可以用软件编程方法实现，数字滤波技术成为目前滤除心电干扰的有效手段8。1巴特沃斯滤波器用MATLAB信号处理工具箱函数设计巴特沃斯滤波器，其中有五种格式，在本次设计任务中，用到的是格式之一的N,WNBUTTERODWP,WS,RP,RS和B,ABUTTERN,WN,TYPE10。其中，N,WNBUTTERODWP,WS,RP,RS用于计算巴特沃斯数字滤波器的阶数N和3DB截止频率WC。调用参数WP，WS分别为数字滤波器的通带、阻带截止频率的归一化值，要求0WP1，0WS1。RP，AS分别为通带最大衰减和组带最小衰减DB。当WSWP时，为高通滤波器；当WP和WS为二元矢量时，为带通或带阻滤波器，这时WC也是二元向量。N，WC作为BUTTER函数的调用参数。2滤波过程1选择工频干扰滤波，选取带阻滤波器，参数设置为FS360HZ，WP149HZ，WP251HZ，WS146HZ，WS254HZ，RP3DB，RS30DB。2选择基线漂移滤波，选取高通滤波器，参数设置为FS360HZ，WP1HZ，WS005HZ，RP3DB，RS30DB。3选择肌电干扰滤波，选取低通滤波器，参数设置为FS360HZ，WP80HZ，WS100HZ，RP3DB，RS30DB。3程序设计与实现（1）高通滤波器对原信号滤波后仿真结果如图31所示。图31去除基线漂移可以看出经过高通滤波后抑制了基线漂移。（2）低通滤波器对原信号滤波后仿真结果图32所示。图32去除肌电干扰可以看出经过低通滤波后抑制了肌电干扰。（3）带阻滤波器对原信号滤波后仿真结果如上图33所示。可以看出经过带阻滤波后抑制了工频干扰。图33去除工频干扰总体的滤波前后信号对比如图34所示。图34滤波后信号与原始信号对比带噪声的心电信号经三次滤波后，噪声被滤除，效果良好。设计的滤波器能够有效地滤除心电信号中的低频干扰、高频干扰及50HZ工频干扰，利用MATLAB强大的计算功能及丰富的工具箱函数可以快速有效地设计滤波器10，大大简化了计算量。334心电信号频谱及功率谱的计算心电图是心电信号的一种时域描述。心电信号经傅立叶变换可用频域表示。心电图对多种心脏病有诊断意义。据而推论频域分析对心功能了解和疾病的诊断可能更有价值。心电信号是心电电压振幅随时间变化的信号，称为时域函数。按照函数傅里叶展开理论。任何时间函数均可分解成许多不同频率的正弦波，各正弦波都有其振幅和位相，按频率的高低。表示此不同频率正弦波振幅和位相大小的图形即为信号频谱。5任何实际信号均可通过一系列不同频率和振幅的正弦波精确叠加而体现出来。其最低频率的正弦波称为基频,较高频率的正弦波谐波均为基频的整数倍的频率,FFT分析基于计算机数学演算,可实现时域与频域的相互转换，在时域与频域之间提供最直接的联系12。对处理过后的信号进行频谱分析，程序运行结果如图45所示。图35心电信号频谱图通过对信号进行频谱分析，能够更好的把握心电信号的特征，从而更有利于对患者情况的判断。335心率的计算由于心电信号的R波算是一个波峰，我们只需要抓住波峰所在位置就能很容易算出心率，下面的程序找到了第一个波峰出现的时间T1，以及第N个波峰出现的时间T2，然后通过公式XINLV60N/T2T1,就可以算出一分钟心跳的数。进行心率的计算是对心电信号一个基本而重要的环节，通过MATLAB编程实现使得该过程简单而精确。336信号的保存和回放MATLAB中用FPRINTF函数可以对数据进行保存；回放过程只需要再次将RESULTTXT进行读取即可。4基于MATLAB的GUI界面设计MATLAB为用户开发图形界面提供了一种方便高效的集成环境MATLAB图形用户界面开发环境GUIDEGRAPHICALUSERINTERFACEDEVELOPMENTENVIRONMENT,因此用户界面所有开发工作都能利用GUIDE方便地实现GUIDE主要是一个界面设计工具集,将所有支持的用户控件都集中在这个环境中并提供界面外观、属性和行为响应方式的设置方法。GUIDE将用户保存设计好的GUIDE界面保存一个FIG资源文件中,同时还能够生成包含GUI初始化和组件界面布局控制代码的M文件。这个M文件为实现回调函数当用户激活GUI某一组件时执行的函数提供了一个参考框架。虽然也可以使用用户自己编写的、包含GUI所有发布命令的M文件也能够实现一个GUI,但是使用GUIDE执行效率更高,不但能够交互式地进行组件界面布局,而且能够生成一个用来保存和发布GUI的文件。41GUI所要实现的功能1心电信号读取。2心电信号滤波。包括去肌电干扰滤波；去工频干扰滤波；去基线漂移滤波。3心电信号的频谱显示。4心电信号心率的计算。5心电信号的保存和回放。42GUI界面的布局设计1设计4个输出窗口分别显示输入信号，输出信号，信号频谱及功率谱，以及滤波器的特性。2设计选择框包括文件的选择，滤波处理过程的选择，以及查看频谱对象的选择。3设计提示框包括当前操作，心率计算结果。4设计保存和回放功能调用保存、读取对话框。具体布局如图41所示。图41GUI界面的布局图从图中可以看出，MATLAB中许多工具采用的是图形用户界面，包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索、浏览器等。这些工具使用户更方便直观地使用MATLAB的函数和文件。43GUI对应函数的编写431信号的读取模块先清空所有图片，然后根据文件选择的情况给变量NAME赋值，并修改提示窗口的显示内容，最后调用DUQU函数。在文件选择模块中选择所要处理的文本文件后点击“信号读取”后程序运行结果如图42所示。此过程实现的是对信号的读取。图42信号读取432信号的滤波模块在读取被保存的数据后根据滤波选择框中选择的情况调用不同的滤波器函数进行滤波，并将滤波前后的波形显示。最后将滤波后的信号数据保存到变量空间。选择滤波过程并点击“滤波”后，程序运行结果如图43、44、45、46所示。图43去除基线漂移图44去除工频干扰图45去除肌电干扰图46整体滤波进行工频干扰滤波时，因其频率一般在50HZ左右，选取带阻滤波器进行滤除；进行肌电干扰滤波时，因其频率一般在5HZ2KHZ之间，选取低通滤波器进行滤除；进行基线漂移滤波时，因其频率一般在005HZ2HZ之间，选取高通滤波器进行滤除。433信号的频谱显示根据选择框的选择，调用PINPU函数对输入或输出信号进行频谱分析，并将频谱显示在窗口中。选择信号后点击“查看”，程序运行结果如图47、48所示。图47输入信号频谱显示图48输出信号频谱显示434信号的心率计算模块直接调用XINLVJISUAN函数即可，计算值在EDIT2中显示。程序运行结果如图49所示。图49心率计算模块435信号保存模块调用FPRINTF函数进行保存即可。程序运行结果如图410所示。图410信号保存点击保存即可保存滤波后的数据。436信号回放模块先用UIGETFILE记录打开文件的位置及名称，然后调用HUIFANG函数进行回放。44GUI界面运行结果打开JIANHUYIMONIFIG