变分模式分解(vmd)算法
vmdp的Python项目详细描述
变分模式分解python包
计算信号变分模式分解的函数(Dragomiretskiy和Zosso,2014年)
原始vmd纸张:
Dragomiretskiy,K.和Zosso,D.(2014)“变分模态分解”,
《信号处理ieee汇刊》,62(3),531–544页。doi:10.1109/tsp.2013.2288675。
原始Matlab代码:https://www.mathworks.com/matlabcentral/fileexchange/44765-variational-mode-decomposition
安装
- 从https://github.com/vrcarva/vmdpy下载项目,然后从项目文件夹运行“python setup.py install”
或
- pip安装vmdpy
引用和联系
如果您觉得这个包有用,我们恳请您在工作中引用它。
Vinicius Carvalho(2019-),Python中的变分模式分解
一篇论文将很快提交并链接到这里。
联系人:vrcarva@ufmg.br
卡瓦略市
巴西贝洛奥里藏特米纳斯格拉斯联邦大学
神经科学学院-米纳斯杰拉斯联邦大学
示例脚本
#%% Simple example importnumpyasnpimportmatplotlib.pyplotaspltfromvmdpyimportVMD#. Time Domain 0 to T T=1000fs=1/Tt=np.arange(1,T+1)/Tfreqs=2*np.pi*(t-0.5-fs)/(fs)#. center frequencies of components f_1=2f_2=24f_3=288#. modes v_1=(np.cos(2*np.pi*f_1*t))v_2=1/4*(np.cos(2*np.pi*f_2*t))v_3=1/16*(np.cos(2*np.pi*f_3*t))f=v_1+v_2+v_3+0.1*np.random.randn(v_1.size)#. some sample parameters for VMD alpha=2000# moderate bandwidth constraint tau=0.# noise-tolerance (no strict fidelity enforcement) K=3# 3 modes DC=0# no DC part imposed init=1# initialize omegas uniformly tol=1e-7#. Run actual VMD code u,u_hat,omega=VMD(f,alpha,tau,K,DC,init,tol)