实测信号VMD分解

% 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