driver_3d_example

cleanup

0001 %cleanup
0002 clearvars
0003 close all
0004 
0005 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0006 %Parameters
0007 param.n_grid_1    = 30;  %Grid size in dimension 1
0008 param.n_grid_2    = 30;  %Grid size in dimension 2
0009 param.n_grid_3    = 30;  %Grid size in dimension 3
0010 param.nimfs       = 3;    %Maximum number of IMFs that can be stored
0011 param.type        = 5;    %type of window size
0012 param.tol         = 0.05; %sifting tolerance
0013 param.nslice      = 10;   %number of slices in volume plot
0014 param.plot        = 'on'; %toggle plotting
0015 
0016 %USE FOR EMD3D3V_parallel_var
0017 % param.xend        = 5; %x-domain length in physical units
0018 % param.yend        = 5; %y-domain length in physical units
0019 % param.zend        = 5; %z-domain length in physical units
0020 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0021 
0022 X = linspace(0,30,param.n_grid_1);
0023 Y = linspace(0,15,param.n_grid_2);
0024 Z = linspace(0,15,param.n_grid_3);
0025 
0026 [x,y,z] = meshgrid(Y,X,Z); 
0027 
0028 uc1 = (sin(3.5*x)+sin(3.5*y)+sin(3.5*z));  uc2 = (sin(x)+sin(y)+sin(z)); ru  = (cos(x/24)+cos(y/24)+cos(z/24));
0029 u   =  uc1 + uc2 + ru; %Signal 'u'
0030 
0031 vc1 = (sin(x)+sin(y)+sin(z));  vc2 = (sin(3.5*x)+sin(3.5*y)+sin(3.5*z)); rv  = (cos(x/24)+cos(y/24)+cos(z/24));
0032 v   =  vc1 ; %Signal 'v'
0033 
0034 wc1 = (cos(x/24)+cos(y/24)+cos(z/24));  wc2 = (sin(x/5)+sin(y/5)+sin(z/5)); rw = (sin(x)+sin(y)+sin(z)) ;
0035 w   =  wc1 ;%+ wc2 + rw; %Signal 'w'
0036 
0037 Results = EMD3D3V(u,v,w,param);
0038 
0039 
0040