傅里叶变换和小波分析

2023-10-30

无论是傅立叶变换还是小波变换,其实质都是一样的,既:将信号在时间域和频率域之间相互转换,从看似复杂的数据中找出一些直观的信息,再对它进行分析。由于信号往往在频域有比在时域更加简单和直观的特性,所以,大部分信号分析的工作是在频域中进行的。音乐——其实就是时/频分析的一个极好例子,乐谱就是音乐在频域的信号分布,而音乐就是将乐谱变换到时域之后的函数。从音乐到乐谱,是一次傅立叶或小波变换;从乐谱到音乐,就是一次傅立叶或小波逆变换。

1、傅立叶变换

可以理解为:任意一条在实数域内有意义的曲线都可以分解为若干个正弦曲线的叠加。傅立叶变换与分形的原理其实是同源的,不要小瞧这个变换,它可能就是宇宙的一个最基本法则。换句通俗的话说就是:无论多么复杂的物质,都可以用几种简单的基本物质通过一定方式的组合来构成。小到分子原子,中到地形地貌,大到银河系宇宙,其实都是这样。

但是,傅立叶变换也有它的缺陷。由于正弦波是无限宽度的,这使得被分析的信号也需要具有从负无穷大到正无穷大都有意义的特性,所以傅立叶变换不能很好的处理一些局部信号。比如,一个在局部范围内有非0值而其余所有地方都等于0的函数,它的频谱会呈现出一幅相当混乱的状况。这时,频域的信号反而不如时域的直观,频谱分析变得很艰难。

2、小波分析

为了克服傅立叶变换的这些缺陷,数学家和工程师们已经开发出若干种使用有限宽度基函数进行变换的方法。这些基函数不仅在频率上而且在位置上是变化的,这些有限宽度的波被称为小波Wavelet)。基于它们的变换被称为小波变换Wavelet transforms)。

小波变换的概念是由法国从事石油信号处理的工程师J.Morlet1974年首先提出的,通过物理的直观和信号处理的实际需要经验的建立了反演公式,当时未能得到数学家的认可。正如1807年法国的热学工程师J.B.J.Fourier提出任一函数都能展开成三角函数的无穷级数的创新概念未能得到著名数学家J.L.LagrangeP.S.Laplace以及A.M.Legendre的认可一样。幸运的是,早在七十年代,A.Calderon表示定理的发现、Hardy空间的原子分解和无条件基的深入研究为小波变换的诞生做了理论上的准备,而且J.O.Stromberg还构造了历史上非常类似于现在的小波基;1986年著名数学家Y.Meyer偶然构造出一个真正的小波基,并与S.Mallat合作建立了构造小波基的同意方法棗多尺度分析之后,小波分析才开始蓬勃发展起来,其中比利时女数学家I.Daubechies撰写的《小波十讲(Ten Lectures on Wavelets》对小波的普及起了重要的推动作用。它与Fourier变换、窗口Fourier变换(Gabor变换)相比,这是一个时间和频率的局域变换,因而能有效的从信号中提取信息,通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析(Multiscale Analysis),解决了Fourier变换不能解决的许多困难问题,从而小波变化被誉为数学显微镜,它是调和分析发展史上里程碑式的进展。

小波分析的应用是与小波分析的理论研究紧密地结合在一起的。现在,它已经在科技信息产业领域取得了令人瞩目的成就。 电子信息技术是六大高新技术中重要的一个领域,它的重要方面是图像和信号处理。现今,信号处理已经成为当代科学技术工作的重要部分,信号处理的目的就是:准确的分析、诊断、编码压缩和量化、快速传递或存储、精确地重构(或恢复)。从数学的角度来看,信号与图像处理可以统一看作是信号处理(图像可以看作是二维信号),在小波分析地许多分析的许多应用中,都可以归结为信号处理问题。现在,对于其性质随实践是稳定不变的信号,处理的理想工具仍然是傅立叶分析。但是在实际应用中的绝大多数信号是非稳定的,而特别适用于非稳定信号的工具就是小波分析

    事实上小波分析的应用领域十分广泛,它包括:数学领域的许多学科;信号分析、图像处理;量子力学、理论物理;军事电子对抗与武器的智能化;计算机分类与识别;音乐与语言的人工合成;医学成像与诊断;地震勘探数据处理;大型机械的故障诊断等方面;例如,在数学方面,它已用于数值分析、构造快速数值方法、曲线曲面构造、微分方程求解、控制论等。在信号分析方面的滤波、去噪声、压缩、传递等。在图象处理方面的图象压缩、分类、识别与诊断,去污等。在医学成像方面的减少B超、CT、核磁共振成像的时间,提高分辨率等。

    (1)小波分析用于信号与图像压缩是小波分析应用的一个重要方面。它的特点是压缩比高,压缩速度快,压缩后能保持信号与图像的特征不变,且在传递中可以抗干扰。基于小波分析的压缩方法很多,比较成功的有小波包最好基方法,小波域纹理模型方法,小波变换零树压缩,小波变换向量压缩等。

    (2)小波在信号分析中的应用也十分广泛。它可以用于边界的处理与滤波、时频分析、信噪分离与提取弱信号、求分形指数、信号的识别与诊断以及多尺度边缘检测等。

    (3)在工程技术等方面的应用。包括计算机视觉、计算机图形学、曲线设计、湍流、远程宇宙的研究与生物医学方面。

 

4Gabor FilterGabor滤波器)

     

