1.有了Fourier,为什么还需要Wavelet?

先来揭揭短:

(1)Fourier分析不能刻画时间域上信号的局部特性。

(2)Fourier分析对突变和非平稳信号的效果不好,没有时频分析。

傅立叶变换将函数投影到正弦波上,将函数分解成了不同频率的正弦波,这不能不说是一个伟大的发现,但是在大量的应用中,傅立叶变换的局限性却日趋明显,事实上在光滑平稳信号的表示中,傅立叶基已经达到了近似最优表示,但是日常生活中的信号却并不是一直光滑的,而且奇异是平凡的,傅立叶在奇异点的表现就着实让人不爽,从方波的傅立叶逼近就可以看出来,用了大量不同频率的三角波去逼近其系数衰减程度相当缓慢,而且会产生Gibbs效应。其内在的原因是其基为全局性基,没有局部优化能力,以至局部一个小小的摆动也会影响全局的系数。实际应用中很需要时频局部化,傅立叶显然缺乏此能力了。即使如此,由于其鲜明的物理意义和快速计算,在很多场合仍然应用广泛。傅立叶变换在从连续到离散的情形是值得借鉴与学习的,大家都知道,时间周期对应频域离散,时间离散对应频域周期,时间离散周期对应频域离散周期,DFT其实是将离散信号做周期延拓然后做傅立叶变换再截取一个周期,反变换同样如此,所以DFT用的是块基的概念,这样如果信号两端的信号连接后不再光滑(即使两边都光滑),同样会在边界上产生大幅值系数(边界效应),延伸到图像中就是块效应。当对信号做对称周期延拓后再做傅立叶变换得到的正弦系数全部为0,也就是任何对称函数可以写成余弦的线性组合,同样按照离散的思路构造得到的是离散块余弦基,即DCT变换,虽然DCT可以通过对称后周期延拓再变换减少了边界效应(两边信号接上了,但不一定平滑),但也不能消除块效应,尤其是图像变换中人为将图像分成8*8处理后块效应更加明显。但是DCT很好的能量聚集效应让人惊奇,加之快速计算方法使它替代DFT成为图像的压缩的标准了很长时间(JPEG)。

闲扯了这么多,不要迷糊,开始上图

第一个就是傅立叶变换是整个时域,所以没有局部特征。这是由基函数决定的,同时如果在时域发生突变,那么在频域就需要大量的正弦波去拟合,这也是傅立叶变换性质决定的。

第二个就是面对非平稳信号,傅立叶变换可以看到由哪些频率组成,但不知道各成分对应的时刻是什么。也就是没有时频分析,看不出来信号频域随着时间变换的情况,反过来说就是,一个的频图对应好几个时域图,不知道是哪个,这个在实际应用中就不方便了,如图:

如上图,最上边的是频率始终不变的平稳信号。而下边两个则是频率随着时间改变的非平稳信号,它们同样包含和最上信号相同频率的四个成分。做FFT后,我们发现这三个时域上有巨大差异的信号,频谱(幅值谱)却非常一致。尤其是下边两个非平稳信号,我们从频谱上无法区分它们,因为它们包含的四个频率的信号的成分确实是一样的,只是出现的先后顺序不同(时间分辨性太差)
可见,傅里叶变换处理非平稳信号有天生缺陷。它只能获取一段信号总体上包含哪些频率的成分,但是对各成分出现的时刻并无所知。因此时域相差很大的两个信号,可能频谱图一样。
然而平稳信号大多是人为制造出来的,自然界的大量信号几乎都是非平稳的,所以在比如生物医学信号分析等领域的论文中,基本看不到单纯傅里叶变换这样naive的方法。

上图所示的是一个正常人的心电相关电位。对于这样的非平稳信号,只知道包含哪些频率成分是不够的,我们还想知道各个成分出现的时间。知道信号频率随时间变化的情况,各个时刻的瞬时频率及其幅值——时频联合分析

2.预备知识:何为基?何为內积?

2.1 基

