diff options
Diffstat (limited to 'xmds2')
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm/pp.m | 35 |
1 files changed, 26 insertions, 9 deletions
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m index ab536be..186d822 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m @@ -4,21 +4,27 @@ Nlevels_no_dopler_with_z z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS figure(1) -subplot(1,3,1); imagesc(z_1, t_1, I1_out_1); colorbar +subplot(2,2,1); imagesc(z_1, t_1, I1_out_1); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_1') title('I_1') -subplot(1,3,2); imagesc(z_1, t_1, I2_out_1); colorbar +subplot(2,2,2); imagesc(z_1, t_1, I2_out_1); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_2') title('I_2') -subplot(1,3,3); imagesc(z_1, t_1, I3_out_1); colorbar +subplot(2,2,3); imagesc(z_1, t_1, I3_out_1); colorbar xlabel('z (cm)') ylabel('t (uS)') zlabel('I_3') title('I_3') +subplot(2,2,4); imagesc(z_1, t_1, I4_out_1); colorbar +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_4') +title('I_4') + print('-color','fields_propagation.eps') @@ -26,7 +32,7 @@ print('-color','fields_propagation.eps') %% fields before and after the cell figure(2) -subplot(1,3,1); +subplot(2,2,1); plot( ... t_1,I1_out_1(:,1),'-;before;', "linewidth", 4, ... t_1,I1_out_1(:,end), '-;after;', "linewidth", 4 ... @@ -34,7 +40,7 @@ plot( ... xlabel('t (uS)') ylabel('I_1 (1/s)^2') title('I_1 before and after cell') -subplot(1,3,2); +subplot(2,2,2); plot( ... t_1,I2_out_1(:,1),'-;before;', "linewidth", 4, ... t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... @@ -42,7 +48,7 @@ plot( ... xlabel('t (uS)') ylabel('I_2 (1/s)^2') title('I_2 before and after cell') -subplot(1,3,3); +subplot(2,2,3); plot( ... t_1,I3_out_1(:,1),'-;before;', "linewidth", 4, ... t_1,I3_out_1(:,end), '-;after;', "linewidth", 4 ... @@ -60,9 +66,20 @@ printf('Second field delay time = %f uS\n',delay_time); print('-color','fields_before_after_cell.eps') +subplot(2,2,4); +plot( ... + t_1,I4_out_1(:,1),'-;before;', "linewidth", 4, ... + t_1,I4_out_1(:,end), '-;after;', "linewidth", 4 ... + ) +xlabel('t (uS)') +ylabel('I_3 (1/s)^2') +title('I_3 before and after cell') + figure(4) -I2_max_in=max(I2_out_1(:,1)); -I2_max_out=max(I2_out_1(:,end)); +settling_time=0.5; %uS +t_good_indx=t_1> min(t_1 + 0.5); +I2_max_in=max(I2_out_1(t_good_indx,1)); +I2_max_out=max(I2_out_1(t_good_indx,end)); I2_in_norm=(I2_out_1(:,1))/I2_max_in; I2_out_norm=(I2_out_1(:,end))/I2_max_out; tmin=-0.05; @@ -78,7 +95,7 @@ ylabel('I_2') title('I_2 before and after cell normalized') print('-color','probe_before_after_cell.eps') -%return; +return; %% all density matrix elements in one plot % diagonal populations, |