%%?-------?Preparation
clear?all;
close?all;
clc;
%%?Time?Domain?0?to?T
T?=?1000;
fs?=?1/T;
t?=?（1:T）/T;
freqs?=?2*pi*（t-0.5-1/T）/（fs）;
%%?center?frequencies?of?components
f_1?=?2;
f_2?=?24;
f_3?=?288;
%%?modes
v_1?=?（cos（2*pi*f_1*t））;
v_2?=?1/4*（cos（2*pi*f_2*t））;
v_3?=?1/16*（cos（2*pi*f_3*t））;
%%??for?visualization?purposes
wsub{1}?=?2*pi*f_1;
wsub{2}?=?2*pi*f_2;
wsub{3}?=?2*pi*f_3;
%%??composite?signal?including?noise
f?=?v_1?+?v_2?+?v_3?+?0.1*randn（size（v_1））;
%%?some?sample?parameters?for?VMD
alpha?=?2000;????????%?moderate?bandwidth?constraint
tau?=?0;????????????%?noise-tolerance?（no?strict?fidelity?enforcement）
K?=?4;??????????????%?4?modes
DC?=?0;?????????????%?no?DC?part?imposed
init?=?1;???????????%?initialize?omegas?uniformly
tol?=?1e-7;
%%?--------