clear all
s=('matlab123.mat')
alpha=2000;
tau=0;
k=3;
DC=0;
init=1;
tol=1e-7;
[u, u_hat, omega] = VMD(s, alpha, tau, K, DC, init, tol);
%--------------- Visualization
% For convenience here: Order omegas incclcreasingly and reindex u/u_hat
[~, sortIndex] = sort(omega(end,:));
omega = omega(:,sortIndex);
u_hat = u_hat(:,sortIndex);
u = u(sortIndex,:);
linestyles = {'b', 'g', 'm', 'c', 'c', 'r', 'k'};
figure('Name', 'Composite input signal' );
plot(t,f, 'k');
set(gca, 'XLim', [0 1]);
for sub = 1:length(fsub)
figure('Name', ['Input signal component ' num2str(sub)] );
plot(t,fsub{sub}, 'k');