plot.m

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

function plot(A,profileonly)

% plot.m Plot method for aixtron objects. Called as
%
%          plot(A,profileonly)
%
% INPUT PARAMETERS
%            A : an aixtron object
%  profileonly : a flag which if set to 1 plots the
%                composition profiles only (optional)

if nargin == 1, profileonly = 0; end

unpack(A)

if ~profileonly

    subplot(2,2,1)
    plot(velocity(A),2)
    xlabel('R (m)')
    ylabel('v (m/s)')
    axis tight
    box off
    grid

    subplot(2,2,2)
    plot(T,2)
    xlabel('R (m)')
    ylabel('T (K)')
    axis tight
    ax = axis;
    axis([ax(1:3) 1.1*ax(4)])
    box off
    grid

    subplot(2,1,2)
else
    subplot(2,1,1)
end

plot(xA,xB,xC,2)
axis tight
ax = axis;
hold on
plot( [Rs-Rw Rs-Rw],[-1 1],'k', ...
    [Rs+Rw Rs+Rw],[-1 1],'k' )
plot([Rs Rs],[-1 1],'k-.')
legend('x_A','x_B','x_C')
xlabel('R (m)')
ylabel('x_i')
grid

for i = ax(1):0.005:ax(2)
    for j = ax(3):0.005:ax(4)
        if i < Rs-Rw | i > Rs+Rw
            plot(i,j,'b.','MarkerSize',2)
        end
    end
end

hold off
drawnow