www.pudn.com > MATLAB.rar > my_cwt_3.m, change:2014-08-30,size:1308b


% ------------------------------------------------------------------------- 
% 
%   Archivo:        my_cwt.m 
% 
%   Descripción:    version 3 del algoritmo cwt propio. Se modificó el 
%   algoritmo sin usar diff ni int 
% 
%   Autor:          ROUX, Federico G. (rouxfederico@gmail.com) 
% 
% 
% ------------------------------------------------------------------------- 
 
 
function [y] = my_cwt_3(x, escalas, wavelet) 
 
N = length(x);                                                              % Longitud de la señal 
 
escalas_N = size(escalas, 2);                                               % Cantidad de escalas a evaluar 
 
y = zeros(escalas_N, N); 
 
li = -8; 
ls = -li; 
 
for i = 1:escalas_N 
     
    a = escalas(i);                                                         % Guardo en 'a' la escala a procesar 
 
    L = 16*a + 1;                                                           % Calculo el largo de la wavelet  
    wsub = wavelet(li, ls, L);                                              % Recalculo la wavelet 
             
    yt = myConv(x, wsub);                                           % Hago la convolución entre la señal y la wavelet     
    y(i, :) = -sqrt(1/a) * wkeep1(yt, N);                             % Me quedo con la parte central 
        
end 
     
end