整理下时频分析变换的方法,遇见好的文章就记录下来了,本篇博客参考知乎https://www.zhihu.com/topic/19621077/top-answers上的一个回答,自己手敲一遍,增强记忆

首先说明这里是连续小波变换,不会涉及离散小波变换,不涉及尺度函数。

对于一个morlet小波变换,其频率成分分别是4、6、10,具体表达式如下:

经小波变换,得到如下时频图:

图一、二维时频图

三维时频图

从图中可以清晰辨认出原始信号含有的频率成分及其个自对应的时间区间。如果你和我一样认真的话,你会产生三个疑问?

(1)为什么小波变换能确定信号频率和其对应的时间区间?
(2)为什么频率小的条纹比频率高的条纹要细?

(3)为什么三个条纹颜色的深浅不一样 ?

其实第(1)个问题可以从上一篇的博客中总结出来。

第一个问题我们把它一分为二,为什么CWT能辨认出信号的频率成分?为什么CWT能确定频率对应的时间区间?

当然先看下公式。来一个morlet小波基函数表达式

morlet小波的基函数是由复三角函数乘上一个指数衰减函数构成,这里表示中心频率。还是先看图(分别令

图三

上图左边是基函数时域图像,右边是傅里叶变换图像。可以看到基函数的频率正是其中心频率的值。这里给出morlet小波的傅里叶变换表达式

从表达式也可以看出当频率等于其中心频率时,取极大值。这里复三角函数可以辨认频率,衰减函数可以保证其时域的有限支撑。只给一个固定的中心频率可不能辨认信号的频率,同样,基函数只在【-2,2】之间也确定不了时间区间。所以这里的小波基函数需要平移和伸缩 。又是一个公式(按传统,这个变大时前还有一个根号a,但我不讨论这个)

继续上图,给出b=5(平移因子),a=2(尺度变量),即平移5个单位,缩小2倍的morlet小波基函数的图像

图四

相对图三,从时域可以明显看出平移,从频域也很明显看出尺度伸缩。当然从傅里叶变换的性质也可以推断出来。

好了,基函数算是告一段落,下面是激动人心的时刻,还是得有公式。

这算是morlet小波变换的公式了。

好,回到第一个问题,为什么小波变换能确定信号频率和其对应时间区域?

从时频图一可以看到,能看出条纹的部分明显是关于频率的积分极大值。

所以,把问题转换一下,为什么小波变换会在原始信号固有的频率上产生极大值,

而这极大值对应的区间正好是该频率的时区。连续小波变换的运算是什么?

是积分,是原始信号与小波基函数乘积后积分的过程,而积分就是一个求和的过程

举个例子:

一个频率为5的正弦函数,和乘积后积分,问当频率为何值时,他们的乘积积分后结果最大。答案你肯定知道,当两个正弦函数频率相等的时候,而就是n=5的时候,他们的乘积积分后最大。有图为证:

好,到这里你应该就明白了为什么小波变换会在原始信号固有的频率上产生极大值。因为基函数里包含复三角函数。

下面把换成小波基函数,换成morlet小波函数,去掉平移参数,并令初始中心频率,令尺度参数

再次求其与y乘积积分后的最大值(这个得取模)。很明显,当a=2的时候,也就是小波基函数频率为5的时候,会取的极大值,如图

到这里,应该就很明朗了。小波小波,顾名思义,既要小又要有波动。morlet小波的波动性可以用复三角函数表达,小则用衰减函数表达,数学上把这种小称为有限支撑。即morlet小波的有限支撑是通过一个指数衰减函数实现的。复三角函数使其能分析频率(和原始信号乘积积分求极大值),衰减函数使其可以定位时间,他们加起来,才使得morlet小波可以用来做时频分析。

CWT就是选一个中心频率然后通过尺度变换得到一大堆中心频率,又通过时移得到一系列不同区间的基函数,分别和原始信号的某一段(对应基函数的区间)乘积再积分,产生的极值对应的频率就是原始信号这一区间含有的频率。

你不想换个角度看看吗?

换个角度,我们来从频域看看,好吗?根据傅里叶变化的性质Parseval公式可知:

两信号乘积后时域的积分等于其 傅里叶变换后频率的积分。好,换成信号和小波基函数

时域里面乘积积分看的不明显,频域里看就一目了然了,还是那个频率为5的原始信号,看看他的频域图像

再看看五个中心频率不同的小波基函数

好了,两个函数再频域乘积再积分,又一次很明显,只有当小波中心频率为5的时候,小波基函数和原始信号乘积的积分才会取得极大值。其他的好像都是0啊

所以,小波可以认为是一个带通滤波器,只允许频率和小波中心频率(经过尺度伸缩后)相近的信号的通过

连续小波变换(CWT)的更多相关文章

  1. 连续小波变换CWT(2)

    如果让你说说连续小波变换最大的特点是什么?多分辨分析肯定是标准答案.所谓多分辨分析即是指小波在不同频率段会有不同的分辨率.具体表现形式,我们回到前一篇文章的第一个图, 图一 对应的信号为 低频时(频率 ...

  2. 小波变换检测信号突变点的MATLAB实现

    之前在不经意间也有接触过求突变点的问题.在我看来,与其说是求突变点,不如说是我们常常玩的"找不同".给你两幅图像,让你找出两个图像中不同的地方,我认为这其实也是找突变点在生活中的应 ...

  3. 对小波变换中DWT和CWT的理解

    转载自 http://blog.sina.com.cn/s/blog_633750d90100hbco.html 连续小波变换的概念.操作.及时间尺度图的显示 最近很多网友问到关于连续小波变换的诸多问 ...

  4. matlab 小波变换

    MATLAB小波变换指令及其功能介绍 1 一维小波变换的 Matlab 实现 (1) dwt函数 功能:一维离散小波变换 格式:[cA,cD]=dwt(X,'wname') [cA,cD]=dwt(X ...

  5. 基于MATLAB的离散小波变换

    申明,本文非笔者原创,原文转载自:  基于Matlab的离散小波变换         http://blog.sina.com.cn/s/blog_725866260100ryh3.html 简介 在 ...

  6. 【DWT笔记】傅里叶变换与小波变换

    [DWT笔记]傅里叶变换与小波变换 一.前言 我们经常接触到的信号,正弦信号,余弦信号,甚至是复杂的心电图.脑电图.地震波信号都是时域上的信号,我们也成为原始信号,但是通常情况下,我们在原始信号中得到 ...

  7. 小波包分解 仿真 matlab

    clc;close all;clear;fs = 100000;t = 1: 100;x = sin(2*pi*4000* t/fs) + sin(2*pi*40000*t/fs); %db8[Lo_ ...

  8. paper 41:正交变换

    正交变换是保持图形形状和大小不变的几何变换,包含旋转,轴对称及上述变换的复合. 几何意义 正交变换是保持图形形状和大小不变的几何变换,包含旋转,轴对称及上述变换的复合. 代数定义 欧几里得空间V的线性 ...

  9. dennis gabor 从傅里叶(Fourier)变换到伽柏(Gabor)变换再到小波(Wavelet)变换(转载)

    dennis gabor 题目:从傅里叶(Fourier)变换到伽柏(Gabor)变换再到小波(Wavelet)变换 本文是边学习边总结和摘抄各参考文献内容而成的,是一篇综述性入门文档,重点在于梳理傅 ...

随机推荐

  1. spring boot 监控与管理(actuator)

    Spring POMs 中提供了一个特殊的依赖模块,即spring-boot-starter-actuator,我们只需要在我们的POM中添加依赖即可 <!-- 监控 管理 --> < ...

  2. Backbone源码风格

         代码风格: 一.自执行匿名函数创建执行环境 var root = this; root保存全局执行环境的指针.浏览器端为window对象 二.依赖库 (1).underscore 如果bac ...

  3. httpclient通过post提交到webapi

    var client = new HttpClient(); var url = BASConfig.Instance.SiteSettingsModule.SyncWorkLogAppUrl; va ...

  4. oracle数据库,表被锁,如何解锁?

    4.批量解锁declare cursor mycur isselect b.sid,b.serial# from v$locked_object a,v$session bwhere a.sessio ...

  5. 使用CMake生成VS2010项目查看OpenCV源代码

    近期项目需要用到OpenCV中的几个函数,但其函数无法全部实现自己需要的功能,故而需要改进部分函数,为安全及效率起见,想参考OpenCV的源码来改进,这样节省时间的同时亦可提供代码的鲁棒性和通用性.那 ...

  6. 《译》准备做一些 AR/增强现实的 翻译

    中文这方面资料实在少之又少. 准备做一些这方面翻译,关注于Vuforia, Unity3d, Hololens等方面. 如有问题.建议,随时联系.Fell free ton contact me.

  7. 转 zigbee学习笔记---Channel、PANID、发射功率及其它参数

    现对z-stack里几个网络参数的设置以及如何获取总结一下.信道配置:Zigbee在3个频段定义了27个物理信道:868MHz频段中定义了1个20Kb/s信道,915MHz频段中定义了10个40Kb/ ...

  8. java Vamei快速教程04 封装和接口

    作者:Vamei 出处:http://www.cnblogs.com/vamei 欢迎转载,也请保留这段声明.谢谢! 总结之前的内容,对象(object)指代某一事物,类(class)指代象的类型.对 ...

  9. linux 命令——49 at (转)

    在windows系统中,windows提供了计划任务这一功能,在控制面板 -> 性能与维护 -> 任务计划, 它的功能就是安排自动运行的任务. 通过'添加任务计划'的一步步引导,则可建立一 ...

  10. ABAP Netweaver和Cloud Foundry上的环境变量Environment Variable

    Netweaver 更准确的说应该是系统变量:结构体sy 设一个断点,调试器里看这些字段的值就能知道每个字段是用来做什么的. sy-dbsys sy-sysid sy-opsys sy-saprl s ...