《基于环境振动响应的桥梁损伤识别》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
⑧申请网济大学工学硕士学位论文y‘78S6《蠢基于环境振动响应的桥梁损伤识别(国家自然科学慕金项目编号:50278067上海市科委辩磷玻关计划重点矮嚣编号:032112060)培养单位:一级学科:二级学科:研究生:指导教师:土木工程学洗土木工程褥粱与隧道王程黄莉萍张启伟教授二oo五年三胃 摘鼹璃要损伤识别技术是目前桥梁健康监测研究的关键问题之一。基于数值模型的损伤识臻硬究已经不笈双真正意义上瀵足绥梅,j、撰巍弱要求。考虑刭壤橥结稳铰多的不确定因索、复杂的工作环境以及出此而引起的结构动力特性的变异性,本文提出了抗噪能力强、环境适应性好、对中、小损伤敏感的损伤识别新方法,并遵过嚣跨遣续梁试验验谣该方法的有效性。本文综述了墓予结鞫响应静损伤谈潮鼢理论磺究躐、欹和发展。放院较线髂、非线性和其它时频方法用于非平稳信号分析的优缺点入手,引出了适用于环境振动响应分析的Hilbert.Huang变换,详细地阐述了Hilbert.Huang变换的基本组成——经验模式分薅秘蓑尔{鑫特交接,劳攒窭Hilbert-Huang变换在≤}乎稳岩号处瑷方面的优越性和有待改进的方面;同时介绍了本文采用的模式识别系统的主鼷组成,如何将模式识别系统引入到土术工程领域,并使其适用于结构的健康浆测也是本文麴一项重要任务。本文摄蹬了一释桥粱缭输损伤识剐静新方法。蕾先对采集静环境振动下的缀构响应进行预处理,由于缝自适应性强的经验模式分解任何复杂的时间序列都可分成小尺度的本征模式函数,从而可利用Hilbert.Huang变换提取各个尺度上稳基予缝撬瓣瓣特毽麓瓣辩旗篷特薤,继嚣形戏瓣辩疆篷特薤群本空藏,对该样本空间设定阈值进行特锻样本的选取以剔除噪声等对结构损伤识别造成影响的样本,最艏对选择后的特征样本进行点到集合的标准欧氏距离度量,实现结构豹二类模式识别飙面诊凝结橡损伤是落发生。同时,零文在损伤位置豹辕宠方面氇敲了籀成麴探讨,秘爝与瓣爵幅缓对应的瞬时频率酶特征释本修正均餐标出了损伤位置。在对试验锯原有的有机玻璃模型梁段改进的基础上,采用试验方法对本文提窭戆方法遂行检验。在试骏室繇凌条锌下,弱霜小车激瓣两踌连续模型粱莰采集梁段结构的加速度响应,进行数据处理验证了本文方法的有效性。同时,讨论了特征选择中不同闯值选取、传感器数量大小和损伤处是否安鼹传感器对该方法弱影哟。搔出整套损谯识澍方法基鸯遥弱毪,霹逶蠲予其它鹣时闽序列损伤特征参数。缀后提出了遮~研究领域的有待进一步开展的研究谦邀。 摘要关键词:键康蓝灏,Hilbert—Huang交换,特征祥本选择,损伤识翔Il ABSTRACTA88TRACTDamagedetectiontechniqueisoneofthekeyissuesforbridgehealthmonitoringresearchrecently.DamageidentificationresearchbasedonhumericatmodelCall’tsatisfylightNructuredamagedetectionfromactualmeaning.Consideringlotsofuncertainfactors,complicatedworkingconditionandnoveltyofdynamiccharacteristicofbridgestructural,thisthesispresentsastructuralnoveltydiagnosisapproach+whichispromisingforpracticalapplication诚氇strongabitityofresistingnoise,self-adaptabilitytoambientvibration,sensitivitytolightdmnages,asillustratedintheeasestudyona2-spancontinuousorganicglassbridge。Thecurrentandfuturetheoreticalresearchinthefieldofdamageidentificationisintroduced。Ontheonehand,Hilbert·Huangtransfomlwhichisapplicablefordealingambientvibrationresponsehasbeenintroducedbylistingouttheadvantagesanddisadvantagesoflineartime—frequencyanalysismethod,nonlineartime—frequencyanalysismethodandothertime-frequencyanalysismethod.Thenthedetailsof‘EmpiricalModeDecomposition’methodandHilbe娃Transformfollowed.ThedeficiencyofHitbe菇-Huangtransforminanalyzingnon-stationarydataispointedoataswellasitssuperiority.Ontheotherhand,theessenceofthepatternrecognitionsystemispresented.Howtointroducethepaaemrecognitionsystemintothecivilengineeringfieldandhowtomakeitaccommodatetostructurehealthmonitoringarealsothetaskofresearch.Anewmethodofbridgedamageidentificationisproposed。FirstlNthedatasequencecollectedfromthestructureunderambientvibrationconditionshouldbepretreated.Secondly,becausewithself-adaptiveEmpiricalModeDecompositionmethod,anycomplicateddatasetcanbedecomposedintoafiniteandoftensmallnumberof‘IntrinsicModeFunction’.Thedamagesensitivefeatureindexofsm.1eturalcorrelatedwithinstantaneouscharacteristicoa鼗beextractedindifferentscale-spacebasedonHilbert*Huangtransform.Theninstantaneousamplitudefeaturesamplescomeintobeingouthethirdstep.Fourthly,thesamplesinfectedbyIn ABSTRACTsuchasnoisesaretobedeletedsincetherealizationoffeatureselectingthroughdealingtheinstantaneousamplitudesampleswithconstantthreshold。Fiflhl弘bridgedamageidentificationbasedontwo—patternrecognitionmethodcanbewelldonerefetingtostandardEuclideandistancemeasurementbetweenapointwithpointsets.Inaddition,theproblemofdamagelocmionhasbeentouchedupon.Withstatisticalanalysisoftheinstantaneousfrequencytrimmeanwhichiscorrespondingtoinstantaneousamplitude,thedamagecanbelocated.Atestingmethodbasedonmodifiedlaboratoryorganicglassbeamshadbeenadoptedtocheckupthenewmethod.Thefeasibilityandeffeetivenessofthenewproposedmethodareinvestigatedviadealingwiththemnbientvibrationresponsemeasuredfromthedynamiccasestudyona2-spancontinuousorganicglassbridgeinlaboratory.Inaddition.theinfluenceofthethresholdduringthecoLn'seoffeatureselecting,thenumberofsensorstotheproposedmethodisanalyzedinthisprogramaswellasthenecessaryofsensoratthedamagelocation.ItispointedoutthatthenewmethodiSpracticableforothersubstitutabletimesequencedamagefeatureindexinthcstructuraldamageidentification.Finally,directionsforfurtherresearchOnbridgedamagedetectionarepointedoutbasedonambientvibration.Keywords:healthmonitoring,Hilbert—Huangtransform,featureindexsamplesselecting,damageidentification 同济大学学位论文原创性声明本人郑重声明:所坚交的学位论文,是本人在导师的指导下,进行囊暑究工作所取褥酶威袋,撰写戏漕±,硬士学位论文竺蓬量蔓蕉燕越响应的撅攘塑鱼迟别==。除文中已经注明引用的内容外,本学位论文戆褒究成暴不包含任凭匏人剞豫憋、基公开发表或者没有公开发表的作品的内容。对本论文所涉及的研究工作做出贡献的其他个人和集髂,均毫在文孛以饔确方式标翻。本学位论文簇剖性声明煞法律责任由本人承担。签名:截臻1降z弼?年冬月菇曩 学位论文版权使用授权书本入完全了解嗣济大学关于收集、保存、使用学位论文鲍规定,阐意如下各项内容:按照学校要求提交学位论文的印蒯本和电子敝本;学校有权保存学位论文的印刷本和电子版,并采用影印、缩印、扫接、数字化或其它手段保存论文;学校有权提供目录检索以及提供本学位论文全文或者部分盼阅览服务:学校有投按有关蕊定向国家有关部门或者机构送交论文的复印件和电子版;猩不以赢利为目的的前提下,学校可以适当复劁论文的部分或全部内容用于学术活动。学位论文作者签名:蕊寡,3辫j器$年弓囊j吾掰经指导教师同意,本学位论文属于保密,在j年解密后适用本授权书。指导教师签名:糊学位论文作者签名:巨茹1’烨v。书年≥月以目j咎轻拿年弓男荔曩 第1章绪论1,1桥梁结构损伤识别l。l。{桥梁缝褥损伤谖剃的意义一方面,随着社会的进步和经济。。。体化的加速推进,建设跨江跨海桥梁工程尤为追切;另~方面,综台国力的增强粒主本工程技术的发展又绘大跨凄桥梁鹣修建提供了飘遥秘撬战。疆溯蓦薷有不少已建、在建或都将建设翁大跨度轿梁项目,如同本的明石海峡大桥(悬索桥,主跨1991m)、中国的江苏苏通大桥(斜拉桥,主跨1088m)、香港的昂船洲大桥(斜拉桥,主跨1018m)等。人们列大鍪}重要援粱夔安全经、瓣久缝与燕露菠瑗功能醚海关注帮重视。虽然合理保:r的设计怒结构安全的搬本保证,但魑限于当前对大型复杂结构的认识程度,许多不可确定预知因素,如超期服役、疲劳、环境腐蚀、材料老化、构件缺陷、越风、地震、爆炸等危害瞧搴{孛,人们并不都能进行有效控制或预测,惩这些生命线工程缡褥又对国民经济、社会稳定帮人民的叠三命财产产生蘩大影噙,一鹰大桥的塌毁将影响到很大范围地区(一个或几个省)。为了确保设计的使用安全性和耐久性达到预期的标准,时时了解其“健康”状态,并及时对存在的损饶徽篷诊鞭势疆警是蒌露耋要豹。耩鬃毽瘫夔测系绫被陵续褒瘸予大壅撬粱中。到目前为止,已有不少大跨度桥梁安装了不同规模的桥梁健康赡测系统“。“1,如:明石海峡大桥、南备赞濑户大桥、大贝尔特东悬索桥、青马大桥、虎门大褥、徐溱大援等。强前援粱在线的毽痰媳测系统中输趣熬结构晌嫩时程是在环境随机振动下采集酶,不嚣激励设备麓不影响正露逶营,环境振动试验是丈登土木结构最实际可行的振动试验方法,遂就给基于环境振动响应的损伤识别方法研究带来了广阔的应用前景。基予缕携翡应懿诊鼗毅零已经在辊辍矮城褥妥了羧必深入兹疲瑗,芽取褥了一定的成果,为土木工程的结构损伤谚{别发展奠定了基础。基于结构的损伤识别涉及结构分折、结构动力学、信号处理、统计学及模式识别等领域,应用研究的发展氇襁戏的促进了这些领域砑究浆逶步。同融,缝梅毽康黢测已经弱裁成为国内外桥梁学术眷和工程界的研究热点,被多个领域国际学术会议嬲为中心议题之一。结构损伤识别也已成为动力学的六大挑战性问题之一o。。因此,基 第1章绪论予环壤振动响应静桥梁损伤识潮研究非常有意义。{.{,2基卡结构鹅瘟戆攫毪谖裂基本暴蕊Housner等定义结构健康监测为:结合无损感测与结构系统特性分析(包括结构响应),侦测结构的交化,议别结构损伤与性能邋优“1。结构损伤即结构木孝料或几何性熊包括边界条件的有懑或无意的改变,其反过来将影响结构现在以及将来的行为和能力。1。臻蕊发生豹形式多耱多群,逶拳穰=l委①累计损伤(如:疲劳、徐变、腐蚀等)<蛩突然损伤侄可碳测(如;蠢蛋攘桥墩等>(D突然损伤且不可预测(如:地震、台风导致的土木结构破坏等)一般情况下,土木工程结构突然损伤(可预测帮不可预测)发生的时间和部位都能被较准确地把握,从恧进霉亍梭测修复:丽绣掏的累计损伤不易被察觉,即通常人们都不知道损伤发生在哪里,也不知道什么时候损伤会影响到结构的正露镬窝。瑞典人Rytter-于1993年提出将损伤识别分为以下四个等级“3:④翔剐损伤是否发生;②确定损伤位置;(D评价损伤程度;(D评佳预测缝麴东常使翅豹剩余毒会。国际学者普遍认同这种分类方式,并认为结构损伤诊断的三个基本问题其难易程度蒎次为:捩伤与否懿判繇、按侮懿定位帮豢伤装窥量。瓣藩对损经识裂的研究基本还处于第一个阶段,即损伤发生与否的诊断。基于结构响陂的损伤识别基本恩想是损伤将改变结构的刚度、质量或耗能能力,进覆弓|起结构动力特性或赡痤臻号的改变,遵过提取结麴嘲应时闽序列媳参数信息或其衍生信息,并比对结构无损状态下的相应倍息,实现结构的损伤渗黼。1.2基予结构晌应的损伤谖别方法渣予主本工糕结稳稳特殊往,魏绪穆型式多棒、自骚役周翘长、影稳缝鞫妻羔笈2 第1豪绪论黝因素复杂多变等,对其邀行监测比一般豹机械系统嶷杂得多。桥梁结构尺寸天,质量重,羧动东平较低,结构酶麓力螭应极容荔受弱蒡|:凌因素戮及:{}篓鞫构件等的影响。损伤诊断方法的局限性归结为下”。:①桥梁结构赘余度大,而损伤一一般为局部现象;②主本绥季盘较多熬不麓定因素、复杂熬工彳乍环缓以爱大型槛罄致结构行为非线性时变,丽基于模型的诊断方法通常都将结构俸为线性时不变系统考虑:(D桥梁结构的响应与环境和运营祭件有关,采鬃信号不可避免具有离散性。目前方法对结构损伤的谈别灵敏度过低,与早期发现损伤这一目标差距较天;④由于实际结构可能的损伤模式及各种损伤模式下的结构响应通常难以预测,并且响应的数据样本匮乏,因而逡用模式识别方法时很难通过监督式学习进行凌态攘式分类。事实上,褥粱蘩梅键瘫整测鹣痘翔中,只§g依靠获褥的避续监测信号运用非监督式模式识剐方法进行结构状态的分类与诊断;⑤桥梁结构体量大,相比之下,监测测点分布稀少,测量信息不足。另外,测量信号噪声通常使得损俊识别更加困难,有时噪声水平甚至与损伤引起的结稳响应交纯稠当。由于影响结构振动的因索较多,目前尚无统⋯和有效的对结构损伤发生与甭、损伤位鬣和损伤程度’卜分敏感的特征。一般情7兄下,测量的数据相对结构懿塞交度数豁楚投乡粒,对大羹工程缡魏攒巍浚爨,不完全测量受楚一令营遮存在的问题。由于测量点的噪声影响和浏量点的密度分布的局限憾,很难得到司+靠的测量方法和高振动模态的识别方法,但充分的商振动模念信息是识别损伤位置的关键。为撂到测鬃和分析之间的一一对应关系,必须应用模态缩减或模态扩震蔽零,这些技术黪应矮蹩菇丧失精度力霞徐酶,瓢雨使寝多基j二完全测量的损伤识别方法的效果大为降低。湛于结构时獠响应的研究溅免了振型、频率的不完全测量性,因此被越来越多地应用到特,饺提取中。臻镑谀别中耱{垂参数戆提取萋霉关键。泰寒磷突发蓑戆特{歪参数澎具有豹特征。1:①对结构低频B瞄瘦敏感;(多抗噪声能力强;(葑适于“在线”环境随桃激励资料的处理;(D直接间接参数的翁溯性;(D对中、小损伤邋应度好;(国对损伤单调。嗣射还可以搬握特定桥型,进行不同损伤敏感特征的提取硬究。基于环境振动时程响应的特征提取骚用蜀狠多处疆方法,蟊经舞静谱分季厅方法、传统的时域波形法、小波变换、时域分布和高阶统计量分析法、统计能量3 第1章绪论分析法、非线性动力分析法,具体分为”1:(D频域法,包括傅立叶谱分析(功率谱、细化谱等)、多谱分析(信号高次矩的傅立叶变化)、倒谱分析等”。【”“”I。(D时域法,包括基于脉动输入的方法(RD函数、ARMA、Kalman滤波等)、基于脉冲反应的方法(IbrahimTimeDommn(ITD)⋯1、EigenrealizationAlgorithms(ERA)、随机子空间法等)“⋯。(D时频法,包括包括线形时频分析方法(短时傅立叶变换、Gabor展开、小波变换“”。“”等)、非线性时频分析方法(维格纳分布、Cohen分布等)以及其它时频分布(Cohen.Posch类正值分布、高阶时频分布、Hilbert—Huang变换等)。时频法将时域和频域组合在一个域内,是较好的信号处理方法。小波变换把信号分解成多重子信号可以更好显示时序信号,损伤处子信号呈显著变化,其中某些信号对损伤灵敏度很高,但小波基难于选取。维格纳分布是一科吲}常有效的手段,它是时间一频率两维联合函数,具有很高的分辨率、能量分布集中性和跟踪瞬时频率的特性,通过比较损伤前后系统的时频特性,可以获得有价值的信息,但是必须解决互谱项的干扰问题“⋯。Hilbert-Huang变换””是上世纪术提出的一种新的针对非平稳非线性信号的时频分析方法。本文采用了此方法,详细的内容在后边的章节中予以深入讨论。结构损伤识别是桥梁健康监测系统的核心部分。目前用于结构特征分离的分类方法主要有监督式学习方法与非监督式学习方法“7。“““⋯。应用较,11泛的监督式学习方法有神经网络法、遗传算法、支持向量机等,非监督学习方法有神经网络法、聚类方法等。模式识别的方法多种多样,主要包括:①聚类方法聚类方法对具有相同统计行为的多个因素进行归类分析,将功能上可能相似或有关联的因素归为一类。使用聚类方法可以评定内部数据的相似性或在模式识别研究的探索阶段想从一堆没有分组的信息中提取有意义的数据。⋯。②人工神经网络法不同的单元计算特性(神经元类型)、连接方式(网络结构)和连接强度调节规律(学习算法)形成了不同的神经网络结构。神经网络在自学习、自组织、自联想及容错等方面都表现出强劲优势,可以实现监督和非监督学习条件下的4 第1章绪论分类和回归工彳乍,能解决缎多非线性问题,并在很多:I:程中取得了成功”“’”⋯。}率经两络爱予维穆损伤谈弱有三个霞点;绸练遥静棒经隧终麓够存德有关过程的知识,根据对象的正常历史数据训练网络,然后将此信息与当前测量数据进行比较,以确定损伤;神缀网络具有在噪声环境中得出正确结论的能力;神经溺终具有势辨擐困及损绕类塑瓣囊力。释港瑾工大学以香港鲍毒马、汀丸和汲水门3座大跨轿为背景运用人工神经网绦方法对大跨桥的健康和安全状态的监测方法进行了探讨”⋯。④统计分褥将不确定的因素与确定性的损伤谈剐方法相结合,并应用于实簖的复杂一l:稷结构中是当前隧待解决的课题。“’“⋯。损伤可以看作是附加激励,作用于结构上褥;|起结构嚷艇的改变,然两在环境条件下由于结构受噪声影响缎大,结构的小损伤弓l起绥梅响应信号黪交证可麓被鞣声信号掰淹没。Trendafilove。1等人稳出了利用统计模式识别的方法进行损伤识别,根据损伤概率的大小来确定损伤的发生。1997年HoonShon嗣Law“”等人根据贝叶斯估计原理,提如了对结构刚瘦交纯遂幸亍矮诗送行摈揍裣溯疆报。HoonSohn”“奁2000年提出了剩羯统诗过程控制进行连续监测损伤的方法。支持向疑机是在研究有限样本情况下的统计模式识别和更广泛的机器学=习问题上建立的新的理论樵架,在解决小样本、非线憋及高维模式议剃方蔼表现出了优越性,是疑前研究蚋热点。模式谖剐方法还有旬法方法,我方法考虑要谈剩对象静各部分之闻酌联系,对模式进行络构描述和分析,该模式识别方法较多地应用于指纹等图像识别中,当前在结构损伤识别中应用还较少。同时,模糊模式识别也被借撩到桥梁损伤渗蔹领域汹】。实嚣应震数戆势是夔上述豹一耱或耍秘方法结台趁寒,互逯畜无,阱实现尽量缩小主观识别和客观实际的纛异而进行模式识别。{.3本文的课题鸷景秘主要研究内蜜本文的课越来自国家融然科学基会项隧“斜拉桥健康监测中的损伤敏感参数特征与损伤识别研究”(课题编号:50278067)和上海市科委科研攻关计划重点项耳“既寄棱粱结构检测译定与加固瓣关键技术”(谍题编号:032112060)。镑对斜拉耩缝构健康蓝溺帮损伤谖嗣阏簇,先进行损伤敏感的褥髹生成、选择和提取,获得适用于结构损伤诊断的最优特征或最优特征组合,蒋运用模式识5 第1章绪论别方法进行结构分类器设计完成损伤的识别,为完善健康监测与安全评估技术建立理论基础。研究项目中的两个关键问题是损伤敏感参数特征提取与损伤识别⋯。针对项目中的两个关键问题,本文在基于环境振动响应的损伤识别理论和方法进行了研究。在特征提取方面,引入一种适合于非平稳非线性的时频分析方法一Hilben.Huang变换处理环境振动信号,该方法由经验模式分解方法和lHilbert变换两部分组成,可以在多尺度上显示子信号在时域和频域内的信息。在模式识别方面,借鉴信控领域的模式识别技术,力求实现损伤与健康特征模式的分类。同时由于结构响应信号具有随机性,对信号的处理必须在统计意义上进行。本文的主要工作包括:①寻找对桥梁结构的损伤敏感特征参数;(D形成特征样本并进行特征样本选择;③就已提取的桥梁结构损伤敏感特征进行桥梁损伤与否的诊断和损伤位置的标定。④设计试验室模型试验验证所提取出的敏感特征和损伤识别方法的有效性和实用性。参考文献[1]张启伟,袁万诚,范立础.火型桥梁结构健康监测的研究现状与发展[J].同济大学学报,1997,25(增刊):76~81[2]张启伟.大型桥梁健康监测概念与监测系统设计[J].同济大学学报,2001,29(1):65468[3]FarrarCR,eta1.Grandchallengesforstructuraldynamics[R].18“IMAC,Feb7-10,2000,SanAntonio,Texas.[4]GWHousner,LABergman,TKCaughey,eta1.Structrualcontrol:past,present.andfuture[J].JournalofEngineeringMechanics,1997,123(9):897~971【5jF.G.Yuan.StructuralHealthMonitoring.NorthCarolinaStateUniversityRaleigh.NC27695[6]Rytter,A.VibrmionBasedInspectionofCivilEngineering[D].Denmark:UniversityofAalborg,1993[7]韩大建,谢峻.大跨度桥梁健康监测技术的近期研究进展[J].桥梁建设,2002,6:69~73[8]韩大建,王文东.基于振动的结构损伤识别方法的近期研究进展[J].华南理r大学学报,2003,31(1):92~966 第1章绪埝[9]DebraGeta1.Identifyingdamagesensitivefeaturesusingnonlineartime-seriesandbispectratanalysis[A],Proceedingsof18”1IMAC[C],2000,1796~1802[10]SammanMVibrationtestingfornondemruetiveevaluationofbri媳es田.JournalofStructuralEngineering,1994,120(1):269-、300【11]lbrahimSR.Randomdecrementtechniqueformodalidentificationofstructures[J].JournalofSpacecraftandRockets,1977。14(11):696~700}i2】Yij鞋。YunCB.CompatativestudyonroodalidentificationmethodsusingOHtout-onlyinformation[J】.StructuralEngineeringandmechanics,2004,17:445~466[13]GurleyK,KareemA,applicationofWavelettransformsinearthquake,windandoceanengineering[J]。JournalofStructuralEngineering,1999,21:149~167Il《』HOHZeta1.Wavelet-basedapproachforstructuraldamagedetection[J]。JournalofMechanicalEngineering,2000,126(7):677~683[15]张贤达.现代信号处理[M],北京:清华大学出版社。1995116jHuangN.E.,ZhengS.,Long,etal.TheempiricalmodedecompositionandtheHill∞rtspectrumfofnonlinearandnonstationarytimesedesanalysis17181920lc】.Proc.R.Soc.London1998,A(454):903~995温熙森,胡茑庆,邱静.模式识别与状淼监控[M].长沙:国防科技大学出版社,1997吴逸飞译.模式识别一原理、方法及麻闷!j[M].北京:清华大学出版社,2003迭肇祺,张学’1:嚣.模式谈捌[醚]。能京:瀵华大学出舨社,2001GodinN,HuguetS,Gae撞nerR.IntegrationoftheKohonen'sself-organisingmapandk-meansalgorithmforthesegmentationoftheAEdatacollectedduringtensiletestsoncross-plycomposites[J】.NDTandEInternational.2005,38“):299~309【21]YamLH,YanYJ,JiangJS.Vibration-baseddamagedetectionforcompositestructuresusingwavelettransformandneuralnetworkidentification}≯1.CompositeStructures+2003,60(4):403~412[22】姜缁飞等.大跨度悬索桥损伤定位的自适应概率神经网络研究[J].土木工稀学报,2003,36(8):74~78[233姜鲻飞.藏+F耱经潮络熊皱秘优纯与换惦捡溯[翻。I£索:科学出敝社,2002124』SohnH,LawK}{.ABayesianprobabilisticapproachforstructuraldamagedetection【J].EarthquakeEngineeringandStructtwalDynamic,1997,26:1259~1281[25]TshilidziMarwala.Probabilisticfaultidentificationusingvibrationdataandneuralnetworks{J】.MechanicalSystemsandSignalProcessing,2001,15(6):1109~l12812翻Trendafiloveleta1.Twostatisticalpatternrecognitionmethodsfordamagelocalization.Proceedingsof17“‘IMAC。1999.1380~1386[27]HoonSohneta1.Abaysianprobabilitstieapproachforstructuraldamagedetection.EarthquakeEngineeringandStructuralDynamics,1997,1259~t281}28】HoonSohnetal。Continuousstructuralhealthmonitoringusingstatisticalprocesscontr01.Proceedingsofl8“‘IMAC.2000。660~667[29】ZhaoZ,ChenC.Concretebridgedeteriorationdiagnosisusingfuzzyinferencesystem.7 第1章绪论AdvancesinEngineeringSoftware,2001,317“-325f30】张痘接。耩粱健凄夔溅中瓣臻馥i特薤提取与买霞诊瑟闭.溺婆丈学学缀。2003,31(3)2弱~2628 第2章环境振动信号分析一般环境振动信号存在以下两个问题:①信号是非平稳的;(D信号是非线性的。首先须对信号非平稳相关的概念有所了解。平稳过程的特点是产生随机现象的主要因素不随时间推移而改变,因此随机过程的统计特性也不随时削推移而改变。严平稳过程是指相隔时间f的两个时刻,与,+f处随机过程所处的状态具有相同的概率分布,即F(x】,z2,⋯,X。;^,t2,⋯,f。)=F(xl,x2∥.,x。;t+f,‘+f,⋯,t。+r)(2.1)其中F(+)代表概率分布。宽平稳过程则只要求随机过程x(f)满足E[x(t1)]-E[x{tl+r)】E[x(t1)x(t2)]_Nx(t】+r)x(t2+f)】(2.2)(2.3)其中日+】是数学期望值。2.2和2.3两式分别表示单矩和联合矩的平移不变性。凡是不具有以上特性的随机过程即为非平稳过程。通常在对事物的认知过程中,我们总假定研究对象是一个线性系统,满足线性叠加规律,然而客观世界里一切运动及规律都是非线性的,线性系统只是在特定状态的近似结果。随机过程的时间序列非线性更是复杂,故一般将信号的非线性分析简化为二次的线性分析,如后面介绍对信号的处理方法——二次时频分析方法。频域法仅可以确定哪些频率存在,时域法可以确定时间信息,时频法可以确定某一时刻频率成分的分布情况。所以对于环境振动这样的非平稳非线性的信号处理一般采用时频分析法。2.1信号的时频分析方法乜钔一维时域信号z0)或频域信号x(国)映射到时间一频率平面上得到j维信号PAt,05')。由此,信号的瞬时能量和功率谱可以表示为㈨12=£只(f,co)dm(2.4)9 第2章环境振动信号分析Ⅳ(珊)I2=£只(f,o)at(2.5)信号在时域内f∈(^,f:),国e(co。,∞:)的能量成分表示为』2r2eat,∞)dtd&,,理想的eat,m)表示信号在时间频率点(f,甜)处的能量密度。对非平稳信号的时频方法分为线形时频和非线性时频两种。2.1.1线性时频分析方法线性时频变换从傅氏谱演化得来,其变换为线性的。若时域信号及其频谱是线性的,如x(t)=alxl(t)+a2x2(f)(2.6)X(co)=alXl(m)+(12Ⅳ2(∞)(2.7)其中d.与日:为常数,则时频表示只.(f,∞)与只,(f,oJ)是线性变换的,有只(f,珊)=a1只.(f,co)+d2只,(f,∞)(2.8)即所有的线性时频表示都满足线性原理。信号的线性时频表示主要有短时傅立叶变换、和小波变换两种形式。】)短时傅立叶变换短时傅立叶变换是采用加窗技术将信号在时间上分成许多段,然后对每--d,段进行傅氏变换,得到对应于不同时刻的信号频谱。短时傅立叶变换可以理解为信号x(f)在分析时间t附近的局部谱。对于信号x(f),可以表示为STFT(t,co)=bo谛’p—t)eⅢ7dr(2.9)其中^(,)为窗函数;“”表示取复共轭。在时间t的STFT是被加上窗函数后信号x(f)的谱,其所有的以t为中心的局部窗间隔内的信号特性都会在时间t的STFT内显示出来。高分辨率要求窗函数较短,同时要求带宽较窄,然而这两者矛盾,设出和Ⅳ分别表示STFT的时间分辨率和频率分辨率,两者应满足Heisenberg不确定性原理△r·6f≥二(2.10)4石这就阻碍了既有任意小的时间间隔又有任意小的带宽的窗函数的存在,也就使1n 第2章环境振动信号分析得STFT时频表示无法在时间和频率两个方向同时获得较高的分辨力。2)小波变换小波变换是80年代中后期逐渐发展起来的一种数学分析方法。该方法克服了传统的傅立叶分析和卷积滤波分析的局限性。由于它在时域和频域的局部化和可变时频窗的特点使得用它分析非稳态信号具有比传统的傅立叶分析更显著的优点。所谓小波,就是指具有零均值,在时域和频域内能量局部化的函数,其波形表现为两端衰减为零的小的波形。设小波函数∥(f)∈L2(C),r(C)为平方可积的复函数空间,则y(r)必须满足:c,:£嵴<。泣⋯其中≯(国)是妒(r)的傅氏变换。上式为容许条件,蕴含着步(o)=0,即函数具有零均值。基小波经尺度伸缩和时间平移得到一个予波簇,其形式为‰。(f):{P仁鱼)(2.12)其中a为尺度参数,b为时间位置参数。a用来调整子波覆盖的频率范围,b用来调整子波的时域位置,系数—;用来实现子波能量的归一化。4a对基小波函数∥(f),如果其时窗宽度为出,经傅氏变换后谱≯(∞)的频窗宽度为Aco;那么对于9t(t/a),其时窗宽度为aAt,经傅氏交换后谱眵(Ⅱ∞)的频窗宽度将为Aco/a。因此,小波变换对低频信号(此时a相对较大)在频域里有很好的分辨率,而对高频信号(此时a相对较小)在时域里也有好的分辨率。如果变动a和b,形成一簇小波函数,然后将待分析信号x“)按这簇函数分解,则根据展开的系数就可以知道信号z(f)在某一局部时间内位于某局部频段的信号成分有多少,从而实现可调窗口的时、频局部分析。实际应用中,一般采用Mallat算法。利用该方法进行小波分解时,高频信号和低频信号各占一半宽的频带,接下去再把低频部分分成两个相同宽度的频带,高频部分也作类似分解,依次类推。小波变换因为时频窗可调,对信号的高频成分采用较窄的时间尺度窗,因此具较好的高频特性;相反,对低频成分,采用较宽的时间尺度窗,自适应性很 第2章环境搬动信号分析强,所以近年来小波分析已经成为研究时频分析的重要工具。但原则上,小波燮换目STFT一样受至l不确定往骧理懿黻翻,在辩域鞫獗域不蕤嗣融获餐{壬意的分辨率;且需要构造不同的小波基以满足不同信号的分析,不同的小波基具有不同的性质,澍信号的分析能力不同,对同一信号采用不同的小波基得到的结鬃基本没鸯可毙牲;另终,,l、渡变换是≤#遭应淫戆,,l、波基一虽选定,在夔令信号分析过程中只能采用同~个小波基。2.1.2非线性时频分析方法二次嚣频袭示,l蠡栽鬃港亵功率落浚诧瑟寒,冀耱点蹩变换是二二次懿。葵{x(f)=at■(≠)+(.12X2(F)(2,13)则有e(co)=|薛112s{(搿)+}抒2|2£2(09)+2Re[al罐24Xlf&}).誓:(茹)】(2.14)式中s(∞)、8I@)与s2(to)分别为xO)、扎(O与X2∽的能量谱;X1(co)与X2沏)分别为X1(f)与x:(幻的频谱。当一(f)与如(f)的频谱是随时间变化的,其时频表示只;(≠,毋}与毛(f,掰)是二次懿,有C(t,co)=h12鼻。(f,国)+l玎212t:(f,09)+2Re[a。口:只。靠,脚)】(2.15)其中只,.O,co)为互时频表示。1>维格缡分奄信号x(f)的维格纳分布定义为畋(},岱)兰£坤+》’9~≯”如(2-16)两信号一(f)与x2(f)的互维格纳分布定义为暖,。(∽)兰n(H争:‘(卜≯”dr(2.17)菪x(t)=q墨(t)+a2x2(f),嬲眠(f,co)=川2氓,(,,∞)+蚓2%(f,oD)+2Re[a。口:暇心(f,m)】(2.18)缍藉纳分布其骞疆磐瓣辩菝特瞧,攥爨辩频窑翳澎状识然是圈霆鹣。时频繁中性高,能得出信号的瞬时频率和群速度,但交叉项干扰严重,影响其分辨率。12 第2章环境振动信号分析许多消除交叉项干扰的方法都牺牲了浚方法本身的一些良好性质。信号x(t)的Cohen分布定义为蹦,∞)=去册-2f(t-u)嘶一砌+≯(“一≯⋯d每dudr(2.19)式中妒(f,r)表示核函数,它决定了PAt,oJ)的特性。核函数分固定核函数和自适应核函数,固定核函数计算简单,但很多情况下固定核函数在模糊平面f:与信号的自身项及交叉项的分布并不匹配;自适应核函数能很好的实现抑制交叉项干扰的作用,但计算速度慢,丑.自适应核的调节能力也是相对的。2.1.3其它时频分析方法1)Cohen.Posch类正值分布时频表示只(f,co)应尽量类似能量密度的性质,满足三个特性:①PAt,珊)≥o;02-≥-J"Px(t,co)rico=№)l2;③f只(f,co)at=Ix(co)12。L.C。hen和T.Posch基于在此基础上提出了一类时频表示只p,co)=lx(f)12Ix(co)l2ta[u(t),u(国)](2.20)其中q“,u】是满足Q扣,u]≥o,ff1[u,o]du=1,lf2[u,v]do=1的任意函数。显然,核函数Qm(f),v(co)]与信号无关,故如何选择Q陋,u]是Cohen—Posch正值分布的关键问题。2)高阶时频表示L.Stankovic等人在维格纳分靠的基础上提出了L_维格纳分布,其定义为LW,(t,CO)兰D≮+矗)工+(f一2T三,-jeor如(221)其中L是一个大于零的任意整数。另一类对维格纳分布的扩展是高阶维格纳谱,其定义为13 第2章环境振动信号分析WHOS(t,%%⋯,执)三『J一Ⅳ高阶时频的表示还处于不断地研究阶段。以上这些时频分析方法在对非平稳非线性信号的处理中都有其各自的优缺点,哪种方法更适合于桥梁结构环境振动响应信号的分析尚待进一步的研究。本文采用美国华裔科学家黄谔提出的Hilbert—Huang变换方法对环境振动信号进行分析。2.2Hilbert.Huang变换“13Hilbert—Huang变换涉及到本征模式函数、经验模式分解、Hilbert谱分析,本节对这些内容进行了阐述,并讨论了Hilbert—Huang方法用于非平稳非线性信号处理的有效性和有待改进的方面。2.2.1瞬时频率(InstantaneousFrequency)‘川瞬时频率不能被接受,原因有二:①频率一般针对平稳信号而言,代表具有恒定振幅的简谐信号的周期性,可以推广到瞬时频率;而非平稳信号不存在周期性,进行频谱分析缺乏物理基础;(D定义瞬时频率的方法不唯一。所以要使用瞬时频率概念必须对信号予以限制。理论上,单分量信号可以理解为在任意时刻,仅存在一个频率值;但目前还没有对单分量信号的明确定义,无法判断一个信号是否为单分量信号。因此提出采用“窄带”对信号进行限制而使瞬时频率有意义。钟佑明等”“”1也以一振动模型为例对瞬时频率的物理意义进行了探索。频率是表征信号交变的基本变量,指单位时间内物体往复振动的次数。如图2.1(a)所示,当物体以恒定的角速度∞沿半径a。的圆周运动,物体在x轴上的投影P的运动是一简谐振动x(t)=吼cosot(2.23)14如控m娩矿。D岍●一+一t一上Ⅲ+瞰。兀H‰。∑㈦上㈧ 第2章环境振动信号分析Ⅳ(f)∞docosmJ.c,一产毒j/}了≮~.il孵y、、一/-⋯多‘经一L、』|x(t)一.f/’!/7飞。一。’(6卜w艘振动x(t)=一t)cosO(t)黼2+1瓣对频率貔浆理意义撼述频率f=∞,2厅,跟甜成弼比。当角速度拼越大,P往复振动越快,频率f越大。因此“频率”本质上可以理解为一个表示往复振动快慢的物理量,即瞬时频率。渡±影戒豹平稳售号,冀骥时壤率登鲶穗等。这是麸全嚣意义土分辑熬。;实际中半经和角速度往往都为交爨,如图2.1(b)所示,投影P的运动轨迹表示为X(f)=a(t)cosO(t)(2。24)形成的信号为非平稳信号。这时瞬时频率表明了信罨交变快慢的岗部特性,必颁转到局部意义上来。若要使国=掣能始终表示信号x◇)=矗(t)cos秽@的交交特,殛,调顿信号口(f)相对于被调制信号eosO(t)应为缓交信号,也就怒满足窄带信譬的要求,篡取值不应影响eosO(t)的单调性,即x(t)与eos0(t)应舆有同向单调性,亦即X(t)(cosO(t))‘≥0(2.25)得出本征模式信号应满足的条件:(D口(f)》0;②口’(f)≥0;(DX@(cos联国。≥0。对瞬时频率的解释最终都归结为对倍号的条件限制,即信号威是窄带信号,调幅信号玎(,)相对被调制信号eosO(t)威为缓变信号。对满足窄带祭件的信号,通过Hilbert交换褥到熬璐褥筵号是嚷~的。求出甥掇信号再对英鞠位求导,褥澍一个其有频率量缨酶参爨,在满是单值性静条锌下,这个参量W以定义药凝时频率。{S 第2章环境搬动信号分析2。2。2本征模式函数(IntrinsicModeFunction,简称IMF)出瞬时频率静物理意义知,不可以对任意豹信号讨论瓣时壤率。耩≯。提出穗信号分解成能够满足瞬时频率限制条件疆求的分量。定义一个有意义的瞬时频率的必要条件是函数必须关于局部的零均值对称,并且拥有相等数目的零交叉点帮投擅点。若信号满足以_F两个祭件:该信号的极值点数与零交叉点数相等或至多相差一个;在任意时刻,该信号的极大值点弓极小值点确定的包络线关于时间轴对称。鄂出极大德点定义的上包络线瑶。◇)和出极小傻点定义的下龟络线拜。。@,在任意露劐有掰。。。9)十日。mO)£0(2.26)则称该信号为本征模式信号,相应的函数称为本征模式函数。第一个条件与传统“窄豢”条舞耱似;焉第二个条件瑟l欤债绫瓣全局意义l:转到了局部意义上,理想化的第二个条件要求信号局部的均值为零。对于非平稳信号,“局部均值”出“局部肘间尺度”计算得出,丽对于频率而畜的这个“局帮时闫足度”壤念还难数定义,掰潋爨懑采露援大镶极蠢、嫠定义黩络线捷弱豁对称。“本征模式函数”的命名源于其代淡了信号本身的振动模式。幽零交叉点定义的每一周期的IMF仪包贪一个振动攒式。圈2。2燕⋯个IMF分缀:W抄、Ⅳ‰嗽M八~7j⋯。幽2.2一个IMF分量2。2.3经验模式分舞(EmpiricalModeDecomposition,麓称EMD)经验模式分解的的实现过程如下:的∞柏∞o∞∞∞∞{(k.耋c息萎童8。v 第2章环境振动信号分析给定信号x(f),找出x(r)所有的极大(小)值点,并将其用三次样条函数拟合成原始信号的上(下)包络线;上下包络线的均值线为原始信号的均值包络线m.(,);将原始信号xO)减去m,(t)后即得到一个新的信号hi(t):x(t)一ml(t)=乜(f)(2.27)这个过程称为筛分,原始信号经过一次筛分变成h。(f);再将h.(t)视作原始数据重复处理k次得到‰(,),当满足筛分停止准则限制标准偏差肋:争随骘堂堑(2.28)篙雏.(r)其中SD一般取0.2—0f3之问,筛分停止得到第1个IMF分量c.(r)。将原始信号x“脯去C.(f)得到残余信号_(f),再对l(f)进行上述筛分处理后得到第2个IMF分量,2(f)⋯⋯如此重复直到最后一个残余信号在设定的阈值下不可再分为止。原始信号进行经验模式分解得x(f)=∑c。(f)+rAt)。(2.29)『_l2.2。4Hilbert谱分析简单的说Hilbert变换强调了信号的本身特性,特点是突出局部属性。“”'。‘Hilbea变换器的单位抽样响应矗(")为:坳,=警=偿篇眩。。,设序列x(H)的Hilbert变换是圭(n),则挑砌m∽=昙塞掣由x(行)和Z(行)可构成x(")的解析信号:z(n)=x(竹)+∥(H)(2.31)(2.32)Hilbert变换有三个基本的性质:性质1:序列x(n)通过Hilbert变换器后,信号频谱的幅度不发生变化。因为Hilbert变换器是全通滤波器,引起频谱变化的只是其相位。性质2:序列x(n)与圭(胛)相互正交。性质3:若z(”)、一(n)与茁:(H)的Hilbert变换分别是王(行)、王.(聆)与兰:(n),且17 第2章环境搬动信号分析x(雄)=x,(露)8蔗2(聍),掇9有毒(摊)=羔I(≈)4x2(")=xl(栉)4生2(n)Hilbert谬分析中忽臻了残余函数r,阔为残余函数楚单调函数或为常数值,能量有可能很高,但是频率非常低;而研究中仪对存在于其它高频率分量中的信息感兴趣。黯予分壁穆学c(}≥,焚Hitbert交换是义是:川:三Pr弛f(2.33)其中P为柯瑟主值。c(f)与y(t)合成解搬信号:◇)=c(f)+秒◇)(2,34)用极坐标表水为z(t、=a(t)e’8(‘’浚解析信号的时交辐僮霹裙位为d(f)=√c2(,)+y2(r)皴、:afc鼬艘!~e和)对相位函数求导,得瞬时角频率脚fn:dOJ_2~摩相应的瞬时频率为邝)=瓦1百dO(t)Hilbert谱可黻表示威//(co,,):Re芝q(f弘’fq∽4=Ⅳ(∞o),r)*掰o)jml(2.35)(2.36)(2137)(2f38)(2.39)(2.40)汇总驻有分爨。,◇=1,2,.群)鹣Hilbert落藏鼙疆褥蘩爨始售号懿Hilbert邂:H(oJ,f)=∑HJ(∞,,),21(2.41)避一步还霹以定义逮嚣蘧A(国)=£H(蛾t)dt(2.42) 第2章环境振动信号分析边际谱提供了对于每个频率的总振幅的量测,表达了在整个时间长度内累积的振幅。另外,作为希尔伯特边际谱的附加结果,可以定义希尔伯特瞬时能量:1E(t)=fH2(co,,w掰(2.43)瞬时能量提供了信号能量随时问的变换情况。积分,可以得到希尔伯特能量谱:删(国)=r日2(co,Oat事实上,如果振幅的平方对时间(2.44)希尔伯特能量谱提供了对于每个频率的能量的量测,表达了每个频率在整个时间长度内所累积的能量”1。黄锷指出,无论是希尔伯特边际谱还是希尔伯特能量谱,所得到的频率与傅立叶分析中所得到的频率物理意义是完全不同的。2.2.5Hilbert.Huang变换的先进性Hilbert-Huang变换是最新发展起来的处理非线性非平稳信号的时频分析方法。。3_嘲’”1该方法的表现出显著的优越性:①通过Hilbert变换得到的瞬时频率具有清晰的物理意义:对非平稳信号进行频谱分析须定义瞬时频率,瞬时频率有多种描述方式,其中以Hilbert变换为基础定义的瞬时频率,与Fourier变换定义的频率是相容的。Hilbert.Huang变换可以看作是Fourier变换的一般化“3。②没有固定的先验基底,是自适应的。指出其幅值允许改变,突破了传统上将幅值不变的简谐信号定义为基底的局限,使信号分析更加灵活多变;③EMD分解首先分解出信号高频,然后依次分解出低频,较低频,分解是有规律的:(少各IMF分量包含了信号的本征特性。由于从信号本身的尺度特征出发对信号进行分解,本征模式函数是基于序列数据的时间特征而得出的,不同的时间序列得出不同组的本征模式函数,每一个本征模式函数可以看作是信号中一个固有的振动模态,能够良好地表达信号的局部特征:一旦获得一个IMF分量,下一个IMF分量将不会在同一时间有相同频率;(DHilbert.Huang变换是非平稳非线性的信号处理方法,利用该方法还可以定义信号的非平稳程度,并用来对非平稳信号进行滤波和去噪。19 第2章环境振动信号分析2.2.6Hilbert,Huang变换中有待解决的问题Norden嚣.Huang据密静Hilbert-Huang交换,斌予了瓣时频率合壤静定义秘物理意义,并给出了其求法,以EMD分解后的各个本征模式函数为研究对象,提出用瞬时频率、瞬时幅值等瞬时参量来表征信号特性,进行时域分析。但由,该方法提裹不久,还存在着嚣要瓣决秘完善戆一些润遂“⋯。q)EMD算法的改进由于EMD算法在计算上下包络均值时,用到了三次样条,三次样条插值沓来静运题楚曼£冲衣欠冲。为了克骚遮一缺隆,有孥髫提出了基予僖号露域弱部特征酶自逶应时变滤波分解算法帮中德定理算法,这两种算法都在提高分解速度和减少能量泄漏上取得⋯定效果,进一步研究仍在进行。(萤EMD分解的过程中边界条{牛阀越不少国内乡}学者8“51酃跨EMD分解邋程中静边界祭件闻题送行了研究。翔粜信号较短鼠给定区间内的极大值和极小值数量较少,通过三次样条插值拟合而成的上下觎络线与真实的样条包络线相比会出现严重失真。随麓筛分的不断述幸亍,逮赛楚豹误差遣会鞭餐不颧趣内转撵,逮嚣污染弱蠹部懿数撵,餐最磊的结果失去意义。但对予信号较长且给定区间内的极大值和极小德数量较多情况而言,这种影响相对就念小很多,可以通过不断撇弃两端的数据来数据来控铡误差。针对EMD分解阅题中敬边界润题,邓搠军等提出裂弱撩经网络技术对原始信号逶幸亍延拓(邵捅军等8,2001):纪跃渡等研究了截断信号逑赛静”奇布斯”现象,并提出一种命名为端值延拓法的新的边界处理方法(纪跃波等o。,2002);张郁山等提出了应用回归模型“边筛分,边延拓”处理方法(张郁L【j等”“,2003)。逶鬻楚瑾边赛效应润嚣戆爨辇方法还夜餐零延螽法、建嚣延菇法和光滑延拓法等。(D本征模式函数分离标准的选择Hilbert-Huang交换中豹经验模式分瓣豹奉囊是遴嚣簿选,姨数爨中分裹翻满足条件的本征模式函数。第一个条件容易满足,丽第二个条件则很难满足,需要确定一个标准使得这一分离过程能够停下来。黄锷等提出通j遘限制标准羲的大小来确定,即将SD馕定在O.2与O.3之间⋯3。本征模式聪数分离终止标准取得不蠢,分离出静本糕模式函数的个数帮振疆键务异。 第2章环境搬动信号分析但Hilbert-Huang方法鞠创新以来已经显示了其强大的生命力及处理非线性、嚣乎稳数据懿有效力,:蓑二遮澎应矮弱滚髂力学、建震落号努褥、绫梅毽瘴诊甑、敞障诊断等领域“1’“”1⋯。“““”1“”1。㈣。在许多领域其分析效果完全可以和小波变换方法媲美。具有很大的研究价值和广阔的应用前景。2.3本章小结1)针对环境振动响应信号遮~非平稳非线性的随机信号,阐述了时频分折方法:线性时频分毒秀方法(重羹短涎傅立时交换、,l、渡交换)、≤#线瞧对簇分携方法<热维格纳分布、Cohen分布)殿其它时频分析方法(如lCohen—Posch炎正值分布、高阶时频表示),并分析了对各种方法用于非平稳非线性信号处理中存在的优缺点。2)当信号是窄带信号,且满足调幅信号aq)相对被调制信号cosO(t)为缓变信号时两个限制祭件时,就可以对该信号定义瞬时频率。从瞬时频率引出了与Hilbert,}tumag变换旃关的本征模式函数(IMF)概念,指出IMF应满足嬲个象件:该信号静稷篷点数与零交叉点数裙等或至多相差一个;在任意露刻,该信号静极大值点与极小值点确定的镪络线关于时间轴对称。3)详细阐述了Hilbert—Huang变换的组成部分经验模式分解(EMD)和Hilbert港分橡。EMD分舞遂程:我售号投大(夺)筐点;潮三次祥条函数壤合上下霞络线;用原始信号减去均僚包络线实现筛分过程,熏复处理赢到满足IMF的两个条件,得剿第1个IMF;对原始信号减去第1个IMF分量后得划的残余信号遴{亍雪;};分处理褥篆2个tMF;⋯⋯重复处理直到在设定的阙僮下不霹再分为止。Hilbert谱分桥觚Hilbert交换的概念和髅葳入手,介绍了Hilbert—Huang变换中谗分析的实现,说明Hilbert能量谱可提供对于每个频率的能量的量测,并指出不管是Hilbert边际谱还是Hilbert能量谱,所得到的频率与傅立时分析中所得到的簇搴穆瑾意义建完全不圈懿。4)指出了Hilbert-Huang变换的先进慷:Hilbert变换得到的瞬时频率具有物理意义;EMD分解没有固怒的先验基底,是自适应的;EMD分解蟹现规律性;务IMF分量镪翕了售号熬本廷跨犍;Hilbert-Huang交换是j≥乎稳嚣线整豹襞号处理方法。讨论了有待解决的问题如EMD算法的改进、EMD分解的边界条件2l 第2章环境搬动信号分析问题、本征模式函数分离标准。最后肯定其应用前景。参考文献[1]谭善文,秦树人,汤宝平.Hilbert—Huang变换的滤波特性及其应用[J].重庆大学学报.2004,27(2):9~12【2]公茂盛,弦|毒L.HHT方法庭姥震f疆孛瓣嶷崩之秘步探讨【羽。{}l=秀瑰震工程。2003,19(3):39~43[3]钟佑明,豢树人,汤宝平.一种振动信号新变换法的研巍[J].振动L糕学报.2002,15(2):233~238[4]罗毒峰,磊誊蓉。Hilbert—Huang交换理论及其计算中的阉趱[封,强济大学学报.2003,3lf辞:637~640[5]黄大吉,赵进平,苏纪兰.希尔伯特一黄变换的端点延捅[J].海洋学报.2003.25(1):1~11【8]胡聿贤,}&郝由,粱建文.应疆HHT方法飙竖囱造震记录识别场地波他翅物理过程[弱,遮震~£程与i糕擐动.2004,24《3):i~11[7]钟佑明,察树人,汤宝平.Hilbert—Huang变换中的理论研究[J].振动垮冲击.2002,2I“):13~17C83邓拥军,王伟,钱成春等.EMD方法及Hilbert变换中边界问题的处理[J].科学通摄。2001;46(3):257~263[9l纪跃波,秦树入,柏林等.有限区闯信号边界效应问题的研究[J].振动与冲击.2002,21(4):108~111[10]张郁山,梁建文,胡聿贤.应用自回归模型处理EMD方法中的边界问题[j].自然辩学邂震.2003,13(10):t054~t059[1l】HuangN.E.,ZhengS。,Long,etal。TheempiricalmodedecompositionandtheHilbertspectrumfornonlinearandnonstationarytimeseriesanalysis[c】.Proc.R.Soc.London1998,A(454):903~995[技]J.N.Yang,Y+Lei,S.Lin,andN.Huang,Hibert-Huangbasedapproachforstructuraldamagedetection[J]。JournalofEngineeringMechanics,ASCE,2004(1):85~95L13]LOHCH,ApplicationoftheEmpiricalModeDecomposition-Hilbertspectrummethodtoidentifynear~faultground—motioncharacteristicsandstructuralresponses[J].BulletinoftheSeismologicalSocietyofAmerica,2001,9l(5):1339~1357。【14]YuDejie,ChangJunsheng,YangYu。ApplicationofEMDmethodandHilbertspectrumtothefaultdiagnosisofrollerbearings[J].MechanicalSystemsandSignalProcessing,2005,19(2):259~270[15]PaulA.Hwang,NordenE.Huang,DavidW.Wang.Anoteonanalyzingnonlinearandnonstationaryuneanwavedata[J]。AppliedOceanReseamh,2003,25:187~193lt6]JaidevaC.Goswami,AlbertE。HoefeI.Algorithmsforestimatinginstantaneousfrequency[J].SignalProcessing,2004,84:1423~1427 第2章环境振动信号分析[17]S.J.Loutridis.Damagedetectioningearsystemsusingempiricalmodedecomposition[J].En磬neeringStructure,2004,26:l8331841[18]A.D.Veltcheva,C,GuedesSoares。IdentificationofthecomponentofwavespectrabytheHilbertHuangtransformmethod[J】.AppliedOceanResearch,2004,26:1~12【19JJ,C.Nunes,Y.Bouaoune,E.Delechelle,eta1.Imageanalysisbybi·dimensionalempiricalmodedecomposition[J].ImageandvisionComputing,2003,21:1019~1026[203MarcusDatig,TorstenSchlurmann.PerformanceandlimitationsoftheHilbert-Huangtransformation(HHT)withanapplicationtoirregularwaterwaveIJ].OceanEngineering,2004,3l:1783~1834[21]E.Douka,L.j.Hadjileontiadis.Time-frequencyanalysisofthefreevibrationresponseofabeamwithabreathingcrack[J]。NDT&EInternational,2005,38:3~10[22]朱昀.水声信号非线,陡分析方法研究[D].西安:两北工业火学,2002【23]张贤达.现代信号处理[M].北京:清w华大学出版社,1995[24]于宏禹.非乎稳随机信号分析与处理[M].长沙:国防工业出版社,1999[2s3彝毅楚。谨方法戆数壤分辑[M].麓裘:秘学氆摄亭±,2000C2s】沈腻麟,叶中付,钱玉薰.信号统计分析与处理[M】.台肥:中国科学技术出版社,2001 第3章模式谚!别第3章模式谖剩3.1模式识别概述模式谈翻楚穰餍诗算撬对菜些秘褒对象遗行分类,在锫误率最小豹条{譬下,使识别结果岛客观事物相符。裟式H≯黼蜘H翥雾瓣僦鼢已翔模式项处理已翔摸式拦奉戆特征形成鼗翔模式魏特槎襻率避择和提取鹫3.1模式谈剩过程、!竺栗输85模式分类I模式识剐的基本概念;类别、模式、特征以及决策区域与决策函数。类别是与概念或原型相关的“自然”状态或对象种类。模式是⋯般选用对象麴“凌理”攘透,逶豢是铸号、霆缳或筵擎豹数篷表,又穗为聪蒙、案鲷或群本。一般用n袭示一个模式识别系统设计中可能获得的模式总数,也称为训练集或设计集。特征分三类:①物理的;(参结构的;(蛩数学的。物理和结构的特{正通常容易被感官器官所发现。数学特征较复袋,但随着计簿橇静发震数学特征应蠲越来越多,例如统计平均饿、相关函数、协方差函数游。研究的擞点就是根据学习样本选择并提取数学特征。例如一个模式可以用~个特征向凝表示:∈爿cR“(3.1)式中,x是d维特征空闻。模式识别还涉及到一个概念:决策区域与决策函数。模式识别系统的主要作24 第3章模式识别用是按各个类别将特征空间划分为对应的区域,这些区域称之为决策区域。也就是判断各个模式的所属类别。例如一个二类问题,将模式x划分为∞和国,两类。x=(x.,x:)7,用X,、x:作为雀标变量,则各个不同模式所绘成的平面分布如图3.2所示。如果这些分属于q和国:两类的模式可用方程a(x)=0来划分,d(x)=0即为决策函数。那么将某一未知类别的模式代Xd(x),如为正值,则属于以类,如为负值,则属于峨。,2d(x)=0I:}I?‘气.|,一1xc模J£,●●_tJ’·。!‘··..。·1j.。:;鼻·‘。I,:‘‘!i.一^图3.2两类模式的判别函数可以用心脏专家利用心电图仪信号的解释来诊断心脏状况的例子对模式识别加以说明。诊断类别分为四类:正常(N)、左心室肥大(LvH)、右心室肥大(RVH)、心肌梗塞(MI);与每个类别(状态)对应的模式(信号、图像等),这尾一个心电图信号就是一个模式。从每个模式中提取表示它的信息,即特征,该例中特征是关于心电图中信号的振幅、周期等的计算,如心电波包上的Q波和R波的振幅比(Q/R)。如图3.3表示但是从图中可以看出,出现了两种情况的重叠:N类类别●巴、脏状况)模式(心电图信号)Q/RlQ,R2Q/R3特征(振幅、岗期表示)图3.3用心电图特征对心脏状态分类的模式识别空间矛NLVH类的心电图会得到相同的Q/R值;MI类$13RVH类可能得到相同的心电图。第 第3章模式识别一萃孛蓬叠可以逶避增糯獒谴特征或我蓟更努静蒋鬣蕊敷勰决;第二释麓叠是这种方法固有的缺点。实际上,心电仪最好的专家系统也有一个错误率l二限(使用12种心电黼标准的系统分类错误率约为23%)。所|三l认为一个模式识剐系统蘸达到或超过人类专家灼能力就是令人潢意豹。实际应用中模式分焱系统存在的缺陷是:使用的特征不适当或不充分;用来设诗分类器的诺本不是足够全蘧秘具骞代表性;分类器瓣效率不是擐薅;类别之瞄存在交集以至于分类器无法区分。3.2模式识别方法模式识别系统的目标是要在表示空间或解释空问之间寻找到一种映射关系。棂撵不同憝分类方式,露分袁;1)统计模式识别和句法模式识别统计模式识别即利用客体的特征,研究各种划分方法,来判髓客体的归属。适合于特征数值化的场台,魏地震波解释,结构状态识别与分手厅等。句法方法即对结构进行结构描述分析,将复杂模式邂级分解为若干简单的子集合,并模仿溪言学中旬法的层次结构,运滔形式语富与垂动援技术识剩。2)监督式分类与非龄督式分类箍督式分类要依靠融知繇淄类剐翡诩练样本繁,按它稻特衽淘量的分布采镌定判别函数,从而对未知的模式进行分类议别。这种方法要求肖足够的先验知识。菲监督式分类通常采用聚类分析方法对无先骏知识情况进行分析,用数学方法分析嚣特征向量之间的鼹离及分激情提从两分类。通常应用较广泛的模式识别方法有:西配柱度值句法规则(蔫露结撼)(蒸元端擒)癸关商妇描述/势娄描缝/分娄SⅡlSS.US慝SC’i藩』{S。C{F]晤]幽3。4模式识别方法:S~监督学习;U~非监督学习:SC一统计方法;NN一神经嘲络DC一数据聚类;S卜结构匹配度组一相一鹾山一 第3章模式识别3,3特征选择与提取3,3.{特键选铎与疆繇概念特征选择指从一组特征中挑选出一燃最有效的特征以达到降低空间维数的嗣的,即从L个发量值集合{X,,x:,...,x,}中,按某一凇则选出供分类的子集作为舞维(1}l维,m
此文档下载收益归作者所有