CL多小波预处理方法在故障数据压缩中的应用

CL多小波预处理方法在故障数据压缩中的应用

摘要：在介绍CL(Chui-Lian)多小波基本理论的基础上，探讨了CL多小波的预处理方法及其对CL多小波原有滤波器响应的影响，分析后认为Haar和平衡预处理方法是CL多小波最有效的预处理方法。文章还做了将具有这两种预处理方法的CL4多小波应用于电力系统正弦信号和故障暂态信号数据压缩的仿真实验。仿真结果表明，基于Haar和平衡预处理方法的CL4多小波具有较GHM多小波和传统db4小波更好的数据压缩效果。

关键字：CL多小波预处理方法数据压缩电力系统故障1引言

多小波分析是一种基于小波理论的近几年发展起来的新理论，多小波可同时具有对称性、正交性、短支撑性、高阶消失矩等属性，而这些属性是传统实系数小波不能同时具有的[1]。多小波有许多构造方法，如Geronimo等人[2]应用分形插值方法构造了具有短支撑、正交性、对称性和二阶消失矩属性的GHM多小波，Chui等人[3]利用多小波的正交性、紧支撑性、对称性和插值性构造了CL(Chui-Lian)多小波，Jiang[4]利用时频分析中的窗函数性质构造了具有最优时频分辨率的Jiang系列多小波，MariantoniaCotronei等人[5]利用Hurwitz块矩阵和Gram矩阵构造了半正交多小波。本文在介绍CL多小波理论的基础上，深入探讨了CL多小波的预处理方法，并将其应用于电力系统正弦信号数据和故障暂态数据的压缩，还比较了GHM多小波与CL多小波的数据压缩效果。

2CL多小波的基本理论

小波分析中的多分辨率即是将平方可积信号f∈L2(R)的逐级逼近视为采用低通平滑函数φ(t)对f(t)作平滑滤波的结果，且逐级逼近时平滑函数φ(t)也作逐级伸缩。一个多分辨率分析由一个尺度函数生成，且包含一个经平移与伸缩构成L2(R)空间基的小波函数。类似地，多小波分析中也存在多分辨率分析，一个多分辨率分析由多个尺度函数生成，且包含多个经平移与伸缩构成L2(R)空间基的小波函数，这些小波函数即称为多小波[6]。

多小波的多尺度函数φ(t)和多小波函数Ψ(t)满足以下二尺度矩阵方程[3]式中

0≤k≤L，Hk和Gk为r×r维系数矩阵；L为多小波滤波器长度；r为多小波维数。根据多小波的多分辨率分析，有如下快速多小波分解与重构公式[7]

式中

为多小波分解和重构的低频系数；为多小波分解和重构的高频系数；分别为Hk和Gk的复共轭矩阵。

位于区间[0,2]上的CL多小波的两个尺度函数φ1(t)、φ2(t)和两个小波函数Ψ1(t)、Ψ2(t)的支撑区间均为[0,2]；位于区间[0,3]上的CL多小波的两个尺度函数φ1(t)、φ2(t)和两个小波函数Ψ1(t)、Ψ2(t)的支撑区间均为[0,3]。为方便起见，作者根据Daubechies系列小波的定义方式，将上述两种CL多小波暂称为CL3多小波和CL4多小波(即它们的滤波器长度分别为3和4)，CL4多小波的滤波器系数矩阵[3]为

式中

H0、H1、H2、H3为低通滤波器系数矩阵；G0、G1、G2、G3为高通滤波器系数矩阵；下标0、1、2、3表示系数矩阵的次序，与传统小波的滤波器系数序列相同；S=diag[1,-1]。根据多小波的尺度矩阵方程绘出CL4多小波的尺度函数φ1(t)、φ2(t)和小波函数Ψ1(t)、Ψ2(t)的时域波形和频域响应波形分别如图1和图2所示。与传统小波相比，CL多小波具有更为优良的属性：CL多小波的尺度函数和小波函数均具有紧支撑属性，使其具有良好的局域性；两个尺度函数分别与两个小波函数对称和反对称，保证其具有线性相位；CL多小波是正交的，使其变换后保持能量衡定；CL3多小波具有二阶逼近阶，CL4多小波具有三阶逼近阶，使其具有良好的逼近性能，GHM多小波与CL4多小波相似，滤波器长度均为4，但其逼近3CL多小波的预处理方法

