vmd去噪MATLAB代码

在emd基础上改进过的vmd去噪方法，可直接使用。有注释

clc
clear all
close all
tic
load 130.mat;%读取外圈故障数据，这个数据在VMD文件夹给你了
toc
tic
x=X130_DE_time（1:4096）;%从读取的数据中取X130_DE_time数据作为要分解的信号
fs=12000;%采样频率，这个是信号处理里面要设定信号采集时候的频率，是由采集信号传感器决定，你不用管
N=4096;%采样点数，这是采集点数，也就是信号的长短
n=0:N-1;
f=n*fs/N;
T=1000;
n=0:N-1;
t=n/fs;

% some sample parameters for VMD
%这个是我要优化的VMD信号处理的算法，包含6个参数，其中alpha和K是要优化的，其他都不用管，取默认值
alpha = 4000;        %这叫惩罚因子，可以取100~2000以内
tau = 0;            % noise-tolerance （no strict fidelity enforcement）
K = 5;              % 这个是信号分解的个数，设定多少就分解成多少个
DC = 0;             % no DC part imposed
init = 1;           % initialize omegas uniformly
tol = 1e-5;
%上面的参数设定好后再进行下面的信号分解，设定K=8就是把信号分解成8个，如果alpha设定不同，也会影响分解的8个分量的频率不一样，所以这两个参数都会影响分解结果
%我们的思路是当alpha取100~2000，K取1~10，为了减少计算量，把alpha取100200300，整百间隔的跳，也就是20个alpha值和10个K值组合

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2928192  2010-03-08 15:43  VMDtest\130.mat

     文件       2968  2018-03-06 12:53  VMDtest\main.m

     文件       4645  2017-03-06 21:32  VMDtest\VMD.m

     文件        852  2017-12-22 15:11  VMDtest\yyshang.m

     目录          0  2018-03-07 19:54  VMDtest

----------- ---------  ---------- -----  ----

              2936657                    5