一种改进的经验小波变换方法(Python环境)

152 篇文章 33 订阅
125 篇文章 46 订阅

经验小波变换EWT是Gilles基于小波分析理论提出的一种新的自适应信号分解方法,该方法主要分为三个步骤:1.根据傅里叶谱的特性自适应划分频谱,获得一组边界;2.根据边界序列和Meyer小波构造滤波器组;3.滤波重构,获得若干个具有紧支撑傅里叶谱的经验模态分量。相比于EMD与EEMD,EWT在故障信息的提取与计算效率方面更有优势。EWT以其计算效率高、严谨的数学推导过程和较小的模态混叠等优点在许多领域得到广泛应用,如滚动轴承、行星齿轮箱、风机轴承的故障诊断、语音信号处理、太阳能电池阵列的冲击振动试验等。

随着经验小波变换在各领域的不断应用,专家学者们逐渐发现EWT的不足之处。在Gilles提出的经验小波变换方法中,需要预先设置大量的参数,从而造成很多时候难以获得合理的边界,同时由于该方法对噪声的敏感性,可能会获得很多无效分量,专家学者对此进行了改进,涌现出了大量经验小波变换的改进方法。目前对EWT的改进主要集中在四个方面:第一个改进角度是优化经验模态,例如模态的分类和利用某些准则来合并经验模态。第二个改进角度是以新的谱替换传统的傅立叶频谱。频谱划分是经验小波变换中的关键步骤之一,但Gilles提出的尺度空间表示等频谱划分方式对噪声比较敏感,容易产生大量的无效边界。第三个改进角度是优化EWT的设定参数,比如利用Otsu方法和K-Mean方法建立一组带通滤波器以自适应分解信号。第四个改进角度是优化频谱划分方式,比如利用顺序统计滤波器估计谱的上包络函数。

鉴于此,采用一种级数展开代替傅里叶频谱来实现频谱划分,Python代码编写,结果如下:

图片

图片

图片

图片

图片

图片

图片

图片

图片

代码可通过知乎学术咨询获得,知乎ID:哥廷根数学学派。

工学博士,担任《Mechanical System and Signal Processing》《中国电机工程学报》《控制与决策》等期刊审稿专家,擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。

  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 这个工具箱组织如下: EWT ?| ?| -1D:1D EWT功能 ?| -2D:2D EWT功能 ?| | - 小波:经验曲线变换 ?| | -Littlewood-Paley:经验的Littlewood-Paley小波变换 ?| | -Ridgelet:经验Ridgelet变换 ?| | - 张量:经验张量小波变换 ? - 边界:用于执行傅里叶支持的函数 ?| | - LocalMaxima:根据当地最大值,中途或当地最小值执行检测的功能 ?| | -MorphoMath:执行形态学操作符对谱进行预处理的功能 ?| | -PowerLaw:通过去除其幂律近似来预处理谱 ?| | -ScaleSpace:基于尺度空间方法执行检测的函数 ?| - 文档:工具箱文档 ?| -Tests ?| | -1D:对几个1D信号执行基本测试的功能 ?| | -2D:用于在不同图像上执行几个2D变换的基本测试的功能 ?| -utilities ?| | -1D:在1D情况下绘制结果的有用函数(时频平面,分量,边界) ?| | -2D:用于在2D情况下绘制结果的有用函数(不同类型的组件,2D边界,...) 如果要运行所有功能,则需要在计算机上正确安装以下Matlab工具箱: - Flandrin的EMD工具箱(在一维变换中需要执行希尔伯特变换并且可视化时频平面) 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱(2D变换除了基于张量变换之外) 可从http://www.cs.technion.ac.il/~elad/software/获得 已经包含了
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

哥廷根数学学派

码字不易,且行且珍惜

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值