五种常用小波基含MATLAB实现

1、1给出五种常用小波基的时域和频域波形图与标准的傅里叶变换相比，小波分析中使用到的小波函数具有不唯一性，即小波函数中(t)具有多样性。小波分析在工程应用中，一个十分重要的问题就是最优小波基的选择问题，因为用不同的小波基分析同一个问题会产生不同的结果。目前我们主要是通过用小波分析方法处理信号的结果与理论结果的误差来判定小波基的好坏，由此决定小波基。常用小波基有 Haar 小波、Daubechies(dbN)小波、MexicanHat(mexh)小波、Morlet 小波、Meyer 小波等 5 种。(1)Haar 小波Haar 函数是小波分析中最早用到的一个具有紧支撑的正交小波函数，也是最简答的一个

2、小波函数，它是支撑域在t三0,1范围内的单个矩形波。Haar10t2函数的定义如下：甲(t)=-12t_110其他Haar 小波在时域上是不连续的，所以作为基本小波性能不是特别好。但它也有自己的优点，如：计算简单；空(t)不但与空(2jt)j亡z正交，而且与自己的整数位移正交。因此，在a=2j的多分辨率系统中 Haar 小波构成一组最简单的正交归一的小波族。中的傅里叶变换是：W(Q)=j-sin2(-)e-jQ211ahaar 时域1.5.c10.5-0-0.5-1-1.5cc:00.511.5ti=20;wav=haar;phi,g1,xval=wavefun(wav,i);subplot(

3、1,2,1);plot(xval,g1,-r,LineWidth,1.5);xlabel(t)title(haar 时域)；g2=fft(g1);g3=abs(g2);subplot(1,2,2);plot(g3);xlabel(f)title(haar 频域)haar 频域x101015f5x10(2)Daubechies(dbN)、波Daubechies 小波是世界著名的小波分析学者 InridDaubechies 构造的小波函数，简写为 dbN,N 是小波的阶数。小波中(t)和尺度函数4(t)中的支撑区为 2N-1,中(t)的消失矩为 N。除N=1外，dbN 不具有对称性(即非线性相位)

4、。dbN 没有明确的表达式(除N=1外)，但转换函数 h 的平方模是明确的。Daubechies 小波系是由法国学者 Daubechies 提出的一系列二进制小波的总称，在 Matlab 中记为 dbN,N 为小波的序号，N 值取 2,3,，10。该小波没有明确的解析表达式，小波函数小与尺度函数的有效支撑长度为 2N-1.当 N 取 1 时便成为 Haar小波。N1,N-1kkN-1k.一.一一.令p(y)=Cky,其中Ck为二项式的系数，则有kz012N式中，mod)=r=hhke，2k4Daubechies 小波具有以下特点：(1)在时域是有限支撑的，即中(t)长度有限。(2)在频域中(s

5、)在=0 处有 N 阶零点。(3)手(t)和它的整数位移正交归一，即*(t)V(t-k)dt=6k。(4)小波函数空(t)可以由所谓“尺度函数”(t)求出来。尺度函数(t)为低通函数，长度有限，支撑域在 t=0(2N-1)范围内。2m。()22(cos-)p(sin2)db4频域10001900-800700-600-500-400300-200-1000111f02000400060008000fi=10;wname=db4;phi,g1,xval=wavefun(wname,i);subplot(1,2,1);plot(xval,g1,-r,LineWidth,1.5);xlabel(t)

6、title(db4 时域)；g2=fft(g1);g3=abs(g2);subplot(1,2,2);plot(g3,-r,LineWidth,1.5);xlabel(f)title(db4 频域)wname=db4;%计算该小波的4个滤波器Lo_D,Hi_D,Lo_R,Hi_R=wfilters(wname);subplot(2,2,1);stem(Lo_D);title(分解低通滤波器);subplot(2,2,2);stem(Hi_D);title(分解高通滤波器,);subplot(2,2,3);stem(Lo_R);title(重构低通滤波器);subplot(2,2,4);stem

7、(Hi_R);title(重构高通滤波器,);注思Daubechies 小波常用来分解和重构信号，作为滤波器使用。波形如下:(3)MexicanHat(mexh)小波22-i-MexicanHat 函数为 Gauss 函数的二阶导数：,(t)=(1-t)e2因为它的形状像墨西哥帽的截面，所以也称为墨西哥帽函数MexicanHat(mexh)小波的时域和频域波形图mexihat频域15r110-.II;I!IoLc050fd=-6;h=6;n=100;g1,x=mexihat(d,h,n);subplot(2,2,1);plot(x,g1,-r,LineWidth,1.5);xlabel(t)t

8、itle(Mexihat 时域)；g2=fft(g1);g3=(abs(g2);subplot(2,2,2);plot(g3);xlabel(f)title(mexihat 频域)；特点：墨西哥帽函数在时间域与频率域都有很好的局部化，并且满足(t)dt=0。由于它不存在尺度函数，所以小波函数不具有正交性。100t它是高斯包络下的单频率副正弦函数：2-t,(t)=Ce，co$5x)C 是重构时的归一化常数。Morlet 小波没有尺度函数4(t),而且是非正交分解Morlet 小波的时域波形图和频域波形图Morlet时域td=-6;h=6;n=100;g1,x=morlet(d,h,n);subp

9、lot(2,2,1);plot(x,g1,-r,LineWidth,1.5);xlabel(t)title(Morlet 时域)；g2=fft(g1);g3=abs(g2);subplot(2,2,2);plot(g3);xlabel(f)title(Morlet 频域)151050050100fMorlet频域iHI1(5)Meyer 小波Meyer 小波的小波函数和尺度函数都是在频率域中进行定义的,其定义为:其中，v(a)为构造 Meyer 小波的辅助函数，具有v(a)=a4(35-84a70a2-20a3)aQ11-(2二厂-2二()=(2二)cos(v(0Meyer 小波不是紧支撑的，

10、但它收敛的速度很快:2-nn(1+|t)中(t)无限可微。Iz1-3I(2n)2e-sin(-v(-1)22njwJI-1jw二3()=(2)2e-cos(-v(22二-1)2二034二034-38二3Meyer 小波的时域和频域波形图15meyer频域020406080100120140f1050d=-6;h=6;n=128;psi,x=meyer(d,h,n,psi);subplot(2,1,1),plot(x,psi,-r,LineWidth,1.5)xlabel(t)title(Meyer 时域)；PSI=fft(psi);PSII=abs(PSI);subplot(2,1,2),pl

11、ot(PSII);xlabel(f)title(meyer 频域)2、在信号x(t)=sin(2兀*30t)+cos(2兀*50t)加上噪首后分别进行FFT和CW变换。解：引入随机噪声 randn(1,N)N=100;fs=1000;n=0:N-1;t=n/fs;x=sin(60*pi*t)+cos(100*pi*t);subplot(3,2,1);plot(x,-r,LineWidth,1.5);xlabel(t)ylabel(x(t)title(原信号 x(t)波形图)F1=fft(x);m1=abs(F1);subplot(3,2,2);plot(m1);xlabel(f)0-2050原信号x(t)波形图50tx(t)加噪声后fft变换图100 x(t)的fft变换图III|.J,EJ：50100fmorlet50f尺度为110020406080100time(orspace)b尺度为264216421aQXFacsaQXFacs-100500-2cE100050100%原信号t 田 e(x的 fft 变换图)x1=randn(1,N);%加入噪声x2=x+x1;F2=fft(x2);m2=abs(F2);subplot(3,2,3);plot(m2);xlabel(f)title(