小波的几个术语及常见的小波基介绍

1、小波的几个术语及常见的小波基介绍本篇是这段时间学习小波变换的一个收尾，了解一下常见的小波函数，混个脸熟，知道一下常见的几个术语，有个印象即可，这里就当是先作一个备忘录，以后若有需要再深入研究。一、小波基选择标准小波变换不同于傅里叶变换，根据小波母函数的不同，小波变换的结果也不尽相同。现实中到底选择使用哪一种小波的标准一般有以下几点：1、支撑长度小波函数W、3(0>)尺度函数&t)和4(的支撑区间，是当时间或频率趋向于无穷大时，w(t)、w(5)Mt)和。(Mt一个有限值收敛到0的长度。支撑长度越长，一般需要耗费更多的计算时间，且产生更多高幅值的小波系数。大部分应用选择支撑长度为5

2、9之间的小波，因为支撑长度太长会产生边界问题，支撑长度太短消失矩太低，不利于信号能量的集中。这里常常见到紧支撑”的概念，通俗来讲，对于函数f(x)，如果自变量x在0附近的取值范围内，f(x)能取到值；而在此之外，f(x)取值为0,那么这个函数f(x)就是紧支撑函数，而这个0附近的取值范围就叫做紧支撑集。总结为一句话就是除在一个很小的区域外，函数为零，即函数有速降性2、对称性具有对称性的小波，在图像处理中可以很有效地避免相位畸变，因为该小波对应的滤波器具有线性相位的特点。3、消失矩在实际中，对基本小波往往不仅要求满足容许条件，对还要施加所谓的消失矩(VanishingMoments)条件，使尽量

3、多的小波系数为零或者产生尽量少的非零小波系数，这样有利于数据压缩和消除噪声。消失矩越大，就使更多的小波系数为零。但在一般情况下，消失矩越高，支撑长度也越长。所以在支撑长度和消失矩上，我们必须要折衷处理。小波的消失矩的定义为，若jfP中df=0其中，W(t)为基本小波，0<=p<N。则称小波函数具有N阶消失矩。从上式还可以得出，同任意n-1阶多项式正交。在频域内表示就是3(在3=0处有高阶零点(一阶零点就是容许条件)。4、正则性在量化或者舍入小波系数时，为了减小重构误差对人眼的影响，我们必须尽量增大小波的光滑性或者连续可微性。因为人眼对不规则”(irregular误差比平滑”误差更加

4、敏感。换句话说，我们需要强加芷则性"(regularity除件。也就是说正则性好的小波，能在信号或图像的重构中获得较好的平滑效果，减小量化或舍入误差的视觉影响。但在一般情况下，正则性好，支撑长度就长，计算时间也就越大。因此正则性和支撑长度上，我们也要有所权衡。消失矩和正则性之间有很大关系，对很多重要的小波(比如，样条小波，Daubechies小波等)来说，随着消失矩的增加，小波的正则性变大，但是，并不能说随着小波消失矩的增加，小波的正则性一定增加，有的反而变小。5、相似性选择和信号波形相似的小波，这对于压缩和消噪是有参考价值的。二、常见的小波基以下列出的15种小波基是Matlab中支

5、持的15种。小波函数HaarDaubechiesBiorthogonalCoifletsSymletsMorletMexicanHatMeyer小波缩写名haardbbiorcoifsymmorlmexhmeyr表小形式haardbNbiorNr.NdcoifNsymNmorlmexhmeyr举例haardb3bior2.4coif3sym2morlmexhmeyr止交性有有无有有无无有双止交性有有有有有无无有紧支撑性有有一有有有无无无连续小波变换可以可以可以可以可以可以可以可以离散小波变换可以可以可以可以可以/、可以/、可以可以但无FWT支撑长度12N-1重构:2Nr+1分解:2Nd+16N

6、-12N-1有限长度有限长度有限长度滤波器长度22NMax(2Nr,2Nd)+26N2N-4,4-5,5-8,8对称性对称近似对称不对称近似对称近似对称对称对称对称小波函数消失矩阶数1NNr-12NN-尺度函数消失矩阶数-2N-1-小波函数GausDmeyerReverseBiorCgauCmorFbspShan小波缩写名gausdmeyrbioNr.Ndcgaucmorfbspshan表小形式gausNdmeyrbioNr.NdcgauNcmorfbspshan举例gaus3dmeyrbio2.4cgau3cmorfbspshan紧支撑止交性无无无无无无无紧支撑双正交性无无有无无无无连续小波

7、变换可以/、可以可以/、可以/、可以/、可以/、可以离散小波变换/、可以可以可以/、可以/、可以/、可以/、可以对称性对称对称对称对称对称对称对称小波函数消失矩阶数-尺度函数消失矩阶数-Nr-1-1、Haar小波Haar,一般音译为哈尔”。也是最简单的一Haar函数是小波分析中最早用到的一个具有紧支撑的正交小波函数,个小波函数，它是支撑域在tC0,1范围内的单个矩形波。Haar小波在时域上是不连续的，所以作为基本小波性能不是特别好。在Matlab中输入命令waveinfo('haar')可得到如下信息:Generalcharacteristics:Compactlysuppor

8、tedwavelet,theoldestandthesimplestwavelet.scalingfunctionphi=1on01and0otherwise.waveletfunctionpsi=1on00.5,=-1on0.51and0otherwise.CompactsupportyesFamilyHaarShortnamehaarExampleshaaristhesameasdb1OrthogonalyesBiorthogonalyesDWTpossibleCWTpossibleSupportwidth1Filterslength2Regularityhaarisnotcontinuo

9、usSymmetryyesNumberofvanishingmomentsforpsi2、Daubechies（dbN）、波（紧支集正交小波）Daubechies,一般音译为多贝西”。Daubechies小波是由世界著明的小波分析学者IngridDaubechies（一般音译为英格丽多贝西）构造的小波函数，我们一般简写成dbN,N是小波的阶数。小波函数邙t）和尺度函数旭）中的支撑区为2N-1,W（t）的消失矩为NodbN小波具有较好的正则性，即该小波作为稀疏基所引入的光滑误差不容易被察觉，使得信号重构过程比较光滑。dbN小波的特点是随着阶次（序列N）的增大消失矩阶数越大，其中消失矩越高光滑性就

10、越好，频域的局部化能力就越强，频带的划分效果越好，但是会使时域紧支撑性减弱，同时计算量大大增加，实时性变差。另外，除N=1外，dbN小波不具有对称性（即非线性相位），即在对信号进行分析和重构时会产生一定的相位失真。dbN没有明确的表达式（除了N=1外，N=1时即为Haar小波）。在Matlab中输入命令waveinfo（'db'）可得到如下信息：Generalcharacteristics:Compactlysupportedwaveletswithextremalphaseandhighestnumberofvanishingmomentsforagivensupportwi

11、dth.Associatedscalingfiltersareminimum-phasefilters.FamilyDaubechiesShortnamedbOrderNNstrictlypositiveintegerExamplesdb1orhaar,db4,db15OrthogonalyesBiorthogonalyesCompactsupportyesDWTpossibleCWTpossibleSupportwidth2N-1Filterslength2NRegularityabout0.2NforlargeNSymmetryfarfromNumberofvanishingmoments

12、forpsiN3、Symlet(symN)/、波(近似对称的紧支集正交小波)Symlet小波函数是IngridDaubechies提出的近似对称的小波函数，它是对db函数的一种改进。Symlet小波系通常表示为symN(N=2,3,8)。symN小波的支撑范围为2N-1,消失矩为N,同时也具备较好的正则性。该小波与dbN小波相比，在连续性、支集长度、滤波器长度等方面与dbN小波一致，但symN小波具有更好的对称性，即一定程度上能够减少对信号进行分析和重构时的相位失真。在Matlab中输入命令waveinfo('sym')可得到如下信息：Generalcharacteristic

13、s:Compactlysupportedwaveletswithleastasymmetryandhighestnumberofvanishingmomentsforagivensupportwidth.CWTAssociatedscalingfiltersarenearlinear-phasefilters.FamilySymletsShortnamesymOrderNN=2,3,.Examplessym2,sym8OrthogonalyesBiorthogonalyesCompactsupportyesDWTpossiblepossibleSupportwidth2N-1Filtersle

14、ngth2NRegularitySymmetrynearfromNumberofvanishingmomentsforpsiN4、Coiflet(coifN)/、波根据R.Coifman的要求，Daubechies构造了Coiflet小波，它具有coifN(N=1,2,3,4,5)这一系列。Coiflet的小波函数W的2N阶矩为零，尺度函数帕)的2N-1阶矩为零。W(t)和蛆)的支撑长度为6N-1。Coiflet的照t)和小具有比dbN更好的对称性。在Matlab中输入命令waveinfo('coif)可得到如下信息：Generalcharacteristics:Compactlysu

15、pportedwaveletswithhighestnumberofvanishingmomentsforbothphiandpsiforagivensupportwidth.FamilyCoifletsShortnamecoifOrderNN=1,2,，5Examplescoif2,coif4OrthogonalyesBiorthogonalyesCompactsupportyesDWTpossibleCWTpossibleSupportwidth6N-1Filterslength6NRegularitySymmetrynearfromNumberofvanishingmomentsforp

16、si2NNumberofvanishingmomentsforphi2N-15、Biorthogonal(biorNr.Nd)小波为了解决对称性和精确信号重构的不相容性，引入了双正交小波，称为对偶的两个小波分别用于信号的分解和重构。双正交小波解决了线性相位和正交性要求的矛盾。由于它有线性相位特性，所以主要应用在信号与图像的重构中。通常的用法是采用一个函数进行分解，用另外一个小波函娄进行重构。双正交小波与正交小波的区别在于正交小波满足<Wj,k,Wl,m>=亚6,m,也就是对小波函数的伸缩和平移构成的基函数完全正交，而双正交小波满足的正交性为<当上，W1m>=也，也就是对

17、不同尺度伸缩下的小波函数之间有正交性，而同尺度之间通过平移得到的小波函数系相应的滤波器也不能由同一之间没有正交性，所以用于分解与重构的小波不是同一个函数,个小波生成。该小波虽然不是正交小波，但却是双正交小波，具备正则性，同时也是紧支撑的，其重构支撑范围为2Nr+1,分解支撑范围为2Nd+1。biorNr.Nd小波的主要特征表现在具有线性相位特性。一般来说为了获得线性相位，需要降低对于正交性的局限，为此该双正交小波降低了对于正交性的要求，保留了正交小波的一部分正交性，使小波攻得了线性相位和较短支集的特性。在Matlab中输入命令waveinfo('bior')可得到如下信息：Ge

18、neralcharacteristics:CompactlysupportedbiorthogonalsplinewaveletsforwhichsymmetryandexactreconstructionarepossiblewithFIRfilters(inorthogonalcaseitisimpossibleexceptforHaar).FamilyBiorthogonalShortnamebiorOrderNr,NdNr=1,Nd=1,3,5rforreconstructionNr=2,Nd=2,4,6,8dfordecompositionNr=3,Nd=1,3,5,7,9Nr=4,

19、Nd=4Nr=5,Nd=5Nr=6,Nd=8Examplesbior3.1,bior5.5SupportwidthFilterslengthOrthogonalnoBiorthogonalyesCompactsupportyesDWTpossibleCWTpossible2Nr+1forrec.,2Nd+1fordec.max(2Nr,2Nd)+2butessentiallybiorNr.NdldlreffectivelengtheffectiveofLo_DofHi_Dbior1.122bior1.362bior1.5102bior2.253bior2.493bior2.6133bior2.

20、8173lengthbior3.144bior3.384bior3.5124bior3.7164bior3.9204bior4.497bior5.5911bior6.81711psirec.RegularityforyesNr-1andNr-2attheknotsSymmetryNumberofvanishingmomentsforpsidec.NrRemark:bior4.4,5.5and6.8aresuchthatreconstructionanddecompositionfunctionsandfiltersarecloseinvalue.6、ReverseBior小波由Biorthog

21、onal而来，因此两者形式很类似。MatlabGeneralcharacteristics:CompactlysupportedbiorthogonalsplinewaveletsforwhichsymmetryandexactreconstructionarepossiblewithFIRfilters(inorthogonalcaseitisimpossibleexceptforHaar).FamilyBiorthogonalShortnamerbioOrderNd,NrNd=1,Nr=1,3,5rforreconstructionNd=2,Nr=2,4,6,8dfordecomposit

22、ionNd=3,Nr=1,3,5,7,9Nd=4,Nr=4Nd=5,Nr=5Nd=6,Nr=8Examplesrbio3.1,rbio5.5OrthogonalnoBiorthogonalyesDWTCompactsupportyespossibleCWTpossibleSupportwidth2Nd+1forrec.,2Nr+1fordec.Filterslengthmax(2Nd,2Nr)+2butessentiallyrbioNd.NrlrldeffectivelengtheffectivelengthofHi_DofLo_Drbio1.122rbio1.362rbio1.5102rbi

23、o2.253rbio2.493rbio2.6133rbio2.8173rbio3.144rbio3.384rbio3.5124rbio3.7164rbio3.9204rbio4.497rbio5.511rbio6.81711Regularityforpsirec.yesMeyerShortnamemeyrOrthogonalyesNd-1andNd-2attheknotsSymmetryNumberofvanishingmomentsforpsidec.NdRemark:rbio4.4,5.5and6.8aresuchthatreconstructionanddecompositionfunc

24、tionsandfiltersarecloseinvalue.7、Meyer小波Meyer小波的小波函数和尺度函数都是在频率域中进行定义的，它不是紧支撑的，但它的收敛速度很快。在Matlab中输入命令waveinfo('meyr')可得到如下信息：Generalcharacteristics:Infinitelyregularorthogonalwavelet.FamilyBiorthogonalyesCompactsupportnoDWTpossiblebutwithoutFWTFIRbasedapproximationprovidesFWTCWTpossibleSuppor

25、twidthinfiniteEffectivesupport-88RegularityindefinitelyderivableSymmetryyes8、Dmeyer小波Dmeyer即离散的Meyer小波，它是Meyer小波基于FIR的近似，用于快速离散小波变换的计算。在Matlab中输入命令waveinfo('dmey')可得到如下信息：Definition:FIRbasedapproximationoftheMeyerWavelet.FamilyDMeyerShortnamedmeyOrthogonalyesBiorthogonalyesCompactsupportyesD

26、WTpossibleCWTpossible9、Gaussian小波Gaussian小波是高斯密度函数的微分形式，它是一种非正交与非双正交的小波，没有尺度函数。在Matlab中输入命令waveinfo('gaus')可得到如下信息：Definition:derivativesoftheGaussianprobabilitydensityfunction.gaus(x,n)=Cn*diff(exp(-xA2),n)wherediffdenotesthesymbolicderivativeandwhereCnissuchthatthe2-normofgaus(x,n)=1.Famil

27、yGaussianShortnamegausWaveletnamegaus"n"OrthogonalnoBiorthogonalnoCompactsupportnoDWTnoCWTpossibleSupportwidthinfiniteEffectivesupport-55Symmetryyesneven=>Symmetrynodd=>Anti-Symmetry10、MexicanHat(mexh)小波MexicanHat函数为Gauss函数的二阶导数。因数它的形状像墨西哥帽的截面，所以我们称这个函数为墨西哥草帽函数。它在时域和频率都有很好的局部化，但不存在尺度

28、函数,所以此小波函数不具有正交性。在Matlab中输入命令waveinfo('mexh')可得到如下信息：Definition:secondderivativeoftheGaussianprobabilitydensityfunctionmexh(x)=c*exp(-xA2/2)*(1-xA2)Familywherec=2/(sqrt(3)*piA1/4)MexicanhatShortnamemexhOrthogonalnoBiorthogonalnoCompactsupportnoDWTnoCWTpossibleSupportwidthinfiniteEffectivesup

29、port-55Symmetryyes11、Morlet小波Morlet小波是高斯包络下的单频率正弦函数，没有尺度函数，是非正交分解。在Matlab中输入命令waveinfo('morl')可得到如下信息：Definition:morl(x)=exp(-xA2/2)*cos(5x)FamilyMorletShortnamemorlOrthogonalnoBiorthogonalnoCompactsupportnoDWTnoCWTpossibleSupportwidthinfiniteEffectivesupport卜44Symmetryyes12、ComplexGaussian/

30、卜波属于一类复小波，没有尺度函数。在Matlab中输入命令waveinfo('cgau')可得到如下信息:Definition:derivativesofthecomplexGaussianfunctioncgau(x)=Cn*diff(exp(-i*x)*exp(-xA2),n)wherediffdenotesthesymbolicderivativeandwhereCnisaconstantFamilyComplexGaussianShortnamecgauWaveletnamecgau"n"OrthogonalnoBiorthogonalnoCompa

31、ctsupportnoDWTnoComplexCWTpossibleSupportwidthinfiniteSymmetryyesneven=>Symmetrynodd=>Anti-Symmetry13、ComplexShannonWaveletsshan在Matlab中输入命令waveinfo('shan')可得到如下信息：Definition:acomplexShannonwaveletisshan(x)=FbA0.5*sinc(Fb*x)*exp(2*i*pi*Fc*x)dependingontwoparameters:FbisabandwidthparameterFcisawaveletcenterfrequencyTheconditionFc>Fb/2issufficienttoensurethatzeroisnotinthefrequencysupportinterval.FamilyComplexShannonShortnameshanWaveletnameshan"Fb"-"Fc"OrthogonalnoBiorthogonalnoCompactsupportnoDWTnocomplexCWTpossibleSupportwidthinfinite14、ComplexFre