说到这里,可能你对scaling function以及多解析度分析已经比较理解了。但是,我们还没有涉及到它们在小波变换中的具体应用,也就是还没有回答刚才那个问题:凭空插了一个scaling function到小波basis组合中干嘛。也就是说,我们希望理解scaling function是怎么和小波函数结合的呢,多解析度能给小波变换带来什么样的好处呢。这其实就是是小波变换中的核心知识。理解了这个,后面的小波变换就是纯数学计算了。

好,我们已经知道,对于子空间V0,basis是scaling function:

对应的小波函数是:

然后子空间V1的basis集合是这俩哥们:

看出什么规律了么?多看几次这三个图,你会惊讶地发现,在V0中的scaling function和wavelet function的组合,其实就是V1中的basis!继续这样推导,V1本来的的basis是:

然后V1中对应的wavelet function是

他们的组合,本质上也就是V2的basis(参考图2)。你继续推导下去,会得到同样的结论:在scale j的wavelet function,可以被用来将Vj的basis扩展到V(j+1)中去!这是一个非常非常关键的性质,因为这代表着,对任何一个子空间Vj,我们现在有两种方法去得到它的orthonormal basis:

1. 一种就是它本来的basis

,对任意k。

2. 第二种就是它上一个子空间的basis

,对任意k,以及上一级子空间的wavelet function

,对任意k。

第二种选择能给我们带来额外的好处,那就是我们可以循环不断地用上一级子空间的scaling function以及wavelet function的组合来作为当前子空间的基。换句话说,如果针对V3这个子空间,它实际上就有四种不同的,但是等价的orthonormal basis:

1. 本级(V3)的scaling function basis set

2. 上一级(V2)的scaling function + wavelet function;

3 . 上上一级(V1)的scaling function + 上上一级(V1)的wavelet function + 上一级(V2)的wavelet function;

4. 上上上一级(V0)的scaling function + 上上上一级(V0)的wavelet function + 上上一级(V1)的wavelet function + 上一级(V2)的wavelet function

好,看看最后一种选取方式,有没有感到眼熟?对了,它就是我们之前提到的“针对此信号space的哈尔小波basis组合”,参见图1。现在我们知道了,这个scaling function不是凭空插进去的,而是通过不断的嵌套迭代出来的:)

那为什么我们最后选定的是这种选取方式呢?实际上,刚才介绍的这个性质已经告诉我们,对于任何的scale j0,我们都可以给我们的signal space找到一组orthonormal basis,这个basis是通过组合scale j0上的scaling function以及所有在scale j,j>j0上的wavelets得到的。这样,基于这个orthonormal basis,所有信号空间中的信号都可以写成组成这个basis的functions的线性组合:

对应的系数的计算和平常一样:

这,就是最终的,也是最核心的,小波变换形式。不管是信号压缩,滤波,还是别的方式处理,只要是用小波变换,都逃不出这个基础流程:

1. 选取合适的wavelet function和scaling function,从已有的信号中,反算出系数c和d。

2. 对系数做对应处理

3. 从处理后的系数中重新构建信号。

这里的系数处理是区别你的应用的重点。比如图像或者视频压缩,就希望选取能将能量聚集到很小一部分系数中的小波,然后抛弃那些能量很小的小波系数,只保留少数的这些大头系数,再反变换回去。这样的话,图像信号的能量并没有怎么丢失,图像体积却大大减小了。

还有一个没有解释的问题是,为什么要强调尺度函数和小波函数组成一个orthonormal basis呢?计算方便是一方面,还有一个原因是,如果他们满足这个性质,就满足瑞利能量定理,也就是说,信号的能量,可以完全用每个频域里面的展开部分的能量,也就是他们的展开系数表示:

到这里,我们对小波变换的形式就讲完了。虽然是用的最简单的哈尔小波为例子,但举一反三即可。我们着重介绍了多解析度分析以及它给小波变换带来的杀手锏:时域频域同时定位。结束之前,再多说几句小波变换的意义。我们拿刚才例子中V3子空间的第二种可选择的orthonormal basis作为例子:

左边这四个basis组成元素,也就是scaling functions,的系数,表征的是信号的local平均(想想它们和信号的内积形式),而右边的这四个basis组成元素,也就是wavelet functions,的系数则表征了在local平均中丢失的信号细节。得益于此,多解析度分析能够对信号在越来越宽的区域上取平均,等同于做低通滤波,而且,它还能保留因为平均而损失的信号细节,等同于做高通滤波!这样,我们终于可以解释了wavelet function和scaling function背后的物理意义了:wavelet function等同于对信号做高通滤波保留变化细节,而scaling function等同于对信号做低通滤波保留平滑的shape!

