对于一个连续的周期信号,可以将其分解为一组频率不同的三角函数信号的线性组合,这就是傅里叶级数的本质,将信号从时域投影到频域中的不同频段上来完成分解。

当这个周期信号的周期趋近于无穷大时,傅里叶级数就变成了傅里叶变换。此时的信号本质上是一个连续非周期信号,傅里叶变换的意义就在于对其进行分解,同样也是以一组三角函数作为正交基,并通过这组三角函数基的线性组合来表示原信号。数学表达为:

由于三角函数是一个无限长的信号,在时域上不具有局部性,因此以其作为正交基对信号进行拟合时,具有以下两个不足:第一,对于突变信号,如阶跃信号或尖峰信号,其需要大量的三角函数基进行组合才能完成较好的信号拟合;第二,由于三角函数不具备在时域上的局部性,因此在对信号进行傅里叶变换时,仅仅只能获取到信号在频域上的分布信息,并不能获取到这些不同频率的信号分量在时域上出现的位置。因此傅里叶变换对于非平稳信号的分解会遗失其在时域上的变化信息。

小波变换就是为了解决对非平稳信号的分解问题而产生的数学方法。相比于傅里叶变换使用一组无限长的三角函数基进行信号拟合,小波变换使用的是一组正交的、迅速衰减的小波函数基进行信号拟合。这种小波函数基可通过其尺度变量和平移变量,获得不同的频率和时间位置。因此在利用这种小波函数基对信号进行分解时,可以用较少的小波函数基就拟合出突变信号(稀疏编码特性),同时也能获得不同频率的信号分量在时域上的出现位置。

用于生成一组不同频率和时移的小波函数的小波函数

,称为基本小波(Basic Wavelet),由其生成的一组小波函数,是该基本小波的一个小波族(Wavelet Family),表示为:

,其中

为尺度参数,通过伸缩控制小波的尺度(频率),

为平移参数,通过移位控制小波在时域中的出现位置。这两个参数的作用顺序是先作平移,再作伸缩。对这一族小波函数进行归一化,即得到一组小波函数基。

可进行小波变换的对象是平方可积的信号,也即是位于希尔伯特空间中的一个函数。希尔伯特空间保证了在空间中的每一个函数,都可由该空间中的其他函数线性组合得到。

小波变换值,也即是用于拟合原始信号的小波函数基的系数,可通过计算原始信号与各小波函数基的内积得到。其意义为原始信号在各个小波函数基上的投影值,投影值越大,说明对应的小波信号所携带的原始信号的特征信息的比例越大。数学表达如下:

其中,

表示

的共轭函数。

在实际应用中,由于计算机的处理对象为离散数据,因此在使用小波变换时,一般指的都是离散小波变换。因此对于上式,可写成离散形式:

其中,

表示

的共轭函数。

从数学的角度理解,在小波变换中,一个位于希尔伯特空间中的函数,可以分解成一个尺度函数和一个小波函数,其中尺度函数对应原始函数中的低频部分,小波函数对应原始函数中的高频部分。通过尺度函数可以构建对原始信号的低通滤波器,通过小波函数可以构建对原始信号的高通滤波器。

从信号处理的角度理解,在小波变换中,信号可通过信号滤波器分解为高频分量(高频子带(subband))和低频分量(低频子带(subband)),高频子带又称为细节(detailed)子带,低频子带又称为近似(approximate)子带。细节子带是由输入信号通过高通滤波器后再进行下采样得到的,近似子带是由输入信号通过低通滤波器后再进行下采样得到的。

小波变换是这样一个过程:首先将原始信号作为输入信号,通过一组正交的小波基分解成高频部分和低频部分,然后将得到的低频部分作为输入信号,又进行小波分解,得到下一级的高频部分和低频部分,以此类推。随着小波分解的级数增加,其在频域上的分辨率就越高。这就是多分辨率分析(MRA,MultiResolution Analysis)。

离散小波变换在逐级分解时,由尺度函数所张成的空间为:

其中

为第

级的尺度函数所张成的空间,

为尺度变量,

为平移变量。

为产生

一族尺度函数的父函数,又称父小波。上式中的

是归一化因子。

离散小波变换在逐级分解时,由小波函数所张成的空间为:

其中

为第

级的小波函数所张成的空间,

