dmakessfigure6.m

beamlet变化的工具箱

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

% MakeSSFig4 -- Draw Figure showing the definitions of lines.
%figure
gcf; subplot(1,1,1)

n = 8;
w = zeros(2,n+1,n);
xi = zeros(2,2,n+1,n);
Total = n^2+1;
for s=1:2,
	first = 3-s;
	secnd = s;
	if s==1, sg = 1; else sg= -1; end
	for l=1:(n+1),
		for m=1:n,
			lt = l-1 -n/2;
			mt = m-1 -n/2;
			if lt==0,
				w(s,l,m) = 1/(2*n);
			elseif abs(lt) < n/2,
				w(s,l,m) = 2*abs(lt)/n;
			elseif abs(mt) < n/2,
				w(s,l,m) = 1/2;
			else
				w(s,l,m) = 1/4;
			end
			xi(first,s,l,m) = 2*pi*lt/n;
			xi(secnd,s,l,m) = sg*2*pi*mt/n * lt*2/n;
		end
	end
end
%
% Display lines indicating plot extent
%
%plot([-pi, pi,pi,-pi,-pi],[-pi, -pi,pi,pi,-pi],'-k'); hold;
plot([-pi pi],[-pi -pi],'r');
axis([-4 4 -4 4]); 
axis('equal');
hold on;
plot([-pi pi],[pi pi],'r');
plot([-pi -pi],[-pi pi],'b');
plot([pi pi],[-pi pi],'b');
for s=1:2, 
	if s==1,
	 	pch = '-.r';
		pc  = 'xr';
	else
	 	pch = ':b';
		pc  = 'xb';
	end
	for m=1:n, 
		plot(xi(1,s,(n+1),m),xi(2,s,(n+1),m),pc); 
        plot(squeeze(xi(1,s,[1 (n+1)],m)), squeeze(xi(2,s,[1 (n+1)],m)),pch); 
	end ; 
end
hold off; 
axis off;
%drawnow; figure(gcf);
%print -depsc RawSSFig4.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]
%