beamlet3dfig09.m

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

function Beamlet3Dfig09(OpenNewWindow,HaveTitle,LoadData)
LoadData=0;
if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end

ind=[16 1;16 5;16 9;16 13;16 17;16 21;16 25;16 29;16 32];
for k=1:9
     i=ind(k,1);
     s=ind(k,2);
     subplot(3,3,k);
     if s<17
         vertices=[0,0,-8;0,16,-8;16,16,-8;16,0,-8;0,0,24;0,16,24;16,16,24;16,0,24;0,0,0;0,16,0;...
              16,16,0;16,0,0;0,0,16;0,16,16;16,16,16;16,0,16];
     else
         vertices=[0,-8,0;0,-8,16;16,-8,16;16,-8,0;0,24,0;0,24,16;16,24,16;16,24,0;0,0,0;0,16,0;...
               16,16,0;16,0,0;0,0,16;0,16,16;16,16,16;16,0,16];
     end

     faces=[1,2,3,4;5,6,7,8;1,2,6,5;3,4,8,7;9,10,11,12;13,14,15,16;9,10,14,13;11,12,16,15];
     p=patch('faces',faces,'vertices',vertices,'facecolor','n','EdgeColor','b');
     axis image
     hold on
     if s<17
         [x,y]=meshgrid(0:16,0:16);
         z=(i/2)+(-1.125+s/8)*(8-y);
         m=mesh(x,y,z,~((z>=16)+(z<0)),'FaceColor','flat','EdgeColor','r');
     else
         [x,z]=meshgrid(0:16,0:16);
         y=(i/2)+(3.125-s/8)*(8-z);
         m=mesh(x,y,z,~((y>=16)+(y<=0)),'FaceColor','flat','EdgeColor','r');

   
     end
     set(gca,'view',[-38,30]);
     colormap summer
     if HaveTitle
         ttlstr=strcat('b(:,:,',num2str(ind(k,1)),',',num2str(ind(k,2)),')');
         title(ttlstr)
     end
     set(gca,'XTick',[],'YTick',[],'ZTick',[]);
     h=xlabel('x');
     set(h,'position',get(h,'position')+[5 8 0]);
     h=ylabel('y');
     set(h,'position',get(h,'position')+[7 9 0]);
     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]
%