From 8549601a679147d5bd7f5416528b927f78a6b459 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Wed, 11 May 2016 15:25:42 -0400 Subject: added MOR vs EIT comparison in postprocessing --- .../pp_Nlevels.m | 49 ++++++++++++++++++++++ 1 file changed, 49 insertions(+) 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 a5d61a8..9de8ccb 100644 --- a/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m +++ b/MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m @@ -6,6 +6,15 @@ IdR_out_1 = abs(EdRaR_1 +1i*EdRaI_1).^2; IpL_out_1 = abs(EpLaR_1 +1i*EpLaI_1).^2; IpR_out_1 = abs(EpRaR_1 +1i*EpRaI_1).^2; +EpR = EpRaR_1 + 1i*EpRaI_1; +EpL = EpLaR_1 + 1i*EpLaI_1; + +Ieit = abs(EpR).^2 + abs(EpL).^2; +Epos = 1/sqrt(2)*(EpR + 1i*EpL); +Eneg = 1/sqrt(2)*(EpR - 1i*EpL); +Imor = abs(Epos).^2 - abs(Eneg).^2; + + %% field propagation z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS @@ -78,6 +87,46 @@ plot( ... xlabel('t (uS)') ylabel('I_{pR} (1/s)^2') title('I_{pR} before and after cell') +legend('before', 'after') + +%% EIT and MOR fields before and after the cell +figure(6) +subplot(2,2,1); +plot( ... + t_1, (Ieit(1,:)'), ... + t_1, (Ieit(end,:)'),'LineWidth', 4 ... + ) +xlabel('t (uS)') +ylabel('I_{eit} (1/s)^2') +title('I_{eit} before and after cell') +legend('before', 'after') + +%% +subplot(2,2,3); +plot( ... + t_1, (Imor(1,:)'), ... + t_1, (Imor(end,:)'), 'linewidth', 4 ... + ) +xlabel('t (uS)') +ylabel('I_{mor} (1/s)^2') +title('I_{mor} before and after cell') +legend('before', 'after') + +%% +subplot(2,2,2); +imagesc(z_1, (t_1), Ieit.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_{eit}') +title('I_{eit}') + +subplot(2,2,4); +imagesc(z_1, (t_1), Imor.'); colorbar +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_{mor}') + +title('I_{mor}') %% return -- cgit v1.2.3