3.1采用预处理方法的必要性

与传统小波相比，多小波在实际应用中必须解决的关键问题是对原始信号的预处理。预处理的关键问题是：①由于多小波的尺度函数和小波函数是多维的，而需处理的信号一般是一维的，必须对所有多小波的原始信号进行预处理；②不同的预处理方法对多小波应用性能的影响非常大，怎样根据应用需要选择相应的最优预处理方法是多小波应用的关键问题。

如不采用预处理方法，简单的方法是将一维信号分解为其多相形式

式中z=ejω。CL4多小波的低通滤波器响应H1(ω)和H2(ω)分别为同理可求出其相应的高通滤波器响应G1(ω)和G2(ω)，如图3所示。由图3可见，CL4多小波的低通滤波器和高通滤波器响应均表现出带通和带阻两种不相同的属性，其他多小波也有类似属性，这种现象易导致多小波的低通和高通频带相互混叠，不利于分解和重构信号，即如不采用有效的预处理方法，多小波不能像传统小波那样对信号很好地进行去噪和压缩等处理。

目前多小波的预处理方法主要分为两类：预滤波(prefilter)法[8]和采用平衡多小波(balancedmultiwavelet)法[9]。

3.2预滤波法

对于预滤波方法的研究主要集中在GHM多小波的预滤波，本文将文[5]提出的odd/even、deriv.、Haar、mod.Haar预滤波法和文[10]提出的预滤波法应用于CL4多小波。经大量仿真分析后，作者认为采用Haar法可取得较好的滤波效果，图4为Haar预滤波方法对CL4多小波原有滤波器响应的影响。由图4可见，Haar预滤波法在一定程度上改善了CL4多小波两个低通滤波器和两个高通滤波器的响应。

3.3平衡多小波法

多小波低通部分如有不同的频谱属性，会使通道不平衡和复杂化矢量化过程，矢量化过程的多相方法产生了混合粗细分辨系数，仅用粗分辨系数重构信号时将产生强烈的振荡，称为不平衡现象[9]。消除多小波不平衡现象的方法是构造平衡多小波，主要有两种构造方法：采用复Daubechies小波系列滤波器构造平衡多小波和平衡已有的不平衡多小波(简称平衡法)。

以下采用平衡方法来平衡CL4多小波，平衡算子为，平衡后CL4多小波的低通和高通滤波器分别为

CL4多小波的低通合成算子LT和高通合成算子TT满足等式

器葛计式填中款

庙博u竞1叙=秋[变.好.军.脏,犯1经,洽1夏,秆1怕,趴1女,羊.锣.许.移]源T费画，众即辽矢扛量钞邮在察低浓通芦支绸路右茧保校持提不辰变鹿，迁在援高影通及支慨路而被撕取伴消券，编这傅时榴C青L弟4笋多惰小隙波味即波是淋平劳衡既的融多吊小凉波熔[捕9虽]摆。搜图身5穷为辣平促衡舟法锋对院C淡L辟4惯多摘小枣波谈原树有京滤桂波冒器跌响县应僚的壳影缓响劝。杂翼穷浅由罪图猪5身可湖见橡，乓对甘C普L斧4掀多谨小蜘波披采舍用叛平已衡况法酸不慢仅欧可渡使扇两遭个写低果通藏滤搬波映器感和挖两刚个访高系通妻滤危波油器芳的蜜响励应约分还别菜重妈合缩，枪还图可伙较拳好甘地猴改向善逼系米统捕本惯身文的列低适通揉、字高帖通鸡滤音波桐器脆的时响怀应怎性厅能轿。

窝4括

啦地C慕L索多做小订波晒的过应色用

扮植胀本巨文朝采春用珠不慌同奸预土处古理肤方缩法取的先C库L酸4三多熄小竿波枕、楼G魂H扰M异多顽小俯波而和翅d吼b划4络小延波而来菊实威现蜓对追数脆据齐的棚压普缩腰，挤火G效H条M渔多勿小庄波懂和迅d任b怀4纲小撇波县与忙C洗L注4透多叔小课波代非某常巩类艺似扮，轧且欢具枯有躲正迅交册性棚和烧紧揉支甩撑闸性猎，镜滤束波峡器每长刺度速均掏为米4虎。趋关朗于私传忽统脾小言波赶在蚁电扇力室系故统秆故吩障糕数派据经压纷缩椅中粥的冒应窄用验可叠参挑见舞文杏[腹1弊1同]祸。

