时间序列信号处理（一）-----变分模态分解（VMD)

1.简述变分模态分解

变分模态分解--vmd，适用于非线性时间序列信号，主要是利用求解变分问题的思想去对信号进行提取，在不丢失原始信号特征的情况下，把一个原始信号分解成多个不同中心频率的信号，即不在同一个调制信号内。

2.以轴承信号为例

安装vmd库，直接pip install vmdpy 即可

import matplotlib.pyplot as plt
import numpy as np
from vmdpy import VMD

读取信号

ball_18_0 = np.loadtxt('E:/12k1/0HP/inner18.txt')
data = ball_18_0[2048*2:2048*3]

设置参数，但分解层数K和二次惩罚系数α需要仔细考虑设定，特别是K值，对分解效果影响巨大。

# 参数设置
alpha = 2000      # moderate bandwidth constraint
tau = 0.            # noise-tolerance (no strict fidelity enforcement)
K = 5            # 3 modes
DC = 0             # no DC part imposed
init = 1           # initialize omegas uniformly
tol = 1e-7u, u_hat, omega = VMD(data, alpha, tau, K, DC, init, tol)

绘制分解信号图

plt.figure(figsize=(10, 8))
for i in range(K):plt.subplot(K+1, 1, 1)plt.plot(data)plt.title("outer")plt.subplot(K+1, 1, i+2)plt.plot(u[i, :], linewidth=0.2, c='r')plt.ylabel('u{}'.format(i + 1))
plt.tight_layout()
plt.show()

最后绘出其频谱图

Fs = 12000
Ts = 1.0/Fs
t = np.arange(N)
k = np.arange(N)
T = N/Fs
frq = k/T
frq1 = frq[range(int(N/2))]lt.figure(figsize=(10, 8))
for i in range(K):plt.subplot(K + 1, 1, 1)data_f = abs(np.fft.fft(data)) / Ndata_f1 = data_f[range(int(N / 2))]plt.plot(frq1, data_f1)plt.title("outer")plt.subplot(K+1, 1, i + 2)data_f2 = abs(np.fft.fft(u[i, :])) / Ndata_f3 = data_f2[range(int(N / 2))]plt.plot(frq1, data_f3, 'red')plt.xlabel('pinlv(hz)')plt.ylabel('u{}'.format(i + 1))
plt.tight_layout()
plt.show()

由上述结果可知，分解出来的信号基本上不在同一个频率段，分解效果很好。

时间序列信号处理（一）-----变分模态分解（VMD)相关推荐

鲸鱼算法优化变分模态分解(VMD)包络熵和参数的特征提取及MATLAB代码实现
目录 1 简介 2 变分模态分解VMD原理 3 鲸鱼优化算法优化VMD原理 3.1. 鲸鱼优化算法优化VMD原理及流程 3.2. 特征提取流程 4 优化效果 4.1. VMD各分量信号时域图 4.2. ...
vmd变分模态分解程序matlab论坛_博士兼职辅导员论坛分享会第三期
新一期经验分享报告会又来了!本次报告满满干货,快来看看都有些什么内容吧? 报告题目齿轮箱关键零部件复合故障特征提取方法研究报告摘要: 基于振动信号的复合故障特征提取技术一直以来都是旋转机械故障诊断 ...
变分模态分解 python_浅谈VMD(变分模态分解)
学号:19011210554 姓名:袁博 [嵌牛导读]:好多人看着VMD看博客最想知道的就是这东西的应用和大概步骤原理,而具体原理算法不太感兴趣,而且也不太容易看懂.本文既然是浅谈,就讲解一下VM ...
matlab中使用VMD(变分模态分解)
最近我们被客户要求撰写关于VMD(变分模态分解)的研究报告,包括一些图形和统计输出. 拨号音信号的变模分解创建一个以4 kHz采样的信号,类似于拨打数字电话的所有键.将信号另存为MATLAB®时间数 ...
利用智能算法优化参数的自适应变分模态分解，VMD实现混合储能系统的分频
关键词:混合储能,VMD,麻雀搜索算法,遗传算法,混合储能容量配置优化,混合储能功率分配,利用智能算法优化参数的自适应变分模态分解,VMD实现混合储能系统的分频,高频分配给超级电容器,低频分配给蓄电池 ...
分解得到的时频域特征_【推荐文章】基于变分模态分解和广义Warblet变换的齿轮故障诊断...
<机械传动>2018年第42卷第7期文章编号:1004-2539(2018)07-0157-05 DOI:10.16578/j.issn.1004.2539.2018.07.0 ...
变分模态分解_Android小部件示例中的模态对话框（弹出）
变分模态分解在此示例中,我们将看到如何在主屏幕中创建一个可以打开弹出对话框的Android小部件. 如您所知,Android Widgets是小型应用程序,基本上可以做两件事. 按下时启动新的活动, ...
linux环境vmd下载,VMD Linux版下载|VMD(变分模态分解程序) V1.9.3 Linux版下载_当下软件园_软件下载...
VMD Linux版是款适用于Linux操作系统的分子运动绘图分析模拟软件.它可以真实的模拟分析的运动场景,提供可视化的图形界面,帮助用户更好的理解分析变化,提高研究效率,操作简单,方便快捷,非常好用 ...
VMD-SSA-LSTM基于变分模态分解和麻雀算法优化长短期记忆网络的时间序列预测MATLAB代码(含LSTM、VMD-LSTM、VMD-SSA-LSTM等模型的对比)。
clc; clear all close all %% VMD-SSA-LSTM预测 tic load vmd_data.mat load lstm.mat disp('--------------- ...
【VMD-SSA-LSSVM】基于变分模态分解与麻雀优化Lssvm的负荷预测【多变量】（Matlab代码实现）

时间序列信号处理（一）-----变分模态分解（VMD)

1.简述变分模态分解

2.以轴承信号为例

时间序列信号处理（一）-----变分模态分解（VMD)相关推荐

最新文章

热门文章