ourphase.m

斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善

function phi = ourphase(g,epsilon);
% phi = ourphase(g);
% calculates the complex phase of vector g ;
% coordinates whose moduli are smaller than a defined 
% threshold epsilon are associated to phi = -pi/2;



phase = rem(imag(log(g)),2*pi.*ones(size(g)));
phi = phase.*(abs(g) > epsilon) - pi * (abs(g) <= epsilon);

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]
    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%