toon0212.m

老外写的小波变换的工具箱

% toon0212 -- 2d Wavelet Families
%
%   Two-dimensional wavelets may be formed by taking the tensor
%   product of two one-dimensional wavelet functions.
%
%   Here we plot the four tensor products:
%
%       Father(x)*Father(y)
%       Mother(x)*Father(y)
%       Father(x)*Mother(y)
%       Mother(x)*Mother(y)
%
%  for the S8 wavelet.
% 
	qmf = MakeONFilter('Symmlet',8);
	n=128;
	x = linspace(0,1.,n); 
	y = linspace(0,1.,n); 
%
	w = zeros(n,n);
	[kx,ky] = quad2ix(2,1,1,0);
	w(kx,ky)=1;
	z = IWT2_PO(w,0,qmf);
	subplot(221);
	mesh(x,y,z,z);
	title(' S8 Wavelet: FxF ');
%
	w = zeros(n,n);
	[kx,ky] = quad2ix(2,2,2,1);
	w(kx,ky) = 1;
	z = IWT2_PO(w,0,qmf);
	subplot(222);
	mesh(x,y,z,z);
	title(' S8 Wavelet: MxF ');
%
	w = zeros(n,n);
	[kx,ky] = quad2ix(2,2,1,2);
	w(kx,ky)=1;
	z = IWT2_PO(w,0,qmf);
	subplot(223);
	mesh(x,y,z,z);
	title(' S8 Wavelet: FxM ');
%
	w = zeros(n,n);
	[kx,ky] = quad2ix(2,2,1,3);
	w(kx,ky)=1;
	z = IWT2_PO(w,0,qmf);
	subplot(224);
	mesh(x,y,z,z);
	title(' S8 Wavelet: MxM ');
    
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:43 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]