对小波变换的基础知识,我们就讲到这里。需要注意的是,这只是小波变换最基本最基本的知识,但也是最核心的知识。掌握了这些,代表你对小波变换的物理意义有了一定的了解。但对于小波变换本身的讲解,一本书都不一定能将讲透,还有很多的基础知识我都没有讲,比如如何构建自己的scaling function,选取合适的系数集h[k],并由此构建自己的wavelet functions。所以,如果有深入下去研究的同学,好好买一本书来看吧。而只是希望用小波变换来服务自己的应用的同学,个人觉得这些知识已经足够让你用来起步了。

小波变换 完美通俗讲解 之 二相关推荐

  1. 二维小波变换_小波变换完美通俗讲解系列之 (一)

    声明:该篇文章转自csdn,原始博主已经找不到了,在这里给出转载博主地址,如有侵权,请私信我删除. https://blog.csdn.net/liusandian/article/details/5 ...

  2. 小波变换完美通俗讲解系列之 (一)

    申明:小波变换 完美通俗解读,是<小波变换和motion信号处理>系列中的第一篇. 原始出处为​windstorm的网站http://www.kunli.info/,并非本站原创,但这位大 ...

  3. 小波变换 完美通俗解读

    申明:小波变换 完美通俗解读,是<小波变换和motion信号处理>系列中的第一篇. 原始出处为​windstorm的网站http://www.kunli.info/,并非本站原创,但这位大 ...

  4. 人工智能算法通俗讲解系列(二):逻辑回归

    2019独角兽企业重金招聘Python工程师标准>>> 今天,我们介绍的机器学习算法叫逻辑回归.它英语名称是Logistic Regression,简称LR. 跟之前一样,介绍这个算 ...

  5. 小波变换 完美通俗解读【转载

    转自: http://www.kunli.info/2011/02/15/fourier-wavelet-motion-signal-1/ 这是<小波变换和motion信号处理>系列的第一 ...

  6. 小波变换和motion信号处理(二)【转载】

    这是<小波变换和motion信号处理>系列的第二篇,深入小波.第一篇我进行了基础知识的铺垫,第三篇主要讲解应用. 在上一篇中讲到,每个小波变换都会有一个mother wavelet,我们称 ...

  7. 网络编程懒人入门(九):通俗讲解,有了IP地址,为何还要用MAC地址?

    1.前言 标题虽然是为了解释有了 IP 地址,为什么还要用 MAC 地址,但是本文的重点在于理解为什么要有 IP 这样的东西.本文对读者的定位是知道 MAC 地址是什么,IP 地址是什么. (本文同步 ...

  8. Hadoop平台K-Means聚类算法分布式实现+MapReduce通俗讲解

        Hadoop平台K-Means聚类算法分布式实现+MapReduce通俗讲解 在Hadoop分布式环境下实现K-Means聚类算法的伪代码如下: 输入:参数0--存储样本数据的文本文件inpu ...

  9. 冲突域、广播域的通俗讲解

    冲突域.广播域的通俗讲解 1.冲突域(物理分段)         连接在同一导线上的所有工作站的集合,或者说是同一物理网段上所有节点的集合或以太网上竞争同一带宽的节点集合.这个域代表了冲突在其中发生并 ...

最新文章

  1. 02-基于注解的入门案例
  2. 小米路由器mini改打印服务器_如何把家里的闲置路由器用起来
  3. 数据结构----归并排序
  4. PowerDesigner(二)-项目和框架矩阵
  5. Linux多线程与Linux多进程混合项目的死锁问题
  6. PTA4-11 先序输出叶结点
  7. 微软发布ASP.NET 5路线图
  8. phoneGap实际开发中的某些雷区
  9. ionic之自动生成app图标和启动页面
  10. html怎么动态背景,HTML 动态背景
  11. 【方法论】从入门到精通是怎样一种体验
  12. LS文法构图算法(3) Hilbert-Peano曲线
  13. 2017上半年计算机教学计划,2017上半年数学教学计划
  14. 蓝桥杯单片机12届第二场
  15. 爱老婆,除了性,还有更价值的东东
  16. #距离#JZOJ 3256 BZOJ 3170 洛谷 3964 松鼠聚会
  17. ubuntu系统中文件解锁方法
  18. golang slice扩容机制
  19. 特斯拉称正调查上海停车场特斯拉车辆自燃事故
  20. 如何删除 MOSS 共享服务提供程序 (SSP, Shared Service Provider)

热门文章

  1. hibernate 悲观锁乐观锁
  2. selenium找不到元素
  3. 把事务封装成类似Serializable用法的特性
  4. nginx日志自动切割
  5. POSTMAN 数据关联
  6. textarea 高度调整
  7. 08-GCD常用方法
  8. crontab+php-cgi/php 定时执行PHP脚本
  9. cocos2d-x学习笔记16:记录存储1:CCUserDefault
  10. LACP/PAGP的定义与差别—Vecloud微云