drpdrawfig01.m

beamlet变化的工具箱

function dRPDrawFig01(OpenNewWindow,HaveTitle,LoadData)

%Plot FIO tiling, Ridgelet Tiling and Wavelet tiling

% if you don't want to process it, set LoadData=1

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end

LoadData=0;
%Plot FIO tiling, Ridgelet Tiling and Wavelet tiling


subplot(1,3,1)
hold on;
J=6;
color='b';
for j=0:J,
	rj = 2^(-j);
	ij = (J-j);
	nij = 2^(ij)*2;
	for l=0:(nij-1),
		thetal = 2*pi*l ./ nij;
		poutx = rj * cos(thetal); pouty = rj * sin(thetal);
		pinx  = poutx/2; piny = pouty/2;
		plot([poutx pinx],[pouty piny],color)
	end
    plotcircle(rj);
end
axis('square')
raxp = axis;
axis('off')
if HaveTitle
    title('Ridgelet Tiling')
end
hold off

subplot(1,3,2)
hold on


for j=0:J,
	rj = 2^(-j);
	
	ij = (J-j);
	nij = 2^(ij);%/2;%*2;
	for l=0:(nij-1),
		thetal = 2*pi*l ./ nij;
		poutx = rj * cos(thetal); pouty = rj * sin(thetal);
		pinx  = poutx/2^(j+2);
        piny = pouty/2^(j+2);
		plot([poutx pinx],[pouty piny],color)
	end
    plotcircle(rj);

end
axis('square')
raxp = axis;
axis('off')

if HaveTitle
    title('FIO Tiling')
end

hold off

subplot(1,3,3)
%figure
hold on
for j=0:J,
	rj = 2^(-j);
    plotcircle(rj);
end

plot([-1 1],[0 0])
plot([0 0],[1 -1])
axis('square')
raxp = axis;
axis('off')
if HaveTitle
    title('Wavelet Like Tiling ')
end
hold off

   
% Copyright (c) 2002 David Donoho and Ana Georgina Flesia

%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%