drpdrawfig07.m

beamlet变化的工具箱

function dRPDrawFig07(OpenNewWindow,HaveTitle,LoadData)

%DrawFig 7:  Make Pseudo-Polar Grid (Fourier Space)


if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


LoadData=0;



%DrawFig 7:  Make Pseudo-Polar Grid (Fourier Space)
	
n = 8;	
m = 2*n;
w = zeros(2,m+1,n);
xi = zeros(2,2,m+1,n);
%
for s=1:2,
		first = 3-s;
		secnd = s;
		if s==1, sg = 1; else sg= -1; end
		for k=1:(m+1),
			for l=1:n,
				kt = k-1 -n;
				lt = l-1 -n/2;
				if kt==0,
					w(s,k,l) = 1/(2*n);
				elseif abs(kt) < n,
					w(s,k,l) = abs(kt)/n;
				elseif abs(lt) < n/2,
					w(s,k,l) = 1/2;
				else
					w(s,k,l) = 1/4;
				end
				xi(first,s,k,l) = pi*kt/n;
				xi(secnd,s,k,l) = sg*2*pi*lt/n * kt/n;
			end
		end
	end
%
%figure
% Display lines indicating plot extent
%
plot([-pi pi],[-pi -pi],'w');
axis([-4 4 -4 4]); 
axis('equal');
hold on;
plot([-pi pi],[pi pi],'w');
plot([-pi -pi],[-pi pi],'w');
plot([pi pi],[-pi pi],'w');

for s=1:2, 
	if s==2, 
		for l=1:n,   % radial lines
			plot(squeeze(xi(1,s,[1 (m+1)],l)),squeeze(xi(2,s,[1 (m+1)],l)),'-r'); 
		end ;
%	else			% horizontal lines
		for k=1:(m+1), 
			plot([xi(1,s,k,1), xi(1,s,k,1)],...
				 [xi(2,s,k,1), (-xi(2,s,k,1))],'-r'); 
		end
	end 
end
for s=1:2, 
	if s==1, 
		for l=1:n, 
			plot(squeeze(xi(1,s,[1 (m+1)],l)),squeeze(xi(2,s,[1 (m+1)],l)),'-r'); 
		end 
%	else
		for k=1:(m+1), 
			plot([xi(1,s,k,1), -xi(1,s,k,1)],...
				 [xi(2,s,k,1), (xi(2,s,k,1))],'-r'); 
		end
	end 
end

if HaveTitle
    title('Pseudo-Polar Grid')
end
hold off;
axis 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]
%