scipy连续小波变换

连续小波变换,使用小波函数对数据执行连续小波变换。CWT使用由宽度参数和长度参数表征的小波函数执行与数据的卷积。允许小波函数是复数形式。

接口:

scipy.signal.cwt(data, wavelet, widths, dtype=None, **kwargs)

“”"
Continuous wavelet transform.

Performs a continuous wavelet transform on data, using the wavelet function. A CWT performs a convolution with data using the wavelet function, which is characterized by a width parameter and length parameter. The wavelet function is allowed to be complex.

Parameters
data(N,) ndarray
data on which to perform the transform.

waveletfunction
Wavelet function, which should take 2 arguments. The first argument is the number of points that the returned vector will have (len(wavelet(length,width)) == length). The second is a width parameter, defining the size of the wavelet (e.g. standard deviation of a gaussian). See ricker, which satisfies these requirements.
(小波函数,它应该带2个参数。第一个参数是返回的向量将具有的点数(len(小波(长度,宽度))==length)。第二个是宽度参数,定义了小波的大小(例如高斯的标准差)。参见ricker,它满足了这些要求。)

widths(M,) sequence
Widths to use for transform.

dtype:data-type, optional
The desired data type of output. Defaults to float64 if the output of wavelet is real and complex128 if it is complex.
所需的输出数据类型。如果小波的输出为实值,则默认为float64,如果它是复数的,则复数128。

kwargs
Keyword arguments passed to wavelet function.

Returns
cwt: (M, N) ndarray
Will have shape of (len(widths), len(data)).

对于非对称复值小波,输入信号与小波数据的时反复共轭卷积[1].
“”"

代码

from scipy import signal
import matplotlib.pyplot as plt
import numpy as np

t = np.linspace(-1, 1, 200, endpoint=False)
sig = np.cos(2 * np.pi * 7 * t) + signal.gausspulse(t - 0.4, fc=2)
widths = np.arange(1, 31)
cwtmatr = signal.cwt(sig, signal.ricker, widths)
plt.imshow(cwtmatr, extent=[-1, 1, 31, 1], cmap=‘PRGn’, aspect=‘auto’,
vmax=abs(cwtmatr).max(), vmin=-abs(cwtmatr).max())
plt.show()

引用

[^1] : S. Mallat, “A Wavelet Tour of Signal Processing (3rd Edition)”, Academic Press, 2009.
et Tour of Signal Processing (3rd Edition)”, Academic Press, 2009.
[^2] : scipy.signal.cwt — SciPy v1.9.3 Manual

scipy连续小波变换相关推荐

  1. matlab 小波变换_matlab小波工具箱实例(二):时频分析和连续小波变换

    本文讲解matlab小波工具箱实例(二):时频分析和连续小波变换.目录如下: 链接:https://www.mathworks.com/help/wavelet/ug/time-frequency-a ...

  2. 小波的秘密5_多分辨率分析和连续小波变换2

    1.如何理解时间分辨率和频率分辨率? In this section we will take a closer look at the resolution properties of the wa ...

  3. matlab 小波变换_连续小波变换实现方法的总结及其程序详解

    在帖子"给大家分享我自己编的程序-连续小波变换" 中,pengzk版友给出了morlet小波变换的源代码,但其中的许多参数和语句意义不够明确,这就给一些希望了解连续小波变换实现方法 ...

  4. 小波分析实验: 实验1 连续小波变换

    实验目的: 在理解连续小波变换原理的基础上,通过编程实现对一维信号进行连续小波变换,(实验中采用的是墨西哥帽小波),从而对连续小波变换增加了理性和感性的认识,并能提高编程能力,为今后的学习和工作奠定基 ...

  5. 基于连续小波变换的厄尔尼诺海平面周期变化数据集分析

    基于连续小波变换的厄尔尼诺海平面周期变化数据集分析 1. 厄尔尼诺海平面周期变化数据集介绍 2. 基于连续小波变换的信号分析 2.1 原始信号读取可视化 2.2 傅里叶变换结果可视化 2.3 小波变换 ...

  6. 【小波分析】连续小波变换的概念swt,cwt,dwt

    连续小波变换的概念.操作.及时间尺度图的显示 最近很多网友问到关于连续小波变换的诸多问题,我用了点时间,写了个底层程序,提供给大家参考. 1.连续小波的概念.就是把一个可以称作小波的函数(从负无穷到正 ...

  7. 使用连续小波变换(cwt)对时间轴上频率发生变化的非平稳信号进行处理

    使用连续小波变换对时间轴上频率发生变化的非平稳信号进行处理 傅里叶变换对时间轴上的平稳信号处理有着很大的优势,但是对非平稳信号就无能为力了. 小波变换是基于傅里叶变换产生的可以用于处理非平稳信号的信号 ...

  8. 小波变换学习笔记(2):连续信号的连续小波变换与离散小波变换

    前言 上次说到小波变换的知识体系,这篇博客就主要说小波变换里的连续信号的连续小波变换与离散小波变换. 连续信号的连续小波变换 话不多说,我们先放公式,如果你是第一次接触小波,你可会有点懵,但是不要怕, ...

  9. 制作凯斯西储轴承故障信号(CWRU)的连续小波变换(CWT)时频数据集

    制作CWRU的连续小波变换的时频数据集 1. 遍历指定文件夹下的mat文件 2. 基于子函数的CWRU信号分割 3. CWT时频提取 3.1 主程序 3.2 批量CWT绘图 4. 提取结果的可视化 5 ...

最新文章

  1. Java基础学习总结(14)——Java对象的序列化和反序列化
  2. vlookup 2张表 显示na_Vlookup函数的这7个应用技巧都不掌握,那就真的Out了
  3. ImportError: No module named protobuf
  4. phpmyadmin #1045 - Access denied for user 'root'@'localhost' (using password: NO)
  5. 快速获取青年大学习完成截图
  6. django中的认证登陆与用户的创建
  7. S/4HANA extension field search的SQL语句是在什么地方生成的
  8. ai无法启动产品_启动AI启动的三个关键教训
  9. 智能包装结构,提高可测性
  10. php console postman,postman内置脚本说明
  11. Mongo 与 SQL 操作对应介绍
  12. java 图片地址,java中图片地址存入数据库
  13. python能制作ppt动画效果吗_你听说过Python可以做动画吗
  14. Apache+tomcat配置集群(1)
  15. directx9 截屏
  16. 如何提高你的工作效率
  17. python在abaqus中的应用光盘文件下载_Python语言在Abaqus中的应用DVD光盘
  18. android从应用到驱动之—camera(2)---cameraHAL的实现
  19. Android display
  20. SVM支持向量机的应用

热门文章

  1. 采用thinkpphp5 采集完成的优惠券采集网站,采集有淘宝,好单库,大淘客网站商品数据
  2. 分享一款精美的后台模板
  3. 曾轶可哥哥太吓人了~~~
  4. iptables实现网络限制下ntp自定义端口同步时间
  5. 【算法】数字读法转汉字
  6. 第一个C++语言,用的VS2019,感觉不错
  7. 用Tkinter做一个每日英语的程序
  8. qt5报错Installed Qt WebEngine locales directory not found at location
  9. 2017c语言 形成性考核,(2017年电大)《c语言程序设计》形成性考核作业()解答.doc...
  10. nested exception is java.lang.ClassCastException: java.lang String cannot be cast to java.util.Map