基于压缩感知的心电信号压缩算法研究

基于压缩感知的心电信号压缩算法研究

ID:34873664

大小:6.33 MB

页数:57页

时间:2019-03-13

上传者:U-24835
基于压缩感知的心电信号压缩算法研究_第1页
基于压缩感知的心电信号压缩算法研究_第2页
基于压缩感知的心电信号压缩算法研究_第3页
基于压缩感知的心电信号压缩算法研究_第4页
基于压缩感知的心电信号压缩算法研究_第5页
资源描述:

《基于压缩感知的心电信号压缩算法研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

来兩?:k參硕壬学位论文基于巧缩感知的也电信号压缩算法研巧专业名觀电路与系统研究生姓名:彭娜导师姓名;堵国模 THERESEA民CHOFEGG別GNALCOMPRESSIONALGORITHMBASEDONCOMPRESSEDSENSINGAThesisSubmittedtoSoutheastUniversityFor化eAcademicDegreeofMasterofEngineeringBYPENG-HuihuiSuervisedbpyProf-.DUGuoliangSchoolofElectronicScienceandEngineeringSoutheastUniversityMa2015y 东南大学学位论文独询性声明本人声明所呈交的学位论文是我个人在导师指导下进斤的研巧工作及取得的研巧成果,,。尽我所知除了文中特别加W标注和致谢的地方外论文中不包含其他人己经发表或撰写过的研究成果,也不包含为获得东南大学或其它教育机构的学位或证书而使用过的材料一。与我同工作的同志对本研究所做的任何贡献均己在论文中作了明确的说明并表示了谢意。:日期z研究生签名與吟碑;a矿言东南大学学位论文使用授权声明东南大学、中国科学技术信息研巧所、国家图书馆有权保留本人所送交学位论文的复印件和电子文档,可W采用影印、缩印或其他复制手段保存论文。本人一电子文档的内容和纸质论文的内容相致,。除在保密期内的保密论文外允许论文被查阅和借阅,可W公布(包括W电子信息形式刊登)论文的全部内容或中、英文摘要等部分内容)。论文的公布(包括W电子信息形式刊登授权东南大学研究生院办理。'^的‘化研巧生签名:聲峰巧产师篇名期: ^巧要也电监护仪的广泛应用将产生大量的也电数据,必电压缩技术在满足临床诊断要求的条件下一,可有效减少传输和存储的数据量。压缩感知理论是近年发展起来的种新的压缩采样理论,压缩感知理论表明信号的采样率不再取决于信号的带宽,而是取决于信息在信号中的结构和内容。本文对基于压缩感知的也电信号压缩算法进行了研究。、信号的稀疏性是实现压缩感知的本质。本文结合屯电信号自身的特点,重点研究了—屯,电信号的稀疏性及其稀疏字典的构造K-SVD。利用常用的字典训练算法算法对历史也电数据进行了字典训练。通过相同稀疏度下的表示误差对比和表示误差相同条件-下的稀疏度对比,证明了屯吨信号在KSVD算法训练的字典下比固定DCT字典下有更K-SVD训练字典的启发下好的稀疏表示性,考虑到必电监护仪的实际应用姜求,。在-D的结合也电持征波识别,SV、特征波分类等技术将K无监督字典训练算法演变为有、监督的字典训练算法。,提出了针对屯电信号的基于特征波匹配的字典训练算法并用相同稀疏度下的表示误差对比实验证明基于特征波匹配训练的字典能更好地稀疏表示屯、电信号,且能更好地符合也电信号的准周期特点。DCT字K-VD训在研巧了稀疏字典的构成方法后,实现了也电信号分别基于典、S练字典、時征波匹配法训练的字典的压缩。基于特征波匹配训练字典的压缩实验表明本文的字典训练算法训练出的原子由也电将征波模板中的典型模板构成,符合医学应用要-IH也求:并且将有损压缩转换为了无损压缩B律失常数据库中全部样本的ML。对MTII导联也电数据基于特征波匹配字典的压缩实验,取得了3.51的无损平均压缩比。,,关巧词:压缩感知,也电信号压缩稀疏字典特征波匹配I AbstractAbstractThewidelyuseofEGGmonitorproducealargeamountofECOdataneeded化becompres化d.Theemergingcompres化d化nsing(CS)也eoryisanewkindofsignalacquisitionandcompressio打paradigm.Thisnewtheoryprovidesanewdirectionforthecomre巧巨CG.pionofsignalThisthesisstudies化e巨CGsignalcompre化ionalgorithmbasedoncomressedsensin.pgThesparsityofthesignalis比〇ftindamental化realizethecompressedsensing.According化化echaracteristicsofECOsignal,化is化esismainlystudi的化esparsityof也eEGG-slddseignaandhowtotrainitssparseictionary.KSVDictionarytrainingalorithmithg’mo巧widelyused出ctionartraininalorithm.ECGssarsedictionaristrainedusintheyggpygK-SVDamlorithTrouennneg.hgh比eexprimetsofsparserepresentatioerrorcontrastithsamesarsiitlevelandsarstcontrastinthesamesarsereresentationerrorconditionsitpypypp,ssesnaasarsundereanedcaranxedihowd化化化eEGGiglhsbeterpity比triditiony化化efidictionarofDCT.Thenthroughanalyzingthedeficienciesofthetraineddictionaryanewy,dictionarytrainingalgorithmisproposedbasedonECGsignalfeaturewaveformmatchingtechnolo.Theexerimentofsarsereresentationerrorcontrastin化esamesarsitlevelgyppppyshows化atthe出ctionarytrainedby化epropo化dalgorithmcansparserepresentationofECG^naeer-wesandcanuaveriglbttmeettheECGsqsiperiodicyll.Analyzin良thecompre巧ionexperimentsofthetrainingdictionarybaseontheECGsignalfeaturewaveformmatchingtechnology,theresultsshowthattheatomsofthedictionarytrainedby化eECGsignalfeaturewaveformmatchingtechnologyarethetypicalmmee*化plateofECGsignalfeaturewaveformtintheieuirementsofmedicalalications.,gqppAtthesamelime,lo巧compre巧ionisconvertedtolo巧lewcompr的sion.Atlast;alltheECGsam-hhmplesofML"leadofMITBIHarrytiadatabasearecomr的sedusingtheroosedpppmethodiinheaveraincomressionratioof3.51.,achevgtgpgKeywords:ComressedSensing,ECGSignalCompression,SparseDictionary,ECGFeaturepWaveformMatchingTechnologyn 15目录摘要IAbstractIISIll录一第胃绪11题背景与意义1.1课、1.2也电信号压缩算法研巧概况11.2.1直接编码压缩...21.2.2变换编码压缩21.2.3参数提取压缩3241..压缩感知的研巧进展概述31.3论文内容与姐织4第二章也电信号和皮缩感知原理介绍52.1也电信号简介52丄1也电信号产生机理52丄2也电信号波形简介62丄3、屯电信号的特点72.2压缩感知理论72.2.1压缩感知的数学模型892.2.2压缩感知的约束条件92.2.3直观例证23压缩感知在也电信号中的应用112.3.1总体框架11211.3.2线性观测2.3.3信号重构口241.本章小结3第兰章基于也电信号特征的稀疏字典构造15、315.也电信号的稀疏性分析131.2常用稀疏字典的构造方法5必电信号的K-3.3SVD稀疏字典训练1633-1KSVD17..算法原理332K-SVD19..算法在也电信号中的应用.213.33实验结果与分析3.4基于特征波匹配的稀疏字典设计233.4.1基于特征波匹配的稀疏字典设计构思243.4.2基于特征波匹配的稀疏字典设计25328.4.3实验结果与分析35.本章小结30第四章也电信号的压缩感知算法314-1KVD训练字典的压缩感知算法实现1.基于S34.2基于持征波匹配字典的压缩算法32m gm4.2.1压缩原理3242233..算法实现4.35实验结果与分析34.31也.电数据压缩性能评价指标354.3.2评估数据集35-4.3KVD训练字典的压缩实.3基于S验结果与分析%4.3.4基于特征波匹配字典的压缩实验结果与分析巧4.4本章小结4243第五章总结与展望5.1总结435.2展望43賴45#考5:献47IV 第一章绪论第一*绪论1.1课趣背巧与意义一必血管疾病正严重威胁人类身体健康,已经成为高死亡率疾病之。国家也血管病、于204、血管病报舍中屯1年8月8日发布的《中国屯2013》显示,20口年我国也血管病在城市居民疾病死亡构成中占41.1%,在农村占38.7%,平均每5例死t中就有2例死于必血管病,死亡率居各种疾病之首。此报告还显示平均每5个成年人中就有1人患有也血管疾病。在此报宵的新闻发布会上,高润霖院±指出,随着社会经济的发展和国民、血管病的危险因素普遍增生活方式的改变,尤其是老龄化及城镇化进程的加速,诱发屯加、,也血管病的发病人数预计仍将持续上涨。也血管病的髙发病率、高患病率、高死亡率和极高的医疗负担,已经成为我国面临的重大公共卫生问题W。实时监测高危也脏病患的也电数据,有利于提早发现病情,及时就医降低发病死亡P一1率必血管病人在发病前般并不能发现较大身体异常、。。屯血管疾病发病往往十分危一急,,旦发作就可能给病人带来极大病痛病情发展迅速且严重,并且很有可能导致病患晕巧,甚至巧死。而高危也脏病患者在发病前屯电图就会出现异常,在晕厥或挣死前W一伴随有相应的也电图表现。因此,实时监测商危也血管患者的也电国,可提早步发现屯、脏异常情况。随着人们生活和文化水平的提高,预防胜于治疗的自我保健意识得到普及,医疗服务需求随之増长。而也脏病、也血管疾病的发病率日益提髙,使人们对于也电监护仪的、需求日益増长。也电监护仪的广泛应用将产生大量的必电数据,这使得需妄传输和存储、的也电数据量同比增长。也电压缩技术在满足临床诊断要求的条件下,可有效减少也电信息间的兀余,减少传输和存储的数据量。一W压缩感知理论是近年发展起来的种在信号采样的同时完成信号压缩的新理论。这一一持性为屯电信号的压缩提供了个新的切入点,本文对也电信号基于压缩感知理论实现压缩进行了探索研究。1.2也电信号压缩巧法研巧概况屯、电监护仪对病患多导联全天候的实时监护会产生大量的也电数据需要传输和存储,将采集到的屯电数据先进行压缩巧传输,可W降低传输功耗,还可减少所需存储空间,。数据压缩技术是W减少描化给定信息集合所需要的信号空间为目的其中数据是指W一对给定信息集合的描述。信息空间般是指数据占用存储介质空间的大小、传输数据所需要的时间抖及传输数据所需带宽,。简而言之对也电信号进行压缩就是要用尽可能减少描述给定也电信号所需的数据量。屯电压缩方法有2种分类法,按压缩技术是否对信号造成失真可分为两类:无损压W缩和有损压缩。无损压缩能够实现完全恢复原始信号,不会对原始信号造成失真,无损压缩的缺点是无法获得较高的压缩比;而有损压缩不能完全恢复原始信号,会造成信一号定的失真,但有损压缩能够取得更高的压缩比。=按压缩时去除冗余信息的方法来分类,也电信号压縮算法分为大类:直接编码压1 东南大学硕±学位论文PWU缩、变换编码压缩、参数提取压缩。1.2.1宜接编码压缩直接编码压缩方法也就是时域压缩法,多是通过分析信号的时域特性,去除信号冗一W余并进行编码,从而达到压缩的目的。。般有多项式巧合、预测编码等处理方法如i"*DPCMDPuk]TPT[]AZTECIAmeiferentialeCodeModulationuminoint,l加d(',(g)>(p4i5i【]’me【]-AZoneTEnAPAAnPoonarLADTneaipochCodig)SScanlolloximationLir,(gygpp)(Wi9lHAroximationDistanceThreshold每方法pp)。直接编码压缩算法大致可W分为S类阳例,一1容差比较压缩技术:该技术使用种简单的样条来近似信号样本点,与预置的误)差值进行比较,如果在误差范围内,则认为该样本点兀余,将其删除,从而达到压缩数据的目的。在重构时则需要用样条来巧似被删除的样本点。LADT算一法是容差比较压缩技术的种,在现有的折线拟合压缩算法中,其性能最佳,161^且己有学者提出了[111LADT的快速实现算法。2差分脉冲编码调制技术差分脉冲编码调制技术与容差比较压缩技术有一):定的相似性,即两者都需要利用邻近样本点的相关性来达到压缩的目的,都使用了预测或插值技术。差。但在实现上采用了完全不同的方法分脉冲编码调制技术W当前样本点之前的样本为先验知识,对当前样本点的值作出预测,然后巧真实值与预测值之间的差值经量化后编码,从而实现压缩。DPCM算法是典型的il差分脉冲编码调制技术tl。3)滴编码技术:摘编码即编码过程按贿原理不丢失任何信息的编码。常见的贿编码有香农编码一、哈夫曼编码和算术编码等编码方法。贿编码技术般是其他数一据压缩方法的最后步,即其他压縮技术需要结含厢编码技术使用。直接编码压缩法的优点是算法原理简单,压缩速度快,但它的缺点是压缩比偏低。直接编码压缩方法在无损压缩或者对压缩比要求不太离,实时性要求较高的应用领域应用广泛。1.2.2巧换编码压缩变换编码压缩是通过某种数学变换将时域中的原始数据变换至某种更适合去兀余的其他变换域上,然后对变换系数进行编码压缩处理,重构时,对变换系数做反变换恢P01复原信号I。变换编码压缩技术的优点是压缩比高且可控是目前主流的压缩技术。但变换编码压缩方法需要做运算复杂度高的数学变换,因此压缩的实时性较差。瓜电信号的常用变换编码压缩算法有基于DCTDiscreteCosineTransform的压缩方法PUPq,基于()u24StK1Pl小波变换的压缩方法,W及基于小波包的压缩方法等。小波变换是时域和频域的局部变换,通过伸缩和平移可W实现对信号的多尺度细化分析,因此,小波变换的时域局部化和多分辨率特性特别适合分析处理非平稳的也电信号。SPIHT牌eSetPartitioningin曲erarchicalTress)算法由EZW(EmbeddedZerotreeWavelet)算法发展而来,首先在图像压缩中取得了不错的压缩效果。文献24将SP町T算法应用于也电数据的压[],2,对各尺度的小波系数进行判决和煽编码缩其原理是采用倍递减的多级口限值,把数据变成比特流,而这些比特流是按其重要性进斤排序的,编码过程可在任意时刻终止,并能根据要求确定压缩比。因此SPIHT算法得到了广泛的关注和应用。针对小波变py--换的计算量大,间相国等人提出了基于整整小波变换的也电数据压缩算法,整整映2 第一草绪论射的整数运算提高了算法的运算速率。1.2.3M担鹏缩参数提取压缩方法是提取信号的特征点或建立信号模型,在信号重构时利用这些特、征点或模型参数来重建波形,。例如李延龙基于隐式马尔可夫模型进行屯电数据压缩、该方法充分利用屯电信号的各个脉动周期的相关性对原始信息进斤处理,在高数据压缩比的情况下仍然能够很好地恢复原始数据。参数提取压缩方法中最具代表性的方法是神经网络压缩法,神经网络中的BP网络的压缩方法是应用最广泛的神经网络压缩算法。W原始信号作为训练样本,输入到BP网络中,经过学习后存储网络节点的状态值作为28压缩结果。例如文献中,整个皮缩系统采用两个结构相同的前馈神经网络,经过预[]处理后一,将个周期的也电图分为P波、QRS波、T波3个部分,分别送到压缩网络的输入端。如果此时压缩网络的输出值的误差满足给定误差容限要求,则只需要记泉该网络的隐含层输出值一,通过利用隐含层的输出值和压缩网络的输出层权重就可W得到个误差小于误差容限的重构波形,这个波形就被压缩为6个数值,可见,用这种方法可W达到巧高的压缩比。但对于训练结果不满足误差容限要求的波形则要将整个波形都储存一,因此,这种算法的压缩性能容易受噪声干扰影响起来。神经网络算法的另个缺点是BP网络的收敛速度较慢,因此训练压缩时间很长。1.2.4达渝S巧的研究进展巧述压缩感知口9]Pa】(compressivesensing,CS)理论是由EmmanuelJ.Cand知和陶哲轩iW及DavP2006一idL.Donoho咕年发表的文献中提出的,并提出了实现压缩感知的般算法流程,。传统的奈奎斯特采样定理指出若要无失真地恢复模拟信号采样速率必须是一信号带宽的两倍,,。W带宽的两倍进行采样获取的信号中定包含了大量的冗余信息一因此信号压缩的种方法就是降低系统采样率,W减少信息总量,CS理论就为此提出了可行的方案,将模拟信号直接转换为数字形式的压缩信号,也就是CS理论在采样的一同时实现了压缩。CS理论指出只要信号在某个变换空间中是稀疏的,那么就可^用一个与变换基不相关的观测矩阵将变换所得的离维信号投影到个低维空间上,并且只要通过解一D个优化问题就可W从这些少量的观测中高概率地重构出原信号,onoho等人证明了这样的投影包含了重构需要的足够信息。在该理论下,信号的采样率不再取决于PS一,信号的带宽,而是取决于信息在信号中的结构和内容。在采集信号的时候采用个能实现CS观测的传感器,从模掛到数字,采样地同时完成对信号压缩的是模拟式的压,称为模拟CS(ana),缩感知logCS这是CS理论的终极目标。但目前的传感采样技术还不能实现,因此研巧者提出了基于传统传感器的数字式店缩传感体系(出italCS),g将经过ADC之后含大量兀余信息的信号用CS理论压缩后再传输,同样也能为系统节33巧存储空间I1、降低汁算和传输功巧。压缩感知理论绘信号处理领域带来了不小的冲击,为各应用领域的数据分析指出了一一一S理条不样的途猛,。C论提出就在各个领域得到了广泛的应用如雷达探测、地P23质勘探、医疗成像、图像压缩、信号去噪、图像超分辨率重建等。一、Mamaann将压缩感知理论应用于屯电信号的压缩也取得了定的成果,ghia等人首先将CS理论引入无线体域网(WirelessBodySensorNetworks,WBSN)中,应用于也电信号的圧缩采样,并通过实验证明用压缩感知进行压缩比数字小波变换压缩方法可W延71%,WBSN的低成本、低功长也电采集终端3.的使用寿命满足了耗、超长时间的无线3 东南大学硕±学位论文、mahan屯电监护需求ian,。Mag构建了也电信号压缩采样的皮缩感知框架并采用巧优化化-方法求解sicPursuitDeNoising(BPDN)问题实现了也电信号的重构。这个研究成果为后续研巧提供了将压缩感知应用于屯、电信号压缩的参考方向。在与基于小波变换的也电一压缩的实验对比结果中,可W看出虽然基于压缩感知的压缩性能稍差,但其个显著的4D人P优点是有助于WBSN实现超低功耗解决方案。ixon等財也电信号进行时域阀值稀疏化I并采用1位伯努利压缩感知观测矩阵,实施压缩感知,指出当信号经过闲值处理后的稀疏性大于等于98%时,可W实现16倍的压缩比。这结果令人激动,但显然是htWECG牺牲信号质量为前提取得的。Zang等人利用的块内相关性,提出基于块稀疏贝叶斯学习(BlockSpareBayesianLearning,BSBL)的重构算法,并成功将该算法运用到了非稀疏的胎儿也电信号的无线传输中。但BSBL算法的缺点是计算量太大,计算时间过长,。这些研巧都成功地将压缩感知理论应用到也电信号的压缩中并且用实验和数据、证明了采用压缩感知理论进斤也电信号的压缩能有效降低信号采集端的能耗,适用于信号采集端的低功耗要求.信号都有。但上述论文中的实验结果中压缩比结果并不太理想一定程度的失真,。压缩感知理论表明信号的采样率不再取决于信号的带宽而是取决于,信息在信号中的结构和内容,结合也电。基于此本文希望能在应用压缩感知的基础上信号的自身特点,研巧提高压缩性能的方法。1.3论文内容与组织本文的主要研究内容是基于压缩感知的也电信号压缩算法研巧。本文在前人的研究基础上进一,,,步结合屯吨信号的自身特点从压缩感知H要素的稀疏性出发研究如何构造更化符合必电信号恃征的稀疏字典,再利用此字典进行压缩,W提高压缩性能。本论文的内容与组织如下:一第章为绪论,介绍课题背景、研巧意义及也电信号压缩算法研巧概况。、第二章对屯电信号和压缩感知的基本原理进行了介绍,为后续研巧奠定理论基础。、电信号中的应用流趕介绍了压缩感知在屯,对压缩感知H要素中的线性观测和非线性优化重构做了简要介绍。=一第三章就压缩感知要素中的稀疏性作了进步研巧,分析了也电信号的稀疏性和可压缩性。指出了兀余字典仍然能应用于压缩感知的理论巧据。分析比较了固定字典与用K-SVD训-VD练的字典的稀疏表示性能,受KS字典训练算法的启发,结合也电信号的实际应用,提出了本文基于持征波匹配的必电稀疏字典训综算法。第四章分别介绍了也、电信号基于K-SVD训综字典的压缩感知实现方法和也电信号基于持征波匹配字典的压缩方法,并分别进行了相关的实验和数据分析。一第五章为总结和展望,在对本文的基本内容进行总结的基础上,指出了本文进步研巧的方向。4 巧二章私电信号和压缩感知原理介绍第二章也电信号和圧缩感知原理介绍、通过了解屯电信号的基本原理,可更好地研巧也电信号的压缩方法。其次,压缩感知理论是近年才提出的全新理论,深刻地理解其基本原理有助于更好地应用这个新理论。2.1也电信号简介、2丄1屯电信号产生机理实语结么屯、房话查结—左束支后支右心室浦觸纤维2-图1也脏示意图人体的也电信号顾名思义就是人体也脏活动产生的电信号。也脏的主要功能是累血、W维持周身血液循环,这主要依靠屯肌细胞实现。也肌细胞的电活动主要由0去极化过程,1快速复极化初期,2缓慢复极期,3快速复极化末期,4恢复期共五个动作电位循pq-、环组成。参看也脏示意图2i,正常情况下,每个屯动周期由窦房结也肌细胞产生电兴奋、、、,这种电兴奋W有序方式通过屯房内的传导束,首先激活右也房,其次是左屯房,电兴奋经过房室结时,稍有延迟,然后进入希氏束,左、右束支到达浦肯野纤维,最后通过浦肯野纤维迅速激动也室壁的普通P^也肌细胞依次产工作性也肌细胞,使得相应的生电兴奋。一、屯肌细胞在去极化和复极化过程中,与尚处于静息状态的邻近细胞膜构成对电偶,使此生物电变化与有序传导结果通过周围组织传递至身体各部位,从而使身体各部分出屯、现有规律的各向异性的电变化.动周期。因此将测量电极放在肢体表面就可W记录到一的节律性电位变化曲线、也动周期所产生的电活动变化图形。为了从体表记录屯脏每,人们设计了相应的导联系统。临床上通常采用标准十二导联系统,该导联系统包括H个、:VR、VL、aVF:双极肢体导联:In、曲:H个肢体加压导联aa:六个单极胸导联VI、V2、V3、V4、V5、V6。肢体导联系统反映了也脏矢状面情况,胸前导联系统反1W应了屯、脏水平面情况。、、采用导联系统来记录也脏电活动的技术称为必电图。也电图可W反映屯脏的生理状一况,是临床最常用的检查之。也电图可^用于监测病人生命体征,还能有助于判断治、、疗药物对屯脏的影响,用途广泛。同时由于病变屯电国与正常也电图波形存在差异,分、肌梗死、析屯吨图可帮助诊断也肌缺血、屯、也脏肥大等也律失常疾病。而临床上采用5 东南大单硕±学位论文、、二导联系统长期腔测,会产生大量右电数据需要存储,标准十病人也电图,因此研巧屯电信号圧缩技术具有重要的实际意义。2.1.2私电信号液形简介一2-2P-RR-T段、T如图所示,个周期的也电波形通常由P波SS、段、Q复合波、、脏跳动过程中的不同阶段PS面对也电波形的主要组波和U波组成,分别反映了屯。下成及其特点进行简要介绍。PPRQRSST一T山波段波巧段,波溢Iii%ft.PRrQT¥誓-闻巧..片關图2-2正常也电图图解、、、P:也叫屯房除极波。正常情况下波,每个屯动周期由窦房结屯肌细胞产生电兴奋2-1起始,如图,窦房结位于右也房与上腔静脉交界处,所W窦房结的电兴奋首先传导至右也房、、,反映为也电图P波的前半部。右屯房的电兴奋经过房间束传导至左屯房,形一成P波的后半部、P,即P波反映了屯房如肌细胞的去极化过程。波持续时间般为、0-、.08011S,025mV屯.波幅不超过.。当屯房肥大时,左右房的传导出现障碍,使得P波变异,可能出现高尖或双峰的P波。P-RP--间期(或称Q间期):是指从P波的起点到QRS波群的起点,PR间期反映了、、屯室开始除极的过程、房开始除极至屯。正常情况下,电兴奋由也房传导至屯室需要经过、--房室结,由于房室结的传导速度缓慢,因此形成了屯电图的PR段PR。正常的间期持一-、、般为012020屯-续时间.S.S。当屯房到室的传导出现障碍时,则表现为PR段延长或P波之后的也室波消失。'、RS复合波:正,电兴奋向下经过希氏束、左右束支Q常情况下,同时传导至左屯室、右也室形成复合波、S,即S复合波室除极过程SQRQR反映了屯。典型的QR复合波由一三个紧密相连的波姐成,第个向下的波称为Q波,紧接着向上、高而尖峭的称为民一S ̄、波,最后是向下的波。Q民S复合波持续时间般为0.06s0.10s。当屯脏的左右束支出现传导阻滞、有也室扩大或肥厚等情况时,QRS复合波则表现为増宽或变形等。一S-T段R:指QS复合波的终点与T波起点之间的线段,理论上,般与零电位基线平齐。在这段时期内,因为也室各部分都已全部进入除极化状态,但尚未开始复极,、室各部分之间没有电位差存在、故屯,屯电曲线恢复到基线水平。但若有冠状动脉供血不一足或‘I、-I肌梗死等情况发生时,ST段常会偏离基线,并超过定的幅度范围。一T波:S之后的TQR波反映了也室复极化过程的电位变化,。方向常与R波致一 ̄T化05s025。必肌缺血表现为T波低平倒置部分导联可倒置.s,高。波持续时间般为血巧或急性屯、肌梗死的超急期则表现为T波高耸。-T间期:指从QRSQ复合波的起点至T波的终点之间的波形,反映了从也室开始6 第二章也电信号巧压缩患巧原理介绍;-除极到全部完成复极化的过程-T间期T。Q与也率、年龄及性别有密切的关系。分析Q、间期的变化,对屯脏疾病的早期诊断有益。U波:在某些导联上的T波之后可见U波,其产生机制目前尚未有定论。2.0也电倍号的特点、屯电信号是通过电极在人体体表采集得来I属于强噪声背景下提取的微弱生物电信-M-号23,ITB阻100。如图所示是标准也律失常数据库的号数据的MLII导联也电信号的前3s波形图。-0.8--0.6■0.40.-.2§貪Q--巧^.4?.-0.6々丑00.5t1.522.53巧闻T肉图2-3屯、电信号波形图也电信号具有W下特点:一、1准周期性:屯电信号随着也脏节拍有定的规律和周巧性,但是仔细观察,每)一个周期并不完全致,且巧时间产生随机变化,因此也电信号是准周期的非平稳信号。、也电信号是低幅度信号,,典型的屯最大峰峰值巧微弱性:电信号峰峰值为ImV—mV般不会超过lO。巧强噪声:也电信号中存在基线漂移、工频干扰、肌电噪声、电极接触噪声等多种噪声。H ̄之间4)低频化也电信号的频率集中在0.05z150Hz。T、hakorNV也..等人研巧了电信号的各个成分的频谱分布,他指出也电信号大部分能 ̄量都集中在0.5Hz45Hz其中QRS复合波的能量在也电信号中占最大比例,这使得QRS,波成为也电信号中最重要的波形,其往往是也电信号分析的出发点和突破点。QRS波的 ̄Hz ̄-020Hz之PS峰值大约在1间,波和T波的频率范围是0.5Hz10Hz,T段的频率范化 ̄2此0,围是.7基线漂移和运动伪迹的频带在7HzW下其中由呼吸引起的基线漂移 ̄005 ̄Hz2出,30z300化,大约在.。肌电干扰的分布范围较广在H是也电信号中常见一PqPS高频噪声,工频干扰般在50Hz(或者60Hz)及其倍频附近。2.2压缩感知理论7 东南大学硕±学位论文^压缩感知理论是信号处理领域里近年来提出的一个全新的采样理论。压缩感知提出可W利用变换空间描述信号,只需要直接采集少数包含信号全局信息的线性观测数据,一就可W通过解个优化问题从这些少数的数据中髙概率恢复出原始信号。可见,与传统的奈奎斯特采样不同,信号的采样速率不再需要大于信号带宽的两倍,而是取决于信号1MI自身的结构与内容。但压缩感知理论要求信号是稀疏的,或者至少在某个变换空间是稀疏的,同时要求观测系统与变换系统是不相关的。只有同时满足这两个条件时才能从Pq低分辨观测数据中恢复出高分辨的信号。2.2.1压巧惑知的巧学模型W原始信号序列X是R空间的一W维个有限列向量。假定信号序列X是稀疏的,或者若X在时域中并不是稀疏的,但是根据经典的变换域压缩方法可知,总能找到某个稀ww=eR巧?,疏变换使信号^稀疏,令稀疏基?*...?且满足:[,2,,1NjqnpT_iTr_qqj(2.1)X在这组正交基下展开的表达式如下:x=)强0(2.21地其中系数向量0是必稀疏的,即其中非零系数的个数K?W,或者说0中系数呈指数衰减,具有很少的K个大系数和许多小系数,这表明原信号是可压缩的。展开系数为:==目[<x,冲i>xf>Jx(2.3)写成矩阵形式,可W得到:x=ve(2).4也就是说0是X在稀疏字典W中的K-稀疏表示。一W将W个采样点沮成原始信号X,采用另个与稀疏字典不相关的观测矩阵TA=xn,对,/?m[化.化扣X进行CS压缩得到M个压缩数据(其中AA〇組成i]一一观测向量《!>。可形的每,_^象地将行都看成个传感器中的每个元素即为相应的传感器采集原始信号得到的结果.这使得J中的每个元素都包含了原始信号X的全局信25)息,数学表达式如式(.。'===x<JC>1...A()巧如斯,,护Q2.5NM£R巨R口其中々i,因此通过CS的线性压缩得到的观测向量可W被表达为=々x(2y.6)^也就是说CS压缩只要经过式(2?^^.6)的压缩观测就可抖将7个原始信号压缩为个信号,一一.其中从数学上来说,巧压缩仅仅做了个矩阵运算JC是个解线。从中恢复一32U11性方程组的问题,进步定义为CS信息算子,记为a,24)结合式子(.,式子(2乂)可W被改写为:=(2)y.7一0重构原始信号的个重要工作就是根据得到的观测向量护选取观测矩阵>和稀疏一一2一字典求解上述方程组(.7)进行重构,般情况下通过求解此方程得到个唯确一定解是不可能的,因为这是个未知数大于方程个数的病态方程,存在无穷多个解。但由于CS压缩要求0是必稀疏的,这样未知数的个数大大减少了,这使得信号重构成为了可能。得到原始信号在稀疏字典W中的必稀疏表示系数0后,利用(2.4)就可重构原始信号。为了确保重构信号的有效性和鲁棒性,观测矩阵和稀疏字典必须严格满足一定的约束条件。8 第二章也电信号和压缩感知原理介绍2.2.2达巧感巧的约束条件一Can沁S和Donoho等人论证了只要信号在某个正交空间具有足馈的稀疏性,就能PWWPU(M?W)0从少数个观测值中高概率重构该信号。因此,只要存在使足够稀疏的稀疏基存在,我们就能从M个观测值中髙概率重构原信号:C。在0是稀疏的条件下来求解式(2.7)的问题可(U通过求解最小心范数的优化问题来精确重构或者近似逼近0,数学表达式如下:CS=in0sA2m.t.0(.8)||y||。其中L〇范数是指向量中非0元素的个数。因此求解(2乂)得到的是原信号的最优0稀疏逼近,。但是为求解其最稀疏解需要穷举中所有非零项位置种可能的线性组合,一-才能得到最优解,这是个NPhard的非凸优化问题andfesonoho。C和D等人证明了在一tW定条件下,最小^范数问题的解可与最小k范数的解等价:CSmin0s=A9).t.0(2.||y||i^范数是指向量中各个元素的绝对值之和,山范数也被称作稀疏规则算子asso(Lreua)。应,glrizatkm使非凸优化问题的求解转换成了凸优化问题的求解问题变成更加便于求解的线性规划问题。同时Cand扣和陶哲轩等人论证了式(2.8)与式(2.9)求得同解的充要条件,即观测矩阵与稀疏字典必须满足有限等距性质(RestrictedIsometryPropert,RIP):y1-各々W0含1+(2.10)(WWl|(成川训2llz||z其中称为有限等距常数,《k必须小于1,且《k越小精确重构的概率越大。0II||2是系数向量0的^范数,k范数是指向量各元素的平方和的平方根。常用的观测矩阵一如独立同分布的高斯矩阵或者伯努利矩阵,有学者已经证明它们与般的标准正交基的33[]RIP。有限等距常数大都非常小,满足性质但采用独立同分布的高斯矩阵或者伯M的一努利矩阵作为观测矩阵时大小有定的要求:M>CKIoW/(2.11)kg(/〇一个很小的常数其中Ck是。U一对于其他测量矩阵和稀疏基,判定a是否满足RIP性质是个组合复杂度问题,U难度很大。文献指出只要保证观测矩阵和稀疏字典的不相关性,则a就高概率满足了RIP性质。因此,实际应用中,我们常常通过计算观测矩阵和稀疏字典的相关性片来£Sa;2是否满足要求.12)验证。;的计算式如式(。屯=/,Fmaxi<>(212)|V.sI,I(ijsN机A^;的值越小,说明稀疏字典与观测矩阵么间越不相关。2.2.3直观巧证zzzzJC由四个频率分别为50H、lOOH、200H、400H的余弦信号叠加构成原始信号,表达式如下:*=*■?扣?xs)0.3cos(27t/its75)+0.6COS口/zts了5)+(7^■I?■=1…0.1cos27ttsT+0.9cos27rtsTrW(2.1)(/(鳥?其中33s)s):9 东南大学硕±学位论文^,,,,式中乂、/2、/、K为采样间隔价j為为余弦信号的四个频率與为采样频率=7;为采样序列。/尸50Hz,尸lOOHz,佑200Hz,/尸400出,采样频率为最高频率的两倍/=-XX/5800。如图24,对做FFT变换,可W看出在频域是稀疏的。0-.908■.070.60-I-S^0乂.0■-30■,2。IlLLLLj:400-300-200-1000100200300400巧率f[Hz>图24构造信号的频谱分析ts)巧抬巧号与S巧任号的时试巧形巧比图40.010.的0.的0化DU5a扣时间T(S)"-WD脚巧始话号与重巧信号的误董曲巧—2111111|--1心^。' ̄■八?^A/wVvW\A/爸^/\^y\/Vy---1壬0.010.02D.030.M0.化0.06时邮(S)图2-5原始信号与重构信号的对比图取W为256,采用随机高斯矩阵作为观测矩阵对X进行线性观测,得个压缩观测值,稀疏基为傅里叶正变换矩阵,用正交匹配追踪算法求解式(2.7)进行信号的重-构,原始信号与重构信号的对比图如图25,其中(a)图是原信号与重构信号前50个一点的时域波形对比图,由图可见,原始信号与重构信号基本重叠在起b)。图(是原始信号与重构信号之间的误差曲线,量级在1〇?14,可见重构误差非常小,基本实现了完全重构。原始信号在时域并不是稀疏的,但是由于原信号是由四个频率不同的余弦信号叠加而成的,因此原信号在频域是稀疏的,因此可W用压缩感知算法对此信号进行压缩观测。按照奈奎斯特采样定理精确重构信号需要信号带宽2倍上的采样频率,也就是本例至少需要256个点才能无失真地恢复模拟信号。但按照压缩感知理论,此例用64个观测值就可W实现精确重构原信号,可见在CS理论下,信号的采样频率不再取决于信号的10 第二章&电信号和压纸感知原理介绍2P3带宽,而是取决于信息在信号中的结构和内容,。总之本例说明了只要信号是稀疏的或者可压缩的,就可W利用压缩感知理论,从少数包含全局信息的观测值中高概率恢复原信号。2.3压缩感知在也电信号中的应用2?义1总体框架一-、压缩感知算法流程图如图26,首先要保证能找到个稀疏字典,使右电信号在该稀疏字典下巧疏。其次需要设计相应的观测矩阵,要使其与稀疏字典满足RIP性质。信一号压缩时经过压缩观测得到的还需要经过去冗余编码,实现进步的压缩。重构恢复,,2)0信号时利用解码得到的J再通过解式(.8得到稀疏系数,已知稀疏字典就可W2、通过式(.4)重构信号得到。将压缩感知用于也电信号实现压缩和重构的流程也与一,可W直接采用的流程其他类信号致。稀疏功!WV1 ̄ ̄^*iyIy6iI重枯传导———絹解码 ̄ ̄^:Hi观巧矩阵J<s>n\’国2>6压缩感知算法流程从的压缩感知流程图可W看出,实现压缩感知有H个关键要素:稀疏性、线性观测P3和非线性优化重构。一信号是稀疏的或者具有可压缩性,是实现压缩感知的本质。稀疏性这块最重要的工作就是构造令信号稀疏的稀疏字典,例如常用的由傅里叶变换,小波变换等正交变换的正交基构成的字典。关于也电信号的稀疏字典的构造的相关内容将在后续章节进行研究。线性观测是实现皮缩感知的关键,也是实现信号压缩的关键步骤。关键是如何设计R一IP与稀疏字典满足性质的观测矩阵。般地,在线性观测后,为得到更好的压缩效果,压缩观测后需要进一步去冗余编码。本文采用的是应用广泛的哈夫曼编码。非线性优化重构则是重构信号的必要途径。如何从低维观测值中稳定恢复出高维信一号是个数学难题,信号的稀疏性使非线性优化重构信号成为了可能。相关重构算法已、趋于成熟,可W直接应用于屯电信号的压缩重构中。2.3.2线性观測4>M个行向量对原信号进行投影线性观测的过程实际就是利用观測矩阵mxn中的I得到M个投影值,实现从W个信号点到M(M?A〇个观测值的压缩。而实现压缩观测一个满足条件的观测矩阵。首先需要设计,保证在重构时能很好地恢复信号设计的主要U一原则是:(1)a妄满足RIP性质;(2)普适性,即该测量系统与般的稀疏字典或者11 东南大学硕±学位论文稀疏表示系统都高概率不相关;(3)易于实现,如计算简单,硬件要求低等。目前常用的测量矩阵是服从独立同分布的高斯随机矩阵,伯努利分布矩阵等。随机矩阵能满足上一般的固定稀疏字典W很高的概率满足R述H个设计原则,独立同分布的随机矩阵与IP*KlAW性质,并且指出0(og能将#维信号的K个最大值较島的概率稳定(g)的采样就PW重构。Candfcs还指出用服从独立同分布的高斯随机变量形成的测量矩阵和任意的超UPiKWtui完备冗余字典构成a,压缩观测信号的稳定重构仍然有很高的概率。2.33倍号■构重构原始信号就是要从低分辨率观测中恢复出高分辨率信号,是压缩感知理论的核一、屯内容之。快速稳定的信号重构算法是压缩感知理论能够走向实际应用的关键,是当一前压缩感知理论的研究热点。目前存在的主要有S大类的重构算法:第类是基于最小L〇范数的优化求解问题的贪婪算法,主要包括匹配追踪(MatchingPurs山t,MP)和正tWh交匹配追踪〇rtoonalMatchinPursuitOMP;第二类是基(gg,)及其改进算法于1>1范数4t3,主要有基追踪算法(BascPursu,B)最小化的凸优化算法iitP梯度投影算法,4"【。t](GradientProectionForSarseGPSR),内点迭代法,总变差算法(TotalVariation,jp,TV),最小角度回归算法(LeastAngleRegression,LARS)等;第;类是组合算法,主要包括傅里叶采样法PqPU,链式追踪算法(ChaininPurs山tCP),和HHS(Heavg,g2口]HitersonSteroids)追踪算法等。贪婪算法是压缩感知的信号重构算法中最经典的算法,其原理简单,重构效果好且容易实现。在小规横信号领域就可W取得较快的运算速度与全局最优解。但若信号是大规模数据,贪婪算法就需要很高的运算成本,这时往往采用基于k范数最小化的凸优一化算法、,W牺牲精度来换取速度。本文的压缩对象屯电信号是维信号,可W采用滑动窗口进行分段压缩,数据规模小。因此本文选用贪婪算法中的正交匹配追踪算法来实现信号的重构,就可W取得较好的效果。正交匹配追踪算法(OMP)是由匹配追踪算法(MP)发展而来,它继承了MP算法的原子选择思想,但又对其进行了优化。即OMP算法在每次迭代时选取原子库中与当前误差向量最匹配的原子加入到逼近原信号的原子集中,此处的匹配是指误差向量与原子库中原子的内积大小。OMP算法的优化则体现在对得到的最佳匹配原子进行了正交化处理,得到信号在正交基上的投影,再更新迭代误差。因此OMP算法最多只需要K次迭代就可W完成运算,而MP算法则无法推断迭代次数,这说明了OMP算法可W实现比MP算法更快地收敛。令观测信号为y,,稀疏度为&用于存放与残差匹配的向量的矩阵称为増量矩阵记为P,用于存放正交基索引的称为索引向量,记为r,OMP算法的具体实现过程如下:=1初始化:残差r,増量矩阵P为空,索引向量r为空。)^UU2算一一一)计r与a中每列的内积,找出内积结果最大的列,将这列从a中U一移到増量矩阵P中,也就是扩充增量矩阵的同材从a中移除这列。将最大内积向量的索引添加到索引向量r中。T-iT=3利用最小二乘法求解近似解sppp;更新残差)j,。()4最大迭代次数一)般为信号稀疏度&如果残差小于迭代阀值或达到迭代次数,则停止迭代。否则转到步骤(2)继续迭代。迭代完成后,由信号在稀疏字典下的稀疏系数和已知的稀疏字典就可W得到原信号的近似解。12 第二章也电信号和压巧感知原理介绍2.4本章小结本章内容主要分为两部分、。首先介绍了也电信号的产生机理、屯电信号典型波形、特征波形和波段间期的基本特征、、分析总结了也电信号的时频特点,为后续的屯电信号的压缩算法研究奠定基础。其次介绍了压缩感知的基本原理,压缩感知的数学模型,及其约束条件一。并用个混合了四种频率的余弦信号直观地展示了压缩感知的应用。最后介绍了将压缩感知应用到屯、电信号进行压缩的算法流程。指出压缩感知主要有H个关键部分:线性观测,、非线性优化重构和稀疏表示简要介绍了线性观测和非线性优化重构一的原理,章中展开而稀疏表示的相关内容将在下。13 第s章基于也电信号恃征的稀疏字典构造第H章基于也电信号特征的稀巧宇典构造Pq必电实现压缩感知有H个关键要素:稀疏性、线性观测和非线性优化重构。实现、信号的稀疏表示,是实现压缩感知的本质,本章将详细分析也电信号的稀疏性,探讨如何构造符合也、电信号特点的稀疏字典。3.1也电信号的稀巧性分析1IbI■■II0.-8[0怎-0■.4?02.-I0.巧-0■-6IIII40IIItI012345678910时间巧S)3-、图1屯电信号波形图-、M-如图31是标准屯律失常数据库ITBIH的100号数据的MLII导联也电信号的前1化波形图,,因为也电信。也电信号在时域并不是稀疏的但也电信号是可W被压缩的主要体现在两个方面一:也电信号是,号存在着相关性,个低频窄带信号其95%的能量、HzHz的-0,1,都集中在.25到35频带中间!也电信号是准周期信号如图3连续屯动周P31期所产生的必电信号波形十分相似。人与人之间的也电信号波形大致相同,但是根据各个人的身体素质和也血管状况会有较大出入一。但同个人的各个正常也动周期的也电一波形基本致,只在也律失常时出现较大变化。在研究也电信号的压缩方法时可充分利用也电信号送些特性来去除也电信号兀余?,构建基于也电信号特征的稀疏字典。3.2常用稀巧字典的构迸旅若信号的稀疏性越好,压缩感知算法将取得更好的压缩效果。前文已经阐明对于服*M?从随机高斯分布或者随机伯努利分布的观测矩阵,只需要〇(Kl〇g(次观测就可W))鳥概率稳定重构信号32指出由统计理论和组合优化理论可知;在满足重构。同时文献[],条件时,通过选择合适的观测方式和重建算法仅仅需要K+1次观测就可将W维空间的K稀疏信号精确地重建。因此,必电信号在某稀疏字典下稀疏度越高,离概率稳定重构所需要的观测数目就越少,重构算法需要的送代次数也越少。,,超完备兀余字典等常见的稀疏字典有正交基字典正交级联字典。正交基字典是15 东南大学硕±学位论文一指由个傅里叶变换,离散余弦变换,小波变换或小波包变换等正交变换构成的字典。正交基级联字典是指由多组正交基级联而成的字典,如Haar函数和Wa沈函数组成的级D-llO联正交字典,aubechies系列小波基dbdb构成的正交级联字典等。而在超完备兀余()字典中,基的正交性不再被保证,此时的基改称为原子atom。根据构造过程来分类,P3常用的稀疏字典构造方法主要分为两类,即无训练字典构造法巧有训练的字典构造法。一无训练字典构造方法,般是根据待分解信号的先验知识■确定字典生成函数,对生成函数做平移、伸缩、旋转、调制等运算来确定字典的原子。例如,常见的DCT字一旦确定也就不再变化,也就是不具有自典。无训练字典构造方法虽实现简单,但字典适应性,,。若要求信号在此字典下稀疏就必须先掌捏该信号大量的先验知识确保生成一一函数能准确反映待分解信号的结构,。每个人的也电信号波形都不致而且同个人的、波形每个周期也都有差异,造成的屯。各种各样的也血管疾病对也电波形也有极大影响PW电波形变异。因此远种固定不变的字典将不能很好地自适应所有病人的也电波形。有训练字典的构造方法相对较复杂一些初始原子构成初始字典原,首先选取,初始子可W由正交基构成,也可抖根据先验知识构造的生成函数确定。然后准备好训练字典的样本,用这些样本来训练字典的过程就是先固定字典,,求取表示系数:再固定系数更新字典,,直到训练结果满足要求为止,交替循环。由此得来的字典与信号特征相近因而保证了分解的稀疏性,,。与无训练的字典不同稀疏字典将不是由先验知识确定但是先验知识有助于字典的训练P81。有训练字典构造方法的关键是选取一个有效快速的字典训练算法一。字典训练算法一般分为两步:(1)固定字典,分解样本得到在此字典下稀疏系数步实际就是稀疏。这,可用己有的MP2,表示、OMP、BP等多种稀疏表示算法完成;()根据样本集更新字典原子,使原子更加符合样本特性,是字典训练的关键步骤。常用的字典更新算法有56最大似然法Pq【](MaximumLikelihood,ML)、最化方向法(MethodofOptimalDirections,口。。"--MOD)、(MaaosterioriMAP)、KSYD(KS最大后验概率逼近法ximuminularp,gaueecomoson---VlDpiti,KSVD)算法等SVDKSinlar。K算法是由均值算法和奇异值分解(guValueDecompositionSVD),是目。,扩展而来前应用最广泛的字典训练算法超完备的冗余字典不仅可W使稀疏表示更加灵活,而且能提高信号表示的稀疏度。一用训练算法训练出的字典般不能保证原子的正交性,然而在压缩感知理论提出之初,均设定稀疏字典为标准正交基字典,压缩感知理论是否适用于兀余相干字典呢?。那么这将决定我们能否采用超完备的冗余字典作为稀疏字典。文献54的研巧将压缩感知系统的稀疏字典由标准正交基字典扩展到了兀余字典。[]文献通过理论和数值实验证明当压缩理论常用的随机测量矩阵与确定性的兀余字典具有很小的有限等距常量时,那么在冗余字典下的稀疏信号依然可W通过BP算法或者OMP算法从少量的随机观测值中高概率稳定重构出来。因此,原压缩感知理论中的标WwK><L准正交基WeR可W扩展为超完备冗余字典DG民(L>N),压缩感知理论的数学表达式(2.7)扩展为:xi=^xnD(31)ynmnxlXlxi.送表明,只要冗余相干字典与观测矩阵满足RIP性质,并且信号在该冗余相干字典下稀疏时,仍然能通过常用的重构算法如OMP实现高概率重构。這为后续的实验提供了理论依据。3-.3也电信号的KSVD稀巧字典训综16 =第章基于屯、电信号特征的稀疏字典构造信号在不同字典下的稀疏度是不同的,本章希望构造出符合屯电信号特征的稀疏字典,使也电信号在该字典下表示更稀疏。通过对历史如电数据训练的超完备兀余字典将-更加符合也电信号特征VD、,选用应用最广泛的KS算法来进行屯电稀疏字典的训练。331K-SVD..算法原理一K-SVD算法是由Aharonad和Bucksen、Elrti提出的种新的字典训综算沾,也---K奇异值分解算法。KSVDK均SVD称算法与值算法有着密切地联系,当K算法中一K-SVDK均要求的每个样本只用个原子表示时,算法就退化为值聚类算法,因此可一-VD心均值算法-VDm兑KS算法就是广义的。KS算法非常灵活,可W与任何种稀疏表示算法结合-,并且KSVD算法有明确的目标函数,是,收敛效果好应用最为广泛的字典学习算法。WxS令字典a>A〇,字典的列向量为字典的原子,训练样本集Xgrj中tx一一ERS每列都是个训练样本,胡I练信号的稀疏表示系数向量集r其中L代表字典原,-子个数,S代表样本的个数。KSVD训练算法的目标函数可表示为;-<K。inXDrsi.2)m.t.VDr,|||yi|,{|口||〇其中K为稀疏表示系数中非零分量个数的上限,尸i为r的第i列的列向量。从误-VD训差逼近的角度KS综算法的目标困数还可W表示为:-<minsXDr.3)y.t.£(3Di,r|l||{||〇}||^一送两个式子实质相同,只是考虑问题的角度不同。式0.2)的求解过程是个迭代过程,D,用MP、OMP或邸算法求解样X的稀疏表示系数矩阵r首先固定字典本;然后根据系数矩阵r,更新字典D。交替进行这两个步骚,直到满足迭代条件。-,也是字典学习算法的关键32字典更新算法比较复杂。可W参看图来理解字典的更新过程。咖匪了円飞」d.^叫Jy(x>jj图3-2字典训练示意图一r和字D字典更新需要逐更新字典的列,将要更新。假设系数矩阵典都是固定的字典的第ri(注意这并不是r中第行的转j列屯令系数矩降中与4相应的第j行为ypj-2。3.3)?:置),如图3则式(中的优化式可y改写为—-=X-rx£(34)Dd!.IllHllI:比II?||III(冲灿j片广片17 东南大学硕±学位论文*心式中乘积Dr被分解为L个矩阵的和,由于除4外的其他列是固定的,则其他1个矩阵是固定的,将要处理的第j个矩阵提出来,其他项合并为E,代表的是去掉原子j和^S个祥本中造成的误差SVD分解更新,SVD能4的成分在所有。如果此时就用4yT找到距离最近的秩为1的矩阵,这能使式(3.4)的误差有效减小。但是,这样得到的y!■将不;*是稀疏向量,也就是说因为X中0的影响,用SVD得到的的更新向量中非零值的位T""置和数量将和原中非零值的位置和数量不同,并且大多数是非零值,称为出现发散。为了解决这个问题,需要将中所有的零值点固定不变,仅考虑非零值的影响再做}iSVD更新和W=0就不会再发散,,i《,<s*。为此定义集合,r为用到44片j{iU〇}-£lX"矩阵巧在的所有样本集合*的索引构成的集合,如图32。定义表示为W{}j|j怖'!)处的值为1,其他点的值都为0。定义:=xk(3.5)X=Xft(3.6)|jjE=Eft(3.7)jjj|?式(3.5)、(3.6)、(3.7)分别为、X、E中去掉零输入后的收缩结果,的长度yjj:Kk!IWW为,X和e是长度为WX的矩阵。此时,式〇4;.)转化为|j|rr|j|0d1=-4(3).8I巧广iI||j知川柄!知38最小化式(.)得到的户K和原rk就会有相同的支乾不会再出现发散。此时再对E做llSVD分解:T=UeAV(3.9){^一U的第一A令4的更新结果a,为V的第(11)后列列和,内诗,的乘积。得到诗用它更新原rk对应的非零元的值。逐列更新完宇典D后得到氏再用B做稀疏分解,并判断是否达到停止条件,没有达到就继续更新字典,达到停止条件,说明字典训练己经完^成。停止条件可1是指定1^的迭代次数或者重构信号与原信号之间的误差率。在保证稀疏的条件下,每次字典更新时,式(3.8)的误差总是减小或者保持不变。送可W保证总误差是单调不増的-,因此KSVD算法能保证收敛到局部(或全局)最小K-SVD值,。同时算法又能与现有的稀疏表示算法灵活结合因此在稀疏表示、压缩感知等领域得到了广泛的应用。18 第兰韋基于也电信号特征的稀疏字典构造3-.3.2KSVD算法在也吨信号中的应用也电监护常常是对某个病人进行长时间的监护、,因此屯电数据的压缩应该考虑对单个病人的长时间数据进行压缩的应用场景。由于也电信号是准周期信号,个人的也电波一形除屯、律失常情况外般各个也动周期的波形差异都较小,具有较大相关性。因此取单个病人的历史必电数据进行字典训练得到的稀疏字典最符合该病人的也电特征,将其应S压缩中将有利于提高压缩性能-VD用于C。本文选用应用广泛的KS算法进行也电信号心1的稀疏字典1练。KxlW"K-SVD算法应用于屯、电信号的稀疏字典训练SR>GR将,令D(王a〇,X(於bsreR分别代表字典、训练样本集、训练信号的稀疏表示系数向量集S,其中代表样本一的个数,//,,。令字典训综次数为W即个字典总共训练次是固定的本文每次训练按-VD、-照KS算法的目标函数(3.2)电信号稀疏字典的训练SVD算法流程进行屯。K图如-图33。初始化D、B、H=l.kk=lj用OMPg法计算I户Tx=〇r?计算误差矩阵Epj巧出用离4的样II本索引集含Uj-r计算误差巧阵Ei—-—-VSVD分巧—」I更新4和I=+ijj=kkkk+l图3-3K-SVD算法流程图瓜、对训练字典的性能分析采用仿真电信号进行定量分析。在仿真也电信号上加入=S个样本组成矩阵X=i=720SNR30dB随机商斯噪声构造,作为训练样本。取W360,,19 东南大学硕±学位论文-护1440按照KSVD算法公式3.2)进行宇典训练,护10,字典训练次数护10。初(,-始字典原子可W从训练样本中选取,或者使用固定字典。关于KSVD算法的具体实现EK-,本文采用Michaellad提供的SVD工具包。由于本文字典的规模并不是特别大PMchaeEadK-VD工具包因此稀疏表示算法就选择了原理简单的OM算法。ill提供的S中有OMP算法的相关函数。为了便于与DCT字典比较,此处采用固定DCT字典作为初始字典。DCT字典的构Kxl造方法是令字典DGR,L为字典原子的个数,N为原子的数据点数。按照公式(3.10)951构结:")。12???^VcosV2cosVcos。;^^!化如?1VI??coscos巧osVIc2養:D=(3.10)...些^AiV2cosVIcosVI扣sg晋豈'口w--如^^^^1V2cosVicosL2L2L乎2lj由于7>W,组成的DCT字典是兀余字典。构造的DCT字典中第1、9、27、81个一3-4原子波形如围其中第个原子为直流分量,其他原子由不同频率的余弦波姐成。,口磯-1〇0510日名1.巧间TW时间T似<c>荣巧个巧子巧形因(d)*81个巧子巧化图imm00100-515巧间TU)时间巧钟图3?4DCT字典其中4个原子波形图阳策1个巧子巧形图CD)?3个用子巧化图0?0之.61|1A0.A.l04I|I?L.-■■■.■一一.ft?々I■!0051日0.-51巧间Ts巧间T【)C糾)278(c*个巧字巧巧图<d)i?1个巧字巧化因00.£.6.1jA0-4^0.4!Il-.-?0■2〇0.51.00051巧间T的巧间Ts【}图3--SV5KD训练字典中4个原子图20 第三章基于也电信号特征的巧疏字典构造-K-VD训-进行字典训练按照图33■图35是S练字典中第9、27、811、个原子波-形图,可见,KSVD算法将初始字典中的不同频率的余弦波都训练为了也电信号波形和高斯噪声曲线。此处的训练样本就是由仿真的无噪也电信号加上随机髙斯噪声构成的。这说明了K-SVD算法能在无监督下训练出符合样本特征的字典原子。实验结果与分析此处将分别从两个方面进行实验对比:相同稀疏度下的表示误差对比和表示误差相同条件下的稀疏度对比。一一般地信号的稀疏性定义是,个长度为W的信号具有稀疏性是指信号仅有少量的,K个非零分量或者信号在某种变换基下的系数仅有少量的K个非零分量,若非零分量?W。,的个数K,则说此信号具有稀疏性更广泛地当信号中的元素或在稀疏字典下的系数呈指数衰减,仅有很少的K个大系数和许多小系数,那么也认为此信号是稀疏的。信号的稀疏性可通过下面的Lp范数来度量:1化=往(3.11)?=1邮)|叫=特别地,当>0化就是L〇范数,表对時元素的个数。常用的如2分别/扩 ̄2指^范数和b范数。P的取值范围为0。1、相同祐51成下的表示误差对比按照L〇范数来度量表示系数的稀疏度,相同稀疏度也就是表示系数中的非零值个数相同。(a)巧始巧号与*巧巧号巧比图_1I*511fI■?1fIII一+?—巧巧化号I1'1一^—i巧化号;ItIIIII>I_flC'00.1020.30.40.50.¥0.70.80.91,巧间T(s>巧始化号与重巧巧号的误爱曲线0.111111111.11-A〇.化I笠也化■??IIIIIIIIJJ1■00.10.20.30.4氏50.£0.70.80.91时间T<?)图3-6基于K-SVD训练字典的稀疏表示结果3.3,采用.2小节同等条件进行字典训练字典训练完成后,在仿真也电信号上加入NR=S30曲随机高斯噪声产生100个测试样本00个样本I。用训练得到的字典去表示这1MP=-VD训采用O算法计算表示系数,取迭代停止条件为K90。在KS练字典下,取这-100个测试样本中表示误差最大的样本,原信号与重构信号的对比围如图36,图(a)为原始信号与重构信号的波形对比图,图(b)是原始信号与重构信号之间的误差曲线?>/化>=)重构信号与原始信号之间的百分比均方误差1.71%(f皮0的定义见4.3节。用DCT字典去表示这100个測试样本,仍然取这100个测试样本中表示误差最大一-7)的个样本,其原信号与重构信号的对比图如图3,图(a为原始信号与重构信号的21 东南大学硕±学位论文()。重波形对化图,图b是原始信号与重构信号之间的误差曲线构信号与原始信号么间的百分比均方误差户化>3.96%。(a)巧巧倍号与*构话号巧比图_5■1.I<■■*—-II巧抬化号'1今I',M±iII是0S.11^'iIIIII1115001.1日之0.30乂D.6060.70.80.9巧间TU)(b)巧抬巧号与童巧巧号的误爱曲巧0.111111I11110.-05SI-0-^-05IIIItIIII_〇100.1日_20.30.40.50.60.70.80.?1巧邮(S)图-的稀疏表示结果37基于DCT字典表3.1两种字典下的差对比^最大PW)平均f/{Z)最小户化)I(%)(%)(%)-基于KSVD...训练字典171%153%134%基于DCT3.96%3.57%3.22%字典K-VD训DCT表3.1给出了这100个样本在S练字典下和字典下的表示误差统计结>==-SVD1J?D14果.53%,/化)1.7%,最小_P.3%。。在K训练字典下的表示误差均值为1最大==而在DCT字典下的表示误差均值3,?.57%,最大巧?03.96%最小/701.34%。也就是K-SVD说在相同稀疏度的条件下,基于训练字典的最大表示误差、最小表示误差、平均表示误差都小于基于DCT字典相应的表示误差K-SVD。这说明用算法训练的字典更屯、符合电信号特征,能更好地表示也电信号。2、表巧误差相同条件下的稀巧度对化2-采用3,VD训练字典。在仿真也电信号上.3.小节同等条件进行字典训练得到KS一=dBX-加入SNR30随机高斯噪声构造个测试样本。采用OMP算法计算JC分别在KSVD训练字典和DCT字典下的表示系数,使两者达到相同迭代误差后停止,OMP迭代阐值—f=X-后5l〇3。运算结束,与原始信号相比,两者戶/?2)致,都是2.14%。在两种、)3-8-不同字典下的迭代误差曲线如图,图(a)为在KSVD训练字典下的误差曲线,图K-VD(b)为在DCT字典下的误差曲线。从图中可W看出X在S训练字典下收敛更快。按照L〇范数来度量,信号在稀疏字典下的稀疏度就是表示系数中非零系数的个数。此处迭代次数等同于稀疏表示需要的非零系数的个数,即此处的迭代次数等于信号在此字-SVD53DCT典下的稀疏度。在K训练字典下的迭代次数为,在字典下的迭代次数为1-32,这说明通过KSVD算法对也电历史数据训练得到的字典更加符合也电信号特化-屯、电信号在KSVD训练字典下的稀疏度离于固定DCT字典。通过相同稀疏度下的表示误差对比和表示误差相同条件下的稀疏度对比,可发现用K-VD屯、、S算法对电信号进行训练得到的稀疏字典比固定DCT字典更符合屯电信号特征,22 第兰章基于也电信号持征的巧疏字典构造-VD训练字典下的稀疏度也更离也电信号在KS。K-SVD(a)奋于W练宇巧的巧爸曲巧如 ̄1111111[15■-W■10巧6?-?'-■--巧II‘,.II0巧)站60邮1邮120巧巧巧巧(b)空于DCTiJII练字巧的巧S曲巧20111,PJis'■\W0.-1\巧\5.■V>>〇!■..‘A0204060如100120迭代次巧图3-8两种字典下的迭代误差曲线义4基于特征波匹巧的稀疏字典设计在与DCT固定字典的对比中-SVD,经过分析得出用K算法训练的字典更符合也电也电信号在K-信号恃征,VD训能更好地表示也电信号;S练字典下的稀疏度高于固定DCT字典。这说明超完备冗余字典确实具有更好的稀疏表示性。0.4■IIII■II■1。3.-I會-IP师嗎i嗎HIp-0--.3IAIIIIIIIII^00.1020.3日-4050.60,70.80.91时间T树-3-9图KSVD训练字典所有原子波形图===、NR取iV3如,1720,护1440S30犯,在仿真屯电信号上加入随机髙斯噪声构造SX-D算法公=个样本組成矩阵,作为训练样本SV.2)进,K10,。按照K式G行字典训练伴-=10,采用DCT字典作为初始字典,得到KSVD训练字典字典训练次数。字典训练一=--完成后,将所有1720个原子重叠展示在同张图上,如图39,结合困35,可W看出23 东南大学硕±学位论文一一K-VD算法用必电S信号训练出的字典中,原子分为两类,类是也电信号波形,另、类是随机小信号。这个结果符合算法原理,因为此处的训练样本是由无噪仿真也电信号加随机高斯噪声构成的。但是训练出的字典中大部分原子是重复的波形,这些过多冗余是没有意义的。且用一K-SVD算法进行字典训练,计算量大,训练时间长,字典旦训练完成就固定不变了。一、、、而也电监护般是针对屯血管患者,其也电波形将经常发生变异,变异的屯电波形又反,字典长时间固定不更新,不能有效反映也律失常情况映了也血管的病变情况。因此。3乂1塞于特征液匹巧的话巧宇典设计构思、、,结合屯电信号的实际应用,希望字典训练算法能将也电信号分为各个恃征波使字典原子由各类特征波典型模板-SVD训练字典、各类疾病的典型特征波模板构成。而K、,电特征波并进行分类方法是无监督的字典学习算法不能在训练字典时自动识别屯。本一一文基于这需求,设计了种新的针对也电信号的字典训练算法。屯、电图波形的形态直接反映了屯、血管的健康状况因此恃征波的识别与特征参数提,stw处理算法常把必电信号分为P波取是必电自动诊断的基础。也电信号、QRS复合波、T一定的电生理学含义波等特征波及其间巧组成。其中的每个特征波或间期都有。人体正常的也电波形具有准周期性,形态类似,若某个也动周期出现了异常波形,则可W将此异常波形作为疾病判断的重要依据IW。本文W此为出发点,希望在训练字典的同时提取出样本中各类特征波及巧变异波形。特征波IS1II」I,[^识别电1iJliiJ[i广II叫而也电监护仪早_」传输1I.一....—.._.一?一一.一—图3-10常规也电分析与压缩流程图也电监护仪不仅要能稳定采样也电信号,还要实现对也电图的动态自动分析,W实^21现监测屯、血管疾病患者并能发出危险预普的功能。也就是说也电监护仪在实现了稳定、采样屯电信号和稳定传输也电信号给远程监护中也后,还要实现动态也电图的自动诊断。、-如图310,常规的屯电分析与压缩流程图,图中虚线框中的是也电监护仪需要完成的工作,图中点线框中的是远程监护中也需要完成的工作,常分为预处。在也电信号处理中理、特征波识别、特征参数提取、自动诊断等几部分。如果将也电信号的压缩重构与特征波识别,可W让、特征参数提取、自动诊断联合起来考虑本文的压缩重构算法更有实际意义。24 第兰章基于屯、电信号特征的稀疏字典构造广品:U^i电igirisir^1!.I"__i_传输'主商翁^^""'画掛5暮' ̄—IIag-M-ierI、图3-11改进的屯电分析与压缩流程图3-、如围11,将屯电信号的压缩重构与也电监护仪上初步的特征波识别、特征参数提取,、自动诊断联合起来考虑直接通过字典训练算法提取也电信号中的各持征波纪成稀疏字典、。压缩时就将也电信号分为各个特征波的线性组合,传输压缩也电信号给远程监护中必的同时I就相当于传输了该波段的特征波组成信息。在重构时远程监护中也就能一、获取每个也动周期的特征波组成成分,直接为进步的高精度屯电自动诊断打好了基础。-与图310相比,这个设计构思更好地利用了系统资源,精简了系统流程,可有效降低系统功耗。3.4.2基匹R的Hit字典设计、-1按照框图31所示,在训练稀疏字典进行压缩化可W获取屯电监护仪为了实现一些也电信号处理结果自动诊断进行的,如特征波的识别、特征波位置标定等。但本文,仅仅是进行必电信号压缩算法的仿真实验,并不是基于某个具体的平台所W实现压缩一前必须先获得这些先验知识:第步,。字典设计的基本原理如下对采集到的也电信号、二步,,,进行预处理和特征波识别;第在定位屯电特征波后计算相应特征参数对特征波模板进行分类,提取各类特征波模板的典型模板构造特征波字典原子。-<■>个巧巧屯电巧形【b)pg#巧.1-1nI言一0.5^0.5^香Am*w0—A00-5100.51时间T(s〉巧间T(s)(cQRS(d)T)巧巧巧巧模巧A1A1>IE一一0.5O.S>>mH3A巧巧/\—0。I,00.5100.61时间T(s)时间T(s)图3-12也电恃征渡模板提取示意團25 东南大学硕±学位论文^必电特征波识别算法己经有大量的研究结果可用可直接采用已有的屯、电、特征波识别算法,电信号中的屯电特征波模板提取出来构。定位各特征波后将待压缩屯Nxs成样本矩阵X巨r,W为每个模板的长度,S为样本中模板的个数。此处仿真也电一-12P信号来说明构造X的方法:如图3,个模板的长度W为RR间期的平均值,波模、动周期提取的P波组成板由该屯,位置仍保持在P波定位点处,其余部分补零;QRS模板和T波模板构造方法类似。获得特征模板W后,需要通过也电信号的波形分类方法筛选特征模板中的典型模板,将之提取到字典中一。也电信号的波形分类方法般有基于时域波形形态结构的分类方法和基于持征参数的分类方法两种、。由于基于屯电信号时域波形的分类方法非常符合临床医生的诊断习惯tW,因此本文将根据也电的时域波形形态进行模板的分类和提取。本文一屯、选用模板匹配法进行模板的形态分类和提取,必电信号的模板匹配方法般是对电特征波构成的矢量按相似性进行分类。对于两个不同的模板和63=13/?.../?片定义为:巧的4,9,,如\表示它们相似程度的相关系数-而-跑脚W>如均巧)(3.12),狂2巧0)甸?=1(邮广巧户)J==式中巧i雄imw试中巧雖i?n均。吉吉采用模板匹配方法选择出的典型模板的含义是在也跳周期中多次出现且相似度高一一一的特征波形。典型横板需要满足两个要求:方面个典型模板代表着类特征波或者一一、跳周期的信号能由少量的典型模板线性组合得到方面类变异特征波,使不同屯,;另不同的典型模板应该尽可能地不同,用模板间的相关系数片小于给定的模板形态分类闲<值片〇来进行分类提取,即片^/町巧公〇。)K-SVD与字典训练算法训练出的原子多是冗余原子不同,本文提出的字典训练算、血管疾病变异波形的法训练出的原子应有实际的医学意义,应当是不同特征波或不同屯PR典型模板,如波典型模板、QS波典型模板、T波典型模板,或者同种特征波的不同一屯、3血管疾病变异波形的典型模板。1^〇民波为例,训练出的字典中除了包含般的正RS一R常Q波典型模板外,还可能包含种或几种也血管疾病变异波形的QS典型模板,64R-S变,313,(a)nodaunc如文献[列出的常见的Q异波形如图图为交界性逸搏(ltional](j)escaebeat,NESC)(b)(rematureventricularontractions,PVC),p,图为室性早搏pc图(c)为左束支传导阻滞(leftbun出ebranchblockbeat,LBBB),图(d)为右束支传导阻滞(rightbundlebranchblockbeat,RBBB),国(e)为起搏也跳(pacedbe化PACE),f、图()为起搏融合屯搏(fusionofacedandnormalbeat,FP)。p训练字典时,根据模板匹配法从模板样本X中提取典型模板组成字典D,并计算出、相应的系数向量,字典更新应使屯电特征波模板集X在减去相应的字典原子后残差的Wxi_^范数最小。设字典DeR,N为模板的长度,L为字典中原子的个数,相应的系数矩^sTGR,K。阵稀疏度用表示本文提出的基于特征波匹配的稀疏字典训练算法,按照压缩感知理论可归结为如下的优化目标函数:X-〇rsi<K.t.V(.1)||,|y|33i{|口||〇26 第H章基于屯、电信号恃征的稀疏字典构造PVCPVCPVC/W.^本洛麥kJII图3-13常见QRS变异波形)(图(a)为交界性逸搏,图(b)为室性早搏,图(C为左束支传导阻滞.图(d)为右束支传导阻tw滞、),困(e)为起巧也跳,困(f)为起搏融合屯搏r的更字典D和新是逐列进行的,由于必电,。同时信号是准周期信号特征波在每个也跳用期的位置会有微小偏移,。在不考虑也电的准周期性提取模板时并不能实现特征波的自对准,这会造成模板匹配偏差。在模板匹配时,将使典型模板进行微小位移来实现与样本的最佳匹配。则第A个原子&可表达为典型模板关于时间的函数:—…d=nitAts.tAt=t(3.14)k的k()n?P根据第i个样本x更新字典时,需要计算最佳匹配系数。若可将也电信号分为波、i、RS复合波、T波H个特征波,按L范数度量,屯电信号在基于持征波匹配的字典下Q〇K应为X中每一一3,稀疏度。但由于此处的模板样本集个模板由个持征波组成故此处一法如下w。fk:系数向量中只有个非零值ri的计算方-三…=…=max;t,mk(tAOs.tfc1,2,LandAtt(3.15)p(也ji伯i^(当字典D中存在与样本模板A匹配的原子时,即当逐个模板更新字典的优化目标函数如下:_=*-…-"化)mtAt12Land么3s.tfct=<:t(.m化k().,,.化|時。2n公式(3.16)表示逐个模板更新字典时,字典更新算法应当寻找出字典中第几个原,子在位移量为多少时与X,4位移偏移量4沁匹配i最匹配并记录下相应的原子序号巧 东南大学硕±学位论文系数化。而当rk<^〇时,则表明字典D中不存在与样本模板JCi匹配的原子,将;Ci作为,新的=原子加入到字典D中,需更新王位移偏移量却,匹配系数八l,i此时x中第/个向量的残差将为o。由于本文提出的算法字典是可(^^动态更新的,因此对应的系数矩阵也需要适时补0。记取得最佳匹配时的位移偏移量为^fmax,将字典原子更新算法伪代码总结如下:基于特征波匹配的字典训练算法Nxs输入:XGR,p〇,位移范围[ti,tn],共n个点。Sxi输出:特征波字典D,系数矩阵r,位移偏移向量JFkER。D=i=1初始化:,0,扛。[]=for/1:1:5=0化;1if王^0()-=…=Xm…=ywmaxpi,k(tAt))S.Lk1,2,LandAt。,【2,tn;(()endifyd<^〇)(i=〇阳,础王=王+1;k=L;AT^,irln〇;]=rkii;-=〇-ri:1:^1i/;j()=-iii-irk::/)(/;jelseAT=AT^,/m^[kax;]=〇=r(h:\L^ky,jJ);endend^-法在进斤字典训练时KSVD算,并不要求获取信号大量的先验知识,是无监督的字典训练算法,。而基于特征波匹配的字典训练算法需要大量的特征波识别、特征参数提取等先验知识,是有监督的字典训练算法。3.4.3实验结果与分析若可将必电信号分为P波、QRS复合波、T波三个特征波,按L〇范数度量,也电信号在基于特征波匹配的字典下稀疏度K应为3,3。在DCT字典下非零系数只有个。将不能有效表示仿真也电信号为了展示基于特征波匹配设计的字典的稀疏表示效果,K-SVD训练字典将分别用、基于特征波匹配法设计的字典,在相同的稀疏度下去表示仿真也电信号。R=SX在仿真必电信号上加入SN30证随机高斯噪声构造个样本组成矩阵,作为训一、练样本。为了仿真屯电信号的准周期性,每个样本都进行个随机的微小位移,偏移量-W==-取值为10ms10ms。取360,正720,护1440照KSVD3.2)[,],按算法公式(进行字28 第s草基于也电信号時征的稀疏字典构造==-典训练,K10,字典训练次数好10,令初始字典为DCT字典,通过训练得到KSVD训练字典。按照基于特征波匹配算法思想从样本X中提取出典型模板,由于坦里的样本是由仿真也电信号构成,不含也律失常变异波形Jl,。巧t提取出的字典原子只有H个分别是PRST、波典型模板、Q典型模板、波典型模板。用此字典去表示仿真屯电信号时,需要将字典原子进行微小位移W寻找最佳匹巧。因此,此处可将H个典型模板分别进行偏移量-为[10ms10ms,组,,]随机位按成冗余字典然后再用此冗余字典通过OMP算法去稀疏表示测试样本一。考虑到也电信号有基值,令冗余字典第个原子为直流。=在仿真必电信号上加入SNR30犯随机高斯噪声产生100个测试样本,并将每个样一--1s1sVD训本都进行个随机的微小位移,偏移量取值为[0m0m。S,分别对K练字典和]基于特征波匹配法训练的字典采用OM=P算法计算表示系数,取迭代停止条件为K4。-VD训在KS练字典下,取这100个测试样本中表示误差最大的样本,原信号与重构信-号的对比图如图314,图(a)为原始信号与重构信号的波形对比图,图(b)是原始信=?£)9号与重构信号之间的误差曲线.的%。。重构信号与原始信号之间的百分比均方误差巧(a)巧始化号与重巧信号巧比图—■1.61111T111II1——B'巧始宿号--1一扫宿号HIf-iIIyy—.…cI?I■?IlIII^'00.10.20.30.4OS0.6化7化B0.31时邮u)bC)巧始宿号与意巧巧号的误餐曲线A王AZq01〇_20.30.40-50怎070.8031时间T(S)3--图14准周期化真信号在KSVD训练字典下的稀疏表示结果(a)巧始巧号与重巧巧号对比图 ̄''''I巧抬k号I^-一 ̄1〇1兵倍号>H|t0-,511—。‘?(4I111111111500.10名0-30.40.50.60.70.80.91时邮(S)(b)巧巧巧号与*巧宙号的误金曲线01.1.....1.1.10--.的^-0--1.05'00.10.20.30.40-508070.8091时邮(S)图3-15准周期仿真信号在特征波匹配字典下的稀疏表示结果29 东南大学硕±学位论文在基于特征波匹配法训练的字典下,取这00,1个测试样本中表示误差最大的样本原信号与重构信号的对比图如图3-15,图(a)为原始信号与重构信号的波形对比图,图(b)是原始信号与重构信号之间的误差曲线。重构信号与原始信号之间的百分比均方误差f/?D=3.39%。表3.2100测试样本在两种典下的表示误差较^^最大平均f/gg(%)最小(%)III-基于KSVD训练字典9.89%4.26%3.21% ̄基于特征波匹配法的字典3...39%314%288%I统计这100个测试样本的表示误差,结果如表3.2,从表中可得知在非零系数个数4K-VD训为的条件下,基于S练字典的最大表示误差、最小表示误差、平均表示误差都大于基于特征波匹配字典相应的表示误差:(1)基于巧征波匹配。这说明了两个结论()算法训练的字典可W提高也电信号稀疏表示的稀疏度;2基于特征波匹配算法训练的、字典更符合屯电信号的准周期特性,能有效表示准周期也电信号。义5本章小结本章分析了必电信号的稀疏性和可压缩性。指出了冗余字典仍然能应用于压缩感知的理论依据、K-VD,介绍了用于屯电信号压缩的稀疏字典的构造方法S。在阐述了字典-训练算法的基本原理后,将KSVD字典训练算法应用于也电稀疏字典构造,通过实验证明用K-SVD字典训练算法构造的也电稀疏字典比固定DCT字典更加符合也电信号特-SVD点,。但K训练的字典原子多是重复的;且字典训练完成后就固定不变了字典长、时间固定不更新,不能有效反映也律失常情况.电信号分析。受此启发在详细分析了屯与压缩的具体流程后,本文提出应将也电信号的压缩重构与也电特征波识别、特征参数提取,直、自动诊断联合起来考虑接通过字典训练算法提取也电信号中的各特征波狙成屯、电字典,。并详细介绍了基于特征波匹配的字典设计构思及其设计原理同时给出了本文的字典训练算法的目标函数及其算法伪代码。30 第四章也电信号的压缩感知笠法第四章也电信号的压缩巧知31法实现压缩感知有H个关键要素:稀疏性、线性观测和非线性优化重构信号具有、稀疏性或者可压缩是实现压缩感知的本质,第H章重点研巧了如何提高屯电信号的稀疏-巧,使基于压缩感知的也电压缩取得更好的效果SVD训。本章将基于K练字典和本文设计的字典分别探讨如何实现也电信号的压缩。在给出压缩实现方法后,进行了相关实验,并对实验结果进行了分析。4-.1基于KSVD训练字典的压缩感知巧法实现K-SVD训3-3基于练字典的压缩感知算法流程图如图,图左边是进行数据压缩的基本流程,右边是重构流程。开始之前,需要准备好进行压缩的也电数据,本文先W仿真的也电数据作为实验对象,与基于固定DCT字典的压缩感知算法实验结果进行比较。再用实际也电信号M-ITB阻也律失常数据库的也电数据作为实验对象,验证实际也电-信号基于KSVD训练字典的CS压缩重构效果-B旧屯、。而MIT律失常数据库中的数据是长时间段殺据,不能直接进行如此大维度的库缩重构处理,因此,本文采用滑动窗口进行分段压缩重构。初始化时,需要确定谓动窗口的长度W,字典原子个数I,CS观测点个数和=AM确定了,也就确定了压缩比0T/M。同时,设定OMP算法的送代停止条件:可设置稀疏度&使每个样本确定由K个原子构成;或者设定迭代误差阔值f,使迭代误差f后才停止迭代低于。开始K-)JflSVD#a训练字典0初始化;I少计算CS算子X-…一—读取必电信号:i—'IOMP算法实现稀疏编码得r一产生观測矩阵oi',*X=Dr]=—Y<>ixT4--图1基于KSVD训练字典的压缩感知算法?一一<!>观测矩阵是由随机髙斯向量构成,个扁矩阵本是是,但般认为从JVXAT的随机高斯方阵中随机抽取M个向量构成?的。将原始也电信号X经过观测炬阵线性观测后就实现了从W个数据点到M个数据点的压缩。用K-SVD算法训练符含也电特征的稀疏字典己经在第H章进行过详细探讨,并给-SVDhae出了則练字典的流程图。关于K算法的具体实现本文采用MklElad提供的ssUiiK-SVD工具包。重构时,首先要计算a,在已知观测矩阵和稀疏字典时,这仅是做UM矩阵运算。有a后就是利用OP算法进行非线性优化重构。第二章已经对OMP算法Pq-ll的基本原理进行了介绍,本文直接采用MichaeEad提供的KSVD工具包中相关的31 东南大学硕±学位论文OMP函数实现。通过OMP得到稀疏系数后就可W通过线性组合重构信号。4.2基于特征波匹祀字典的压缩算法K-SVD训练的字典原子多是重复的,;且字典則练完成后就固定不变了字典长时K-VD训间固定不更新,不能有效反映必律失常情况S。针对练字典的不足结合对常规也电分析与压缩流程的分析,本文提出了基于特征波匹配的字典设计方法。4.2.1压缩原理根据压缩感知的理论,利用本文提出的字典训练算法,仍然可W实现基于压缩感知、4-2的屯电信号压缩,皮缩重构流程如图。按照库缩感知理论,用非线性优化重构算法进行信号重构的条件是:信号是K稀疏的且应满足K<M?£。因此,要基于特征波匹配、D必须是超完备的兀余字典法训练的字典实现基于压缩感知的右电压缩,训练出的字典。而基于特征波匹配法训练出的字典原子剔除了过多的冗余,针对单个病人训练得出的字,不能构成冗余字典典仅仅包含了此病人的也电特征波的典型模板及其变异模板。为了训练出超完备冗余字典,必须对大量病人的大量也电数据进斤字典训练,训练出的字典包括了各种病人的也电特征波的典型模板及其变异模板,。这种方法虽然理论上可行但实际工作量太大,太复杂。''C开始)」用特征波匹配法 ̄ ̄训练字典0___^初始化;I^^计算CS算子读取屯、电信号X1■iOMP算法实现I'巧疏编码得。产生观測矩阵本I.i占,*X=o*r1Y=C>X ̄i图4-2基于恃征波匹配字典的压缩感知算法一进步探讨压缩感知理论,压缩感知理论中在重构时需要利用稀疏表示算法求取压U缩观测值在a中的稀疏表示。而基于恃征波匹配的字典训练算法在进行字典训练的时候,就计算出了也电信号模板由哪些特征波原子组成及其相应的系数。也就是说,基于特征波匹配进行稀疏字典训练的同时完成了线性观测和稀疏表示,将相应信息传输给远,重构时就不再需要进行非线性优化重构汁算了,在实现基于特征波程监护中必。因此一匹配字典的屯、电信号压缩时,不需要再生硬地套巧压缩感知的般流程,可W在字典训练的同时完成屯、电信号的压缩。4-3基于特征波匹配的训练字典进行也电信号压缩的桓图如图,将字典训练得到的ED、字典D、残差向量rrata、系数向量中的非零值及其索引经过贿编码后传输给远屯电,必电远程监护中也就可W实现信号的重构,程监护中也。基于压缩感知实现压缩时都一*是有损的压缩x-Dr,重构信号有定失真。本文将稀疏表示后的残差进行了烟编码32 第四章屯、电信号的压缩感知算法后再传输给远程监护中屯,将皮缩感知方法的有损压缩转换为了无损压缩。X去基值梓征波传输字典D、r重构_巧__?、」.紹拉se=aBa识别更新ErrData编码支Dr+ErxDat+化se图4-3基于特征波匹配字典的压缩桓图本文采用的煽编码算法是常用的哈夫曼编码算法。哈夫曼编码使用变长编码表对数据进巧编码,需要首先评估信源中信号出现的概率,出现概率高的字母使用较短的编码,而出现概率低的则使用较长的编码。这使得编码之后的字符串平均长度最低,具有接近理论极限的压缩比、。但是变长编码表给实际屯电信号的编码和解码实施造成了困难。本6一5文采用文献[提出的基于阶差分的H字长哈夫曼编码,将码长度固定在凡个固定的]长度上,当然这是W牺牲压缩率为代价,换取应用实施的便捷性。基于一阶差分的H字长哈夫曼编码的基本原理是预先通过实验得出信号的概率统计表,将妈长固定在四个固定长度上,出现概率高的数据给予短码,出现概率低的数据给予长码一一。为保证编码的唯性可译性,在固定码长前面加上唯的标志码。分段数值编码表如表4.1。4.表1H字长哈夫曼编码方案^ii编码长度 ̄^I‘标志位方式数据值(2进制)数据位码长I,---21910a,之[)---?1。1。1_*__22210b,[川分段a-ib-ia-】2+2,)[2c-i-a_-_数值_2_2bi-2i-2bi-2aiU110c,[)-1*-1C-1-i3-1编码6b口+2222,++)b-i-—-322^U11112c-ib-a-+i+iU口22,2)从表4.:,ac.1中可W看出完整的编码长度为标志位妈长加数据位码长其中6为,,一3组编码的码长组对超过上述H段范围的数据不压缩,直接采用原有的码长进。最后-、,,MITBIH12,6行编码保证信号不失真的也电数据是位的码长最高位为符号位,。AC的值需要根据对信号的数组出现概率的统计结果进行确定。4.2.2*法实现ThakorN乂等人硏巧了也电信号的各个成分的频谱分布,他指出也电信号大部分能量都集中在0-、.5出45HZ,这S,其中QRS复合波的能量在屯电信号中占最大比例使得QRtwi波成为也电信号中最重要的波形,其往往是也电信号分析的出发点和突破点。因此,本文也将QRS波作为基于特征波匹配的字典设计的出发点和突破点,算法实现时仅仅、RS研究了提取屯电的Q复合波作为字典原子,为基于恃征波匹配的字典设计做出初步、的探索,为后续完备的屯电特征波字典研巧奠定基础。一进行压缩的第、民波可直接使用前人的研巧成果步是要识别出屯电信号的,例如,出的基于联合阐值的QRS复合波识别算法,在加入回溯处理后正确识别率高达99.74%。为了便于压缩和计算,将提取字典原子算法进巧了简化:(1)不将QRS模板提取、出来组成特征模板矩阵,而是直接利用滑动窗口方法实现逐个R波的巧缩。在定位屯电33 东南大学硕±学位论文一R-民波后,由于QRS复合波持续时间般为0.060.10S,最大化2s,因此1取波前后各一RR60ms波形作为个QS特征波窗口,令S等于波个数2)。(由于民波的定位采用的是M-ITBIH数据库给出的专家标注,定位非常准确,因此计算相似系数时没有将恃征模板进行微小位移,寻找最大相似系数3)压。(缩开始时从初始时间段中的模板中选择一一个与其他模板相关系数最大模板,作为第个原子加入字典中。开始取出与Xi相关系) ̄Jri数最大的原子化 ̄=初始^il!jxik<pI,d()^^读取也电信号]V ̄±""""*去除基线N1n将XI加入D中NS111—,也电*-R波识别n==JLblkL‘,1N^s-ErrDataixi()夺J典D是否\、dkY('‘叫空—保存k、丫、\(i)ITki.ErrDataiY)(()计算初i信号段^;++中各个QRS横板i的相关系数矩阵"1zt取该段信号中与其他各弧S模板;最相关的作为「飛柿ItTi迸行视始模扳一阶差分的H字I__^I长啥夫曼编码Sj.L=L+1 ̄结束J)图冬4基于特征波匹配字典的压缩算法流程图-按照本文设计的字典更新算法,进行也电信号压缩的流程图如图44。本文的实验M-、、对象ITB阻屯律失常数据库中的屯电数据都提供了民波的医学专家标注与分类,所W可W直接根据数据库中的R波标注进行压缩。压缩前,首先需要去掉信号的基线。初始原子由信号的前1化数据中与其他各QRS模板最相关的模板构成。在获取初始原子口进。口之后,按窗行特征模板的提取和也电信号的压缩窗外的部分都作为残差,不作处理:对/S,。字典和系数更新方法是第个QR波Xi先取出D中与A相关系数最大的.1)原子成,按照式(32)计算,U,成,若公(斬成)小于给定的模板形态分类阐一一值片〇,则将jc,ki作为个新的字典原子加入D中最后列更新为现在D中原子个数。*、屯、令^产片(A,成)i,同,将当前屯拍减去n成得到残差时需要记录当前拍所用到的字典原子序号k,系数rw,和当前也拍的时间坐标。当所有也拍都完成了特征模板的匹一配后、,对残差进行阶差分哈夫曼编码,至此就完成了本文基于特征波匹配字典的屯电34 第四韋也电信号的压缩感知算法信号压缩。^^;重构时,己知字典〇和系数矩阵11及,将残差1每个如拍的时间坐标解码后叠加上*Dr和基值就完成了信号的恢复。4.3魏结果与分析:本文的实验主要在电脑上进行算法仿真测试,平台为宏碁笔记本电脑,其配置如下CPU为Intel?Core(TM)2,双核CPU,主频都是2.10GHz;内存为2GB;32位Windows7旗舰版操作系统。所用的软件平台是MATLABR2013b。4.3.1私电巧拒压巧性化勝巧巧一、电压缩算法个性能卓越的屯,不仅要求算法能取得较高的压缩比,而且要求算法能实现较高的重构精度。但送二者是矛盾的,提島压缩比的同时会降低重构精度,加大、re波形失真电压缩算法通常采用两个指标:压缩比(ComssionRatioCR)和百。评价屯p,…eenaeoo-Mean-SuareDerencePRD)分比均方误差(PrctgRtqif,。压缩比是最常用的压缩效率评价指标,能直观的反映算法的压缩效率。压缩比的定41I义如式(.)。但在实践应用中需要注意,即使采用相同的压缩算法C及并不是固定不变,它会随着待压缩信号的待压缩信号的信号质量、采样率和量化精度等因素的变化而一变化,。因此单看某个信号或某段时间的压缩比并不能总体反映算法的压缩效率般给出的是统计平均压缩比。原始巧据比特致^—(41)^岳缩后极据比特巧?是用来评价重构信号与原信号之间的相似性的指标,?的定义如式(4.2)。2X一班i(nrg阳rec做pRn_J理品品式中表示原始信号,式中Xi表示对应的重构信号,式中N表示样本点数。recf)4.3.2船致巧集针对也电信号研发了各种分析处理算法后需要评估数据集来验证算法性能,本文的一一实验数据分为两种,种是便于进行定量分析的仿真必电数据,种是标准也电数据库6^M1^的数据SI/AAI£〇13:200268。本文参照美国国家标准AN和文献[构造的仿真也电波]、民-形如图45,图(a)是仿真的无噪也电波形,有典型的P波、QS波、T波,但忽略了=不常出现的U波,图(b)是在无噪波形上加了SNR30dB的高斯白噪声后的波形。标准也电数据库是权威医学专家标注和确认的也电数据集,它被广泛应用于临床研究中,除了用于测试&电算法性能外,还可用于评估市场上各种也电图机的性能。目前、-B常用的标准也电数据库为MITIH也律失常数据库,大多数研究者都采用它来评估算-法性能,本文也采用了MITBIH也律失常数据库的数据来进行算法性能评估实验。MIT-B町也律失常数据库中的数据是麻省理工学院的BethIsrael医院也律失常实验室在7519巧4000Holter4819年到年间采集于多条数据。该数据库总共有个记录,每个记录时长30分钟且包含两个导联的同步数据。这48个记录取自47个个体35 东南大学煩±学位论文一2睹从23(其中201和202来自于同个体),包括年齡从32岁到89岁的5个男性和年22一岁到的岁的个女性,约60%的数据取自于住院病人。第组数据,100到124之间的23个数据从4000多条数据中随机选取;第二组数据,200到234之间的25个数据包含不常见但具有非常重要的临床现象的数据。1111111..i45!,11.00034.05060708091.10立.0.....巧间T<s)(b>巧巧巧■必电巧巧I:_:II<I?II5II00,10.20..0.50.60.70.80.91304时间T<s)图4-5仿真屯电波形-MITB出的数据采样率为360Hz112,,采样精度为位。每条记录包當个通道数据一-二般第通道采用的是校正的肢体I)1II导换(即MLI导联第通道采用校正的VI、VV、2或5导联数据,。对屯电数据进行压缩不需要转化为实际意义的也电数据反而需要的是将-I11位的二进制数据直接转为十进制的数据,因此,我们需要了解MITBH的数-B据格式。MITIH的数据存储采用212数据格式。212格式用3个字节(24位)来存储2个数,即每个数据占12位,最离位为符号位。212数据格式如表4.2所示,第I字节和第3字节分别存储信号1、信号2的低8位数据244。第字节的低位、高位分别存124储信号、信号的高位数据,最高位为符号位。-表4.2MITB田的212数据格式示意胃巧一字节第二字节第H字节 ̄[I信号1的低8位信号2的由4位信号I的高4位信号2的化g位 ̄ ̄^^I胃胃76s432l〇巧号l〇98巧号109876s432l〇|||||j||||III|j|||||也电数据的压缩的应用对象是从临床上的病人身上实时采集来的也电数据,不仅包括正常信号,,还包括大量也律失常信号这要求算法不仅要能很好的压缩正常屯电信号,还要能有效压缩也律失常信号。由此可见IT-BIH也律失常数据库采集于临床病人也,M电数据,非常适合用于评估也电压缩算法的性能。4-.3.3基于KSVD训综宇典的压缩实離果与分析-将則练字典与固定DCT字典作对比实验,进行分析,展示用KSVD算法训练的字典的优点-SVD。基于DCT字典的压缩感知算法流程图与基于K训练字典的压缩感知算--法流程类似,如图33,只是不需要用KSVD算法去训练字典,只需按照式(3.10)构造DCT字典即可。本小节将从三个方面进行分析;不同压缩比下的重构误差分析;不同巧声下的重构、误差分析;实际屯电信号基于压缩感知的压缩重构效果分析。乂 第四章也电信号的压缩盛知算法L不同压缩比下的重构误差分析压缩观测值M的个数不同,压缩比就不同,压缩重构后的误差也不同。这里将测-VD训练字典和固定DCT字典下不同压缩比的重构误差试在KS。===首先取W360,王720,护1440在仿真也电信号上加入SNR30dB随机离斯噪声构,-SVD算法公式造S个样本组成矩阵X,作为训综样本。按照K(3.2)进行字典训练,护=DCT-1〇|字典训练次数巧10,采用字典作为初始字典,得到KSVD训练字典。=360=f=X-仍令W,/;720,OMP算法迭代巧值l10、3,测试信号固定为仿真无噪)=屯、SNR30dB/及-电信号加上高斯随机噪声。当A依次减少,C依次増加时SVD,在K训练字典和固定DCT字典下,测得加噪仿真必电信号在不同伍缩比下的百分比均方误4*6差如图。如图可知,当C;?増加时,在DCT字典下的CS压缩重构误差快速増加,K-SVD训S压但在练字典下的C缩重构误差増加庭势缓慢,且都小于DCT字典情况。--仙VT::厂7;r:1*—^DC!:I::巧巧 ̄巧…………………一…………+Ksvi卻傑宇iIi4;i^iii!MMT/i-?■?0*3!、t中!t^}产-*…--——^——25i\i^j\jj\jJIiiiIMIk?"■…"?5乐?i4,;1户;}證……H-…-1IiiiiIiiiiiq1141五1822224之6283.>-压结比供图4-6不同压缩比下的百分比均方误差由本实验可得出两个结论;(1)压缩观测点数M越少,即C及越大,压缩感知实现-重构的难度越大,重构误差越大D算法训练的稀疏字典性能;(2)对如电信号用KSV优于固定DCT字典,更加符合也电信号特征。2、不同SNR下的重构误差分析沉……-…-r:;:T:一-:i;iI6c巧巧:I——iiHK4VD抑练车巧i:If—————r沒TiJJjJIj—……一—…---20j41jI片—--r5度15rt?iMiIIIIIji/SfjI ̄-."ytji1IJ\iiIi置100抑抑70饼幼40如进巧巧比SNR图4-7不同信噪比下的百分比均方误差巧 东南大学硕±学位论文、实际的屯电信号通常都包含各种各祥的噪声,主要有基线漂移、工频干扰、肌电干扰、运动干化等。本实验測试有噪声的必电信号在压缩感知算法下的性能。==护=V360,£720,1440也电信号上加入SNR30dBS取i,在仿真随机窝斯噪声构造=-。D32),,个样本组成矩阵X,作为训练样本按照KSV算法公式(.进行字典训练K10=-,K。字典則练次数W10,采巧DCT字典作为初始字典得到SVD训练字典==e=A=-360Xl〇,令W,A仁180,即C及2,王720,迭代闭值l(3给仿真也电信号加)-DCT字典上SNR依次降低的高斯噪声,在不同信,测试基于KSVD训练字典和固定-噪比下基于压缩感知理论的压缩重构情况,47测试结果如图。由图中可W看出,随着一,髙斯噪声的増加,压缩感知方法的重构误差将増加,且当噪声定大W后将严重影响-VD训练DCT字典压缩感知的重构性能于KS,;基字典的压缩重构的误差始终低于固定一‘。电稀疏字典的训综、这进步说明了用历史也电数据进行,得出的稀疏字典更加符合屯电信号持征。3、实际也电信号基于压缩感知的压缩重构效果分析MIT-。BIH本实验针对实际也电信号,进行基于压缩感知的也电信号压缩重构实验一口进行分段压缩的条数据是30分钟的也电信号,压编时可W采用滑动窗。每个窗口一的压缩参数和方法都样,这里仅展示前10s数捉的压缩重构效果。=AT360=720=144,取,1,50初始字典为DCT字典字典训练样本为标准也律失常,M--VD数据库ITB旧的100号数据的ML日导联也电信号的前1440s也电信号。按照KS==-算法公式(3.2)进行字典训练,K10,字典训练次数好10,得到实际也电信号的KSVD训练字典。A=====elX1-3俯,A^180,即(:^?2,王720,巧代停止条件设为迭代闽值03令成(),也、-输入信号为标准律失常数据库MITBIH的100号数据的MLII导联也电信号的前1化-也电信号,观测矩阵仍为随机高斯分布阵,分别在DCT字典和KSVD训练字典下进行4-8,,b)压缩重构实验,重构波形对比图如图图(a)为原始也电波形图图(为基于-DCT字典实现重构的波形图,图(C)为基于KSVD训练字典实现重构的波形图。在两-种字典下的重构误差对比图如图49,图(a)为基于DCT字典的重构误差曲线,其原P=-D训练字典的重构误差曲JU>2.16%(b)始信号与重构信号之间的,图为基于KSVZ=1线,其原始信号与重构信号之间的Wf).21%。(a)巧抬宙号1.C£tIII1tJ1Ittt-■fQ12345678910巧邮(s)(b)丢于DCT夺巧的重巧巧化。l>i" ̄IIIIIIItI^*1012347891056巧间T(8)(b)K-SVD里于^I重巧巧形I蘇字巧的。四!■-356789f012410时间T(a)图4-8实际也电信号在两种字典下的重构波形对比38 第四章也电信号的压缩巧知算法<a》圣于Dcr?:巧的S巧巧艺1.11111111j1A比占■-*'IIIII|I*1〇1234S910567时间T(8)b)i^svDil(S于jl谋字巧的重构谋gI1I0-■-550訓"片曲.1I片Mi晰置"O■'012467358910时间T<i)-图49实际亡、电信号在两种字典下的蛋构误差对比-SVD-对比在DCT字典和K训练字典下的压缩重构效果,可知在KSVD字典下重构波形更接近原始信号,原始信号与重构信号之间的误差更小。可见经过对历史也电数据的训练构成的也电稀疏字典更加符合也电信号特征,用于也电巧缩感知,可W取得更好的效果。-分析重构波形与原始波形可知,基于压缩感知在KSVD训练字典下的R波重构效11£)果较化但除去波外的波段重构效果不太理想.21%,比。此处户/^仅为1仿真也吨信号的压缩重构实验中的f及Z)要化,但实际的重构波形与原始波形差异却很明显,分析其原因,可能是本例实验数据的R波远离于P波和Q波,因此R波所占的能量比仿真信号的R波所占能量高得多,而民波处重构效果好,使民波处的残差能量占比小,使本例变小。4玉4基匹E字典的达巧实验结呆与分祈-W屯、B出中=律失常数据库MIT102号作为测试对象,给定模板形态分类闭值00.9,0=2护=3=5*。4402,,c按照流程图对1中的MLII导联也电数据进行压缩巧重构。对残差一阶差分王字K哈夫曼编码后r坐标所占进行,算上字典D,稀疏系数,每个也拍时间的字节数后,与原始数据每个数据占12字节相比.得到压缩比为3.巧,原始信号与重4-构信号之间的误差曲线如图10,1说化的数据都是没有误差的,可见本文的提出的压缩算法确实可W实现无损压缩。最后对本文提出的字典训综方法训练出的字典进行分析。按本文提出的字典训练方法对M-B出中IT102号的MLII导联也电数据进行字典原子提取和压缩后,训练出的字一D中N-典有个normalbeat(lPVC57,4411),个,个PACE,个FP总共1个原子。图分别是送四种QRS模板的波形图。可见本文提出的基于特征波匹配的字典原子设计方法确实能提取出屯、电信号中不同种类的特征波作为字典原子,使字典原子具有实际医用价值PVC是。例如常见也血管疾病中的室性早搏,由于采用本文的压缩方法,定位了每个也跳周期所用的特征波模板类型,知道了巧些也跳周期由PVC类原子构成,也就是定位了也电信号中发生室性早搏的位置。由此可见,本文提出的压缩方法在实现特征波提取和压缩的同时,也实现了屯、血管疾病的分类和定位,为也电图的自动诊断打好了基础。39 东南大学硕±学位论文111111111|0-60■-.60?.42‘-AD-看一■0^.-04-0.6-0.8II?IIIIii■1020040D说0800100012001400160018002000巧间T<s)图4>始信号与重构信号的误差曲线10原(a)normiybeM(b》prematurewntrioilarcontraction—004.511°A■"""I里A/A415〇0.的0.10.化0之〇0化0.10.巧0.2巧间T(s)巧间T<?)(c)pacedbeA(d)fustonofpacedandnormalbeAIf1|1""iIK。。I!J\.■.J-0!.■50O.OS0.10.巧0.200.姑0.10.巧0.2巧间T<s)时间T(s)4-n图102号数据训练出的字典原子波形图在本例训练出的字典中,有重复的PACE和FP模板.分析其原因有二:本文采用的模板匹配的也电特征波分类方法进行的特征波分类并不完全准确的:同种也血管疾病的变异波形类型并不是唯一的。因此,本文提出的基于持征波匹配的字典设计方法是有可能会提取重复的同类型特征模板的,但绝不是无意义的兀余,这些同类模板之间形态-SVD字典训上也是有差异的。且与K练算法相比,利用本文提出的算法针对实际也电,信号训练出的字典原子都是特征波模板中的典型模板,且能实现自适应更新更加符合医学应用要求。-WTBIH也48个样、MI律失常数据库中全部本的MLII导联屯电数据作为评估数据-465集,按照压缩流程图4进行压缩,将压缩结果与文献[中对也电信号直接进行基于]一阶差分的H字长哈夫曼编码的压缩结果进行比较,。根据对残差中的数值概率统计结果===23c5-o6io,I确定,给定0.9TB阳也律失常数据库中每个,,模板形态分类闽值尸按M-,压缩比结果如表43。表中实际特征波匹配数据占口位计算.哈夫曼编码压缩比是指除算上残差所占的字节外,还要加上字典D,稀疏系数r,各必拍的时间坐标所占的字节48-33数后的皮缩比。其中个数据的直接差分哈夫曼编码压缩比平均值为.1,模板匹配40 第四章必电信号的压缩感知算法3-1夫曼编码压缩比的平均值为.巧,实际模板匹配哈夫曼编码压缩比的平巧值为3.5。因此按照实际压缩比,本文基于特征波匹配字典的压缩比比文献[65压缩方法压缩比]的提髙了6.04%。显然,与直接差分字长哈夫曼编码压缩相比,本文提出的压缩方法可一政进歩减少必电信号的冗余,提高压缩比。表4.3基于特征波匹配字典的压缩实验结果C民CR ̄ ̄数据直接差分-轉征波匹巧-实际特征波匹数据品差分持征波匹馬-实际梯征妓[IIP序号哈夫曼编码哈夫曼编码配-哈夫曼编巧序号-哈夫曼编码哈夫曼编码匹巧-哈夫曼压缩比皮缩比压缩比压缩比压缩比编码压缩比1003.653.953.85203.874.化4.051...巧2023...1013诚3823343653571023.693.813.722032.812.962.871033.503.773.692053.744.174.05......巧10434235534420734736431053.143.483.巧2083.033.巧3.14......1063193273202093083403301072.532.%2.71203.463.803.6811083.巧3.443.%222.903083011......巧2..1093033563472132斯288..巧3.1113,353.613.巧2143073311123.巧3.953.化22.863..1517307..1133.253.423.巧2172843103.03...1143533703.612193453.843.751153.巧4.033.952203.543.923.831163.053.巧3.262213.293.623.巧1173.巧3.如3.巧2223.503.603.511182.903.093.032233.%3.783.671193.203.巧3.452283.283.443.%1213.954.214112303.301723.63.1223.143.623.巧2313.693.%3.881233.5..13803.742323773.933.861243.513.953.882巧2.933.393.巧2州3.103.203.122343.433.933.82虽然,在考虑字典D,稀疏系数r,各私拍的时间坐标所占的字节数后,实际的压缩比有所下降,但是,由于也电信号是准周期信号,因此,字典D中的原子个数并不会随、MIT-BIH屯着信号数据量加大而快速增加,而是相对保持稳定。同时由于律失常数据库中的每个样本仅包含了30分钟左右的数据,而也电监测都是长时间持续的,所若一是也电数据持续时间不断增加-,表中的实际特征波匹配哈夫曼编码压缩比将进步提-升,趋向于表中的特征波匹配哈夫曼编码压缩比。综上所述,针对也电信号,基于特征波匹配字典实现私电压缩时有四大优点:(1)^可处剔除过多的冗余原子,训练的字典原子具有生理医学意义2)D;(实现了字典的动、态増长,在出现屯律失常变异波形时,本文提出的算法可是及时捕捉到变异模板将其加入到字典中;(3)将模板匹配后的残差进巧了哈夫曼编码和传输,使有损压缩转换为无41 东南大学硕±学位论文、损压缩(4),在实现字典训练的同时,;结合屯电压缩与特征波识别、特征波分类与提取、可完成也电特征波分类,、屯律失常点定位和也电压缩为瓜电图自动诊断奠定基础。4.4本章小结-、VD训练字典的压缩感知实现方法本章首先介绍了屯电信号基于KS,阐述了在、K-SVD训练字典下的压缩过程,给出了流程图。然后介绍了屯电信号在本文设计的基于特征波匹配字典下的压缩算法的基本原理和实现方法,并给出了算法流程图。、在介绍了压缩原理后,在进行实验前对实验平台、屯电圧缩性能评价指标、实验的评估数据集作了简要的介绍。通过在不同压缩比下的重构误差分析和在不同SNR下的、-重构误差分析,得出了利用历史也电数据通过KSVD算法训练得到的稀疏字典比固定DCT字典更加符合也电信号特征,可W降低重构误差。在指出了用本文提出的算法训练,的字典基于压缩感知框架进行也电信号皮縮的实施难处后,提出了新的压缩思路并将压缩感知的有损压缩转变为了无损压缩、电。通过实验证明了基于特征波匹配法的进行屯M-、ITBIH屯信号的皮缩提出的字典原于具有实际的医学意义。最后给出了律失常数据1。库中全部样本的压缩结果.取得了3.5的无损平均压缩比42 第五章总结与展望第五章总结与展望5.1总结利用也电监护仪实时监测也脏病患的也电数据,有利于提早发现病情,降低发病死亡率,但这会产生大量的也电数据需要传输和存储。将也电数据压缩后再传输和存储将sPi有效降低系统传输功耗,节约存储空间hanian。Mamag等人首先将压缩感知理论应用于也电信号的压缩,但只是简单地应巧圧缩感知理论进行也电信号的压缩,并没有根据也电的特点进行深入优化压缩重构效果一。基于此,本文进步研究了基于压缩感知的也电信号皮缩算法。本文进行的主要工作和成果有:、、1分析了也电信号的自身特点,指出利用历史屯电数据进行字典训练,得到的字)典将更加符合也电信号特征-VD算法训练的字典比固定。用实验证明了用KSDCT字典有更好的稀疏表示性。2-SVD-分析了巧K算法训练的字典的不足:KSVD训练的字典原子多是重复的):且字典训练完成后就固定了,字典长时间固定不更新,不能有效反映也律失常情况。3屯、分析了常规也电信号分析与压缩的具体流程,指出应将电信号的压缩重构与)屯、电监护仪上初步的特征波识别,直、特征参数提取、自动诊断联合起来考虑接通过字典训练算法提取屯、电信号中的各特征波姐成稀疏字典。这个设计构思精简了系统流程,可有效利用系统资源。-4、SVD训练字典的压缩感知压缩实验)实现了也电信号基于K。与固定DCT字典屯、,相比,经过对历史必电数据的训练构成的也电稀疏字典更加符合电信号特征用于也电压缩感知,可W取得更好的压缩效果。5)结合2)、3)两点,提出了基于特征波匹配的稀疏字典设计方法,给出了基于特征波匹配的稀疏字典训练算法的目标函数及其算法的伪代码。分析了直接根据压缩感知流程使用本文提出的字典训练算法进行压缩的不便,阐述了本文的压缩方法框图,将信号经过特征波匹巧后的残差进行了编码并传输给信号重构-、端,将压缩感知的有损压缩转换为无损压缩。观测统计了MITBIH屯律失常数102据库中号样本用基于特征波匹配的稀疏字典训练算法训练出的字典原子,M-都是由私电特征波的典型模板沮成,具有实际的医学意义ITB^^也。对律失ML、常数据库中全部样本的II导联屯电数据进行了基于特征波匹配字典的压缩一1实验,取得的平均压缩比为3.5,比阶差分H字长哈夫曼编码方法提窩了6.04%。5.2展望虽然采用特征波匹配法训练的稀疏字典进行也电信号的压缩,取得了平均压缩比达一3.51的无损压缩,但在进行实际足缩时,进行了定的简化。本文仅为基于特征波匹配一的字典设计做出了初步地探索,为后续研究奠定了定基础。后续的研巧可料从W下几43 东南大学硕±学位论文个方面进行拓展;1)本文在进行字典训练时,仅考虑了屯吨信号中能量占比最大的QRS特征波,后续研究应该考虑将也电信号中所有特征波的典型模板加入到字典中,使也电信3号在本文训练字典下的稀疏表示的稀疏度为,即每个也动周期的必电波形由P波、QRS复合波、T波的典型模板组成,完成完整意义上的基于特征波匹配的字典训练算法。并在此基础上可W使用与同类压缩算法的比较实验。2i十、在算特征波匹配度时,应该充分考虑到屯电信号的准周巧性W及特征点定位)的不准确性,计算相似度系数时应将特征模板进巧微小位移,W寻找最佳匹配。3)基于特征波匹配字典进巧压缩实验村,采用的是滑动窗曰法进行QRS特征模板的提取,且窗口长度都固定为120ms。这并不能很好地适应所有人的必电数据,显然本文的特征模板提取方法有待改进。特征模板提取需要采用更精准的特征波定位算法。4本文采用的模板形态分类和提取方法是模板匹配法,模板匹配法虽然计算简单,)、但分类准确度仍然有待提島。因此,后续的研巧可W将可靠稳定的也电特征波分类方法应用于本算法中,提高算法也能。44 致谢=年研究生生活即将结束,在此,我对在此期间给予我指导、支持、关必与照顾的所有老师、校友、同学、朋友和家人表示最诚擎的谢意。首先我要感谢我的导师堵国楼老师,还有我们教研室的朱为老师和董志芳老师,他们都给我三年的学习生活给予了许多的指导和帮助。堵国操老师和朱为老师给我指明了课题研究方向和方法,并旦化我根据我个人的兴趣爱好选择课题,让我在研巧生期间学到了很多东西,并提离了实践能力,让。感谢董老师为我们组织的小组讨论会我和其他人能够相互学习。最后感谢几位老师为我的论文提出了许多宝贵意见。其次-我要感谢我的校友徐拥军老师和王丽艳老师,他们是我的校友,也是与我的教研室有合作的纳龙公司的工作人员。我研究生期间的最后两个项目都与他们有洽作,他们对我的项目研巧和本论文的完成都给予了很多帮助和指导,在此我再次对他们表示深深地感谢。还要感谢课题組的蒋富贵、彭鹏、李梓横、李婷同学对我完成本课题的帮助。与蒋一富贵的讨论交流使我收获颇丰,感谢彭鹏同学对本文阶H字长哈夫曼编码算法做出的贡献,感谢李子恒同学提供的硬件平台支持.感谢李婷同学为我修订论文。最后感谢所有在这H年期间帮助、关也、照顾过我的同口师兄师弟、朋友和家人,谢谢你们!45 参考文献参考文巧、I.[刘志学.《20B》血管疾病己成我国重大公共,胡洋中国也血管病报告在京发布也]卫生问题J20-[.中国医药导报,14,26:117117.]()口].T慎平,王应海,吴卫荣,邓玲黎,吕继东.便携式必电监护仪设计进展[J].生物医学-工程学杂志201403J08713.,,、3.林加锋张建华陈晓曙等.也脏性晕厥或巧死的若干高危也电表现J.也电学杂[],,,[]志04-:1951.,2的,098仰戴琼海,付长军,季向阳.压缩感知研巧町计算机学化2011,34(3):425434.[5].赵安,吴宝明,高上凯.也电数据压缩技术的新进展[J].国际生物医学工程杂志,20052-:8993.,().JalaleddineSMSHuchensCGSaan民DECGdami间,t,trttetal.tacoresson,ptechn-iuesaunifiedaroachJ.SiocalEnineeringIEEETrani[]me出gsactonson1990qpp,,374329-343);.(7.吴丽华.也电信号压缩方法的比较化中国科技信息,2009,(24).[][8].高力蠢.ECG信号的自适应无损压缩[D]:[硕±学位论文].西安:西安电子科技大学,2010.ECG赵志华.信号处理与数据压缩:.[外[巧[硕±学位论文]哈尔滨:哈尔滨工程大学,2006.10.杨胜天.生物医学信号无损/近无损压缩技术的研充[D.浙:[]]:[硕±学位论文]江浙江大学,2002.II].梅靈.ECG信号的数据压缩:硕±学位论文.:2003.[[D][]南京东南大学,12.AlamSSAAlamSetal,ADPCMbasedElectrocardioramcoderwith[],,,gthresholdingforrealtime化lemonitoringalicationsC.CommunicationsandSinalpp[]gPro-cessinICCSP2014InternationalConferenceon.fiEESOMinGlSO.g(),i[13].MuellerWC,arrhythmiadetectionprogramforanambulatoryECG卫-monitorJiomedSciInstruinent^voL14,1978.8185-[],pp"W.CoxJ氏No化FM,FozzardHA,巧化AZTEC,aPreprocessingProgramforReakTimeECGRhytiimAnalysis[J].BiomedicalEngineering,IEEETransactionson,-1%82:128129.,()5IMSh-A1.shiimaasainetal.ScanlPolnalAritioforDa[]j,,,ongygoppoxmantaCompressionofElectrocardiograms[J.BiomedicalEnineerinIEEETransactionson]gg,,1983ll723-72乂:,()一16.费敏.ECG压缩算法J.[],胡广书,周礼果种新的瓜电()数据实时[]中国医疗器-械杂志1991:2492%.,,口)17ADT屯、?李刚.L.[],冯静电数据压缩算法的快速实现阴中国生物医学工程学报,1993-5:236240.,()18.谭新?可.[],刘虔城,徐彬锋穿戴式多参数监护装置信号处理平台的设计与实现阴 ̄中国组织工程研究2012如:4化74871.,,()。刊?GangL,JingF,LingL巧al.FastrealizationoftheLADTECGdatacompre巧ionmennndMa-thodicineandioloazneIEEE199413:225.EgieeriginMeBgyi558.…,,,g[20].王春光,刘金江,孙即禅.基于稀疏分解的也电数据压缩算法[J].中国生物医学47 参考文献2008-1工程学报:1317.,,()21不详、T压[.王培康?屯.994],电图数据的分段DC缩方法阴数据采集与处理,1,4-(:273277.)22.尚宇..1[],梅翠松离散余弦变换在必电信号压缩中的应巧研究[J]微处理机203.,23G.毛云甲?基于小波变换和也搏间相关性的EC信号压缩[D:.[]][硕±学位论文]云南:云南大学,2013.[24.LuZZLLuZ,约al.Wavektcomressio打ofECGsinalsbthesetartitioninin],,pgypghierarchicaltreesalgorithmJ.BiomedicalEnineerin,IEEETransactionson,2000,[]gg477-:849856.()2-5.阁相国..西[,郑祟就,伍晓宇等基于整整小波变换的必电数据压缩算法[J]安]交通大学学报1999巧12.,,()[26].HiltonM,LHM,HiltonML约al.Wavel幻andwaveletpacketcompressionofelectrocardiogranis[J].BiomedicalE打gineering,lEE巨Transactionson,1997,44口):394-402.口7].李延妃吕文红,郭银景等?基于HMM的也电数据压缩方法研究阴.计算机工程与应用,2006,(25:183-184+2%.)28[]?林泉宋文强彭承琳?基于BP神经网络的也电数据压缩研究饥?医疗玉生装,,备-,20089:1921.,()9CanTear-Oimal窃nalRecoverFromRandomProions:Universal.desEJTaoNtect口],pgyjnraees?-Encod.rmaonornaonon:.igSttgiNfoHy,2006,巧540654巧^Ra]口0]?CandesEJ,Rombe巧J,TaoT.Robustuncertaintyprinciples:exactsignalreconstruction-fromhighlyincompletefrequencyinformalionJ.RanaononNformaonHor,6,:489[]y20052509.班-31.DonohoDLCompressedsensinJ.NformaonHoryRanaonon2006:.[]g[],,1289130631[斗焦李成,杨淑媛刘芳等.皮缩感知回顾与展望阴.电子学报201,,,7-1於(:16516.)3MHAtiD-!.amahanianKhaledNenzaetal.Comressedseilime口,,,nsngforreat]gpener-eficECGnonwegyientcompressioirelessbodsensornodsJ.Biomedicaly[]EnEEn2015824%-ineerin按圧Transactionso1:2466.g,,巧)[3叫.DixonA,RDAM,DixonAMR,etal.CompressedSensingSystemConsiderationsfbrECGandEMGWirelessBiosensorsJ.BiomedicalCircuitsandSystems,圧EE[]Transactionson,2012,6:156-166.ZhT-S.ZhilinanzPinJunMakeietal.Gomreinforpyg,yygg,g,p巧edSensgEn-nVergyEficientWirel的sTelemoitoringofNoninvasiveFetalECGiaBlockSarsepBayesianLeaming[J].BiomedicalEngineering,IEEETransactionson,2013,60口):300-309.屯、3:.电信号智能分析关键技术研究[D.:[化姚成][博古学位论文]吉林吉林大学,2012.7:,李振廷.也电信号检测及压缩算法的研巧[D]硕±学位论文].哈尔滨口:哈尔滨工][业大学,2008.38、、,王.电信号特征波检测及屯:[]春光基于稀疏分解的也电数据压缩[D][博古学位论文.长沙:国2010.]防科学技术大学,[39].ThakorNV,WebsterJG,TompkinsWJ.EstimationofQRSComplexPower48 参考文巧"SectraforDesignofa民SFilter?BiomedicaErnIEEETpQ口]lngineeig,ransactionson,-19841102706八.,(40.Tsaig义onohoDL.Extensionsofcom巧巧edsensinJ.SinalProcessin2006D,,[]pg[]gg^86:549571.41BDLuDBouP.oufounosifounosTlAlectureoncomressivesensinJ,IEEE[],,pg[]roc的snaazne4-SinalPigMgi200724:118121.g,,()42.C加化sEJEldarYCNllDtl.Comressedsensinwihcoherenandredundan[],,eede,eapgtttd.tt-ictionariesJAliedandComuaionalHarmonicAnalsis2011313:5973.f]pppy(),,43MSGZanZMaturuwme-uencrs.allath.chingsitsith巧fredictionakJ.IEEE[];gpqy[]T-ransactionsilrocessin:715onSgnaPg,19巧,41口3巧34.()44.TroJAGilbertAC?幻nalRecoverFromRandomMeasurementsVi泣[]pp,gyOhonamaon-rtlMatchinPursuitNforHorRanaonon20075346gog.y:巧。,,〇巧4666.[4y.ChenSS,DonohoDkSaundersMA,etahAtomicDecompositionbyBasisPursuJSm-IAMJoumaonScenit.litificCoutin1998l:;336L[]pg,,()[4y.FigueiredoM,Nowak民D,WrightSJ.GradientProjectionforSparseReconstruction:ApplicationtoCompressedSensingandOtherInverseProblemsJ.SelectedToicsin[]p圧EE-SinalProcessinJournalof200714586巧7gg:.,,,()47ansenooser-Mlak巧alInriPoinlfoLinr[.JBRCTT.teortethodoorea],,y,gyProgramming:Duality,SensitivityAnalysisandComputationalAspectsJ.[]OptimizationinPlanninandOerationofElectricPowersems13.gSt,99pyK--4imSKohKLusiMeal.AnIneriorPoinMehodforLarecale?t,ttttS,,gg[巧-圧距II民eularizedLeastSuaresJ.SelectedToicsin^nalProcessinurnalofgq,Jo,[]pgg20074-1606617.,();4Has*?EfronBtieTJohnstoneIetal.LeastanleieressionJ.Ann.Statist2002[刊,,,gg[],,2407-4532:1.()-ma0?GilbertAMuthukrishnanSStraussM.ImprovedTimeBoundsforNearOtil口],,pSparseFourierRepresentations饥.inProc.SPIEWaveletsXI,2004.口。.GilbertA,Strau巧M,TroppJ,巧al.Algorithmiclineardimensionreductionmthenormforsarsevectors?erintarXiv乂s/06080792006.p…p,口巧,GilbertA,Strau巧M,Vershynin民.Onesketchforall:fastalgorithmsforcompressedsensingJ,InProc.39thACMSm.TheorfComtin,2007.[]ypyopug53、电信号压缩算法阴.赵泽文??北京生物医学工程,方滨,王普基于稀疏分解的屯,[]20092848-4851.,:口叫.RauhutH,Schnass1C,VandcrghcynstP.Compres化dSensing抑dHcdundantD-ictionariesJ.InformationTheorIEEETransactionson2008545:22102219.[]y,,,()ausenBAearsecodnwtovercomeeassse:泣srae.Olsh巧WtDJ.Spiihanltbittt口3,gpgyemployedbyVI[J].VisionResearch,1997.Enan-wme1C民aoBKreutzDeFOCUSS口巧?glgadoK.Framedesignusingiththodof,,optimaldirections(MOD).Proc.NORSIG,1999.99.JFKreuz-De-MurratladoK.AnimrovedFOCUSSbasedlearninalorithmfor.y口可,gpggsolvingsparselinearinverseproblemsC.SinalsSstemsandComuters2001.[]g,yp,nferencer-圧EE-CoirFftAsilrnferen1:5Recodof化e孔tyihomaCoceon.20034731?,49 参考文献AronMEadMBrucAK-VDAnADehalrithmforsininof?lkstein?化eS:o口糾,,gggOvercompleteDictionari巧forSparseHep化化ntation,IE££lYansactionsonSignalfflP-rocessn200654ll:43114ig322.,,()59.M:,,胡广书数字信号处理理论、算法与实现(第三版).北京清华大学出版社[][]-2012;328333.60屯、:..郑臣.实时动态电图认知与自动诊断系统的研巧[D]硕±学位论文]北京:北[[]京邮电大学,2012.、D6博±学位论文:1.季虎.屯电信号自动分析关键技术研究:.长沙国防科学技[][][]术大学,2006.62.段会龙,冯靖讳,洪巧.基于模板匹配和特征识别相结合的也室期前收缩波形分类[]200212-102算法[J].航天医学与医学工程598.,,():一63.马杰.1998,姜向中种综合的也电数据压缩方法[化复旦学报:自然科学版,,[]-6:759764(),E-64lendiMFastSetectionwithanotimizedknowledebasedmethod:[,.d]gQRpgevauaonon11srECG.Fl201389乂73557.ltitandaddatabases.osOne,,()!;。65.张跃徐廷松.单片机也电数据无损压缩方法研究与实现J?医疗卫生装备[],[],-200829:2427.,[6句?IvayloIChristov.RealtimeelectrocardiogramQRSdetectionusingcombinedadaptivethre化old[J].Biomed.Eng.Online,2004,3:28(1):28.67.Cardiacmonitors,heartratemeters,andalarms.ANSI/AAMIEC13^002/(R)2007?[]巧巧.McSharryPE,ClifordGD,Tarassenkoetal.Adynamicalmodelfbrgenerating-lIEB2003sntheticeectrocardioramsinals.J.EETransiomedEn503289ygg[]g,,():294.6-.Phsionet.MITBIHArrhthmiaDatabaseDB/〇.htt://hsion巧.or/hsiobank[种yy[mppygpy-aabase2001-/dt/mitdb/1031.,〇MoodyGM民GTMIT-BMArriDJIE.Bark.heimactof化ehytfmiaatabase.EE,[7]p[]Enn-M-MedandBo204550June2001PM瓜:114462giil(3):(ay).(09).71oererAI>AmaraANasseaPhsobanksoooandsone.Gldbg!GlLtlihyitlkithyi:,,,,y^[]ppComponentsofaNewResearchResourceforComplexPhysiologicSignals[J].Crcuaonvo-101e2152000iltil..e220.,,pp50

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
大家都在看
近期热门
关闭