为尺度变量,

为平移变量。

为产生

一族小波函数的母函数,又称母小波。上式中的

是归一化因子。

关于

的正交空间补集,两者存在关系:

,根据该关系,可以递归展开,得到:

由上式可知,分解级数越高,信号在时域和频域的分辨率就越高,包含的信息也就越多。

又根据上式,对于

中的任一函数

,可分解为以下形式:

其中

中的函数,

中的函数。

在小波变换中,若令尺度参数为

,即对尺度按幂级数作离散化,同时对平移保持连续变化,则此类小波变换称为二进小波变换(Dyadic Wavelet Transform)。

在二进小波变换中,各级小波分解时,相邻级数的尺度函数之间满足关系:

,

即是:

相邻级数的小波函数和尺度函数之间满足关系:

,

即是:

其中,

为第$i$级的尺度函数,

为第

级的尺度函数,

为第

级的小波函数。

为低通滤波器,定义为:

为高通滤波器,定义为:

在小波变换中,紧支撑小波基是性质较好的一类小波基,紧支撑(Compact Support)函数是指这样的一类函数:其自变量仅在0附近的取值范围内能得到非零函数值,而在其他区间取值,则得到的函数值全为零。能得到非零函数值的自变量取值区间被称为该函数的支撑区间。一个函数的支撑区间长度主要由其尺度参数决定。支撑区间越大,计算复杂度越高,边界拖尾效应越明显。不仅如此,支撑区间越大,会产生更多的高幅值小波系数,关于这个结论的解释,可参考傅里叶变换使用无限长(支撑区间大)的三角函数基进行信号拟合的情况,相比于使用信号迅速衰减(支撑区间小)的小波基,三角函数基拟合信号时需要更多的数量。因此在选择小波基时,以支撑长度较短的小波基为宜。此外,小波基的正交性也是一类重要的性质,它确保了信号的分解没有冗余(最优分解)。

在小波变换中,另一个重要的概念是消失矩。小波函数的消失矩定义如下:

若:

则称该小波函数具有N阶消失矩。其中,

为基本小波,

。消失矩越大,在进行小波分解时,会产生更多数值为0的小波系数,使得信号的分解更为稀疏,便于进行数据压缩。但与此同时,更大的消失矩也会产生更大的支撑区间,这是一个权衡关系。

然而小波变换仍存在着不足之处,由于小波变换在每级信号分解时仅仅对低频子带进行分解,因此无法对高频子带的信息进行同样高分辨率的提取。

小波包分解(Wavelet Packet Decomposition),又称为最优子带树结构(Optimal Subband Tree Structuring)正是对小波变换的进一步优化。其主要的算法思想是:在小波变换的基础上,在每一级信号分解时,除了对低频子带进行进一步分解,也对高频子带进行进一步分解。最后通过最小化一个代价函数,计算出最优的信号分解路径,并以此分解路径对原始信号进行分解。

同样地,在二进小波包变换中,各级小波包分解时,相邻级数的尺度函数和小波函数之间也具有递推关系。

记小波包变换中的父小波

,母小波

,其中的上标表示该小波包所在的分解级数,下标表示该小波包在其级里的位置。

于是上述的递推关系可以表述如下:

即是:

其中

的定义同小波变换。

即称为小波包。

在小波包变换中,常用的代价函数为信息熵函数。最小化代价函数,也即是最大化逐级信号分解的信息熵。

无论是小波变换还是小波包变化,一个重要的环节是选择合适的小波函数基进行信号分解。常用的几类小波函数基包括:

Haar小波:Haar函数是最简单的一个小波函数,其具有紧支撑性和正交性,函数图像为在支撑区间[0,1)上的单个矩形波。Haar小波在时域上不连续,作为基本小波时性能不是很好。尺度函数 小波函数

Daubechies小波:通常简称为

,其中

为小波分解的级数。该小波的尺度函数和小波函数的支撑区间长度为

,小波函数的消失矩为

小波具有较好的正则性(本身的函数较为光滑,其拟合出的信号也较为光滑),且随着小波分解级数$N$的增加,其小波函数的消失矩也增大,小波更加光滑,在时域的紧支撑性降低(支撑区间变大),在频域的局部性增加(频域可分性更强)。此外,由于支撑区间的变大,计算量也将增加,使得计算实时性下降。当

