diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m | 25 |
1 files changed, 19 insertions, 6 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 579925d..e107cb2 100644 --- a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m +++ b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m @@ -14,27 +14,40 @@ Epos = 1/sqrt(2)*(EpR + 1i*EpL); Eneg = 1/sqrt(2)*(EpR - 1i*EpL); Imor = abs(Epos).^2 - abs(Eneg).^2; +%% Maps are not so useful with too many pixels in the image +% Here we limit number of image points +Npx = 800; +Npy=600; + +Ndx = length(z_1); +Ndy = length(t_1); + +skip_xp = max(1, floor(Ndx/Npx) ); +skip_yp = max(1, floor(Ndy/Npy) ); +ind_x = 1:skip_xp:Ndx; +ind_y = 1:skip_yp:Ndy; + %% field propagation 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(z_1, (t_1), IdL_out_1.'); colorbar +subplot(2,2,1); imagesc(z_1(ind_x), t_1(ind_y), IdL_out_1(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{dL}') title('I_{dL}') -subplot(2,2,3); imagesc(z_1, (t_1), IpL_out_1.'); colorbar +subplot(2,2,3); imagesc(z_1(ind_x), t_1(ind_y), IpL_out_1(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{pL}') title('I_{pL}') -subplot(2,2,2); imagesc(z_1, (t_1), IdR_out_1.'); colorbar +subplot(2,2,2); imagesc(z_1(ind_x), t_1(ind_y), IdR_out_1(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{dR}') title('I_{dR}') -subplot(2,2,4); imagesc(z_1, (t_1), IpR_out_1.'); colorbar +subplot(2,2,4); imagesc(z_1(ind_x), t_1(ind_y), IpR_out_1(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{pR}') @@ -114,14 +127,14 @@ legend('before', 'after') %% subplot(2,2,2); -imagesc(z_1, (t_1), Ieit.'); colorbar +imagesc(z_1(ind_x), t_1(ind_y), Ieit(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{eit}') title('I_{eit}') subplot(2,2,4); -imagesc(z_1, (t_1), Imor.'); colorbar +imagesc(z_1(ind_x), t_1(ind_y), Imor(ind_x, ind_y).'); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_{mor}') |