beamlet3dfig05.m

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

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

subplot (1,2,2)

r1=2; c1=2; r2=3; c2=2; 
dd=0.5;
n=16;
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];
n=n/2;
v=[v;0 0 n;0 n n;n n n;n 0 n;0 0 2*n;0 n 2*n;n n 2*n;n 0 2*n];
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];
f2=[6 7 8 16 15 14;1 4 8 16 12 9;1 2 6 14 10 9];
   f=[f;f+8];
   p=patch('vertices',v,'faces',f2,'facecolor',[0.7 0.7 0.7]);
   p=patch('vertices',v,'faces',f([7 10 11],:),'facecolor',[.9 .9 .9]);
   
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];
n=n/2;
v=[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];

set(gca,'view',[-30,30],'Projection','orthographic');
axis off
axis image
hold on
plot3([8 r2],[c1 8],[r1 c2]+8,'r')

plot3([8 8],[c1-dd c1+dd],[r1 r1]+8,'b');
plot3([8 8],[c1 c1],[r1-dd r1+dd]+8,'b');
plot3([r2-dd r2+dd],[8 8],[c2 c2]+8,'b');
plot3([r2 r2],[8 8],[c2-dd c2+dd]+8,'b');
if HaveTitle
    title('1st scale beamlet')
end

subplot(1,2,1)
r1=4; c1=4; r2=12; c2=13; 
dd=.5;
n=16;
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];
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];
p=patch('vertices',v,'faces',f,'facecolor','n');
set(gca,'view',[-30,30],'Projection','orthographic');
axis off
axis image
hold on

plot3([0 r2],[c1 c2],[r1 16],'r')

plot3([0 0],[c1-dd c1+dd],[r1 r1],'b');
plot3([0 0],[c1 c1],[r1-dd r1+dd],'b');
plot3([r2-dd r2+dd],[c2 c2],[16 16],'b')
plot3([r2 r2],[c2-dd c2+dd],[16 16],'b')
if HaveTitle
    title('0th scale 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]
%