- 1、本文档共13页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
五种常用小波含MATLAB实现
给出五种常用小波基的时域和频域波形图。
与标准的傅里叶变换相比,小波分析中使用到的小波函数具有不唯一性,即小波函数 具有多样性。小波分析在工程应用中,一个十分重要的问题就是最优小波基的选择问题,因为用不同的小波基分析同一个问题会产生不同的结果。目前我们主要是通过用小波分析方法处理信号的结果与理论结果的误差来判定小波基的好坏,由此决定小波基。常用小波基有Haar小波、Daubechies(dbN)小波、Mexican Hat(mexh)小波、Morlet小波、Meyer小波等5种。
Haar小波
Haar函数是小波分析中最早用到的一个具有紧支撑的正交小波函数,也是最简答的一个小波函数,它是支撑域在范围内的单个矩形波。 Haar函数的定义如下:
Haar小波在时域上是不连续的,所以作为基本小波性能不是特别好。但它也有自己的优点,如:
?计算简单;
?不但与正交,而且与自己的整数位移正交。
因此,在的多分辨率系统中Haar小波构成一组最简单的正交归一的小波族。
的傅里叶变换是:
Haar小波的时域和频域波形图
i=20;
wav = haar;
[phi,g1,xval] = wavefun(wav,i);
subplot(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 频域)
Daubechies(dbN)小波
Daubechies小波是世界著名的小波分析学者Inrid Daubechies构造的小波函数,简写为dbN,N是小波的阶数。小波和尺度函数中的支撑区为,的消失矩为N。除外,dbN不具有对称性(即非线性相位)。dbN没有明确的表达式(除外),但转换函数h的平方模是明确的。
Daubechies小波系是由法国学者Daubechies提出的一系列二进制小波的总称,在Matlab中记为dbN,N为小波的序号,N值取2,3,…,10。该小波没有明确的解析表达式,小波函数φ与尺度函数Φ的有效支撑长度为2N-1.当N取1时便成为Haar小波。
令,其中为二项式的系数,则有
式中,。
Daubechies小波具有以下特点:
在时域是有限支撑的,即长度有限。
在频域在=0处有N阶零点。
和它的整数位移正交归一,即。
小波函数可以由所谓“尺度函数”求出来。尺度函数为低通函数,长度有限,支撑域在t=0~(2N-1)范围内。
Daubechies小波的时域和频域波形图
i=10;
wname = db4;
[phi,g1,xval] = wavefun(wname,i);
subplot(1,2,1);
plot(xval,g1,-r,LineWidth,1.5);
xlabel(t)
title(db4 时域);
g2=fft(g1);
g3=abs(g2);
subplot(1,2,2);plot(g3,-r,LineWidth,1.5);
xlabel(f)
title(db4 频域)
注意
Daubechies小波常用来分解和重构信号,作为滤波器使用。波形如下:
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(Hi_R);
title(重构高通滤波器);
Mexican Hat(mexh)小波
Mexican Hat函数为Gauss函数的二阶导数:
因为它的形状像墨西哥帽的截面,所以也称为墨西哥帽函数。
Mexican Hat(mexh)小波的时域和频域波形图
d=-6;
h=6;
n=100;
[g1,x]=mexihat(d,h,n);
subplot(2,2,1);plot(x,g1,-r,LineWidth,1.5);
xlabel(t)
title(Mexihat 时域);
g2=fft(g1);
g3=(abs(g2));
subplot(2,2,2);plot(g3);
xlabel(f)
title(mexihat 频域);
特点:墨西哥帽函数在时间域与频率域都有很好的局部化
文档评论(0)