小波的

=1时,

小波即为Haar小波;当

小波的

时,其不具有对称性(即其相位为非线性),该特性会使其在进行信号重构时产生一定程度的失真。尺度函数 小波函数

Biorthogonal小波:双正交小波,使用对偶的两个小波分别用于信号的分解和重构。双正交小波与正交小波的区别在于,在正交小波中,由基本小波的伸缩和平移所产生的一族小波基函数之间完全正交。而在双正交小波中,由基本小波的伸缩所产生的一族小波基函数完全正交,由基本小波的平移所产生的一族小波基函数则不正交。双正交小波借由牺牲一部分的正交性,来减少基本小波平移产生小波基时需满足的约束(扩大函数空间),从而使小波可以通过平移来获取线性相位(在函数空间中可取到具有线性相位的函数),以得到精确重构信号所需的对称性。双正交小波的作用过程可以描述为:首先对输入信号使用一个双正交小波进行正交分解,然后在该双正交小波的对偶空间中寻找一个具有线性相位的对偶小波进行对输入信号的非正交重构。双正交小波具备正则性和紧支撑性,其重构支撑区间长度为

,分解支撑区间长度为

分别代表信号重构的级数和信号分解的级数。左上:用于信号分解的尺度函数 右上:用于信号分解的小波函数 左下:用于信号重构的尺度函数 右下:用于信号重构的小波函数

Symlets小波:通常简称为

,其中

为小波分解的级数。

小波是对

小波的一种改进,其除了具备

小波所拥有的特性外,还具备了更好的对称性,是一类近似对称的紧支撑正交小波。因此,

小波能在一定程度上减少信号重构时产生的相位失真。

小波的支撑区间长度为

,消失矩为

上:尺度函数 下:小波函数

Mexican Hat小波:Mexican函数为Gaussian函数的二阶导数,形似墨西哥帽的截面,它在时域和频域上都具有很好的局部性,但不存在尺度函数,因此该小波函数不具备正交性。小波函数

在选择合适的小波函数基时,需要考虑的因素除了小波的支撑区间大小、小波的消失矩之外,还需考虑其对称性、正则性和相似性。小波的对称性主要体现在保证信号重构时不会产生相位畸变,即是不会产生重构信号的相位失真。小波的正则性保证了信号的光滑和可微性,对于大部分小波而言(非全部),其与消失矩存在关系:小波的消失矩越大,正则性也就越大。最后,选择与输入信号的波形相似性高的小波,意义在于使数据压缩和降噪变得更容易(信号的拟合和分解都更容易)。

关于小波分解的应用:

由于小波分解的稀疏编码特性,其可用于数据压缩,主要做法为:将信号进行小波分解,并将较小的小波系数置零。相当于将不重要(特征不明显)的信息分量去除,达到数据精简的效果。

小波分解也可用于信号滤波,主要做法为:将信号进行小波分解,并将特定级数以上的小波系数置零。相当于将高分辨率的信息分量去除,达到数据平滑的效果。

小波分解还可用于信号降噪,主要做法为:将信号进行小波分解,并通过设置一个阈值,将其中低于阈值的小波系数置零。相当于将信号中占成分比例较低的噪声部分去除。

小波包分解的应用与小波分解同理。

