|
发表于 2021-12-13 22:00:34
|
显示全部楼层
%% Figure(1)
figure(1)
if h == 1 || h == 2
if h == 1
subplot(2,2,h)
else
subplot(2,2,h+1)
end
hold on
p1 = plot(time,cbsim./csat,'-','LineWidth',1.0);
p2 = plot(X,Y./csat,'LineWidth',1.0,'LineStyle',':');
p3 = plot(t,cb./csat,'LineWidth',1.0);
legend([p1 p2 p3],{'c_{b} (PSS)','c_{b} (LT)','c_{b} (OM)'},'Location','southeast')
axis([0 t(end) 0 cbsim(end)./csat*1.05])
xlabel('$$ Time (s) $$','interpreter','latex')
ylabel('$$ \frac{c_{b}}{c_{sat}}\ \left( \frac{g\ m^{-3}}{g\ m^{-3}} \right) $$','interpreter','latex')
grid on
grid minor
box on
if h == 1
title({'k = 10^{-7.0} s^{-1}';['Ha = ',num2str(Ha)]})
elseif h == 2
title({'k = 10^{-6.0} s^{-1}';['Ha = ',num2str(Ha)]})
end
if h == 1
subplot(2,2,h+1)
else
subplot(2,2,h+2)
end
p1 = plot(time,cbsim./csat,'-','LineWidth',1.0);
p2 = plot(X,Y./csat,'LineWidth',1.0,'LineStyle',':');
p3 = plot(t,cb./csat,'LineWidth',1.0);
legend([p1 p2 p3],{'c_{b} (PSS)','c_{b} (LT)','c_{b} (OM)'},'Location','southeast')
axis([0 t(end) 0 1.05])
xlabel('$$ Time (s) $$','interpreter','latex')
ylabel('$$ \frac{c_{b}}{c_{sat}}\ \left( \frac{g\ m^{-3}}{g\ m^{-3}} \right) $$','interpreter','latex')
if h == 1
title({'k = 10^{-7.0} s^{-1}';['Ha = ',num2str(Ha)]})
elseif h == 2
title({'k = 10^{-6.0} s^{-1}';['Ha = ',num2str(Ha)]})
end
hfig1 = figure(1);
set(hfig1, 'position', get(0,'ScreenSize'));
print('0 k = 10^(-7 -6) dimensionless','-djpeg','-r1200')
print('0 k = 10^(-7 -6) dimensionless','-dpdf','-fillpage') %% Figure(1)
figure(1)
if h == 1 || h == 2
if h == 1
subplot(2,2,h)
else
subplot(2,2,h+1)
end
hold on
p1 = plot(time,cbsim./csat,'-','LineWidth',1.0);
p2 = plot(X,Y./csat,'LineWidth',1.0,'LineStyle',':');
p3 = plot(t,cb./csat,'LineWidth',1.0);
legend([p1 p2 p3],{'c_{b} (PSS)','c_{b} (LT)','c_{b} (OM)'},'Location','southeast')
axis([0 t(end) 0 cbsim(end)./csat*1.05])
xlabel('$$ Time (s) $$','interpreter','latex')
ylabel('$$ \frac{c_{b}}{c_{sat}}\ \left( \frac{g\ m^{-3}}{g\ m^{-3}} \right) $$','interpreter','latex')
grid on
grid minor
box on
if h == 1
title({'k = 10^{-7.0} s^{-1}';['Ha = ',num2str(Ha)]})
elseif h == 2
title({'k = 10^{-6.0} s^{-1}';['Ha = ',num2str(Ha)]})
end
if h == 1
subplot(2,2,h+1)
else
subplot(2,2,h+2)
end
p1 = plot(time,cbsim./csat,'-','LineWidth',1.0);
p2 = plot(X,Y./csat,'LineWidth',1.0,'LineStyle',':');
p3 = plot(t,cb./csat,'LineWidth',1.0);
legend([p1 p2 p3],{'c_{b} (PSS)','c_{b} (LT)','c_{b} (OM)'},'Location','southeast')
axis([0 t(end) 0 1.05])
xlabel('$$ Time (s) $$','interpreter','latex')
ylabel('$$ \frac{c_{b}}{c_{sat}}\ \left( \frac{g\ m^{-3}}{g\ m^{-3}} \right) $$','interpreter','latex')
if h == 1
title({'k = 10^{-7.0} s^{-1}';['Ha = ',num2str(Ha)]})
elseif h == 2
title({'k = 10^{-6.0} s^{-1}';['Ha = ',num2str(Ha)]})
end
hfig1 = figure(1);
set(hfig1, 'position', get(0,'ScreenSize'));
print('0 k = 10^(-7 -6) dimensionless','-djpeg','-r1200')
print('0 k = 10^(-7 -6) dimensionless','-dpdf','-fillpage') |
|