beamlet3dfig02.m

beamlet变化的工具箱

function Beamlet3Dfig02(OpenNewWindow,HaveTitle,LoadData)
LoadData=0;
if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end
n=16; dx=0; dy=0; dz=0;
v=[0 0 0;0 n 0;n n 0;n 0 0;0 0 n;0 n n;n n n;n 0 n];
v=v+[dx*ones(8,1),dy*ones(8,1),dz*ones(8,1)];
f=[1 2 3 4;5 6 7 8;1 5 6 2;2 6 7 3;3 7 8 4;4 8 5 1];
set(gca,'view',[-23,24],'Projection','orthographic');
axis off
axis image
hold on
[x,y]=meshgrid(0:n,0:n);
m2=mesh(x+dx,y+dy,n*ones(n+1)+dz,ones(n+1),'facecolor','n','Edgecolor','k');
m3=mesh(x+dx,zeros(n+1)+dy,y+dz,ones(n+1),'facecolor','n','Edgecolor','k');
m5=mesh(zeros(n+1)+dx,x+dy,y+dz,ones(n+1),'facecolor','n','Edgecolor','k');

r1=4; c1=4; r2=12; c2=13; 
dd=.5;

p=plot3([0 r2],[c1 c2],[r1 16],'r','lineWidth',2);

plot3([0 0],[c1-dd c1+dd],[r1 r1],'g','lineWidth',2);
plot3([0 0],[c1 c1],[r1-dd r1+dd],'g','lineWidth',2);
plot3([r2-dd r2+dd],[c2 c2],[16 16],'g','lineWidth',2)
plot3([r2 r2],[c2-dd c2+dd],[16 16],'g','lineWidth',2)

if HaveTitle
    title('A Beamlet')
end

%
% 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]
%