plot.m

来自「MDPSAS工具箱是马里兰大学开发的」· M 代码 · 共 41 行

function plot(A)

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

unpack(A)

% Reconstruct solution

vx    = a*Bvx;
vy    = b*Bvy + vyBC;
press = d*Bp;

S = get(Bvx,'pd');

Vx = get(vx,'val');
Vy = get(vy,'val');
x  = get(S,'qp','x');
y  = get(S,'qp','y');
skip = 3;

figure(1), clf

plot(press)
xlabel('x','FontSize',14)
ylabel('y','FontSize',14)
zlabel('P','FontSize',14)

figure(2), clf

quiver( y(1:skip:end),x(1:skip:end),Vy(1:skip:end,1:skip:end),...
                                    Vx(1:skip:end,1:skip:end),2,'b');
axis([0 1 0 1])
xlabel('y','FontSize',14)
ylabel('x','FontSize',14)
title('Moving Wall -->','FontSize',14)