1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
|
figure(2)
clf;
set(gcf, 'color', 'white')
set(0,'defaultaxesfontname','Arial');
n=10;
tx = 0.54;
ty = 0.084;
title(['a = ',num2str(Rai),'\delta = ',num2str(delta1*180/pi),'° and \Delta x = ',num2str(dx)])
subplot(3,1,1)
k=1;time=round(T(k+1)*100);
[AX,H1,H2]=plotyy(Axx,R(ns:nx,n,k),Axx,U(n,ns:nx,k));
set(get(AX(1),'Ylabel'),'String','Flow thickness (m)') ;
set(get(AX(2),'Ylabel'),'String','Velocity (m.s^{-1})') ;
set(AX(1), 'ylim', [0 .1],'YTick',0:.05:.1,'xlim', [0 xmax]);
set(AX(2), 'ylim', [0 2],'YTick',0:1:2,'xlim', [0 xmax]);
set(get(AX(2),'Xlabel'),'String','Distance (m)') ;
axis equal
text(tx,ty,['Time: ',num2str(time),' ms.']);
subplot(3,1,2)
k=4;time=round(T(k+1)*100);
[AX,H1,H2]=plotyy(Axx,R(ns:nx,n,k),Axx,U(n,ns:nx,k));
set(get(AX(1),'Ylabel'),'String','Flow thickness (m)') ;
set(get(AX(2),'Ylabel'),'String','Velocity (m.s^{-1})') ;
set(AX(1), 'ylim', [0 .1],'YTick',0:.05:.1,'xlim', [0 xmax]);
set(AX(2), 'ylim', [0 2],'YTick',[0:1:2], 'xlim', [0 xmax]);
set(get(AX(2),'Xlabel'),'String','Distance (m)') ;
axis equal
text(tx,ty,['Time: ',num2str(time),' ms.']);
subplot(3,1,3)
k=10;time=round(T(k+1)*100);
[AX,H1,H2]=plotyy(Axx,R(ns:nx,n,k),Axx,U(n,ns:nx,k));
set(get(AX(1),'Ylabel'),'String','Flow thickness (m)') ;
set(get(AX(2),'Ylabel'),'String','Velocity (m.s^{-1})') ;
set(get(AX(2),'Xlabel'),'String','Distance (m)') ;
set(AX(2), 'ylim', [0 2],'YTick',[0:1:2],'xlim', [0 xmax]);
set(AX(1), 'ylim', [0 .1],'YTick',[0:.05:.1],'xlim', [0 xmax]);
axis equal
text(tx,ty,['Time: ',num2str(time),' ms.']); |
Partager