糠壶京采皆用户电商力腥系如统膨的熊正但弦纷信荡号才和屡5流0违0芒k法V炎高句压挨输住电恭线街路腔单券相因接础地喉短投路羡故雄障丝相冠的即电扎压柿信必号寨(愉故油障呢点讯距火输咱电追线纳首套端哀的贝长民度惧与熄输愚电活线路总话长僚度自的绢百猎分钳比或分勾别纽为敞2影.驴5型%朵、旷1否2伟.事5慰%奶、迈2魔5婆%辫、旦3遣7阴.师5犹%尖、档5龙0塘%商、孙7眯5垄%句、衣8质7哑.串5勾%芦)牲作服为竖原移始庆信受号酱，疤采播样狐点称为笔1狼0独2握4璃个擦，策分局解贩层天数挖为辆6缝层舅。居对健C欠L蒸4歌多盾小谋波称采爬用俭H搏a厦a羊r影法捐和胜平圆衡亲法垒两拳种由预亚处遵理团方化法补；邀对才G袖H铃M臂多时小估波煌采存用踩G郑H耳M战.鞋i作n戚i婚t蚁预需处樱理冷方辈法陈(屯G脑H挂M骡.舟i棚n贺i耽t带是蹈G肿H荡M之多驴小轨波雕的握所率有厚已婆知跟预恼处医理镇方侵法筐中碌综掀合误效桌果尤最掩好练的也预磨处婆理额方卫法嚼之岸一锡)渴；赛对罗d姐b茶4厉传顽统挺小支波度采辩用箱直抖接者处久理征采油样黎点伪的吊方毒法性。忽本叔文筒采福用四保置留真多丈小究波偶分孤解额后防的屋若笔干知最厘大赞系弓数扣的稻数辫据爬压便缩柜方头法鹊。其该嘴方纱法符的厦原缘理止是怎：膨对量原暮始勿信真号艳进霞行食小忠波柏尺柔度陕的零扩跃展命，肝保错留竭绝没对挖值键最择大特的惩系灰数歪，同这翼种跑情且况辞下首，舱可话仅抵用惜全麻局摊阈蕉值泻来泪压鉴缩跃信讨号伟，坚以拐实畏现膝信必号量的展压欣缩悄或璃相改对冠均烘方辟差浇赋蹄范变信征号扒的艳恢梨复袍，守该并方轰法欣可难先贼确袍定它数辉据脏的丢压战缩币比警。扣定贫义庄信团号某重盯构宏后铃的盆赋乐范出均测方怒误卡差急为

朗式杰中伯

舰闲f环(虑n叔)汪、谱扛分倡别蔬为迹原智始唐信婶号项与失恢匀复嫌信币号法。猴述表鞋1俊为葵采副用朗该倾方彩法汽时计，纷采项用最不苍同咱多电小杂波究对投相孩应卵信荷号旱数辽据沈的掌压期缩刚结弃果罗。捕沙弦上计保算添出厌其嫁他圣故围障奶点斤信梯号卧数悲据肾压闸缩铺的芦赋劝范均误毯差秘、项平零均好误横差伟和半均中方规差割，陕得惕到惰压股缩钩比当为他2兴0赚:消1合时堂随赛不肢同禾故隙障乌地众点皱变触化壮的餐赋强范涌误陷差境曲岁线唇如裹图户6绕所闷示增。六5艰0棚%兴处裤故施障少的副暂法态贝信击号轻随沈不虎同习压辈缩锤比扰变案化担的近赋触范摇误洋差厅曲杆线凡如庆图雄7顾所镇示袋。高叼励般由似表菜1毯、堪图行6拒和妈图警7躺中慢压层缩哄效盼果斯参昆数绕可间知辱：

清占追(幸1名)慰G现H上M侵多街小端波磁及熄C职L冻多购小焦波膊对卸平缝稳希信贼号厌(殖正牛弦汽信育号烤)沃或疗暂杰态悬信湿号柔(搅故饼障臂信避号倦)君的陪数驻据储压茄缩森效干果花均袭比扣d婶b为4营传彻统矛小吸波哀的破好秋；

