plot.m

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

function plot(A)

% plot.m Plot method for radialrxrgeom objects. Called as
%
%          plot(A)
%
% INPUT PARAMETERS
%         A : an radialrxrgeom object

unpack(A.param)

th = [0:0.01:2*pi 2*pi]';

xplaten = Rp*cos(th);
yplaten = Rp*sin(th);

xinlet = Ro*cos(th);
yinlet = Ro*sin(th);

xRs = Rs*cos(th);
yRs = Rs*sin(th);

plot(yplaten,xplaten,yinlet,xinlet,'b',yRs,xRs,'b--', ...
    0,0,'k.','MarkerSize',15,'LineWidth',2)
grid

hold on

for rcontour = Rs-Rw:Rw/4:Rs+Rw
    
    xrcontour = rcontour*cos(th);
    yrcontour = rcontour*sin(th);
    
    plot(yrcontour,xrcontour,'r:')
    
end

xRw = Rw*sin(th);
yRw = Rw*cos(th);

nowafer = 5;

for iwafer = 1:nowafer
    
    xcenter = Rs*sin(2*pi*(iwafer-1)/nowafer);
    ycenter = Rs*cos(2*pi*(iwafer-1)/nowafer);
    plot(yRw+ycenter,xRw+xcenter,'g',0+ycenter,0+xcenter,'k.', ...
        'MarkerSize',15,'LineWidth',2)
    
end

% arrows and text

arrow([0 0],[Rp*cos(pi/6) -Rp*sin(pi/6)],7,'FaceColor','y','EdgeColor','k')
text(0.05,-0.05,'R','FontWeight','Bold','FontSize',14)
text(0.005,0.005,'R_o','FontWeight','Bold','FontSize',14)
text(0.01+Rp*cos(pi/6),-0.01-Rp*sin(pi/6),'R_p','FontWeight','Bold','FontSize',14)
text(-0.02+Rs,0.065,'R_s','FontWeight','Bold','FontSize',14)

arrow([0+Rs 0],[-0.0311+Rs 0.0251],7,'FaceColor','y','EdgeColor','k')
text(-0.02+Rs,0.005,'r','FontWeight','Bold','FontSize',14)
text(0.02+Rs,-0.02,'R_w','FontWeight','Bold','FontSize',14)

hold off

axis equal
% axis([-0.75*Rp+Rs 0.5*Rp+Rs -0.75*Rp 0.75*Rp])
xlabel('y')
ylabel('x')

box off
drawnow