dmakessfigure11_4.m

beamlet变化的工具箱

function MakeSSFigure11_4(varargin)
% MakeSSFigure11_4.m -- Figure 11_4 of Slant Stack Paper
% Usage: MAKESSFIGURE11_4 will start a new figure window and plot Figure 11_4 of Slant Stack Paper.
%        MAKESSFIGURE11_4(n) for any argument type n will plot the figure on the current figure 
%        window and plot Figure 11_4 of Slant Stack Paper.
if nargin == 0
	 figure;
else
     gcf;
end

% MakeOnion0 -- Show Sector Conversion Problem
%figure
	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
%
% 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==1, 
		for l=1:n, 
			plot(squeeze(xi(1,s,4:(m-2),l)),squeeze(xi(2,s,4:(m-2),l)),'xb'); 
%			plot(xi(1,s,[1 (m+1)],l),xi(2,s,[1 (m+1)],l),':b'); 
		end ;
	elseif s==2, 
		for l=1:n, 
			plot(squeeze(xi(1,s,3:(m-1),l)),squeeze(xi(2,s,3:(m-1),l)),'xb'); 
%			plot(xi(1,s,[1 (m+1)],l),xi(2,s,[1 (m+1)],l),':b'); 
		end ;

% 	else
% 		for k=1:(m+1), 
% % 			plot([xi(1,s,(n+1),m), (-xi(1,s,(n+1),m))],...
% % 				 [xi(2,s,(n+1),m), (-xi(2,s,(n+1),m))],'x');
% 			plot([xi(1,s,k,1), (xi(1,s,k,1))],...
% 				 [pi, (-pi)],'-r');
% 		end
	end 
end
for k1=-8:8,
  for k2 =-8:8,
  	if abs(k2) == abs(k1),
		plot(pi*k1/8,pi*k2/8,'ro');
	elseif abs(k1) == 8 | abs(k2)==8,
		plot(pi*k1/8,pi*k2/8,'ro');
	elseif abs(k1)==7,
		plot(pi*k1/8,pi*k2/8,'ro');
	elseif abs(k2)==7 | abs(k2)==6,
		plot(pi*k1/8,pi*k2/8,'ro');		
	end	
  end
end
w1 = 6.5*pi/8; w2 = 5.5*pi/8; ed = 3.3;
plot([-w1 -w2 -w2 -w1 -w1],[-ed -ed ed ed -ed],'-k')
plot([w1 w2 w2 w1 w1],[-ed -ed ed ed -ed],'-k')
axis([-4 4 -4 4])
hold off;
axis off;
%drawnow; figure(gcf);
%print -depsc SSOnion3.eps
%
% 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]
%