傅立叶变换和小波变换,都会听到分解和重构,其中这个就是根本,因为他们的变化都是将信号看成由若干个东西组成的,而且这些东西能够处理还原成比原来更好的信号。那怎么分解呢?那就需要一个分解的量,也就是常说的基,基的了解可以类比向量,向量空间的一个向量可以分解在x,y方向,同时在各个方向定义单位向量e1、e2,这样任意一个向量都可以表示为a=xe1+ye2,这个是二维空间的基:

而FT的基是不同频率的正弦曲线(整个时间),所以FT是把信号波分解成不同频率的正弦波的叠加和:而对于小波变换就是把一个信号分解成一系列的小波(短时间),也许就会问,小波变换的小波是什么啊,定义中就是告诉我们小波,因为这个小波实在是太多,一个是种类多,还有就是同一种小波还可以尺度变换,但是小波在整个时间范围的幅度平均值是0,具有有限的持续时间和突变的频率和振幅,可以是不规则,也可以是不对称,很明显正弦波就不是小波,什么的是呢,看下面几个图就是:

有了基,有什么用呢?下面看一个傅立叶变换的实例:
对于一个信号的表达式为x=sin(2*pi*t)+0.5*sin(2*pi*5*t); 下面看图形表示,感受一下频域变换给人的一目了然:

基具有非冗余性,即使基不是正交的,有相关性,但若去掉其中任何一个,则不成为基,这一点也叫完备性;基的表示有唯一性,即给定一族基对一个函数的表达是唯一的;一般情况下基非正交,也称为为exact frame(Resize basis),这个时候要表示信号可以将基正交化成唯一的正交基(对偶为其自身);也可以求其对偶框架(dual frame),其对应了小波变换中的双正交情形!信号可以依框架分解,然后用对偶框架重构。若在基集里添加一些新的向量,并随意调整空间位置,则有可能成为框架。把函数与基或框架作内积,也可以说成是一种函数空间到系数空间的变换。若某种变换后的能量(内积的平方和度量)仍然有一个大于0的上下界,才可以成为框架,由于框架的冗余性,所以系数的表达也不具有唯一性。若上下界相等,则为紧框架,且界表示冗余度。若上下界相等为且为1,称为pasval identity frame,此时不一定为正交基,此时若加上基的长度均为一的条件,则框架退化为正交基。可能你会问我们用基来表示信号就行了啊,为什么还要框架呢?其实很多信号表示方法不能构成基,却能构成框架,如短时傅立叶变换中如要求窗函数满足基条件,则可推出该函数有很差的时频局部化性质(事实上退化为了傅立叶变换)

2.2 內积

如果两个向量的内积为0 ,就说他们是正交的。
如果一个向量序列相互对偶正交,并且长度都为1,那么就说他们是正交归一化的。

3.小波诞生的前一个晚上:短时傅里叶变换(STFT)

有了缺点当然就想着改进了,这就出来了短时傅立叶变换,也叫加窗傅立叶变换,顾名思义,就是因为傅立叶变换的时域太长了,所以要弄短一点,这样就有了局部性。
定义:把整个时域过程分解成无数个等长的小过程,每个小过程近似平稳,再傅里叶变换,就知道在哪个时间点上出现了什么频率了。”这就是短时傅里叶变换。下面就是示意图

时域上分成一段一段做FFT,不就知道频率成分随着时间的变化情况了吗!
可能理解这一点最好的方式是举例子。首先,因为我们的变换是对时间和频率的函数(不像傅立叶变换,仅仅是对频率的函数),它是二维的(如果加上幅度则是三维)。以下图所示的非平稳信号为例:

在这个信号中,在不同时刻有四个频率分量。0-250ms内信号的频率为100Hz,其余每个250ms的间隔的信号频率分别为50Hz,25Hz和10Hz。很明显,这是一个非平稳信号,让我们看一看它的短时傅立叶变换:用这样的方法,可以得到一个信号的时频图了:

既能看到10Hz, 25 Hz, 50 Hz, 100 Hz四个频域成分,还能看到出现的时间。两排峰是对称的,只用看一排就行。这貌似解决了问题,好像有了加窗技术,整个世界都亮了。没有那么简单的,面对一个随机的非平稳信号,那么这个窗要多大了呢?(此时加窗傅里叶变换一脸的懵逼)

