diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m | 45 |
1 files changed, 22 insertions, 23 deletions
diff --git a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m index 8bdf885..5575956 100644 --- a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m +++ b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m @@ -4,24 +4,24 @@ Nlevels_with_MOR z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS figure(1) -subplot(2,2,1); imagesc(t_1, z_1, IdL_out_1); colorbar -ylabel('z (cm)') -xlabel('t (uS)') +subplot(2,2,1); imagesc(z_1, flip(t_1), IdL_out_1.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') zlabel('I_{dL}') title('I_{dL}') -subplot(2,2,2); imagesc(t_1, z_1, IpL_out_1); colorbar -ylabel('z (cm)') -xlabel('t (uS)') +subplot(2,2,2); imagesc(z_1, flip(t_1), IpL_out_1.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') zlabel('I_{pL}') title('I_{pL}') -subplot(2,2,3); imagesc(t_1, z_1, IdR_out_1); colorbar -ylabel('z (cm)') -xlabel('t (uS)') +subplot(2,2,3); imagesc(z_1, flip(t_1), IdR_out_1.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') zlabel('I_{dR}') title('I_{dR}') -subplot(2,2,4); imagesc(t_1, z_1, IpR_out_1); colorbar -ylabel('z (cm)') -xlabel('t (uS)') +subplot(2,2,4); imagesc(z_1, flip(t_1), IpR_out_1.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') zlabel('I_{pR}') title('I_{pR}') @@ -34,8 +34,8 @@ title('I_{pR}') figure(5) subplot(2,2,1); plot( ... - t_1,IdL_out_1(1,:)', ... - t_1,IdL_out_1(end,:)','LineWidth', 4 ... + t_1, flip(IdL_out_1(1,:)'), ... + t_1, flip(IdL_out_1(end,:)'),'LineWidth', 4 ... ) xlabel('t (uS)') ylabel('I_{dL} (1/s)^2') @@ -45,8 +45,8 @@ legend('before', 'after') %% subplot(2,2,2); plot( ... - t_1,IpL_out_1(1,:)', ... - t_1,IpL_out_1(end,:)', 'linewidth', 4 ... + t_1, flip(IpL_out_1(1,:)'), ... + t_1, flip(IpL_out_1(end,:)'), 'linewidth', 4 ... ) xlabel('t (uS)') ylabel('I_{pL} (1/s)^2') @@ -56,8 +56,8 @@ legend('before', 'after') %% subplot(2,2,3); plot( ... - t_1,IdR_out_1(1,:)', ... - t_1,IdR_out_1(end,:)', 'linewidth', 4 ... + t_1, flip(IdR_out_1(1,:)'), ... + t_1, flip(IdR_out_1(end,:)'), 'linewidth', 4 ... ) xlabel('t (uS)') ylabel('I_{dR} (1/s)^2') @@ -91,8 +91,8 @@ display( strcat('Second field delay time = ', num2str(delay_time), ' uS/n')); subplot(2,2,4); plot( ... - t_1,IpR_out_1(1,:)', ... - t_1,IpR_out_1(end,:)', 'linewidth', 4 ... + t_1, flip(IpR_out_1(1,:)'), ... + t_1, flip(IpR_out_1(end,:)'), 'linewidth', 4 ... ) xlabel('t (uS)') ylabel('I_{pR} (1/s)^2') @@ -111,8 +111,8 @@ tmin=-0.05; tmax=0.05; indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region plot( ... - t_1(indx),IpR_in_norm(indx), ... - t_1(indx),IpR_out_norm(indx), 'linewidth', 4 ... + t_1(indx),flip(IpR_in_norm(indx)), ... + t_1(indx),flip(IpR_out_norm(indx)), 'linewidth', 4 ... ) xlim([tmin,tmax]); xlabel('t (uS)') @@ -184,4 +184,3 @@ xlabel('z (cm)') ylabel('t (uS)') zlabel('Imag(rho_{ca})') title('Imag(rho_{ca})') - |