驳晶妥(燥2项)拍采朱用做H弊a伪a服r废和地平睡衡息法谎预木处宗理老方盒法闲的顷C田L调4疾多领小鸭波牲的刷压谋缩得效华果摘均对明搬显旷比方采聋用缓G蜜H引M巧.轮i纲n物i婶t贴预浴处州理著方鉴法督的球G须H氧M后多仆小顶波黄的沈压住缩坛效棵果岔好显；

屠休乘(类3短)拴C煎L沙4喘多芦小解波饶采扩用燥平择衡弓预想处雀理形方云法导时酸，术压脖缩斧效火果贷总毙体夸上淋比沃H嘱a楚a驼r典预失处皇理亮方傻法塑的意好正一党些换。

板5说甩结奥论惕锯悟本岗文咳介何绍临和帅分穿析株了给C仍L视多煌小母波瓦及论其豪不乘同螺的涌预躺处狡理购方优法愈，椒采逢用皮保稿留澡分恼解差后捞若喇干窃最熄大于系薯数酷的忌数醋据活压景缩交方裕法罗，旷对稀采皂用独不彻同赢预使处白理误方限法轿的财C药L庄多煮小巨波业在斧故粗障签数棋据光压膊缩词中常的收效鸭果咐与惭G独H漂M稻多治小奇波黎及炉传吐统气d凯b冷4注小非波慢进眨行杨了匙比企较疮。容综虚合介各盗种划压乱缩距性驰能攻指注标既后胳认邻为棒，赴基钟于源H善a踩a徐r梳法绍和减平老衡丝法罩的至C缸L到4虫多晚小岗波就比雹基此于军G籍H慌M殖.和i潮n烫i愉t跟法伤的响G任H昼M割多怀小倒波察压踢缩躁效胁果赌优解越仅，捐更膜适舱用预于槐对订电汽力崖系奶统艰信苗号辛的散数迹据插压摔缩倍。参考文献

[1]Daubechies．Tenlecturesonwavelets[M]．Philadelphia：SIAM，1992．

[2]GeronimoJS，HardinDP，MassopustPR．Constructionoforthogonalwaveletsusingfractalinterpolationfunctions[J]．SIAMJournalMathematicAnalysis，1996，27：1158-1192．

[3]ChuiCK，LianJA．Astudyoforthonormalmulti-wavelets[J]．AppliedNumericalMathematics，1996，20(3)：273-298．

[4]JiangQT．Orthogonalmultiwaveletwithoptimumtime-frequencyresolution[J]．IEEETransonSignalProcessing，1998，46(4)：830-844．

[5]CotroneiM，MontefuscoLB，PuccioL．Multiwaveletanalysisandsignalprocessing[J]．IEEETransonCircuitsandSystems-II：AnalogandDigitalSignalProcessing，1998，45(8)：970-987．

[6]GoodmanTNT，LeeSL．Waveletsofmultiplicity[J]．TransAmericanMathematicalSociety，1994，342：307-324．

[7]XiaXG，GeronimoJS，HardinDPetal．Designofprefiltersfordiscretemultiwavelettransforms[J]．IEEETransonSignalProcessing，1996，44(1)：25-34．

[8]StrelaV，HellerPN，TopiwalaPetal．Theapplicationofmultiwaveletfilterbankstoimageprocessing[J]．IEEETransonImageProcessing，1999，8(4)：548-563．

[9]LebrunJ，VetterliM．Balancedmultiwavelettheoryanddesign[J果]奸．却I幼E悠E获E岗悟T杆r才a乖n凶s叙鲜o聋n驼约S至i渔g鞠n岗a拔l拨仿P浴r炸o秧c芽e两s北s肿i忆n恶g扛，息1尿9帖9评8揭，诸4鸡6授(夹4暴)洞：脂1笔1鲜1荒9眼-低1旱1蚀2悦4址．

膨[侮1速0治]钩旷T绕h卖a秆m蝶溜J通扣Y蒙，则S轨h轻e斧n伸畜L殃识X蛮，惊L扒e女e鹊逼S耽挤L奇踏e每t旱哪a茶l皇．骑A则赚g飘e岂n贯e植