%% ====================================================
%% FUNCIÓN ANALITICA Y ENVOLVENTE
%%
%% !OBS: El directorio debe contener el código
%%       hilbert_ricker.m que genera un ondícula de
%%       Ricker y calcula su transformada de Hilbert.
%%       El directorio debe contener también la función 
%%       ricker.m
%% ====================================================

clear all; close;

% ---------------------------------------------------
% 1. Corro el código hilbert.m que genera la ondícula
%    de Ricker (dato) y su transformada de Hilbert.
% ---------------------------------------------------

hilbert_ricker ; % rw: ondícula de Ricker
                 % hw; transformada de Hilbert de rw

% --------------------
% 2. Función analítica
% --------------------
 
gw = rw - 1i*hw ;

% --------------------------------------------
% 3. Envolvente (módulo de la señal analítica)
% --------------------------------------------

env = abs(gw) ;

% -----------
% 4. GRÁFICOS
% -----------

t=(-M:M)*dt; % vector de tiempo

figure(1)
clf;
hold on ;
plot(t,rw,';Ondicula de Ricker;b-') ;
plot(t,hw,';Transformada de Hilbert;g-') ;
plot(t,env,';Envolvente;r-','linewidth',2) ;
plot(t,-env,'r-','linewidth',2) ;
axis([-M*dt,M*dt,-1.1,1.1])
hold off

% Guardo el gráfico en un postscript
print('fun_analitica_time.eps','-depsc')
