This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: on plot(Mone, Tratio, 'k' , 'LineWidth' ,2) plot(Mone, rratio, '.k' , 'LineWidth' ,2) grid on legend( 'p_2/p_1' , 'T_2/T_1' , '\rho_2/\rho_1' ) set(gca, 'YLim' ,[0,20]) xlabel( 'Upstream Mach number, M_1' ) ylabel( 'Property Ratio' ) hold off 1 1.5 2 2.5 3 3.5 4 4.5 5 2 4 6 8 10 12 14 16 18 20 Upstream Mach number, M 1 Property Ratio p 2 /p 1 T 2 /T 1 2 / 1 Mone = [1:0.01:5]; pratio = NS_p2p1(Mone); Tratio = NS_T2T1(Mone); rratio = NS_r2r1(Mone); % clf semilogy(Mone, pratio, 'k' , 'LineWidth' ,2) hold on semilogy(Mone, Tratio, 'k' , 'LineWidth' ,2) semilogy(Mone, rratio, '.k' , 'LineWidth' ,2) grid on legend( 'p_2/p_1' , 'T_2/T_1' , '\rho_2/\rho_1' ) set(gca, 'YLim' ,[0,20]) xlabel( 'Upstream Mach number, M_1' ) ylabel( 'Property Ratio' ) hold off set(gca, 'YTick' ,[1:10 20]) 1 1.5 2 2.5 3 3.5 4 4.5 5 1 2 3 4 5 6 7 8 9 10 20 Upstream Mach number, M 1 Property Ratio p 2 /p 1 T 2 /T 1 2 / 1...
View Full
Document
 Fall '09
 COLLICOTT

Click to edit the document details