Home > FA-MVEMD > 3D > driver_3d_example.m

# driver_3d_example

## PURPOSE

cleanup

## SYNOPSIS

This is a script file.

## DESCRIPTION

    cleanup

## CROSS-REFERENCE INFORMATION

This function calls:

- EMD3D3V Purpose:

This function is called by:

## SOURCE CODE

    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