python小波包分解_小波包变换(Wavelet Packet Transform)的学习笔记相关推荐

  1. python小波包分解_小波包获得某个节点信号的几个细节问题

    小波包分解主要用了这四个函数  wpdec,wprec,wprcoef,wpcoef.1)Wpcoef: 求解某个节点的小波包系数,数据长度是L/(2^n)(n表示分解的层数) 比如:原始信号长度为1 ...

  2. 小波分解与小波包分解代码_分解的功能参数和代码可维护性

    小波分解与小波包分解代码 Code keeps changing, there's no doubt about that. We always do our best to set some roc ...

  3. db9小波包分解matlab,小波包分解、重构 matlab代码

    clear all close all clc % 加载信号 load leleccum; % 采样点 x = leleccum(300:2348); % x = importdata('V1.txt ...

  4. python祝福祖国代码_“小程序 大梦想”之创意编程校园邀请赛------53信息技术学科周...

    原标题:"小程序 大梦想"之创意编程校园邀请赛------53信息技术学科周 如今"人工智能"."创意编程"这些新热词已经深入人们的学习和生 ...

  5. 小甲鱼python课后题答案_小甲鱼Python课后习题-仅学习 - 随笔分类 - 小红帽爱吃大灰狼 - 博客园...

    摘要:1.字典是Python中唯一的映射类型 2.字典包含两个要素:键(key)和值(value)他们是成对出现的,用大括号括起来,多对存在时用逗号隔开. 3.可以用dict()直接创建字典,如dic ...

  6. 用python画小仓鼠教程_小仓鼠简笔画教程

    近期啊,在刷视频的时候啊,看见一个特别可爱的小松鼠,在吃东西!简直被萌坏啦!以前小时候啊,住在乡下,晚上的时候呢,经常听见小松鼠的叫声,经常就会带着一些小零食去喂松鼠,有一次,哥哥们抓住一只松鼠,激动 ...

  7. android小程序案例_小程序案例赏析:高质量的小程序怎么做

    很多新手想做小程序,但却不知道好的小程序应该做成什么样子.下面就跟大家分享几个做得比较好的微信小程序案例,你可以从这些案例中学习一下,然后再做自己的小程序. 1.商城小程序案例 商城小程序如今是比较常 ...

  8. 小尼机器人_小尼被机器人嫌弃“唱歌难听,长相一般”?

    我们如今所处的时代,科技创新的速度日新月异,生活方式多彩多姿.人人都说:科技改变了生活.今晚<开门大吉>也迎来了三大改变生活的神奇黑科技! 智能且生态的"移动城堡" 在 ...

  9. audio 小程序 放大_小程序之争:支付宝、百度“放大招”,春节欲打突围战

    2017年1月9日,微信小程序正式上线,上线后引发舆论热议.一年后,小程序再次凭借着小游戏大火,成功收获大量新用户. 如今,微信小程序已成为众多微信用户日常高频应用.腾讯之外,阿里.百度.今日头条等也 ...

  10. java小组队徽_小组队徽设计.ppt

    小组队徽设计 小组队徽设计 制作团队:F7 小组商讨设计方案 小组各成员队徽设计草图 小组队徽设计问题 初期遇到的问题:对于色彩的搭配比较迷茫,关于色彩的问题一直无法确定,还有关于色彩的明暗问题如何给 ...

最新文章

  1. PyTorch关键算法疑似侵权,Facebook被起诉
  2. SpringBoot实战(四)之使用JDBC和Spring访问数据库
  3. python基于条件、规则构建已有字典的子集
  4. android 开发框架 怎么使用,Android快速开发框架dyh详解(二)---控件层的使用
  5. KeUserModeCallback用法详解
  6. markdown转html
  7. 【light 1341Aladdin and the Flying Carpet】
  8. Java偏向锁、轻量级锁、重量级锁
  9. oracle入库的速度能到多少_多线程能提高Oracle的入库速度吗
  10. 隐藏或者修改nginx信息
  11. Netty工作笔记0011---Channel应用案例2
  12. 微软职位内部推荐-Senior Network Engineer
  13. radius认证服务器无响应,squid radius认证“RADIUS服务器没有响应”
  14. python贝叶斯网络预测模型_概率图模型之:贝叶斯网络
  15. c语言程序设计在未来前景,转型背景下C程序设计课程的实施方案
  16. 计算机组成体系结构复习笔记
  17. 软件项目管理 第六章 成本计划
  18. 入手评测 锐龙r7 5825u和i5 1240p选哪个好
  19. 视频直播:职场如何升职加薪?
  20. Sa函数 与 sinc函数

热门文章

  1. 分享一个NI软件卸载工具
  2. Linux快捷键的使用
  3. 波场地址解析易语言代码
  4. 机顶盒系统升级服务器,服务器及其软件升级方法、IP机顶盒及其软件升级方法...
  5. android wifi信号检测工具,推荐4个专业又实用的WiFi检测工具,了解一下
  6. OD教程(去除NAG窗口--PE文件结构)
  7. bandgap带隙基准电路设计,cadence,
  8. IDE、SATA、SCSI、SAS、FC、SSD 硬盘类型
  9. 【Fltk】Fltk1.3.3+VS2015 编程
  10. zbrush史上最全笔刷下载43G(1200个)