Gabor小波是小波集中的一种,公式见(下面的链接)。从外形上看,Gabor小波被封装在一个Gaussion分布的形状中,而且它的积分为零。如果要将Gabor filter在三维中显示出来,应该是下面这个样子,左边是实部(偶函数),右边是虚部(奇函数)。http://cs.ccnu.edu.cn/wjylwjyl/kjqy/paper/Recognition.htm

通过改变k的相位和波长,可以得到一组不同的Gabor滤波器。出于速度和效果的综合考虑,一般使用8个方向和5种频率,这样,一共可以产生5x8=40个不同的Gabor滤波器。对于大小为128x128的图象,最小和最大频率的波长分别为16和4个象素。

4、Jet

将原始图像分别与每一个Gabor滤波器做卷积,会得到40个结果(注意,实部和虚部必须分开才能做卷积)。对于每一个输入的像素而言,则会产生40个输出的复数,我们将这40个复数按滤波器的顺序排列好,就是一个Jet.

如何做卷积?对于每个象素,将两个图像错开一定的距离放在一起,将重叠的像素相乘后累加。当然别忘了还有个经典公式:时域卷积 = 频域相乘。只要事先将图象变换到频域,相乘后,再反变换回时域即可。好在老外早在几十年前就发明了快速傅立叶变换算法(FFT),将原本是N*N的计算量减少到了N/2*log2N,使得信号在频域和时域中的转换变得十分快速。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

傅里叶变换和小波分析 的相关文章

  • sqli-labs靶场Less-15

    1 访问首页 Less 15 index php 这里的传参点是表单中的uname passwd 判断是否是数字形传参 判断是否有单引号闭合 判断是否是双引号闭合 判断是否单引号 括号闭合 判断是否是双引号 括号闭合 如果以上都不适 哪以上

