cwtex.m

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

%Script file that computes the  CWT of a rectangular pulse using 
%Haar wavelet. Demonstration of Example 1.1 in the book.
%
%Author: Raghuveer M. Rao 
%Copyright (c) 1998 by Addison Wesley Longman, Inc.
%

for a=-5:0.1:5,
  for b=-5:0.1:5,
	if a~=0
	CWT(round(a*10+51),round(b*10+51))=sign(a)*(2*tri(b+a/2)-tri(b)-tri(b+a))/sqrt(abs(a));
	end
  end
end
CWT(51,:)=zeros(1,101);
a=-5:0.1:5;
b=-5:0.1:5;
figure;contour(b,a,CWT.^2,20)