窗太窄,窗内的信号太短,会导致频率分析不够精准,频率分辨率差
窗太宽,时域上又不够精细,时间分辨率低

下面通过一组实验,我们深入探索这种加窗技术,问题出现在哪里?

上图对同一个信号(4个频率成分)采用不同宽度的窗做STFT,结果如右图。用窄窗,时频图在时间轴上分辨率很高,几个峰基本成矩形,而用宽窗则变成了绵延的矮山。但是频率轴上,窄窗明显不如下边两个宽窗精确。所以窄窗口时间分辨率高、频率分辨率低,宽窗口时间分辨率低、频率分辨率高。对于时变的非稳态信号,高频适合小窗口,低频适合大窗口。然而STFT的窗口是固定的,在一次STFT中宽度不会变化,所以STFT还是无法满足非稳态信号变化的频率的需求。

4.小波来了

时势造英雄,小波开始一展拳脚了!
对于加窗傅立叶变换让人头疼的就是窗口的大小问题,如果我们让窗口的大小可以改变,不就完美了吗?答案是肯定的,小波就是基于这个思路,但是不同的是。STFT是给信号加窗,分段做FFT;而小波变换并没有采用窗的思想,更没有做傅里叶变换。小波直接把傅里叶变换的基给换了——将无限长的三角函数基换成了有限长的会衰减的小波基。这样不仅能够获取频率,还可以定位到时间了~
就又回到了最开始的基了。
小波变换采用的这些基函数会伸缩、会平移(其实是两个正交基的分解)。缩得窄,对应高频;伸得宽,对应低频。然后这个基函数不断和信号做相乘。某一个尺度(宽窄)下乘出来的结果,就可以理解成信号所包含的当前尺度对应频率成分有多少。于是,基函数会在某些尺度下,与信号相乘得到一个很大的值,因为此时二者有一种重合关系。那么我们就知道信号包含该频率的成分的多少。如前文所述,小波做的改变就在于,将无限长的正弦函数基换成了有限长的会衰减的小波基。效果如下图:

从公式可以看出,不同于傅里叶变换,变量只有频率ω,小波变换有两个变量:尺度a(scale)和平移量 τ(translation)。尺度a控制小波函数的伸缩,平移量 τ控制小波函数的平移。尺度就对应于频率(反比),平移量 τ就对应于时间。如下图

当伸缩、平移到这么一种重合情况时,也会相乘得到一个大的值。这时候和傅里叶变换不同的是,这不仅可以知道信号有这样频率的成分,而且知道它在时域上存在的具体位置。而当我们在每个尺度下都平移着和信号乘过一遍后,我们就知道信号在每个位置都包含哪些频率成分。看到了吗?有了小波,我们从此再也不害怕非稳定信号啦!从此可以做时频分析啦!

3.1解决了局部性

3.2解决时频分析

时域信号                                            FFT变换                                小波分析

有了这些,小波分析也就入门了。

感谢大连理工大学机械工程硕士小木匠的出色工作。

感谢杜克大学方圆之中对STFT的现状所做的科普工作。

此外,本文参考了:http://blog.sina.com.cn/s/blog_13bb711fd0102w9x1.html 一并致谢

