VMD变分模态分解算法

根据单个或多个宽频偶数长度电磁信号，根据VMD计算分量数据，展示含噪曲线、去噪曲线、分量时域、分量频谱、多测道曲线的图形绘制和数据生成。实现单条信号的独立运算，多条信号可选参数的独立运算并实现合并。 源码中包含一组测试数据，结果与matlab一致。

import numpy as np
import numpy.fft as nf
import matplotlib.pyplot as py
#import pylab as py
alpha = 2000        # moderate bandwidth constraint 惩罚因子
tau = 0             # noise-tolerance （no strict fidelity enforcement） 0 噪声容忍度
K =5                # VMD分解个数，这里需要人工设置
DC = 0              # no DC part imposed
init = 1            # initialize omegas uniformly中心频率初始化
tol = 1e-10
time=[1.00E-051.01E-051.02E-051.04E-051.05E-051.06E-051.07E-051.08E-051.10E-051.11E-051.12E-051.14E-051.15E-051.16E-051.18E-051.19E-051.20E-051.22E-051.23E-051.24E-051.26E-051.27E-051.29E-051.30E-051.32E-051.33E-051.35E-051.37E-051.38E-051.40E-051.41E-051.43E-051.45E-051.46E-051.48E-051.50E-051.51E-051.53E-051.55E-051.57E-051.59E-051.60E-051.62E-051.64E-051.66E-051.68E-051.70E-051.72E-051.74E-051.76E-051.78E-051.80E-051.82E-051.84E-051.86E-051.88E-051.91E-051.93E-051.95E-051.97E-052.00E-052.02E-052.04E-052.07E-052.09E-052.12E-052.14E-052.16E-052.19E-052.21E-052.24E-052.27E-052.29E-052.32E-052.35E-052.37E-052.40E-052.43E-052.46E-052.49E-052.51E-052.54E-052.57E-052.60E-052.63E-052.66E-052.69E-052.73E-052.76E-052.79E-052.82E-052.85E-052.89E-052.92E-052.95E-052.99E-053.02E-053.06E-053.09E-053.13E-053.17E-053.20E-053.24E-053.28E-053.32E-053.35E-053.39E-053.43E-053.47E-053.51E-053.55E-053.59E-053.64E-053.68E-053.72E-053.76E-053.81E-053.85E-053.90E-053.94E-053.99E-054.03E-054.08E-054.13E-054.17E-054.22E-054.27E-054.32E-054.37E-054.42E-054.47E-054.53E-054.58E-054.63E-054.68E-054.74E-054.79E-054.85E-054.91E-054.96E-055.02E-055.08E-055.14E-055.20E-055.26E-055.32E-055.38E-055.44E-055.50E-055.57E-055.63E-055.70E-055.76E-055.83E-055.90E-055.97E-056.04E-056.11E-056.18E-056.25E-056.32E-056.39E-056.47E-056.54E-056.62E-056.70E-056.77E-056.85E-056.93E-057.01E-057.09E-057.18E-057.26E-057.34E-057.43E-057.51E-057.60E-057.69E-057.78E-057.87E-057.96E-058.05E-058.15E-058.24E-058.34E-058.43E-058.53E-058.63E-058.73E-058.83E-058.93E-059.04E-059.14E-059.25E-059.35E-059.46E-059.57E-059.68E-059.79E-059.91E-050.0001002310.0001013930.0001025680.0001037570.0001049590.0001061760.0001074070.0001086520.0001099110.0001111850.0001124740.0001137770.0001150960.000116430.000117780.0001191450.0001205260.0001219230.0001233360.0001247660.0001262120.0001276750.0001291550.0001306520.0001321660.0001336980.0001352480.0001368160.0001384020.0001400060.0001416290.000143270.0001449310.0001466110.000148310.0001500290.0001517680.0001535280.0001553070.0001571070.0001589280.000160770.0001626340.0001645190.0001664260.0001683550.0001703070.0001722810.0001742770.0001762980.0001783410.0001804080.0001824990.0001846150.0001867550.0001889190.0001911090.0001933240.0001955650.0001978320.0002001250.0002024450.0002047910.0002071650.0002095660.0002119950.0002144530.00021