clear
c0=(1+sqrt(3))/4*sqrt(2);
c1=(3+sqrt(3))/4*sqrt(2);
c2=(3-sqrt(3))/4*sqrt(2);
c3=(1-sqrt(3))/4*sqrt(2);
N=64;
I=eye(N);
for k=1:N,
f=[f;daubidwt(c0,c1,c2,c3,I(:,k),2)'];
end

figure(1) 
imshow(mat2gray(f),256)
title('Daubechies wavelet basis functions for N=64')

%f=daubidwt(c0,c1,c2,c3,[0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;zeros(48,1)],2);
%plot (f,':');
%hold
%f=daubidwt(c0,c1,c2,c3,[0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;zeros(48,1)],2);
%plot (f);

%f=daubidwt(c0,c1,c2,c3,[0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;zeros(48,1)],2);
%plot(f,'-.');
%title('Scaled and Translated Daubechies wavelets for 4th, 5th, and 6th basis')
%hold
