%?test-script?for?VMD
%?authors:?Dominique?Zosso?and?Konstantin?Dragomiretskiy
%?zosso@math.ucla.edu?---?http://www.math.ucla.edu/~zosso
%?Initial?release?2013-12-12?（c）?2013
%
%?When?using?this?code?please?do?cite?our?paper:
%?-----------------------------------------------
%?K.?Dragomiretskiy?D.?Zosso?Variational?Mode?Decomposition?IEEE?Trans.
%?on?Signal?Processing?（in?press）
%?please?check?here?for?update?reference:?
%??????????http://dx.doi.org/10.1109/TSP.2013.2288675

%---------------?Preparation
clear?all;
close?all;
clc;

%?Time?Domain?0?to?TT?=1000;fs?=?1/T;t?=?（1:T）/T;freqs?=?2*pi*（t-0.5-1/T）/（fs）;
fs=20480;
T=600;
nn=fs*T;
t=1/fs:1/fs:T;
freqs?=?（（1/T:1/T:fs）-fs/2-1/T）;
%?center?frequencies?of?components
%f1=load（‘LIFEa.mat‘）;
%dlmwrite（‘f2.txt‘f1‘delimiter‘?‘\t‘‘precision‘‘%6.2f‘）;
load（‘LIFEa.mat‘）
a=reshape（data20480*9841）;
f=a‘;
f_hat?=fftshift（fft（f））;grid?on;

%?some?sample?parameters?for?VMD
alpha?=?500;????????%?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?