plotold.m

MDPSAS工具箱是马里兰大学开发的

function plot(A,R)

% plot.m Displays a relquad object graphically in relation
%        to a quadgrid object. Called as
%
%        plot(A,R)
%
%        Input Parameters:
%
%        A : the base quadgrid object
%        R : the relquad object
%

if isa(A,'relquad')
   plot(A.quadgrid)  
else
   plot(A);
   end

hold on

geom = get(R,'geom');
qp   = get(R,'qp');
name = get(R,'name');

baseGeom = get(A,'geom');

nodir = size(qp,2);

switch nodir

case 1

   qp{1} = qp{1}+R.offset(1);

   plot(ones(size(qp{1})),qp{1},'ro','MarkerSize',5)
   ylabel(name{1})

case 2

   if strcmp(baseGeom{2},'peri') % translate (base) polar coordinates to Cartesian

      x0 = R.offset(1)*cos(R.offset(2));
      y0 = R.offset(1)*sin(R.offset(2));
      R.offset(1) = x0; R.offset(2) = y0;
      end

   if strcmp(geom{1},'slab')

      qp{1} = qp{1}+R.offset(1);
      qp{2} = qp{2}+R.offset(2);

      plot(qp{2}(1)*ones(size(qp{1})),qp{1},'ro','MarkerSize',5)
      hold on
      for i = 2:size(qp{2},1)
         plot(qp{2}(i)*ones(size(qp{1})),qp{1},'ro','MarkerSize',5)
         end
      hold off
      xlabel(name{2})
      ylabel(name{1})
      box on
      end

   if strcmp(geom{1},'cyln') & strcmp(geom{2},'peri')
      for i = 1:size(qp{1},1)
         for j = 1:size(qp{2},1)
            x = qp{1}(i)*cos(qp{2}(j))+R.offset(1);
            y = qp{1}(i)*sin(qp{2}(j))+R.offset(2);
            plot(y,x,'ro','MarkerSize',5)
            if i == 1 & j == 1, hold on, end
            end
         end
      hold off
      xlabel('2nd dimension')
      ylabel('1st dimension')
      box on
      end

case 3

   qp{1} = qp{1}+R.offset(1);
   qp{2} = qp{2}+R.offset(2);
   qp{3} = qp{3}+R.offset(3);

   axis([0 1 0 1 0 1])
   hold on

   for i = 1:size(qp{1},1)
      plot3( [qp{1}(i),   qp{1}(i)], ...
             [qp{2}(1),   qp{2}(1)], ...
             [qp{3}(1),   qp{3}(end)], 'b-')
      plot3( [qp{1}(i),   qp{1}(i)], ...
             [qp{2}(end), qp{2}(end)], ...
             [qp{3}(1),   qp{3}(end)], 'b-')
      plot3( [qp{1}(i),   qp{1}(i)], ...
             [qp{2}(1),   qp{2}(end)], ...
             [qp{3}(1),   qp{3}(1)], 'b-')
      plot3( [qp{1}(i),   qp{1}(i)], ...
             [qp{2}(1),   qp{2}(end)], ...
             [qp{3}(end), qp{3}(end)], 'b-')
      end

   for i = 1:size(qp{2},1)
      plot3( [qp{1}(1),   qp{1}(1)], ...
             [qp{2}(i),   qp{2}(i)], ...
             [qp{3}(1),   qp{3}(end)], 'b-')
      plot3( [qp{1}(end), qp{1}(end)], ...
             [qp{2}(i),   qp{2}(i)], ...
             [qp{3}(1),   qp{3}(end)], 'b-')
      plot3( [qp{1}(1),   qp{1}(end)], ...
             [qp{2}(i),   qp{2}(i)], ...
             [qp{3}(1),   qp{3}(1)], 'b-')
      plot3( [qp{1}(1),   qp{1}(end)], ...
             [qp{2}(i),   qp{2}(i)], ...
             [qp{3}(end), qp{3}(end)], 'b-')
      end

   for i = 1:size(qp{3},1)
      plot3( [qp{1}(1),   qp{1}(1)], ...
             [qp{2}(1),   qp{2}(end)], ...
             [qp{3}(i),   qp{3}(i)], 'b-')
      plot3( [qp{1}(end), qp{1}(end)], ...
             [qp{2}(1),   qp{2}(end)], ...
             [qp{3}(i),   qp{3}(i)], 'b-')
      plot3( [qp{1}(1),   qp{1}(end)], ...
             [qp{2}(1),   qp{2}(1)], ...
             [qp{3}(i),   qp{3}(i)], 'b-')
      plot3( [qp{1}(1),   qp{1}(end)], ...
             [qp{2}(end), qp{2}(end)], ...
             [qp{3}(i),   qp{3}(i)], 'b-')
      end

   hold off
   xlabel(name{1})
   ylabel(name{2})
   zlabel(name{3})
   box on

otherwise
   disp('Too many dimensions')
   end

hold off