五种常见小波基函数及其matlab实现.doc

上传人:hw****26 文档编号:3820123 上传时间:2019-07-31 格式:DOC 页数:8 大小:214.92KB
下载 相关 举报
五种常见小波基函数及其matlab实现.doc_第1页
第1页 / 共8页
五种常见小波基函数及其matlab实现.doc_第2页
第2页 / 共8页
五种常见小波基函数及其matlab实现.doc_第3页
第3页 / 共8页
五种常见小波基函数及其matlab实现.doc_第4页
第4页 / 共8页
五种常见小波基函数及其matlab实现.doc_第5页
第5页 / 共8页
点击查看更多>>
资源描述

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

2、函数0,1t的定义如下: 21(t)-0t其 他Haar 小波在时域上是不连续的,所以作为基本小波性能不是特别好。但它也有自己的优点:1. 计算简单。2. 不但与 正交,而且与自己的整数位移正交,因(t)j(t)z2此,在 的多分辨率系统中,Haar 小波构成一组最简单的正交ja归一的小波族。的傅里叶变换是:()t2/24=sin()jea( ) jHaar 小波的时域和频域波形phi,g1,xval = wavefun(haar,20);subplot(2,1,1);plot(xval,g1,LineWidth,2);xlabel(t)title(haar 时域);g2=fft(g1);g3

3、=abs(g2);subplot(2,1,2);plot(g3,LineWidth,2);xlabel(f)title(haar 频域) Daubechies(dbN)小波Daubechies 小波是世界著名的小波分析学者 InridDaubechies 构造的小波函数,简写为 dbN,N 是小波的阶数。小波 和尺度函数 中的支撑区为(t)(t), 的消失矩为 。除 (Harr 小波)外,dbN 不具有对称12(t)1性(即非线性相位) 。除 (Harr 小波)外,dbN 没有明确的表达式,但转换函数 h 的平方模是明确的:令 ,其中 为二项式的系数,则有kNkypC10-()kN1- )2)

4、p(sin(co)220 m其中:ehjkNk-1200)(mDaubechies 小波具有以下特点:1. 在时域是有限支撑的,即 长度有限。t2. 在频域 在 处有 N 阶零点。)(=03. 和它的整数位移正交归一,即 。tk)dt-(t4. 小波函数 可以由所谓“尺度函数” 求出来。尺度函数为低通函数,长度有限,支撑域在 的范围内。(t)t=02N-1db4 的时域和频域波形:phi,g1,xval = wavefun(db4,10);subplot(2,1,1);plot(xval,g1,LineWidth,2);xlabel(t)title(db4 时域);g2=fft(g1);g3=

5、abs(g2);subplot(2,1,2);plot(g3,LineWidth,2);xlabel(f)title(db4 频域)Daubechies 小波常用来分解和重构信号,作为滤波器使用:Lo_D,Hi_D,Lo_R,Hi_R = wfilters(db4); %计算该小波的 4个滤波器subplot(2,2,1); stem(Lo_D,LineWidth,2);title(分解低通 滤波器 );subplot(2,2,2); stem(Hi_D,LineWidth,2);title(分解高通 滤波器 );subplot(2,2,3); stem(Lo_R,LineWidth,2);t

6、itle(重构低通 滤波器 );subplot(2,2,4); stem(Hi_R,LineWidth,2);title(重构高通 滤波器 ); Mexican Hat(mexh)小波Mexican Hat 函数为 Gauss 函数的二阶导数: 2(t)1t)et2因为它的形状像墨西哥帽的截面,所以也称为墨西哥帽函数。Mexihat 小波的时域和频域波形:d=-6; h=6; n=100;g1,x=mexihat(d,h,n);subplot(2,1,1);plot(x,g1,LineWidth,2);xlabel(t);title(Mexihat 时域);g2=fft(g1);g3=(abs

7、(g2);subplot(2,1,2);plot(g3,LineWidth,2);xlabel(f);title(mexihat 频域);Mexihat 小波的特点:1. 在时间域与频率域都有很好的局部化,并且满足 。(t)d0R2. 不存在尺度函数,所以 Mexihat 小波函数不具有正交性。 Morlet 小波它是高斯包络下的单频率副正弦函数: 2(t)Cecos(5)tt其中 C 是重构时的归一化常数。Morlet 小波没有尺度函数,而且是非正交分解。Morlet 小波的时域和频域波形图:d=-6; h=6; n=100;g1,x=morlet(d,h,n);subplot(2,1,1)

8、;plot(x,g1,LineWidth,2);xlabel(t);title(morlet 时域);g2=fft(g1);g3=(abs(g2);subplot(2,1,2);plot(g3,LineWidth,2);xlabel(f);title(morlet 频域); Meyer 小波1. Meyer 小波不是紧支撑的,但它收敛的速度很快2. 无限可微(t)Meyer 小波的时域和频域波形图:d=-6; h=6; n=128;psi,x=meyer(d,h,n,psi);subplot(2,1,1);plot(x,psi,LineWidth,2);xlabel(t);title(meyer 时域);g2=fft(psi);g3=(abs(g2);subplot(2,1,2);plot(g3,LineWidth,2);xlabel(f);title(meyer 频域);

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 实用文档资料库 > 策划方案

Copyright © 2018-2021 Wenke99.com All rights reserved

工信部备案号浙ICP备20026746号-2  

公安局备案号:浙公网安备33038302330469号

本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。