beamlet3dfig16.m

来自「beamlet变化的工具箱」· M 代码 · 共 46 行

function Beamlet3Dfig16(OpenNewWindow,HaveTitle,LoadData)

%Auxillary functions: 3D beamlets

LoadData=0;
if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end
n=16;
ind1=[8,8,14,8;7,7,14,8;8,5,14,8;5,8,14,8;8,8,13,7;8,8,12,6;];
ind2=[8,8,14,8;9,9,14,8;8,11,14,8;11,8,14,8;8,8,15,9;8,8,16,10;];
for i=1:6
    subplot(2,3,i)
    B=zeros(n,n,n,n);
    B(ind1(i,1):ind2(i,1),ind1(i,2):ind2(i,2),ind1(i,3):ind2(i,3),ind1(i,4):ind2(i,4))=1;
    vol=adj_newbt3aa(B);
    p=patch(isosurface(vol,0),'facecolor','g');
    p2=patch(isocaps(vol,0),'facecolor','g');
    set(gca,'view',[70 12],'Xlim',[0 n],'Ylim',[0 n],'Zlim',[0 n]);
    set(gca,'Xtick',[],'Ytick',[],'Ztick',[])
    axis square
    light
    box on
    zoom(1.1)
    xlabel('x');ylabel('y');zlabel('z');
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]
%