diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-11 15:25:42 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-11 15:25:42 -0400 |
| commit | 8549601a679147d5bd7f5416528b927f78a6b459 (patch) | |
| tree | a47793295131be8e4fe35e91654bdf67d2a2abc5 | |
| parent | 23b06abc5c734357cbd40ba4fcc61bb4f32f736d (diff) | |
| download | noisy_eit_xmds-8549601a679147d5bd7f5416528b927f78a6b459.tar.gz noisy_eit_xmds-8549601a679147d5bd7f5416528b927f78a6b459.zip | |
added MOR vs EIT comparison in postprocessing
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/pp_Nlevels.m | 49 |
1 files changed, 49 insertions, 0 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 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 |