haar.m

这是伯克里wavelet transforms一书中的例子的代码

function x=haar(a,ts)
%ha = haar(a,ts) computes the Haar wavelet for a given dilation factor, a, and 
%sampled at a rate of ts. 
%This routine is used by the routine, ctwt.
%
%Author: Ajit S. Bopardikar
%Copyright (c) 1998 by Addison Wesley Longman, Inc.
%

   if(a~=0)
    x = abs(a)/2 - (0:ts:abs(a));
    x = sign(a)*sign(x)/sqrt(abs(a));
  else
   fprintf('Dilation factor of 0 encountered. Cant compute wavelet');
  end %endif