summaryrefslogtreecommitdiff
path: root/MOR_5_levels_with_doppler_and_propagation
diff options
context:
space:
mode:
authorMelissa Guidry <maguidry@email.wm.edu>2016-04-27 15:41:46 -0400
committerMelissa Guidry <maguidry@email.wm.edu>2016-04-27 15:41:46 -0400
commit5e2a11994f8c3c3be6d9e0d394eeda51025ef96f (patch)
tree884b8fee257ef145194f126e688c2ad4b3418851 /MOR_5_levels_with_doppler_and_propagation
parent870e37ab2f8cc8c5c1d26a5ed294784431275371 (diff)
downloadnoisy_eit_xmds-5e2a11994f8c3c3be6d9e0d394eeda51025ef96f.tar.gz
noisy_eit_xmds-5e2a11994f8c3c3be6d9e0d394eeda51025ef96f.zip
Changed post-processing file to have proper time on y-axis.
Diffstat (limited to 'MOR_5_levels_with_doppler_and_propagation')
-rw-r--r--MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m45
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})')
-