随机推荐

  • type object ‘Path‘ has no attribute ‘home‘

    报错 根据堆栈定位 是matplotlib init py文件里面562行位置报错的问题 讲道理 应该没有问题来着 进入到matplotlib init py文件 然后悬停Path和home 发现 emmm 文件居然不一样 好的 在我尝试更
  • OpenCV缺陷检测

    新手上路 欢迎高手指点 原图 OpenCV缺陷检测代码 An highlighted block include head h int main Mat img imread image 1204 bmp assert im data if
  • python是面向对象还是面向过程的

    Python虽然是解释型语言 但从设计之初就已经是一门面向对象的语言 对于Python来说一切皆为对象 正因为如此 在Python中创建一个类和对象是很容易的 当然如果习惯面向过程或者函数的写法也是可以的 Python并不做硬性的限制 Py
  • python连接sqlsever_使用Python连接到Microsoft SQL Server

    我试图通过python连接到SQL以在Microsoft SQL服务器上的某些SQL数据库上运行一些查询 从我的在线研究和在这个论坛上 最有前途的图书馆似乎是pyodbc 所以我做了以下代码 import pyodbc conn pyodb
  • 【C#学习笔记】Hello World

    using System namespace ConsoleApplication class Program static void Main string args Console Write Hello World Console R
  • Class.forName()用法详解

    一 CLASS类概念 Class也是一个Java类 保存的是与之对应Java类的meta信息 元信息 用来描述这个类的结构 比如描述一个类有哪些成员 有哪些方法等 一般在反射中使用 详细解释 Java 源程序 java 文件 在经过 Jav
  • mysql 查询不出结果_mysql 执行查询SQL 一直执行不出结果

    今天执行一个mysql 语句 一直在 执行 执行了5分钟了 还是没有出来结果 每个组织下包括 同级或者下级的 注册店铺数 查询 历史每一天的每个组织下的当天存在的门店数 select dt time startDate o brand id
  • Spring Cloud 序列化和反序列化过程定制(Jackson)

    现在都是基于Spring Cloud Feign进行微服务的调用 并且序列化的过程都封装完成的 只是自己可以定制序列化的方式 但是为了调用的时候能方便的找到问题所在等 基本都会使用json Jackson等 方式的序列化 虽然性能比较差 但
  • 学习笔记-弗洛伊德算法

    弗洛伊德算法解决最短路径问题 弗洛伊德算法和迪杰斯特拉算法都可以求解最短路径的问题 但区别是迪杰斯特拉会求出某个顶点到其他顶点的最短路径 而弗洛伊德算法会求出各个顶点到各个顶点的最短路径 最终结果采用一个二维表表示 而迪杰斯特拉算法只需要用
  • 文件读写的并发操作分析

    前言 涉及到多进程 线程间对文件的并发读写 首先说明一下比较常见的多进程读写方法是在保证数据不混乱的前提下 让某一个进程专门负责写该文件 其它进程负责往该进程发消息 通常在日志系统中 开启一个专门的进程 线程进行文件的写操作 其他进程 线程
  • Windows下用pip安装lib时报错的简单解决思路

    今天在Windows下用pip尝试安装Python爬虫库Scrapy 但是安装的时候报错 不过看命令行里也能运行scrapy命令 以为没问题 结果写了个简单爬虫 不停地报各种错 怀疑可能是这个系统中的pip版本问题 想了想 可以先重装pip
  • 操作系统中的硬链接和符号链接的不同

    硬链接和符号链接 这里没有介绍它们的功能和优点缺点 太多博客已经写过了 主要是解释了它们的原理 出现硬链接和符号链接的原因 我们需要共享文件 如果一个共享文件同时出现在属于不同用户的不同目录下 工作起来就会很方便 但是 怎么解决不同用户下对
  • 在未来的十年无疑将是刷脸支付的世界

    在未来的十年 交易场景的闭环 无疑将是刷脸支付的世界 最近比较热的词就是 刷脸支付 相对扫码支付 刷脸支付的方式更加安全 更加便捷 支付宝和微信同时推出刷脸支付 微信叫 青蛙 用户只需站在屏幕前 输入支付金额即可完成刷脸支付 按住屏幕下方的
  • vue项目性能优化详解汇总

    提起性能优化 很多人眼前浮现的面试经验是不是历历在目呢 反正 性能优化在我看来他永远是前端领域的热度之王 先说一下性能优化的方案 一 基础优化 代码以及编码规范 1 v if 和 v show 区分使用场景 v if false时不渲染DO
  • python输入姓名_Python基础篇--输入与输出

    在任何语言中 输入和输出都是代码最基础的开始 so 先来聊一聊输入和输出 输出 在python中 我们一般用print 输出 在括号里输入你想输出的信息 用引号包裹起来 单双三都可以 例如我们来输出一个 hello python gt gt
  • 【endnote】利用endnote批量修改参考文献格式

    1 情况描述 本来是在word中用交叉引用插入了参考文献 没有用endnote 但是格式没有统一 现在需要把全部参考文献统一成gb7714的格式 大概三百多篇 2 方法 1 在谷歌学术中挨个搜索参考文献 点击引用 点击 endnote 下载
  • Ubuntu14.04桥接网络设置与SSH登陆

    操作系统 Unbuntu14 04 虚拟机 VMware10 一 网络设置 1 设置vmware Bridge Protocol 本地链接 gt 属性 gt vmware Bridge Protocol打钩 2 然后主机设定静态ip 如果已
  • Matlab学习6-图像处理之直方图处理、灰度变换

    1 直方图均衡化 代码 直方图均衡化 img imread img rice png 显示 subplot 2 2 1 imshow img xlabel 原图 subplot 2 2 2 imhist img xlabel 原图的灰度直方
  • qt中漂亮的几款 qss 样式

    Qt中漂亮的几款QSS Shared QStackedWidget QLabel QPushButton QRadioButton QCheckBox QGroupBox QStatusBar QToolButton QComboBox Q
  • 傅里叶变换和小波分析

    无论是傅立叶变换还是小波变换 其实质都是一样的 既 将信号在时间域和频率域之间相互转换 从看似复杂的数据中找出一些直观的信息 再对它进行分析 由于信号往往在频域有比在时域更加简单和直观的特性 所以 大部分信号分析的工作是在频域中进行的 音乐