vmd算法代码

 概述
变分模态分解由Konstantin Dragomiretskiy于2014年提出，可以很好抑制EMD方法的模态混叠现象（通过控制带宽来避免混叠现象）。

与EMD原理不同，VMD分解方式是利用迭代搜索变分模型最优解来确定每个分解的分量中心频率及带宽，属于完全非递归模型，该模型寻找模态分量的集合及其各自的中心频率，而每个模态在解调成基带之后是平滑的。

Konstantin Dragomiretskiy通过实验结果证明：对于采样和噪声方面，该方法更具有鲁棒性。

模态变分分解是一种自适应、完全非递归的模态变分和信号处理的方法。该技术具有可以确定模态分解个数的优点，其自适应性表现在根据实际情况确定所给序列的模态分解个数.

随后的搜索和求解过程中可以自适应地匹配每种模态的最佳中心频率和有限带宽，并且可以实现固有模态分量（IMF）的有效分离、信号的频域划分、进而得到给定信号的有效分解成分，最终获得变分问题的最优解。

它克服了EMD方法存在端点效应和模态分量混叠的问题，并且具有更坚实的数学理论基础，可以降低复杂度高和非线性强的时间序列非平稳性，分解获得包含多个不同频率尺度且相对平稳的子序列，适用于非平稳性的序列，VMD的核心思想是构建和求解变分问题。

构造变分问题
假设原始信号f被分解为k个分量，保证分解序列为具有中心频率的有限带宽的模态分量，同时各模态的估计带宽之和最小，约束条件为所有模态之和与原始信号 相等，则VMD约束变分模型如下：

其中，
为各模态函数，
为各模态中心频率。

变分求解问题(引入拉格朗日)
解决上述的约束最优化问题，将约束变分问题转变为非约束变分问题，利用二次惩罚项和拉格朗日乘子法的优势，引入了增广拉格朗日(Lagrangian)函数，如式所示：

其中，α \alphaα是罚参数，λ \lambdaλ是Lagrangian乘子。

对于所有ω \omegaω>=0，更新泛函

更新泛函ω \omegaωk：

所有ω \omegaω>=0，进行双重提升：

γ \gammaγ表示噪声，当信号含有强噪声时，可以设定γ \gammaγ=0达到更好的去噪效果。
不断重复更新泛函，知道满足下面的迭代约束条件：

对于所有