小波的秘密1_小波变换概况与综述相关推荐

  1. 小波的秘密3_连续、离散小波变换定义

    1.前言 小波变换主要包括连续小波变换和离散小波变换.本篇博客主要想弄清楚连续小波变换.离散小波变换.高维小波连续变换的意义. 2.连续小波变换 2.1 连续小波变换的定义 将任意L2(R)空间中的函 ...

  2. 小波的秘密10_图像处理应用:图像增强

    1.前言: 图像融合是信息融合的一个重要分支,广泛地用于目标识别.机器视觉.智能系统.医学图像处理等领域.传统图像融合方法方法主要是在时间域通过算术运算实现融合,具有算法简单直观,适合实时处理等优点, ...

  3. 小波的秘密9_图像处理应用:图像增强

    1.前言: 图像增强问题的基本目标是对图像进行一定的处理,使其结果比原图更适合用于特定的应用领域. 在图像增强领域,图像增强问题主要通过时域和频域两种方法进行处理.时域方法通过直接在图像上作用算子来解 ...

  4. 小波的秘密8_图像处理应用:图像降噪

    1.前言: 图像去噪是信号处理的一个经典问题,传统的去噪方法多采用平均或线性法法进行,最常用到的就是维纳滤波,但是他的降噪效果并不是很明显.小波分析法开辟了非线性降噪的先河,小波能够降噪得益于小波变换 ...

  5. 小波的秘密6_小波包的理解

    1.各种变换的适合处理对象 小波变换 加窗Fourier变换   Fourier变换 突变信号或奇异性函数 自适应信号处理 处理渐变信号 实时信号处理 稳定和渐变信号 实时信号处理 2.小波包分解概述 ...

  6. 基于小波的图像边缘检测,小波变换边缘检测原理

    1.什么是"小波神经网络"?能干什么用呀 小波神经网络(Wavelet Neural Network, WNN)是在小波分析研究获得突破的基础上提出的一种人工神经网络.它是基于小波 ...

  7. 谐波小波matlab,基于MATLAB小波变换在电网谐波检测的仿真与研究

    <工业控制计算机>2012 年第 25 卷第 10 期 谐波分析是当前国内外电能质量检测研究的一个热点,其研究的难点是对突变的. 暂态的. 非平稳扰动信号的检测与分析,只利用传统的傅立叶变 ...

  8. 小波的秘密10_小波包的数学支撑

    1.前言: 小波多分辨率分析可以对信号进行有效的时频分解,但由于其尺度函数是按照二进制变化的,所以在高频段其频率分辨率较差.只能对信号的频段进行指数等间隔划分.小波包分析还能够为信号提供一种更精细的分 ...

  9. 小波的秘密2_MATLAB图像多尺度分析

    1.MATLAB程序编写步骤 1.把小波w(t)和原函数f(t)的开始部分进行比较,计算系数C.系数C表示该部分函数与小波的相似程度. 2.把小波向右移k单位,得到小波w(t-k),重复1.重复该部知 ...

最新文章

  1. php采集分析,PHP采集程序原理分析篇
  2. Recom Airbnb embedding
  3. ZeroC Ice启用SSL通讯的配置
  4. qaxobject控制word到某一行_太简单了!Excel信息导出为每页打印一行,且批量添加同一表头...
  5. 柏林噪声实践 水与火,顶点纹理拾取
  6. regexp_like 方法
  7. Windows下配置QGIS和Python
  8. Leetcode200岛屿数量(深搜)
  9. 算法洗脑系列(8篇)——第四篇 枚举思想
  10. thinkpadt410接口介绍_ThinkPadT410全程拆解
  11. w ndows7如何清理垃圾,Win7系统清理:如何清理Win7系统盘垃圾
  12. 【SOLIDWORKS学习笔记】工程图基础操作
  13. 圆柱体积怎么算立方公式_圆柱的体积换算立方怎么算
  14. SVG_16_defs标签_use标签_style标签_红绿灯效果
  15. WeUI源码学习笔记一
  16. 01背包,完全背包,多重背包
  17. 抽象工厂模式(三):抽象工厂模式概述
  18. 他整整复习了三个月,金三成功上岸字节,面试题及答案分享
  19. 人生第一笔试加面试 SAP胜略
  20. IT计算机前后端相关专业——20个常用的学习教程及在线编程工具

热门文章

  1. oracle 定时任务 job 调用存储过程有回到输出参数(含out参数)
  2. iOS 各种编译错误汇总
  3. 利用Inotify和Rsync将webproject文件自己主动同步到多台应用server
  4. AS5 安装 JAVA 1.6.0 用于 TOTOplayer的启动
  5. hdu 5067(状态压缩dp)
  6. leetcode 11. Container With Most Water
  7. CodeForces 459C(构造题)
  8. NYOJ 595 乱七八糟
  9. Java反射之反射权限和静态属性的赋值
  10. 3ds Max制作客厅场景实例教程