driver_2d_example

cleanup

0001 %cleanup
0002 clearvars
0003 close all
0004 
0005 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0006 %Parameters
0007 param.n_grid_1    = 256;  %Grid size in dimension 1
0008 param.n_grid_2    = 256;  %Grid size in dimension 2
0009 param.nimfs       = 3;    %Maximum number of IMFs that can be stored
0010 param.type        = 6;    %type of window size
0011 param.tol         = 0.01; %sifting tolerance
0012 param.plot        = 'on'; %toggle plotting
0013 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0014 
0015 X = linspace(0,30,param.n_grid_1);
0016 Y = linspace(0,30,param.n_grid_2);
0017 
0018 [x,y] = meshgrid(Y,X); 
0019 
0020 uc1 = (sin(4*x)+sin(4*y));  uc2 = (sin(x)+sin(y)); ru  = (cos(x/24)+cos(y/24));
0021 u   =  uc1 + uc2 + ru; %Signal 'u'
0022 
0023 vc1 = (sin(x)+sin(y));  vc2 = (sin(3.5*x)+sin(3.5*y)); rv  = (cos(x/24)+cos(y/24));
0024 v   =  vc1 ; %Signal 'v'
0025 
0026 wc1 = (cos(x/24)+cos(y/24));  wc2 = (sin(x/5)+sin(y/5)); rw = (sin(x)+sin(y)) ;
0027 w   =  wc2 ;%+ wc2 + rw; %Signal 'w'
0028 
0029 Results = EMD2D3V(u,v,w,param);
0030 
0031 
0032