beamlet3dfig07.m

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

function Beamlet3Dfig07(OpenNewWindow,HaveTitle,LoadData)

%Auxillary function: dividing_planes

LoadData=0;
if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end

subplot(2,2,1);
dividing_planes
dd=.03;
p1=[0 .75 .75];
p2=[.4 1 .9];
p=plot3([p1(1) p2(1)],[p1(2) p2(2)],[p1(3) p2(3)],'r','linewidth',2);

plot3([p1(1) p1(1)],[p1(2)-dd p1(2)+dd],[p1(3) p1(3)],'b');
plot3([p1(1) p1(1)],[p1(2) p1(2)],[p1(3)-dd p1(3)+dd],'b');
plot3([p2(1)-dd p2(1)+dd],[p2(2) p2(2)],[p2(3) p2(3)],'b');
plot3([p2(1) p2(1)],[p2(2) p2(2)],[p2(3)-dd p2(3)+dd],'b');
zoom(0.8)

if HaveTitle
    t=title('a beamlet belongs to a single child');
    %set(t,'fontunits','normalized','fontsize',8);
end

subplot(2,2,2)
dividing_planes
p1=[.4 1 .9];
p2=[.2 0 .9];
p=plot3([p1(1) p2(1)],[p1(2) p2(2)],[p1(3) p2(3)],'r','linewidth',2);
plot3([p1(1)-dd p1(1)+dd],[p1(2) p1(2)],[p1(3) p1(3)],'b');
plot3([p1(1) p1(1)],[p1(2) p1(2)],[p1(3)-dd p1(3)+dd],'b');
plot3([p2(1)-dd p2(1)+dd],[p2(2) p2(2)],[p2(3) p2(3)],'b');
plot3([p2(1) p2(1)],[p2(2) p2(2)],[p2(3)-dd p2(3)+dd],'b');
zoom(0.8)

if HaveTitle
    title('decomposition into 2 continuum beamlets')
end

subplot(2,2,3)
dividing_planes
p1=[0 .75 .9];
p2=[1 0.1 .7];
p=plot3([p1(1) p2(1)],[p1(2) p2(2)],[p1(3) p2(3)],'r','linewidth',2);

plot3([p1(1) p1(1)],[p1(2)-dd p1(2)+dd],[p1(3) p1(3)],'b');
plot3([p1(1) p1(1)],[p1(2) p1(2)],[p1(3)-dd p1(3)+dd],'b');
plot3([p2(1) p2(1)],[p2(2)-dd p2(2)+dd],[p2(3) p2(3)],'b');
plot3([p2(1) p2(1)],[p2(2) p2(2)],[p2(3)-dd p2(3)+dd],'b');
zoom(0.8)

if HaveTitle
    title('decomposition into 3 continuum beamlets')
end

subplot(2,2,4)
dividing_planes
p1=[0 .75 .9];
p2=[1 0.1 .35];
p=plot3([p1(1) p2(1)],[p1(2) p2(2)],[p1(3) p2(3)],'r','linewidth',2);

plot3([p1(1) p1(1)],[p1(2)-dd p1(2)+dd],[p1(3) p1(3)],'b');
plot3([p1(1) p1(1)],[p1(2) p1(2)],[p1(3)-dd p1(3)+dd],'b');
plot3([p2(1) p2(1)],[p2(2)-dd p2(2)+dd],[p2(3) p2(3)],'b');
plot3([p2(1) p2(1)],[p2(2) p2(2)],[p2(3)-dd p2(3)+dd],'b');
zoom(0.8)

if HaveTitle
    title('decomposition into 4 continuum beamlets')
end
set(GCF, 'Renderer','Zbuffer');

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